Original Research ARTICLE
Temporal Expression Dynamics of Plant Biomass-Degrading Enzymes by a Synthetic Bacterial Consortium Growing on Sugarcane Bagasse
- 1Microbial Ecology Cluster, Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, Netherlands
- 2Department of Biological Sciences, Universidad de los Andes, Bogotá, Colombia
Plant biomass (PB) is an important source of sugars useful for biofuel production, whose degradation efficiency depends on synergistic and dynamic interactions of different enzymes. Here, using a metatranscriptomics-based approach, we explored the expression of PB-degrading enzymes in a five-species synthetic bacterial consortium during cultivation on sugarcane bagasse as a unique carbon source. By analyzing the temporal expression dynamics of a selection of enzymes we revealed the functional role of each consortium member and disentangled the potential interactions between them. Based on normalized expression values and the taxonomic affiliation of all the transcripts within thirty carbohydrate-active enzyme (CAZy) families, we observed a successional profile. For instance, endo-glucanases/-xylanases (e.g., GH8, GH10, and GH16) were significantly expressed at 12 h, whereas exo-glucanases (e.g., GH6 and GH48) and α-arabinosidases/β-xylosidases (e.g., GH43) were highly expressed at 48 h. Indeed, a significant peak of extracellular β-xylosidase activity was observed at this stage. Moreover, we observed a higher expression of several CAZy families at 12–48 h, suggesting easy access to the main plant polysaccharides. Based on this evidence, we predicted that the highest level of collaboration between strains takes place at the initial stages of growth. Here, Paenibacillus, Brevundimonas, and Chryseobacterium were the most important contributors, whereas Stenotrophomonas was highly active at the end of the culture (96–192 h) without contributing to a large extent to the expression of lignocellulolytic enzymes. Our results contribute to the understanding of enzymatic and ecological mechanisms within PB-degrading microbial consortia, yielding new perspectives to improve the PB saccharification processes.
Sugarcane bagasse (SCB) has been considered a useful bioenergy feedstock due to its low cost and huge availability (Cardona et al., 2010). As the majority of agricultural residues, SCB is composed of about 40% cellulose (β-1,4-linked D-glucose units), 25% hemicellulose, 20% lignin and a small percentage of pectin (Kim and Day, 2011; Szczerbowski et al., 2014). In agricultural residues, three main polysaccharides constitute the hemicellulosic fraction, i.e., xylan/arabinoxylan, xyloglucan, and galacto(gluco)mannan. These polymers are usually classified according to the sugar residues present in the backbone. For instance, xylan/arabinoxylan is composed of β-1,4-linked D-xylose units, which may be substituted by different side groups (e.g., D-galactose, L-arabinose, and glucuronic acid) (van den Brink and de Vries, 2011; de Souza, 2013). Thus, due to the high SCB complexity, an efficient release of the monosaccharides requires the synergic interaction of several enzymes (e.g., endo/exo -glucanases, endo-xylanases, β-glucosidases, α-arabinosidases, lytic polysaccharide monooxygenases, β-xylosidases, and α-glucuronidases) (Gao et al., 2011; Rytioja et al., 2014; Kim et al., 2017). In the production of biofuels, an efficient SCB saccharification process still represents a bottleneck due to the recalcitrant nature of the most plant biomass (PB) (Himmel et al., 2007; Canilha et al., 2012; Visser et al., 2015).
Recently, constructed PB-degrading microbial consortia have proven to be a excellent sources of (hemi)cellulolytic enzyme cocktails (Park et al., 2012; Jiménez et al., 2015b). A common strategy to obtain these types of consortia is by the dilution-to-stimulation method (Lee et al., 2013). In this approach, a microbial community (inoculum) is set to grow on agricultural residues in liquid batch cultures that are sequentially diluted across several transfers (Zhang et al., 2013; de Lima Brossi et al., 2015; Jiménez et al., 2016). Throughout the enrichment, the inoculum gradually changes in composition, yielding a specialized microbial consortium with progressively reduced richness. In this process, two types of successional dynamics may occur: the first one is the succession along the sequential transfers – the selection phase – and the other one occurs within the relative ‘stable’ consortia, where microbial species with complementary activities synergistically degrade the PB substrate (Jiménez et al., 2017). These PB-degrading microbial consortia are usually composed of hundreds of microbial species, from which small core sets of organisms – i.e., the minimal effective consortium – contribute to most of the degradation process whereas the remaining species could act mostly as ‘cheaters,’ consuming the sugars released by the effective core. For instance, in our previous studies, Sphingobacterium and Klebsiella species were the most relevant degraders in microbial consortia cultivated on wheat straw (Jiménez et al., 2014a, 2015a,b; Maruthamuthu et al., 2016). These types of microbial enrichments have served as sources for the isolation of efficient PB-degrading strains, which can be further used either as source of new enzymes or as active strains useful for the design of synthetic microbial communities (Evans et al., 2017).
The efficiency of PB-degrading microbial consortia relies probably on the level of complementarity and synergism that is achieved, either in terms of enzymatic pool (or the timing of enzymatic release) and/or number of active species. The temporal aspects of microbial PB degradation are of high relevance, because enzyme-enzyme synergism can occur in a temporal scale, as is the case for the conversion of cellulose to cellobiose and the subsequent release of glucose (Lynd et al., 2002). Two recent studies explored the temporal expression dynamics of lignocellulolytic enzymes in PB-degrading microbial consortia. In these, the expression profile of some glycoside hydrolases (GHs) was assessed by metaproteomics (Zhu et al., 2016) as well as through the hybridization of total consortium RNA against a new CAZychip (Abot et al., 2016). For instance, Zhu et al. (2016) have evaluated the secreted proteins in a corn stover-adapted consortium along 1, 3, and 7 days of growth. They detected the CAZy families GH1, GH3, GH9, GH10, GH11, GH13, GH43 and GH94 after 24 h of cultivation. However, an in-depth understanding of the enzymatic mechanisms and ecological interactions in such consortia is still missing. The modulation of the species diversity/composition in a synthetic microbial consortium provides an experimental tool to unravel the enzymatic mechanisms that can potentially lead to a significant improvement in PB saccharification processes (Lindemann et al., 2016; Cavaliere et al., 2017; Jiménez et al., 2017). In addition, here we posit that to understand interactions, it is crucial to disentangle the temporally explicit expression dynamics of lignocellulolytic enzymes in these microbial systems.
In this study, we analyzed the temporal expression dynamics of thirty PB-degrading enzymes in a five-species synthetic bacterial consortium – previously selected for complementarity and high degradation potential – along a single culture batch on SCB. A successional expression profile was observed to occur in the consortium, where Paenibacillus, Brevundimonas, and Chryseobacterium were the most relevant contributors to the expression of lignocellulolytic enzymes. Ecological interactions and successional dynamics between the consortium members were discussed based on a recently developed conceptual framework (Jiménez et al., 2017). Our metatranscriptomics analysis portrays the temporal enzymatic processes underlying a synthetic PB-degrading microbial consortium, improving our ecological understanding on these systems and giving new perspectives to improve PB saccharification processes.
Materials and Methods
Growth of the Synthetic Bacterial Consortium on Different Substrates
Five phylogenetically different bacterial strains isolated from a soil-derived lignocellulolytic microbial consortium were selected. This selection was based on parameters established by Puentes-Téllez and Salles (2018), after screening of 60 synthetic five-species consortia varying in their functional metabolic diversity and drawn from a pool of 18 bacterial species. The selected strains (mixture 48) when grown together on PB (after 96 hrs of cultivation), could degrade up to 50% of the total lignocellulose fraction, with specific degradation percentages of 43.01, 39.65, 69.18% of lignin, cellulose and hemi-cellulose, respectively (Puentes-Téllez and Salles, 2018). Based on previous 16S rRNA gene sequencing, these strains were identified, at the genus level, as Stenotrophomonas, Paenibacillus, Microbacterium, Chryseobacterium, and Brevundimonas.
In order to analyze the temporal growth dynamics of the synthetic bacterial consortium, we cultivated the bacterial mixture (five strains) on four different substrates as sole carbon source varying in their complexity (glucose -GLC, carboxymethylcellulose -CMC, xylan -XYL and sugarcane bagasse -SCB). Reagents like GLC, CMC and XYL were obtained from Sigma–Aldrich. The SCB was previously milled into small pieces (≤1 mm), washed twice with 70% of ethanol and distilled water, and then dried at 55°C for 48 h. Firstly, each bacterial strain was cultivated in 5 ml of LB broth for 48 h. The cultures were centrifuged (6,500 g × 10 min) and the cells were washed with 1 ml of 0.9% NaCl. Absorbance (OD600 nm) for all the strains was adjusted with 0.9% NaCl at 1.0 in order to ensure similar starting densities. Subsequently, 25 μl of each strain (relation 1:1) were inoculated into 25 ml of sterile mineral salt medium containing 0.5% of each substrate, vitamins and trace element solutions (Jiménez et al., 2014b). Triplicate flasks per time point were incubated at 28°C in aerobic conditions (180 rpm) and samples were taken at 12, 48, 96, and 192 h for subsequent analysis. Microbial growth was evaluated by absorbance (OD600 nm) for each time point. Absorbance in the SCB-derived cultures was measured directly at the upper layer, after substrate sedimentation. Two controls, i.e., one without substrate and the other without microbial source were set up (Supplementary Figure S1).
Quantification of Sugars, Proteins and Enzymatic Activities in the Extracellular Fractions from the Synthetic Bacterial Consortium Cultivated on Xylan and Sugarcane Bagasse
Sugar and protein concentrations, as well as the enzymatic activities of the consortium cultivated on XYL and SCB were measured in the extracellular fraction. To obtain these fractions, 10 ml of each culture at each time point were centrifuged (6,500 g × 8 min). The supernatant was clarified by passing through a 0.22 μm filter. Quantity of sugars and proteins was obtained by 3,5-dinitrosalicilic acid (DNS) method (Miller, 1959) and the Quick StartTM Bradford Protein Assay (Biorad, Hercules, CA, United States), respectively. Enzymatic activities related with PB degradation were evaluated using pNP-labeled substrates, i.e., p-nitrophenyl α-D-glucopyranoside (pNPαGlu), p-nitrophenyl β-D-xylopyranoside (pNPβXyl), p-nitrophenyl β-D-galactopyranoside (pNPβGal) and p-nitrophenyl α-D-mannopyranoside (pNPαMan). The reaction mixtures consisted of 200 μl of 7.5 mM of each p-nitrophenol-glycoside (diluted in 50 mM of Tris-HCL pH 7.5) and 150 μl of each supernatant. The mixtures were incubated at 37°C for 4 h, after which the reactions were stopped on ice. Two negative controls were used for all assays: (i) reaction mixture without pNP-substrate; (ii) reaction mixture using sterile water. Enzymatic activities were determined from the measured absorbance units using a standard calibration curve. Activity was detected by the presence of yellow color and the amount of para-nitrophenol liberated was measured by absorbance at 410 nm (Maruthamuthu et al., 2017). All the analyses were performed using three biological replicates and standard deviations (σ) were obtained.
Extraction of Total RNA and Metatranscriptomic Sequencing
Total microbial RNA extraction from the consortium was carried out along the batch of culture on XYL and SCB. For each time point (12, 48, 96, and 192 h) three biological replicates were taken. Briefly, 10 ml of RNAprotect Bacteria Reagent (Qiagen, Hilden, Germany) were added to 5 ml of microbial culture. After 5 min (at room temperature) the mixture was centrifuged at 5,000 g for 10 min. The cell pellets were lysed using proteinase K (50 mg/ml) and lysozyme (15 mg/ml), and the RNA extraction was performed using the RNeasy Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The quality and quantity of the extracted total RNA was checked by non-denaturing agarose gel electrophoresis and absorbance ratios 260/280 nm. Depletion of ribosomal RNA (using Ribo-ZeroTM Kit), mRNA libraries preparation and sequencing (Illumina NextSeq 2 × 150 bp; 1 lane) were performed at LGC Genomics (Berlin, Germany).
Analysis of Metatranscriptomic Data
A total of 116,205,704 quality filtered paired-end reads (>20 bp) from 21 samples were used as input for de novo transcript assembly using Trinity (Haas et al., 2013). Transcript abundance estimation and normalization was performed using scripts included in the Trinity package. The script abundance_estimates_to_matrix.pl was used to create a matrix with the estimated RNA-Seq fragment raw counts that were used for differential expression analysis with edgeR (Robinson et al., 2010). To enable direct in-between sample comparison of transcript abundances, raw abundances were converted to TPM (transcripts per million transcripts) and TMM (trimmed mean of M-values) normalized in RSEM (RNA-Seq by Expectation Maximization) (Dillies et al., 2012). The script analyze_diff_expr.pl was used to extract those transcripts that were at least fourfold differentially expressed at a FDR (false discovery rate) corrected P-value cut off of 0.001 in any of the pairwise sample comparisons, followed by hierarchical clustering based on the Pearson correlation matrix or pairwise comparisons of samples using edgeR (Robinson et al., 2010). Moreover, total transcript contigs were converted to potential coding regions using the TransDecoder function developed within Trinity. For functional classification, the proteins were annotated using the GhostKOALA (Kanehisa et al., 2016) and dbCAN web-platforms (Yin et al., 2012). In order to unveil the relative expression dynamics of lignocellulolytic proteins, and in accordance with previous omics studies (Simmons et al., 2014; Zhou et al., 2014; Berlemont and Martiny, 2015; Jiménez et al., 2015a; 2016; Abot et al., 2016; Antunes et al., 2016; Zhu et al., 2016), we selected thirty carbohydrate-active enzyme (CAZy) families (Lombard et al., 2014) that are highly relevant in the deconstruction of PB. All transcripts affiliated within those CAZy families were manually annotated by BLASTp against the NCBI-nr database. The taxonomic affiliation of each transcript was performed by manually taking the most frequent BLASTp hit between the top-50. Finally, the TMM-normalized TPM values were added to each transcript classification and summed up to obtain the values per each CAZy family (Pfreundt et al., 2016). Relative expression of the CAZy families that segregated significantly between the time points were identified using random forest analysis with 1000 trees followed by the Boruta algorithm for feature selection (average z-scores of 1000 runs > 3) (Cardenas et al., 2015). All metatranscriptomic data are publically accessible on the MG-RAST (Meyer et al., 2008) server under the IDs mgm4734165.3 to mgm4734173.3; and mgm4735448.3 to mgm4735459.3.
Growth of the Synthetic Bacterial Consortium on Different Substrates
The growth of the consortium was evaluated on four carbon sources with different structural complexity. The results showed that the consortium is able to grow using GLC, XYL, and SCB, but not CMC as a sole source of carbon (Figure 1A). In the GLC culture, a maximal growth was observed at 48 h, gradually decreasing toward the end of the batch (192 h). Regarding the XYL culture, we observed the exponential growth phase after 12 h, followed by the start of the stationary phase at 96 h. In the case of SCB, the exponential phase started at an early growth stage (between 0 to 12 h), and apparently was maintained until 96 h of the culture. Subsequent analyses were performed on the cultures with bacterial growth in substrates with higher complexity (XYL and SCB cultures).
FIGURE 1. (A) Microbial growth (OD600 nm) of the synthetic bacterial consortium on glucose –GLC (), carboxymethylcellulose –CMC (), xylan -XYL () and sugarcane bagasse –SCB (). (B) Total amount of reducing sugars (mg/ml) produced by the consortium on XYL and SCB along a batch of culture (12–192 h). Relative enzymatic activities (β-xylosidase: pNPβXyl; α-glucosidase: ; pNPαGlu; β-galactosidase: pNPβGal; and α-mannosidases: pNPαMan) of the consortium along a batch of culture on (C) XYL and (D) SCB. All the analyses were performed using three biological replicates and the error bars correspond to standard deviations.
Quantification of Extracellular Sugars and Enzymatic Activities
In the XYL and SCB cultures, the concentration of the total proteins secreted by the consortium was less than 0.1 mg/ml. Regarding release of reducing sugars, maximum values (0.30–0.34 mg/ml) for the XYL culture were found between 48 and 192 h, during the exponential and stationary phases. In the SCB, we observed an increase of sugars during the entire growth period, with a maximum value at the end of the culture (0.28 ± 0.01 mg/ml; mean ± SD) (Figure 1B). Along the incubation time, we evaluated four types of enzymatic activities (target enzymes that can act on terminal ends of substrates) related with hemicellulose (β-xylosidase, β-galactosidases, α-mannosidase) and pectin (α-glucosidase) degradation processes. Thus, the consortium cultivated on XYL displayed maximal β-xylosidase and α-glucosidase activity values at 48 and 192 h, which was consistent with the sugar accumulation data (Figure 1C). Concerning the consortium cultivated on SCB, we observed a high β-xylosidase, α-glucosidase and β-galactosidases activity at 48 h versus a low α-mannosidase activity. The activities were relatively stable until the last measurement, except for β-xylosidase where a significant peak (p < 0,001, t-test of pairwise comparisons) was observed at 48 h (Figure 1D).
Metatranscriptomic Overview of the Synthetic Bacterial Consortium Cultivated on Xylan and Sugarcane Bagasse
Total RNA extraction from the consortium cultivated on XYL and SCB was performed on triplicate systems along four time points of incubation. However, it was not possible to isolate sufficient RNA at the initial stage (12 h) of the XYL culture due to the low microbial activity and the presence of undigested xylan. Thus, a total of 21 samples were subjected to mRNA sequencing (9 from XYL and 12 from SCB cultures). Approximately 5.7 and 7.1 Gb of total mRNA sequencing data were obtained for the XYL and SCB samples, respectively. The approximately 116 million reads were assembled into 59,734 contigs. Given that N50 values discard read coverage information, we report the E90N50 value (1.1 kb) as a more useful indicator of transcriptome assembly quality than the N50 (665 bp) (Haas et al., 2013). Also note that the E90 number of transcripts (23,364) for which the E90N50 value is computed is roughly a third of the total number of transcripts assembled (68,738) and for which the N50 statistic was calculated. Based on the RNA-Seq fragment raw counts, we observed that 4,213 transcripts were differentially expressed in all pairwise combinations. As expected, the clustering of these transcripts showed that the XYL and SCB cultures had different expression profiles. For each, the profiles clustered by triplicates divided by treatment. Only the SCB time points 96 and 192 h formed an intermingled cluster based on the differential expression profile (Supplementary Figure S2), suggesting that after 96 h of growth on SCB the consortial expression patterns were ‘stable.’
Relative Expression Profile of Plant Biomass-Degrading Enzymes in the Synthetic Bacterial Consortium Cultivated on Sugarcane Bagasse
Thirty CAZy families were selected in order to evaluate the expression of transcripts related with PB-degrading enzymes in the consortium. Here, we focus our analyses on the SCB cultures. A heat map using TMM-normalized TPM values showed a trend in the expression profiles of enzymes from the consortium along the incubation time (Figure 2 and Supplementary Table S1). In general terms, most of the selected CAZy families were highly expressed at 12 and 48 h compared with 96–192 h. For instance, xyloglucanases/endo-glucanases (GH16), endo-xylanases (families GH8 and GH10) and α-glucuronidases (GH67) were significantly expressed at 12 hrs. Moreover, lytic polysaccharide monooxygenases -LPMOs (AA10), exo-glucanases (GH6 and GH48), α-arabinosidase/β-xylosidase (GH43), next to other endo-glucanases (GH5 and GH9) and endo-xylanases (GH11) were highly expressed at 48 h compared with all other time points. In addition, enzyme families involved in pectin degradation (e.g., GH28 and PL10) and β-glucosidases (specifically GH3) were mostly expressed at 48 hrs. Regarding lignin degradation, the family AA2 (catalases-peroxidases) showed higher values of expression at 12 (117 ± 42) and 192 hrs (118 ± 26). In terms of relative expression values (i.e., TMM-normalized TPM) between the CAZy families in each time point, we observed that transcripts of enzymes of families GH16 (1115 ± 447) and GH10 (1005 ± 207) were highly expressed at 12 h (Figure 3A). Moreover, those from families GH11 (439 ± 42), GH10 (339 ± 2.5), GH43 (414 ± 44), GH3 (236 ± 31) and GH9 (163 ± 20) were the most expressed at 48 hrs. At the final stages of the cultures (96–192 h), the relative expression values of enzymes belonging to families GH43, GH13, GH11, GH3 and CE3 were higher in comparison with those of the other families.
FIGURE 2. Heat map displaying the relative expression values (row z-score of TMM-normalized TPM) of thirty CAZy families (mostly involved in plant biomass degradation) in the synthetic bacterial consortium along a batch of culture (12–192 h) on sugarcane bagasse (SCB). Asterisks represent families differentially segregated along the incubation time that were identified by Boruta featured selection.
Taxonomic Affiliation of Expressed Transcripts Related with Thirty Plant Biomass-Degrading Enzymes in the Synthetic Bacterial Consortium
The relative expression of the rpoA gene (housekeeping single-copy gene throughout bacteria) was used as a marker to evaluate the metabolic activity of each member the consortium along the incubation time on SCB, thus providing a depiction of the relative contribution of each strain in the CAZy expression profile described above. The taxonomic affiliation and the TMM-normalized TPM values of the rpoA transcripts showed that Brevundimonas and Paenibacillus were the most active members at 12 and 48 h, respectively. In contrast, Chryseobacterium and Stenotrophomonas were highly active at 96–192 h compared with the other three consortium members. Microbacterium showed low activity along the entire incubation time (Figure 3B). In addition, we characterized all transcripts encoding enzymes of the thirty CAZy families to species (Figure 4). This analysis indicated that Brevundimonas and Paenibacillus were major contributors to the expression of enzymes involved in cellulose degradation. Specifically, Brevundimonas was relevant in the expression of β-glucosidases/galactosidases (GH1 and GH3) at 12 and 48 h (transcript IDs 14 and 15). The expression of enzymes of CAZy families GH10, GH11, GH30, GH43, GH51, and GH67, Paenibacillus and Chryseobacterium were the most important drivers of the degradation of xylan/arabinoxylan along incubation time (transcript IDs 17–26). In addition, endo-xylanases (GH8), endo-xyloglucanases (GH74), β-galactosidases (GH2), xyloglucanases/mannosidases (GH44) and β-endomannases (GH26) were mainly expressed by Paenibacillus at 48 h. Pectin degradation was mainly associated with the expression of enzymes of families GH28 and PL10 from Brevundimonas at 48 h (transcript IDs 47 and 50). The expression of catalases-peroxidases (AA2) was carried out by Paenibacillus (at 12 h), Chryseobacterium and Brevundimonas (at 192 h). Taxonomic affiliation, TMM-normalized TPM values (average and standard deviations) and the functional annotation of the 50 most important transcripts within the 30 CAZy families can be found in the Supplementary Table S2.
FIGURE 3. (A) Relative expression values (TMM-normalized TPM) of thirty CAZy families in the synthetic bacterial consortium, cultivated on sugarcane bagasse (SCB), within each incubation time. Error bars correspond to standard deviations. (B) Relative expression values and taxonomic affiliation of rpoA gene transcripts (RNA polymerase alpha subunit; KEGG orthology number K03040) within the consortium members along a batch of culture. In time zero, we assumed a low metabolic activity and probably undetected expression of rpoA gene.
FIGURE 4. Taxonomic affiliation and TMM-normalized TPM values (average) of individual transcripts within thirty CAZy families. Axis X shows the incubation time of the synthetic microbial consortium cultivated on sugarcane bagasse (SCB). In time zero, we assumed a low metabolic activity and probably undetected expression of these CAZy families. Average transcription values lines affiliated to Paenibacillus, Brevundimonas, Chryseobacterium, Stenotrophomonas, and Microbacterium are colored in red, green, blue, purple and black, respectively. Numbers beside the highest expression lines within each panel correspond to the transcript ID in the Supplementary Table S2.
In this study, we applied a metatranscriptomics-based approach to unveil the temporal expression dynamics of thirty PB-degrading enzymes (or CAZy families) in a synthetic consortium composed of five bacterial strains that was previously selected for its high degradation capacity. Overall, the data suggest that a specific successional expression profile of lignocellulolytic enzymes occurred in the consortium cultivated on SCB. For instance, the expression of specific endo-glucanases and endo-xylanases becomes highly relevant at 12–48 h, whereas enzymes acting on the external side linkages of the hemicellulose (e.g., α-arabinosidases/β-xylosidases) and cellulose (exo-glucanases and β-glucosidases) fractions are relatively highly expressed only after 48 h of cultivation (Figures 2, 5). It is important to highlight that all predicted functions for the CAZy families were based on sequence annotation and that we observed high variability in the expression of some CAZy families, indicating that the different flasks (biological replicates) might show slight variation in the relative abundance of each strain. Nevertheless, the observed profile fits a conceptual model that was recently proposed (Jiménez et al., 2017) and suggests the successive stages of PB deconstruction. The main findings of this study are discussed based on three stages: 12, 48, and 96–192 h.
FIGURE 5. Cartoon of the successional profile representing the expression of lignocellulolytic enzymes by the synthetic bacterial consortium, along a batch of culture on sugarcane bagasse (SCB), in relation to the degradation process in the (hemi)cellulosic fraction of the plant biomass. The model shows how the selected CAZy families can degrade hemicellulose based on the reported evidence and also the taxonomic affiliation to genus level of the proteins expressed by members of the consortium. The size of the enzymes (semi-circles shape) is relative proportional to the relative expression values. Enzyme shapes colored in red, green, blue and purple correspond to enzymes expressed mostly by Paenibacillus, Brevundimonas, Chryseobacterium, and Stenotrophomonas, respectively. Specific molecules conforming each polymer follow the shape and color code shown in the top right corner.
First Stage (12 h)
Out of the thirty CAZy families evaluated, a total of 23 were highly expressed at this stage (row Z-score > 0 in at least one biological replicate). By using a random forest analysis and the Boruta statistical test, eleven CAZy families were significantly more expressed at 12 h (Figure 2), suggesting that this is the most functionally diverse phase compared with the other time points. In this stage, transcripts encoding enzymes of CAZy families GH8, GH10 and GH16 were significantly expressed compared with other time points. The family GH16 contains enzymes that can be assigned to five subgroups according to their substrate specificities, including xyloglucan transglucosylases/hydrolases (XTHs), (1,3)-β-galactanases, (1,4)-β-galactanases/κ-carrageenases, “non-specific” (1,3/1,3; 1,4)-β-D-glucan endohydrolases and (1,3; 1,4)-β-D-glucan endohydrolases (Strohmeier et al., 2004). Specifically, in our study, the high expression of this family at 12 h was due to one transcript that was affiliated to a 1,3-1,4-endoglucanase from Paenibacillus (78% identity) (Supplementary Table S2). The relevance of GH16 proteins has been acknowledged in previous studies. For instance, proteins from this family have been found in high relative abundances in the metagenomes of three aerobic PB-degrading microbial consortia (Jiménez et al., 2016). In addition, Zhou et al. (2014) reported that one of the microbial members (Clostridium -FC2) of the cellulolytic consortium F1RT has the capacity to produce cellulosomal enzymes belonging to family GH16. In Trichoderma reesei and T. harzianum these types of enzymes have been found highly expressed (or induced) in presence of cellulosic substrates (e.g., cellulose, avicel and sophorose) (Häkkinen et al., 2012; Ferreira Filho et al., 2017). On the other hand, proteins from GH16 have a putative structural and evolutionary connection with endoxylanases from family GH11 (Strohmeier et al., 2004). Based on the lack of experimental support on the cellulolytic activity of GH16 transcripts, we can speculate that, in our system, these enzymes may have an important role in transformation of hemicellulose (e.g., xyloglucan) present in the SCB. Based on the ability of our consortium to grow on GLC, XYL and SCB, but not on CMC, we suggest that the expression of enzymes involved in cellulose degradation (e.g., endo-glucanases) could be regulated by the presence natural cellulose of other type of PB polysaccharides present in the SCB. CMC could be used as an inductor of endoglucanases (Xiao et al., 2005). However, this polymer has a different structure compared to the natural cellulose found in PB.
In our study, two transcripts (IDs 17–18) contributed to a large extent to the high expression level of enzymes from family GH10 at 12 h (Figure 4). These were affiliated to endo-1,4-β-xylanases from Paenibacillus (>96% identity). Proteins from the family GH10 and GH16 were also strongly expressed and secreted in a corn stover-adapted microbial consortium at 24 and 72 h, respectively (Zhu et al., 2016). The high expression level of endo-xylanases in the consortium at 12 h suggests that xylan/arabinoxylan is attacked at early growth stages. One of the most important features of this stage was the high metabolic activity of Brevundimonas (Figure 3B). This organism may participate in the expression of endo-glucanases and endo-xyloglucanases from families GH5 and GH74, respectively. It is key to note that proteins from family GH5 can breakdown different plant polymers [e.g., cellulose, xyloglucan, and galacto(gluco)mannan]. Proteins derived from organisms affiliated with Brevundimonas species have been found in high relative abundance in two consortia cultivated with once-used PB (Jiménez et al., 2016). These results suggest that Brevundimonas may play an important role in xylan, cellulose and lignin degradation. Moreover, the expression of α-glucuronidases – which cleave 4-O-methyl-α-D-glucuronic acid (MeGlcA) from xylose residues when they are attached to the terminal, non-reducing end of xylo-oligosaccharides (Dodd and Cann, 2009) – by Chryseobacterium and Paenibacillus was significantly higher at 12 h, indicating that the consortium can carry out this process at the early stages of cultivation. Interestingly, the removal of MeGlcA could expose the xylan to attack by endo-xylanases and β-xylosidases (Rogowski et al., 2014). Two transcripts (ID 41 and 44) affiliated to AA2 and AA6 (and derived from Paenibacillus and Brevundimonas, respectively) were comparatively highly expressed at 12 h (Figure 4). We speculate that the expression of these transcripts is probably related to cell protection from oxidative damage and the transfer of electrons.
At this stage, the low cell densities and the presence of accessible PB polysaccharides could enhance the fast growth of Paenibacillus and Brevundimonas. Here, competition and niche partitioning events are likely to occur (Jiménez et al., 2017). Thus, based on the high expression of Brevundimonas-derived transcripts affiliated to α-fucosidases (GH29) and β-galactosidases (GH35) at this stage (Figure 4), we suggest that this species may consume preferably fucose and galactose, whereas Paenibacillus could have more access to arabinose and mannose. Competition for the most abundant sugars (glucose and xylose) could be a valid explanation for the lower activity of Brevundimonas in the next stages. At an early phase, Microbacterium could consume easily available oligo- or monosaccharides present in the SCB. However, we suggest that the low metabolic activity of Microbacterium could be directly related with an intrinsic characteristic of the organism (e.g., oligotrophy), niche competition and possibly a low metabolic versatility. Probably, in the consortium, Microbacterium could confer functional stability, redundancy and robustness, thus increasing the metabolic diversity and providing a potential buffer against external disturbance (Awasthi et al., 2014). In addition, Microbacterium species could secrete molecules that synergistically stimulate the degradation capacity of other PB-degrading bacterial species.
Second Stage (48 h)
This is the key stage for the degradation of the main PB polysaccharides. Indeed, we observed raised extracellular enzymatic activities at 48 h, with a significant peak of β-xylosidase (Figure 1D). Here, efficient hydrolysis of cellulose can be carried out by a synergistic action of non-processive endo-glucanases (GH5 and GH9), which produce new ends at random within the polysaccharide chain, and processive exo-glucanases (GH6 and GH48), which remain attached to the substrate and split off cellobiose from such free ends. Then, cellobiose can be converted to glucose by the action of β-glucosidases expressed by Brevundimonas. An interesting synergistic interaction is carried out at this stage, where consumption of cellobiose by Brevundimonas potentially prevents metabolite repression of endo/exo-glucanases produced by Paenibacillus (Lynd et al., 2002; Berlemont and Martiny, 2013). Moreover, LPMOs from Paenibacillus may play a relevant role in the cellulose degradation at this stage. These LPMOs are metalloenzymes (copper-dependent) that catalyze the oxidative cleavage of (1,4)-linked glycosidic bonds of crystalline PB polysaccharide surfaces (Horn et al., 2012; Agger et al., 2014).
Concerning hemicellulose degradation, we found that transcripts affiliated to Paenibacillus from the families GH11 (endo-xylanases) and GH43 (α-arabinosidases/β-xylosidases) were highly expressed when compared to the other time points. We found a shift in the expression of endo-xylanases – i.e., from GH10 at 12 h to GH11 at 48 h – which display an endo-action and thus hydrolyze β-(1,4)-linkages between adjacent xylose residues within the xylan/arabinoxylan chain. The enzymes of these families essentially differ in their structure and substrate specificity. For instance, GH10 xylanases have a (β/α)8-fold and are much more versatile. In contrast, GH11 enzymes are the smallest xylanases with a β-jelly-roll architecture, are highly substrate-specific and do not tolerate high substitutions on the backbone (Biely et al., 1997; Paës et al., 2012). This suggests that a high level of enzymatic complementarity can occur in Paenibacillus within the consortium along the growth period. Regarding GH43, it can display dual action (Maruthamuthu et al., 2017), thus removing arabinose from the external-side chains attached to the xylan backbone or hydrolyze xylobiose into its monomeric units, releasing D-xylose (Rytioja et al., 2014). The family GH43 enzymes seem to play a key role in lignocellulose degradation, being highly enriched in two soil-derived microbial consortia cultivated on wheat straw (Jiménez et al., 2015a) and highly expressed and secreted between 72 to 168 h in PB-degrading microbial consortia cultivated under aerobic and anaerobic conditions (Abot et al., 2016; Zhu et al., 2016).
Although Paenibacillus was the most metabolically active organism and possibly the biggest contributor to the (hemi)cellulose deconstruction, Brevundimonas may play a key complementary role through the expression of pectinolytic enzymes (e.g., GH28 and PL10). Pectins comprise a diverse and complex group of polymers that function as structural components in PB (Mohnen, 2008). Thus, the action of pectinolityc enzymes could increase the sugar release from the (hemi)cellulosic fraction (Xiao and Anderson, 2013). Theoretically, different consortium members can contribute in distinctive ways and proportions to the SCB degradation process. However, we suggest that organisms that contribute more to the substrate transformation may profit from the most abundant sugars (glucose or xylose) (Jiménez et al., 2017). Based on this conceptualization, at this stage, Paenibacillus might indeed be the most abundant organism in terms of cell density. In other words, Paenibacillus invests a lot in the production of lignocellulolytic enzymes, and it might grow very well at the expense of the large proportion of sugars.
Third Stage (96–192 h)
At this final stage, the expression of the majority of the CAZy families decreased. A similar trend occurs with the metabolic activity of Paenibacillus and Brevundimonas. These two organisms may be acting as biggest contributors to the degradation process at early stages (12 and 48 h) of the culture, At this stage, Paenibacillus probably entered stationary phase, which is reflected in the decrease of GHs gene expression. We suggest that at these time points the major fraction of the available (hemi)cellulose is already degraded. Here, not only carbon sources but also other essential nutrients, such as nitrogen, can become progressively limiting. Although production of sugars (by PB polysaccharide degradation) and their consumption occurs simultaneously in the consortium, we observed an increase in the concentration of reducing sugars along the time of incubation (Figure 1B). This suggests that, at the end of the culture, the environment becomes relatively more homogeneous. Thus, the increasing release of easily assimilable carbon sources is expected to enhance competition among the bacterial strains. Stenotrophomonas and Chryseobacterium showed higher normalized expression values of their rpoA gene at the end of the culture (192 h) (Figure 3B), suggesting that they become highly metabolically active when the concentration of monosaccharides is higher. We considered Stenotrophomonas to be a ‘sugar cheater’ in the consortium. That means that it profits from the sugars made available by the other community members, without effectively contributing to the PB degradation process (Jiménez et al., 2017). This is the case of Thermus thermophilus, which lacks the metabolic potential for PB polysaccharide hydrolysis, in a switchgrass-degrading consortium (D’haeseleer et al., 2013). Moreover, the exact functional roles of Microbacterium and Stenotrophomonas in the consortia remained unknown and further studies need to be performed in order to understand their impact in this system. However, we speculate that they could be involved in feedback regulation, cross-feeding, maintenance of functional stability and inhibition of metabolite repression. Finally, at this stage, lignin bioconversion was thought to take place by the expression of catalases-peroxidases (CAZy family AA2) from Chryseobacterium and Brevundimonas. Similar to patterns observed in composting (Antunes et al., 2016), where high expression of ligninolytic enzymes occurs at the end of the process.
Conclusion and Final Remarks
In this study, we provide evidence for a temporally explicit dynamics of expression of PB-degrading enzymes in a five-species synthetic bacterial consortium growing on a complex substrate such as SCB. Initially, Paenibacillus and Brevundimonas were found to become active, presumably synergistically, in the degradation of cellulose and hemicellulose (between 12 to 48 h), whereas Chryseobacterium and Stenotrophomonas were most active at the end of the culture. Based on our results, we suggest that the consortium may follow an a ‘endogenous heterotrophic succession’ trajectory (Fierer et al., 2010; Jiménez et al., 2017), which is characterized by a rapid initial activity of so-called pioneer populations that initiate the deconstruction of PB (in the consortium, Paenibacillus and Brevundimonas), followed by the relative decline of these initial populations. Other species that outcompete the pioneer populations in sugar uptake then start to proliferate (in the consortium, Chryseobacterium and Stenotrophomonas). Thus, we posit that a shift from synergism to competition is an important trait in these systems. This dichotomy could be correlated with the complexity of the available carbon sources (Deng and Wang, 2016).
Understanding the ecological interactions and the temporal expression dynamics of PB-degrading enzymes in a synthetic microbial consortium can be used to improve aerobic saccharification processes for biofuel production. Thus, we propose a dynamic saccharification process where two types of enzyme cocktails could be applied in a sequential manner. The first cocktail could contain enzymes acting mainly in the backbone of hemicellulose (e.g., xyloglucanases/endo-glucanases from CAZy family GH16; and endo-xylanases from families GH8 and GH10). The second cocktail could be composed by a diverse array of endo-enzymes, exo-enzymes and LPMOs where proteins from CAZy families GH1, GH3 (β-glucosidases/galactosidases), GH5, GH9 (endo-glucanases), GH6, GH48 (exo-glucanases), GH11, GH30 (endo-xylanases) and GH43 (α-arabinosidase/β-xylosidases) need to be considered as essential components. Although lignin is an important polymer of PB, it is still unclear which type of lignolytic enzymes could increase the release of sugars from the (hemi)cellulosic fraction. This last subject will require more attention in the future.
As a cautionary note, mRNA abundance is not necessarily indicative of protein levels and enzymatic activities, and so it is possible that actual activities last longer than the time of (instantaneous) mRNA detection. One may argue here that the successional expression dynamics of PB-degrading enzymes on the basis of mRNA provides an overview of the inferred activity dynamics, and so still improves the ecological understanding of PB-degrading microbial consortia. Metaproteomics approaches as well as the relative abundances (e.g., cell population densities and/or 16S rRNA gene copy numbers) of each organism – growing alone and in a consortium – could improve our understanding of the successional dynamics in these systems. In addition, the growth of the individual strains on different carbon sources would unveil substrate preferences, interactions, and niche partitioning and competition events. Moreover, we suggest that further enzymatic tests, substrate composition and chromatographic analysis could improve the evidence that the bacterial consortium starts breaking the internal bonds of the plant polysaccharides, moving to the external branches as succession takes place.
DJJ performed most of the experiments, participated in the bioinformatic analyses and drafted the manuscript. MCDM performed all the bioinformatic procedures. JFS conceived the study, participated in its design, coordination, and contributed to the drafting of the manuscript.
This work was financed by the Microwaste Project (No. 729004006; FAPESP-NWO).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Pilar Puentes-Tellez for providing the bacterial strains; to Francisco Dini-Andreote and Maria Julia de Lima Brossi for their help in the statistical analysis. In addition, to Jan Dirk van Elsas for his comments and support.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00299/full#supplementary-material
FIGURE S1 | Flasks containing the synthetic bacterial consortium at 96 h of growth in (A) glucose -GLC, (B) carboxymethylcellulose -CMC, (C) xylan -XYL and (D) sugarcane bagasse –SCB. Flasks marked with WI are without inoculum.
FIGURE S2 | Heat map showing the expression values (Log2_TMM) of all significantly differentially expressed transcripts of the synthetic bacterial consortium across all samples along the batch of culture in xylan –XYL (48, 96, and 192 h) and sugarcane bagasse –SCB (12, 48, 96, and 192 h).
TABLE S1 | Relative expression values (TMM-normalized TPM) used to construct the heat map in the Figure 2.
TABLE S2 | Taxonomic affiliation, TMM-normalized TPM values (average and standard deviations-SD) and the functional annotation (by BLASTp) of the 50 most important transcripts, pointed out in the Figure 4, within the thirty CAZy families.
Abot, A., Arnal, G., Auer, L., Lazuka, A., Labourdette, D., Lamarre, S., et al. (2016). CAZyChip: dynamic assessment of exploration of glycoside hydrolases in microbial ecosystems. BMC Genomics 17:671. doi: 10.1186/s12864-016-2988-4
Agger, J. W., Isaksen, T., Várnai, A., Vidal-Melgosa, S., Willats, W. G., Ludwig, R., et al. (2014). Discovery of LPMO activity on hemicelluloses shows the importance of oxidative processes in plant cell wall degradation. Proc. Natl. Acad. Sci. U.S.A. 111, 6287–6292. doi: 10.1073/pnas.1323629111
Antunes, L. P., Martins, L. F., Pereira, R. V., Thomas, A. M., Barbosa, D., Lemos, L. N., et al. (2016). Microbial community structure and dynamics in thermophilic composting viewed through metagenomics and metatranscriptomics. Sci. Rep. 6:38915. doi: 10.1038/srep38915
Awasthi, A., Singh, M., Soni, S. K., Singh, R., and Kalra, A. (2014). Biodiversity acts as insurance of productivity of bacterial communities under abiotic perturbations. ISME J. 8, 2445–2452. doi: 10.1038/ismej.2014.91
Canilha, L., Kumar Chandel, A., dos Santos, Milessi, T. S., Fernandes Antunes, F. A., da Costa, et al. (2012). Bioconversion of sugarcane biomass into ethanol: an overview about composition, pretreatment methods, detoxification of hydrolysates, enzymatic saccharification, and ethanol fermentation. J. Biomed. Biotechnol. 2012:989572. doi: 10.1155/2012/989572
Cardenas, E., Kranabetter, J. M., Hope, G., Maas, K. R., Hallam, S., and Mohn, W. W. (2015). Forest harvesting reduces the soil metagenomic potential for biomass decomposition. ISME J. 9, 2465–2476. doi: 10.1038/ismej.2015.57
Cardona, C. A., Quintero, J. A., and Paz, I. C. (2010). Production of bioethanol from sugarcane bagasse: status and perspectives. Bioresour. Technol. 101, 4754–4766. doi: 10.1016/j.biortech.2009.10.097
Cavaliere, M., Feng, S., Soyer, O., and Jimenez, J. I. (2017). Cooperation in microbial communities and their biotechnological applications. Environ. Microbiol. 19, 2949–2963. doi: 10.1111/1462-2920.13767
de Lima Brossi, M. J., Jiménez, D. J., Cortes-Tolalpa, L., and van Elsas, J. D. (2015). Soil-derived microbial consortia enriched with different plant biomass reveal distinct players acting in lignocellulose degradation. Microb. Ecol. 71, 616–627. doi: 10.1007/s00248-015-0683-7
de Souza, W. R. (2013). “Microbial degradation of lignocellulosic biomass,” in Sustainable Degradation of Lignocellulosic Biomass -Techniques, Applications and Commercialization, eds A. Chandel and S. Da Silva (Rijeka: InTech), 208–209.
D’haeseleer, P., Gladden, J. M., Allgaier, M., Chain, P. S., Tringe, S. G., Malfatti, S. A., et al. (2013). Proteogenomic analysis of a thermophilic bacterial consortium adapted to deconstruct switchgrass. PLoS One 8:e68465. doi: 10.1371/journal.pone.0068465
Dillies, M. A., Rau, A., Aubert, J., Hennequet-Antier, C., Jeanmougin, M., Servant, N., et al. (2012). A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. Brief. Bioinform. 14, 671–683. doi: 10.1093/bib/bbs046
Evans, R., Alessi, A. M., Bird, S., McQueen-Mason, S. J., Bruce, N. C., and Brockhurst, M. A. (2017). Defining the functional traits that drive bacterial decomposer community productivity. ISME J. 11, 1680–1687. doi: 10.1038/ismej.2017.22
Ferreira Filho, J. A., Horta, M. A. C., Beloti, L. L., Dos Santos, C. A., and de Souza, A. P. (2017). Carbohydrate-active enzymes in Trichoderma harzianum: a bioinformatic analysis bioprospecting for key enzymes for the biofuels industry. BMC Genomics 18:779. doi: 10.1186/s12864-017-4181-9
Fierer, N., Nemergut, D., Knight, R., and Craine, J. M. (2010). Changes through time: integrating microorganisms into the study of succession. Res. Microbiol. 161, 635–642. doi: 10.1016/j.resmic.2010.06.002
Gao, D., Uppugundla, N., Chundawat, S. P., Yu, X., Hermanson, S., Gowda, K., et al. (2011). Hemicellulases and auxiliary enzymes for improved conversion of lignocellulosic biomass to monosaccharides. Biotechnol. Biofuels 4:5. doi: 10.1186/1754-6834-4-5
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Häkkinen, M., Arvas, M., Oja, M., Aro, N., Penttilä, M., Saloheimo, M., et al. (2012). Re-annotation of the CAZy genes of Trichoderma reesei and transcription in the presence of lignocellulosic substrates. Microb. Cell Fact. 11:134. doi: 10.1186/1475-2859-11-134
Himmel, M. E., Ding, S. Y., Johnson, D. K., Adney, W. S., Nimlos, M. R., Brady, J. W., et al. (2007). Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science 315, 804–807. doi: 10.1126/science.1137016
Jiménez, D. J., Chaves-Moreno, D., and van Elsas, J. D. (2015a). Unveiling the metabolic potential of two soil-derived microbial consortia selected on wheat straw. Sci. Rep. 5:13845. doi: 10.1038/srep13845
Jiménez, D. J., de Lima Brossi, M. J., Schückel, J., Kračun, S. K., Willats, W. G., and van Elsas, J. D. (2016). Characterization of three plant biomass-degrading microbial consortia by metagenomics- and metasecretomics-based approaches. Appl. Microbiol. Biotechnol. 100, 10463–10477. doi: 10.1007/s00253-016-7713-3
Jiménez, D. J., Dini-Andreote, F., DeAngelis, M. K., Singer, S. W., Salles, J. F., and van Elsas, J. D. (2017). Ecological insights in the dynamics of plant biomass-degrading microbial consortia. Trends Microbiol. 25, 788–796. doi: 10.1016/j.tim.2017.05.012
Jiménez, D. J., Dini-Andreote, F., and van Elsas, J. D. (2014a). Metataxonomic profiling and prediction of functional behaviour of wheat straw degrading microbial consortia. Biotechnol. Biofuels 7:92. doi: 10.1186/1754-6834-7-92
Jiménez, D. J., Korenblum, E., and van Elsas, J. D. (2014b). Novel multispecies microbial consortia involved in lignocellulose and 5-hydroxymethylfurfural bioconversion. Appl. Microbiol. Biotechnol. 98, 2789–2803. doi: 10.1007/s00253-013-5253-7
Jiménez, D. J., Maruthamuthu, M., and van Elsas, J. D. (2015b). Metasecretome analysis of a lignocellulolytic microbial consortium grown on wheat straw, xylan and xylose. Biotechnol. Biofuels 8:199. doi: 10.1186/s13068-015-0387-8
Kanehisa, M., Sato, Y., and Morishima, K. (2016). BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol. 428, 726–731. doi: 10.1016/j.jmb.2015.11.006
Kim, I. J., Lee, H. J., and Kim, K. H. (2017). Pure enzyme cocktails tailored for the saccharification of sugarcane bagasse pretreated by using different methods. Process Biochem. 57, 167–174. doi: 10.1016/j.procbio.2017.04.006
Kim, M., and Day, D. F. (2011). Composition of sugar cane, energy cane, and sweet sorghum suitable for ethanol production at Louisiana sugar mills. J. Ind. Microbiol. Biotechnol. 38, 803–807. doi: 10.1007/s10295-010-0812-8
Lee, D. J., Show, K. Y., and Wang, A. (2013). Unconventional approaches to isolation and enrichment of functional microbial consortium-a review. Bioresour. Technol. 136, 697–706. doi: 10.1016/j.biortech.2013.02.075
Lindemann, S. R., Bernstein, H. C., Song, H. S., Fredrickson, J. K., Fields, M. W., Shou, W., et al. (2016). Engineering microbial consortia for controllable outputs. ISME J. 10, 2077–2084. doi: 10.1038/ismej.2016.26
Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42, D490–D495. doi: 10.1093/nar/gkt1178
Lynd, L. R., Weimer, P. J., van Zyl, W. H., and Pretorius, I. S. (2002). Microbial cellulose utilization: fundamentals and biotechnology. Microbiol. Mol. Biol. Rev. 66, 506–577. doi: 10.1128/MMBR.66.3.506-577.2002
Maruthamuthu, M., Jiménez, D. J., Stevens, P., and van Elsas, J. D. (2016). A multi-substrate approach for functional metagenomics-based screening for (hemi)cellulases in two wheat straw-degrading microbial consortia unveils novel thermoalkaliphilic enzymes. BMC Genomics 17:86. doi: 10.1186/s12864-016-2404-0
Maruthamuthu, M., Jiménez, D. J., and van Elsas, J. D. (2017). Characterization of a furan aldehyde-tolerant β-xylosidase/α-arabinosidase obtained through a synthetic metagenomics approach. J. Appl. Microbiol. doi: 10.1111/jam.13484 [Epub ahead of print].
Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386
Park, J. I., Steen, E. J., Burd, H., Evans, S. S., Redding-Johnson, A. M., Batth, T., et al. (2012). A thermophilic ionic liquid-tolerant cellulase cocktail for the production of cellulosic biofuels. PLoS One 7:e37010. doi: 10.1371/journal.pone.0037010
Pfreundt, U., Spungin, D., Bonnet, S., Berman-Frank, I., and Hess, W. R. (2016). Global analysis of gene expression dynamics within the marine microbial community during the VAHINE mesocosm experiment in the southwest Pacific. Biogeosciences 13, 4135–4149. doi: 10.5194/bg-13-4135-2016
Puentes-Téllez, P., and Salles, J. F. (2018). Construction of effective minimal active microbial consortia for lignocellulose degradation. Microb. Ecol. doi: 10.1007/s00248-017-1141-5 [Epub ahead of print].
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Rogowski, A., Baslé, A., Farinas, C. S., Solovyova, A., Mortimer, J. C., Dupree, P., et al. (2014). Evidence that GH115 α-glucuronidase activity, which is required to degrade plant biomass, is dependent on conformational flexibility. J. Biol. Chem. 289, 53–64. doi: 10.1074/jbc.M113.525295
Rytioja, J., Hildén, K., Yuzon, J., Hatakka, A., de Vries, R. P., and Mäkelä, M. R. (2014). Plant-polysaccharide-degrading enzymes from Basidiomycetes. Microbiol. Mol. Biol. Rev. 78, 614–649. doi: 10.1128/MMBR.00035-14
Simmons, C. W., Reddy, A. P., D’haeseleer, P., Khudyakov, J., Billis, K., Pati, A., et al. (2014). Metatranscriptomic analysis of lignocellulolytic microbial communities involved in high-solids decomposition of rice straw. Biotechnol. Biofuels 7:495. doi: 10.1186/s13068-014-0180-0
Strohmeier, M., Hrmova, M., Fischer, M., Harvey, A. J., Fincher, G. B., and Pleiss, J. (2004). Molecular modeling of family GH16 glycoside hydrolases: potential roles for xyloglucan transglucosylases/hydrolases in cell wall modification in the poaceae. Protein Sci. 13, 3200–3213. doi: 10.1110/ps.04828404
Szczerbowski, D., Pitarelo, A. P., Zandoná Filho, A., and Ramos, L. P. (2014). Sugarcane biomass for biorefineries: comparative composition of carbohydrate and non-carbohydrate components of bagasse and straw. Carbohydr. Polym. 114, 95–101. doi: 10.1016/j.carbpol.2014.07.052
Visser, E. M., Leal, T. F., de Almeida, M. N., and Guimarães, V. M. (2015). Increased enzymatic hydrolysis of sugarcane bagasse from enzyme recycling. Biotechnol. Biofuels 8:5. doi: 10.1186/s13068-014-0185-8
Zhang, Q., Tian, M., Tang, L., Li, H., Li, W., Zhang, J., et al. (2013). Exploration of the key microbes involved in the cellulolytic activity of a microbial consortium by serial dilution. Bioresour. Technol. 132, 395–400. doi: 10.1016/j.biortech.2012.11.097
Zhou, Y., Pope, P. B., Li, S., Wen, B., Tan, F., Cheng, S., et al. (2014). Omics-based interpretation of synergism in a soil-derived cellulose-degrading microbial community. Sci. Rep. 4:5288. doi: 10.1038/srep05288
Zhu, N., Yang, J., Ji, L., Liu, J., Yang, Y., and Yuan, H. (2016). Metagenomic and metaproteomic analyses of a corn stover-adapted microbial consortium EMSD5 reveal its taxonomic and enzymatic basis for degrading lignocellulose. Biotechnol. Biofuels 9:243. doi: 10.1186/s13068-016-0658-z
Keywords: bacterial consortium, expression dynamics, lignocellulolytic enzymes, metatranscriptomics, sugarcane bagasse, synergism
Citation: Jiménez DJ, Chaib De Mares M and Salles JF (2018) Temporal Expression Dynamics of Plant Biomass-Degrading Enzymes by a Synthetic Bacterial Consortium Growing on Sugarcane Bagasse. Front. Microbiol. 9:299. doi: 10.3389/fmicb.2018.00299
Received: 21 September 2017; Accepted: 08 February 2018;
Published: 26 February 2018.
Edited by:Suhelen Egan, University of New South Wales, Australia
Reviewed by:Johan Larsbrink, Chalmers University of Technology, Sweden
François Thomas, UMR8227 Laboratoire de Biologie Intégrative des Modèles Marins, France
Guillermina Hernandez-Raquet, Institut National de la Recherche Agronomique (INRA), France
Fabio M. Squina, Universidade de Sorocaba, Brazil
Copyright © 2018 Jiménez, Chaib De Mares and Salles. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Diego Javier Jiménez, firstname.lastname@example.org