- Open Access
Unique organization and unprecedented diversity of the Bacteroides (Pseudobacteroides) cellulosolvens cellulosome system
Biotechnology for Biofuels volume 10, Article number: 211 (2017)
(Pseudo) Bacteroides cellulosolvens is an anaerobic, mesophilic, cellulolytic, cellulosome-producing clostridial bacterium capable of utilizing cellulose and cellobiose as carbon sources. Recently, we sequenced the B. cellulosolvens genome, and subsequent comprehensive bioinformatic analysis, herein reported, revealed an unprecedented number of cellulosome-related components, including 78 cohesin modules scattered among 31 scaffoldins and more than 200 dockerin-bearing ORFs. In terms of numbers, the B. cellulosolvens cellulosome system represents the most intricate, compositionally diverse cellulosome system yet known in nature.
The organization of the B. cellulosolvens cellulosome is unique compared to previously described cellulosome systems. In contrast to all other known cellulosomes, the cohesin types are reversed for all scaffoldins i.e., the type II cohesins are located on the enzyme-integrating primary scaffoldin, whereas the type I cohesins are located on the anchoring scaffoldins. Many of the type II dockerin-bearing ORFs include X60 modules, which are known to stabilize type II cohesin–dockerin interactions. In the present work, we focused on revealing the architectural arrangement of cellulosome structure in this bacterium by examining numerous interactions between the various cohesin and dockerin modules. In total, we cloned and expressed 43 representative cohesins and 27 dockerins. The results revealed various possible architectures of cell-anchored and cell-free cellulosomes, which serve to assemble distinctive cellulosome types via three distinct cohesin–dockerin specificities: type I, type II, and a novel-type designated R (distinct from type III interactions, predominant in ruminococcal cellulosomes).
The results of this study provide novel insight into the architecture and function of the most intricate and extensive cellulosomal system known today, thereby extending significantly our overall knowledge base of cellulosome systems and their components. The robust cellulosome system of B. cellulosolvens, with its unique binding specificities and reversal of cohesin–dockerin types, has served to amend our view of the cellulosome paradigm. Revealing new cellulosomal interactions and arrangements is critical for designing high-efficiency artificial cellulosomes for conversion of plant-derived cellulosic biomass towards improved production of biofuels.
Cellulosic biomass and waste are raw materials of great abundance, and its deconstruction conversion to soluble sugars is an important resource within the context of production of biofuels and valuable chemicals [1, 2]. Some anaerobic cellulolytic bacterial strains have developed the cellulosome, an efficient enzymatic strategy to utilize cellulosic biomass as a major carbon source. One of the major advantages of cellulosome-producing bacteria is their ability to degrade different types of carbohydrates present in various types of biomass . The organization of enzymes into a cellulosome serves to concentrate them physically and position them in suitable orientation, both with respect to each other and to the cellulosic substrate, for efficient decomposition of the recalcitrant insoluble substrate . Moreover, the fact that the complex is both attached to the substrate and to the cell results in minimal diffusion loss of enzymes and hydrolytic products, and precludes product-mediated feedback inhibition of the cellulolytic enzymes. The cellulosomal enzymes are incorporated into the complex via their resident dockerin module and interact specifically with the cohesin modules of a structural scaffoldin subunit [3,4,5,6]. The scaffoldin subunit can selectively integrate enzymes or additional scaffoldin subunits into the cohesive complex via specific and high-affinity cohesin–dockerin interactions, which thus determine overall cellulosome architecture [7,8,9,10].
Cohesins and dockerins have been classified traditionally into types (I, II, and III) based on sequence similarity [3, 11, 12]. Primary scaffoldins, the backbone of the cellulosomal complex, have thus far been demonstrated to contain multiple type I cohesins, each of which interacts with a type I dockerin harbored by each cellulosomal enzyme [13,14,15]. The primary scaffoldin may contain a dockerin module that interacts with the cohesins of an adaptor and/or anchoring scaffoldin, thereby allowing the expansion of the cellulosome complex by integration of multiple enzymes and/or the attachment of the cellulosome to the cell surface [16,17,18]. These scaffoldin assemblies are generally mediated by type II cohesins and dockerins that are located on the adaptor or anchoring scaffoldins [9, 18]. The anchoring scaffoldins can contain one or more cohesins and anchor the cellulosome complexes to the cell surface via a surface-layer homology (SLH) domain [19, 20]. The cohesin–dockerin interactions are considered to be species- and/or type-specific, although some cross-species interactions have been observed .
In this context, (Pseudo) Bacteroides cellulosolvens is an anaerobic, mesophilic, cellulolytic bacterium that was isolated from a methanogenic cellulose-enrichment culture of municipal sewage sludge . This bacterium produces an extracellular multi-enzyme cellulosome complex for efficient degradation of plant cell wall polysaccharides and cellulosic wastes  and is capable of utilizing cellulose or cellobiose as a sole carbon source . Originally termed Bacteroides cellulosolvens, the bacterium was subsequently found to be phylogenetically related to the clostridial assemblage  and more recently reclassified as Pseudobacteroides cellulosolvens . Earlier work reported two major scaffoldins in B. cellulosolvens  and the cellulolytic potential of the bacterium [22, 27]. The two proteins, a primary scaffoldin and an anchoring scaffoldin, were the largest yet described, bearing 11 and 10 cohesins, respectively [28, 29]. Recently, the B. cellulosolvens genome was sequenced to near-completion  allowing comprehensive bioinformatic studies that will represent a milestone in current research on this bacterium. Therefore, in this work, we explored the architectural and functional aspects of the cellulosome of B. cellulosolvens, and in particular the cohesin–dockerin specificities of interactions between different scaffoldin and enzymatic modules. Its large range of cellulosomal components was revealed, and we demonstrated binding activity and specificity of selected cohesin and dockerin modules, thus revealing overall cellulosome architecture in this intriguing cellulosome-producing bacterium.
Anaerobic fermentation of Bacteroides cellulosolvens
Bacteroides cellulosolvens ATCC 35603 was grown under anaerobic conditions essentially as described by Murray et al.  with either cellobiose (CB, Sigma Chem. Co. St. Louis, MO) or microcrystalline cellulose (MCC, Avicel, E. Merck, Darmstadt, Germany) as carbon and energy source. B. cellulosolvens cell lysates were prepared using PopCulture Reagent (Novagen Inc, Darmstadt, Germany), as described by Slutzki et al. .
Fractionation of high-molecular-weight complexes
The spent growth medium of B. cellulosolvens cells, grown on either CB or MCC, was concentrated 100-fold and subjected to gel-filtration chromatography on a Superose 6 gel-filtration column (GE Healthcare) as described earlier . The two resultant peaks (I and II) were pooled and concentrated using a Vivaspin concentrator (50-kDa cutoff; Sartorius Stedim Biotech GmbH, Göttingen, Germany).
Blastp searches were performed against predicted B. cellulosolvens proteins, using deduced amino acid sequences of the known cohesin and dockerin modules as queries [16, 17, 33]. Hits above an E-value of 10−4 were examined individually, by searching for characteristic sequence features. For example, for dockerin modules, we searched for two Ca2+-binding repeats, putative helices and linker regions. Multiple sequence alignments were created using the Clustal Omega server [http://www.ebi.ac.uk/Tools/msa/clustalo/]. Phylogenetic trees were generated by iTOL version 3 [http://itol.embl.de/] according to the “One Click” Phylogeny analysis tool [http://www.phylogeny.fr/simple_phylogeny.cgi]. Signal peptide sequences were predicted using the SignalP server [http://www.cbs.dtu.dk/services/SignalP/]. Amino acid sequence logos were performed using the WebLogo3 application, version 3.5.
Annotation of dockerin-containing enzymes
The proteins were annotated using the carbohydrate-active enzymes database (CAZy) http://www.cazy.org/ . The analysis was based on sequence conservation between catalytic modules, and the different catalytic modules were sorted into different families.
Cloning and expression plasmid cassettes
The XynDoc gene cassette consists of xylanase T6 from Geobacillus stearothermophilus with an N-terminal His-tag cloned into plasmid pET9d (Novagen Inc., Madison, WI, USA), into which a dockerin-encoding sequence was introduced between the KpnI and BamHI restriction sites of the plasmid . The CBM-Coh gene cassette consists of a family CBM3 (family 3 carbohydrate-binding module) from the Clostridium thermocellum CipA scaffoldin cloned into plasmid pET28a (Novagen Inc., Madison, WI, USA), into which a cohesin gene was introduced between BamHI and XhoI restriction sites of the plasmid [35, 36].
Polymerase chain reaction (PCR)
An expanded high-fidelity PCR system (Boehringer Mannheim) was used in all PCRs. PCR was performed using a Mastercycler personal instrument (Eppendorf, Hamburg, Germany), programed as follows: a 3-min predenaturation step at 95 °C was followed by 30 cycles comprising a 45-s denaturation step at 94 °C, an annealing step of 30 s at 50–60 °C (depending on the primer), and an extension step at 72 °C for 1 min. The primers used for the cloning of 43 cohesins and 27 dockerins are listed in Additional file 1: Table S1.
PCR products were purified and double digested at 37 °C for 15–30 min with FastDigest restriction enzymes (Thermo Scientific) and ligated into the desired plasmid. Positive clones were verified by sequencing.
The pET28a cassette containing the CBM-Coh fusion proteins and the pET9d cassette containing the XynDoc fusion proteins were transformed into Escherichia coli BL21 (DE3) strains and plated onto LB-kanamycin plates. For each plate, 4–5 ml of Luria–Bertani broth (LB) were added in order to resuspend the cells. The resuspended cells were added to 1 l of LB with 50 µg/ml kanamycin and 2 mM CaCl2 and were grown for 2 h at 37 °C to A600 ≈ 0.8–1. Protein expression was induced by adding isopropyl-1-thio-β-D-galactoside (IPTG) (Fermentas UAB, Vilnius, Lithuania) in a final concentration of 0.2 mM, and the growth was continued in 16 °C for 16 h. Cells were harvested by centrifugation at 5000 rpm for 15 min.
Purification of CBM-containing proteins
The supernatant fluids of the cohesin-containing proteins (fused to a CBM tag, both for increased solubility and for affinity purification) were added to 2 g of preswollen cellulose gel macroporous beads (IONTOSORB, Usti nad Labem, Czech Republic) and incubated for 1 h with rotation at 4 °C. The mixture was then loaded onto a column, and washed with 100 ml of Tris-buffered saline (TBS: 13.7 mM NaCl, 0.27 mM KCl, 2.5 mM Tris, pH 7.4) brought to 1 M NaCl, and then washed with 100 ml TBS. Three 5-ml elutions of 1% triethanolamine (TEA) were then collected, protease-inhibitor cocktail was added. The fractions were subjected to SDS-PAGE in order to assess protein purity.
Purification of Xyn-containing and His-tagged proteins
The supernatant fluids containing the dockerin-bearing proteins were mixed with 4 ml Ni–NTA, for 1 h on a 20-ml Econo-pack column, on a rotator at 4 °C (batch purification system). The column was then washed by gravity flow with 100 ml wash buffer (TBS, 15 mM imidazole). Elution was performed first using 100 mM imidazole, followed by 250 mM imidazole. Fractions (2 ml) were collected and were run on SDS-PAGE. The fractions containing relatively pure proteins were pooled, and CaCl2 (10 mM), as well as protease-inhibitor cocktail was added.
Protein concentration and storage
Protein concentration was evaluated by absorbance at 280 nm, based on the extinction coefficients derived from the known composition of amino acids of each protein. Extinction coefficients were calculated using the ExPASy ProtParam tool http://web.expasy.org/protparam/. Some proteins were concentrated by Amicon ultra concentrators (Millipore, Ireland), and stored at −20 °C in 50% (vol/vol) glycerol.
ELISA-based affinity assay
The standard ELISA procedure was performed as described previously . Representative cohesin and dockerin modules were selected and expressed using one of the two cassettes described above. In this manner, we cloned 43 CBM-fused cohesins and 27 as Xyn-fused dockerins (13 from the scaffoldins and 14 from the putative enzymes). The 96-well ELISA plates (Nunc, A/S, Roskilde, Denmark) were coated with the fusion proteins CBM-Cohs or full-length scaffoldins at a concentration of 1–10 µg/ml, and variable concentrations of Xyn–Docs (0.001–1000 ng/ml) were used to detect specific cohesin–dockerin interactions. Interactions with the Xyn–Doc fusion proteins were examined immunochemically by using anti-xylanase primary antibody and HRP-labeled secondary antibody. The experiments were performed three times in duplicate.
For cell lysate-based ELISA, the 96-well ELISA plates were coated with cellobiose-grown B. cellulosolvens cell lysate, and graded concentrations of the desired Xyn–Docs were used to examine cohesin–dockerin interactions.
Absorbance was plotted as a function of Xyn–Doc fusion proteins concentration. For comparative purposes, the reference concentration of a Xyn–Doc standard that generates a maximum response was used in order to normalize the data as a relative binding of maximum response, as described earlier . The results were presented as a heatmap (iTOL version 3, http://itol.embl.de/), whereby each node is associated with multiple numerical values, which are displayed as a set of colored boxes. Dataset values are mapped to a color gradient corresponding to the binding strength.
Xylan activity assay was performed in triplicate in a total volume of 500 µl, containing 50 mM citrate buffer (pH 6.5), 12 mM CaCl2, 2 mM EDTA, and 25 µg/ml of purified cellulosome complex from B. cellulosolvens. Xylan degradation was assayed at a final concentration of 1% beechwood xylan (Sigma-Aldrich, Rehovot, Israel), for 1 h at 42 °C (according to predetermined optimal conditions for B. cellulosolvens cellulosome activity). The assay performed for the purified C. thermocellum cellulosome was incubated at 70 °C (the optimal temperature for C. thermocellum cellulosome activity). The tubes were incubated under shaking (400 rpm), and the reaction was terminated by flash-cooling the tubes on ice. The tubes were centrifuged (14,000 rpm, 5 min), and 100 µl of the supernatant was transferred into 150 µl dinitrosalicylic acid (DNS) solution. The tubes were boiled for 10 min at 100 °C, and absorbance was measured at 540 nm in a plate reader. A glucose standard curve served to determine the amount of reducing sugars (in mM).
The elaborate cellulosomal system of B. cellulosolvens revealed by bioinformatics
We have recently sequenced the near-complete genome of B. cellulosolvens DSM 2933 (ATCC 35603), which appears to be the largest among the currently known cellulolytic bacteria (~6.9 Mbp) (Fig. 1) . Detailed bioinformatics analysis revealed multiple cellulosomal components. In fact, this bacterium contains the largest number of cellulosomal components currently known. To delineate the mode by which the components may assemble into cellulosomal complexes remains an intriguing assignment . We herein revealed 78 cohesin modules, scattered among 31 scaffoldins, and 212 dockerin-bearing ORFs, representing 197 putative carbohydrate-degrading enzymes [including assorted glycoside hydrolases (GHs), carbohydrate-binding modules (CBMs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), and defined X-modules], and 15 dockerin-bearing scaffoldins (Fig. 2a). Almost half of the enzyme-borne type II dockerins (92 out of 212) possess an X60 module upstream of the dockerin sequence. As noted earlier for the then-discovered isolated B. cellulosolvens components [28, 29], in comparison to previously described cellulosome systems, the apparent roles of the B. cellulosolvens cohesins are curiously reversed, compared to all previously described cellulosomal components, in that the type II cohesins are located on the enzyme-binding primary scaffoldin, whereas the type I cohesins are located on the anchoring scaffoldin. In addition, significant numbers (17) of scaffoldin genes were found to be arranged in genomic clusters (Fig. 2b), whereas dockerin-containing genes were scattered more evenly throughout the genome (Fig. 1).
Diversity of CAZy-associated cellulosomal enzymes
Bacteroides cellulosolvens was found to contain three times more dockerin-bearing proteins, as compared to other clostridia, such as Clostridium cellulolyticum (~60 dockerins), C. thermocellum (>70 dockerins), or Clostridium clariflavum (79 dockerins) [17, 37, 38]. The Acidothermus cellulolyticus genome contains 143 dockerin-containing ORFs . The number (212) of dockerin-bearing ORFs in the B. cellulosolvens genome, however, is more comparable to those of Ruminococcus flavefaciens strains FD-1, 17, and 007c, which contain between 180 and 223 dockerins [33, 39, 40]. Table 1 presents the abundance of CAZy-associated modules (cellulosomal and non-cellulosomal) in the B. cellulosolvens genome. In general, about 50% of the B. cellulosolvens dockerins are associated with carbohydrate-active enzymes (GH, PL, CE). About 85 out of 212 dockerin-containing proteins were not associated with a defined CAZy module (Table 2).
The GH48 enzymes are known to be the definitive exoglucanase and quantitatively most abundant enzyme type, in all known cellulosomes [41, 42]. Remarkably, the B. cellulosolvens cellulosome contains three distinct GH48 enzymes in contrast to A. cellulolyticus, C. thermocellum, and C. clariflavum that contain only a single GH48 cellulosomal enzyme [17, 41]. As opposed to other known clostridial species that possess only one GH48 cellulosomal enzyme, in B. cellulosolvens all three GH48 possess a dockerin. The cellulase systems of other complex cellulosome-producing clostridia mentioned above also contain a second GH48 enzyme, but it bears a CBM3 rather than a dockerin and is thus not cellulosomal. One of the B. cellulosolvens GH48 enzymes (WP_050753099.1) was shown previously to bind a ScaA1 cohesin .
As reported earlier [22, 43, 44], B. cellulosolvens grows on cellulose and cellobiose as sole carbon sources but is also able to degrade xylan and has high xylanase activity in secreted and cell-associated fractions (Additional file 2: Figure S1). Recently, it was shown that B. cellulosolvens is a highly active lignocellulolytic microorganism able to efficiently digest cellulose, hemicellulose, and lignin together with Clostridium stercorarium . Here, CAZy analysis revealed 147 GH modules, with a wide array of cellulolytic and hemicellulolytic enzymes, either cellulosomal (about 60% containing dockerins) or free enzymes (40%). Many of the cellulosomal and free enzymes (94 and 60, respectively) possess a CBM in addition to the definitive catalytic module(s), and in some cases more than one, thereby enabling extensive interaction with the lignocellulosic substrate.
Thirteen GH families in the genome are non-cellulosomal, suggesting that they could support biomass degradation of cellulose in free form or contribute to the degradation of distant or concealed carbohydrates. Interestingly, the percentage of dockerin-containing GHs (60%) in the B. cellulosolvens genome is similar to that of C. thermocellum and A. cellulolyticus, despite the superior number of enzymes in B. cellulosolvens. The GH9 family includes the highest number of enzymes, similar to the other cellulosome-producing bacteria. Among the 40 GH9s, 33 contain a dockerin. This is the most abundant GH family known today among the clostridia. In addition to dockerins, most of the cellulosomal GH9 enzymes contain CBM3 modules: fifteen CBM3c-possessing enzymes (two of them containing two CBM3s), fourteen CBM4-possessing enzymes. and two with CBM30s, in accordance with known modular architectures the CBM modules, would be expected to provide a significant contribution to enzyme action [46,47,48]. The wealth of the GH9 family in this bacterium indicates its important role in biomass degradation by the ability of its members to bind and hydrolyze cellulosic and xylan/xyloglucan substrates .
GH10 is the second most abundant family with 15 enzymes (Table 1), four of which are multifunctional enzymes together with an additional GH motif and a dockerin (Table 3). Four non-cellulosomal GH10 members are associated with CE4 and CBM22/CBM9 elements, and three of them contain a triple SLH repeat, suggesting that these enzymes are attached to the cell surface. Four of seven cellulosomal GH10 enzymes contain a CBM6, suggesting strong cellulose binding .
Eleven enzymes were revealed as containing GH5 and GH43 modules, most of them cellulosomal (Table 1). In C. thermocellum and A. cellulolyticus, GH5 is the second most abundant family, but in B. cellulosolvens the representation is somewhat different. Similar to other bacteria, CBM3 is prevalent in the B. cellulosolvens genome with 30 representative proteins, 16 of which possess a dockerin module, mostly associated with GH and CE enzymes. A. cellulolyticus has 24 CBM3 members, 19 of which are associated with a GH9 . In B. cellulosolvens, out of 30 CBM3s, 16 are associated with GH9 enzymes and 15 of them appear to be cellulosomal. CBM4 and CBM6 are less abundant, but represent a considerable part of the CBM family, with 23 and 20 members, respectively; 14 of these proteins in each group are associated with a dockerin. In general, we observed a very wide array of enzymatic and structural modules, which may collectively assemble into a robust machinery of both cellulosomal and free biomass-degrading components.
The variety of B. cellulosolvens GH catalytic modular representatives emphasizes the robustness of its cellulosome system. Another intriguing feature is the presence of 17 multifunctional enzymes (Table 3), which harbor a combination of at least two catalytic modules in the same polypeptide. Seven of these enzymes include two different GH families, and two have two catalytic modules from the same family—GH16 (GH16-GH16-Doc, KNY27855.1) and GH43 (GH43-Doc-CBM42-GH43, KNY29222.1), respectively (Table 3). There are also several mixed bifunctional hemicellulase/carbohydrate esterases, a dual carbohydrate esterase and a bifunctional polysaccharide lyase. Similar types of multifunctional protein architectures have been reported in C. thermocellum , A. cellulolyticus , R. flavefaciens , Ruminococcus champanellensis , and other bacteria , indicating that multifunctional enzymes are a common component in cellulosomal systems.
A high number of dockerins was associated with putative peptidases (71 proteins in total), suggesting a broader role of the cellulosomal complex in parallel with fiber degradation. Predicted peptidase modules were also found in scaffoldins ScaO and ScaP of B. cellulosolvens (associated with multiple PPC modules in addition to a single cohesin and dockerin), as well as in ScaO and ScaP of A. cellulolyticus . The role of cellulosomal peptidases has not been defined experimentally, but recent studies have suggested the presence of peptidase modules associated with dockerins in the metagenome of the bovine rumen . Similarly, in R. flavefaciens FD-1, numerous dockerins are associated with putative peptidase modules . One putative cysteine peptidase associated with a C-terminal X-dockerin modular dyad from R. flavefaciens exhibited functional binding to the surface-anchoring ScaE cohesin .
Characterization of the numerous scaffoldins and cohesins
Only two scaffoldins were previously reported in B. cellulosolvens [26, 28]. The present work revealed an unprecedented number (29) of additional cohesin-containing scaffoldins for a total of 31 B. cellulosolvens scaffoldins. Figure 2a presents the modular architecture of all the putative B. cellulosolvens scaffoldin proteins and their diverse types of cohesin and dockerin components. All scaffoldins (except ScaA2) contain a predicted signal peptide , suggesting that these proteins are secreted. Proteomics experiments indicated the presence of ScaA2 in the spent supernatant fluids of B. cellulosolvens growth cultures (data not shown), indicating that this scaffoldin was also secreted, despite the apparent lack of a signal peptide.
In naming the different B. cellulosolvens scaffoldins, we tried to compare their predicted architecture with those of previously described scaffoldins from other cellulosome-producing bacteria. Scaffoldins A to P (19 scaffoldins) contain cohesins and dockerins possessing modular arrangements similar to those of other known bacteria, particularly to those of Acetivibrio cellulolyticus, but with one important difference—the cohesin types are always reversed, i.e., if the primary cohesins of the homologous scaffoldins of A. cellulolyticus (and other species) are of type I, then those of B. cellulosolvens will be of type II and vice versa . Remarkably, we observe this pattern in all B. cellulosolvens cellulosomal proteins that have orthologues in other cellulosome-producing bacteria.
The cohesin modules within the scaffoldins exhibit a variety of intriguing sequence features. This bacterium also has some unique cohesin sequences which are somewhat different from the canonical type I or type II classification, according to the majority of known cellulosomal systems [21, 53]. Multiple sequence alignment of the cohesins can be found in Additional file 3: Figure S2 (the file includes scaffoldin accession numbers). Of the various B. cellulosolvens cohesins, 75 are classified into the two main types: type I (33 modules) and type II (42 modules). In addition to the canonical cohesin types I and II (and type III of the ruminococci), three B. cellulosolvens scaffoldins (ScaR1, ScaR2, and ScaR3) represented significantly different cohesin and dockerin sequences that exhibited only weak similarity to the main types and were therefore classified as ‘group R.’ We then examined the conservation patterns of the cohesin sequences, both within and among the different scaffoldins (Fig. 3). Clustered scaffoldins (Fig. 2b) may share homologous cohesins of similar types (ScaQ, ScaX1, ScaX2, and ScaD), although more distant cohesins may share some similarity as well. Two adjacent ORFs (scaffoldins ScaA1 and ScaB, Fig. 2) include different cohesin types, similar to the ScaA and ScaB pairs observed in Clostridium thermocellum, A. cellulolyticus, and C. clariflavum , albeit, as noted above, reversed in type.
Scaffoldins Q to X were found to share less similarity to the known scaffoldins from other bacteria and were named alphabetically taking into account intraspecies similarity (Fig. 2a). Overall, the modular organization of the proteins in B. cellulosolvens bear similarities to A. cellulolyticus and C. thermocellum, but the bacterium also contains new types of scaffoldins which were not described before (ScaQ–ScaX2). Most of the multiple cohesins within a scaffoldin are very similar, e.g., ScaA1 and ScaM1, but the phylogenetic tree also reveals variability among cohesin sequences, even within a single scaffoldin (Fig. 3). For instance, ScaA2 cohesins 6 and 10 are significantly distant from each other, although they are all classified as type II.
The numbers of the cohesins on the scaffoldins vary from a single cohesin (20 different scaffoldins) to 11 cohesins on ScaA1, the largest number of cohesin modules found on a single scaffoldin to date. ScaB, the adjacent downstream ORF of ScaA1, is an anchoring scaffoldin with 10 type I cohesins and an S-layer homology (SLH) domain, which is believed to form a non-covalent interaction with peptidoglycan-associated polymers to attach the protein to the cell surface . These two largest ORFs are clustered on the genome, resembling the clusters described in other cellulosome-producing species, notably C. thermocellum , A. cellulolyticus [15, 18], and R. flavefaciens . According to the bioinformatics analysis of regulatory regions flanking the scaA and scaB genes , it is not likely that scaA1 and scaB of B. cellulosolvens are transcribed together. SLH domains and cell surface-binding modules (CSBMs) enable attachment of the anchoring scaffoldin to the bacterial cell surface [20, 56] and are present in nine scaffoldins (Fig. 2a). For some B. cellulosolvens scaffoldins (ScaB, ScaF2, and ScaJ), we observed the presence of an interesting SLH domain structure that includes a unique type of X-module at the N terminus, which differs from the X60 module of the ScaA subunit, known to stabilize dockerin interactions in C. thermocellum, A. cellulolyticus, and C. clariflavum . The X-SLH modular dyad in scaffoldin proteins thus far seems to be unique to B. cellulosolvens. In addition to SLH and CSBM domains, we have identified other domains that may participate in anchoring of the scaffoldin to the cell wall or interactions with the substrate, including a PA14 domain of ScaU, a cadherin domain in ScaD (in addition to the SLH), VCBS in ScaH3, and a PPC (pre-peptidase C-terminal) domain in ScaO and ScaP (Fig. 2a).
The B. cellulosolvens genome contains genes for seven scaffoldins with no dockerin, CSBM, or SLH domain, which implies that they may serve as cell-free scaffoldins, which contain either type I or type II cohesins (Fig. 2a). Two examples of this type of scaffoldin, ScaE and ScaS, bear type I cohesins, potentially forming cell-free cellulosomes with up to 77 dockerin-bearing enzymes (Fig. 4). Intriguingly, C. thermocellum, A. cellulolyticus, and C. clariflavum all produce ScaE homologues, bearing seven type II cohesins in these species (for binding primary scaffoldins, reverse in type, compared to B. cellulosolvens). In B. cellulosolvens, enzyme-binding type II cohesins from the cell-free primary scaffoldins ScaM1 and ScaM2 have a CBM3, which would bind the substrate and thus allow targeted degradation by the enzymes. Enzyme-integrating CBM-bearing ScaM homologues (with CBM2s) have been detected in the cellulosome systems of A. cellulolyticus and C. clariflavum but not C. thermocellum.
Classification of the dockerins into types based on sequence homology
In B. cellulosolvens all of the type I dockerins are associated with primary scaffoldins (Fig. 2a; Additional file 4: Figure S3). Scaffoldins ScaR1, ScaR2, and ScaR3 possess a unique type of dockerin that did not fit the main types and were collectively termed ‘group R.’ Similarly, the ScaP dockerin with its unique sequence was not included into any of the main dockerin types (Additional file 4: Figure S3).
Dockerin modules are characterized by two duplicated segments, consisting of a calcium-binding loop that precedes an α helix (Fig. 5), connected by a linker sequence, with distinctive N- and C-terminal stretches . Dockerins usually display a conserved pattern within the given type. Type I dockerins would presumably bind type I cohesins within the same bacterial species, and the same applies to types II .
Bacteroides cellulosolvens dockerins were herein classified into two major previously defined groups (type I and type II) and a new group (group R that includes ScaR1-3 dockerins). Almost half of the dockerins are located downstream of an X-module and have distinctive sequence features compared to the rest of the B. cellulosolvens dockerins. Their X-modules belong to family X60, which displays significant sequence identity (30–57%) with the X-module at the C-terminus of the C. thermocellum CipA scaffoldin . As mentioned above, these X-modules are known to stabilize their adjacent dockerin and render them more soluble . In addition, X60-modules were described at the C-terminus of the primary scaffoldin of A. cellulolyticus and C. clariflavum, all related to type II dockerins [16, 17]. The high number of X60 modules in B. cellulosolvens may indicate their crucial role in the dockerin interactions. Intriguingly, all of the type II dockerins described previously were accompanied by a neighboring X-module [16, 17, 57], but here, in B. cellulosolvens, the presence of many proteins containing type II dockerins (111) lacking the adjacent X-module remains an enigma. The sequence alignment of the type II dockerins is presented in Additional file 5: Figure S4 and Additional file 6: Figure S5.
Prior to the sequencing of the B. cellulosolvens genome, the scientific community was cognizant of only a few type II dockerin and cohesin sequences. Only the type II dockerins of the primary scaffoldins of C. thermocellum [11, 61], A. cellulolyticus , and C. clariflavum , and the early discovery of the B. cellulosolvens Cel48 dockerin  had been reported. The 197 type II dockerins revealed by this genome has thus significantly enriched our understanding of this very basic dockerin type.
The characteristic sequence conservation profile of the B. cellulosolvens type II dockerin is shown in Fig. 5 [21, 53, 59]. Examination of the putative recognition residues revealed the highly conserved calcium-coordinating residues: Asp in positions 1 and 12; mostly Asn/Asp in position 3; Asn/Asp in position 5; and Asn in position 9. Positions 3, 5, and 9 are infrequently replaced by Ser, Thr, and sometimes Lys. Other variations have also been observed [24, 25, 31]. The putative calcium-binding residues are consistent with those of type II dockerins known from the literature , and the predicted recognition residues show Met and Ala dominating in positions 10 and 17, respectively. Interestingly, Phe dominates at position 20 of the helix and Asn and Gly are prevalent at position 21, consistent with the few previously known type II dockerin sequences.
Selection of cohesins and dockerins
To shed light on cellulosome assembly in this unique bacterium, multiple cohesins and dockerins were selected for further experimental investigation. A total of 43 cohesins and 27 dockerins were cloned as fusion proteins with a solubility/stability tag, expressed and purified. The solubility tags, CBM3 from C. thermocellum CipA and xylanase T6 from G. stearothermophilus for the cohesin and dockerin modules, respectively, also served as general affinity tags for semiquantitative detection of the cohesin–dockerin interaction by specific antibodies directed against them using an ELISA-based system .
The initial alignment of the cohesin (Additional file 3: Figure S2) and dockerin sequences (Additional files 4, 5, 6: Figures S3–S5) served to determine their type distribution and was used for the selection of cohesins and dockerins for biochemical characterization. Despite the very large number of modules in this bacterium, we tested experimentally at least one cohesin from nearly all of the scaffoldins as well as the most divergent cohesins within a given multiple-cohesin scaffoldin (Fig. 2a). Only three scaffoldins: ScaW1, ScaW2, and ScaG that were discovered by scrutinizing the genome at a later stage of the study were not tested. Almost all of the scaffoldin-borne dockerins (except ScaO and ScaH2) were selected, owing to their important role in cellulosome assembly. In addition, we selected 17 dockerins that diverged in their sequences from the most prevalent GH modules (e.g. GH5, GH8, GH9, GH10, GH30, GH43, and GH48). When an X-module was adjacent to the selected dockerin, it was included in the cloned sequence.
Following expression and purification of these 70 proteins, SDS-PAGE analysis revealed single major protein bands in each case, in agreement with their calculated molecular mass.
Identification of cohesin–dockerin interactions
In total, 103 positive interactions were detected (Fig. 4). The specificity of the various cohesin and dockerin counterparts revealed in this study served to determine the theoretical supramolecular organization of its known cellulosomal components. The cellulosomal architectures of B. cellulosolvens cellulosome are represented in Fig. 6. Our analyses underscore the highly heterogeneous and diverse supramolecular architecture of this cellulosome system.
In accordance with previous reports [28, 29], our results indicated that the type I ScaB cohesins bind selectively to the ScaA1 dockerin, whereas the GH48 (WP_050753099.1) dockerin binds specifically to the type II ScaA1 cohesins (Fig. 4). The cohesins of ScaA1 also bound to various cellulosomal enzymes, in particular to GH10 (WP_036936763.1), GH43 (WP_081926929.1), and GH30 (WP_081927211.1). In addition, cohesins of other primary scaffoldins (namely, ScaA2, ScaH1, ScaH2, ScaJ, ScaL2 ScaM2, and ScaO) shared the same binding specificities as ScaA1 (Fig. 4). Type II cohesins from scaffoldins ScaL1 and ScaV showed clear preference to bind type II dockerins but with lower intensity. The type II dockerin of ScaI failed to show any interaction, but its type II cohesin showed low levels of interaction with the ScaH1 dockerin.
Some of the tested enzyme-borne dockerins, i.e., two dockerin modules from GH9 enzymes (WP_050753192.1; KNY25939.1) and a GH5 with an X60-dockerin modular dyad (WP_050753119.1), failed to show binding specificity to ScaA1 cohesins or those of other primary scaffoldins. In addition, six of the selected dockerins failed to interact with any of the selected cohesins: one originated from an ORF containing only two dockerin modules (WP_036940956.1), another from an ORF containing a putative peptidase and two similar dockerins (WP_036945116.1), one GH30-associated dockerin along with its X60 module (KNY28903.1), dockerin modules from scaffoldins ScaI and ScaT along with their adjacent X60 and an X60-dockerin pair originating from a GH43 module (KNY26505.1). The dockerin sequence of the latter is similar to another dockerin from a GH43 enzyme (WP_038290784.1) that exhibited high binding interaction (Fig. 4). In this particular case, we designed two Xyn–Doc fusions with or without the X-module and expressed the full-length enzyme but none of those constructs exhibited any binding activity.
The recent sequencing of the B. cellulosolvens genome enabled comprehensive bioinformatics identification of the numerous cellulosomal components and cell-anchoring modules. The high quality of the genome sequence  allowed us to identify an unprecedented number of scaffoldins in this bacterium. The cohesin and dockerin modules contain some unique and intriguing sequences, which were separated on the basis of bioinformatics and complementary biochemical analysis into the two major conventional types and one novel group—group R, which contain only a few members.
The cellulosomal system of B. cellulosolvens represents the most complex so far discovered in nature, by virtue of its 31 different scaffoldins—nearly four times that of the C. thermocellum standard (i.e., initially discovered cellulosome system) and double that of A. cellulolyticus, the next most extensive known system. Theoretically, B. cellulosolvens can assemble up to 110 enzymes in a single, cell-associated complex (Fig. 6), consisting of only two interacting scaffoldins, ScaA1 and ScaB, without the aid of an adaptor scaffoldin . In addition, its genome contains numerous free scaffoldins (which lack any apparent cell surface-binding modules) that interacted specifically with selected dockerins. This putative cell-free cellulosome system may serve an important role of degrading carbohydrates distant from the bacterial cell, as hypothesized earlier . Free cellulosomes have been described before for several cellulosome-producing bacteria and are believed to contribute to a more efficient carbohydrate degradation [17, 19]. In this context, ScaE (which contains seven type I cohesins) and ScaS (which contains three type I cohesins) were shown to bind primary scaffoldins (Fig. 4) and could assemble into comparatively large free cellulosome complexes (Fig. 6). The additional cell-free scaffoldins, ScaM1 and ScaM2, contain type II cohesins along with a CBM3 module that would target the attached enzymes to the substrate. In contrast, monovalent (single cohesin) free scaffoldins (ScaW1, ScaW2) may either serve as molecular shuttles, stabilize enzyme activity through cohesin–dockerin interaction, or serve a regulatory function [63,64,65]. Multiple monovalent scaffoldins are widespread among complex cellulosome-producing bacteria [16, 19, 66] and could be part of a regulatory mechanism for cellulosomal composition.
Previous research suggested that a dual-binding mode of the dockerins would result in increased flexibility characteristics of the catalytic subunits [67, 68], it was reported previously that the type I ScaA dockerin has a dual-binding mode in the recognition of the B. cellulosolvens type I cohesins . The additional type I dockerins reported here share similar binding preferences and sequence similarity. It can therefore be assumed they also have a dual mode of binding. In this context, the meaning of “reversed” cohesin types in B. cellulosolvens would indicate higher requirement in flexibility of scaffoldin assembly (type I interactions) rather than in enzyme integration (type II interactions).
As has been reported for other cellulosome-producing bacteria, some of the B. cellulosolvens scaffoldin genes are assembled in sca gene clusters (Fig. 2b) [16, 70]. Here, most of the clusters (with the exception of Cluster B) are only composed of scaffoldin genes without any GH genes. Cluster C consisted of genes coding for the two major multiple cohesin-carrying proteins (ScaA1 and ScaB) and represents the classic sca cluster . However, we also find additional scaffoldin clusters in the B. cellulosolvens genome: cluster E, formed from genes for ScaD, ScaX1, ScaX2, and ScaQ, also includes genes encoding two ABC (ATP-binding cassette) transporters: one in between scaX2 and scaX and another (including a C-terminal dockerin) downstream of scaX1. The ABC transporter may assist in uptake of degraded carbohydrates by the cell and the dockerin-bearing ABC transporter seems to result from a fusion of two ORFs. ScaD, ScaX1, and ScaX2 all include cell-anchoring elements (an SLH and CSBMs, respectively). Interestingly, ScaQ has an unknown protein module, which might serve a similar cell-anchoring function, considering its location on the same cluster. Moreover, the type I cohesins of the four genes are all very similar in sequence composition (Fig. 3). Only one cluster that contains genes encoding ScaI, ScaS, and ScaR3 (Fig. 2b), includes downstream genes for a GH9 and a two-component sensor, which may participate in their regulation.
The binding specificities of some of the expressed cohesins and dockerins remain unknown (data not shown). Possible reasons could be inappropriate folding of the recombinant fusion proteins or the fact that a relevant cohesin partner was not among those selected in our study. Intriguingly, none of the selected X60-linked dockerins (both scaffoldin- and enzyme-associated) bound to any of the cohesins. Since the X60 module is widely represented in the B. cellulosolvens genome, we expected to observe positive interactions for X60-linked dockerins. Indeed, a previous study reported that X60-dockerin modular dyads from other bacteria did in fact exhibit binding interactions with appropriate type II cohesins . For some cohesins and dockerins, the expression of longer scaffoldin sequences that included linkers or additional cohesin(s) improved significantly the binding capacities. This suggests, that not only the specific cohesin sequence is important for dockerin binding, but the adjacent protein sequence and modular structure could impact the stabilization of the interaction. Previous research also emphasized the importance of linker length and specific position of a module in a given scaffoldin .
Interestingly, the ScaL1 dockerin, that failed to show any binding activity for the selected cohesins, exhibited high-affinity binding with the lysate of B. cellulosolvens grown on cellobiose (Additional file 7: Figure S6). Despite the fact that we failed to demonstrate an in vitro interaction (for the reasons stated above) between the ScaL1 dockerin and any of the cohesins tested, this result confirms that the module presumably interacts with its partner in vivo.
In addition to the intraspecies interactions described in this work, we revealed inter-species interaction between the type II B. cellulosolvens cellulosomal components and those derived from two other cellulosomal bacteria: A. cellulolyticus and C. clariflavum (Additional file 8: Figure S7). These results raise the possibility of inter-species cross-reactivity, which may reflect diversification and increased cellulosomal degradation capacities for efficient carbohydrate degradation in nature. It was particularly interesting to examine the interaction between the same types of cohesins and dockerins, which are reversed in all other known bacteria in comparison to B. cellulosolvens. Thus, in the case of inter-species interaction, the cohesins from primary scaffoldin ScaA1 of B. cellulosolvens (type II) successfully bound a dockerin harbored by an adaptor scaffoldin (type II in other species). It is interesting to note that both B. cellulosolvens and A. cellulolyticus were originally isolated from sewage sludge, and C. clariflavum was first isolated from an anaerobic thermophilic methanogenic sludge. In this context, cross-species interaction among the different type II components has indeed been observed previously , and the ability of the B. cellulosolvens primary scaffoldin ScaA1 to bind strongly to the primary ScaA scaffoldins of either the thermophilic C. clariflavum or the mesophilic A. cellulolyticus via type II cohesin–dockerin interactions may indicate its complex relationship with other cellulose-degrading microbes within a specific ecological niche.
The present work has revealed the binding properties of a large number of cellulosomal elements and described a multiplicity of potential cell-free or cell-associated elaborate cellulosomal arrangements in B. cellulosolvens. These cell-free or cell-associated cellulosome complexes could be targeted to the polysaccharides substrate and include an extremely large variety of different plant cell wall-degrading enzymes and proteases via multiple scaffoldin assemblies. The accumulated knowledge of the cellulosomal components in newly discovered cellulosome-producing bacteria enables comparative evaluation of the variety of possible cellulosome architectures and/or cohesin-dockerin functions in as-yet-undescribed and/or uncharacterized cellulosome-producing bacteria. Moreover, the extensive cellulosomal system of B. cellulosolvens bears potential to provide a significant reservoir of novel components for subsequent cellulosomal research thus promoting future application of designer cellulosomes and other types of biotechnological assemblies [72, 73].
cell surface-binding module
fibronectin type III domain
enzyme-linked immunosorbent assay
splicing and auto-cleavage bacterial intein-like domain
open reading frame
anthrax protective antigen domain
polymerase and histidinol phosphatase domain
bacterial pre-peptidase C-terminal domain
repeat domain in Vibrio, Colwellia, Bradyrhizobium, and Shewanella
X-module coupled with a type II dockerin
Liao JC, Mi L, Pontrelli S, Luo S. Fuelling the future: microbial engineering for the production of sustainable biofuels. Nat Rev Microbiol. 2016;14:288–304.
Bayer EA, Lamed R, Himmel ME. The potential of cellulases and cellulosomes for cellulosic waste management. Curr Opin Biotechnol. 2007;18:237–45.
Bayer EA, Belaich JP, Shoham Y, Lamed R. The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Annu Rev Microbiol. 2004;58:521–54.
Fontes CM, Gilbert HJ. Cellulosomes: highly efficient nanomachines designed to deconstruct plant cell wall complex carbohydrates. Annu Rev Biochem. 2010;79:655–81.
Bayer EA, Morag E, Lamed R. The cellulosome-a treasure-trove for biotechnology. Trends Biotechnol. 1994;12:379–86.
Artzi L, Bayer EA, Moraïs S. Cellulosomes: bacterial nanomachines for dismantling plant polysaccharides. Nat Rev Microbiol. 2016;15:83–95.
Brown SD, Lamed R, Morag E, Borovok I, Shoham Y, Klingeman DM, et al. Draft genome sequences for Clostridium thermocellum wild-type strain YS and derived cellulose adhesion-defective mutant strain AD2. J Bacteriol. 2012;194:3290–1.
Morrison M, Daugherty SC, Nelson WC, Davidsen T, Nelson KE. The FibRumBa database: a resource for biologists with interests in gastrointestinal microbial ecology, plant biomass degradation, and anaerobic microbiology. Microb Ecol. 2010;59:212–3.
Ding SY, Lamed R, Bayer EA, Himmel ME. The bacterial scaffoldin: structure, function and potential applications in the nanosciences. Genet Eng. 2003;25:209–25.
Shimon LJ, Bayer EA, Morag E, Lamed R, Yaron S, Shoham Y, et al. A cohesin domain from Clostridium thermocellum: the crystal structure provides new insights into cellulosome assembly. Structure. 1997;5:381–90.
Leibovitz E, Béguin P. A new type of cohesin domain that specifically binds the dockerin domain of the Clostridium thermocellum cellulosome-integrating protein CipA. J Bacteriol. 1996;178:3077–84.
Ding SY, Rincon MT, Lamed R, Martin JC, McCrae SI, Aurilia V, et al. Cellulosomal scaffoldin-like proteins from Ruminococcus flavefaciens. J Bacteriol. 2001;183:1945–53.
Gerngross UT, Romaniec MP, Kobayashi T, Huskisson NS, Demain AL. Sequencing of a Clostridium thermocellum gene (cipA) encoding the cellulosomal SL-protein reveals an unusual degree of internal homology. Mol Microbiol. 1993;8:325–34.
Kakiuchi M, Isui A, Suzuki K, Fujino T, Fujino E, Kimura T, et al. Cloning and DNA sequencing of the genes encoding Clostridium josui scaffolding protein CipA and cellulase CelD and identification of their gene products as major components of the cellulosome. J Bacteriol. 1998;180:4303–8.
Ding SY, Bayer EA, Steiner D, Shoham Y, Lamed R. A novel cellulosomal scaffoldin from Acetivibrio cellulolyticus that contains a family 9 glycosyl hydrolase. J Bacteriol. 1999;181:6720–9.
Dassa B, Borovok I, Lamed R, Henrissat B, Coutinho P, Hemme CL, et al. Genome-wide analysis of Acetivibrio cellulolyticus provides a blueprint of an elaborate cellulosome system. BMC Genom. 2012;13:210.
Artzi L, Dassa B, Borovok I, Shamshoum M, Lamed R, Bayer EA. Cellulosomics of the cellulolytic thermophile Clostridium clariflavum. Biotechnol Biofuels. 2014;7:100.
Xu Q, Gao W, Ding SY, Kenig R, Shoham Y, Bayer EA, et al. The cellulosome system of Acetivibrio cellulolyticus includes a novel type of adaptor protein and a cell surface anchoring protein. J Bacteriol. 2003;185:4548–57.
Ben David Y, Dassa B, Borovok I, Lamed R, Koropatkin NM, Martens EC, et al. Ruminococcal cellulosome systems from rumen to human. Environ Microbiol. 2015;17:3407–26.
Lemaire M, Ohayon H, Gounon P, Fujino T, Béguin P. OlpB, a new outer layer protein of Clostridium thermocellum, and binding of its S-layer-like domains to components of the cell envelope. J Bacteriol. 1995;177:2451–9.
Pagès S, Bélaïch A, Bélaïch JP, Morag E, Lamed R, Shoham Y, et al. Species-specificity of the cohesin-dockerin interaction between Clostridium thermocellum and Clostridium cellulolyticum: prediction of specificity determinants of the dockerin domain. Proteins Struct Funct Genet. 1997;29:517–27.
Murray WD, Sowden LC, Colvin JR. Bacteroides cellulosolvens sp. nov., a cellulolytic species from sewage sludge. Int J Syst Bacteriol. 1984;34:185–7.
Lamed R, Morag E, Mor-Yosef O, Bayer EA. Cellulosome-like entities in Bacteroides cellulosolvens. Curr Microbiol. 1991;22:27–33.
Lin C, Urbance JW, Stahl DA. Acetivibrio cellulolyticus and Bacteroides cellulosolvens are members of the greater clostridial assemblage. FEMS Microbiol Lett. 1994;124:151–5.
Horino H, Fujita T, Tonouchi A, et al. Description of Anaerobacterium chartisolvens gen. nov., sp. nov., an obligately anaerobic bacterium from Clostridium rRNA cluster III isolated from soil of a Japanese rice field, and reclassification of Bacteroides cellulosolvens Murray et al. 1984 as Pse. Int J Syst Evol Microbiol. 1984;2014(64):1296–303.
Haimovitz R, Barak Y, Morag E, Voronov-Goldman M, Shoham Y, Lamed R, et al. Cohesin-dockerin microarray: diverse specificities between two complementary families of interacting protein modules. Proteomics. 2008;8:968–79.
Murray WD. Increased cellulose hydrolysis by Bacteroides cellulosolvens in a simplified synthetic medium. J Biotechnol. 1985;3:131–40.
Ding SY, Bayer EA, Steiner D, Shoham Y, Lamed R. A scaffoldin of the Bacteroides cellulosolvens cellulosome that contains 11 type II cohesins. J Bacteriol. 2000;182:4915–25.
Xu Q, Bayer EA, Goldman M, Kenig R, Shoham Y, Lamed R. Architecture of the Bacteroides cellulosolvens cellulosome: description of a cell surface-anchoring scaffoldin and a family 48 cellulase. J Bacteriol. 2004;186:968–77.
Dassa B, Utturkar S, Hurt RA, Klingeman DM, Keller M, Xu J, et al. Near-complete genome sequence of the cellulolytic bacterium Bacteroides (Pseudobacteroides) cellulosolvens ATCC 35603. Genome Announc. 2015;3:e01022.
Slutzki M, Ruimy V, Morag E, Barak Y, Haimovitz R, Lamed R, et al. High-throughput screening of cohesin mutant libraries on cellulose microarrays. Methods Enzymol. 2012;510:453–63.
Artzi L, Morag E, Barak Y, Lamed R, Bayer EA. Clostridium clariflavum: key cellulosome players are revealed by proteomic analysis. MBio. 2015;6:e00411–5.
Rincon MT, Dassa B, Flint HJ, Travis AJ, Jindou S, Borovok I, et al. Abundance and diversity of dockerin-containing proteins in the fiber-degrading rumen bacterium, Ruminococcus flavefaciens FD-1. PLoS ONE. 2010;5:e12476.
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42:D490–5.
Barak Y, Handelsman T, Nakar D, Mechaly A, Lamed R, Shoham Y, et al. Matching fusion protein systems for affinity analysis of two interacting families of proteins: the cohesin-dockerin interaction. J Mol Recognit. 2005;18:491–501.
Morag E, Lapidot A, Govorko D, Lamed R, Wilchek M, Bayer EA, et al. Expression, purification, and characterization of the cellulose-binding domain of the scaffoldin subunit from the cellulosome of Clostridium thermocellum. Appl Environ Microbiol. 1995;61:1980–6.
Blouzard JC, Coutinho PM, Fierobe HP, Henrissat B, Lignon S, Tardif C, et al. Modulation of cellulosome composition in Clostridium cellulolyticum: adaptation to the polysaccharide environment revealed by proteomic and carbohydrate-active enzyme analyses. Proteomics. 2010;10:541–54.
Fendri I, Tardif C, Fierobe HP, Lignon S, Valette O, Pagès S, et al. The cellulosomes from Clostridium cellulolyticum: identification of new components and synergies between complexes. FEBS J. 2009;276:3076–86.
Berg Miller ME, Antonopoulos DA, Rincon MT, Band M, Bari A, Akraiko T, et al. Diversity and strain specificity of plant cell wall degrading enzymes revealed by the draft genome of Ruminococcus flavefaciens FD-1. PLoS ONE. 2009;4:e6650.
Dassa B, Borovok I, Ruimy-Israeli V, Lamed R, Flint HJ, Duncan SH, et al. Rumen cellulosomics: divergent fiber-degrading strategies revealed by comparative genome-wide analysis of six ruminococcal strains. PLoS ONE. 2014;9:e99221.
Lamed R, Setter E, Bayer EA. Characterization of a cellulose-binding, cellulase-containing complex in Clostridium thermocellum. J Bacteriol. 1983;156:828–36.
Wu JHD, Orme-Johnson WH, Demain AL. Two components of an extracellular protein aggregate of Clostridium thermocellum together degrade crystalline cellulose. Biochemistry. 1988;27:1703–9.
Giuliano C, Khan AW. Cellulase and sugar formation by Bacteroides cellulosolvens, a newly isolated cellulolytic anaerobe. Appl Environ Microbiol. 1984;48:446–8.
Giuliano C, Khan AW. Conversion of cellulose to sugars by resting cells of a mesophilic anaerobe, Bacteriodes cellulosolvens. Biotechnol Bioeng. 1985;27:980–3.
Hu Y, Hao X, Wang J, Cao Y. Enhancing anaerobic digestion of lignocellulosic materials in excess sludge by bioaugmentation and pre-treatment. Waste Manag. 2016;49:55–63.
Bayer EA, Shoham Y, Lamed R. Lignocellulose-decomposing bacteria and their enzyme systems. The Prokaryotes. 2013:215–66.
Jindou S, Xu Q, Kenig R, Shulman M, Shoham Y, Bayer EA, et al. Novel architecture of family-9 glycoside hydrolases identified in cellulosomal enzymes of Acetivibrio cellulolyticus and Clostridium thermocellum. FEMS Microbiol Lett. 2006;254:308–16.
Ravachol J, Borne R, Tardif C, de Philip P, Fierobe HP. Characterization of all family-9 glycoside hydrolases synthesized by the cellulosome-producing bacterium Clostridium cellulolyticum. J Biol Chem. 2014;289:7335–48.
Bayer EA, Shimon LJ, Shoham Y, Lamed R. Cellulosomes-structure and ultrastructure. J Struct Biol. 1998;124:221–34.
Bensoussan L, Moraïs S, Dassa B, Friedman N, Henrissat B, Lombard V, et al. Broad phylogeny and functionality of cellulosomal components in the bovine rumen microbiome. Environ Microbiol. 2016;19:185–97.
Levy-Assaraf M, Voronov-Goldman M, Rozman Grinberg I, Weiserman G, Shimon LJ, Jindou S, et al. Crystal structure of an uncommon cellulosome-related protein module from Ruminococcus flavefaciens that resembles papain-like cysteine peptidases. PLoS ONE. 2013;8:e56138.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Mechaly A, Yaron S, Lamed R, Fierobe HP, Belaich A, Belaich JP, et al. Cohesin-dockerin recognition in cellulosome assembly: experiment versus hypothesis. Proteins Struct Funct Genet. 2000;39:170–7.
Zhao G, Ali E, Sakka M, Kimura T, Sakka K. Binding of S-layer homology modules from Clostridium thermocellum SdbA to peptidoglycans. Appl Microbiol Biotechnol. 2006;70:464–9.
Fujino T, Béguin P, Aubert JP. Organization of a Clostridium thermocellum gene cluster encoding the cellulosomal scaffolding protein CipA and a protein possibly involved in attachment of the cellulosome to the cell surface. J Bacteriol. 1993;175:1891–9.
Lupas A, Engelhardt H, Peters J, Santarius U, Volker S, Baumeister W. Domain structure of the Acetogenium kivui surface layer revealed by electron crystallography and sequence analysis. J Bacteriol. 1994;176:1224–33.
Adams JJ, Jang CJ, Spencer HL, Elliott M, Smith SP. Expression, purification and structural characterization of the scaffoldin hydrophilic X-module from the cellulosome of Clostridium thermocellum. Protein Expr Purif. 2004;38:258–63.
Slutzki M, Jobby MK, Chitayat S, Karpol A, Dassa B, Barak Y, et al. Intramolecular clasp of the cellulosomal Ruminococcus flavefaciens ScaA dockerin module confers structural stability. FEBS Open Bio. 2013;3:398–405.
Mechaly A, Fierobe HP, Belaich A, Belaich JP, Lamed R, Shoham Y, et al. Cohesin-dockerin interaction in cellulosome assembly: a single hydroxyl group of a dockerin domain distinguishes between nonrecognition and high affinity recognition. J Biol Chem. 2001;276:9883–8.
Adams JJ, Webb BA, Spencer HL, Smith SP. Structural characterization of type II dockerin module from the cellulosome of Clostridium thermocellum: calcium-induced effects on conformation and target recognition. Biochemistry. 2005;44:2173–82.
Salamitou S, Raynaud O, Lemaire M, Coughlan M, Béguin P, Aubert JP. Recognition specificity of the duplicated segments present in Clostridium thermocellum endoglucanase CelD and in the cellulosome-integrating protein CipA. J Bacteriol. 1994;176:2822–7.
Xu Q, Resch MG, Podkaminer K, Yang S, Baker JO, Donohoe BS, et al. Dramatic performance of Clostridium thermocellum explained by its wide range of cellulase modalities. Sci. Adv. 2016;2:e1501254.
Pagès S, Gal L, Bélaïch A, Gaudin C, Tardif C, Bélaïch JP. Role of scaffolding protein CipC of Clostridium cellulolyticum in cellulose degradation. J Bacteriol. 1997;179:2810–6.
Pinheiro BA, Gilbert HJ, Sakka K, Sakka K, Fernandes VO, Prates JAM, et al. Functional insights into the role of novel type I cohesin and dockerin domains from Clostridium thermocellum. Biochem J. 2009;424:375–84.
Voronov-Goldman M, Yaniv O, Gul O, Yoffe H, Salama-Alber O, Slutzki M, et al. Standalone cohesin as a molecular shuttle in cellulosome assembly. FEBS Lett. 2015;589:1569–76.
Moraïs S, Ben David Y, Bensoussan L, Duncan SH, Koropatkin NM, Martens EC, et al. Enzymatic profiling of cellulosomal enzymes from the human gut bacterium, Ruminococcus champanellensis, reveals a fine-tuned system for cohesin-dockerin recognition. Environ Microbiol. 2016;18:542–56.
Carvalho AL, Dias FMV, Nagy T, Prates JAM, Proctor MR, Smith N, et al. Evidence for a dual binding mode of dockerin modules to cohesins. Proc Natl Acad Sci USA. 2007;104:3089–94.
Pinheiro BA, Proctor MR, Martinez-Fleites C, Prates JA, Money VA, Davies GJ, et al. The Clostridium cellulolyticum dockerin displays a dual binding mode for its cohesin partner. J Biol Chem. 2008;283:18422–30.
Cameron K, Weinstein JY, Zhivin O, Bule P, Fleishman SJ, Alves VD, et al. Combined crystal structure of a type-I cohesin, mutation and affinity-binding studies reveal structural determinants of cohesin-dockerin specificity. J Biol Chem. 2015;290:16215–25.
Jindou S, Brulc JM, Levy-Assaraf M, Rincon MT, Flint HJ, Berg ME, et al. Cellulosome gene cluster analysis for gauging the diversity of the ruminal cellulolytic bacterium Ruminococcus flavefaciens. FEMS Microbiol Lett. 2008;285:188–94.
Jindou S, Borovok I, Rincon MT, Flint HJ, Antonopoulos DA, Berg ME, et al. Conservation and divergence in cellulosome architecture between two strains of Ruminococcus flavefaciens. J Bacteriol. 2006;188:7971–6.
Vazana Y, Barak Y, Unger T, Peleg Y, Shamshoum M, Ben-Yehezkel T, et al. A synthetic biology approach for evaluating the functional contribution of designer cellulosome components to deconstruction of cellulosic substrates. Biotechnol Biofuels. 2013;6:182.
Vazana Y, Moraïs S, Barak Y, Lamed R, Bayer EA. Designer cellulosomes for enhanced hydrolysis of cellulosic substrates. Methods Enzymol. 2012;510:429–52.
OZ and EAB designed the research. OZ performed the experiments. OZ, SM, RL, and EAB analyzed the results. BD and BH analyzed the genome data. SMU and SDB performed the genome sequencing. OZ, BD, SM, and EAB wrote the manuscript. All authors read and approved the final manuscript.
The authors appreciate the assistance of Dr. Yoav Barak for his expert professional contributions to this work.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
This research was supported by Grant No. 1349 from the Israel Science Foundation (ISF) and by the United States-Israel Binational Science Foundation (BSF), Jerusalem, Israel. In addition, a joint Research Grant (No. 2566/16) from the Israel Science Foundation (ISF)-National Natural Science Foundation of China (NSFC) entitled “Decoding the Pseudobacteroides cellulosolvens biomass-sensing mega-regulon” is gratefully acknowledged. The authors appreciate the support of the European Union, Area NMP.2013.1.1-2: Self-assembly of naturally occurring nanosystems: CellulosomePlus Project number: 604530 and European Union Horizon 2020 contract: Sustainable production of next generation biofuels from waste streams: Waste2Fuels. This research was funded in part by the Bioenergy Science Center (BESC), which is a U.S. Department of Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research in the DOE Office of Science. ORNL is managed by UT-Battelle, LLC, Oak Ridge, TN, USA, for the DOE under contract DE-AC05-00OR22725. EAB is the incumbent of The Maynard I. and Elaine Wishner Chair of Bio-organic Chemistry.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. List of primers for the Bacteroides cellulosolvens cohesin and dockerin modules that were cloned in this study. Restriction enzyme sites are shown in bold.
Additional file 2: Figure S1. Hydrolysis of beechwood xylan by cellulosome fractions of Bacteroides cellulosolvens. The two cellulosomal complexes (high-molecular-weight complex and low-molecular-weight complex) isolated from two different growth media (CB and MCC) were tested for their catalytic activity on beechwood xylan in order to demonstrate its ability to degrade it. The Clostridium thermocellum cellulosome (kindly provided by CelDezyner LTD, Rehovot, Israel) was also tested for catalytic activity as a positive control of the catalytic activity of the B. cellulosolvens cellulosomes. CB, cellobiose; MCC, microcrystalline cellulose; Ct, Clostridium thermocellum.
Additional file 3: Figure S2. Multiple sequence alignment of 87 cohesin sequences, originating from the genomes of Bacteroides cellulosolvens, Acetivibrio cellulolyticus (Ac), Clostridium thermocellum (Ct), Ruminococcus flavefaciens (Rf) and Ruminococcus champanellensis (Rc). Alignment length: 175; Strongly similar (:): 1 residue =0.57%; Weakly similar (.): 0 residue =0.57%.
Additional file 4: Figure S3. Multiple sequence alignment of the five Bacteroides cellulosolvens type I and miscellaneous dockerin modules. The alignment shows two internal dockerin repeats of B. cellulosolvens type I and miscellaneous dockerins that contain unique sequences. The left part of the sequence represents duplicated sequence 1, and the right sequence part represents duplicated sequence 2. Cyan highlight indicates putative calcium-binding residues. Yellow highlight indicates putative recognition residues. A: Alignment length: 71. Identity (*): 30 residues = 42.3 %. Strongly similar (:): 11 residues = 15.5 %. Weakly similar (.): 9 residues = 12.7 %. C: Alignment length: 65. Identity (*): 28 residues = 43.1 %. Strongly similar (:): 17 residues = 26.2 %. Weakly similar (.): 7 residues = 10.8 %. D: Alignment length: 78. Identity (*): 7 identical residue = 9 %. Strongly similar (:): 13 residues = 16.7 %. Weakly similar (.): 5 residues = 6.4 %. E: Fragmented dockerins: Two ORFs that resemble a dockerin sequence were found in the genome. One of the two ORFs was not annotated and presented with its ordinal number – ORF1413.
Additional file 5: Figure S4. Multiple sequence alignment of the 146 Bacteroides cellulosolvens type II dockerin modules. The alignment shows two internal dockerin repeats of B. cellulosolvens and was used to create Figure 6 representing the Weblogo of the dockerin repeats. The left part of the sequence (before the hyphen) represents duplicated sequence 1 and the right part (after the hyphen) represents duplicated sequence 2. Cyan highlight indicates putative calcium-binding residues. Yellow highlight indicates putative recognition residues. Alignment length: 65. Identity (*): 3 residues = 4.6 %. Strongly similar (:): 3 residues = 4.6 %. Weakly similar (.): 4 residues = 6.2 %.
Additional file 6: Figure S5. Multiple sequence alignment of the 49 miscellaneous Bacteroides cellulosolvens type II dockerin modules. The alignment shows two internal dockerin repeats of the B. cellulosolvens type II dockerins that contain unique sequences (particularly in the first calcium-binding loop) but remain type II dockerins. The left part of the sequence (before the internal hyphens) represents duplicated sequence 1 and the right part (after the internal hyphens) represents duplicated sequence 2. Cyan highlight indicates putative calcium-binding residues. Yellow highlight indicates putative recognition residues. Alignment length: 77. Identity (*): 5 residues = 6.5%. Strongly similar (:): 9 residues = 11.7%. Weakly similar (.): 3 residues = 3.9%.
Additional file 7: Figure S6. Determination of dockerin specificity to Bacteroides cellulosolvens cell lysate components by affinity-based ELISA. The 96-well ELISA plates were coated with B. cellulosolvens cell lysate (grown on cellobiose), and various concentrations of Xyn-Docs were used to detect cohesin-dockerin interactions. Abbreviations: Doc, dockerin.
Additional file 8: Figure S7. Determination of the inter-species interactions of Bacteroides cellulosolvens cell lysate by affinity-based ELISA. The 96-well ELISA plates were coated with B. cellulosolvens cell lysate (grown on cellobiose) and various concentrations of Xyn-Docs from three different bacteria were used to detect cohesin-dockerin interactions. Abbreviations: Doc, dockerin; Cc, Clostridium clariflavum; Rf, Ruminococcus flavefaciens; Ac, Acetivibrio cellulolyticus; GH9, Glycoside hydrolases of family 9. The dockerins were chosen to include the three previously defined types: Ac-GH9-Doc is a representative of type I dockerins; Cc-ScaA-Doc and Ac-ScaA-Doc represent type II dockerins; and Rf-ScaA-Doc represents type III dockerins. Here we show that B. cellulosolvens cell lysate is capable of crossreaction with type II dockerins from the primary scaffoldin ScaA from two different bacteria but not type I or III.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zhivin, O., Dassa, B., Moraïs, S. et al. Unique organization and unprecedented diversity of the Bacteroides (Pseudobacteroides) cellulosolvens cellulosome system. Biotechnol Biofuels 10, 211 (2017). https://doi.org/10.1186/s13068-017-0898-6
- Cellulolytic bacteria
- Glycoside hydrolases