- Research article
- Open Access
CAZyChip: dynamic assessment of exploration of glycoside hydrolases in microbial ecosystems
BMC Genomics volume 17, Article number: 671 (2016)
Microorganisms constitute a reservoir of enzymes involved in environmental carbon cycling and degradation of plant polysaccharides through their production of a vast variety of Glycoside Hydrolases (GH). The CAZyChip was developed to allow a rapid characterization at transcriptomic level of these GHs and to identify enzymes acting on hydrolysis of polysaccharides or glycans.
This DNA biochip contains the signature of 55,220 bacterial GHs available in the CAZy database. Probes were designed using two softwares, and microarrays were directly synthesized using the in situ ink-jet technology. CAZyChip specificity and reproducibility was validated by hybridization of known GHs RNA extracted from recombinant E. coli strains, which were previously identified by a functional metagenomic approach. The GHs arsenal was also studied in bioprocess conditions using rumen derived microbiota.
The CAZyChip appears to be a user friendly tool for profiling the expression of a large variety of GHs. It can be used to study temporal variations of functional diversity, thereby facilitating the identification of new efficient candidates for enzymatic conversions from various ecosystems.
The degradation of polysaccharides such as cellulose, chitin, starch and glycogen is an essential feature of carbon cycle in the biosphere, a process that requires the contribution of various microorganisms that together deploy an arsenal of carbohydrate-degrading enzymes. Plant cell walls (PCWs) are composed of a composite network of macromolecules, including polysaccharides and lignin. The major polysaccharide in most plant cell walls is cellulose, which is composed of β-1,4 linked glucose polymers that interconnect through strong hydrogen bonds, forming crystalline microfibrils that are very stable. Cellulose is further embedded in a 3 D matrix composed of hemicelluloses, pectin and lignin  resistant to degradation. Compared to cellulose, hemicelluloses are heteropolymers that are variable in both chemical composition and structure, with heteroxylans and mannans being the two major categories of hemicelluloses in PCWs . The exact compositional and structural features of hemicelluloses are dependent on a number of determinants, including the botanical origin of the plant, and also the pedoclimatic conditions prevailing at the time of growth [13, 14, 62]. Therefore, microorganisms that are responsible for biomass degradation are faced with a formidable task, which they achieve through the deployment of complex arsenals of enzymes .
Among the key PCW-degrading enzymes that are produced by microorganisms, the glycoside hydrolases (GH) and the carbohydrate esterases (CE) belong to a wide class of enzymes that modify, synthesize or hydrolyze carbohydrates: Carbohydrate Active enZymes, or CAZymes (ref CAZy). The CAZymes are prominent and highly diverse and have been identified in all taxa, representing typically 1–5 % of the predicted coding sequences in their genomes . These proteins are expressed by microorganisms inhabiting almost all ecological niches (e.g., soil, marine environment and digestive tracts), where they participate in carbon cycling. The strategies of carbohydrate-degradation are often different at both the level of the microbial community and of individual microorganisms .
GH and CE can be encoded by multigenic operon-like clusters , such as Sus system [15, 51], that have been designated as Polysaccharide Utilization Loci in Bacteroidetes species [41, 44]. Evidence so far reveals that the proteins produced by such clusters display functional interplay with CAZyme components, displaying synergy on complex substrates [1, 48, 53]. In some anaerobic biomass-degrading bacteria, CAZymes, such as cellulases and hemicellulases, are arranged on cellulosomes, which are extracellular, cell-bound multi-enzyme complexes. In cellulosomes, the enzyme components are brought into close physical proximity, thus optimizing their synergistic actions and enhancing their biomass-degrading ability [3, 20].
GH and CE, and particularly those that are active on PCWs, are sought after for a wide range of industrial applications, including biorefining. In this field, the enzymes that are of particular interest include those active on cellulose (e.g., endoglucanases, EC 22.214.171.124, exoglucanases, EC 126.96.36.199 and EC 188.8.131.52) and on heteroxylans (e.g., endoxylanases, EC 184.108.40.206, β-D-xylosidases, EC 220.127.116.11 and α-L-arabinofuranosidases, EC; 18.104.22.168). Cellulose and hemicellulose yield monomeric sugars readily fermentable to produce alcohols, organic acids, or alkenes. The exploration of glycoside hydrolase (GH) diversity, and to a lesser extent CE can provide efficient biocatalysts and new insight into the different enzyme mechanisms that are used by microorganisms in biomass degradation. GHs have been used in many industries such as in paper production, textiles, detergents, feed and food [4, 33] as well as to promote healthy human nutrition and prevent diseases . In the last decade, cellulases and more recently hemicellulases have been considered for biorefining [23, 30]. The discovery of GHs has been considerably accelerated with the metagenomic and metatranscriptomic approaches, which allow the identification of new enzymes in an unprecedented manner.
GH exploration is largely facilitated by the existence of the CAZy database (CAZy; www.cazy.org). This database describes the families of enzymes that catalyze the breakdown, biosynthesis or modification of carbohydrates and glycoconjugates. In the CAZy database, GHs are classified into families based on amino acid sequence similarities and others conserved features [7, 25, 26, 39]. GH- are classified in 135 families and represent approximately 47 % of the entire database. (April 2016) . The vast majority of currently known GH are from bacterial origin.
DNA microarrays are widely used to profile gene expression and represent a relevant tool to study expression of key enzymes and monitor physiological changes of pure cultures or microbial communities [12, 18, 28, 42, 46, 50, 68]. This approach can also be useful to link microbial diversity to ecosystem processes and functions [22, 29, 67].
In this study, we developed the first microarray tool, termed CAZyChip, to quickly and accurately explore, at transcriptomic level, the GH composition of environmental samples. The CAZyChip provides snapshot views of the enzymes expressed by a single microorganism or more interestingly by microbial consortia derived from complex and various ecosystems. The biochip gives an opportunity to highlight enzyme cooperation along with the plant biomass degradation pathway. The present study demonstrates that the CazyChip represents a unique, robust and yet generic tool to dynamically analyze the expression of a large variety of GHs in parallel. The current version of this biochip allows the detection of 55,220 bacterial annotated GHs and contains the signatures of all bacterial GH in all families available to date in the CAZy database in addition to 53 CE sequences. The CAZy chip was validated using characterized enzymes from gut metagenomic libraries of different species, which were chosen for their known abilities to degrade plant cell walls. The encoding sequences of the enzymes of interest were recovered from microbiome of worm (Pontoscolex corethrurus), human, rumen, and termites these latter include fungus-growing (Pseudacanthotermes militaris), wood-feeding (Nasutitermes corniger), or soil-wood feeding (Termes hispaniolae). Furthermore, the developed biochip was tested to highlight the GH functional diversity of complex lignocellulolytic microbial communities, using a cow rumen-derived microbial consortium. The resulting biochip is able to test the GH functional diversity of complex microbial communities that present high metabolic and taxonomic diversity.
Custom microarray design
The design of oligonucleotides for the microarray was performed using either the Agilent e-Array online portal (https://earray.chem.agilent.com/earray/) or, when sequences were rejected by eArray, the ROSO software [16, 52]. When the design of 60-mers were impossible, a 40-mer or a pair of 25-mers associated with inert nucleotidic linkers was generated. For each targeted CAZyme gene (GH and CEs), three different 60-mer probes were designed and for each probe. The Agilent probe design algorithm assigned a BC score, which reflects uniqueness, secondary structure considerations, GC content and thermodynamic parameters, that predicts hybridization quality on the basis of their nucleotidic composition . Five grades of BC scores were defined and indicated the quality of the designed probes. These different scores were, from the best to the worst: BC_1, BC_2, BC_3, BC_4 and BC_Poor. A total of 180,000 probes, including 4848 Agilent internal positive or negative control probes, were selected and synthesized in situ, on a glass slide using Agilent SurePrint technology to obtain a high-density DNA microarray tool on 4x180 K format (Agilent Technologies, Massy, France) . The full description of the CAZyChip microarray has been deposited in the Gene Expression Omnibus (GEO) public database (GSE80173 study is at: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE80173).
Strains and growth conditions
Different GH cloned in plasmid or fosmid (pDest vector) were expressed by recombinant E. coli strains as previously described, [1, 2, 10, 34, 57, 59, 66]. Briefly, cultures were stopped at OD600nm between 0.4 and 0.6, and cells were harvested by centrifugation for 10 min at 5000 rpm at 4 °C. The supernatant was then discarded and the bacterial pellet immediately frozen at −80 °C before RNA extraction.
Microbial consortia analysis were performed on an anaerobic rumen-derived consortium RWS, which efficiently degrades lignocellulose, as reported by Lazuka et al. .
Availability of materials section
The GH gene sequences used in this study were deposited under the GenBank accession number: TxAbf CAA76421; THSAbf ABZ10760; CfXyn AEA30147; TM1225 AAD36300.1; Abn43a and Pm08 CCO20984.1; Abn43b CCO20993.1; Abf51b CCO20994.1; Pm06 HF548274; Pm13 CCO21046.1, Pm14 CCO21057.1, Pm15 CCO21059.1; Pm21 CCO21105.1; Pm25 CCO21110.1; Pm31 CCO21136.1; Pm41 CCO21355.1; Pm43 CCO21392.1;Pm55 CCO21443.1; Pm65 CCO21487.1; Pm66 CCO21489.1; Pm69 CCO21492.1; Pm80 CCO21560.1; Pm81 CCO21564.1; Pm83 CCO21640.1; Pm85 CCO21658.1; and Pm87 CCO21793.1.
Bacterial pellets were lysed with 1 mg/ml lysozyme (Sigma-Aldrich, Isle d’Abeau Chesnes, France) for 5 min at 25 °C, followed by Total RNA extraction using the RNeasy Mini Kit (Qiagen, Courtaboeuf, France) according to the manufacturer’s recommendations. RNA concentration and purity was evaluated by measuring the absorbance ratio at 260/280 nm and 260/230 nm using a Nanodrop spectrophotometer (Labtech, Palaiseau, France). The Ratio Integrity Number (RIN) was evaluated using 2100 Bioanalyzer® (Agilent Technologies, Massy, France) and only samples with a RIN greater than 8 were hybridized on the microarray.
Total RNA of rumen derived consortium was extracted in two steps from nitrogen frozen samples using the PowerMicrobiome RNA isolation kit (MoBio Laboratories, Carlsbad, CA, USA) . RNA purification was performed using AllPrep DNA/RNA minikit (Qiagen), according to the manufacturer’s recommendations.
Labelling and amplification of total mRNA
The One-Color Low Input Quick Amp WT Labeling Kit™ (Agilent Technologies, Massy, France) was used to amplify and label 100 ng of RNA according to the manufacturer’s recommendations. The labelling efficiency was checked using a NanoDrop spectrophotometer operating at 260 nm to quantify cRNA and at 550 or 660 nm to measure cyanine 3 (Cy3) and cyanine 5 (Cy5) dye incorporation, respectively. Labeling efficiency was calculated as indicated by the manufacturer’s protocol (ratio cyanine quantity / amount of RNA) and was above 6.
Microarray hybridization, washing and scanning
For each sample, 1650 ng of labeled and amplified cRNA was used for hybridization. The hybridization master mix was prepared according to manufacturer’s protocol (Agilent Technologies, Massy, France) and 100 μl were deposited onto a gasket slide, according to the Agilent Microarray Hybridization Chamber User Guide. Next, the active side of the microarray slide was placed on top of the gasket to form a properly aligned “sandwich slide pair”. The microarray slides were inserted into an Agilent Technology hybridization chamber then placed at 65 °C for 17 h with rotation at 10 rpm. After hybridization, the microarray was washed over a 1-min period, first using Gene Expression Wash Buffer 1 and then Gene Expression Wash Buffer 2 (Agilent Technologies, Massy, France) pre-warmed at 37 °C. After washing, the arrays were immediately scanned using an MS200 scanner (NimbleGen Roche Diagnostics, Meylan, France) with NimbleGen MS200 software v1.2 at 2 micron resolution.
The median signal of each spot in the hybridized arrays were determined and quantified using Feature Extraction software v22.214.171.124. The data from all the microarrays were normalized using the “limma” package function “normalizeQuantiles” and the “quantile” method [5, 56]. Normalization and statistical analyses of the data were performed using the Bioconductor packages (http://www.bioconductor.org) and R software v3.1.3. For each sample, the normalized fluorescence intensities of the three experimental replicates were analyzed and the mean values, standard deviations and correlation coefficients (%CV) were calculated. To determine whether probes were specific and target genes present, limma one way ANOVA test was carried out with False Discovery Rate adjusted p value < 0.05. Limma t test using “limma” package, was conducted to know in which comparison(s) this gene is differentially expressed (DE).
Analysis of mRNA levels by qRT-PCR
One microgram of RNA was used as template to generate cDNA using the High Capacity cDNA reverse transcriptase kit (Applied Biosystems, Life Technologies, Saint Aubin, France). The reverse transcription reaction (20 μl final volume) was performed for 10 min at 25 °C, and then 2 h at 37 °C. Quantitative real-time PCR (qRT-PCR) assays were performed using SsoFast EvaGreen Supermix (Bio-Rad, Marnes-La-Coquette, France) on the StepOne instrument (Applied Biosystems, Life Technologies, Saint Aubin, France). Primers were validated by testing qRT-PCR efficiency using standard curves (95 % efficiency 105 %) as described previously . Gene expression was quantified using the comparative Ct (threshold cycle) method. The RNA polymerase sigma S (rpoS) gene encoding the sigma factor sigma-38 was used as a reference to normalize the expression level of the targeted genes. Gene-specific primers sequences are described in Additional file 1: Table S1.
To design a generic microarray for the high-throughput detection of bacterial CAZymes mainly composed of GH’s, all of the bacterial GH protein sequences referenced in the CAZy database (www.cazy.org) up to January 2015 (133 families), were selected and their nucleotide sequences downloaded from the National Center for Biotechnology Information database (www.ncbi.nlm.nih.gov). We also selected sequences of interest obtained from human or termite guts and cow rumen metagenomic libraries created in our laboratory [2, 10]. The initial dataset used for probe design contained a total of 55,220 sequences and for each gene we designed three non-overlapping probes, with the aim to validate at least one probe per GH for use in a future prototype. With the e-array software, probe design has been possible on 55,012 sequences with a BC score attribution. This score reflects several criteria including the predicted hybridization quality, GC content and steric hindrances (Additional file 2: Table S2). A total of 56 % of probes displayed a BC_score of BC_1, 22 % of BC_2 reflecting the highest quality of predicted hybridization and a stable and consistent duplex with their targets. Only a small fraction of the probes were scored as BC_3 (11 %), BC_4 (11 %) and no BC_Poor were detected. Using the ROSO software we designed probes for the 208 of the remaining sequences.
The final CAZyChip was constructed using 180,000 probes, targeting 55,220 GHs able to detect 117 GH families on the 133 available in the CAZy database (www.cazy.org). We included 4848 positive and negative control probes. Non-bacterial families and GH7, 22 and 133, for which thermodynamical parameters did not provide specific probes, were not represented on the CAZyChip.
Regarding the high score of BC_1 and BC_2, we considered our CAZyChip as a promising high-density oligo-DNA microarray, which allows high throughput exploration of bacterial GHs.
Validation of the CAZyChip
The specificity of the CAZychip probes was first evaluated using a set of plasmid bearing GH-encoding sequences, some of which encode well-characterized enzymes [1, 2, 10, 34, 57, 59]. To achieve this, 26 RNA samples from plasmid-bearing bacteria were labeled and hybridized with the probes on the CAZyChip. Figure 1 shows the heatmap (relative signal intensities) for this experiment and illustrates the fact that the vast majority of the samples hybridized quite specifically to the probes on the chip. Pm83 specific probes 2 and 3 not only hybridized with their target RNA, but also to a lesser extent with RNA from Pm85. This cross-hybridization can be easily explained by the fact that both Pm83 and 85 belong to GH8 family and share 81 % nucleotide sequence identity. Regarding probes specific for Pm65 (probe 2), Pm06 (probes 2 and 3), CfXyn (probes 1 and 2), and Pm15 (probe 3), these mostly failed to properly detect their target RNA in the test set (weak signals or no signal). Nevertheless, for each of these targets at least one probe proved to be adequate to properly hybridize to the target RNA and provide unambiguous detection.
To further validate the CAZyChip, RNA from 23 metagenomic clones derived from different gut microbial communities were used ([1, 2, 59, 61]; Table 1). These clones are all characterized by the fact that they bear more than one GH-encoding sequence, with at least one metagenomic clone containing up to 9 GH-encoding sequences (Additional file 3: Table S3). Upon hybridization with the CAZychip, the 23 metagenomic clones resulted in 69 positive signals (Fig. 2), which corresponds to a high detection rate. Most of the GH-encoding sequences were detected by at least one probe, but in some cases by two or three specific probes (Table 1). All genes were expressed in Rum33M21, or Cor367 whereas in Cor28 or Hum5 only a few genes were expressed (sequences GH3- and GH95- from the metagenomic clones Hum5 and Cor28 respectively were not detected), allowing identification of the gene responsible for the activity of each clone (Table 1 and Fig. 2a and b).
Validation of the CAZyChip using individual GH-encoding sequences borne on multi-copy plasmids provided large amounts of RNA that procured strong, saturated hybridization signals for most of the specific probes. However, in the case of fosmid born sequences (metagenomics clones) the intensity of the different hybridization signals was variable, allowing us to determine an accurate minimal detection threshold. This threshold is defined as the minimum signal necessary to differentiate between positive and negative hits in a significant way. As in standard DNA Chip protocols, our samples were labeled with either Cy3 or Cy5. The minimal detection threshold was 8.00 (log 2 of intensity) for Cy3-labelled RNA and 6.70 (log base 2 of intensity) for Cy5-labeled samples. Calculation of the median of variation coefficients (CV) for all experimental probes revealed that this value lies in a narrow range from 1.43 and 4.75 % (Additional file 4: Figure S1), underlining the robustness of the CAZyChip. In addition, 14 GH-encoding sequences cloned either in plasmids (Uhbg_MP, TM1225, XylB, CfXyn and TxXyn) or in fosmids (Cor428 and Hum10), were randomly chosen to be analyzed by qRT-PCR. The results of this analysis were consistent with those obtained using the CAZyChip (Additional file 5: Figure S2).
Exploration of GH diversity evolution in microbial consortium from cow rumen
The CAZyChip was used to investigate the dynamic evolution of stable rumen-derived microbial community displaying good wheat straw degrading ability and a reduced complexity when compared to the parental inoculum . Culture of this stable rumen-derived microbial community presented a 3-phase dynamic behavior over a 15 day period. The initial lag phase was characterized by stable, low-level enzyme activity and very little biomass degradation. The second phase (day 3 to 7), was characterized by an exponential burst of enzyme activities and the third phase was characterized by a stabilized level of enzyme activity . The CAZyChip was used to compare two points that characterize the second phase of the culture, in order to highlight and identify what enzymes are the key players of the wheat straw degradation. The first point corresponded to the beginning of phase 2 (day 3), the second point was in the middle of the phase 2 (day 5), where enzymatic activities were high (Fig. 3).
A limma t test revealed that 2567 GHs were expressed in the two time points: day 3 and day 5 (Additional file 6: Table S4). Both samples displayed a common group of 257 expressed GHs. The two sample points also displayed GH expression unique to the specific time point, with the day 3 sample containing the expression of an additional GH belonging to the GH66 family (accession number AFH61494), and the day 5 sample containing expression of 2309 additional GH’s. Among the total 2566 GHs that were expressed at day 5, only 2 were down-regulated on day 5 compared to day 3 (Additional file 6: Table S4). The weighted differentially expressed genes, and those present at day 5, belong to 96 GH families and are displayed on Fig. 3.
Most of the differentially expressed genes encoding GHs are found in families that are correlated with either cellulose (e.g., GH1, GH3, GH5, and GH8) or hemicellulose (notably heteroxylan) hydrolysis (e.g., GH5, GH10, GH30, GH39, GH43, GH51) (in green Fig. 3b), is consistent with the known chemical composition of wheat straw [21, 35, 54]. CAZyChip analysis also revealed that GH arsenal deployed by the microorganisms in the rumen-derived microbial community contains an extensive range of GH families, including those related to starch hydrolysis (e.g., GH13) and others related to bacterial cell wall degradation (e.g., GH23; Fig. 3b), enzyme activities that are known to be highly represented in all kingdoms.
Using CAZyChip, we are able to explore expression of specific GH families implicated in the targeted functions of plant cell wall polysaccharide degradation. While focusing on GH families involved in enzymatic activities necessary to reach 25 % of wheat straw degradation , we observed an increase of the genes differentially expressed between day 3 and day 5, from GH families containing cellulase, xylanase, exoglucanase and beta-glucosidase activities in accordance with  (Table 2). We observed an enhanced expression of GH1, GH3 and GH5, which according to CAZy, some members of these families are beta-glucosidases and exoglucanases (for GH1 and GH3) or cellulases (for GH5) (Table 2). However, Lazuka et al. have previously shown enhanced cellulase and exoglucanase activities with a constant beta-glucanase activity . Our results strongly suggest that enhanced GH5’s were implicated in efficient cellulase activity and that GH1 and GH3 explained the increased of exoglucanase activity. Our tool allows evaluation of the genetic potential of microbial consortium and highlights complementarity between GHs to contribute to these mechanisms of degradation of plant cell walls.
DNA microarray is one of the most popular technologies for gene expression profiling used in the past 15 years [28, 42, 46, 50, 68]. In this study we presented the development and the validation of the microarray CAZyChip dedicated to analyze the bacterial glycoside hydrolase expression. This is the first high throughput tool, based on DNA microarray technology, allowing the rapid characterization and exploration of the GHs arsenal of complex microbiota at the transcriptomic level. For design purposes, we first collected all sequences of bacterial GHs available in the CAZy data base, belonging to cultivated species, as well as some metagenomic sequences issued from uncultivated species. We then performed a probe bioinformatic design using eArray and ROSO softwares, which took into account the thermodynamics and specificity regardless of the secondary structures that probes can adopt. We validated probe specificity and the robustness of the biochip with different RNAs obtained from well characterized GHs cloned in plasmids and expressed in E. coli. For each GH, we validated at least one specific probe on the three designed per gene. For the great majority of GHs tested, the three probes gave a positive and specific hybridization signal, meaning that our probe design was highly effective.
Following this first validation step with unique GH overexpressed in bacteria, we studied the hybridization behavior of a series of metagenomic clones obtained from different metagenomic libraries. Metagenomic clones were selected for their enzymatic activity and can express up to 9 identified GHs. The CAZyChip allowed for the identification of genes responsible for the activity detected in each metagenomic clone. The multi-genic hybridization step allowed us to validate probes to identify 69 GHs. As an example, His28, which showed arabinofuranosidase activity, encodes two GH51 typical arabinofuranosidases, F but only one was expressed. 96 % of tested GHs had at least one validated probe. Previous studies have demonstrated that the use of multiple probes per target sequence is not essential for in situ synthesized 60mer oligonucleotides in bacterial Agilent’s arrays . Our results demonstrate the robustness of the CAZyChip for GHs detection at transcriptomic level with experimental reproducibility.
Among naturally-occurring biomass-degrading systems, cow rumen represents a natural bioreactor. It is colonized by large communities of symbiotic microorganisms that produce an impressive arsenal of biomass-degrading enzymes, usually including cellulases and hemicellulases. With the CAZyChip, GH expression profiles at two different time points (day 3 and day 5) characterized by an exponential burst of enzyme activities were analyzed. At day 5, we identified overexpression of the GH families associated with cellulase (GH5, GH6, GH8, GH9 and GH48), xylanase (GH8, GH10), and exoglucanase (GH1, GH3) activities, which is in agreement with previous results . The most common activities of GH3 include glucosidases, arabinofuranosidases, xylosidases and glucosaminidases and GH43 shows xylosidase, arabinofuranosidase, arabinanase, xylanase and galactosidase activities. Thus, these two families are implicated in degradation of arabinoxylan, the most abundant hemicellulose component in wheat straw  which explains the great number of genes overexpressed at day 5 in these GH families. An over representation of members of family GH13 and GH23 was seen, as they are implicated in common bacterial physiological processes and known to possess one of the broadest distributions among the gut microbiota [8, 17, 18]. It is the first time that such a generic tool is developed for GH detection from complex microbial ecosystems, although a custom microarray has been previously developed by El Kaoutari et al., to explore partial CAZome of specific human microbiota [17, 18]. This microarray contained probes targeting approximately 7000 genes encoding glycoside hydrolases and selected from 174 reference genomes from specific bacteria present in the human feces.
Our new CAZyChip tool allows the identification of an unprecedented amount of bacterial GHs (55,220) and few CEs, offering opportunities to study expression of a variety of GHs and combinations of enzymes in non-cultivable microorganisms found in any environment. The CAZychip provides an efficient method to explore complex environments, to analyze enriched niches for lignocellulose degradation, and to perform comparative studies. This transcriptomic screening approach (microarrays), reveals the genes that are being actively expressed by lignocellulolytic communities. This in turn allows us to consider the stability and/or performance of target enzymes, enabling the design of new enzyme cocktails and engineering microbial mixed cultures for an optimized lignocellulose bioconversion. Technologies for the rapid screening of GH’s activities are currently in development for high throughput analysis [6, 11, 38, 65]. Functional metagenomic has been proven to be useful tool to achieve this screening of GH’s activities (see review [27, 58]). However, like any screening related technology they face the paradigm “you get what you screened for”. In this context, the CAZyChip allows the observation of the enzymatic arsenal developed by microbial consortia on complex substrate, and could represent a decisive support before choosing a sample for further analysis.
As few CE sequences were included in the CAZyChip, in the near future, others CAZymes (i. e. glycosyltransferases, polysaccharide lyases, or auxiliary activities) could be detectable on the CAZyChip with the same approach. Thanks to its flexible design, this biochip will be able to accommodate additional probes  and could be upgradable, taking into account the regular updates of the CAZy database. Minty et al. have previously proven that fungal-bacterial consortia are efficient for the biosynthesis of valuable products from lignocellulosic feedstocks . Probes for detection of this kind of CAZymes could easily be added on the CAZyChip, in order to highlight a large number of enzymes that work synergistically for cellulose and hemicelluloses breakdown [31, 40]. Understanding the biological process used by bacteria for carbohydrates depolymerization and metabolization is a considerable biotechnological interest not only for biorefineries but also to appreciate carbon flow in the environment, or to promote healthy human nutrition and prevent diseases [17, 18, 40]. The CAZyChip has been developed in a context of lignocellulosic biomass degradation but this biochip represents an excellent tool for other applications in the field of health and nutrition and more widely in any field interested in carbohydrate metabolism. Indeed, GHs are widely characterized in many biological systems such as human intestinal microbiota [17, 18] and the GHs profile are modified depending on eating habits and evolutionary plasticity of the human gut microbiome, playing a major role in nutrition and maintaining human health . Modifications of their expression induce a number of diseases like colon cancer, Crohn’s disease, lactose malabsorption, food allergies, metabolic syndrome, type II diabetes, mucopolysaccharidoses [49, 55, 60, 64]. Thus applications referred for diagnostic or preventive health and nutrition could be explored, if considering GHs as biomarkers. Following glycosyltranferase expression could be also of great interest as they play an important role in the human antigenic system .
In conclusion, the CAZychip developed in this study is a user-friendly, high-throughput, and reliable method to quickly explore GHs expression from complex environmental samples. It can be used to explore functional and ecological dynamics of the enzymatic machinery used by microbes for carbohydrate degradation. This approach can enhance the understanding of how the microbes metabolize polysaccharides and optimize polysaccharide or glycan deconstruction. The CAZyChip could guide the design of enzyme cocktails or the engineering of microbial mixed cultures for many applications.
CE, carbohydrate esterase; DE, differentially expressed; GH, glycoside hydrolase; PCR, polymerase chain reaction; PCW, plant cell wall; qRT-PCR, quantitative real-time PCR; RIN, Ratio Integrity Number
Arnal G, Bastien G, Monties N, Abot A, Anton Leberre V, Bozonnet S, et al. Investigating the function of an arabinan utilization locus isolated from a termite gut community. Appl Environ Microbiol. 2015;81:31–9.
Bastien G, Arnal G, Bozonnet S, Laguerre S, Ferreira F, Fauré R, et al. Mining for hemicellulases in the fungus-growing termite Pseudacanthotermes militaris using functional metagenomics. Biotechnol Biofuels. 2013;6:78.
Bayer EA, Belaich J-P, Shoham Y, Lamed R. The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Annu Rev Microbiol. 2004;58:521–54.
Beg QK, Kapoor M, Mahajan L, Hoondal GS. Microbial xylanases and their industrial applications: a review. Appl Microbiol Biotechnol. 2001;56:326–38.
Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics Oxf Engl. 2003;19:185–93.
Boutard M, Cerisy T, Nogue P-Y, Alberti A, Weissenbach J, Salanoubat M, et al. Functional diversity of carbohydrate-active enzymes enabling a bacterium to ferment plant biomass. PLoS Genet. 2014;10:e1004773.
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.
Cantarel BL, Lombard V, Henrissat B. Complex carbohydrate utilization by the healthy human microbiome. PLoS One. 2012;7:e28742.
Carpita NC. Progress in the biological synthesis of the plant cell wall: new ideas for improving biomass for bioenergy. Curr Opin Biotechnol. 2012;23:330–7.
Cecchini DA, Laville E, Laguerre S, Robe P, Leclerc M, Doré J, et al. Functional metagenomics reveals novel pathways of prebiotic breakdown by human gut bacteria. PLoS One. 2013;8:e72766.
Chauvigné-Hines LM, Anderson LN, Weaver HM, Brown JN, Koech PK, Nicora CD, et al. Suite of activity-based probes for cellulose-degrading enzymes. J Am Chem Soc. 2012;134:20521–32.
Chen X, Luo Y, Yu H, Sun Y, Wu H, Song S, et al. Transcriptional profiling of biomass degradation-related genes during Trichoderma reesei growth on different carbon sources. J Biotechnol. 2014;173:59–64.
Chundawat SPS, Beckham GT, Himmel ME, Dale BE. Deconstruction of lignocellulosic biomass to fuels and chemicals. Annu Rev Chem Biomol Eng. 2011;2:121–45.
Cosgrove DJ. Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005;6:850–61.
D’Elia JN, Salyers AA. Effect of regulatory protein levels on utilization of starch by Bacteroides thetaiotaomicron. J Bacteriol. 1996;178:7180–6.
Dugat-Bony E, Peyretaillade E, Parisot N, Biderre-Petit C, Jaziri F, Hill D, et al. Detecting unknown sequences with DNA microarrays: explorative probe design strategies. Environ Microbiol. 2012;14:356–71.
El Kaoutari A, Armougom F, Gordon JI, Raoult D, Henrissat B. The abundance and variety of carbohydrate-active enzymes in the human gut microbiota. Nat Rev Microbiol. 2013;11:497–504.
El Kaoutari A, Armougom F, Leroy Q, Vialettes B, Million M, Raoult D, et al. Development and validation of a microarray for the investigation of the CAZymes encoded by the human gut microbiome. PLoS One. 2013;8:e84033.
Ferraresso S, Vitulo N, Mininni AN, Romualdi C, Cardazzo B, Negrisolo E, et al. Development and validation of a gene expression oligo microarray for the gilthead sea bream (Sparus aurata). BMC Genomics. 2008;9:580.
Fontes CMGA, Gilbert HJ. Cellulosomes: highly efficient nanomachines designed to deconstruct plant cell wall complex carbohydrates. Annu Rev Biochem. 2010;79:655–81.
Gilbert HJ. The biochemistry and structural biology of plant cell wall deconstruction. Plant Physiol. 2010;153:444–55.
Häkkinen M, Valkonen MJ, Westerholm-Parvinen A, Aro N, Arvas M, Vitikainen M, et al. Screening of candidate regulators for cellulase and hemicellulase production in Trichoderma reesei and identification of a factor essential for cellulase production. Biotechnol Biofuels. 2014;7:14.
Hasunuma T, Okazaki F, Okai N, Hara KY, Ishii J, Kondo A. A review of enzymes and microbes for lignocellulosic biorefinery and the possibility of their application to consolidated bioprocessing technology. Bioresour Technol. 2013;135:513–22.
Hehemann J-H, Kelly AG, Pudlo NA, Martens EC, Boraston AB. Bacteria of the human gut microbiome catabolize red seaweed glycans with carbohydrate-active enzyme updates from extrinsic microbes. Proc Natl Acad Sci U S A. 2012;109:19786–91.
Henrissat B. A classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem J. 1991;280(Pt 2):309–16.
Henrissat B, Romeu A. Families, superfamilies and subfamilies of glycosyl hydrolases. Biochem J. 1995;311(Pt 1):350–1.
Heux S, Meynial-Salles I, O’Donohue MJ, Dumon C. White biotechnology: state of the art strategies for the development of biocatalysts for biorefining. Biotechnol Adv. 2015;33:1653–70.
He Z, Deng Y, Van Nostrand JD, Tu Q, Xu M, Hemme CL, et al. GeoChip 3.0 as a high-throughput tool for analyzing microbial community composition, structure and functional activity. ISME J. 2010;4:1167–79.
He Z, Gentry TJ, Schadt CW, Wu L, Liebich J, Chong SC, et al. GeoChip: a comprehensive microarray for investigating biogeochemical, ecological and environmental processes. ISME J. 2007;1:67–77.
Himmel ME, Bayer EA. Lignocellulose conversion to biofuels: current challenges, global perspectives. Curr Opin Biotechnol. 2009;20:316–7.
Himmel ME, Ding S-Y, Johnson DK, Adney WS, Nimlos MR, Brady JW, et al. Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science. 2007;315:804–7.
Hughes TR, Mao M, Jones AR, Burchard J, Marton MJ, Shannon KW, et al. Expression profiling using microarrays fabricated by an ink-jet oligonucleotide synthesizer. Nat Biotechnol. 2001;19:342–7.
Kirk O, Borchert TV, Fuglsang CC. Industrial enzyme applications. Curr Opin Biotechnol. 2002;13:345–51.
Ladevèze S, Tarquis L, Cecchini DA, Bercovici J, André I, Topham CM, et al. Role of glycoside phosphorylases in mannose foraging by human gut bacteria. J Biol Chem. 2013;288:32370–83.
Lagaert S, Pollet A, Courtin CM, Volckaert G. β-xylosidases and α-L-arabinofuranosidases: accessory enzymes for arabinoxylan degradation. Biotechnol Adv. 2014;32:316–32.
Lazuka A, Auer L, Bozonnet S, Morgavi DP, O’Donohue M, Hernandez-Raquet G. Efficient anaerobic transformation of raw wheat straw by a robust cow rumen-derived microbial consortium. Bioresour Technol. 2015;196:241–9.
Leiske DL, Karimpour-Fard A, Hume PS, Fairbanks BD, Gill RT. A comparison of alternative 60-mer probe designs in an in-situ synthesized oligonucleotide microarray. BMC Genomics. 2006;7:72.
Liu G, Qin Y, Li Z, Qu Y. Development of highly efficient, low-cost lignocellulolytic enzyme systems in the post-genomic era. Biotechnol Adv. 2013;31:962–75.
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.
Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS. Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev MMBR. 2002;66:506–77. table of contents.
Martens EC, Koropatkin NM, Smith TJ, Gordon JI. Complex glycan catabolism by the human gut microbiota: the Bacteroidetes Sus-like paradigm. J Biol Chem. 2009;284:24673–7.
Maruyama K, Yamaguchi-Shinozaki K, Shinozaki K. Gene expression profiling using DNA microarrays. Methods Mol Biol Clifton NJ. 2014;1062:381–91.
Minty JJ, Singer ME, Scholz SA, Bae C-H, Ahn J-H, Foster CE, et al. Design and characterization of synthetic fungal-bacterial consortia for direct production of isobutanol from cellulosic biomass. Proc Natl Acad Sci. 2013;110:14592–7.
Musso G, Gambino R, Cassader M. Interactions between gut microbiota and host metabolism predisposing to obesity and diabetes. Annu Rev Med. 2011;62:361–80.
Park BH, Karpinets TV, Syed MH, Leuze MR, Uberbacher EC. CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. Glycobiology. 2010;20:1574–84.
Patro JN, Ramachandran P, Lewis JL, Mammel MK, Barnaba T, Pfeiler EA, et al. Development and utility of the FDA ‘GutProbe’ DNA microarray for identification, genotyping and metagenomic analysis of commercially available probiotics. J Appl Microbiol. 2015;118:1478–88.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Purushe J, Fouts DE, Morrison M, White BA, Mackie RI, North American Consortium for Rumen Bacteria, et al. Comparative genome analysis of Prevotella ruminicola and Prevotella bryantii: insights into their environmental niche. Microb Ecol. 2010;60:721–9.
Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490:55–60.
Rasooly A, Herold KE. Food microbial pathogen detection and analysis using DNA microarray technologies. Foodborne Pathog Dis. 2008;5:531–50.
Reeves AR, Wang GR, Salyers AA. Characterization of four outer membrane proteins that play a role in utilization of starch by Bacteroides thetaiotaomicron. J Bacteriol. 1997;179:643–9.
Reymond N, Charles H, Duret L, Calevro F, Beslon G, Fayard J-M. ROSO: optimizing oligonucleotide probes for microarrays. Bioinformatics Oxf Engl. 2004;20:271–3.
Rogowski A, Briggs JA, Mortimer JC, Tryfona T, Terrapon N, Lowe EC, Baslé A, Morland C, Day AM, Zheng H, Rogers TE, Thompson P, Hawkins AR, Yadav MP, Henrissat B, Martens EC, Dupree P, Gilbert HJ, Bolam DN. Glycan complexity dictates microbial resource allocation in the large intestine. Nat Commun. 2015;6:7481. doi:10.1038/ncomms8481.
Scheller HV, Ulvskov P. Hemicelluloses. Annu Rev Plant Biol. 2010;61:263–89.
Sheng YH, Hasnain SZ, Florin THJ, McGuckin MA. Mucins in inflammatory bowel diseases and colorectal cancer. J Gastroenterol Hepatol. 2012;27:28–38.
Smyth GK, Michaud J, Scott HS. Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics Oxf Engl. 2005;21:2067–75.
Song L, Siguier B, Dumon C, Bozonnet S, O’Donohue MJ. Engineering better biomass-degrading ability into a GH11 xylanase using a directed evolution strategy. Biotechnol Biofuels. 2012;5:3.
Steele HL, Jaeger K-E, Daniel R, Streit WR. Advances in recovery of novel biocatalysts from metagenomes. J Mol Microbiol Biotechnol. 2009;16:25–37.
Tasse L, Bercovici J, Pizzut-Serin S, Robe P, Tap J, Klopp C, et al. Functional metagenomics to mine the human gut microbiome for dietary fiber catabolic enzymes. Genome Res. 2010;20:1605–12.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.
Ufarté L, Laville É, Duquesne S, Potocki-Veronese G. Metagenomics for the discovery of pollutant degrading enzymes. Biotechnol Adv. 2015;33:1845–54.
Van Dyk JS, Pletschke BI. A review of lignocellulose bioconversion using enzymatic hydrolysis and synergistic cooperation between enzymes–factors affecting enzymes, conversion and synergy. Biotechnol Adv. 2012;30:1458–80.
Vasconcelos-Dos-Santos A, Oliveira IA, Lucena MC, Mantuano NR, Whelan SA, Dias WB, et al. Biosynthetic machinery involved in aberrant glycosylation: promising targets for developing of drugs against cancer. Front Oncol. 2015;5:138.
Veneault-Fourrey C, Commun C, Kohler A, Morin E, Balestrini R, Plett J, et al. Genomic and transcriptomic analysis of Laccaria bicolor CAZome reveals insights into polysaccharides remodelling during symbiosis establishment. Fungal Genet Biol. 2014;72:168–81.
Vidal-Melgosa S, Pedersen HL, Schückel J, Arnal G, Dumon C, Amby DB, et al. A new versatile microarray-based method for high throughput screening of carbohydrate-active enzymes. J Biol Chem. 2015;290:9020–36.
Vincentelli R, Cimino A, Geerlof A, Kubo A, Satou Y, Cambillau C. High-throughput protein expression screening and purification in Escherichia coli. Methods San Diego Calif. 2011;55:65–72.
Wu L, Thompson DK, Li G, Hurt RA, Tiedje JM, Zhou J. Development and evaluation of functional gene arrays for detection of selected genes in the environment. Appl Environ Microbiol. 2001;67:5780–90.
Zhou A, He Z, Qin Y, Lu Z, Deng Y, Tu Q, et al. StressChip as a high-throughput tool for assessing microbial community responses to environmental stresses. Environ Sci Technol. 2013;47:9841–9.
Acknowledgements and funding
This work was supported by the IDEX Transversality program from Toulouse University (grant number: 2014–628) and by the French National Institute for Agricultural Research (INRA, ‘Meta-omics of Microbial Ecosystems’ research program, and AIP Biotechnology 2011). A. Abot, was funded by a grant from the Conseil Régional Midi-Pyrénées (grant number 12053333). We thank S. Bozonnet for helpful discussion, our colleagues L Ufarté, P. Alvira, B. Guyez, M. Abadie, S. Ladevèze, and S. Comtet-Marre (INRA, UR454 Microbiology, Centre Auvergne-Rhône-Alpes, Saint Genes Champanelle; Université d’Auvergne, EA 4678 CIDAM, Clermont-Ferrand; Lallemand Animal Nutrition, Blagnac, France) for technical assistance, and fruitful discussion. The GeT-TQ and GeT-Biochip Platform from Toulouse Genopole are acknowledged for skillful technical assistance. Roxanne Diaz for her careful rereading of English. BH wishes to thank funding from Programmes d’Avenir BIP: BIP and AMIDEX (programme Microbio-E).
Availability of data and materials
All relevant data are available within the manuscript and its additional files.
AA (post-doctoral position) employed on this project, did all microarrays experiments and interpretations. She carried out the molecular genetic studies, participated in the sequence alignment and drafted the manuscript. GA did all the gene cloning and prepared the termite and rumen libraries. LA and AL did the chemical analyses, enzyme activity assays and characterization of wheat straw degradation by consortium. DL and SL are engineers that performed the statistical analysis, oligonucleotide design, sequence alignments and helped AA to deposit the data in the Gene Expression Omnibus (GEO) public database. LT is the engineer who trained AB to prepare all the RNA extraction and microarrays experiments. EL and GPV did the high-throughput screening and sequencing of human metagenomic libraires and gave final approval of the version to be published, VL and BH made the Cazymes classification within CAZy and extracted oligonucleotides sequences for AA and DL. MO’D and GHR revised the manuscript critically for important intellectual content. GHR is the supervisor of LA and AL. CD and VAL are the project coordinators, conceived the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Management of experimental animals followed the guidelines for animal research of the French Ministry of Agriculture and other applicable guidelines and regulations for animal experimentation in the European Union, in particular the European Directive 2010/63/EU on the protection of animals used for scientific purposes. The experimental farm has a license for keeping rumen-cannulated cattle in their premises. Sampling of rumen contents is considered a minor manipulation that is neither painful nor stressful to the animal. The sampling procedure was performed only once and a specific examination by the regional ethical committee was not considered necessary.
Primers sequences used for real-time qPCR quantification. (XLSX 12 kb)
Summary of probes sequences on the CAZyChip, BC_scores and the accession number of associated gene. (XLSX 8992 kb)
Nucleic sequences of fosmids from metagenomic study. (XLSX 57 kb)
Boxplots of coefficient of variation for specific probes of targeted GHs cloned in fosmids. (JPG 90 kb)
mRNA levels from (A) GHs cloned in plasmid (B-C) or in fosmids were quantified by real-time qPCR and normalized to rpoS mRNA levels. (JPG 43 kb)
GHs differentially expressed in enriched microbial consortium from cattle rumen microbiome between day 3 and day 5. (XLSX 349 kb)
About this article
- CAZymes detection
- Glycoside hydrolase
- Microbial functional diversity
- Plant cell wall degradation
- Transcriptomic analysis