Isolation and characterization of a novel endo-β-1,4-glucanase from a metagenomic library of the black-goat rumen

ABSTRACT

The various types of lignocellulosic biomass found in plants comprise the most abundant renewable bioresources on Earth. In this study, the ruminal microbial ecosystem of black goats was explored because of their strong ability to digest lignocellulosic forage. A metagenomic fosmid library containing 115,200 clones was prepared from the black-goat rumen and screened for a novel cellulolytic enzyme. The KG35 gene, containing a novel glycosyl hydrolase family 5 cellulase domain, was isolated and functionally characterized. The novel glycosyl hydrolase family 5 cellulase gene is composed of a 963-bp open reading frame encoding a protein of 320 amino acid residues (35.1 kDa). The deduced amino acid sequence showed the highest sequence identity (58%) for sequences from the glycosyl hydrolase family 5 cellulases. The novel glycosyl hydrolase family 5 cellulase gene was overexpressed in Escherichia coli. Substrate specificity analysis revealed that this recombinant glycosyl hydrolase family 5 cellulase functions as an endo-β-1,4-glucanase. The recombinant KG35 endo-β-1,4-glucanase showed optimal activity within the range of 30-50 °C at a pH of 6-7. The thermostability was retained and the pH was stable in the range of 30-50 °C at a pH of 5-7.

Keywords: Black goats; Endo-β-1,4-glucanase; Glycosyl hydrolase family 5 (GH5); Metagenomic library; Rumen

Introduction

Cellulose is a major polysaccharide compound in the plant cell wall and is composed of repeating cellobiose units linked via β-1,4-glucanase-glycosidic bonds. It is one of the most abundant renewable bioresources in nature.¹ The enzyme cellulase can be categorized into three major classes based on its catalytic action: endo-β-1,4-glucanases (EC 3.2.1.4), cellobiohydrolases (EC 3.2.1.91), or ^TM-glucosidases (EC 3.2.1.21). Endo-β-1,4-glucanases randomly attack internal amorphous sites in cellulose polymers, generating new chain ends. Cellobiohydrolases remove cellobiose from reducing or non-reducing ends of cellulose polymers. ^TM-Glucosidases hydrolyze cellodextrins and cellobioses to glucoses.² Cellulases have generated commercial interest in various sectors, including the food, energy, pulp, and textile industries.³ Among the three types of cellulases, endo-β-1,4-glucanase can release smaller cellulose fragments of random length.⁴ Consequently, endo-β-1,4-glucanase has a biotechnological potential in various industrial applications. This enzyme has been industrially applied in biomass waste management, pulp and paper deinking, and textile biopolishing.⁵ In addition, this enzyme can be practically used to produce better animal feeds, improve beer brewing, decrease the viscosity of ^TM-glucan solutions, and improve biofinishing in the textile industry.⁵

Cellulases can be produced using a wide range of microorganisms, plants, and animals. Symbiotic microorganisms in the rumen of herbivores can hydrolyze and ferment cellulosic polymers, thereby allowing the host to obtain energy from the indigestible polymers.⁶ Therefore, rumen microbial communities in ruminants are attractive sources of cellulases because these communities have adapted to utilization of lignocellulosic plant biomass.^6-9

Although the ruminal microbial ecosystem is highly diverse, the majority of rumen microorganisms are unculturable.⁶ It is estimated that less than 1% of microbial species in nature have ever been cultured.¹⁰ A culture-independent method has been developed to overcome this problem. This approach involves isolating metagenomic DNA from the environment and cloning this DNA into a vector to generate a metagenomic library. The library is subsequently screened for genes of interest via the detection of genotypic or phenotypic biomarkers using DNA-DNA hybridization, polymerase chain reaction (PCR), or enzymatic assay techniques.¹¹

The Korean black goat is a grazing herbivore with a dietary preference for herbaceous and woody dicots, such as forbs, shrub leaves, and stems.¹² Black goats are classified as small ruminants, and their foregut digestion allows feed particles to move more rapidly through the alimentary tract, promoting higher herb intake in comparison with other larger herbivores.¹³ Therefore, the objectives of this study were to isolate and functionally characterize a novel cellulase metagenomically derived from the rumen contents of black goats.

Materials and methods

Construction and functional screening of a metagenomic fosmid library

