Metagenomic mining pectinolytic microbes and enzymes from an apple pomace-adapted compost microbial community
Biotechnology for Biofuels volume 10, Article number: 198 (2017)
Degradation of pectin in lignocellulosic materials is one of the key steps for biofuel production. Biological hydrolysis of pectin, i.e., degradation by pectinolytic microbes and enzymes, is an attractive paradigm because of its obvious advantages, such as environmentally friendly procedures, low in energy demand for lignin removal, and the possibility to be integrated in consolidated process. In this study, a metagenomics sequence-guided strategy coupled with enrichment culture technique was used to facilitate targeted discovery of pectinolytic microbes and enzymes. An apple pomace-adapted compost (APAC) habitat was constructed to boost the enrichment of pectinolytic microorganisms.
Analyses of 16S rDNA high-throughput sequencing revealed that microbial communities changed dramatically during composting with some bacterial populations being greatly enriched. Metagenomics data showed that apple pomace-adapted compost microbial community (APACMC) was dominated by Proteobacteria and Bacteroidetes. Functional analysis and carbohydrate-active enzyme profiles confirmed that APACMC had been successfully enriched for the targeted functions. Among the 1756 putative genes encoding pectinolytic enzymes, 129 were predicted as novel (with an identity <30% to any CAZy database entry) and only 1.92% were more than 75% identical with proteins in NCBI environmental database, demonstrating that they have not been observed in previous metagenome projects. Phylogenetic analysis showed that APACMC harbored a broad range of pectinolytic bacteria and many of them were previously unrecognized.
The immensely diverse pectinolytic microbes and enzymes found in our study will expand the arsenal of proficient degraders and enzymes for lignocellulosic biofuel production. Our study provides a powerful approach for targeted mining microbes and enzymes in numerous industries.
High worldwide demand for energy and increasing concerns over global climate change have prompted the development of sustainable and environmentally friendly energy [16, 45]. Lignocellulosic biofuel, which derived from the most abundant renewable organic material on our planet, represents a promising alternative to fossil fuels [12, 15]. However, the major obstacles to industrial-scale production of biofuel from lignocellulosic feedstocks lie in the recalcitrant nature of biomass toward enzymatic breakdown and the relatively low activity of currently available hydrolytic enzymes [15, 44].
Pectin is one of the plant cell wall components. It is abundant in the middle lamella and primary cell walls, though presents at low levels in secondary walls . For the cell walls of pectin-rich biomass, for example apple pomace, it contains 12–35% pectin on a dry weight basis . In plant biomass, pectin embeds in the cellulose–hemicellulose network of the cell wall and regulates intercellular adhesion like glues . It is the complex matrix of pectin that masks cellulose and/or hemicellulose through hydrogen bonding interactions , and blocks their accessibility to degradative enzymes , thus resulting in plant biomass that is less susceptible to degradation and more recalcitrant to deconstruction . As a result, degradation of pectin in lignocellulosic materials has been established as essential for efficient bioconversion of lignocellulose . Recently, the reduction of bulk percentage of pectin through genetic manipulation or enzymatic means has been proved to reduce the recalcitrance and accelerate the lignocellulose saccharification of herbaceous plants , Arabidopsis , switchgrass , and woody biomass [3, 4].
Removal of pectin can be achieved in physical, chemical, or biological manner. Biological hydrolysis of pectin by pectinolytic microbes and enzymes is favored as it is environmentally benign and energy efficient . Pectinolytic enzymes have multiple benefits in the efficient hydrolysis of lignocellulosic materials: first, yield of fermentable sugars by hydrolysis of pectin itself ; second, facilitation of sugar release by disrupting the pectin network around cellulose and lignin , and exposure of other polymers to degradation by hemicellulases and cellulases ; third, improvement of cell wall porosity  and reduction of mechanical strength because of its crosslinking and water complexation features . Especially for pectin-rich lignocellulosic biomass which also could serve as the feedstock for lignocellulosic biofuel , for instance apple pomace, pectinolytic enzymes will play a more prominent role.
Despite pectinolytic enzymes playing a crucial part in the lignocellulosic biofuel production, most of the currently available pectinolytic enzymes are costly, inefficient, and susceptive to fluctuations in feedstock . In consequence, search for microbes and enzymes from naturally evolved pectinolytic microbial communities offers a promising strategy for the discovery of new pectinolytic enzymes. Given the unique features of compost habitat, there is tremendous potential to discover robust organisms and novel enzymes which tolerate harsh pretreatment scenarios under industrial conditions . Thus, compost is considered as one of the most attractive DNA pools for target gene discovery [2, 24].
Metagenomics, which directly analyzes the total DNA from environmental samples, provides a powerful strategy in unveiling the novel microbes and enzymes in microbial communities without the technical challenges of cultivation [15, 37]. However, as environmental samples generally hold a huge reservoir of extensive microbes and enzymes, it is unfeasible to characterize them accurately. Hence, to reduce the complexity of metagenomic datasets, render further assembly more amenable, and more importantly, improve the specificity of the sample’s DNA, the oriented enrichment culture technique is essential to be employed for the establishment of microbial consortia with desired functionality . In this manner, the enzyme repertoire of enriched consortia can be tailored to degrade specific feedstock . Since apple pomace is an pectin-rich lignocellulosic biomass , pectinolytic enzymes could be exploited from the established pectinolytic microbes which are selectively enriched in abundance from compost communities by cultivation with apple pomace as the sole carbon source.
In this study, a metagenomic sequence-guided strategy combined with enrichment culture technique was used to targetedly discover the pectinolytic microbes and enzymes from an apple pomace-adapted compost microbial community (APACMC). The pipeline of this strategy is shown in Fig. 1. Firstly, the unique APACMC was constructed from the cow manure compost habitat to boost the enrichment of pectinolytic microorganisms. The dynamic microbial changes of APACMC were characterized by 16S rDNA high-throughput sequencing. Secondly, a targeted metagenomic approach was applied to facilitate the identification of pectinolytic microbes and enzymes. A more accurate microbial taxonomic analysis and function characterization were conducted. Thirdly, the metagenome sequences were annotated and phylogenetically affiliated against carbohydrate-active enzymes (CAZymes) database. Finally, after the specific investigation of genes related to pectinolytic CAZymes and their taxonomic affiliations, the robust microorganisms and novel enzymes processing the degradation of pectin were identified.
Results and discussion
Changes in physicochemical properties during composting
The variations of physicochemical properties during composting are strongly associated with the biological reactions involving organic matter, and thus, these changes reflect the microbial activity and progress of the composting process . The dynamic changes in physicochemical properties (i.e., temperature, water content, and pH) during the apple pomace-adapted compost (APAC) process are illustrated in Fig. 2. During the 30-day enrichment period (Fig. 2a), the temperature of the APAC pile maintained at 25–35 °C for 24 h to allow the compost microbes to establish, then it rapidly reached 60–70 °C to trigger the thermophilic phase. After the temperature reached 68 °C at day 15, the temperature declined gradually back to the ambient temperature over the rest 15 days to trigger the cooling and maturation phase. As shown in Fig. 2b, the initial pH of APAC was in the range of 3.8–4.0 as the acid–base nature of apple pomace. Eventually, the pH value of APAC gradually rose to approximately 8.5. The escalating pH during composting may be attributed to the release of ammonia, methanol, and the decomposition of organic acids of apple pomace . The water content of APAC dropped fast at the early stages and then declined slowly. The small variations of pH and water content at the end of enrichment indicated that the microorganisms were still active and the degradation of apple pomace continued.
Changes in bacterial community structure during composting
To characterize the changes of microbiota structure during composting, 16S rDNA sequencing on the representative samples of different phases, i.e., CM0 (day 0), Mes5 (day 5 of Mesophilic), The15 (day 15 of Thermophilic), and Mat30 (day 30 of Maturation), was performed. As expected, the sample CM0 had the highest α diversity (OUT numbers and Chao1 estimator) while The15 had the lowest (Additional file 1: Table S1). Although the main phyla throughout the entire composting process did not change greatly (Fig. 2c), namely Actinobacteria, Bacteroidetes, Proteobacteria, Firmicutes, and Cyanobacteria, the main genera varied dramatically with some bacterial populations being greatly enriched (Fig. 2d). At the genus level, the bacterial community profiles of the main genera were clustered into two groups, which the bacterial community structures in The15 and Mat30 differed remarkably from CM0. The genera Acinetobacter (21.9%), Planococcaceae unclassified (8.0%), and Ruminococcaceae UCG-005 (6.9%) were the dominant in the CM0, whereas they declined to a very low level or completely disappeared after day 15 (Fig. 2d). By contrast, the genera Pseudoxanthomonas (36.7%), Parapedobacter (16.8%), Chelatococcus (7.2%), Olivibacter (4.7%), and Sphingobacterium (2.4%) were enriched in the composting process. The evolution of specific populations reflected that APAC had adapted to apple pomace degradation. Furthermore, most of the species abundant in Mat30 have been detected showing highly positive correlations on lignocellulose/pectin-degrading activities [9, 36, 44, 52]. The predominance of lignocellulolytic or pectinolytic species suggests that APACMC has potential to degrade lignocellulose and pectin effectively.
Moreover, the pectinolytic activities of APAC increased dramatically after the composting process, which preliminarily proved the effectiveness of the microbial enrichment (Additional file 2: Fig. S1). Further, it was showed that 83.25% pectin was degraded through 30-day composting. Consequently, APACMC was supposed to be successfully established with pectinolytic capability by means of the enrichment culture technique we adopted.
Microbial diversity in APACMC metagenome
To obtain more detailed information on the diversities of pectinolytic microbes and genes encoding pectinolytic enzymes in APACMC, shotgun sequence of Mat30 was performed by using Illumina HiSeq4000 platform. Metagenomic sequencing of APACMC yielded 89,623,103 reads after quality filtering. After assembly, 272,516 predicted ORFs with the average length of 668 bp were obtained (Additional file 3: Table S2). The analysis of metagenomic datasets showed that APACMC was predominately composed of bacterial members (~99.7%), along with very few archaea, eukarya, and uncharacterized organisms. This could be explained by the fact that fungal activities were precluded as APACMC sustained high temperatures between 55 and 68 °C (Fig. 2a).
To estimate the microbial diversity more accurately, various taxonomic protocols such as MEGAN, MetaPhlAn, and MG-RAST were used, while minor differences in the rank abundance order were observed. Taxonomic analysis revealed that APACMC was primarily consisted of members from phyla Proteobacteria, Bacteroidetes, Actinobacteria, and Firmicutes (Fig. 3a), which agreed with the result of 16S rDNA sequencing (Mat30 in Fig. 2c). Several previous studies have reported that thermophilic compost communities contain high-abundance genera within these phyla [18, 36]. Other phyla, such as Verrucomicrobia, Cyanobacteria, and Planctomycetes were presented at very low abundances, which together accounted for only 1.04% of the total sequences. Meanwhile, around 623 predicted genes could not be assigned to a definite bacteria phylum, which may belong to yet uncharacterized bacteria.
The prevalence of genus Sphingobacterium (17.11%), Pseudoxanthomonas (13.46%), Bordetella (4.83%), Microbacterium (2.23%), Pedobacter (2.13%), Niabella (2.10%), and Thermobacillus (1.33%) was in accordance with previous studies, which have found these species to be major components of bacterial consortia with high cellulolytic or pectinolytic activity in compost habitats [9, 18, 36, 40]. The taxonomic compositions at genus level were also consistent with the result of 16S rDNA sequencing (Mat30 in Fig. 2d). In addition, it was found some Enterobacteriaceae species also had high capacity for metabolizing pectin, such as genera Yersinia, Klebsiella, Dickeya, and Pectobacterium . Moreover, the thermophilic pectinolytic bacteria, such as Actinomadura , Bacillus , Geobacillus , Streptomyces , Thermomonospora, and Thermobifida  were also detected.
To gain further insight into the diversity of APACMC, the metagenomic dataset was taxonomically profiled at the species level (Fig. 3b). The result showed that Nannocystis exedens (30.44%, marked with the highest bar, see Fig. 3b) was the most prevalent species accompanying with Sphingobacterium sp. 21 (23.33%), Parapedobacter composti (18.97%), Cytophagaceae bacterium SCN 52-12 (13.69%), Thermobacillus composti (18.24%) and Saccharomonospora glauca (11.91%). It was reported that most of them (marked with asterisk, see Fig. 3b) are proficient degraders of lignocellulose [18, 36, 40]. Simultaneously, a number of thermophilic pectinolytic bacteria (marked with solid circle, see Fig. 3b) including Thermobispora bispora , Thermomonospora curvata , and Thermobifida fusca  were found. Collectively, large proportion of lignocellulolytic and pectinolytic microorganisms further confirms that APACMC has possessed the targeted functions.
Functional profiles of predicated genes in APACMC metagenome
Annotations by the MG-RAST pipeline revealed that 99.98% of the predicted genes were protein coding, among which 81.01% had been assigned a putative function. The COG and KEGG repertoire of the predicted genes was analyzed to assess the primary functions of these genes in APACMC. The COG categories analysis (Fig. 4a) showed that APACMC was enriched for amino acid metabolism (8.8% in all COG functional categories), general function (8.7%), inorganic ion metabolism (7.9%), carbohydrate transport and metabolism (7.1%), energy production and conversion (6.6%), and cell wall/membrane/envelope biogenesis (6.1%). The comparative COG analysis of the APACMC with another four well-known lignocellulose-degrading consortia from rain forest compost, switchgrass-adapted compost , Sao Paulo zoo park compost , and rice straw-adapted compost  revealed that they shared similar metabolic patterns, particularly associated with carbohydrates and amino acids transport and metabolism (Additional file 4: Fig. S2). The KEGG ontology exhibited analogous patterns (Fig. 4b), where carbohydrate metabolism (17.0%), amino acid metabolism (16.0%), energy metabolism (9.4%), nucleotide metabolism (6.6%), and membrane transport (6.2%) were abundant. Generally, these observations indicate that APACMC has successfully enriched several desired functional capacities, especially, for carbohydrate metabolism.
In order to get more detailed information about the decomposition of pectin, specific COGs involved in pectin transport and metabolism were further analyzed (Additional file 5: Table S3). APACMC harbored a broad spectrum of genes involved in the metabolism of different monosaccharide building blocks of pectin (e.g., arabinose, fucose, galactose, mannose, rhamnose, xylose, etc.), all of which accounted for 25.8% of the COG subcategory G (Carbohydrate transport and metabolism). Additionally, the genes associated with carbohydrate transporters and phosphotransferase systems were also very plentiful. For example, ABC-type sugar transport system, permease, TonB, and phosphotransferase system, which are responsible for the uptake, transport, and phosphorylation of sugars [18, 44], took up 7.2, 10.6, 2.6, and 1.0% of the COG subcategory G, respectively. In summary, the rich diversity of gene functions in carbohydrate transport and metabolism indicates that APACMC has enriched a great potential for the degradation of pectin.
The diversity, abundance, and phylogenetic distribution of CAZymes in APACMC metagenome
It is well established that the plant biomass-degrading capacities of microbial consortia are closely related to genes encoding CAZymes . To gain an overview of microbial degradation of main polymers in apple pomace, we screened APACMC metagenome for the discovery of microorganisms and genes encoding CAZymes. All the candidate genes of APACMC metagenome were searched against the CAZy database using dbCAN  for the presence of at least one relevant catalytic domain or carbohydrate-binding module, rather than overall sequence similarity to known CAZymes. The results showed that APACMC harbored a total of 9274 different CAZyme genes, which distributed heterogeneously among glycoside hydrolases (GHs, 35.6%), glycosyltransferases (GTs, 26.9%), carbohydrate esterases (CEs, 17.5%), carbohydrate-binding modules (CBMs, 13.3%), auxiliary activities (AAs, 4.6%), and polysaccharide lyases (PLs, 2.1%) (Additional file 6: Table S4).
To link metabolic functions to abundant consortia members, all the CAZyme genes of APACMC were taxonomically classified. The phylogenetic distribution of CAZyme genes showed that the abundance of CAZyme genes varied across bacterial phyla and most of them were derived from Bacteroidetes (51.86%), Proteobacteria (30.82%), Actinobacteria (13.52%), and Firmicutes (2.31%) (Fig. 5a). Notably, members of Bacteroidetes were predominant in the CAZymes classes (GHs, GTs, CEs, PLs, and CBMs) of APACMC, which was remarkably different from those of rice straw-adapted compost (Actinobacteria)  and corn stover-adapted compost (Proteobacteria) .
As the CAZyme genes were unevenly distributed within each phylum, the extensive phylogenetic distributions of CAZyme genes at lower taxonomic levels were further investigated. In addition to Sphingobacterium and Niabella of the phylum Bacteroidetes, CAZyme genes were also abundant in Pseudoxanthomonas and Chelatococcus of the phylum Proteobacteria, Microbacterium of the phylum Actinobacteria, as well as Thermobacillus of the phylum Firmicutes (Fig. 5b). At the species level, six different species of phylum Bacteroidetes, which accounted for 25.06% of the total CAZyme genes, were present in the top-10 richest members harboring CAZyme genes (Fig. 5c). This finding indicates that members of Bacteroidetes possess a much abundant and wider range of CAZyme catalog in APACMC. Besides, the CAZyme genes were also detected in uncharacterized species of Chelatococcus (5.04%), Pseudoxanthomonas sp. GW2 (3.62%), and Pseudoxanthomonas sp. J31 (4.48%) of the phylum Proteobacteria and Thermobacillus composti (1.57%) of the phylum of Firmicutes.
The phylogenetic distributions of CAZyme genes corresponded well to the structure of the ecologically dominant species in APACMC, which confirms the assumption that the functional traits of consortia have a direct correlation with their taxonomic profiles . In conclusion, the CAZymes profile reveals that polysaccharides of apple pomace are decomposed by the predominant Bacteroidetes in cooperation with Proteobacteria, Actinobacteria, and Firmicutes. Together with the COG profiles for glycan degradation, the diverse repertoire of CAZymes provides a basis for a collaborative system tailored to the processing and metabolizing of apple pomace in the compost habitat.
Mining for pectinolytic enzymes
Pectin is the major composition of apple pomace and is an extremely structurally complex polysaccharide, which is constituted of as many as 17 different monosaccharides and more than 20 different linkages . The representative structure of pectin is schematically shown in Fig. 6. It is basically composed of homogalacturonan (HG), rhamnogalacturonan I (RG-I), the substituted galacturonans rhamnogalacturonan II (RG-II), and xylogalacturonan (XGA) . Due to its complex and heterogeneous structure, the efficient and complete degradation of pectin involves a battery of enzymes which act specifically and synergistically. These pectin-degrading enzymes are classified as de-polymerases (hydrolases and lyases), pectinesterases, and de-branching enzymes based on the action mode and site.
The different types of pectinolytic enzymes and their cleavage sites are depicted in Fig. 6. The degradation of pectin is caused by the de-esterification of methoxyl groups, affecting the texture and rigidity of the cell wall . Pectin methylesterases (PMEs, EC 188.8.131.52) remove the methyl groups from the HG backbone to give access to de-polymerases, while pectin acetylesterases (PAEs, EC 184.108.40.206) remove acetyl groups from acetylated HG and RG . Hydrolases (polygalacturonases PGs, EC 220.127.116.11, 67 and 82) and lyases (pectin lyases: PNLs, EC 18.104.22.168 and pectate lyases: PELs, EC 22.214.171.124 and 9) preferentially degrade the α-1,4-glycosidic bonds of HG/XGA backbones by hydrolysis and β-elimination, respectively. Similarly, in the initial deconstruction of pectin, PNLs also play an essential role since it is the only enzyme that can cleave the α-1, 4 bonds of highly esterified pectin without prior actions of other enzymes. De-branching enzymes are responsible for the cleavage of the backbone or lateral chains of RG-I and RG-II . Rhamnogalacturonases (EC 126.96.36.199, 173 and 174) and rhamnogalacturonan lyases (RGLs: EC 188.8.131.52 and 24) are involved in the cleavage of the RG-I backbone; α-l-rhamnosidases (EC 184.108.40.206) act on hydrolytic cleavage of the RG chain at non-reducing end releasing rhamnose; arabinofuranosidases (EC 220.127.116.11) attack on α-l-arabinofuranosides, α-l-arabinans, arabinoxylans, and arabinogalactans; arabinogalactanases (EC 18.104.22.168) and β-galactosidases (EC 22.214.171.124) act randomly on the galactan core of AGs; and feruloyl esterases (EC 126.96.36.199) are needed to release the ferulic acid attached to C-2 of arabinose and C-6 of galactose.
According to the results of CAZymes annotation, we found an extremely abundant of genes associated with the complete degradation of pectin (Additional file 7: Table S5). A total of 1756 entries were identified as encoding pectinolytic enzymes, which took up to 18.93% of the total CAZy genes. As summarized in Table 1, these entries contained 105 PLs candidates from 6 families, 881 GHs mainly from 17 families, 537 CEs from 5 families, and 233 CBMs from 5 families. Compared to another two compost habitats (Table 1), the catalog of pectinolytic enzymes in APACMC was much more abundant and diverse than RSA (843 candidates from 31 families) and EMSD5 (398 candidates from 28 families), indicating that APACMC has a better potential for pectin degradation based on CAZyme inventory. As shown in Table 1, a large panel of pectinolytic enzymes was found, such as PGs, PNLs, PELs, RGLs, PMEs, PAEs, α-l-rhamnosidases, arabinofuranosidases, arabinogalactanases, and β-galactosidases. Furthermore, we also detected a wealth of CBMs which possibly associate to pectin degradation. For example, some members of family CBM32 have been found to bind oligogalacturonides to counteract the loss of binding affinity between thermophilic pectinases and their substrates at elevated temperature [1, 52]. These findings indicate that APACMC exhibits a collaborative enzymatic system efficient in the complete degradation of pectin.
To assess the identity of these possible pectinolytic enzymes with known proteins, these amino acid sequences of 1756 putative pectinolytic genes were searched against NCBI non-redundant (NCBI-NR), CAZy, NCBI environmental (NCBI-ENV), and Swiss-Prot databases by DIAMOND  (Fig. 7; Additional file 8: Table S6). Firstly, the results based on NCBI-NR showed that the amino acid sequence identity of these 1756 genes ranged from 25 to 100%, with an average of 76.95% (Additional file 8: Table S6). And 23.29% of these sequences were most similar to proteins annotated as “hypothetical/predicted protein” or “proteins of unknown function” in NCBI-NR. Secondly, only 9.83% of these sequences were highly similar (>95% sequence identity) to any CAZy database entry, indicating that most of these sequences had not been previously deposited in CAZy . And 129 sequences were considered as novel with less than 30% identity . Thirdly, only 1.92% of these putative pectinolytic genes are more than 75% identical to sequences deposited in the NCBI-ENV database, demonstrating that these enzymes also have not been observed in previous metagenome projects . Lastly, 145 sequences had less than 30% identity to any known proteins deposited in Swiss-Prot, indicating that their assigned activity has not been verified biochemically. Summarily, the large amount of relatively low identity sequences indicate that the strategy we adopted has great potential in mining novel enzymes from environmental sources.
Mining for pectinolytic microbes
To explore the phylogenetic origins of these pectinolytic enzymes, we examined the top BLASTX hit organism of each identified enzyme at species level, deciphered the role of individual pectinolytic microbe and their potential synergistic action in the process of pectin degradation. Of the 1756 sequences encoding pectinolytic enzymes, most of their phylogenetic affiliations predicted by BLASTX were consistent with the predicted source organisms of APACMC metagenomic bins (Additional file 7: Table S5). Many of these genes are homologous to those found in the top-10 abundant community members (Fig. 5c), such as Sphingobacterium sp. 21, Sphingobacterium spiritivorum, Thermobacillus composti, and Cytophagaceae bacterium SCN 52-12, which further verified that APACMC was successfully target-enriched for pectin degradation.
To provide a systematic overview of pectin degradation by individual member of APACMC, the specific taxonomic assignments of key pectinolytic enzymes, i.e., PGs, PMEs, PNLs, PELs, α-l-rhamnosidases, and arabinofuranosidases, were illustrated in Fig. 6. Clearly, pectinolytic species in APACMC were considerably diverse. The majority of candidate PGs was mainly originated from a variety of Bacteroidetes species, which consisted of Parapedobacter composti, Sphingobacterium sp. 21, Parapedobacter indicus, Sphingobacterium thalpophilum, as well as Opitutus terrae and Verrucomicrobia bacterium IMCC26134 from the Verrucomicrobia. It is generally known that the bacterial sources of PELs and PNLs are some specific bacteria such as Bacillus sp. and Pseudomonas sp. . However, in our study, a broad range of other bacteria were the major producers, including Parapedobacter composti, Sphingobacterium sp. 21 and Sphingobacterium thalpophilum, Opitutaceae bacterium TSB47, Thermobifida fusca TM 51, Saccharomonospora glauca, Nannocystis exedens, and Mesorhizobium sp. LC 103. The taxonomic classifications of genes encoding PMEs and PAEs revealed that they were dominantly from Sphingobacterium sp. 21, Parapedobacter composti, Sphingobacterium thalpophilum, Cytophagaceae bacterium SCN 52-12 of Bacteroidetes, Pseudoxanthomonas sp. J35, Cellvibrio mixtus, Nannocystis exedens of Proteobacteria, Saccharomonospora glauca, and Ruania albidiflava of Actinobacteria along with Thermobacillus composti of Firmicutes. The sequences encoding α-l-rhamnosidases and arabinofuranosidases were predominantly affiliated with organisms from Parapedobacter composti, Sphingobacterium sp. 21 and Parapedobacter indicus of Bacteroidetes, Streptomyces caeruleatus and Ruania albidiflava of Actinobacteria, and Thermobacillus composti of Firmicutes. These collective data indicate a synergistic action of multiple members derived from Bacteroidetes, Proteobacteria, Actinobacteria, Firmicutes, and Verrucomicrobia in the degradation of pectin.
Most of the currently available pectinolytic enzymes are reported from filamentous fungal species (e.g., Aspergillus sp. and Penicillium sp.). However, as the wide functional diversity, broad array of terminal electron acceptors, high ability to degrade lignin , as well as more amenable to genetic manipulation, pectinolytic bacteria are likely to play important roles in future biotechnology strategies. Surprisingly, a variety of bacteria were identified as the major producers of various pectinolytic enzymes, such as Cellvibrio mixtus, Cytophagaceae bacterium SCN 52-12, Nannocystis exedens, Opitutaceae bacterium TSB47, Parapedobacter composti, Parapedobacter indicus, Ruania albidiflava, Saccharomonospora glauca, Sphingobacterium sp. 21, Sphingobacterium thalpophilum, and Thermobacillus composti. Many of these species were initially described to degrade pectin. Strikingly, Parapedobacter composti, Sphingobacterium sp. 21, and Sphingobacterium thalpophilum, which were identified as the top-10 richest members, each of them harbors a great number of genes encoding various pectinolytic enzymes, indicating that they are well equipped with systematic pectinolytic enzymes. This make them to be promising bacterial sources of pectinolytic enzymes and potential efficient degraders of pectin in the future. Although putative pectinolytic enzymes have been annotated in the genomes of these type-strains , so far the information available on the pectinolytic enzymes of these strains is very limited. Our data provide insight into their potentials and highlight their importance in the complex degradation of pectin.
Novel pectinolytic microbes and enzymes have potential application in numerous industrial processes. Here, we adopted a strategy which combined metagenomics sequencing with enrichment culture technique to rapidly discover efficient pectin degraders and novel pectinolytic enzyme sequences. The immensely diverse pectinolytic microbes and enzymes found in our study will not only shed light on the current understanding of microbial interaction and enzymatic synergism in pectin degradation, but also expand the arsenal of proficient degraders and enzymes for lignocellulosic biofuel production. When combined with high-throughput strategies, such as cell-free protein expression system, droplet-based microfluidics, fluorescence-activated cell sorting (FCAS), and nanostructure-initiator mass spectrometry (NIMS), the efficiency of this strategy for obtaining novel enzymes may meet the ever-growing demand from various industries.
Enrichment of apple pomace-adapted microbial community in compost habitat
The composting materials were composed of apple pomace and fresh cattle manure. Apple pomace was kindly provided by Shaanxi Haisheng fresh fruit juice Co. Ltd., China . The cattle manure was collected from the Northwest A&F University farm located in Yangling, China. The composting experiment was conducted from September 15 to October 29, 2016 according to Sun et al.  with slight modifications. Briefly, the cattle manure was mixed with apple pomace to adjust the C/N ratio to 30:1 and the moisture content to around 60% and then the mixture was placed in rectangular foam containers as described in Sun et al. . The piles of compost were turned and sampled daily. Samples were pooled from the top, middle, and bottom of the composting, and then mixed completely. When the pile temperature dropped to ambient temperature at the end of the maturing stage, the composting process was considered as completed. The each sample was split into two parts: one part was stored at 4 °C for subsequent physicochemical analysis and the other was stored at −80 °C for high-throughput sequencing.
The pile temperature was monitored every 24 h by inserting a mercury thermometer in the center of the composting material. The moisture content was measured gravimetrically after drying samples at 105 °C for 24 h. The pH values of the samples were tested in water (solid-to-water ratio of 1:10, w/v) with a pH meter . The contents of pectin were determined by modified carbazole method .
16S rDNA sequencing and phylogenetic classification
The 16S rDNA sequencing was performed at the Frasergen Genoimcs Institute (Wuhan, China) using the Illumina MiSeq platform . The 16S V3-V4 region was amplified using the primers 338F and 806R. After the processing of raw data, the high-quality sequences were subjected to filter singletons, remove chimeras, and cluster into operational taxonomic units (OTUs) at a 97% identity using UPARSE . A representative sequence of each OTU was assigned to a taxonomic level in the SILVA database  using the RDP classifier. Microbial diversity and richness measurements were performed using MOTHUR . The microbial diversity was estimated by Shannon and Simpson, and the richness was determined by Chao and Ace estimators.
Metagenome sequencing, de novo assembly, and Open Reading Frames (ORFs) prediction
Metagenome sequencing, de novo assembly, and ORFs prediction were performed by Frasergen Genoimcs Institute (Wuhan, China) according to Qin with slight modifications . Briefly, a library with 400-bp clone insert size was constructed and sequenced on Illumina HiSeq4000 platform. Sequence reads were quality trimmed to an accuracy of 98.0% and duplicate reads were identified and removed prior to assembly. Nearly 89.6 million high-quality reads were generated (16.2 Gb). High-quality short reads of the DNA sample were assembled by the SOAPdenovo assembler  with a k-mer length of 39–47. The assembled contigs longer than 500 bp were subject to ORFs prediction using the MetaGene  with default parameters. The redundant ORFs were removed by CD-HIT  from the non-redundant gene catalog and the abundances were annotated by SOAPaligner .
Taxonomic assignment and functional classification in metagenomic database
Taxonomic annotation of predicated genes was performed by BLASTP against the NCBI-NR database with an E value of 1e−5. The APACMC metagenomic dataset was also taxonomically profiled at species level by MetaPhlAn2 , MG-RAST  and MEGAN 6 . The phylogenetic tree was generated using iTOL software . Functional classification of predicted gene was performed by BLASTP against eggNOG database and by KOBAS 2.0 (a Orthology Based Annotation System)  against KEGG database. The “function comparison” module of integrated microbial genomes with microbiome samples (IMG/M)  were applied to compare the COG category of APACMC against another four well-known lignocelluloses-degrading microbiomes available on IMG/M, including rain forest compost (IMG Submission ID 5968), switchgrass-adapted compost , sao paulo zoo park compost , and rice straw-adapted compost .
Carbohydrate-active enzymes (CAZymes): annotation and phylogenetic analysis
Searches for CAZymes were performed as described by Wang and coworkers . Briefly, the amino acid sequences of the predicted ORFs in the APACMC metagenome were annotated by dbCAN, an automated CAZyme signature domain-based annotation method based on family-specific HMMs  by MAFFT and HMMER. After identification, these sequences were searched against NCBI non-redundant (NCBI-NR), CAZy database, NCBI environmental database (NCBI-ENV), and Swiss-Prot database by DIAMOND  with a cutoff of E value <1e−5. The phylogenetic distributions in the top ten abundant members possessing CAZymes were visualized via software Circos  at the level of phylum, genus, and specie.
Specific pectin-degrading genes: annotation and phylogenetic analysis
The predicted sequences encoding pectinolytic enzymes were re-annotated and verified using DIAMOND using a sensitive setting  against the proteins deposited in NCBI-NR database. The phylogenetic origins of candidate genes were determined by MEGAN 6 .
Sequence data submission
The assembled metagenome datasets were submitted to IMG/M and Metagenomics RAST server (MG-RAST) under the project ID 117466 and mgs566360, respectively.
apple pomace-adapted compost
apple pomace-adapted compost microbial community
clusters of orthologous groups
open reading frames
operational taxonomic unit
galacturonans rhamnogalacturonan II
Abbott DW, Boraston AB. Structural biology of pectin degradation by Enterobacteriaceae. Microbiol Mol Biol Rev. 2008;72(2):301–16.
Allgaier M, Reddy A, Park JI, Ivanova N, D’Haeseleer P, Lowry S, Sapra R, Hazen TC, Simmons BA, VanderGheynst JS, Hugenholtz P. Targeted discovery of glycoside hydrolases from a switchgrass-adapted compost community. PLoS ONE. 2010;5(1):e8812.
Biswal AK, Hao Z, Pattathil S, Yang X, Winkeler K, Collins C, Mohanty SS, Richardson EA, Gelineo-Albersheim I, Hunt K, Ryno D, Sykes RW, Turner GB, Ziebell A, Gjersing E, Lukowitz W, Davis MF, Decker SR, Hahn MG, Mohnen D. Downregulation of GAUT12 in Populus deltoides by RNA silencing results in reduced recalcitrance, increased growth and reduced xylan and pectin in a woody biofuel feedstock. Biotechnol Biofuels. 2015;8:41.
Biswal AK, Soeno K, Gandla ML, Immerzeel P, Pattathil S, Lucenius J, Serimaa R, Hahn MG, Moritz T, Jönsson LJ. Aspen pectate lyase Ptxt PL1-27 mobilizes matrix polysaccharides from woody tissues and improves saccharification yield. Biotechnol Biofuels. 2014;7(1):1.
Bonnin E, Garnier C, Ralet MC. Pectin-modifying enzymes and pectin-derived materials: applications and impacts. Appl Microbiol Biotechnol. 2014;98(2):519–32.
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.
Chen IA, Markowitz VM, Chu K, Palaniappan K, Szeto E, Pillay M, Ratner A, Huang J, Andersen E, Huntemann M, Varghese N, Hadjithomas M, Tennessen K, Nielsen T, Ivanova NN, Kyrpides NC. IMG/M: integrated genome and metagenome comparative data analysis system. Nucleic Acids Res. 2017;45(D1):D507–16.
Chung Daehwan, Pattathil Sivakumar, Biswal Ajaya K, Hahn Michael G, Mohnen Debra, Westpheling J. Deletion of a gene cluster encoding pectin degrading enzymes in Caldicellulosiruptor bescii reveals an important role for pectin in plant biomass recalcitrance. Biotechnol Biofuels. 2014;7(1):0.
DeAngelis KM, Gladden JM, Allgaier M, D’haeseleer P, Fortney JL, Reddy A, Hugenholtz P, Singer SW, Vander Gheynst JS, Silver WL, Simmons BA, Hazen TC. Strategies for Enhancing the Effectiveness of Metagenomic-based Enzyme Discovery in Lignocellulolytic Microbial Communities. Bioenergy Research. 2010;3(2):146–58.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.
Edwards MC, Doran-Peterson J. Pectin-rich biomass as feedstock for fuel ethanol production. Appl Microbiol Biotechnol. 2012;95(3):565–75.
Francocci F, Bastianelli E, Lionetti V, Ferrari S, De Lorenzo G, Bellincampi D, Cervone F. Analysis of pectin mutants and natural accessions of Arabidopsis highlights the impact of de-methyl-esterified homogalacturonan on tissue saccharification. Biotechnol Biofuels. 2013;6(1):163.
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.
Hahnke RL, Meier-Kolthoff JP, Garcia-Lopez M, Mukherjee S, Huntemann M, Ivanova NN, Woyke T, Kyrpides NC, Klenk HP, Goker M. Genome-Based Taxonomic Classification of Bacteroidetes. Front Microbiol. 2016;7:2003.
Hess M, Sczyrba A, Egan R, Kim T-W, Chokhawala H, Schroth G, Luo S, Clark DS, Chen F, Zhang T. Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science. 2011;331(6016):463–7.
Himmel ME, Ding SY, Johnson DK, Adney WS, Nimlos MR, Brady JW, Foust TD. Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science. 2007;315(5813):804–7.
Huson DH, Beier S, Flade I, Gorska A, El-Hadidi M, Mitra S, Ruscheweyh HJ, Tappu R. MEGAN community edition—interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comput Biol. 2016;12(6):e1004957.
Jimenez DJ, Chaves-Moreno D, van Elsas JD. Unveiling the metabolic potential of two soil-derived microbial consortia selected on wheat straw. Sci Rep. 2015;5:13845.
Khan M, Nakkeeran E, Umesh-Kumar S. Potential application of pectinase in developing functional foods. Annu Rev Food Sci Technol. 2013;4:21–34.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.
Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, Li S, Yang H, Wang J, Wang J. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010;20(2):265–72.
Lionetti V, Francocci F, Ferrari S, Volpi C, Bellincampi D, Galletti R, D’Ovidio R, De Lorenzo G, Cervone F. Engineering the cell wall by reducing de-methyl-esterified homogalacturonan improves saccharification of plant tissues for bioconversion. Proc Natl Acad Sci. 2010;107(2):616–21.
Martins LF, Antunes LP, Pascon RC, de Oliveira JC, Digiampietri LA, Barbosa D, Peixoto BM, Vallim MA, Viana-Niero C, Ostroski EH, Telles GP, Dias Z, da Cruz JB, Juliano L, Verjovski-Almeida S, da Silva AM, Setubal JC. Metagenomic analysis of a tropical composting operation at the sao paulo zoo park reveals diversity of biomass degradation functions and organisms. PLoS ONE. 2013;8(4):e61928.
Mei Y, Chen Y, Zhai R, Liu Y. Cloning, purification and biochemical properties of a thermostable pectinase from Bacillus halodurans M29. J Mol Catal B Enzym. 2013;94:77–81.
Menon V, Rao M. Trends in bioconversion of lignocellulose: biofuels, platform chemicals & biorefinery concept. Prog Energy Combust Sci. 2012;38(4):522–50.
Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A. The metagenomics RAST server–a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinform. 2008;9(1):386.
Montella S, Amore A, Faraco V. Metagenomics for the development of new biocatalysts to advance lignocellulose saccharification for bioeconomic development. Crit Rev Biotechnol 2015: 1–12.
Noguchi H, Park J, Takagi T. MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res. 2006;34(19):5623–30.
Parmar I, Rupasinghe HP. Bio-conversion of apple pomace into ethanol and acetic acid: enzymatic hydrolysis and fermentation. Bioresour Technol. 2013;130(2013):613–20.
Pearson WR. An introduction to sequence similarity (“homology”) searching. Curr Protoc Bioinform. 2013:3.1.1–3.1.8.
Pelloux J, Rusterucci C, J ME. New insights into pectin methylesterase structure and function. Trends Plant Sci. 2007;12(6):267–77.
Qian X, Sun W, Gu J, Wang XJ, Zhang YJ, Duan ML, Li HC, Zhang RR. Reducing antibiotic resistance genes, integrons, and pathogens in dairy manure by continuous thermophilic composting. Bioresour Technol. 2016;220:425–32.
Qin J, Li R, Raes J, Et A. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464(7285):59–65.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–5966.
Reddy AP, Simmons CW, D’haeseleer P, Khudyakov J, Burd H, Hadi M, Simmons BA, Singer SW, Thelen MP, VanderGheynst JS. Discovery of microorganisms and enzymes involved in high-solids decomposition of rice straw using metagenomic analyses. PLoS ONE. 2013;8(10):e77985.
Rosnow JJ, Anderson LN, Nair RN, Baker ES, Wright AT. Profiling microbial lignocellulose degradation and utilization by emergent omics technologies. Crit Rev Biotechnol. 2016:1–15.
Saoudi B, Habbeche A, Kerouaz B, Haberra S, Ben Romdhane Z, Tichati L, Boudelaa M, Belghith H, Gargouri A, Ladjama A. Purification and characterization of a new thermoalkaliphilic pectate lyase from Actinomadura keratinilytica Cpt20. Process Biochem. 2015;50(12):2259–66.
Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS ONE. 2011;6(12):e27310.
Simmons CW, Reddy AP, D’haeseleer P, Khudyakov J, Billis K, Pati A, Simmons BA, Singer SW, Thelen MP, VanderGheynst JS. Metatranscriptomic analysis of lignocellulolytic microbial communities involved in high-solids decomposition of rice straw. Biotechnol Biofuels. 2014;7(1):495.
Sun J, Qian X, Gu J, Wang X, Gao H. Effects of oxytetracycline on the abundance and community structure of nitrogen-fixing bacteria during cattle manure composting. Bioresour Technol. 2016;216:801–7.
Truong DT, Franzosa EA, Tickle TL, Scholz M, Weingart G, Pasolli E, Tett A, Huttenhower C, Segata N. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods. 2015;12(10):902–3.
Valladares Juarez AG, Dreyer J, Gopel PK, Koschke N, Frank D, Markl H, Muller R. Characterisation of a new thermoalkaliphilic bacterium for the production of high-quality hemp fibres, Geobacillus thermoglucosidasius strain PB94A. Appl Microbiol Biotechnol. 2009;83(3):521–7.
Wang C, Dong D, Wang H, Muller K, Qin Y, Wang H, Wu W. Metagenomic analysis of microbial consortia enriched from compost: new insights into the role of Actinobacteria in lignocellulose decomposition. Biotechnol Biofuels. 2016;9:22.
Wang W, Yang H, Zhang Y, Xu J. IoT-enabled real-time energy efficiency optimisation method for energy-intensive manufacturing enterprises. Int J Computer Inter Manuf. 2017;1–18.
Wang X, Chen Q, Lü X. Pectin extracted from apple pomace and citrus peel by subcritical water. Food Hydrocolloids. 2014;38:129–37.
Xiao C, Anderson CT. Roles of pectin in biomass yield and processing for biofuels. Front Plant Sci. 2013;4:67.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.
Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40:W445–51.
Yuan P, Meng K, Shi P, Luo H, Huang H, Tu T, Yang P, Yao B. An alkaline-active and alkali-stable pectate lyase from Streptomyces sp. S27 with potential in textile industry. J Ind Microbiol Biotechnol. 2012;39(6):909–15.
Zhang Y, Li H, Gu J, Qian X, Yin Y, Li Y, Zhang R, Wang X. Effects of adding different surfactants on antibiotic resistance genes and intI1 during chicken manure composting. Bioresour Technol. 2016;219:545–51.
Zhu N, Yang J, Ji L, Liu J, Yang Y, Yuan H. Metagenomic and metaproteomic analyses of a corn stover-adapted microbial consortium EMSD5 reveal its taxonomic and enzymatic basis for degrading lignocellulose. Biotechnol Biofuels. 2016;9:243.
Zykwinska A, Thibault JF, Ralet MC. Organization of pectic arabinan and galactan side chains in association with cellulose microfibrils in primary cell walls and related models envisaged. J Exp Bot. 2007;58(7):1795–802.
MZ conceived and designed the experiments. MZ performed the majority of experimental work, analyzed the results, and drafted this manuscript. PG participated in the analysis of bioinformatics. TW, LG, and HY assisted in the enrichment of APACMC and the chemical composition test. JG and CC helped to guide the composting process and process bioinformatics, respectively. MZ performed all data analyses. XL led and coordinated the overall project. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of supporting data
All data are included in the manuscript and additional information will be provided upon request.
Consent for publication
Not applicable. All the authors agree to the publication of this manuscript.
Ethics approval and consent to participate
This study was supported by the Special Fund for Agro-scientific Research in the Public Interest (No. 201503135), Ministry of Agriculture, China.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1. Red agar test.tif. The effectiveness of enrichment cultures (post- and pre-) by grown on ruthenium red agar plates.
Additional file 2: Table S1. Diversity and OUT distributions.xlsx. Diversity and OTU distribution of CM0, Mes5, The15 and Mat30.
Additional file 3: Table S2. De novo assembly results.docx. Illumnia reads and de novo assembly results of APACMC metagenome.
Additional file 4: Figure S2. The comparison of COG category.tiff. The COG comparison of APAMC with other well-known lignocellulosic metagenomes.
Additional file 8: Table S6. Distributions of identities.xlsx. Distributions of identities of predicated pectinolytic enzyme sequences.
About this article
Cite this article
Zhou, M., Guo, P., Wang, T. et al. Metagenomic mining pectinolytic microbes and enzymes from an apple pomace-adapted compost microbial community. Biotechnol Biofuels 10, 198 (2017). https://doi.org/10.1186/s13068-017-0885-y