- 1College of Environmental and Chemical Engineering, Heilongjiang University of Science and Technology, Harbin, China
- 2State Key Lab of Urban Water Resource and Environment, Harbin Institute of Technology, Harbin, China
The cellulosome of Ruminiclostridium thermocellum was one of the most efficient cellulase systems in nature. However, the product of cellulose degradation by R. thermocellum is cellobiose, which leads to the feedback inhibition of cellulosome, and it limits the R. thermocellum application in the field of cellulosic biomass consolidated bioprocessing (CBP) industry. In a previous study, R. thermocellum M3, which can hydrolyze cellulosic feedstocks into monosaccharides, was isolated from horse manure. In this study, the complete genome of R. thermocellum M3 was sequenced and assembled. The genome of R. thermocellum M3 was compared with the other R. thermocellum to reveal the mechanism of cellulosic saccharification by R. thermocellum M3. In addition, we predicted the key genes for the elimination of feedback inhibition of cellobiose in R. thermocellum. The results indicated that the whole genome sequence of R. thermocellum M3 consisted of 3.6 Mb of chromosomes with a 38.9% of GC%. To be specific, eight gene islands and 271 carbohydrate-active enzyme-encoded proteins were detected. Moreover, the results of gene function annotation showed that 2,071, 2,120, and 1,246 genes were annotated into the Clusters of Orthologous Groups (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively, and most of the genes were involved in carbohydrate metabolism and enzymatic catalysis. Different from other R. thermocellum, strain M3 has three proteins related to β-glucosidase, and the cellobiose hydrolysis was enhanced by the synergy of gene BglA and BglX. Meanwhile, the GH42 family, CBM36 family, and AA8 family might participate in cellobiose degradation.
1. Introduction
Lignocellulosic biomass was considered an ideal sustainable resource for potential feedstock for various value-added chemicals and biofuels (Khare et al., 2015; Usmani et al., 2021). As a result, the rational use of lignocellulosic biomass not only relieves the pressure of fossil energy shortage but also mitigates natural environment damage caused by improper treatment (Staples et al., 2017). Converting lignocellulosic feedstocks into high-value chemical products includes three steps: pretreatment, saccharification, and fermentation (Yadav et al., 2020). Saccharification refers to the hydrolysis of holocellulose into monosaccharides/oligosaccharides, which facilitates the downstream process. The saccharification of lignocellulose feedstocks was one of the bottlenecks of lignocellulosic biomass utilization (Guo et al., 2018; Usmani et al., 2020). Cellulosic biomass can be saccharified by acids or cellulases. Cellulose hydrolysis by cellulase is considered more environmentally friendly than cellulose hydrolysis by acid when conducted under mild conditions. The degradation of cellulose is achieved by the synergistic action of endoglucosidase, extranosidase, and β-glucosidase (Sheng et al., 2016). Both bacteria and fungi can synthesize cellulase in nature; Trichoderma and Aspergillus sp. are known for their potential to produce cellulases, while these fungi lack a complete cellulase system, which leads to a decrease in the catalytic efficiency of cellulase (Srivastava et al., 2018). Anaerobic bacteria degrade cellulose by synthesizing cellulosomes, which gather different cellulases in a narrow space and anchor them on the cell surface, the cellulosome of R. thermocellum is the most efficient cellulase system found at present (Mazzoli and Olson, 2020).
However, the hydrolysis of lignocellulosic biomass by R. thermocellum has not been industrialized for the low activity of β-glucosidase, which is insufficient for the lignocellulosic feedstocks saccharification (Sheng et al., 2016). Moreover, the feedback inhibition of exoglucanosidases induced by the accumulation of cellobiose severely reduced the catalytic efficiency of the cellulosome (Lamed et al., 1985; Tian et al., 2016; Haldar and Purkait, 2020). The addition of exogenous β-glucosidase is a way to improve the hydrolysis efficiency of cellulosome, but additional β-glucosidase directly increases the cost and complexity of the saccharification process. Accordingly, enhancing the activity of β-glucosidase activity of wild R. thermocellum by building recombinant strains that secrete large amounts of β-glucosidase was a promising solution. Meki et al. fused the E. coli plasmid containing cloned bglA into wild R. thermocellum ATCC 27405 to construct a recombinant strain R. thermocellum ATCC 27405 (+McbglA), the result indicated that the β-glucosidase activity expressed by recombinant strain was 2.3 times higher than that of the wild strain at the late logarithmic growth stage (Maki et al., 2013). Waeonukul et al. found that the addition of BglB from R. thermocellum S14 to cellulosome was observed to increase the saccharification rate of cellulose compared to that of Novazyme-188 and cellulosome alone (Waeonukul et al., 2012). Nevertheless, the activity and stability of β-glucosidase secreted by recombinant strain decreased during the hydrolysis process (Zhang et al., 2017). Isiam et al. found that when cellobiose was used as a carbon source, seven cellulosome structural proteins, 31 cellulosome-related glycosidases, and 19 non-cellulosome glycoside hydrolases were expressed in R. thermocellum, which suggests that the degradation of cellobiose by R. thermocellum was not only related to the β-glucosidase gene but also associated with genes other than β-glucosidase (Islam et al., 2006). Therefore, it is particularly important to find strains with a stable ability to degrade cellobiose and reveal the genes involved in stable cellobiose degradation by R. thermocellum.
In previous studies, we isolated an R. thermocellum M3 that can efficiently degrade lignocellulosic biomass from horse manure. Different from other R. thermocellum, 97% of the cellulosic saccharification products of R. thermocellum M3 were monosaccharides. More importantly, R. thermocellum M3 inherited the ability of cellobiose degradation stably that conducts R. thermocellum M3 being an excellent sample for stable expression of exogenous β-glucosidase in the genus of R. thermocellum. In this study, we report the whole genome sequence of R. thermocellum M3 and compare the high-quality complete genome sequence of R. thermocellum M3 with both intra- and inter-generically to those of its close or distant phylogenetic relatives. Moreover, genes related to cellobiose degradation were comprehensively analyzed.
2. Materials and methods
2.1. Bacterial strain and cultivation
Ruminiclostridium thermocellum M3 strain was isolated and enriched from horse manure by Sheng et al. (2016) and deposited in the Microbiology Laboratory of the School of Environment and Chemical Engineering, Heilongjiang University of Science and Technology. The seed was stored in a constant temperature incubator at −20°C and cultivated in an anaerobic bottle (filled with nitrogen) containing modified ATCC 1191 (MA) medium. The main components of culture medium were K2HPO4, 1.5 g/L; MgSO4·7H2O, 0.2 g/L; (NH4)2SO4, 1.0 g KCl, 0.2 g/L; L-cysteine, 0.5 g/L; KH2PO4, 3.0 g/L; CaCl2•2H2O, 0.025 g/L; NaCl, 1.0 g/L; Yeast, 1.5 g/L; and Avicel, 5.0 g/L. The temperature of the culture was maintained at 60°C at 120 rpm.
2.2. DNA extraction and whole genome sequencing
The R. thermocellum samples for whole genome sequencing analysis were cultured in MA medium for 24 h at 60°C, centrifuged for 5 min at 4°C, 12,000 × g, then the cell pellet was washed twice with normal saline (NS), and crushed with the FastPrep-24 instrument in lysing matrix B tubes (MP Biomedical) for 40 s to release the genomic DNA from cells. The extraction of genomic DNA was performed with the Tiangen bacterial DNA mini kit (Tiangen Biotech Co. Ltd., Beijing, China) according to the manufacturer’s protocol. The harvested DNA was detected using agarose gel electrophoresis and then quantified by Qubit 4.0 (ThermoFisher, Q33226). The genomic sequencing was conducted by the Pacbio sequencing platform with de novo assembly (SMRT portal; Berlin et al., 2015).
2.3. Genome annotation and component prediction
The genome annotation and gene function were predicted in the Gene Ontology (GO) database, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, the Clusters of Orthologous Groups (COG) database, and the Non-Redundant Protein databases (NR). The whole genome BLAST search (E-value below 1e−5, minimal alignment length percentage above 40%) was performed with the above four databases. The prediction of carbohydrate-active enzymes was conducted with the Carbohydrate-Active enZYmes Database (Lombard et al., 2013).
Genome component prediction of R. thermocellum M3 strains included the coding gene, repetitive sequences, signal peptide, genomic islands, prophage, lipoprotein, prophage, and clustered regularly interspaced short palindromic repeat sequences (CRISPR). The GeneMarkS program was conducted to retrieve the coding genes (http://topaz.gatech.edu/; Yang et al., 2020). The tRNA was predicted by the Aragorn program (Laslett and Canback, 2004), the rRNA was predicted by the RNAmmer program (Lagesen et al., 2007), and the miscRNA was predicted by Infernal (v1.1.2; Mostajo Berrospi et al., 2019). Repeat Modeler was used to predict the repeat sequence de novo of the assembly results, and the RepeatMasker program was used for identifying repetitive elements in nucleotide sequences.
Genomic islands were predicted by the Island Path-DIOMB program (Hsiao et al., 2003). The CRISPR was predicted by the CRISPR recognition tool (CRT; Grissa et al., 2007). Prophages were predicted using PhiSpy (Hähnke et al., 2009). RepeatMasker was used to identify the location and frequency of repeats on the genome (Saha et al., 2008). NCBI Blast+ was used to compare the protein sequences with CDD, KOG, COG, NR, NT, PFAM, Swissprot, TrEMBL, and other databases to obtain the functional annotation information.
2.4. Whole genome-based comparative genomic analysis
The core genes, specific genes, the gene family phylogenetic tree, single nucleotide polymorphism (SNP), and genome visualization were analyzed to reveal the result of comparative genomics (Zhong et al., 2018). The genome sequences of R. thermocellum DSM 2360, R. thermocellum ATCC 27405, R. thermocellum DSM 1313, and R. thermocellum AD2 were obtained from the NCBI database. Genomic alignments among R. thermocellum M3 and other genomes of R. thermocellum were performed using the MUMmer and LASTZ tools (Kurtz et al., 2004). Core genes and specific genes were analyzed using the CDHIT rapid clustering of similar proteins software with a threshold of 50% pairwise identity and a 0.7 length difference cutoff for amino acids (Li et al., 2001, 2002; Li and Godzik, 2006). The relationships between five R. thermocellum strains were analyzed, and the results were represented in a Venn diagram. NCBI Blast+ was used to compare the predicted 16S rRNA sequence with the NCBI 16S database to obtain its homologous strain information and a phylogenetic tree was constructed using mega software (Hall, 2013).
3. Results
3.1. Feature of the whole genome of Ruminiclostridium thermocellum M3 genome
The whole genome of R. thermocellum M3 was sequenced and analyzed with regard to the predictions of coding genes. The identified total size of the genome R. thermocellum M3 was 3,602,270 bp with 39% GC content using SPAdes, and the number of coding genes was 3,195 with an average gene length of 973.56 bp (Supplementary Figure S1). The major characteristics of the genome of the R. thermocellum M3 and four R. thermocellum are summarized to acquire more perceptions concerning the genetic information (Table 1). The final assemblies indicated that the genomes of five R. thermocellum were similar in size and G + C contents, genome annotation yielded 3,062, 3,196, 2,949, 2,959, and 3,077 genes for strains DSM 2360, ATCC 27405, DSM 1313, AD2, and M3, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-t001.jpg)
Table 1. General features of Ruminiclostridium thermocellum M3 genome and comparison with other closely related species.
The result of genome annotation indicated that a series of genes encoding virulence factors (VFDB), antibiotic resistance (CARD), and pathogen-host interactions (PHI-base) were present in R. thermocellum (Supplementary Table S1). Meanwhile, R. thermocellum strains had several genes related to carbohydrate degradation, biosynthesis, and modification enzymes. The whole-genome visualization map clearly identified the composition, location, and function of the R. thermocellum M3 genome, which is demonstrated in Figure 1, with 100% coverage of sequencing. Detailed information includes analysis of CDS, non-coding RNA, COG, functional classification of genes, and the size of genome and GC%. Additionally, among the 3,077 predicted genes of R. thermocellum M3, only 2,071 CDSs were assigned any COG, which accounted for 67.3% of the predicted genes, and an additional 171 CDSs were assigned to group S with an unknown function. There are 137 genes involved in amino acid production and transformation (group E) in the aligned COG database genes, accounting for 6.62% of the COG annotations; the number of carbohydrate transport and metabolism-related genes (group G) was 130, accounting for 6.28% of the COG-annotated genes; the quantity of energy production and conversion-related genes (group C) was 107, taking up 5.17% of the COG-annotated genes; and only 16 genes fall into group Q (secondary metabolites biosynthesis, transport, and catabolism; see Supplementary Figure S2).
3.2. Genome functional studies of Ruminiclostridium thermocellum M3
Gene Ontology and KEGG databases were used to obtain the elucidation of the “character” of the coding gene in the bacteria from a macroscopic perspective. In the genome of R. thermocellum M3, 2021 genes were annotated in the GO database (Table 2). In the biological process, the metabolic and cellular processes account for the highest proportion; in molecular function, it is mainly related to catalytic activity and binding; in the cellular component, it is closely compared to cells and cell parts, which indicated that strain M3 has more proteins involved in metabolism, cell composition, and enzymatic catalysis. These results are in accordance with the biochemical characteristics of M3, which possesses an outstanding catalytic capacity for cellulose substrate. In the KEGG database, 1,246 genes were annotated, which can be divided into five branches according to the metabolic pathways involved in the genes: cellular processes, environmental information processing, genetic information processing, metabolism, and organismal systems. As shown in Supplementary Figure S3, aligning these annotated genes into metabolic pathways yielded a total of 147 metabolic pathways. It can be clearly seen that the predominant pathway was carbohydrate metabolism and overview with 329 and 229 unigenes, followed by amino acid metabolism and energy metabolism with 197 and 153 unigenes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-t002.jpg)
Table 2. Gene function analysis of Ruminiclostridium thermocellum M3 based on Gene Ontology (GO) annotation.
Meanwhile, the predicted numbers of SignaIP-TM and SignaIP-noTM were 24 and 106, respectively, among the total Signa proteins. A total of four CRISPR arrays were found by CRT (Supplementary Table S2). The number of repeat counts in four arrays was 51, 90, 134, and 144. Moreover, several cas or cas-like genes were found in their neighborhoods. This suggests that R. thermocellum M3 has a defense against phage contamination, as CRISPR is very important in prokaryotes and is involved in resisting foreign phages and plasmids and recognizing and silencing invading functional elements. The results of gene-island (GI) prediction obtained eight GI with an average G + C content of about 36.4%, which is slightly lower than the G + C content of the M3 genome. It is worth noting that the G + C content of GI5 (33.7%) is significantly different than that of other GI, which indicates GI5 may be an exogenous sequence by horizontal transfer. The main components of the exogenous sequence are the transposase protein and the hypothetical protein (Supplementary Table S3).
3.3. Comparative genomics of Ruminiclostridium thermocellum
To reveal the genetic and evolutionary relationships between R. thermocellum M3 and other typical R. thermocellum, the analysis in view of the core-pan gene of the whole genome sequence was conducted. The Gene family boxplot (Supplementary Figure S4) revealed that the pan-genome trend of R. thermocellum is open with the increasing number of R. thermocellum strains sequenced. The open pan-genome indicated that there is a significant capacity for discovering novel genes with the evolution and development of strains. Among the five R. thermocellum strains, the number of genes in the pan-genome is 3,042 and includes a core gene set (2,544 genes; Figure 2). Dispensable genes and unique genes existed across all genomes of five R. thermocellum strains. Among them, strain ATCC 27405 had the largest number of specific genes (267 genes), followed by strain M3 (192 genes). The high resemblance of five R. thermocellum strains was reflected by the large proportion of core genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-g002.jpg)
Figure 2. Venn diagram of core and specific genes among five R. thermocellum strains. Each circle represents an Ruminiclostridium thermocellum strain. The number of orthologous coding sequences (core genome) shared by all strains is shown in the center circle, and the number of specific genes is shown in non-overlapping portions of each oval.
To identify the functional classes of the R. thermocellum pan-genome, the COG database was used to classify the functional genes. In the pan-genome of R. thermocellum, the most core, dispensable, and specific gene clusters fell in the metabolism category (Figure 3). Meanwhile, the results indicated that both gene clusters for general function prediction and unknown functions were also more abundant. Compared with the dispensable and specific gene clusters, the majority of genes in the core gene clusters were involved in translation, ribosomal structure, biogenesis (J), cell wall/membrane/envelope biogenesis (M), carbohydrate transport and metabolism (G), amino acid transport, and metabolism (E). By comparison, the majority of genes in the no-core gene clusters was concerned with housekeeping functions, for example, replication, recombination and repair (L), and defense mechanisms (V). Single nucleotide polymorphism (SNP) represents the variation situation of bacteria in the evolutionary process. We constructed the phylogenetic tree based on SNP, which revealed the similarity in bacterial strain variation in adaptation to the natural environment. Compared with R. thermocellum ATCC 27405, M3 showed stronger evolutionary relationships with the other three R. thermocellum strains (Figure 4). Orthologous protein linear analysis based on MCScanX software was performed to further understand the differences in protein homology and amino acid arrangement between M3 and the other four R. thermocellum strains (Figure 5). Between M3 and other three R. thermocellum strains (DSM 1313, DSM 2360 and AD2), it was found that the proteins are not only homologous but also have a good linear relationship in sequence, while R. thermocellum ATCC 27405 was related to M3 but showed a large number of inversions. At the same time, we performed a genome-linear analysis of M3 and R. thermocellum ATCC 27405 (Supplementary Figure S5). The results were similar to the results of the orthologous protein analysis, and it was more clearly evident that nearly half of the genes were inversions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-g003.jpg)
Figure 3. Distribution of core, dispensable, and specific genes on the Clusters of Orthologous Groups (COG) category of Ruminiclostridium thermocellum M3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-g004.jpg)
Figure 4. Phylogenetic relationship of Ruminiclostridium thermocellum M3 and other four R. thermocellum strains. Numbers along branches indicate bootstrap values with 1,000 times.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-g005.jpg)
Figure 5. Plot of protein linear analysis between Ruminiclostridium thermocellum M3 and other four R. thermocellum.
3.4. Protein-encoding genes related to the CAZyme system of Ruminiclostridium thermocellum M3
The results of the CAZyme system analysis indicated that the multi-modular enzyme system of R. thermocellum M3 consisted of 75 dockerin and eight cohesin, which construct scaffolding structural proteins offering a large number of binding sites for cellulase. The quantity of enzyme protein and composition across dissimilar CAZy families in R. thermocellum M3 were analyzed and compared to those in the other four R. thermocellum to evaluate the inclination for lignocellulose saccharification. In the genome of R. thermocellum, most genes fell into glycoside hydrolases (GHs), carbohydrate binding molecules (CBMs), and glycosyltransferases (GTs), whereas a few genes were annotated in carbohydrate esterases (CEs), polysaccharide lyases (PLs), and auxiliary activities (AAs). To be specific, 271 cellulase proteins were detected in R. thermocellum M3 using dbCAN2 (DIAMOND algorithm; Table 3). Most proteins were detected to be GHs (108 candidates), with GH124 (n = 71), GH9 (n = 17), and GH5 (n = 13) being the most abundant families. In addition, there were 75 carbohydrate-binding module (CBM) proteins, 55 glycosyl transferases (GTs) proteins, 22 carbohydrate esterases (CEs) proteins, seven polysaccharide lyases (PLs) proteins, and four auxiliary activities (AAs) proteins. It is worth noting that there were nine CBM genes and 24 GH genes that had not been identified in R. thermocellum before (Supplementary Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-t003.jpg)
Table 3. Comparison of the number of enzyme protein in Ruminiclostridium thermocellum M3 and other R. thermocellum.
Four encoded proteins associated with cellobiose degradation were found among the cellulase-encoded proteins of strain M3 by Conserved Domain Database (CDD) analysis (Supplementary Figure S6), including BglA (Protein: PROKKA_02182), BglB (Protein: PROKKA_00862), and BglX (Protein: PROKKA_01050 and PROKKA_02059). BglA encoded glycoside hydrolase family 1 (GH1) protein (β-glucosidase). BglB encoded glycoside hydrolase family 1 (GH1) protein (6-phospho-β-glucosidase). BglX (PROKKA_01050) encoded glycoside hydrolase family 3 (GH3) protein (β-glucosidase). BglX (PROKKA_02059) encoded glycoside hydrolase family 3 (GH3) protein (β-glucosidase). The result of protein sequence alignment indicated that the ratio of similarity of pairwise comparison of proteins encoded by three genes (PROKKA_02182, PROKKA_01050, and PROKKA_02059), which correlated with β-glucosidase, were lower than 28.39%. To further demonstrate the diversity among the three Bgls, we constructed phylogenetic trees of Bgls from different strains by MEGA software. Not surprisingly, the phylogenetic tree based on β-glucosidase showed that each of the three Bgls was located in one of three different clusters (Figure 6). Furthermore, a gene encoding β-galactosidase (PROKKA_03012) was present in the M3 genome, which was not found in the other wild R. thermocellum before.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1079279/fmicb-13-1079279-HTML/image_m/fmicb-13-1079279-g006.jpg)
Figure 6. Phylogenetic tree of β-glucosidase from Ruminiclostridium thermocellum M3. Numbers along branches indicate bootstrap values with 1,000 times.
In addition, 11 coding proteins about ABC sugar transport protein and two coding proteins (PROKKA_01086 and PROKKA_02112) about cellobiose phosphorylase coming from glycoside hydrolase family 94 (GH94) were detected in the genome of R. thermocellum M3 (Supplementary Table S5). Two group coding proteins (MglA and UpgA) relating polysaccharide transportation (cellobiose, fructose, and arabinose) ware found in 11 ABC sugar transport proteins by CDD. MglA (PROKKA_00080 and PROKKA_01980) is the ATPase component of ABC sugar transporter, which mainly promotes sugar transport; the primary function of UpgA (PROKKA_01330 and PROKKA_02933) is the transmembrane of sugar substrates, which is a permease component of the ABC transporter. Meanwhile, two unique ABC sugar transport proteins (PROKKA_02931 and PROKKA_02932), which had not been reported in other R. thermocellum, were found in strain M3.
4. Discussion
4.1. Genome specificity of cellobiose hydrolysis by Ruminiclostridium thermocellum M3
Ruminiclostridium thermocellum, as a typical genus of thermophagic cellulosic-degrading bacteria, was a high-utility candidate in lignocellulosic biomass refinement by means of a powerful cellulosome system (Akinosho et al., 2014; Sheng et al., 2016). However, the activity of β-glucosidase (BglA) was low in wild-type R. thermocellum strains (Shinoda et al., 2019), and due to the adsorption of cellulosome to cellulose, only a small fraction of BglA was available on cellulosome. In nature, the efficient decomposition of cellulosic biomass strongly depends on the fast adaptation of R. thermocellum to enable the regulation of its cellulosomal enzymes for a specific substrate composition. Compared to other R. thermocellum, such as genus R. thermocellum ATCC 27045 and AD 2, genus R. thermocellum M3 not only has high cellulosic saccharification ability but also possesses the properties of cellubiose-resistance, which contributes to M3 being a potential industrial option for cellulosic biofuels refine.
In general, the wild type R. thermocellum containing β-glucosidase A (BglA) and β-glucosidase B (BglB), gene BglX (β-glucosidase) was first found in R. thermocellum in this study. In general, β-glucosidase A was sensitive to temperature and thus inactivated under the optimum temperature of R. thermocellum (Yoav et al., 2019). Different from β-glucosidase A, BglB encodes a novel thermostable β-glucosidase B, which is more thermally stable than β-glucosidase A, whereas the biosynthesis of β-glucosidase B is repressed by cellobiose. Moreover, the β-glucosidase B of R. thermocellum was inhibited by low glucose concentration (Waeonukul et al., 2012). Different from BglA and BglB, the gene BglX (Protein:PROKKA_01050) may encode β-xylosidase or β-glucosidase, since R. thermocellum cannot utilize xylose as a carbon source as the previous report (Sheng et al., 2016) suggests that the BglX (Protein: PROKKA_01050) gene mainly plays a role in the hydrolysis of cellobiose in strain M3 and that the BglX (PROKKA_01050) gene encodes an enzyme with β-glucosidase activity predominantly. It is reported that the β-glucosidase encoded by BglX was a periplasmic cellulase that hydrolysis cellobiose into glucose (Souto et al., 2021). In addition, a signal peptide that anchors β-glucosidase BglX to the periplasmatic space was found in the conserved domain of BglX, which might help to reduce the cellobiose concentration in a restricted area of the cell surface (Souto et al., 2021). Kim also found that the β-glucosidase of Aspergillus aculeatus can bind to the yeast surface without any modification (Kim et al., 2013), which suggests that the β-glucosidase of strain M3 could be directly connected to the cell membrane to hydrolyze cellobiose without the involvement of the CBM module.
It is believed that the genome of R. thermocellum does not encode any β-galactosidase (GH42; Ravachol et al., 2016). Surprisingly, GH42 β-galactosidase (PROKKA_03012) was found in the genome of R. thermocellum M3. Some studies found that in addition to β-glucosidase, β-1,4 glucosidic can also be cleaved by β-galactosidase (Morita et al., 2008; Yang et al., 2018). Therefore, GH42 might be another key factor in cellobiose hydrolysis by R. thermocellum M3. Similar to GH42, GH116 was not found in R. thermocellum in the previous study; the GH116 family was reported as a thermophilic β-d-glucosidase, which was found in animals, plants, archaea, and bacteria (Sansenya et al., 2015; Rohman et al., 2019). The GH116 protein from thermophilic bacterial was reported for high hydrolytic activity toward β-1,3- and β-1,4-linked gluco-oligosaccharides and 4-nitrophenyl β-D-glucopyranoside (4NPGlc) artificial substrate (Sansenya et al., 2015); therefore, GH116 protein might be another key factor involved in the hydrolysis of cellobiose by R. thermocellum M3.
The combination of ABC transport protein and cellobiose phosphorylase is a common strategy concerning cellobiose transportation in R. thermocellum (Parisutham et al., 2017). Five putative cellodextrin-specific ABC transporters, labeled as CbpA-D and Lbp, had been identified in R. thermocellum DSM 1313 (Nataf et al., 2009); Yan et al. found that only CbpB plays a key role in cellobiose transport using the functional verification of CbpA-D by genetic inactivation (Yan et al., 2022). We identified four ABC sugar transporters that specifically transport polysaccharides by gene annotation and conserved domain database analysis, but the high affinity for binding to cellobiose needs further proof in future works.
4.2. Genome specificity of carbohydrate binding module of Ruminiclostridium thermocellum
It is believed that the cellulase system of R. thermocellum-cellulosomes contains cellulose-binding modules (CBMs), which leads to catalytic activity varying greatly in different regions. As a result, there are higher local cellobiose concentrations at particular sites. However, the β-glucosidase of wild type R. thermocellum was not bound to the CBM which led to the ineffective hydrolysis of local cellobiose (Yoav et al., 2019). Therefore, the construction of an artificial chimeric cohesion-containing scaffold in which binding the β-glucosidase to the cellulosome and mimicking the enzymatic synergism of native cellulosome systems was feasible to enhance the hydrolysis efficiency of cellobiose (Fierobe et al., 2002, 2005; Morais et al., 2010). Different from other R. thermocellum strains, nine CBM family genes were unique to strain M3, including CBM23, CBM36, CBM37, CBM40, CBM47, CBM53, CBM61, CBM70, and CBM75. It is worth noting that the CBM37s were initially discovered in R. albus on the basis of adhesion-defective strains that lacked specific surface proteins. The CBM37s were collectively shown to exhibit a broad specificity pattern, which indicated a mechanism for binding the parent enzymes to cellulosic substrates (Himmel et al., 2010). It is reported that CBM37 was responsible for anchoring substrate with enzymes to the cell surface. For example, CBM37 was reported as the mode that binds cellobiohydrolase and endoglucanase to the bacterial cell surface by the C terminus of glycoside hydrolases (Ezer et al., 2008). In addition, binding the substrate with enzymes (Himmel et al., 2010) might enhance the linkage between the cellulosome and enzymes related to the cellobiose degradation of R. thermocellum M3.
Meanwhile, CBM13 and CBM15 were also found in R. thermocellum M3. CBM13s acquire a larger variety of carbohydrate binding specificities including endo-β-glucanase (EC 3.2.1.6; Ferrer et al., 1996), α-galactosidase (EC 3.2.1.22; Holan et al., 1993), and some other glycoside hydrolases (White et al., 1995). Different from CBM13, CBM35 was reported to have conserved ligand specificity, which is often appended to plant cell wall-degrading enzymes (Montanier et al., 2009) and xylan-degrading enzymes (Chen et al., 1995), but it is often found in β-galactosidase (EC 3.2.1.23; Ali-Ahmad et al., 2017). Therefore, it is suggested that the β-galactosidase of strain M3 is mainly related to CBM13 and CBM35 modules during the degradation of cellulosic feedstocks.
4.3. Genome specificity of auxiliary activities of Ruminiclostridium thermocellum
AA6 and AA8 families were first identified in R. thermocellum in this study. In general, members of the AA8 family (cellobiose dehydrogenase, CDH) can be isolated or appended to a CBM. Proteins contain iron reductase domains and may generate reactive oxygen species that could contribute to the non-enzymatic degradation of cellulose chains by the generation of highly reactive hydroxyl radicals (OH•) via Fenton’s reaction (Levasseur et al., 2013). It is reported that the poor cellobiose availability of the substrate was a limiting factor to CDH activity; on the surface of R. thermocellum, cellobiose was accumulated in a restricted area, which was a sufficient substrate for the CDH. Good contact with CDH might play an important role in reducing the concentration of cellobiose in the environment (Kittl et al., 2012) and potentially releasing the strain M3 from the feedback inhibition of cellobiose.
5. Conclusion
The genome of R. thermocellum M3 harbored a high level of genomic uniqueness compared to other wild R. thermocellum. The majority of genes of R. thermocellum M3 fell into GHs and CBMs. Moreover, some unique genes were found in R. thermocellum M3 which were not found in R. thermocellum, which belong to the Auxiliary Activity Family (AA), cellulose-binding modules Family (CBM), Carbohydrate Esterase Family (CE), Glycosyl Transferases Family (GT), and Glycoside Hydrolase Family (GH). The hydrolysis of cellobiose by R. thermocellum M3 was conducted by the synergy of BglA and BglX, which not only hydrolyze cellobiose but also increase the affinity of Bgl for cellulosic feedstocks, increasing catalytic activity. Meanwhile, the GH42, GH116, CBM37, and AA8 families might participate in the cellobiose degradation that released the R. thermocellum M3 from the feedback inhibition of cellobiose.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
ST, MQ, LZ, SC, LiL, and LiuL contributed jointly to all aspects of the work reported in the manuscript. ST and LZ designed the experiment. MQ performed the experiments. ST, MQ, LiL, SC, and LiuL contributed to the data analysis. ST and MQ drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (No. 51908200), the Heilongjiang Provincial College Youth Innovation Talent Project (No. UNPYSCT-2020028), and the Basic Scientific Research Operating Expenses (No. 2021-KYYWF-1465).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1079279/full#supplementary-material
References
Akinosho, H., Yee, K., Close, D., and Ragauskas, A. (2014). The emergence of clostridium thermocellum as a high utility candidate for consolidated bioprocessing applications. Front. Chem. 2:66. doi: 10.3389/fchem.2014.00066
Ali-Ahmad, A., Garron, M.-L., Zamboni, V., Lenfant, N., Nurizzo, D., Henrissat, B., et al. (2017). Structural insights into a family 39 glycoside hydrolase from the gut symbiont Bacteroides cellulosilyticus WH2. J. Struct. Biol. 197, 227–235. doi: 10.1016/j.jsb.2016.11.004
Berlin, K., Koren, S., Chin, C.-S., Drake, J. P., Landolin, J. M., and Phillippy, A. M. (2015). Assembling large genomes with single-molecule sequencing and locality-sensitive hashing. Nat. Biotechnol. 33, 623–630. doi: 10.1038/nbt.3238
Chen, H., Leipprandt, J. R., Traviss, C. E., Sopher, B. L., Jones, M. Z., Cavanagh, K. T., et al. (1995). Molecular cloning and characterization of bovine β-mannosidase. J. Biol. Chem. 270, 3841–3848. doi: 10.1074/jbc.270.8.3841
Ezer, A., Matalon, E., Jindou, S., Borovok, I., Atamna, N., Yu, Z., et al. (2008). Cell surface enzyme attachment is mediated by family 37 carbohydrate-binding modules, unique to Ruminococcus albus. J. Bacteriol. 190, 8220–8222. doi: 10.1128/JB.00609-08
Ferrer, P., Halkier, T., Hedegaard, L., Savva, D., Diers, I., and Asenjo, J. A. (1996). Nucleotide sequence of a beta-1,3-glucanase isoenzyme IIA gene of Oerskovia xanthineolytica LL G109 (Cellulomonas cellulans) and initial characterization of the recombinant enzyme expressed in Bacillus subtilis. J. Bacteriol. 178, 4751–4757. doi: 10.1128/jb.178.15.4751-4757.1996
Fierobe, H. P., Bayer, E. A., Tardif, C., Czjzek, M., Mechaly, A., Belaich, A., et al. (2002). Degradation of cellulose substrates by cellulosome chimeras - substrate targeting versus proximity of enzyme components. J. Biol. Chem. 277, 49621–49630. doi: 10.1074/jbc.M207672200
Fierobe, H. P., Mingardon, F., Mechaly, A., Belaich, A., Rincon, M. T., Pages, S., et al. (2005). Action of designer cellulosomes on homogeneous versus complex substrates - controlled incorporation of three distinct enzymes into a defined trifunctional scaffoldin. J. Biol. Chem. 280, 16325–16334. doi: 10.1074/jbc.M414449200
Grissa, I., Vergnaud, G., and Pourcel, C. (2007). CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 35, W52–W57. doi: 10.1093/nar/gkm360
Guo, H., Chang, Y., and Lee, D.-J. (2018). Enzymatic saccharification of lignocellulosic biorefinery: research focuses. Bioresour. Technol. 252, 198–215. doi: 10.1016/j.biortech.2017.12.062
Hähnke, V., Hofmann, B., Proschak, E., Steinhilber, D., and Schneider, G. (2009). PhAST: pharmacophore alignment search tool. Chem. Cent. J. 3:P67. doi: 10.1186/1752-153X-3-S1-P67
Haldar, D., and Purkait, M. K. (2020). Lignocellulosic conversion into value-added products: a review. Process Biochem. 89, 110–133. doi: 10.1016/j.procbio.2019.10.001
Hall, B. G. (2013). Building phylogenetic trees from molecular data with MEGA. Mol. Biol. Evol. 30, 1229–1235. doi: 10.1093/molbev/mst012
Himmel, M. E., Xu, Q., Luo, Y., Ding, S.-Y., Lamed, R., and Bayer, E. A. (2010). Microbial enzyme systems for biomass conversion: emerging paradigms. Biofuels 1, 323–341. doi: 10.4155/bfs.09.25
Holan, Z. R., Volesky, B., and Prasetyo, I. (1993). Biosorption of cadmium by biomass of marine algae. Biotechnol. Bioeng. 41, 819–825. doi: 10.1002/bit.260410808
Hsiao, W., Wan, I., Jones, S. J., and Brinkman, F. S. L. (2003). IslandPath: aiding detection of genomic islands in prokaryotes. Bioinformatics 19, 418–420. doi: 10.1093/bioinformatics/btg004
Islam, R., Cicek, N., Sparling, R., and Levin, D. (2006). Effect of substrate loading on hydrogen production during anaerobic fermentation by clostridium thermocellum 27405. Appl. Microbiol. Biotechnol. 72, 576–583. doi: 10.1007/s00253-006-0316-7
Khare, S. K., Pandey, A., and Larroche, C. (2015). Current perspectives in enzymatic saccharification of lignocellulosic biomass. Biochem. Eng. J. 102, 38–44. doi: 10.1016/j.bej.2015.02.033
Kim, S., Baek, S.-H., Lee, K., and Hahn, J.-S. (2013). Cellulosic ethanol production using a yeast consortium displaying a minicellulosome and β-glucosidase. Microb. Cell Factories 12:14. doi: 10.1186/1475-2859-12-14
Kittl, R., Kracher, D., Burgstaller, D., Haltrich, D., and Ludwig, R. (2012). Production of four Neurospora crassa lytic polysaccharide monooxygenases in Pichia pastoris monitored by a fluorimetric assay. Biotechnol. Biofuels 5:79. doi: 10.1186/1754-6834-5-79
Kurtz, S., Phillippy, A., Delcher, A. L., Smoot, M., Shumway, M., Antonescu, C., et al. (2004). Versatile and open software for comparing large genomes. Genome Biol. 5:R12. doi: 10.1186/gb-2004-5-2-r12
Lagesen, K., Hallin, P., Rødland, E. A., Stærfeldt, H.-H., Rognes, T., and Ussery, D. W. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108. doi: 10.1093/nar/gkm160
Lamed, R., Kenig, R., Setter, E., and Bayer, E. A. (1985). Major characteristics of the cellulolytic system of clostridium thermocellum coincide with those of the purified cellulosome. Enzym. Microb. Technol. 7, 37–41. doi: 10.1016/0141-0229(85)90008-0
Laslett, D., and Canback, B. (2004). ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 32, 11–16. doi: 10.1093/nar/gkh152
Levasseur, A., Drula, E., Lombard, V., Coutinho, P. M., and Henrissat, B. (2013). Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol. Biofuels 6:41. doi: 10.1186/1754-6834-6-41
Li, W., and Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. doi: 10.1093/bioinformatics/btl158
Li, W., Jaroszewski, L., and Godzik, A. (2001). Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics 17, 282–283. doi: 10.1093/bioinformatics/17.3.282
Li, W., Jaroszewski, L., and Godzik, A. (2002). Tolerating some redundancy significantly speeds up clustering of large protein databases. Bioinformatics 18, 77–82. doi: 10.1093/bioinformatics/18.1.77
Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2013). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42, D490–D495. doi: 10.1093/nar/gkt1178
Maki, M. L., Armstrong, L., Leung, K. T., and Qin, W. (2013). Increased expression of beta-glucosidase a in clostridium thermocellum 27405 significantly increases cellulase activity. Bioengineered 4, 15–20. doi: 10.4161/bioe.21951
Mazzoli, R., and Olson, D. G. (2020). “Chapter three—clostridium thermocellum: a microbial platform for high-value chemical production from lignocellulose” in Advances in Applied Microbiology. eds. G. M. Gadd and S. Sariaslani (San Diego, Calif: Academic Press), 111–161.
Montanier, C., van Bueren, A. L., Dumon, C., Flint, J. E., Correia, M. A., Prates, J. A., et al. (2009). Evidence that family 35 carbohydrate binding modules display conserved specificity but divergent function. Proc. Natl. Acad. Sci. 106, 3065–3070. doi: 10.1073/pnas.0808972106
Morais, S., Barak, Y., Caspi, J., Hadar, Y., Lamed, R., Shoham, Y., et al. (2010). Cellulase-xylanase synergy in designer Cellulosomes for enhanced degradation of a complex cellulosic substrate. MBio 1:e00285-10. doi: 10.1128/mBio.00285-10
Morita, T., Ozawa, M., Ito, H., Kimio, S., and Kiriyama, S. (2008). Cellobiose is extensively digested in the small intestine by beta-galactosidase in rats. Nutrition 24, 1199–1204. doi: 10.1016/j.nut.2008.06.029
Mostajo Berrospi, N., Lataretu, M., Krautwurst, S., Mock, F., Desirò, D., Lamkiewicz, K., et al. (2019). A comprehensive annotation and differential expression analysis of short and long non-coding RNAs in 16 bat genomes. bioRxiv [Preprint]. doi: 10.1101/738526
Nataf, Y., Yaron, S., Stahl, F., Lamed, R., Bayer Edward, A., Scheper, T.-H., et al. (2009). Cellodextrin and Laminaribiose ABC transporters in clostridium thermocellum. J. Bacteriol. 191, 203–209. doi: 10.1128/JB.01190-08
Parisutham, V., Chandran, S. P., Mukhopadhyay, A., Lee, S. K., and Keasling, J. D. (2017). Intracellular cellobiose metabolism and its applications in lignocellulose-based biorefineries. Bioresour. Technol. 239, 496–506. doi: 10.1016/j.biortech.2017.05.001
Ravachol, J., de Philip, P., Borne, R., Mansuelle, P., Mate, M. J., Perret, S., et al. (2016). Mechanisms involved in xyloglucan catabolism by the cellulosome-producing bacterium Ruminiclostridium cellulolyticum. Sci. Rep. 6:22770. doi: 10.1038/srep22770
Rohman, A., Dijkstra, B. W., and Puspaningsih, N. N. T. (2019). β-Xylosidases: structural diversity, catalytic mechanism, and inhibition by monosaccharides. Int. J. Mol. Sci. 20:5524. doi: 10.3390/ijms20225524
Saha, S., Bridges, S., Magbanua, Z. V., and Peterson, D. G. (2008). Empirical comparison of ab initio repeat finding programs. Nucleic Acids Res. 36, 2284–2294. doi: 10.1093/nar/gkn064
Sansenya, S., Mutoh, R., Charoenwattanasatien, R., Kurisu, G., and Ketudat Cairns, J. R. (2015). Expression and crystallization of a bacterial glycoside hydrolase family 116 [beta]-glucosidase from Thermoanaerobacterium xylanolyticum. Acta Crystallogr. Sect. F 71, 41–44. doi: 10.1107/S2053230X14025461
Sheng, T., Zhao, L., Gao, L. F., Liu, W. Z., Cui, M. H., Guo, Z. C., et al. (2016). Lignocellulosic saccharification by a newly isolated bacterium, Ruminiclostridium thermocellum M3 and cellular cellulase activities for high ratio of glucose to cellobiose. Biotechnol. Biofuels 9:172. doi: 10.1186/s13068-016-0585-z
Shinoda, S., Kurosaki, M., Kokuzawa, T., Hirano, K., Takano, H., Ueda, K., et al. (2019). Comparative biochemical analysis of Cellulosomes isolated from clostridium clariflavum DSM 19732 and clostridium thermocellum ATCC 27405 grown on plant biomass. Appl. Biochem. Biotechnol. 187, 994–1010. doi: 10.1007/s12010-018-2864-6
Souto, B. M., de Araujo, A. C. B., Hamann, P. R. V., Bastos, A. R., Cunha, I. S., Peixoto, J., et al. (2021). Functional screening of a Caatinga goat (Capra hircus) rumen metagenomic library reveals a novel GH3 beta-xylosidase. PLoS One 16:e0245118. doi: 10.1371/journal.pone.0245118
Srivastava, N., Srivastava, M., Mishra, P. K., Gupta, V. K., Molina, G., Rodriguez-Couto, S., et al. (2018). Applications of fungal cellulases in biofuel production: advances and limitations. Renew. Sust. Energ. Rev. 82, 2379–2386. doi: 10.1016/j.rser.2017.08.074
Staples, M. D., Malina, R., and Barrett, S. R. H. (2017). The limits of bioenergy for mitigating global life-cycle greenhouse gas emissions from fossil fuels. Nat. Energy 2:16202. doi: 10.1038/nenergy.2016.202
Tian, L., Papanek, B., Olson, D. G., Rydzak, T., Holwerda, E. K., Zheng, T., et al. (2016). Simultaneous achievement of high ethanol yield and titer in clostridium thermocellum. Biotechnol. Biofuels 9:116. doi: 10.1186/s13068-016-0528-8
Usmani, Z., Sharma, M., Awasthi, A. K., Lukk, T., Tuohy, M. G., Gong, L., et al. (2021). Lignocellulosic biorefineries: the current state of challenges and strategies for efficient commercialization. Renew. Sust. Energ. Rev. 148:111258. doi: 10.1016/j.rser.2021.111258
Usmani, Z., Sharma, M., Karpichev, Y., Pandey, A., Chander Kuhad, R., Bhat, R., et al. (2020). Advancement in valorization technologies to improve utilization of bio-based waste in bioeconomy context. Renew. Sust. Energ. Rev. 131:109965. doi: 10.1016/j.rser.2020.109965
Waeonukul, R., Kosugi, A., Tachaapaikoon, C., Pason, P., Ratanakhanokchai, K., Prawitwong, P., et al. (2012). Efficient saccharification of ammonia soaked rice straw by combination of clostridium thermocellum cellulosome and Thermoanaerobacter brockii β-glucosidase. Bioresour. Technol. 107, 352–357. doi: 10.1016/j.biortech.2011.12.126
White, T., Bennett, E. P., Takio, K., Sorensen, T., Bonding, N., and Clausen, H. (1995). Purification and cDNA cloning of a human UDP-N-acetyl-alpha- D-galactosamine: polypeptide N-acetylgalactosaminyltransferase. J. Biol. Chem. 270, 24156–24165. doi: 10.1074/jbc.270.41.24156
Yadav, M., Paritosh, K., and Vivekanand, V. (2020). Lignocellulose to bio-hydrogen: an overview on recent developments. Int. J. Hydrog. Energy 45, 18195–18210. doi: 10.1016/j.ijhydene.2019.10.027
Yan, F., Dong, S., Liu, Y.-J., Yao, X., Chen, C., Xiao, Y., et al. (2022). Deciphering Cellodextrin and glucose uptake in clostridium thermocellum. MBio 13:e0147622. doi: 10.1128/mbio.01476-22
Yang, Y., Li, X., Zhang, Y., Liu, J., Hu, X., Nie, T., et al. (2020). Characterization of a hypervirulent multidrug-resistant ST23 Klebsiella pneumoniae carrying a Bla CTX-M-24 IncFII plasmid and a pK2044-like plasmid. J. Glob. Antimicrob. Resist. 22, 674–679. doi: 10.1016/j.jgar.2020.05.004
Yang, X., Liu, Z., Jiang, C., Sun, J., Xue, C., and Mao, X. (2018). A novel Agaro-oligosaccharide-lytic β-galactosidase from Agarivorans gilvus WH 0801. Appl. Microbiol. Biotechnol. 102, 5165–5172. doi: 10.1007/s00253-018-8999-0
Yoav, S., Stern, J., Salama-Alber, O., Frolow, F., Anbar, M., Karpol, A., et al. (2019). Directed evolution of clostridium thermocellum β-glucosidase a towards enhanced Thermostability. Int. J. Mol. Sci. 20:4701. doi: 10.3390/ijms20194701
Zhang, J., Liu, S., Li, R., Hong, W., Xiao, Y., Feng, Y., et al. (2017). Efficient whole-cell-catalyzing cellulose saccharification using engineered clostridium thermocellum. Biotechnol. Biofuels 10:124. doi: 10.1186/s13068-017-0796-y
Keywords: thermocellum, genome, cellobiose, β-glucosidase, CAZyme
Citation: Tao S, Qingbin M, Zhiling L, Caiyu S, Lixin L and Lilai L (2023) Comparative genomics reveals cellobiose hydrolysis mechanism of Ruminiclostridium thermocellum M3, a cellulosic saccharification bacterium. Front. Microbiol. 13:1079279. doi: 10.3389/fmicb.2022.1079279
Edited by:
Mamoru Yamada, Yamaguchi University, JapanReviewed by:
Yejun Han, Institute of Process Engineering (CAS), ChinaRosa Estela Quiroz Castañeda, Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), Mexico
Copyright © 2023 Tao, Qingbin, Zhiling, Caiyu, Lixin and Lilai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) 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: Sheng Tao, tsheng@usth.edu.cn; Li Zhiling, lzlhit@163.com