All animal studies and standard operating procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of the National Institute of Animal Science, Rural Development Administration, Suwon, Korea (No. 2009-007, D-grade, surgery). Three 18-month-old Korean black goats were raised at the National Institute of Animal Science and were freely fed rice straw and mineral supplements 30 d prior to the experiment to maximize cellulolytic adaptation of microorganisms in their rumens. Rumen contents of the three black goats were collected immediately after slaughter. Metagenomic DNA was extracted from the contents using hexadecyl trimethyl ammonium bromide and sodium dodecyl sulfate (CTAB-SDS).¹⁴ DNA fragments ranging from 10 to 50 kb were obtained by partial digestion with the restriction enzyme HindIII to construct a metagenomic fosmid library by means of the CopyControl^TM Fosmid Library Construction Kit (Epicentre, Madison, WI). The technique was based on the pCC1FOS vector, according to the manufacturer's instructions. Escherichia coli DH5α^TM (Life Technologies, Carlsbad, CA) transformants, carrying various recombinant fosmids, were transferred to 384 well plates. To determine an average insert size of the fosmid library, a total of 36 clones were randomly selected and digested with the restriction enzyme NotI.

A fraction of the library was spread on Luria-Bertani (LB) agar plates supplemented with chloramphenicol (15 µg/mL) and 0.2% carboxymethyl cellulose (CMC; Sigma-Aldrich, St. Louis, MO), which is a substrate for screening for endo-β-1,4-glucanase activity. After incubated and stained as described previously,¹⁵ fosmid clones with cellulolytic activity were detected upon the formation of clear zones around the colonies. Activity-based metagenomic screening yielded a single fosmid clone (Ad125D08) that formed one of the largest clear zones.

Shotgun sequencing of the enzyme-positive fosmid clone

Shotgun sequencing and molecular analyses were performed to localize the cellulase gene within the enzyme-positive fosmid clone and to subclone the gene for effective biochemical characterization. Bacteria carrying the cellulolytic fosmid clone (Ad125D08) were cultured in 100 mL of the LB medium containing 25 µg/mL of chloramphenicol at 37 °C for 20 h. After the cells were harvested, the DNA was extracted using a Qiagen Plasmid Midi Kit (Qiagen, Valencia, CA) and subjected to shotgun sequencing that comprised two approaches: conventional Sanger sequencing and pyrosequencing. For shotgun sequencing with Sanger methodology, 15 µg of DNA from the cellulolytic fosmid clone was randomly fragmented into sections measuring 2-3 kb in size. Fragmentation was performed using the HydroShear DNA shearing device (Genomic Solution, Waltham, MA) under the following processing conditions: sample volume 200 µL, speed code 11, and 20 shearing cycles. Small fragments were removed and were subsequently repaired with DNA polymerase and the polynucleotide kinase method (BKL Kit; TaKaRa, Shiga, Japan). The prepared DNA fragments were ligated into the SmaI site of pUC19. The ligation products were introduced into E. coli DH10B cells by electroporation. Approximately 192 recombinant plasmids in the shotgun DNA library were randomly selected for sequencing. Each plasmid DNA was bidirectionally sequenced as described previously.¹⁶ The sequencing data was assembled using Phred and Phrap software (University of Washington, Seattle, WA). One pyrosequencing run was conducted on GS Junior sequencing system (Roche Diagnostics, Oakland, CA) as described previously.¹⁷ The pyrosequencing data were assembled using GS De Novo Assembler, version 2.7.

In silico sequence analyses of the enzyme-positive fosmid clone

