- Open Access
Discovery of new cellulases from the metagenome by a metagenomics-guided strategy
Biotechnology for Biofuels volume 9, Article number: 138 (2016)
Energy shortage has become a global problem. Production of biofuels from renewable biomass resources is an inevitable trend of sustainable development. Cellulose is the most abundant and renewable resource in nature. Lack of new cellulases with unique properties has become the bottleneck of the efficient utilization of cellulose. Environmental metagenomes are regarded as huge reservoirs for a variety of cellulases. However, new cellulases cannot be obtained easily by functional screening of metagenomic libraries.
In this work, a metagenomics-guided strategy for obtaining new cellulases from the metagenome was proposed. Metagenomic sequences of DNA extracted from the anaerobic beer lees converting consortium enriched at thermophilic conditions were assembled, and 23 glycoside hydrolase (GH) sequences affiliated with the GH family 5 were identified. Among the 23 GH sequences, three target sequences (designated as cel7482, cel3623 and cel36) showing low identity with those known GHs were chosen as the putative cellulase genes to be functionally expressed in Escherichia coli after PCR cloning. The three cellulases were classified into endo-β-1,4-glucanases by product pattern analysis. The recombinant cellulases were more active at pH 5.5 and within a temperature range of 60–70 °C. Computer-assisted 3D structure modeling indicated that the active residues in the active site of the recombinant cellulases were more similar to each other compared with non-active site residues. The recombinant cel7482 was extremely tolerant to 2 M NaCl, suggesting that cel7482 may be a halotolerant cellulase. Moreover, the recombinant cel7482 was shown to have an ability to resist three ionic liquids (ILs), which are widely used for cellulose pretreatment. Furthermore, active cel7482 was secreted by the twin-arginine translocation (Tat) pathway of Bacillus subtilis 168 into the culture medium, which facilitates the subsequent purification and reduces the formation of inclusion body in the context of overexpression.
This study demonstrated a simple and efficient method for direct cloning of new cellulase genes from environmental metagenomes. In the future, the metagenomics-guided strategy may be applied to the high-throughput screening of new cellulases from environmental metagenomes.
Cellulose, a renewable biopolymer composed of D-glucopyranose units linked by β-1,4-glucosidic bonds, is commonly used as raw material for the production of important industrial chemicals such as soluble sugars and biofuels . Cellulases involved in hydrolyzing cellulose are composed of endoglucanases, exoglucanases and β-glucosidases. Endoglucanases (EC 220.127.116.11) hydrolyze the internal bonds randomly in the cellulose chain into cellobiose or cello-oligosaccharides. Exoglucanases (EC 18.104.22.168) release cellobiose from either the reducing or the non-reducing ends of the cellulose chain. β-Glucosidases (EC 22.214.171.124) hydrolyze cellobiose to glucose .
Discovery of novel enzymes through metagenomics has recently been shown to have enormous potential for obtaining a wide variety of useful biocatalysts [2, 3]. Functional metagenomics has become a routine method for the discovery of industrially relevant enzymes from natural and artificially engineered ecosystems. However, very few active clones could be obtained from huge quantities of clones tested using functional screening protocols .
Although some cellulase genes had been obtained from the metagenome by functional screening [4–6], the efficiency of the strategy for the discovery of new cellulases from environmental metagenomes can hardly meet the increasing industrial demand . The number of new metagenome-derived cellulases could be significantly increased when novel cellulase sequences obtained by high-throughput sequencing are selected for subsequent expression.
Many researches have focused on the annotation of carbohydrate-active enzymes (CAZymes) in cellulose-degrading consortia using metagenomic analysis [8, 9]. CAZyme families play a crucial role in breakdown of complex carbohydrates. The CAZy database (http://www.cazy.org) defined 135 families of glycoside hydrolases (GHs) based on amino acid sequence similarities and structural features, which provides a powerful tool to annotate functions of obtained GH genes . GHs are a widespread category of enzymes capable of cleaving the glycosidic bond in the polysaccharide chains. These GHs possess different substrate specificities, and some GH families (e.g., GH5 and GH9) were found to have cellulase activity.
It is difficult to obtain new cellulases by functional screening of metagenomic libraries because the strategy is based on cellulase activity rather than sequence similarity. Metagenomic sequencing allows targeted selection of novel cellulase sequences for high-throughput expression. In this work, a metagenomics-guided strategy was applied to clone three new cellulase genes from environmental metagenomes. Furthermore, the recombinant cellulases were functionally expressed in Escherichia coli and their enzymatic characteristics were investigated. A pipeline of the metagenomics-guided strategy used in this study is shown in Fig. 1.
Cloning of three cellulase genes from the metagenome by a metagenomics-guided strategy
In this work, we have focused on the GH family 5 (GH5), mainly consisting of a variety of endoglucanases and exoglucanases that are involved in the hydrolysis of cellulose. In total, we obtained 23 GH sequences belonging to the GH5 family by metagenomic sequencing and ORFs annotation. The nucleotide sequences of 23 GHs are shown in Additional file 1: Figure S1. Among 23 GH sequences, we chose three target sequences (id_7482, 1026; id_3623, 1035; id_36, 1548 bp), which showed low similarity with those GH sequences deposited in GenBank. Using the NCBI BlastP search, id_7482 showed 54 % identity with a putative cellulase from uncultured bacterium in a laboratory biogas digester treating rice straws (GenBank accession no. AEV59734); id_3623 showed 51 % identity with an endoglucanase from Anaerolinea thermolimosa from sludge from a thermophilic UASB reactor (GenBank accession no. GAP08306); id_36 showed 48 % identity with a cellulase from Acetivibrio cellulolyticus (GenBank accession no.WP_010249757).
We designed three pairs of specific primers (F7482c/R7482c, F3623c/R3623c and F36c/R36c) to amplify the full-length target sequences. As a result, three specific DNA fragments with about 1.0, 1.0 and 1.5 kb were amplified by PCR from the metagenome of the anaerobic beer lees converting consortium (Fig. 2). Subsequently, the amplified fragments were sequenced, and the results showed that the nucleotide sequences of the amplified fragments were consistent with those of the assembled target sequences (id_7482, id_3623 and id_36). The nucleotide sequences of cel7482, cel3623 and cel36 have been deposited in GenBank under accession nos. KU168144, KU168145 and KU168146.
Expression and purification of three recombinant cellulases in E. coli
Expression of the recombinant cel7482, cel3623 and cel36 in E. coli BL21 (DE3) was demonstrated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) analysis. After induction with 0.5 mM isopropyl-β-d-thiogalactopyranoside (IPTG) for 4 h at 30 °C, clear bands corresponding to 41, 41 and 58 kDa were observed in whole-cell lysates of E. coli BL21 in SDS-PAGE (Fig. 3), which matched well with the molecular masses estimated from the deduced amino acid sequences of the recombinant cel7482, cel3623 and cel36.
Three recombinant cellulases were purified from 500 ml of E. coli cultures by immobilized metal ion affinity chromatography (IMAC) using Ni Sepharose 6 Fast Flow. The purified target proteins were observed at the position of 41, 41 and 58 kDa in SDS-PAGE (Fig. 3), which were identical to the molecular masses of the recombinant cel7482, cel3623 and cel36. The data on the purification of the recombinant cel7482, cel3623 and cel36 are summarized in Table 1.
Three purified His6-tagged cellulases were shown to have obvious hydrolytic activity toward carboxymethylcellulose (CMC). Cel7482 showed 1.6- and 2-fold higher specific activity than cel3623 and cel36, respectively. However, the recombinant cellulases had no activity against p-nitrophenylcellobioside (pNPC) and cellobiose.
To determine final products of hydrolysis of CMC by cel7482, cel3623 and cel36, HPAEC analysis was performed with enzyme assay mixtures. Cellobiose and cellotriose accumulated as products when each purified enzyme (50 μg/ml) was incubated with CMC for 1 h under optimal conditions, whereas larger cello-oligosaccharides were not visible (Fig. 4). Moreover, no glucose was detected. In contrast, no hydrolysis products were detected in reaction without enzyme.
Enzymatic characteristics of the recombinant cel7482, cel3623 and cel36
The recombinant cel7482 was stable over an acidic pH range of 4.5–6.5, maintaining more than 80 % of the maximum activity (Fig. 5a). The optimal temperature range for activity was 65–75 °C, within which the enzymatic activity was maintained more than 90 % of the maximum activity (Fig. 5b). The recombinant cel7482 retained about 87 % of its initial activity after 1 h of incubation at 70 °C (Fig. 6a). However, after incubation at 80 °C for 1 h, the enzymatic activity dropped to 20 % of its initial activity.
The optimal pH and temperature for activity of the recombinant cel3623 were 5.5 and 65 °C, respectively. The recombinant cel3623 retained over 80 % of the maximum activity between pH 4.5 and 6.5 (Fig. 5). After incubation at 65 °C for 1 h, the recombinant cel3623 maintained about 75 % of its initial activity. However, the enzymatic activity decreased by 88 % after incubation at 75 °C for 1 h (Fig. 6b).
The recombinant cel36 showed high activity in the pH range from 4.5–6.5. The optimal pH and temperature for the enzymatic activity were 5.5 and 60 °C, respectively (Fig. 5). After incubation at 60 °C for 1.5 h, the recombinant cel36 retained about 57 % of its initial activity. However, the recombinant cel36 was almost completely inactive after incubation at 70 °C for 1.5 h (Fig. 6c).
The recombinant cel7482 still retained 80 % of the maximum activity in the presence of 4 M NaCl. In contrast, the recombinant cel3623 only maintained 62 % of the maximum activity in the presence of 2 M NaCl. Activity of the recombinant cel36 was almost completely inhibited by the presence of 2 M NaCl (Fig. 7a). After pre-incubation in 2 M NaCl for 96 h, the recombinant cel7482 still maintained 50 % of its initial activity (Fig. 7b).
We further tested the ability of the recombinant cel7482 to resist three ionic liquids (ILs) (i.e., [Emim]Cl, [Bmim]Cl and [Amim]Cl). The enzymatic activity remained unchanged in the presence of 20 % of one of these ILs (Additional file 1: Figure S2), suggesting that these ILs have no inhibitory effects on the enzymatic activity.
3D structures of cel7482, cel3623 and cel36
Sequence analysis showed that cel7482 shares 64.2 % amino acid identity with cel3623 (Additional file 1: Figure S3). From the high sequence similarity, the two proteins can be deduced to have a high structural similarity. Since there are no actual X-ray crystallographic structures of cel7482 and cel3623, we built 3D structures of cel7482 and cel3623 based on the crystal structure of a family 5 endoglucanase (PDB: 1ceo) that shares 46.0 and 47.1 % amino acid identity with cel7482 and cel3623, respectively.
Figure 8 showed the predicted 3D structure of cel7482 (residues 7–340) with the active site surrounded by the inner β sheets and the outer α helixes. The predicted 3D structure of cel3623 was similar with cel7482. The important functional amino acids R49, H93, LNEL(141–144), H201, Y203, E280, and W313 in the active site of cel7482 as well as α helixes were completely conserved between cel7482, cel3623 and the template. Some of the amino acid substitutions between cel7482 and cel3623 occurred in the entryway of the active pocket of cel7482 and the substitutions contained bulky (G21W, F102S, L178F) and charged (E105S, T107K, T110E, E204L, P217K, Q245R, R246S) residues, which might affect the entry of substrate into active pocket and lead to the variation of basal and ligand-induced activities between cel7482 and cel3623.
The 3D structure of cel36 was generated based on the crystal structure of endo-1,4-beta-glucanase from B. subtilis 168 (PDB: 3pzt) that shares 33.8 % amino acid identity with cel36 (Additional file 1: Figure S4). Figure 8 showed the predicted 3D structure of cel36 (residues 76–370) with the active site surrounded by the inner β sheets and the outer α helixes. 3D structure of cel36 revealed high amino acid identity between cel36 and the template in the catalytic core. The important functional amino acids R129, N206, E207, W246, H268, Y270, S295, E296, and W330 in the active site of cel36 as well as α helixes were completely conserved between cel36 and the template.
Secretion of cel7482 by the twin-arginine translocation (Tat) pathway of B. subtilis 168
To secrete cel7482 into the culture medium, in this study, the twin-arginine signal peptide of YwbN (a strict Tat substrate in B. subtilis) was fused with the N-terminus of cel7482 and used to target cel7482 to the Tat pathway of B. subtilis 168. To demonstrate whether active cel7482 was secreted by the Tat pathway of B. subtilis 168 into the culture medium, extracellular cellulase activity was detected with B. subtilis 168 and its tat mutant strains. As a result, cellulase activity (0.266 U/ml) was detected in the culture medium of B. subtilis 168 expressing YwbN–cel7482 fusion protein. In contrast, cellulase activity (0.028 U/ml) was found in the culture medium of the total-tat 2 mutant strain lacking all Tat translocases. Additionally, cellulase activity (0.062 U/ml) was detected in the culture supernatant of the tatAyCy mutant strain lacking functional TatAyCy translocase. In contrast, cellulase activity (0.192 U/ml) was found in the culture supernatant of the tatAdCd mutant strain lacking functional TatAdCd translocase (Fig. 9).
Complete hydrolysis of CMC by the recombinant cel7482, cel3623 and cel36 was shown to accumulate cellobiose and cellotriose as products. The recombinant cellulases could not hydrolyze pNPC and cellobiose, indicating that they lacked exoglucanase or β-glucosidase activity. Results from the NCBI BlastP search indicated that the three cellulases had highest identity with those endoglucanases deposited in GenBank. Based on product pattern and sequence alignment, the three cellulases were affiliated with endo-β-1,4-glucanases. 3D structure modeling showed that the active residues in the active site of cel7482, cel3623 and cel36 were more similar to each other compared with non-active site residues, which suggested that the three proteins should have a similar function.
In this study, the recombinant cellulases were more active in the acidic pH range (4.5–6.5). The recombinant cellulases showed high activity at high temperatures (60–70 °C) and were thermostable, which might be due to the fact that the cellulases came from a thermophilic consortium. To date, several thermostable cellulases have previously been reported [11, 12].
It is hard to obtain halotolerant cellulases by functional screening of metagenomic libraries . Cel7482, which was obtained in this study by a metagenomics-oriented approach, showed significant resistance to high concentrations of NaCl. Cel7482 was extremely tolerant to 2 M NaCl and was still active in the presence of 5 M NaCl. Cel7482 may be a halotolerant cellulase due to its tolerance to high salinity environments.
Lignocellulose needs to be pretreated prior to enzymatic hydrolysis due to its inherent recalcitrance towards degradation. Although ILs have been shown to be very effective solvents for lignocellulose pretreatment, ILs strongly inhibit cellulase activity . Activity of cellulases is especially influenced by the presence of chloride ions in ILs. Therefore, halophilic cellulases are regarded as promising candidates for screening IL-tolerant cellulases. In this study, a halotolerant cellulase (cel7482) was found to be resistant to three ILs (i.e., [Emim]Cl, [Bmim]Cl and [Amim]Cl), which are widely used for cellulose dissolution. Therefore, cel7482 with ILs resistance has potential for use in a one-pot process (i.e., coupling of ILs pretreatment with enzymatic hydrolysis) in which enzymatic hydrolysis is carried out in aqueous solutions of cellulose-dissolving ILs.
The Tat pathway is utilized by bacteria for the transport of folded proteins across the cytoplasmic membrane. Because B. subtilis lacks the outer membrane, proteins that are exported to the periplasm by the Tat pathway can be secreted directly into the culture medium. In B. subtilis 168, two active Tat translocases with different substrate preferences have previously been identified as TatAdCd and TatAyCy . In this study, when fused to the YwbN signal peptide, cel7482 was secreted by the Tat pathway of B. subtilis 168 into the culture medium as judged from the ratio of extracellular activity (0.266 U/ml) to cell lysate activity (0.12 U/ml). In contrast, secretion of cel7482 was almost completely blocked in a total-tat 2 mutant, indicating that the extracellular secretion of cel7482 was strictly dependent on the Tat pathway in B. subtilis 168. In addition, no obvious secretion of cel7482 was observed in a tatAyCy mutant. However, the effective transport of cel7482 still occurred in the tatAdCd mutant. We assumed that the TatAyCy translocase may play key roles in the export of cel7482 to the extracellular milieu. In the future, the efficiency of Tat-dependent translocation of cel7482 needs to be further improved for industrial production.
In function-based screening of cellulases from environmental metagenomes, very few active clones could be obtained from metagenomic libraries. Activity of specific cellulases may not be detected by functional screening of libraries, resulting from incorrect protein folding in E. coli. More importantly, the strategy cannot select new cellulase sequences prior to cloning.
Metagenomic sequencing has some advantages compared to functional screening. Diverse GH genes have been discovered by metagenomic sequencing in termite guts , cow rumens  and biogas digesters , and they could be utilized as a resource for screening new cellulases with industrial value. In this study, a metagenomics-guided strategy for biomining new cellulases was proposed. Cloning of cellulase genes can be performed by a simple PCR protocol without the need for construction of metagenomic libraries. New cellulase sequences obtained by metagenomic sequencing were chosen for functional expression, which significantly improved the hit rates for new cellulases.
In this work, search for conserved domains of GHs by hidden Markov models (HMMs)  significantly improved the accuracy of GH annotation. In total, we obtained 23 GH sequences belonging to the GH5 family by metagenomic sequencing of the thermophilic anaerobic enrichment. Among 23 GH sequences, we selected three novel GH sequences (id_7482, id_3623 and id_36) as the putative cellulases. Based on the assembled GH sequences, the three putative cellulase genes were obtained by PCR directly from the metagenome of the thermophilic anaerobic enrichment. Furthermore, functionality of the putative cellulases was verified by expression and activity assays. In the future, the rapid detection of activity of a large number of candidate cellulases can be performed using high-throughput screening systems [4, 19]. Compared to function-based screening, the metagenomics-guided strategy is particularly applicable to the high-throughput screening of new cellulases from environmental metagenomes.
New cellulases with unique properties have potential for use in industrial processes. In this work, we propose a metagenomics-guided strategy to rapidly acquire novel cellulase sequences from the metagenome. When combined with high-throughput expression, the efficiency of this strategy for obtaining new cellulases may meet the increasing demand from industrial community. Cel7482 obtained in this study, which has superior enzymatic characteristics, may be a promising candidate for degradation of cellulosic biomass under harsh conditions. Secretion of active cel7482 into the culture medium simplifies the purification procedure and improves the stability of enzyme, which should lay a foundation for large-scale production of the enzyme.
Metagenomic sequencing and GH annotation
Anaerobic digestion sludge (ADS) collected from a local wastewater treatment plant (Shek Wu Hui Sewage Treatment Works, Hong Kong) was applied as seed sludge in this experiment. The enrichment of the beer lees fermenting consortium was carried out in a sequential batch mode in serum bottles with a working volume of 200 ml. Beer lees were applied as the primary substrate, and the temperature was controlled at 55 °C.
After 45, 75 and 120 days of enrichment, genomic DNA was extracted from 4 ml of sludge slurry from the thermophilic consortium with a FastDNA SPIN Kit for Soil (MP Biomedicals). Library size of 300 bp and reading length of 125 bp were applied for illumina high-throughput sequencing of the extracted DNA samples. The sequencing depth for the metagenomic library was 3.0 Gb.
Quality control of metagenomic raw reads derived from the illumina Hiseq 2000 platform was performed as described previously . The trimmed reads were firstly assembled using MetaVelvet (version 1.1.01) [20, 21] with kmer length of 51. The assembled contigs longer than 1000 bp  were chosen for ORFs prediction using MetaGeneMark (version 2.8)  with default parameters. Next, the amino acid sequences of the predicted ORFs were screened against HMMs collected at dbCAN  using hmmscan  with E-value cut-off of 1E−4  for particular GH families classified by the CAZy database .
Cloning of metagenome-derived cellulase genes and their expression in E. coli
Three putative cellulase genes, designated as cel7482, cel3623 and cel36, were amplified by PCR from the metagenome of the anaerobic beer lees converting consortium using three pairs of specific primers (F7482c/R7482c, F3623c/R3623c and F36c/R36c). The 50-μl PCR mixture contained 2 μl of each of forward and reverse primer, 25 μl of 2 × Ex Taq PCR MasterMix (Takara), 1 μl of template DNA, and 20 μl of ddH2O. The PCR protocols were set as below: an initial denaturation at 95 °C for 5 min, followed by 30 cycles of 94 °C for 45 s, annealing at 62 °C for 45 s, and elongation at 72 °C for 90 s, with a final extension at 72 °C for 8 min. The PCR products were run on 0.7 % agarose gel, and DNA bands with the correct size were recovered using a DNA gel purification kit (Tiangen). The purified PCR products were cloned into the pMD19-T simple vector using a Takara TA cloning kit and then the ligation products were transformed into E. coli DH5α. The positive recombinants were identified by colony PCR and sequenced.
Three cellulase genes were re-amplified by PCR from the recombinant pMD19-T vectors using three pairs of primers (F7482e/R7482e, F3623e/R3623e and F36e/R36e), digested with NdeI and XhoI and subcloned into the similarly digested expression vector pET30a (Novagen). The recombinant plasmids, designated as pET-c7482, pET-c3623 and pET-c36, were transformed into E. coli BL21 (DE3). Strains, plasmids and primers used in this study are listed in Table 2.
To achieve high-level expression of cellulase genes in E. coli BL21 (DE3), cells were grown to an OD600 of 0.6 and then induced with 0.5 mM IPTG for 4 h at 30 °C. To check the expression of target proteins, total proteins from IPTG-induced E. coli cells were subjected to 12 % SDS-PAGE analysis .
Purification of three recombinant His6-tagged cellulases
For purification of the recombinant cellulases from E. coli BL21 (DE3), 10 ml preculture was inoculated into 500 ml of LB medium and grown to an OD600 of 0.6 at 37 °C. Subsequently, 0.5 mM IPTG was added to induce recombinant protein expression. After induction at 30 °C for 4 h, cells were harvested by centrifugation at 4 °C and 6000 rpm for 5 min. Cells were resuspended in 20 ml of 10 mM Tris–HCl (pH 8.0) and disrupted by sonication on ice (500 W for 25 min with cycles of sonication of 10 s each and 15 s pause). The crude cell extracts were centrifuged at 4 °C and 10,000 rpm for 15 min to remove cell debris and unbroken cells.
The recombinant His6-tagged cellulases were purified using the IMAC according to the standard procedure  with minor modifications. In brief, the cell-free extracts were loaded onto a Ni Sepharose 6 Fast Flow (GE Healthcare), which had been equilibrated with 100 ml of binding buffer (10 mM Tris–HCl, 0.5 M NaCl, pH 8.0). After washing with 100 ml of washing buffer (10 mM Tris–HCl, 0.5 M NaCl, pH 8.0), His6-tagged cellulases were eluted from the column with elution buffer (10 mM Tris–HCl, 0.5 M NaCl, 300 mM imidazole, pH 8.0). The purity of the purified cellulases was examined by 12 % SDS-PAGE. Proteins were quantified using the Bradford method  with bovine serum albumin (BSA) as standard.
Cellulase activity assays
Enzymatic activity was assayed by measuring the amount of reducing sugar released from CMC using the 3,5-dinitrosalicylic acid (DNS) method . The assay was performed in 50 mM citrate–phosphate buffer (pH 5.5) containing 5 μg/ml purified enzyme and 1 % CMC. After incubation for 30 min at optimal temperature, 200 μl of DNS was added to stop the reaction, followed by boiling for 5 min in water. The absorbance at 540 nm (A540) was measured using a microplate reader (Thermo Scientific). Specific activity is expressed as units (1 μmol of reducing sugars released per minute) per milligram protein.
The effect of pH on activity of the recombinant cellulases was studied by incubating the purified enzyme in 50 mM citrate–phosphate buffer (pH 4.0–8.0) with CMC as substrate at optimal temperature. The effect of temperature on activity of the recombinant cellulases was also investigated by incubating the purified enzyme in 50 mM citrate–phosphate buffer (pH 5.5) with CMC as substrate at different temperatures from 50 to 80 °C.
To evaluate thermostability of the recombinant cellulases, each purified enzyme was incubated at different temperatures for different periods of time, and the residual activity was determined at optimal temperature in 50 mM citrate–phosphate buffer (pH 5.5) with CMC as substrate.
To investigate the effect of NaCl on activity of the recombinant cellulases, each purified enzyme was incubated with CMC at optimal temperature for 30 min in 50 mM citrate–phosphate buffer (pH 5.5) in the presence of 0.5–5 M NaCl. Halotolerance of the recombinant cel7482 was evaluated by incubating the enzyme in 0.5, 2 and 5 M NaCl for different periods of time, and the residual activity was determined at 70 °C in 50 mM citrate–phosphate buffer (pH 5.5) with CMC as substrate.
Resistance of the recombinant cel7482 to ILs was evaluated. The purified enzyme was incubated with CMC at 37 °C for 30 min in 50 mM citrate–phosphate buffer (pH 7.0) in the presence of 20 % of [Emim]Cl, [Bmim]Cl or [Amim]Cl.
Final products of hydrolysis of CMC by the recombinant cellulases were determined by high performance anion exchange chromatography (HPAEC) on a Dionex ICS5000 system equipped with a pulsed amperometric detector and a CarboPac PA200 column (Dionex). The column was equilibrated with 100 mM NaOH and elution was performed at a column temperature of 30 °C using a linear gradient of 0.02–0.5 M NaOH at a flow rate of 0.45 ml/min in 25 min. The enzyme assay mixtures (10 μl) were withdrawn after 1 h of incubation and used for HPAEC analysis.
The peptide sequences of cel7482 and cel3623 were submitted to the automated comparative protein modeling server SwissModel [28–30] (http://swissmodel.expasy.org/) to build up the 3D structures using the crystal structure of a family 5 endoglucanase (PDB: 1ceo)  as modeling template. The 3D structure of cel36 was generated using the crystal structure of endo-1,4-β-glucanase from Bacillus subtilis 168 (PDB: 3pzt)  as modeling template. CLUSTALW software was used to align the peptide sequences (http://www.ebi.ac.uk/Tools/msa/clustalw2/). 3D figures were created with PyMOL (http://www.pymol.org/).
Construction of recombinant B. subtilis 168 for secretory expression of cel7482
The nucleotide sequence encoding the twin-arginine signal peptide of B. subtilis YwbN  and codon-optimized cel7482 was chemically synthesized by BGI Inc., Beijing, China. The synthetic sequence was digested with KpnI and SphI and subcloned into similarly digested pWH1520, an E. coli-B. subtilis shuttle vector, to create the secretory expression vector, pWYC7. Transformation of plasmid into B. subtilis 168 was carried out using the high-osmolarity electroporation method . When B. subtilis cultures reached an OD600 of 0.6, 0.5 % xylose was added to induce expression of YwbN-cel7482 fusion protein. After induction for 24 h at 37 °C, cells and culture supernatant were separated by centrifugation and used for cellulase activity assays.
carbohydrate active enzymes
sodium dodecyl sulfate–polyacrylamide gel electrophoresis
immobilized metal ion affinity chromatography
hidden Markov models
anaerobic digestion sludge
bovine serum albumin
high performance anion exchange chromatography
Lynd LR, Weimer PJ, Van Zyl WH, Pretorius IS. Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev. 2002;66:506–77.
Fernández-Arrojo L, Guazzaroni ME, López-Cortés N, Beloqui A, Ferrer M. Metagenomic era for biocatalyst identification. Curr Opin Biotechnol. 2010;21:725–33.
Ferrer M, Martínez-Martínez M, Bargiela R, Streit WR, Golyshina OV, Golyshin PN. Estimating the success of enzyme bioprospecting through metagenomics: current status and future trends. Microb Biotechnol. 2016;9:22–34.
Ko KC, Lee JH, Han Y, Choi JH, Song JJ. A novel multifunctional cellulolytic enzyme screened from metagenomic resources representing ruminal bacteria. Biochem Biophys Res Commun. 2013;441:567–72.
Mewis K, Armstrong Z, Song YC, Baldwin SA, Withers SG, Hallam SJ. Biomining active cellulases from a mining bioremediation system. J Biotechnol. 2013;167:462–71.
Yan X, Geng A, Zhang J, Wei Y, Zhang L, Qian C, Wang Q, Wang S, Zhou Z. Discovery of (hemi-) cellulase genes in a metagenomic library from a biogas digester using 454 pyrosequencing. Appl Microbiol Biotechnol. 2013;97:8173–82.
Uchiyama T, Miyazaki K. Functional metagenomics for enzyme discovery: challenges to efficient screening. Curr Opin Biotechnol. 2009;20:616–22.
Xia Y, Ju F, Fang HH, Zhang T. Mining of novel thermo-stable cellulolytic genes from a thermophilic cellulose-degrading consortium by metagenomics. PLoS One. 2013;8:e53779.
Xia Y, Wang Y, Fang HH, Jin T, Zhong H, Zhang T. Thermophilic microbial cellulose decomposition and methanogenesis pathways recharacterized by metatranscriptomic and metagenomic analysis. Sci Rep. 2014;4:6708.
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:D233–8.
Leis B, Heinze S, Angelov A, Pham VT, Thürmer A, Jebbar M, Golyshin PN, Streit WR, Daniel R, Liebl W. Functional screening of hydrolytic activities reveals an extremely thermostable cellulase from a deep-sea archaeon. Front Bioeng Biotechnol. 2015;3:95.
Peng X, Qiao W, Mi S, Jia X, Su H, Han Y. Characterization of hemicellulase and cellulase from the extremely thermophilic bacterium Caldicellulosiruptor owensensis and their potential application for bioconversion of lignocellulosic biomass without pretreatment. Biotechnol Biofuels. 2015;8:131.
Voget S, Steele HL, Streit WR. Characterization of a metagenome-derived halotolerant cellulase. J Biotechnol. 2006;126:26–36.
Wahlström RM, Suurnäkki A. Enzymatic hydrolysis of lignocellulosic polysaccharides in the presence of ionic liquids. Green Chem. 2015;17:694–714.
Jongbloed JD, Grieger U, Antelmann H, Hecker M, Nijland R, Bron S, van Dijl JM. Two minimal Tat translocases in Bacillus. Mol Microbiol. 2004;54:1319–25.
Warnecke F, Luginbühl P, Ivanova N, Ghassemian M, Richardson TH, Stege JT, et al. Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite. Nature. 2007;450:560–5.
Hess M, Sczyrba A, Egan R, Kim TW, Chokhawala H, Schroth G, et al. Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science. 2011;331:463–7.
Eddy SR. Accelerated Profile HMM Searches. PLoS Comput Biol. 2011;7:e1002195.
Ko KC, Han Y, Cheong DE, Choi JH, Song JJ. Strategy for screening metagenomic resources for exocellulase activity using a robotic, high-throughput screening system. J Microbiol Methods. 2013;94:311–6.
Namiki T, Hachiya T, Tanaka H, Sakakibara Y. MetaVelvet: an extension of Velvet assembler to de novo metagenome assembly from short sequence reads. Nucleic Acids Res. 2012;40:e155.
Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.
Albertsen M, Hansen LBS, Saunders AM, Nielsen PH, Nielsen KL. A metagenome of a full-scale microbial community carrying out enhanced biological phosphorus removal. ISME J. 2011;6:1094–106.
Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010;38:e132.
Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40:W445–51.
Sambrook J, Russell DW. Molecular cloning: a laboratory manual. 3rd ed. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 2001.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.
Miller GL. Use of dinitrosalicylic acid reagent for determination of reducing sugars. Anal Chem. 1959;31:426–8.
Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22:195–201.
Benkert P, Biasini M, Schwede T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011;27:343–50.
Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, Kiefer F, Cassarino TG, Bertoni M, Bordoli L, Schwede T. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42:W252–8.
Domínguez R, Souchon H, Lascombe M, Alzari PM. The crystal structure of a family 5 endoglucanase mutant in complexed and uncomplexed forms reveals an induced fit activation mechanism. J Mol Biol. 1996;257:1042–51.
Santos CR, Paiva JH, Sforça ML, Neves JL, Navarro RZ, Cota J, et al. Dissecting structure-function-stability relationships of a thermostable GH5-CBM3 cellulase from Bacillus subtilis 168. Biochem J. 2012;441:95–104.
Kolkman MA, van der Ploeg R, Bertels M, van Dijk M, van der Laan J, van Dijl JM, Ferrari E. The twin-arginine signal peptide of Bacillus subtilis YwbN can direct either Tat- or Sec-dependent secretion of different cargo proteins: secretion of active subtilisin via the B. subtilis Tat pathway. Appl Environ Microbiol. 2008;74:7507–13.
Xue GP, Johnson JS, Dalrymple BP. High osmolarity improves the electro-transformation efficiency of the gram-positive bacteria Bacillus subtilis and Bacillus licheniformis. J Microbiol Methods. 1999;34:183–91.
CY and TZ designed the research. CY, YX, HQ and RHL performed the research. CY, YX, HQ, ADL, YBW and TZ analyzed the data. CY, YX and TZ wrote the paper. All the authors read and approved the final manuscript.
The authors would like to thank Hong Kong GRF (HKU 172057/15E). Dr. Chao Yang would like to thank the Hong Kong Scholars Program. Dr. Yu Xia would like to thank the University of Hong Kong for the Postdoctoral Fellowship.
The authors declare that they have no competing interests.
Availability of supporting data
Supporting data could be found in Additional file 1.
Consent for publication
All the authors consented on the publication of this work.
Hong Kong GRF (HKU 172057/15E).
Chao Yang and Yu Xia contributed equally to this work
The nucleotide sequences of 23 glycoside hydrolases. Three target sequences selected in this study are colored in red. Figure S2. Resistance of the recombinant cel7482 to high concentrations of ILs. The recombinant cel7482 was incubated with CMC at 37 °C for 30 min in 50 mM citrate–phosphate buffer (pH 7.0) supplemented with 20 % of [Emim]Cl, [Bmim]Cl or [Amim]Cl. The activity in reaction without ILs was set as 100 %. Figure S3. Alignment of cel7482 and cel3623 proteins. The amino acid sequences of cel7482 and cel3623 were aligned with ClustalX2.0.12. The identity or similarity of the residues is represented by (*), (:), and (.). The residues in the active site are colored in red. The different residues in the entryway of active site between cel7482 and cel3623 are colored in green and underline. Figure S4. Alignment of cel36 and 3PZT proteins. The amino acid sequences of cel36 and 3PZT were aligned with ClustalX2.0.12. The identity or similarity of the residues is represented by (*), (:), and (.). The residues in the active site are colored in red.
About this article
Cite this article
Yang, C., Xia, Y., Qu, H. et al. Discovery of new cellulases from the metagenome by a metagenomics-guided strategy. Biotechnol Biofuels 9, 138 (2016). https://doi.org/10.1186/s13068-016-0557-3
- Glycoside hydrolase