Leandro N. Lemos1,2
Roberta V. Pereira1
Ronaldo B. Quaggio1Layla F. Martins1Livia M. S. Moura1,2Amanda R. da Silva1,2Luciana P. Antunes1
Aline M. da Silva1*†
João C. Setubal1,3*†- 1Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
- 2Programa de Pós-Graduação Interunidades em Bioinformática, Universidade de São Paulo, São Paulo, Brazil
- 3Biocomplexity Institute, Virginia Tech, Blacksburg, VA, USA
Microbial consortia selected from complex lignocellulolytic microbial communities are promising alternatives to deconstruct plant waste, since synergistic action of different enzymes is required for full degradation of plant biomass in biorefining applications. Culture enrichment also facilitates the study of interactions among consortium members, and can be a good source of novel microbial species. Here, we used a sample from a plant waste composting operation in the São Paulo Zoo (Brazil) as inoculum to obtain a thermophilic aerobic consortium enriched through multiple passages at 60°C in carboxymethylcellulose as sole carbon source. The microbial community composition of this consortium was investigated by shotgun metagenomics and genome-centric analysis. Six near-complete (over 90%) genomes were reconstructed. Similarity and phylogenetic analyses show that four of these six genomes are novel, with the following hypothesized identifications: a new Thermobacillus species; the first Bacillus thermozeamaize genome (for which currently only 16S sequences are available) or else the first representative of a new family in the Bacillales order; the first representative of a new genus in the Paenibacillaceae family; and the first representative of a new deep-branching family in the Clostridia class. The reconstructed genomes from known species were identified as Geobacillus thermoglucosidasius and Caldibacillus debilis. The metabolic potential of these recovered genomes based on COG and CAZy analyses show that these genomes encode several glycoside hydrolases (GHs) as well as other genes related to lignocellulose breakdown. The new Thermobacillus species stands out for being the richest in diversity and abundance of GHs, possessing the greatest potential for biomass degradation among the six recovered genomes. We also investigated the presence and activity of the organisms corresponding to these genomes in the composting operation from which the consortium was built, using compost metagenome and metatranscriptome datasets generated in a previous study. We obtained strong evidence that five of the six recovered genomes are indeed present and active in that composting process. We have thus discovered three (perhaps four) new thermophillic bacterial species that add to the increasing repertoire of known lignocellulose degraders, whose biotechnological potential can now be investigated in further studies.
Introduction
Plant biomass can be decomposed by complex lignocellulolytic microbial communities present in natural environments, such as forest soil (Eichorst and Kuske, 2012) and cow rumen (Hess et al., 2011), or in engineered ecosystems, such as composting (Neher et al., 2013) or biogas fermenters (Güllert et al., 2016). For biomass degradation in these environments lignocellulolytic fungal and bacterial species employ hydrolytic and oxidative enzymes, which act synergistically to depolymerize cellulose, hemicellulose, and lignin (Allgaier et al., 2010; Koeck et al., 2014; Hemsworth et al., 2015; López-Mondéjar et al., 2016). These studies have thus shown that taxonomically diverse members within a lignocellulolytic microbial community work in cooperation to fully deconstruct plant biomass.
Microbial consortia selected from complex lignocellulolytic microbial communities are promising alternatives to deconstruct plant waste rather than microbial monocultures (Zuroff and Curtis, 2012; Peng et al., 2016). Indeed, culture enrichment provides a good strategy for studying the interactions among their microbial members as well as their lignocellulolytic enzyme systems (Wongwilaiwalin et al., 2010; D'Haeseleer et al., 2013; Kinet et al., 2015; de Lima Brossi et al., 2016; Zhu et al., 2016). Moreover, natural microbial consortia can be a method of culture enrichment for unculturable microbial species (Vartoukian et al., 2010; D'Haeseleer et al., 2013).
Several studies have succeeded in establishing lignocellulolytic consortia from diverse sources such as soil (Feng et al., 2011; Gao et al., 2014; Jiménez et al., 2014; de Lima Brossi et al., 2016; Cortes-Tolalpa et al., 2016), biogas-producing digester (Yan et al., 2012) and composting (Haruta et al., 2002; Wongwilaiwalin et al., 2010; Gladden et al., 2011; Eichorst et al., 2013; Kinet et al., 2015; Wang et al., 2016; Zhu et al., 2016). Even though all these consortia exhibit lignocellulose degradation capabilities, their microbial community structure, and composition are distinct. The composition of the final microbial consortia appears to be strongly driven by the substrate used as carbon source (Eichorst et al., 2014; Simmons et al., 2014b; de Lima Brossi et al., 2016), the nutrient availability (complex or minimal medium; Mello et al., 2016) and by the inoculum source (Cortes-Tolalpa et al., 2016).
Shotgun metagenomics is a powerful approach to investigate the lignocellulose-degrading potential of microbial communities and has substantially expanded the repertoire of genes and genomes related to plant biomass decomposition (D'Haeseleer et al., 2013; Jimenez et al., 2015; Mhuantong et al., 2015; Antunes et al., 2016; Wang et al., 2016). Furthermore, advances in multi-omics approaches are improving our capacity to explore and mine new biomass-degrading genes (Mhuantong et al., 2015). This approach has also allowed bacterial genome reconstruction of novel microbial species from natural and engineered ecosystems, favoring detailed genome-centric exploration of microbial communities (Albertsen et al., 2013; D'Haeseleer et al., 2013; Rosewarne et al., 2013; Delmont et al., 2015; Nelson et al., 2015; Antunes et al., 2016; Gupta et al., 2016; Sangwan et al., 2016; Stolze et al., 2016; Vanwonterghem et al., 2016). A genome-centric approach is one that emphasizes the analysis of complete or near-complete genomes, as opposed to catalogs of genes based on sequencing read analysis. The advantages of using a genome-centric analysis are that it provides more extensive data on which to base phylogenetic identifications (Hug et al., 2016) and the establishment and quantification of entire gene repertories in single organisms (Prosser, 2015), thereby allowing the possibility of a fuller understanding of individual microbial physiologies as well as their relationship to other organisms sharing the same ecological niche (Albertsen et al., 2016; Li et al., 2016). Based on a literature survey (Albertsen et al., 2013; D'Haeseleer et al., 2013; Delmont et al., 2015; Nelson et al., 2015; Antunes et al., 2016; Stolze et al., 2016; Vanwonterghem et al., 2016) we estimate that the number of reconstructed genomes with at least 90% completeness obtained from biomass degrading environments currently available is about 60. Given the compositional complexity of these environments, it is probably the case that thousands more additional novel species await discovery through genome reconstruction techniques applied to shotgun metagenome sequencing data.
Here, we have used shotgun metagenomics to investigate the microbial community composition of a thermophilic aerobic consortium enriched through multiple passages at 60°C in carboxymethylcellulose (CMC) as exogenous carbon source, using São Paulo Zoo Park compost as the inoculum. In previous studies (Martins et al., 2013; Antunes et al., 2016) we showed that this thermophilic composting operation harbors an impressive variety of bacterial species and metabolic functions related to biomass degradation; that result was the primary motivation for this work. The metagenomic sequences derived from the bacterial consortium under study allowed us to recover six near-complete genomes. We provide putative identifications for these six genomes, showing that at least three of them appear to be novel species. We also provide a detailed analysis of the biomass-degrading gene content in each genome.
Materials and Methods
Preparation of a Composting-Derived Microbial Consortium
The inoculum for consortium preparation derived from a composting facility of São Paulo Zoo Park, São Paulo, Brazil. The sample was collected on October/2013, after 30 days from the beginning of the composting process. Composting pile temperature was in the range of 70–72°C. The São Paulo Zoo composting facility is designed to compost 4 tons/day of all organic waste produced in the park comprising mainly shredded tree branches, leaves and grass from the maintenance of park green areas, plus manure, beddings, and food residues from about 400 species of zoo animals (mammals, avian, and reptiles). The composting process last ~100 days and is performed according a well-defined procedure (Martins et al., 2013; Antunes et al., 2016).
The compost sample was transported to the laboratory in a thermic box and immediately processed. A suspension was prepared by adding 10 g of the compost sample to 50 mL Falcon tube containing 25 mL of NaCl 0.9%. The tube was statically incubated at 60°C for 48 h. The mixture was filtered through layers of sterile cotton gauze and 150 μL of the filtrate were inoculated in 10 mL of M9 minimal medium supplemented with 1% CMC. The tube was statically incubated at 60°C for ~7–10 days, until loss of medium viscosity could be visually inspected as an indication of CMC consumption. The microbial cells were pelleted by centrifugation (8,000 × g) and ressuspended in 1 mL of fresh M9 minimal medium and 100 μL aliquot of this cell suspension were transferred to 10 mL of fresh medium supplemented with 1% CMC. The procedure was repeated for seven times, resulting in eight sequential passages to generate the final thermophilic consortium. M9 minimal medium without CMC was used as control in the sequential passages, and no detectable microbial growth was observed. Aliquots (1 mL) of the final consortium that was named ZCTH02 were stored at −80°C in 20% glycerol or further processed for total DNA extraction and shotgun metagenomic sequencing.
Shotgun Metagenomic Sequencing
The DNA extraction from the ZCTH02 culture (1 mL-aliquot) was performed with MoBio DNA Power Soil kit (MoBio Laboratories, Carlsbad, USA) according to the manufacturer's instructions. Purity and concentration of DNA were evaluated on a ND-1000 spectrophotometer (Nano Drop Technologies, Wilmington, USA). Additional quantification was performed with Quant-iT Picogreen dsDNA assay kit (Life Technologies, Grand Island, USA). DNA integrity was verified using 2100 Bioanalyzer DNA 7500 chip (Agilent Technologies, Santa Clara, USA). Shotgun metagenomic library was generated using the Illumina Nextera DNA library preparation kit (Illumina, Inc., San Diego, USA) as recommended by the manufacturer using 35 ng of total DNA. Size and quality of DNA fragment library were assessed using 2100 Bioanalyzer Agilent High Sensitivity DNA chip. Library quantification with KAPA Library Quantification Kit (Kapa Biosystems, Wilmington, USA), normalization, dilution and pooling were performed following standard protocols for sequencing in the Illumina MiSeq platform (Illumina, Inc., CA). Sequencing run was performed the MiSeq Reagent kit v2 (500-cycle format, paired-end (PE) reads). Illumina PE read1 and read2 presented, respectively, >75 and >70% of bases with quality score above phred 30 (Q30).
Processing, Quality Control, and Assembly of the Metagenomic Sequences
Raw PE sequencing reads were quality-filtered to remove reads shorter than 150 bp and reads with average quality score lower than phred 30 using SICKLE and default parameters (Joshi and Fass, 2011). Although these parameters are restrictive, they warrant better accuracy in the metagenomic assembly process (Sharon et al., 2015). Assembly of merged quality-filtered paired-end reads was performed with SPAdes using Illumina paired end reads and multi-cell data set parameters (Nurk et al., 2013). Contigs smaller than 1,600 bp were removed from downstream analyses. The coverage of each contig was calculated by mapping all high quality reads back to the final assembly by using Bowtie2 (Langmead and Salzberg, 2012).
Genome Reconstruction Based on Binning Methods
The reconstruction of individual genomes was performed with MaxBin software (Wu et al., 2014). This approach combines tetranucleotide frequencies and contig coverage to cluster metagenomic contigs into individual bins, which were validated by the identification of marker genes for each bin. The quality control of the reconstructed genomes was evaluated with CheckM (Parks et al., 2015).
Taxonomic Assignment and Phylogenetic Analysis
For taxonomic assignment and phylogenetic analysis of each reconstructed genome, a strategy based on two steps was employed. The first step was based on the comparative analysis of 16S rRNA and a single gene marker (Albertsen et al., 2013) between nr database using BLAST (Altschul et al., 1997). After the identification of similar microorganisms with each individual genome, the second step was to select phylogenetically close microbial genomes and apply the phylogenetic analysis. Phylogenetic reconstruction was conducted using PhyloPhlAn software (Segata et al., 2013) and trees were generated using iTol (Letunic and Bork, 2016). Digital DNA-DNA hybridization (dDDH) was used to compare genome-to-genome similarity and each reconstructed genome was compared to their closest genome using GGDC tool (Auch et al., 2010).
Functional Annotation
Genomes were annotated using an upgraded version of the NCBI Prokaryotic Genome Automatic Annotation Pipeline (Angiuoli et al., 2008). Protein sequences were compared against Clusters of Orthologous Groups (COGs) database (Galperin et al., 2015) using rpsblast+ (Altschul et al., 1997), with a cut-off e < 10−2. COG categories were assigned for the best hits with cdd2cog script1.
The amino acid sequences of the predicted coding sequences CDSs were screened using the dbCAN HMM-based database (Yin et al., 2012) for carbohydrate-active enzymes (CAZymes; Lombard et al., 2014) and protein domains implicated in lignocellulose degradation with HMMER package (Eddy, 2009). If alignment length >80 amino acids, the cut-off e-value used was 10−5, otherwise 10−3 was used. Data visualization was performed using Circos software (Krzywinski et al., 2009).
Mapping of Reconstructed Genomes against Metagenomic and Metatranscriptomic Datasets
Sequencing reads from previous metagenomics and metatranscriptomics studies of the thermophilic composting cell ZC4 (Antunes et al., 2016) were used for mapping. These datasets consist of nine time-series samples from composting cell ZC4 (days 1, 3, 7, 15, 30, 64, 67, 78, and 99). These sequencing reads were mapped against the reconstructed genomes to investigate their presence and activity in ZC4. For the purposes of determining the presence of reconstructed genomes in ZC4 data, we computed genome coverage under the following definition. A given base pair position in a contig that is part of a reconstructed genome was considered covered if at least one ZC4 sequencing read was mapped to that position. This computation was done using program genomecov in package BEDTools (Quinlan and Hall, 2010). For the purposes of determining relative abundance of reconstructed genomes in the metagenomic datasets we mapped the sequencing reads against the reconstructed contig sequences. For the purposes of determining relative abundance of reconstructed genomes in the metatranscriptomic datasets we mapped the sequencing reads against coding sequence regions annotated in each reconstructed genome. Abundance was measured in terms of number of reads mapped. Mapping was performed with BBMap2. Only reads with length greater or equal to 100 bp were taken into account. A mapping was considered valid if the alignment had at least 95% identity.
Accession Numbers
The recovered genomes were deposited in GenBank under the following accession numbers: LZRU00000000 (BZ1), LZRT00000000 (BZ2), LZRV00000000 (BZ3), LZRS00000000 (BZ4), SAMN05223358 (BZ5), and LZRQ00000000 (BZ6). Unassembled paired-end sequence reads of ZCTH02 metagenome are available at MG-RAST under the under the accession mgm4570098.3 (read1) and mgm4570099.3 (read2).
Results
Reconstruction of Six Genomes from a Compost-Derived Bacterial Consortium
A total of 3,046,968 sequence paired-end reads was obtained from a thermophilic and cellulolytic consortium (ZCTH02; Table 1). Taxonomic analysis of all protein-coding sequences of the ZCTH02 metagenome at the MG-RAST webserver (Keegan et al., 2016) show that it is dominated by Bacteria (99% of total sequences), most of which (86%) are classified as Firmicutes. At the order level, the consortium metagenome protein-coding sequences were identified as belonging primarily to Bacillales (72%) and Clostridiales (10%).
Table 1. Assembly metrics of shotgun metagenomics of the ZCTH02 consortium.
The metagenome shotgun paired-end reads were assembled into 1,468 contigs ≥1,600 bp, comprising 1,263,585 paired-end reads totalizing 0.206 Gbp (Table 1). Using 97% of these contigs we were able to reconstruct six nearly-complete bacterial genomes, which were named BZ1, BZ2, BZ3, BZ4, BZ5, and BZ6 (Table 2). We remark that the length-filtering parameter used during the assembly steps was quite stringent, in order to minimize the presence of chimeric contigs. Quality control (Parks et al., 2015) on the resulting genomes yielded the following values: at least 90% completeness and <5% contamination in all cases.
Table 2. Genomic features of bacterial genomes reconstructed from the ZCTH02 consortium shotgun metagenome.
The reconstructed genomes range in size from 2.7 to 4.4 Mbp, in %GC from 43.4 to 65.5%, and in coding sequences (CDSs) from 2,969 to 4,376 (Table 2). All recovered genomes belong to the Firmicutes phylum. Measuring abundance in terms of reads actually used to assemble the 1,468 binned contigs, BZ1 was the most abundant of all six reconstructed genomes (33%), followed by BZ2 (22%).
Taxonomic and Phylogenetic Assignment of Reconstructed Genomes
For each reconstructed genome, we sought to determine their taxonomic classification and phylogenetic placement using (i) BLAST similarity analysis based on 16S rRNA and DNA primase nucleotide sequences; and (ii) phylogenetic tree reconstruction based on the amino acid sequences of a small set of housekeeping genes. The BLAST results are shown in Table 2 and the phylogenetic trees are shown in Figure 1. Taxonomic and phylogenetic analyses suggest that three of the six reconstructed genomes belong to novel bacteria, with the following preliminary identifications: a new Thermobacillus species (BZ1); a member of a new genus of the Paenibacillaceae family (BZ3), similar to Paenibacillus and Cohnella; and a member of a new deep-branching family in the Clostridia class (BZ6). In addition, we provide a tentative identification of BZ2 as Bacillus thermozeamaize. If this identification is correct this reconstructed genome would be the first for this species, for which currently only 16S rRNA sequences are available. Except for BZ6, the other reconstructed genomes are members of the Bacilli class. We provide more details of these results in what follows.
Figure 1. Phylogenetic analysis of the six reconstructed bacterial genomes. The analyses were based on ~400 conserved single-copy protein sequences, selected among microbial type strain genomes phylogenetically close to BZ1 and BZ3 (A), BZ2 and BZ5 (B), BZ4 (C), and BZ6 (D). Black dots indicate bootstrap values of ≥80%.
For BZ1 the best 16S rRNA hit (99% identity) was to a sequence from a non-culturable microorganism found in an autothermal thermophilic aerobic digester in Australia (Hayes et al., 2011). Thermobacillus composti was the best hit (89% identity) for the DNA primase gene sequence. Previous studies based on 16S rRNA (Ash et al., 1993) and housekeeping genes (Zhang and Lu, 2015) have shown that the Thermobacillus genus is close to the Paenibacillus genus. Based on these results we carried out the phylogenetic reconstruction using members of the Paenabacillaceae family (Figure 1A). Comparative genome analysis based on GGDC between BZ1 and T. composti KWC4 reference genome suggests that they are distinct species (Difference in % G+C = 3.77). The inferred phylogenetic tree and these results indicate that BZ1 is likely a new species in the Thermobacillus genus.
For BZ2, the best 16S rRNA hit (99% identity) was to a sequence from Bacillus thermozeamaize strain L-10997 (accession AY288912.1). This strain was isolated from batch fermentations samples of thermophilic and hyperthermophilic food processing facilities (Mak, 2003). A closely related B. thermozeamaize isolate was also found among thermophilic cellulose-degrading cultivable bacteria from a deep subsurface gold mine (Rastogi et al., 2009). Our BLAST searches using the DNA primase gene sequence did not yield any good matches. In GenBank there are no reports of B. thermozeamaize genomes. Phylogenetic analysis (Figure 1B) shows that BZ2 is closer to Effusibacillus and Thermoactinomyces than to Bacillus. The BLAST analysis and the phylogenetic analysis point to three possible conclusions: (1) BZ2 is a B. thermozeamaize strain but the classification of B. thermozeamaize as belonging to the genus Bacillus is wrong; or (2) BZ2 is not a B. thermozeamaize strain, and the 16S rRNA similarity is misleading; or (3) the classification of the 16S rRNA sequence in record AY288912.1 is incorrect. If BZ2 is not a B. thermozeamaize strain then it would likely be the first representative of a new family in the Bacillales order.
For BZ3, the best hit for the 16S rRNA gene was with an uncultured bacterium detected in asparagus straw compost in YongJi, ShanXi, China (accession JQ775380.1); however, the query was not a complete 16S rRNA gene sequence. DNA primase analysis shows a distant relationship to genes from members of the Paenibacillus genus (79% identity/8% coverage). Phylogenetic analysis (Figure 1A) indicates that BZ3 is a member of the Paenibacillaceae family and might be a member of a new genus.
BZ4 and BZ5 appear to be new strains of known species. 16S rRNA and DNA primase analyses (Table 2) show that BZ4 and BZ5 are very similar to Geobacillus thermoglucosidasius (99% identity) and Caldibacillus debilis (100%), respectively. Phylogenetic analyses (Figures 1C,B) and GGDC analysis confirm this result (Table 2). The genomes for three strains of G. thermoglucosidasius (DSM2542, C56-YS93, and TNO-09.020) are available (Zhao et al., 2012; Brumm et al., 2015; Chen et al., 2015). C. debilis has a genome draft published (Berendsen et al., 2016).
Figure 2. Functional profile of the six reconstructed genomes based on COG categories and CAZy families. (A) Abundance of COGs of each COG functional category based on relative abundance of genes annotated per genome. (B) Abundance of CAZy families based on relative abundance of genes annotated per genome. The figures were drawn based on data shown in Table S1 (A) and Table 4 (B).
For BZ6, the best 16S rRNA hit (98% identity) was to a partial sequence (accession FN667161.1) from an uncultured bacterium in a composting study in Finland (Partanen et al., 2010). We did not get any hits for the BLAST search using the DNA primase nucleotide sequence; using the translated sequence the best hit (66% identity) was to a Clostridiales bacterium (accession KKM09071.1). Phylogenetic analysis (Figure 1D) shows that BZ6 belongs to a divergent group of Clostridia: Clostridiales families XVII (Thermaerobacter marianensis) and XVIII (Symbiobacterium thermophilum) Incertae sedis. Currently these two families are classified as a sister group to other families in the order Clostridiales (Zhang and Lu, 2015). We hypothesize that the organism corresponding to BZ6 is a member of a new deep-branching family within the Clostridia class.
Metabolic Potential Encoded by the Reconstructed Genomes
We analyzed the six reconstructed genomes using Cluster of Orthogous Groups, or COGs (Galperin et al., 2015; Figure 2A and Table S1). BZ4 presents the largest genome and therefore the highest number of CDSs (3,771) classified in COGs. CDSs classified as COG category G (metabolism and transport of carbohydrates) were more numerous in the BZ1 and BZ3 genomes, while CDSs classified as category E (amino acid metabolism and transport) were more numerous in the BZ4 and BZ6 genomes. Among the CDSs assigned to category G in BZ1 and BZ3, the most numerous had the following functional descriptions: xylanase, β-xylosidase, glycosidase, and arabinofuranosidase (Table 3). Genes encoding these predicted functions are likely playing a role in plant biomass biodegradation and some of them were also found in the BZ4 (xylanases and β-xylosidases) and BZ5 genomes (β-xylosidases). We remark that β-xylosidases are a key enzyme for the degradation of the main hemicellulose constituent, xylan, and relatively few xylosidases from thermophilic bacteria have been reported (Shao et al., 2011; Anand et al., 2013; Bhalla et al., 2014). Related to cellulose degradation, CDSs for two key enzymes were annotated: endo-1,4-beta-glucanase (cellulase M) in BZ2, BZ4, BZ5, and BZ6 and beta-glucosidase in BZ1, BZ3, BZ4, BZ5, and BZ6. CDSs encoding multicopper oxidases, which are involved in lignin breakdown, were also identified in BZ1, BZ2, BZ3, BZ4, and BZ6 (Table 3).
Table 3. Number of CDSs assigned to COGs related to lignocellulose metabolism in reconstructed genomes.
Lignocellulose Breakdown Potential Viewed through Cazy Database
We screened the six reconstructed genomes for genes encoding lignocellulose-degrading enzymes using the CAZy database (Figure 2B). We found 691 different CAZyme genes, encompassing all six CAZy families, as follows: 33% glycoside hydrolases (GHs), 24% glycosyltransferases (GTs), 20% carbohydrate esterases (CEs), 1% polysaccharide lyases (PLs), 5% auxiliary activities (AAs), and 18% carbohydrate-binding modules (CBMs; Table 4).
Table 4. Number of CDSs from reconstructed genomes mapped to enzymes in the CAZy database.
The BZ1 genome shows the highest number of CDSs classified as CAZymes (196 CDSs, 6%) followed by BZ3 (128 CDSs, 4%; Figure 2B and Table 4). These numbers are within the range of CAZymes encoding-genes estimated for the genomes of thermophilic lignocellulose degraders such as Thermobispora bispora DSM 43833 (181 CDSs, 5%), Rhodothermus marinus DSM 4252 (182 CDSs, 6.5%), and Clostridium clariflavum str. 4-2a (478, 12%; D'Haeseleer et al., 2013; Rooney et al., 2015; Hiras et al., 2016). BZ1 also stands out for being the richest in diversity and abundance of GHs compared to the other reconstructed genomes (80 CDSs, 35%). For instance, BZ1 has relatively more GHs than BZ4, which is the largest among the reconstructed genomes (BZ1: 0.02 per kbp; BZ4: 0.007 per kbp).
The other five species contribute collectively with 14 GH families not encoded in the BZ1 genome (GH5, GH26, GH38, GH39, GH52, GH57, GH65, GH74, GH76, GH78, GH88, GH108 and GH115, and GH120). Among GH families present in BZ1, GH109 is the most numerous (13 CDSs, 16%), followed by GH43 (9 CDSs, 11%). The GH109 family contains only the α-N-acetylgalactosaminidase enzyme. The GH43 family includes enzymes related to hemicellulose degradation such as β-xylosidase, α-L-arabinofuranosidase, arabinase, and xylanase. The BZ1 genome also presents three GH10 family enzymes, including a potential thermostable celloxylanase that can be involved with xylan degradation. CDSs assigned to cellulase families were found only in the BZ4 (GH5) and BZ1/BZ3 (GH9) genomes.
Overall, nine of the GH families (GH5, GH26, GH57, GH65, GH74, GH88, GH108, GH115, and GH120) are specific, occurring only in one of the five BZ genomes other than BZ1. For instance, the GH57 family (α-amylase) was identified only in the BZ6 genome (Table 4). Among these specific GH families, GH26 (in BZ3) is worth noting, because it includes exo-acting β-mannanase and β-1,3-xylanase, which are involved in hemicellulose degradation (Araki et al., 2000; Cartmell et al., 2008).
The BZ1 genome also encodes the greatest number of CBMs (25%), followed by BZ3 (22%), as computed with respect to all CDSs classified in this family, in the six BZ genomes. The most abundant family in all BZ genomes is CBM50 (80 CDSs, 65%; Table 4). CBM50 is a binding domain associated with chitin or peptidoglycan cleavage, present in several GH families such as GH18, GH23, and GH73. Accordingly, BZ4 and BZ5 encode one CDS classified as GH18 that contains a CBM50 domain. Three CBM families with members encoded in the reconstructed genomes are cellulose-binding (CBM9, CBM16, and CBM30) and one is xylan-binding (CBM22).
We found five AA families (AA1, AA2, AA4, AA6, and AA7) in the six reconstructed genomes. BZ4 possesses the highest number of CDSs assigned to this family (31%). AA4 and AA6 are the most abundant, representing 37 and 28%, respectively, of the total CDSs classified as an AA family. BZ3 and BZ6 encode a cellulose oxidase (AA7), and BZ1, BZ3, BZ4, and BZ6 have the highest number of lignin oxidases (AA1, AA2, AA4, and AA6). BZ2 and BZ5 have the smallest repertoire of genes encoding enzymes for cellulose and lignin oxidation.
Consortium in Action as Revealed by Thermophilic Composting Metagenomics and Metatranscriptomics Data
In a previous work, we studied the microbial community of the composting operation at the Sao Paulo Zoo Park (Antunes et al., 2016). Since the inoculum used to prepare the ZCTH02 consortium originated from a sample of that same composting operation, we had the possibility of verifying whether the genes we studied in this work were present during the actual composting process. First, we verified the presence of the BZ genomes in the metagenomic datasets obtained from the ZC4 composting cell. All BZ genomes except BZ2 had genome coverage values exceeding 70% in at least one sample; for BZ2 the maximum value observed was 38% (Table S2). We interpret these results as providing strong evidence for the presence of the organisms corresponding to the BZ genomes except BZ2 in composting cell ZC4. For BZ2 the evidence is weaker.
Next, we determined the relative abundance of the BZ genomes in the same metagenomic datasets as well as in the metatranscriptomic datasets (Table S3). The results are plotted in Figure 3. The variation between DNA and mRNA levels follows a similar pattern for each one of the reconstructed genomes, providing additional evidence for the presence of these organisms and evidence that they are active throughout the composting process. The number of metatranscriptomic reads mapping to BZ1 shows peaks in days 1 and 15 of composting. BZ2 has a peak in the last day (99), when the compost is considered mature. BZ3 appears more active in later stages, and BZ4 and BZ5 appear to be more active during initial stages (day 1 and until day 15, respectively). BZ6 seems to be present at moderate levels throughout the process.
Figure 3. Variation of metagenome and metatranscriptome reads mapped to the six reconstructed genomes over days of composting. Relative abundance of reads (%) was calculated using total reads of each indicated genome per total reads in the metagenomic (DNA) or metatranscriptomic (RNA) sequences in samples collected from the Sao Paulo Zoo composting process and described previously (Antunes et al., 2016). Shaded bars indicate days of composting where the respective genomes were more abundant.
We observed interesting correlations between the metabolic potential of the reconstructed genomes and the variation in their activity given by the time-series metatranscriptome mapping (Figure 4 and Table S4). On day 15, a point at which composting of easily degradable organic material is in full swing, a 1,4-beta xylanase is the fourth most abundant transcript mapped to the BZ1 genome. The mapping also shows BZ3 to be more active toward the middle and end phases of composting (days 30, 78, and 99). The fact that BZ3 encodes a multicopper oxidase, which is a lignin-degrading enzyme, and the fact that lignin becomes more accessible toward the latter stages of the composting might explain this correlation.
Figure 4. Relative abundance of time-series composting metatranscriptome sequence reads mapped to lignocellulose degrading genes in the reconstructed genomes. (A) The numbers in each day column refer to relative abundance of CDSs representing different enzymatic functions expressed in per thousand to which metatranscriptomic reads were mapped. (B) Colored pie charts show the amount of normalized reads mapped to lignocellulose-related enzyme genes from the six reconstructed genomes over days of thermophilic composting indicated in colored boxes. Red arrows indicate turning step (aeration of the composting pile).
Discussion
In this study, we used shotgun metagenomics to investigate a new thermophilic compost-derived consortium adapted for CMC digestion. This consortium was dominated by Firmicutes. The relative abundance of Firmicutes in diverse lignocellulose-degrading environments (Eichorst et al., 2014; Antunes et al., 2016; Güllert et al., 2016; Stolze et al., 2016; Zhu et al., 2016) underscores their importance and unique adaptation abilities. Previous studies on thermophilic aerobic lignocellulolytic consortia derived from green-waste composting have also shown the presence of Firmicutes, but in addition reveal the presence of other phyla, such as Bacteroidetes, Deinococcus-Thermus, and Actinobacteria, in varied proportions depending on the lignocellulosic substrate used for the enrichment (D'Haeseleer et al., 2013; Eichorst et al., 2013, 2014; Simmons et al., 2014a).
We were able to reconstruct six nearly-complete genomes from sequencing data obtained from this consortium, all of which were determined to be Firmicutes species. Four of these genomes are novel (new Thermobacillus species; the first B. thermozeamaize genome or else the first genome of a representative of a novel family in order Bacillales; a new genus of the Paenibacillaceae family and one new deep-branching family in the Clostridia class), and two are possibly new strains of G. thermoglucosidasius and C. debilis. In some cases, the reconstructed genomes belong to organism groups containing well-known biomass degraders or that play other roles in the biodegrading process. Members of the Thermobacillus genus are aerobic, and are known as plant biomass degraders, producing robust glycoside hydrolases and esterases that are of special interest for industrial applications (Wongwilaiwalin et al., 2010; Rakotoarivonina et al., 2012). G. thermoglucosidasius C56-YS93 is known as a biomass degrader and was isolated from the Obsidian hot spring in Yellowstone National Park (Brumm et al., 2015). A noncellulolytic facultative anaerobe C. debilis strain (GB1) isolated from an air-tolerant cellulolytic consortium has been studied for its ability to supply respiratory protection for cellulolytic bacteria, such as Clostridium thermocellum, when co-cultured (Wushke et al., 2015).
Regarding metabolic potential for lignocellulose degradation a diversity of cellulases (in families GH5 and GH9), endohemicellulases (GH8 and GH11), debranching enzymes (GH51 and GH67), and oligosaccharide-degrading enzymes (GH1, GH2, GH3, GH42, and GH43) were found in the reconstructed genomes. GH5 and GH9 include several enzymes related to cellulose and hemicellulose breakdown and have been described in biomass degrading metagenomes such as sugarcane bagasse (Mhuantong et al., 2015) and thermophilic composting (Antunes et al., 2016), as well as in some consortia enriched from composting (Allgaier et al., 2010; Wang et al., 2016). However, typical cellulases in the GH6 and GH48 families (endoglucanases and cellobiohydrolases) were not found to be encoded in the reconstructed genomes, a result consistent with a study of the bagasse microbiome (Mhuantong et al., 2015), which did not detect these families.
We identified BZ1 as a Thermobacillus. In the literature, the best studied Thermobacillus is T. xylanilyticus (Watanabe et al., 2007). T. xylanilyticus is a thermophilic and highly xylanolytic bacterium isolated from farm soil (Touzel et al., 2000; Paës and O'Donohue, 2006; Paës et al., 2008; Rakotoarivonina et al., 2011, 2012). In these reports, two enzymes were characterized: xylanase (GH11; Paës and O'Donohue, 2006) and α-L-arabinofuranosidase (GH51; Paës et al., 2008), both of which are present in the BZ1 genome. Another study (Rakotoarivonina et al., 2011) identified a hemicellulase enzyme (feruloyl-esterase) in T. xylanilyticus, adding evidence of the lignocellulose biomass degradation capabilities of this species. The most recent study reported that this species secretes a diverse arsenal of lignocellulolytic enzymes depending on the biomass composition used as carbon source (Rakotoarivonina et al., 2012). These studies were conducted with isolated strains. Therefore, this work is the first report of Thermobacillus in a thermophilic consortium.
We found that the BZ1 genome encodes by far the largest number of biomass degrading-related genes as compared to the other genomes. For example, BZ1 encodes more than twice as many GH-classified CDSs as any single other genome (Figure 2B). This suggests that BZ1 has a primary role in lignocellulose breakdown in this consortium. On the other hand, there are 14 GH families for which the other genomes do have representatives while BZ1 has none. Based on this we speculate that BZ1 is a driver organism while the others are auxiliary in the lignocellulose breakdown process, in this consortium. Previous studies have reported that a variety of enzymes from different organisms can work together to improve lignocellulose degradation and access to energy sources (Takasuka et al., 2013; Hiras et al., 2016). Figure 5 summarizes these observations, showing an overview of lignocellulose degradation in the ZCTH02 consortium, depicting the distribution of GH and AA families among consortium component organisms.
Figure 5. An overview of lignocellulose degradation in the ZCTH02 consortium. The diagram shows the main constituents of lignocellulose (hemicellulose, cellulose and lignin). Each colored square represents a GH family or an AA family that contains one or more CDSs that were annotated in a given BZ genome. The color of each square corresponds to the BZ genome where those CDSs were annotated, according to the key in the figure.
We presented evidence that the organisms corresponding to the reconstructed genomes are present and active in a thermophilic composting process. However, none of the genomes reconstructed here were particularly abundant in the composting cell ZC4 (Antunes et al., 2016). This most likely has to do with the particular features of the consortium construction process, with its many selection steps. Aside from high temperature, most or all of the remaining selection features of the consortium construction process are very different from those of composting. This fact, coupled with the high compositional diversity of ZC4 (Antunes et al., 2016) probably explains why our consortium did not capture the most abundant species identified in ZC4.
Overall, our results illustrate the importance and plasticity of Firmicutes members in the bioconversion of lignocellulose in compost-derived consortium selected using CMC as carbon source. The large number and diversity of CAZyme genes investigated indicate that this consortium is composed by organisms that can complement each other with different enzymes relevant for lignocellulose degradation, being therefore a promising potential source for thermostable enzymes in industrial applications.
Author Contributions
Conceived the study and designed the experiments: AMdS, JS, RQ. Performed the experiments: LL, RQ. Processed the samples and performed DNA sequencing: LFM, LA, RQ. Analyzed the data: AMdS, ARdS, LFM, LL, LMSM, RQ, RP, JS. Wrote the manuscript: AMdS, LFM, LL, RP, JS.
Funding
Funding for this research was provided by grant 2011/50870-6 from the São Paulo Research Foundation (FAPESP). LL, LA, and RP were supported by fellowships from FAPESP. LL received fellowship from FAPESP (2013/05325-5). ARdS and LMSM were supported by fellowships from the Coordination for the Improvement of Higher Education Personnel (CAPES). AMdS and JS received Research Fellowship Awards from National Council for Scientific and Technological Development (CNPq). The funders had no role in study design, data collection, analysis, decision to publish or preparation of the manuscript.
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.
Acknowledgments
The authors are grateful to Dr. João Batista da Cruz and Dr. Paulo Bressan for support to this project and for allowing sampling at the São Paulo Zoo Compost Facility.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00644/full#supplementary-material
Footnotes
1. ^https://github.com/aleimba/bac-genomics-scripts/tree/master/cdd2cog.
2. ^https://www.biostarhandbook.com/tools/bbmap/bbmap-help.html.
References
Albertsen, M., Hugenholtz, P., Skarshewski, A., Nielsen, K. L., Tyson, G. W., and Nielsen, P. H. (2013). Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat. Biotechnol. 31, 533–538. doi: 10.1038/nbt.2579
Albertsen, M., McIlroy, S. J., Stokholm-Bjerregaard, M., Karst, S. M., and Nielsen, P. H. (2016). Candidatus Propionivibrio aalborgensis: a novel glycogen accumulating organism abundant in full-scale enhanced biological phosphorus removal plants. Front. Microbiol. 7:1033. doi: 10.3389/fmicb.2016.01033
Allgaier, M., Reddy, A., Park, J. I., Ivanova, N., D'haeseleer, P., Lowry, S., et al. (2010). Targeted discovery of glycoside hydrolases from a switchgrass-adapted compost community. PLoS ONE 5:e8812. doi: 10.1371/journal.pone.0008812
Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Anand, A., Kumar, V., and Satyanarayana, T. (2013). Characteristics of thermostable endoxylanase and beta-xylosidase of the extremely thermophilic bacterium Geobacillus thermodenitrificans TSAA1 and its applicability in generating xylooligosaccharides and xylose from agro-residues. Extremophiles 17, 357–366. doi: 10.1007/s00792-013-0524-x
Angiuoli, S. V., Gussman, A., Klimke, W., Cochrane, G., Field, D., Garrity, G., et al. (2008). Toward an online repository of standard operating procedures (SOPs) for (Meta) genomic annotation. OMICS 12, 137–141. doi: 10.1089/omi.2008.0017
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
Araki, T., Hashikawa, S., and Morishita, T. (2000). Cloning, sequencing, and expression in Escherichia coli of the new gene encoding beta-1,3-xylanase from a marine bacterium, Vibrio sp strain XY-214. Appl. Environ. Microbiol. 66, 1741–1743. doi: 10.1128/AEM.66.4.1741-1743.2000
Ash, C., Priest, F. G., and Collins, M. D. (1993). Molecular identification of ribosomal RNA Group 3 Bacilli (Ash, Farrow, Wallbanks and Collins) using a PCR probe test - proposal for the creation of a new genus Paenibacillus. Antonie Van Leeuwenhoek 64, 253–260. doi: 10.1007/BF00873085
Auch, A. F., Klenk, H. P., and Goker, M. (2010). Standard operating procedure for calculating genome-to-genome distances based on high-scoring segment pairs. Stand. Genomic Sci. 2, 142–148. doi: 10.4056/sigs.541628
Berendsen, E. M., Wells-Bennik, M. H., Krawczyk, A. O., de Jong, A., van Heel, A., Holsappel, S., et al. (2016). Draft genome sequences of seven thermophilic spore-forming bacteria isolated from foods that produce highly heat-resistant spores, comprising Geobacillus spp., Caldibacillus debilis, and Anoxybacillus flavithermus. Genome Announc. 4:e00105–16. doi: 10.1128/genomeA.00105-16
Bhalla, A., Bischoff, K. M., and Sani, R. K. (2014). Highly thermostable GH39 beta-xylosidase from a Geobacillus sp strain WSUCF1. BMC Biotechnol. 14, 963. doi: 10.1186/s12896-014-0106-8
Brumm, P. J., Land, M. L., and Mead, D. A. (2015). Complete genome sequence of Geobacillus thermoglucosidasius C56-YS93, a novel biomass degrader isolated from obsidian hot spring in Yellowstone National Park. Stand. Genomic Sci. 10:73. doi: 10.1186/s40793-015-0031-z
Cartmell, A., Topakas, E., Ducros, V. M., Suits, M. D., Davies, G. J., and Gilbert, H. J. (2008). The Cellvibrio japonicus mannanase CjMan26C displays a unique exo-mode of action that is conferred by subtle changes to the distal region of the active site. J. Biol. Chem. 283, 34403–34413. doi: 10.1074/jbc.M804053200
Chen, J. Y., Zhang, Z. Z., Zhang, C. L., and Yu, B. (2015). Genome sequence of Geobacillus thermoglucosidasius DSM2542, a platform hosts for biotechnological applications with industrial potential. J. Biotechnol. 216, 98–99. doi: 10.1016/j.jbiotec.2015.10.002
Cortes-Tolalpa, L., Jiménez, D. J., Brossi, M. J. D., Salles, J. F., and van Elsas, J. D. (2016). Different inocula produce distinctive microbial consortia with similar lignocellulose degradation capacity. Appl. Microbiol. Biotechnol. 100, 7713–7725. doi: 10.1007/s00253-016-7516-6
de Lima Brossi, M. J., Jiménez, D. J., Cortes-Tolalpa, L., and van Elsas, J. D. (2016). 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
Delmont, T. O., Eren, A. M., Maccario, L., Prestat, E., Esen, O. C., Pelletier, E., et al. (2015). Reconstructing rare soil microbial genomes using in situ enrichments and metagenomics. Front. Microbiol. 6:358. doi: 10.3389/fmicb.2015.00358
D'Haeseleer, P., Gladden, J. M., Allgaier, M., Chain, P. S. G., 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
Eddy, S. R. (2009). A new generation of homology search tools based on probabilistic inference. Genome Inform. 23, 205–211. doi: 10.1142/9781848165632_0019
Eichorst, S. A., Joshua, C., Sathitsuksanoh, N., Singh, S., Simmons, B. A., and Singer, S. W. (2014). Substrate-specific development of thermophilic bacterial consortia by using chemically pretreated switchgrass. Appl. Environ. Microbiol. 80, 7423–7432. doi: 10.1128/AEM.02795-14
Eichorst, S. A., and Kuske, C. R. (2012). Identification of cellulose-responsive bacterial and fungal communities in geographically and edaphically different soils by using stable isotope probing. Appl. Environ. Microbiol. 78, 2316–2327. doi: 10.1128/AEM.07313-11
Eichorst, S. A., Varanasi, P., Stavila, V., Zemla, M., Auer, M., Singh, S., et al. (2013). Community dynamics of cellulose-adapted thermophilic bacterial consortia. Environ. Microbiol. 15, 2573–2587. doi: 10.1111/1462-2920.12159
Feng, Y. J., Yu, Y. L., Wang, X., Qu, Y. P., Li, D. M., He, W. H., et al. (2011). Degradation of raw corn stover powder (RCSP) by an enriched microbial consortium and its community structure. Bioresour. Technol. 102, 742–747. doi: 10.1016/j.biortech.2010.08.074
Galperin, M. Y., Makarova, K. S., Wolf, Y. I., and Koonin, E. V. (2015). Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 43, D261–D269. doi: 10.1093/nar/gku1223
Gao, Z. M., Xu, X., and Ruan, L. W. (2014). Enrichment and characterization of an anaerobic cellulolytic microbial consortium SQD-1.1 from mangrove soil. Appl. Microbiol. Biotechnol. 98, 465–474. doi: 10.1007/s00253-013-4857-2
Gladden, J. M., Allgaier, M., Miller, C. S., Hazen, T. C., VanderGheynst, J. S., Hugenholtz, P., et al. (2011). Glycoside hydrolase activities of thermophilic bacterial consortia adapted to switchgrass. Appl. Environ. Microbiol. 77, 5804–5812. doi: 10.1128/AEM.00032-11
Güllert, S., Fischer, M. A., Turaev, D., Noebauer, B., Ilmberger, N., Wemheuer, B., et al. (2016). Deep metagenome and metatranscriptome analyses of microbial communities affiliated with an industrial biogas fermenter, a cow rumen, and elephant feces reveal major differences in carbohydrate hydrolysis strategies. Biotechnol. Biofuels 9:121. doi: 10.1186/s13068-016-0534-x
Gupta, A., Kumar, S., Prasoodanan, V. P. K., Harish, K., Sharma, A. K., and Sharma, V. K. (2016). Reconstruction of bacterial and viral genomes from multiple metagenomes. Front. Microbiol. 7:469. doi: 10.3389/fmicb.2016.00469
Haruta, S., Cui, Z., Huang, Z., Li, M., Ishii, M., and Igarashi, Y. (2002). Construction of a stable microbial community with high cellulose-degradation ability. Appl. Microbiol. Biotechnol. 59, 529–534. doi: 10.1007/s00253-002-1026-4
Hayes, D., Izzard, L., and Seviour, R. (2011). Microbial ecology of autothermal thermophilic aerobic digester (ATAD) systems for treating waste activated sludge. Syst. Appl. Microbiol. 34, 127–138. doi: 10.1016/j.syapm.2010.11.017
Hemsworth, G. R., Johnston, E. M., Davies, G. J., and Walton, P. H. (2015). Lytic polysaccharide monooxygenases in biomass conversion. Trends Biotechnol. 33, 747–761. doi: 10.1016/j.tibtech.2015.09.006
Hess, M., Sczyrba, A., Egan, R., Kim, T. W., Chokhawala, H., Schroth, G., et al. (2011). Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science 331, 463–467. doi: 10.1126/science.1200387
Hiras, J., Wu, Y. W., Deng, K., Nicora, C. D., Aldrich, J. T., Frey, D., et al. (2016). Comparative community proteomics demonstrates the unexpected importance of actinobacterial glycoside hydrolase family 12 protein for crystalline cellulose hydrolysis. Mbio 7:e01106–16. doi: 10.1128/mBio.01106-16
Hug, L. A., Baker, B. J., Anantharaman, K., Brown, C. T., Probst, A. J., Castelle, C. J., et al. (2016). A new view of the tree of life. Nat. Microbiol. 1:16048. doi: 10.1038/nmicrobiol.2016.48
Jimenez, D. J., Chaves-Moreno, D., and van Elsas, J. D. (2015). 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., Korenblum, E., and van Elsas, J. D. (2014). Novel multispecies microbial consortia involved in lignocellulose and 5-hydroxymethylfurfural bioconversion. Appl. Microbiol. Biotechnol. 98, 2789–2803. doi: 10.1007/s00253-013-5253-7
Joshi, N. A., and Fass, J. N. (2011). Sickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ Files. GITHUB. Available online at https://github.com/najoshi/sickle
Keegan, K. P., Glass, E. M., and Meyer, F. (2016). MG-RAST, a Metagenomics service for analysis of microbial community structure and function. Methods Mol. Biol. 1399, 207–233. doi: 10.1007/978-1-4939-3369-3_13
Kinet, R., Destain, J., Hiligsmann, S., Thonart, P., Delhalle, L., Taminiau, B., et al. (2015). Thermophilic and cellulolytic consortium isolated from composting plants improves anaerobic digestion of cellulosic biomass: toward a microbial resource management approach. Bioresour. Technol. 189, 138–144. doi: 10.1016/j.biortech.2015.04.010
Koeck, D. E., Pechtl, A., Zverlov, V. V., and Schwarz, W. H. (2014). Genomics of cellulolytic bacteria. Curr. Opin. Biotechnol. 29, 171–183. doi: 10.1016/j.copbio.2014.07.002
Krzywinski, M., Schein, J., Birol, I., Connors, J., Gascoyne, R., Horsman, D., et al. (2009). Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639–1645. doi: 10.1101/gr.092759.109
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Letunic, I., and Bork, P. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245. doi: 10.1093/nar/gkw290
Li, M., Jain, S., and Dick, G. J. (2016). Genomic and transcriptomic resolution of organic matter utilization among deep-sea bacteria in guaymas basin hydrothermal plumes. Front. Microbiol. 7:1125. doi: 10.3389/fmicb.2016.01125
Lombard, V., Ramulu, H. G., 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
López-Mondéjar, R., Zuhlke, D., Becher, D., Riedel, K., and Baldrian, P. (2016). Cellulose and hemicellulose decomposition by forest soil bacteria proceeds by the action of structurally variable enzymatic systems. Sci. Rep. 6:25279. doi: 10.1038/srep25279
Mak, P. S. (2003). Isolation and Characterization of Thermophilic and Hyperthermophilic Microorganisms from Food Processing Facilities. Master's Thesis, Iowa State University, Ames.
Martins, L. F., Antunes, L. P., Pascon, R. C., de Oliveira, J. C., Digiampietri, L. A., Barbosa, D., et al. (2013). Metagenomic analysis of a tropical composting operation at the Sao Paulo Zoo park reveals diversity of biomass degradation functions and organisms. PLoS ONE 8:e61928. doi: 10.1371/journal.pone.0061928
Mello, B. L., Alessi, A. M., McQueen-Mason, S., Bruce, N. C., and Polikarpov, I. (2016). Nutrient availability shapes the microbial community structure in sugarcane bagasse compost-derived consortia. Sci. Rep. 6:38781. doi: 10.1038/srep38781
Mhuantong, W., Charoensawan, V., Kanokratana, P., Tangphatsornruang, S., and Champreda, V. (2015). Comparative analysis of sugarcane bagasse metagenome reveals unique and conserved biomass-degrading enzymes among lignocellulolytic microbial communities. Biotechnol. Biofuels 8:16. doi: 10.1186/s13068-015-0200-8
Neher, D. A., Weicht, T. R., Bates, S. T., Leff, J. W., and Fierer, N. (2013). Changes in bacterial and fungal communities across compost recipes, preparation methods, and composting times. PLoS ONE 8:e79512. doi: 10.1371/journal.pone.0079512
Nelson, W. C., Maezato, Y., Wu, Y. W., Romine, M. F., and Lindemann, S. R. (2015). Identification and resolution of microdiversity through metagenomic sequencing of parallel consortia. Appl. Environ. Microbiol. 82, 255–267. doi: 10.1128/AEM.02274-15
Nurk, S., Bankevich, A., Antipov, D., Gurevich, A. A., Korobeynikov, A., Lapidus, A., et al. (2013). Assembling single-cell genomes and mini-metagenomes from chimeric MDA products. J. Comput. Biol. 20, 714–737. doi: 10.1089/cmb.2013.0084
Paës, G., and O'Donohue, M. J. (2006). Engineering increased thermostability in the thermostable GH-11 xylanase from Thermobacillus xylanilyticus. J. Biotechnol. 125, 338–350. doi: 10.1016/j.jbiotec.2006.03.025
Paës, G., Skov, L. K., O'Donohue, M. J., Remond, C., Kastrup, J. S., Gajhede, M., et al. (2008). The structure of the complex between a branched pentasaccharide and Thermobacillus xylanilyticus GH-51 arabinofuranosidase reveals xylan-binding determinants and induced fit. Biochemistry 47, 7441–7451. doi: 10.1021/bi800424e
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114
Partanen, P., Hultman, J., Paulin, L., Auvinen, P., and Romantschuk, M. (2010). Bacterial diversity at different stages of the composting process. BMC Microbiol. 10:94. doi: 10.1186/1471-2180-10-94
Peng, X. F., Gilmore, S. P., and O'Malley, M. A. (2016). Microbial communities for bioprocessing: lessons learned from nature. Curr. Opin. Chem. Eng. 14, 103–109. doi: 10.1016/j.coche.2016.09.003
Prosser, J. I. (2015). Dispersing misconceptions and identifying opportunities for the use of 'omics' in soil microbial ecology. Nat. Rev. Microbiol. 13, 439–446. doi: 10.1038/nrmicro3468
Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033
Rakotoarivonina, H., Hermant, B., Chabbert, B., Touzel, J. P., and Remond, C. (2011). A thermostable feruloyl-esterase from the hemicellulolytic bacterium Thermobacillus xylanilyticus releases phenolic acids from non-pretreated plant cell walls. Appl. Microbiol. Biotechnol. 90, 541–552. doi: 10.1007/s00253-011-3103-z
Rakotoarivonina, H., Hermant, B., Monthe, N., and Rémond, C. (2012). The hemicellulolytic enzyme arsenal of Thermobacillus xylanilyticus depends on the composition of biomass used for growth. Microb. Cell Fact. 11:159 doi: 10.1186/1475-2859-11-159
Rastogi, G., Muppidi, G., Gurram, R., Adhikari, A., Bischoff, K., Hughes, S., et al. (2009). Isolation and characterization of cellulose-degrading bacteria from the deep subsurface of the Homestake gold mine, Lead, South Dakota, USA. J. Indus. Microbiol. Biotechnol. 36, 585–598. doi: 10.1007/s10295-009-0528-9
Rooney, E. A., Rowe, K. T., Guseva, A., Huntemann, M., Han, J. K., Chen, A., et al. (2015). Draft genome sequence of the cellulolytic and xylanolytic thermophile Clostridium clariflavum strain 4-2a. Genome Announc. 3:e00797–15. doi: 10.1128/genomeA.00797-15
Rosewarne, C. P., Greenfield, P., Li, D., Tran-Dinh, N., Midgley, D. J., and Hendry, P. (2013). Draft genome sequence of Methanobacterium sp. maddingley, reconstructed from metagenomic sequencing of a methanogenic microbial consortium enriched from coal-seam gas formation water. Genome Announc. 1:e00082–12. doi: 10.1128/genomea.00082-12
Sangwan, N., Xia, F. F., and Gilbert, J. A. (2016). Recovering complete and draft population genomes from metagenome datasets. Microbiome 4:8. doi: 10.1186/s40168-016-0154-5
Segata, N., Börnigen, D., Morgan, X. C., and Huttenhower, C. (2013). PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat. Commun. 4:2304. doi: 10.1038/ncomms3304
Shao, W. L., Xue, Y. M., Wu, A. L., Kataeva, I., Pei, J. J., Wu, H. W., et al. (2011). Characterization of a Novel beta-Xylosidase, XylC, from Thermoanaerobacterium saccharolyticum JW/SL-YS485. Appl. Environ. Microbiol. 77, 719–726. doi: 10.1128/AEM.01511-10
Sharon, I., Kertesz, M., Hug, L. A., Pushkarev, D., Blauwkamp, T. A., Castelle, C. J., et al. (2015). Accurate, multi-kb reads resolve complex populations and detect rare microorganisms. Genome Res. 25, 534–543. doi: 10.1101/gr.183012.114
Simmons, C. W., Reddy, A. P., D'Haeseleer, P., Khudyakov, J., Billis, K., Pati, A., et al. (2014a). 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
Simmons, C. W., Reddy, A. P., Simmons, B. A., Singer, S. W., and VanderGheynst, J. S. (2014b). Effect of inoculum source on the enrichment of microbial communities on two lignocellulosic bioenergy crops under thermophilic and high-solids conditions. J. Appl. Microbiol. 117, 1025–1034. doi: 10.1111/jam.12609
Stolze, Y., Bremges, A., Rumming, M., Henke, C., Maus, I., Pühler, A., et al. (2016). Identification and genome reconstruction of abundant distinct taxa in microbiomes from one thermophilic and three mesophilic production-scale biogas plants. Biotechnol. Biofuels 9:156. doi: 10.1186/s13068-016-0565-3
Takasuka, T. E., Book, A. J., Lewin, G. R., Currie, C. R., and Fox, B. G. (2013). Aerobic deconstruction of cellulosic biomass by an insect-associated Streptomyces. Sci. Rep. 3:1030. doi: 10.1038/srep01030
Touzel, J. P., O'Donohue, M., Debeire, P., Samain, E., and Breton, C. (2000). Thermobacillus xylanilyticus gen. nov., sp nov., a new aerobic thermophilic xylan-degrading bacterium isolated from farm soil. Int. J. Syst. Evol. Microbiol. 50, 315–320. doi: 10.1099/00207713-50-1-315
Vanwonterghem, I., Jensen, P. D., Rabaey, K., and Tyson, G. W. (2016). Genome-centric resolution of microbial diversity, metabolism and interactions in anaerobic digestion. Environ. Microbiol. 18, 3144–3158. doi: 10.1111/1462-2920.13382
Vartoukian, S. R., Palmer, R. M., and Wade, W. G. (2010). Strategies for culture of ‘unculturable’ bacteria. FEMS Microbiol. Lett. 309, 1–7. doi: 10.1111/j.1574-6968.2010.02000.x
Wang, C., Dong, D., Wang, H. S., Müller, K., Qin, Y., Wang, H. L., et al. (2016). Metagenomic analysis of microbial consortia enriched from compost: new insights into the role of Actinobacteria in lignocellulose decomposition. Biotechnol. Biofuels 9:22. doi: 10.1186/s13068-016-0440-2
Watanabe, Y., Pinsirodom, P., Nagao, T., Yamauchi, A., Kobayashi, T., Nishida, Y., et al. (2007). Conversion of acid oil by-produced in vegetable oil refining to biodiesel fuel by immobilized Candida antarctica lipase. J. Mol. Catal. B 44, 99–105. doi: 10.1016/j.molcatb.2006.09.007
Wongwilaiwalin, S., Rattanachomsri, U., Laothanachareon, T., Eurwilaichitr, L., Igarashi, Y., and Champreda, V. (2010). Analysis of a thermophilic lignocellulose degrading microbial consortium and multi-species lignocellulolytic enzyme system. Enzyme Microb. Technol. 47, 283–290. doi: 10.1016/j.enzmictec.2010.07.013
Wu, Y. W., Tang, Y. H., Tringe, S. G., Simmons, B. A., and Singer, S. W. (2014). MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome 2:26. doi: 10.1186/2049-2618-2-26
Wushke, S., Levin, D. B., Cicek, N., and Sparling, R. (2015). Facultative Anaerobe Caldibacillus debilis GB1: characterization and use in a designed aerotolerant, cellulose-degrading coculture with Clostridium thermocellum. Appl. Environ. Microbiol. 81, 5567–5573. doi: 10.1128/AEM.00735-15
Yan, L., Gao, Y. M., Wang, Y. J., Liu, Q., Sun, Z. Y., Fu, B. R., et al. (2012). Diversity of a mesophilic lignocellulolytic microbial consortium which is useful for enhancement of biogas production. Bioresour. Technol. 111, 49–54. doi: 10.1016/j.biortech.2012.01.173
Yin, Y. B., Mao, X. Z., Yang, J. C., Chen, X., Mao, F. L., and Xu, Y. (2012). dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 40, W445–W451. doi: 10.1093/nar/gks479
Zhang, W., and Lu, Z. (2015). Phylogenomic evaluation of members above the species level within the phylum Firmicutes based on conserved proteins. Environ. Microbiol. Rep. 7, 273–281. doi: 10.1111/1758-2229.12241
Zhao, Y., Caspers, M. P., Abee, T., Siezen, R. J., and Kort, R. (2012). Complete genome sequence of Geobacillus thermoglucosidans TNO-09.020, a thermophilic sporeformer associated with a dairy-processing environment. J. Bacteriol. 194, 4118–4118. doi: 10.1128/JB.00318-12
Zhu, N., Yang, J. S., Ji, L., Liu, J. W., Yang, Y., and Yuan, H. L. (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 genome reconstruction, consortium, thermophilic, cellulolytic, glycoside hydrolases, composting, metagenome
Citation: Lemos LN, Pereira RV, Quaggio RB, Martins LF, Moura LMS, da Silva AR, Antunes LP, da Silva AM and Setubal JC (2017) Genome-Centric Analysis of a Thermophilic and Cellulolytic Bacterial Consortium Derived from Composting. Front. Microbiol. 8:644. doi: 10.3389/fmicb.2017.00644
Received: 20 January 2017; Accepted: 29 March 2017;
Published: 19 April 2017.
Edited by:
Eric Altermann, AgResearch, New ZealandReviewed by:
Dimitris Tsaltas, Cyprus University of Technology, CyprusAlinne Castro, Universidade Católica Dom Bosco, Brazil
Copyright © 2017 Lemos, Pereira, Quaggio, Martins, Moura, da Silva, Antunes, da Silva and Setubal. 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) or licensor 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: Aline M. da Silva, almsilva@iq.usp.br
João C. Setubal, setubal@iq.usp.br
†These authors shared senior authorship.