All contigs from conventional shotgun sequencing with Sanger methodology and from the pyrosequencing assembly were combined to construct a full metagenomic-insert-DNA sequence using the SeqMan software in Lasergene version 7 (DNASTAR, Madison, WI). Gene prediction in each contig was performed using MetaGeneMark with default parameters.¹⁸ Functional annotation was performed by searching Pfam domains on the Pfam website (http://pfam.sanger.ac.uk).¹⁹ Amino acid sequences were deduced from the nucleotide sequences of the enzyme-positive fosmid clone (Ad125D08) and were analyzed using BlastP on the National Center for Biotechnology Information (NCBI) database to search for proteins with similarity to the amino acid sequences derived from the identified open reading frames (ORFs). Multiple alignments among amino acid sequences were deduced from the ORFs, and the sequences of other proteins were performed using ClustalW software (http://www.ebi.ac.uk/Tools/msa/clustalw2/). The PROSITE database was used to predict the module structures and signature sequences of the enzyme (http://prosite.expasy.org/). The N-terminal regions were analyzed using Signal IP software for a signal peptide sequence (http://www.cbs.dtu.dk/services/SignalP/) and by HMMTOP 2.0 software for transmembrane regions (http://www.enzim.hu/hmmtop/html/submit.html).

Heterologous expression and purification of the recombinant KG35 endo-β-1,4-glucanase

The KG35 gene (the complete ORF of the putative cellulase gene) was amplified from the enzyme-positive fosmid clone (Ad125D08) by PCR. The amplicon was ligated into the expression vector pET24 (Novagen, Darmstadt, Germany), and the resulting recombinant vector was introduced into E. coli Rosetta-gami^TM (Novagen, Darmstadt, Germany). Additional screening for cellulolytic activity of the enzyme produced by the transformed E. coli cells was performed using the same CMC agar plate method as described previously.¹⁵ The cellulolytic E. coli cells were grown at 37 °C and the transgene expression was induced with 0.2 mM isopropyl ^TM-d-1-thiogalactopyranoside (IPTG) as described previously.⁷ After incubation for 6 h at 26 °C, the cells were harvested and disrupted by sonication as described previously.⁹ The overexpressed protein was purified from the resulting cell-free extracts using the Chelating Excellose^® Spin Kit (Bioprogen, Daejeon, South Korea). A 2.0-µL aliquot of the purified protein was inoculated into the CMC agar plates to confirm the enzymatic activity. The proteins were separated by SDS polyacrylamide gel electrophoresis (SDS-PAGE), according to Laemmli (1970),²⁰ and the protein concentrations were determined by the Bradford Protein Assay Kit (Bio-Rad, Hercules, CA).

Zymogram analysis of the recombinant KG35 endo-β-1,4-glucanase was performed using a 10% polyacrylamide gel containing 0.1% CMC under denaturing conditions as described previously.²¹

Enzymatic characterization of the recombinant KG35 endo-β-1,4-glucanase

Specific activities of the purified recombinant KG35 endo-β-1,4-glucanase toward various substrates were evaluated under the following standard conditions. The substrates used in this study were CMC, barley glucan, laminarin, avicel, xylan from birchwood, ^TM-mannan, p-nitrophenyl-d-glucopyranoside, and p-nitrophenyl-d-cellobioside. In all the assays, 10 µL of a properly diluted enzyme solution were added to 200 µL of 0.1 M sodium phosphate buffer (pH 7), containing 1% (w/v) of each substrate. The enzymatic reactions were allowed to proceed at 50 °C (different reaction duration): 15 min for CMC, glucan, laminarin, and avicel; 10 min for xylan from birchwood and ^TM-mannan; and 20 min for p-nitrophenyl-d-glucopyranoside and p-nitrophenyl-d-cellobioside. Specific activity was defined as the number of activity units per milligram of protein. One unit of enzymatic activity was defined as the amount of enzyme necessary to release 1 µmol of reducing sugar per minute. Additionally, 1 unit of enzymatic activity producing chromogenic sugar derivatives (p-nitrophenyl-d-glucopyranoside and p-nitrophenyl-d-cellobioside) as the substrates was defined as the amount of the enzyme required to release 1 µmol of p-nitrophenol per minute. The amounts of released reducing sugar and p-nitrophenol were determined by measuring the absorbances at 540 or 400 nm, respectively. All assays were performed in triplicate.

Optimal pH was determined by assaying the enzymatic activity of the purified recombinant, KG35 endo-β-1,4-glucanase, at 50 °C for 15 min at pH ranging from 4 to 10. To determine the optimal temperature, enzymatic activity was analyzed in 0.1 M sodium phosphate buffer (pH 7) for 15 min at a temperature ranging from 20 to 70 °C. All assays were performed in triplicate. The pH stability was assessed after the enzyme was pre-incubated at 4 °C for 24 h at pH ranging from 4 to 10. The thermal stability was assessed after the enzyme was pre-incubated for 1 h in 0.1 M sodium phosphate buffer (pH 7) at various temperatures mentioned above. All assays were performed in triplicate.

Nucleotide sequence accession number

The nucleotide sequence of the enzyme-positive fosmid clone (Ad125D08), including the complete ORF of the KG35 gene (a novel GH5 cellulase gene), was deposited in the GenBank database under the accession number KJ631389.

Results

Construction and functional screening of the metagenomic fosmid library

A metagenomic fosmid library containing 115,200 clones was constructed from the rumen microorganisms of black goats. NotI restriction enzyme analyses revealed that the average insert size was typically 31 kb (data not shown). A total of 155 independent fosmid clones, displaying carboxymethyl cellulase (CMCase) activity, were identified. Among them, a single fosmid clone (Ad125D08) formed one of the largest clear zones on the CMC-containing screening plates and was selected for further analysis.

Conventional shotgun sequencing with the Sanger method was performed under 2.7-fold coverage, and ∼300 kb of pyrosequencing reads were produced. Consequently, one contig of 42.453 kb was obtained without the pCC1FOS vector sequence. Thirty-four ORFs were predicted in the contig sequence of fosmid clone Ad125D08, and an ORF containing a cellulase domain (KG35 gene) was identified and localized within the clone.

In silico sequence analysis of the putative KG35 endo-β-1,4-glucanase

Sequence analysis of the putative KG35 gene revealed a single ORF of 963 bp encoding a polypeptide of 320 amino acid residues, with an estimated molecular mass of 35.1 kDa and pI of 5. Signal IP and HMMTOP 2.0 software programs could not detect any known signal peptide sequence or transmembrane region in the KG35 protein. The module structure search in the PROSITE database predicted that the KG35 protein contains a GH5 cellulase catalytic domain with a conserved (amino acid) signature sequence VLYEICNEPN, and its second Glu residue was assumed to play a catalytic role (Fig. 1).^22,23

Fig. 1
Multiple alignment of the KG35 endo-β-1,4-glucanase GH5 domain with other homologous GH5 protein domains.

The deduced amino acid sequence of a GH5 family domain in KG35 endo-β-1,4-glucanase was aligned with those of selected homologous enzymes from the following microorganisms: Eubacterium cellulosolvens (WP_051527028.1); Agathobacter rectalis (WP_055222932.1); Roseburia inulinivorans (WP_055169535.1); A. rectalis (WP_012741522.1); Ruminococcus torques (WP_055163516.1); Eubacterium rectale CAG:36 (CDC71888.1); R. inulinivorans (WP_055040325.1); and Lachnospiraceae bacterium C6A11 (WP_051646390.1). Signature sequences of GH5 domains are shaded in gray. The Glu residues under the star symbols indicate the catalytic sites.

A BlastP search, comparing the amino acid sequence of the KG35 protein with those of other proteins registered in the GenBank database, revealed that this protein has the closest (58%) sequence identity to a hypothetical protein from Eubacterium cellulosolvens (WP_051527028.1) and shares even less sequence identity with other GH5 family protein homologs from various microbial species (Table 1). The amino acid sequences of these homologous proteins possess characteristic features conserved among the GH5 family of enzymes (Fig. 1).

Table 1
Amino acid sequence identity between the novel KG35 endo-β-1,4-glucanase and other (Homologous) proteins.

Protein name (Organism name)	GeneBank accession number	Query coverage (%)a	Identity (%)a
Hypothetical protein (Eubacterium cellulosolvens)	WP_051527028.1	98	58
Glycosyl hydrolase family 5 protein (Agathobacter rectalis)	WP_055222932.1	98	52
Glycosyl hydrolase family 5 protein (Roseburia inulinivorans)	WP_055169535.1	98	51
Endo-1,4-β-glucanase (A. rectalis)	WP_012741522.1	98	51
Glycosyl hydrolase family 5 protein (Ruminococcus torques)	WP_055163516.1	98	51
Endoglucanase (Eubacterium rectale CAG:36)	CDC71888.1	98	51
Glycosyl hydrolase family 5 protein (R. inulinivorans)	WP_055040325.1	98	50
Hypothetical protein (Lachnospiraceae bacterium C6A11)	WP_051646390.1	98	50

• ^a
  The novel KG35 endo-β-1,4-glucanase gene was used for BlastP searches using the 98% query coverage cutoff value and 50% identity cutoff value.

Purification and SDS-PAGE analyses of the recombinant KG35 endo-β-1,4-glucanase

The putative endo-β-1,4-glucanase gene (KG35) was cloned and overexpressed in E. coli and the resulting recombinant enzyme was purified, as described in the materials and methods section. SDS-PAGE analysis of the purified recombinant KG35 endo-β-1,4-glucanase showed a single band on the gel, with the apparent molecular mass of approximately 35 kDa (Fig. 2). Zymogram analysis revealed a single band corresponding to an active endo-β-1,4-glucanase of the same molecular mass as the band (Fig. 2).

Fig. 2
SDS-PAGE analyses of recombinant KG35 endo-β-1,4-glucanase produced in E. coli Rosetta-gami^TM.

Lane M: protein molecular weight markers; lane 1: total cell-free extracts before induction; lane 2: the soluble fraction after the induction; lane 3: purified KG35 endo-β-1,4-glucanase; lane 4: zymogram analysis of KG35 endo-β-1,4-glucanase.

Characterization of the purified recombinant KG35 endo-β-1,4-glucanase

Substrate specificity toward various substrates was determined (Table 2). The endo-β-1,4-glucanase showed strong specific activities toward CMC (8.40 U/mg) and barley glucan (5.40 U/mg). In contrast, this enzyme showed only trace enzymatic activity toward laminarin (0.48 U/mg), p-nitrophenyl-d-cellobioside (0.51 U/mg), and avicel (0.12 U/mg). In addition, the endo-β-1,4-glucanase could not hydrolyze p-nitrophenyl-d-glucopyranoside, xylan from birchwood, or ^TM-mannan.

Table 2
Substrate specificity of the novel KG35 endo-β-1,4-glucanase toward various substrates.

Substrate	Main linkage types	Specific activity (U/mg)
Carboxylmethylcellulose (CMC)	β-1,4-Glucan	8.40 ± 0.03
Barley glucan	β-1,3/1,4-Glucan	5.40 ± 0.19
P-nitrophenyl-d-cellobioside	β-1,4-Glucan	0.51 ± 0.01
Laminarin	β-1,3/1,6-Glucan	0.48 ± 0.04
Avicel	β-1,4-Glucan	0.12 ± 0.02
Xylan from birchwood	β-1,4-Xylan	0.00 ± 0.00
β-Mannan	β-1,4-Mannan	0.00 ± 0.00
P-nitrophenyl-d-glucopyranoside	β-1,4-Glucan	0.00 ± 0.00

The optimal pH and temperature of the recombinant enzyme were determined under various pH and temperature conditions (Fig. 3). The optimal pH and temperature of the recombinant KG35 endo-β-1,4-glucanase toward CMC was approximately 6-7 and 30-50 °C, respectively.

Fig. 3
Effects of pH and temperature on the activity of recombinant KG35 endo-β-1,4-glucanase. (A) The optimal enzyme activity was measured at a pH range of 4-10. The buffers were 0.1 M sodium acetate (pH 4-6, ●), 0.1 M sodium phosphate (pH 6-8, ○), 0.1 M Tris-HCl (pH 8-9, ▼), and 0.1 M glycine-NaOH (pH 9-10, ▵). (B) The optimal enzyme activity was measured at a temperature range of 20-70 °C. The error bars represent the standard deviation of triplicate measurements.

In addition, the enzymatic stability was tested under various pH and temperature conditions (Fig. 4). The endo-β-1,4-glucanase retained more than 70% of activity after 24 h of incubation at pH 5-7. The same enzyme showed thermostability after a 1-h incubation at 30-50 °C; however, it declined in activity by greater than 60% after a 1-h incubation above 50 °C. Our empirical data collectively indicates that the recombinant enzyme retained its pH and temperature stability approximately within its optimal pH and temperature ranges, respectively.

Fig. 4
Effects of pH and temperature on the stability of the recombinant KG35 endo-β-1,4-glucanase. (A) The pH stability as determined by measuring the residual activity after the enzyme was preincubated at 4 °C for 24 h at the indicated pH levels from 4 to 10. The buffers were 0.1 M sodium acetate (pH 4-6, ●), 0.1 M sodium phosphate (pH 6-8, ○), 0.1 M Tris-HCl (pH 8-9, ▼), and 0.1 M glycine-NaOH (pH 9-10, ▵). (B) The temperature stability as determined by measuring the residual activity after the enzyme was pre-incubated at pH 7 for 1 h over the temperature range of 20-70 °C. The error bars represent the standard deviation of triplicate measurements.

Discussion

It has been reported that the GH5 family of cellulases contain the conserved amino acid signature sequence pattern: [LIV]-[LIVMFYWGA](2)-[DNEQG]-[LIVMGST]-{SENR}-N-E-[PV]-[RHDNSTLIVFY].²⁴ The conserved signature sequence of the KG35 protein could be located in the above-mentioned GH5 family signature sequence pattern. This result indicated that the KG35 protein exhibited a strong relation to the GH5 family of cellulases. Furthermore, the KG35 endo-β-1,4-glucanase protein exhibited the closest (58%) amino acid sequence similarity to its GH5 homologous enzymes, implying that the KG35 endo-β-1,4-glucanase protein has never been identified or characterized, thus is a novel amino acid sequence.

The CMC and barley glucan, which are primarily based on β-1,4-glucanase-glycosidic bonds and have amorphous structures, were found to be the most favorable substrates for the recombinant KG35 endo-β-1,4-glucanase. Amorphous laminarin, containing ^TM-1,3/1,6-glycosidic bonds, and hemicelluloses, including xylan from birchwood and ^TM-mannan, were not specific substrates of the recombinant KG35 endo-β-1,4-glucanase. Although avicel mostly contains β-1,4-glucanase-glycosidic bonds, it was not efficiently hydrolyzed by this enzyme. Avicel contains β-1,4-glucanase-glycosidic bonds akin to other soluble cellulosic substrates, and its crystalline structure can undergo efficient hydrolysis when a cellulose-binding domain is involved in the hydrolysis reaction.^25,26 Nevertheless, sequence analysis could not detect any cellulose-binding domain in the KG35 endo-β-1,4-glucanase; therefore the absence of this domain may explain the weak enzymatic activity toward avicel. Two synthetic substrates used in this study, p-nitrophenyl-d-cellobioside and p-nitrophenyl-d-glucopyranoside, are synthetic substrates for exoglucanase and ^TM-glucosidase, respectively.^27,28 Therefore, the recombinant KG35 endo-β-1,4-glucanase exerted slight exoglucanase activity and no ^TM-glucosidase activity, respectively. All of these results pertaining to substrate specificity indicate that KG35 cellulase is most likely to hydrolyze the β-1,4-glucanase-glycosidic bonds of the amorphous regions of cellulose and consequently function as an endo-β-1,4-glucanase.

The recombinant KG35 cellulase was maximally active and enzymatically stable at a pH range of 6-7 and at a temperature range of 30-50 °C. The optimum and stability characteristics of the recombinant enzyme were similar to those of other cellulose enzymes from rumenal microorganisms of other ruminants. Other cellulases from rumenal microorganisms of the cow were reported to have the optimal pH range of 5-7 and an optimal temperature of approximately 50 °C.^9,29,30 The cellulase from rumen-resident bacteria was reported to retain 70% of its activity in the pH range from 5 to 7 and in a temperature range from 30 to 50 °C.³⁰

The ruminal pH of small ruminants ranges from 5 to 7, under normal dietary conditions, and the ruminal temperature can vary from 38 to 41 °C.³¹ The cellulolytic reactions of ruminal-bacteria cellulase have been reported to be relatively more resistant at the pH of 6.5-7 in culture, but the reactions were weakened as the culture pH decreased.³² KG35 endo-β-1,4-glucanase may be well-adapted to the ruminal environments of black goats because the pH and temperature conditions under which the recombinant KG35 endo-β-1,4-glucanase showed relatively higher activity and stability were in agreement with the ruminal pH and temperature conditions that occurs naturally in ruminants, including goats. The above observation can explain why the recombinant KG35 endo-β-1,4-glucanase functions better at mildly acidic/neutral pH and a moderate temperature than at an acidic/alkaline pH or strongly unphysiological temperature.

Cellulase can effectively convert cellulose, an anti-nutritional component, in animal feedstocks into its hydrolysates, nutritive ingredients, yielding improved performance of animals.⁵ Therefore, applications of different cellulases from various sources have attracted considerable attentions in the feed industry. KG35 endo-β-1,4-glucanase may be exploited as a feed additive to improve feed utilization of ruminants, because the enzyme may be well-adapted to the ruminal environments. Moreover, the enzyme can be potentially used to make cellulosic ethanol, second generation biofuel, from cellulosic feedstocks such as grass, wood, and crop residues.

This study presents novel results pertaining to the construction and screening of a metagenomic library from the rumenal microorganisms of black goats for identification and characterization of a novel cellulolytic enzyme. A total of 155 cellulolytic fosmid clones were selected from the metagenomic library. Among them, the KG35 gene was successfully isolated as a clone with the strongest cellulose-degrading activity. The KG35 gene was found to encode a novel member of the GH5 family of cellulases, and this novel gene was efficiently translated and overexpressed in E. coli Rosetta-gami^TM. Substrate specificity analysis revealed that the recombinant KG35 cellulase functions as an endo-β-1,4-glucanase. The recombinant KG35 endo-β-1,4-glucanase optimally functions in the pH range of 6-7 and temperature range of 30-50 °C. The endo-β-1,4-glucanase retained enzymatic stability in the pH range of 5-7 and the temperature range of 30-50 °C.

• 1
  These authors contributed equally to this work.

Acknowledgments

This work was carried out with the support of "Cooperative Research Programs for Agriculture Science and Technology Development (Project Nos. PJ006649, PJ008701 and PJ011163)", Rural Development Administration, Republic of Korea. J.-S. Lee (for in silico experiment design) and K.-S. Kim (for other experiment design) are co-corresponding authors of this work.

References

• ¹ Bhat MK, Bhat S. Cellulose degrading enzymes and their potential industrial applications. Biotechnol Adv. 1997;15:583-620.
• ² Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS. Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev. 2002;66(3):506-577.
• ³ Bhat MK. Cellulases and related enzymes in biotechnology. Biotechnol Adv. 2000;18(5):355-383.
• ⁴ Sweeney MD, Xu F. Biomass converting enzymes as industrial biocatalysts for fuels and chemicals: recent developments. Catalysts 2012;2:244-263.
• ⁵ Kuhad RC, Gupta R, Singh A. Microbial cellulases and their industrial applications. Enzyme Res. 2011;2011:1-10.
• ⁶ Krause DO, Denman SE, Mackie RI, et al. Opportunities to improve fiber degradation in the rumen: microbiology, ecology, and genomics. FEMS Microbiol Rev. 2003;27(5):663-693.
• ⁷ Feng Y, Duan CJ, Pang H, et al. Cloning and identification of novel cellulase genes from uncultured microorganisms in rabbit cecum and characterization of the expressed cellulases. Appl Microbiol Biotechnol. 2007;75(2):319-328.
• ⁸ Duan CJ, Xian L, Zhao GC, et al. Isolation and partial characterization of novel genes encoding acidic cellulases from metagenomes of buffalo rumens. J Appl Microbiol 2009;107(1):245-256.
• ⁹ Nguyen NH, Maruset L, Uengwetwanit T, et al. Identification and characterization of a cellulase-encoding gene from the buffalo rumen metagenomic library. Biosci Biotechnol Biochem. 2012;76(6):1075-1084.
• ¹⁰ Amann RI, Ludwig W, Schleifer KH. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol Rev. 1995;59(1):143-169.
• ¹¹ Mirete S, Morgante V, González-Pastor JE. Functional metagenomics of extreme environments. Curr Opin Biotechnol 2016;38:143-149.
• ¹² Lee ID, Lee HS. Study on the food habits of Korean native goat (Capra hircus) fed with various roughage sources. J Kor Soc Grassl Forage Sci 2008;38:445-452.
• ¹³ Hofmann RR, Stewart DRM. Grazer or browser: a classification based on the stomach structure and feeding habits of East African ruminants. Mammalia 1972;36:226-240.
• ¹⁴ Zhou J, Bruns MA, Tiedje JM. DNA recovery from soils of diverse composition. Appl Environ Microbiol 1996;62(2):316-322.
• ¹⁵ Geng A, Zou G, Yan X, et al. Expression and characterization of a novel metagenome-derived cellulase Exo2b and its application to improve cellulase activity in Trichoderma reesei. Appl Microbiol Biotechnol. 2012;96(4):951-962.
• ¹⁶ Lee KT, Byun MJ, Kang KS, et al. Neuronal genes for subcutaneous fat thickness in human and pig are identified by local genomic sequencing and combined SNP association study. PLoS ONE. 2011;6(2):e16356.
• ¹⁷ Choi JW, Liao X, Stothard P, et al. Whole-genome analyses of Korean native and Holstein cattle breeds by massively parallel sequencing. PLOS ONE. 2014;9(7):e101127.
• ¹⁸ Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010;38(12):e132.
• ¹⁹ Punta M, Coggill PC, Eberhardt RY, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40:D290-D301.
• ²⁰ Laemmli UK. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970;227(5259):680-685.
• ²¹ Lee CM, Lee YS, Seo SH, et al. Screening and characterization of a novel cellulase gene from the gut microflora of Hermetia illucens using metagenomic library. J Microbiol Biotechnol. 2014;24(9):1196-1206.
• ²² Henrissat B, Callebaut I, Fabrega S, Lehn P, Mornon JP, Davies G. Conserved catalytic machinery and the prediction of a common fold for several families of glycosyl hydrolases. Proc Natl Acad Sci U S A. 1995;92(15):7090-7094.
• ²³ Wang Q, Tull D, Meinke A, et al. Glu280 is the nucleophile in the active site of Clostridium thermocellum CelC, a family A endo-beta-1,4-glucanase. J Biol Chem. 1993;268(19):14096-14102.
• ²⁴ Cheema TA, Jirajaroenrat K, Sirinarumitr T, Rakshit SK. Isolation of a gene encoding a cellulolytic enzyme from swamp buffalo rumen metagenomes and its cloning and expression in Escherichia coli. Anim Biotechnol. 2012;23(4):261-277.
• ²⁵ Le Costaouëc T, Pakarinen A, Várnai A, Puranen T, Viikari L. The role of carbohydrate binding module (CBM) at high substrate consistency: comparison of Trichoderma reesei and Thermoascus aurantiacus Cel7A (CBHI) and Cel5A (EGII). Bioresour Technol 2013;143:196-203.
• ²⁶ Zhao L, Pang Q, Xie J, Pei J, Wang F, Fan S. Enzymatic properties of Thermoanaerobacterium thermosaccharolyticum ^TM-glucosidase fused to Clostridium cellulovorans cellulose binding domain and its application in hydrolysis of microcrystalline cellulose. BMC Biotechnol. 2013;13:101.
• ²⁷ Mallerman J, Papinutti L, Levin L. Characterization of ^TM-glucosidase produced by the white rot fungus Flammulina velutipes. J Microbiol Biotechnol. 2015;25(1):57-65.
• ²⁸ Jing L, Zhao S, Xue J, et al. Isolation and characterization of a novel Penicillium oxalicum strain Z1-3 with enhanced cellobiohydrolase production using cellulase-hydrolyzed sugarcane bagasse as carbon source. Ind Crops Prod. 2015;77:666-675.
• ²⁹ Zhao S, Wang J, Bu D, et al. Novel glycoside hydrolases identified by screening a Chinese Holstein dairy cow rumen-derived metagenome library. Appl Environ Microbiol. 2010;76(19):6701-6705.
• ³⁰ Gong X, Gruninger RJ, Qi M, et al. Cloning and identification of novel hydrolase genes from a dairy cow rumen metagenomic library and characterization of a cellulase gene. BMC Res Notes. 2012;5:566-576.
• ³¹ Mackie RI, McSweeney CS, Aminov RI. Rumen. In: Encyclopedia of Life Sciences Chichester, England: John Wiley and Sons; 2013.
• ³² Russell JB, Dombrowski DB. Effect of pH on the efficiency of growth by pure cultures of rumen bacteria in continuous culture. Appl Environ Microbiol. 1980;39(3):604-610.

Publication Dates

• Publication in this collection
  Oct-Dec 2017

History

• Received
  10 Aug 2016
• Accepted
  1 Mar 2017