- Open Access
Genome sequencing of gut symbiotic Bacillus velezensis LC1 for bioethanol production from bamboo shoots
Biotechnology for Biofuels volume 13, Article number: 34 (2020)
Bamboo, a lignocellulosic feedstock, is considered as a potentially excellent raw material and evaluated for lignocellulose degradation and bioethanol production, with a focus on using physical and chemical pre-treatment. However, studies reporting the biodegradation of bamboo lignocellulose using microbes such as bacteria and fungi are scarce.
In the present study, Bacillus velezensis LC1 was isolated from Cyrtotrachelus buqueti, in which the symbiotic bacteria exhibited lignocellulose degradation ability and cellulase activities. We performed genome sequencing of B. velezensis LC1, which has a 3929,782-bp ring chromosome and 46.5% GC content. The total gene length was 3,502,596 bp using gene prediction, and the GC contents were 47.29% and 40.04% in the gene and intergene regions, respectively. The genome contains 4018 coding DNA sequences, and all have been assigned predicted functions. Carbohydrate-active enzyme annotation identified 136 genes annotated to CAZy families, including GH, GTs, CEs, PLs, AAs and CBMs. Genes involved in lignocellulose degradation were identified. After a 6-day treatment, the bamboo shoot cellulose degradation efficiency reached 39.32%, and the hydrolysate was subjected to ethanol fermentation with Saccharomyces cerevisiae and Escherichia coli KO11, yielding 7.2 g/L of ethanol at 96 h.
These findings provide an insight for B. velezensis strains in converting lignocellulose into ethanol. B. velezensis LC1, a symbiotic bacteria, can potentially degrade bamboo lignocellulose components and further transformation to ethanol, and expand the bamboo lignocellulosic bioethanol production.
Lignocellulose, a widely distributed, renewable and enormous biomass resource, is one of the most important raw materials for bioethanol production . Bamboo, a lignocellulosic feedstock, is a regenerated biomass material with abundant resources, short growth cycle, high yield and similar chemical composition as wood, and it is considered as a potentially excellent raw material [2, 3]. Many studies have evaluated bamboo lignocellulose degradation and bioethanol production, with a focus on using physical and chemical pre-treatment [4, 5]. However, studies reporting the biodegradation of bamboo lignocellulose using microbes such as bacteria and fungi are scarce.
Lignocellulose hydrolysis, especially cellulose degradation, remains a considerable challenge in lignocellulosic bioethanol production . In nature, numerous examples for lignocellulose degradation are present; of these, phytophagous insects are considered the most notable. In these insects, intestinal symbiotic microbes played important roles in lignocellulose degradation . Therefore, the intestines of phytophagous insects were considered as important locations for isolating lignocellulolytic microbes .
Microbial degradation of lignocellulose is a green biological refining method with advantages over physical and chemical methods . The bacterial genus Bacillus is an excellent degrader that exhibits various abilities for degrading lignocellulose biomass, including cellulose, hemicellulose and lignin [10, 11]. Furthermore, genome sequencing, considered an efficient method for investigation of function, has been utilized in lignocellulose degradation research. However, the lignocellulose degradation of Bacillus is still unclear . Dunlap et al.  reported that B. oryzicola and B. methylotrophicus, were classified into the B. velezensis group. Recently, the complete genome and genes associated with lignocellulose degradation of several B. velezensis strains were sequenced and are enriched in the genome [14,15,16]. However, its potential application in converting lignocellulose into bioethanol has received little attention.
In the present study, we isolated an endophytic bacteria from the gut of Cyrtotrachelus buqueti that showed a bamboo lignocellulose-degrading ability  and sequenced the whole genome of the bacteria B. velezensis LC1, determined the cellulase activities and analysed the ethanol production of bamboo shoot. CAZy genes involved in degradation of lignocellulose were identified through genomic analysis. The chemical changes of the cell wall components were investigated, as well as the hydrolytic and ethanol-fermenting properties of bamboo shoots.
Results and discussion
Identification and cellulose-degrading potential of Bacillus velezensis LC1
Five cellulolytic strains, including PX9, PX10, PX11, PX12 and PX13, which produced clear zones around the colonies after Congo red staining, were isolated from the intestine of C. buqueti on CMC agar. Among the five strains, PX12 exhibited the highest cellulose hydrolysis capacity, with a higher hydrolysis capacity ratio (HCR: 4.71) than PX9 (HCR: 2.41), PX10 (HCR: 1.92), PX11 (HCR: 2.12) or PX13 (HCR: 2.56), as determined using the cellulose hydrolysis assay (Fig. 1a, b; Additional file 1: Figure S1). Based on the HCR ratio, many potent cellulolytic bacteria were previously screened from various regions, such as Geobacillus sp. from a hot spring and Paenibacillus lautus BHU3 from a landfill site [18, 19]. Similarly, PX12 was considered as a good cellulolytic bacterium and was used for further study.
PX12 was identified and confirmed by 16S rRNA sequencing. Phylogenetic analysis of the obtained PX12 sequence revealed a 99% resemblance with B. velezensis strain JTYP2 (NZ_CP020375.1) (Fig. 1c). Moreover, a house-keeping gene was used for phylogenetic analysis; it showed that the sequence of PX12 had a 98.6% similarity with B. velezensis strain S3-1 (NZ_CP016371.1) (Fig. 1d). Overall, the PX12 was identified as B. velezensis and was named as B. velezensis LC1.
Several B. velezensis strains have been noted for their lignocellulose-degrading abilities [14,15,16]. To investigate the cellulose-degrading ability, B. velezensis LC1 was cultured on CMC agar to determine cellulase activities by the dinitrosalicylic acid spectrophotometric (DNS) method for 6 days (Fig. 1e–g) . The cellulase activities of strain LC1 were then determined. The endoglucanase activity was 0.689 ± 0.011 U/ml at day 1 and increased to 0.752 ± 0.013 U/ml at day 6, which was in accordance with the exoglucanase activity (from 0.359 ± 0.016 U/ml to 0.385 ± 0.022 U/ml), whereas the β-glucosidase activity decreased from day 6 to day 1. Previous studies have reported the cellulase activities of other lignocellulolytic Bacillus strains. For example, Bacillus sp. 275, Bacillus sp. R2, B. velezensis 157 and B. velezensis ZY-1-1 [10, 11, 21, 22] showed similar results as those achieved in our study. This indicated that the strain played a potential role in cellulose degradation.
Genome sequencing and assembly of Bacillus velezensis LC1
The genome contributes to a clear understanding of bacterial decomposition mechanisms of cellulose; thus, the genome of B. velezensis LC1 was analysed to decipher the genetic code involved in cellulose degradation. The complete genome sequence of B. velezensis LC1 was assembled into a ring chromosome with 3,929,782 bp and had a GC content of 46.5% (Fig. 2). A length of 3,502,596-bp genes was found based on gene prediction, and the ratio of gene length/genome was 89.13%. The intergene region/whole genome ratio was 7.19%, and the GC contents of the gene and intergene region were 47.29% and 40.04%, respectively. Furthermore, 4018 CDSs were contained in the genome, and all were assigned functions. CDSs were further annotated in NR, Swiss-Prot, COGs, KEGGs, GO and Pfam, and their numbers were 4018, 3520, 2996, 2186, 2718 and 3315, respectively (Table 1).
COGs involved in carbohydrate metabolism
In total, 3046 genes were classified into 2996 COGs, of which carbohydrate transport and metabolism, amino acid transport and metabolism and transcription were the most enriched COGs, which represented 9.46%, 7.62% and 7.29%, respectively (Fig. 3a). To elucidate the function of B. velezensis LC1 in cellulose degradation at the genetic level, specific COGs involved in carbohydrate metabolism were analysed. A total of 222 genes were annotated into carbohydrate metabolism, including 130 COGs, of which the most abundant COGs were COG0366 (alpha-amylase), COG0524 (pfkb domain protein), COG2814 (Major facilitator), COG0726 (4-amino-4-deoxy-alpha-l-arabinopyranosyl undecaprenyl phosphate biosynthetic process), COG1263 (PTS System), COG0477 (major facilitator superfamily), COG1455 (pts system), COG1940 (ROK family) and COG2723 (beta-glucosidase) (Additional file 2: Table S1). COG366 encodes an alpha-amylase that acts on a bond between starch and glycogen, hydrolysing polysaccharides into glucose and maltose . COG2814 is involved in cellular transport of some complexes, such as carbohydrates and amino acids. COG0477, a secondary active transporter, helps to catalyse the transport of various substrates [24, 25]. Moreover, other important COGs in carbohydrate metabolism were annotated, e.g. COG0395 was reported to participate in carbohydrate uptake  and COG1109 catalysed the conversion of glucosamine-6-phosphate . The high diversity of function annotations indicated that B. velezensis LC1 had a potent capability in lignocellulose degradation.
GO terms annotations
To explain the relevance of the genome of B. velezensis LC1, GO analysis was used to categorize genes into three categories according to matches with known sequences. In three categories, molecular function contained most numerous GO terms and gene number (3730), followed by biological process (Gene number: 3025) and cellular component (Gene number: 1637) (Fig. 3b). In molecular function, the most five pathway was ATP binding (GO:0005524; 314 genes), DNA binding (GO:0003677; 254 genes), transcription factor activity (GO:0003700; 115 genes), metal ion binding (GO:0046872; 114 genes) and hydrolase activity (GO:0016787; 84 genes). Oxidation–reduction process (GO:0055114) and regulation of transcription (GO:0006355) were most pathways in the biological process, and integral component of membrane (GO:0016021), cytoplasm (GO:0005737) and plasma membrane (GO:0005886) pathways in cellular component.
Furthermore, we analysed the GOs that are associated with carbohydrate metabolism. We identified 114 GO items associated with carbohydrate metabolism, including GO:0004553 (hydrolase activity that hydrolyses O-glycosyl compounds), GO:0005975 (carbohydrate metabolic processes) and GO:0016787 (hydrolase activity) (Additional file 3: Table S2).
The CDSs of B. velezensis LC1 were submitted to KAAS and KEGG pathways to identify metabolism pathways (Additional file 4: Table S3). As shown in Fig. 3c, of the six classification of KEGG pathways, metabolism contained the most numbers of genes, followed by environmental information processing. In KEGG metabolism annotations of B. velezensis LC1, carbohydrate metabolism and amino acid metabolism, which are considered its main functions, contained 392 and 285 genes, respectively. For these metabolisms, some pathways were dominant, such as sucrose and starch metabolism (ko00500), glycolysis/gluconeogenesis (ko00010), and amino and nucleotide sugar metabolism (ko00520). Forty-one genes were related to ko00500, and common enzyme endoglucanase (EC.126.96.36.199), present in ko00500, was involved in cellulose degradation (Fig. 3c). Starch and sucrose metabolic pathways occurred in B. velezensis LC1, indicating that cellulose could be hydrolysed into cellobiose and ultimately, β-d-glucose. In the genome, 39 genes were found in ko00010, in which d-glucose was phosphorylated into d-glucose-6-phosphate. In addition, ko00010 was linked with other pathways. For example, d-glucose-6-phosphate could be converted to pyruvate, which can be oxidized to acetyl-CoA, having an ability to enter the citrate cycle. Furthermore, ko00520 indicated that glucose from the ko00010 finally entered other pathways under various catalytic reactions. α-d-galactose can be transferred and isomerized in ko00520 and then entered into the ascorbate and aldarate metabolism pathways. Additionally, fructose, 1,4-β-d-xylan, and extracellular mannose were metabolized in ko00520.
The annotation involved in the degradation of lignin or aromatic compounds has also been identified. We identified various enzymes associated with lignin degradation, including oxidoreductase, reductases, dehydrogenases, esterases, thioesterases, transferases and hydrolases. Moreover, 13 monooxygenases, 12 dioxygenases, 2 peroxidases (including one DyP-type peroxidase) and 1 laccase were been identified (Additional file 5: Table S4).
Carbohydrate-active enzyme (CAZyme) annotation
The bCAN carbohydrate-active enzymes (CAZy) annotation algorithm was used to analyse CAZy annotations to identify genes involved in lignocellulose degradation. The results showed that 136 genes were identified from CAZy families and distributed into five subfamilies. In B. velezensis LC1, glycoside hydrolases (GHs), which play key roles in carbohydrates degradation, contained 44 members . Additionally, 38 glycosyl transferases (GTs), 30 carbohydrate esterases (CEs), 3 polysaccharide lyases (PLs), 6 enzymes with auxiliary activities (AAs) and 15 carbohydrate-binding modules (CBMs) were identified (Fig. 3d; Additional file 6: Table S5).
The GH family contained various hydrolases that acted on the glycosidic bond. For example, endoglucanases (EC 188.8.131.52) from Bacillus that have cellulose degradation function usually belonged to the GH5 families  (Table 2). In the genome, six GH13s were obliged to hydrolyze starch, such as α-amylase, α-glucosidase and α-glycosidase . Four GH4s, three GH1s and one GH16 exhibited potential cellulose degradation because of enzyme function. Three GH32s were found in the genomic annotations, which contained some hydrolases and levanases, thereby revealing their ability to hydrolyse sucrose . Additionally, Four GH43s, two GH51s and one GH30, responsible for xylan degradation, were considered as other important members for hemicellulose degradation (Table 2). GH43 is an important component of xylan degradation system . One GH5, GH30 and GH1 each were annotated as potential β-glucosidases to utilize cellobiose. One GH53, which hydrolyses (1 → 4)-β-d-galactosidic linkages, was reported as an endo-1,4-β-galactosidase . β-mannosidase was classified into GH26, hydrolysing parts of mannan polysaccharides . Additionally, maltose phosphorylase, which belongs to GH65, was reportedly involved in trehalose degradation .
CEs contributing to the decomposition of xylans were also identified in the genome, including two CE3s, one CE7s, three CE10s, and seven CE4s. CE3 as a potential acetyl xylan esterase enhanced xylan solubilization . The acetylxylan esterase CE7 was considered as a capable xylan-degrading enzyme . CE10 previously exhibited carboxylesterase and xylanase activities involved in hemicellulose degradation . Polysaccharide deacetylases, which play a role in degrading polysaccharides and are classified as a CE4, were also identified. CE4 contained not only highly specific acetylxylan esterases, but also peptidoglycan N-deacetylates involved in chitin degradation .
Moreover, two PL1s and one PL9 were annotated to degrade pectin. Pectate lyase (EC 184.108.40.206) usually has (1 → 4)-α-d-galacturonan cleavage function, causing oligosaccharides present at the end  (Table 2). AA4, AA6, AA7 and AA10 were included in the genome. AA4 included vanillyl-alcohol oxidases, which could transform some phenols . Additionally, AA7 enzymes were involved in biotransformation or detoxification of lignocelluloses . Finally, we proposed a hypothetical cellulose-degrading and ethanol-producing pathway for B. velezensis LC1 (Table 3).
Comparative genomic analysis of CAZymes with other B. velezensis strains
The assembled genome of B. velezensis LC1 was compared to the genomes of other 10 B. velezensis strains, including, B. velezensis S3-1, B. velezensis LS69, B. velezensis JTYP2, B. velezensis DR-08, B. velezensis FZB42, B. velezensis LPL-K103, B. velezensis TB1501, B. velezensis UCMB5036, B. velezensis LB002 and B. velezensis 157. The result showed that GH, GT and PL family numbers were the same in theses strains, while the strain LC1 contained more CE and AA family members and less CBM members (Table 4). The coexistence of these genes suggests that they play important roles in the enzymatic degradation of cellulose and hemicellulose. We consider these degradation enzymes in B. velezensis to have potential use for bioethanol production.
Gene expression analysis for B. velezensis LC1 cultured in bamboo shoot powder or glucose medium
Expression profiles of selected B. velezensis LC1 lignocellulolytic enzymes genes were monitored in cultures grown on bamboo shoot powder compared to cultures utilizing glucose. Quantitative real-time PCR (qRT-PCR) investigation results of genes encoding endoglucanase (gene1950), beta-glucanase (gene3931), 6-phospho-beta-galactosidase (gene1280), glucan endo-1,6-β-glucosidase (gene2084), alpha-glucosidase (gene2945), acetyl xylan esterase (gene0351), xylan 1,4-beta-xylosidase (gene1870), arabinoxylan arabinofuranohydrolase (gene1955), alpha-N-arabinofuranosidase (gene2785), and arabinan endo-1,5-alpha-l-arabinosidase (gene2795) are shown in Fig. 4. Transcript levels of all the genes of interest were significantly up-regulated (at least P < 0.05) in BSP cultures compared to glucose cultures. It indicated that these genes involved in BSP degradation by B. velezensis LC1.
Cellulose degradation efficiency and fermentation efficiency of bamboo shoots by B. velezensis LC1
Several previously reported genomes of B. velezensis strains contained genes encoding enzymes having lignocellulose-degrading potential [14,15,16]. However, lignocellulose-degrading abilities have not been completely verified in B. velezensis . B. velezensis LC1, isolated from the intestine of C. buqueti, which was reported to be a microflora with lignocellulose-degrading ability, was prepared in a bamboo shoot powder (BSP) degradation assay . We first determined the cellulose degradation efficiency to be 39.32% (Fig. 5a). Furthermore, we determined the glucose and xylose content in the degradation products to be 55.30 ± 1.40 mg/L and 488.81 ± 45.06 mg/L, respectively (Fig. 5b). The reducing sugars in the culture medium was mainly derived from the hydrolysis of cellulose and hemicellulose in BSPs, the reducing sugar content was determined to reflect the degree of conversion of lignocellulose. It indicated that the cellulose and hemicellulose of BSP were degraded with incubation with B. velezensis LC1.
Shimokawa et al.  reported that bamboo shoot was an excellent biomass stock for ethanol production due to its high saccharification efficiency. To determine the ethanol production capacity of B. velezensis LC1, a 6-day-treated hydrolysate was prepared for ethanol fermentation using glucose-fermenting S. cerevisiae and xylose-fermenting E. coli KO11. Ethanol production from BSP hydrolysate using B. velezensis LC1 is illustrated in Fig. 5c. The result showed that the ethanol yield continually increased during 48–96 h, and reached 7.21 ± 0.24 g/L at 96 h, while the reducing sugar continually decreased after incubation 48 h. This indicated that B. velezensis LC1 had a potent ability to convert lignocelluloses in bamboo shoot to ethanol.
The genome of a gut symbiotic bacteria, B. velezensis LC1, was sequenced and analysed. The genome comprises a ring chromosome of 3,929,782 bp and has a GC content of 46.5%. A total of 136 CAZyme genes involved in lignocellulose degradation were annotated in the genome, and cellulose-degrading and ethanol-producing pathways were proposed. Moreover, the expression level of some CAZyme genes involved in cellulose and hemicellulose degradation, were up-regulated in bamboo shoot powder substrate. Moreover, a transcriptomic study would perform to identify the role of CAZyme genes in lignocellulose degradation as a future study. The cellulose degradation of bamboo shoot by the strain was 39.32%, and the hydrolysate was subjected to ethanol fermentation; the ethanol yield was 7.2 g/L at 96 h. This study suggests that B. velezensis LC1 could be used for bamboo lignocellulose degradation and bioconversion of lignocelluloses to ethanol.
Insect sample, isolation and identification of cellulolytic bacteria
Cyrtotrachelus buqueti specimens were sampled from the Muchuan County (E 103° 98′, N 28° 96′), China. Gut was extracted from individual insects and stored at 4 °C for isolation of bacteria. The gut was blended, homogenized and serially diluted (10−1 to 10−9), and inoculums of 10−7 to 10−9 dilution were plated on carboxymethyl cellulose (CMC) agar  for cellulolytic bacteria screening. Congo red dye was used to screen the cellulose-degrading bacteria as described by Teather and Wood . Hydrolysis zone = clearance zone/colony diameter.
Molecular characterization of bacterial isolate
The 16sRNA V3–V4 region was considered for amplification for bacterial identification and amplified using bacterial primers (27F 5′-AGAGTTTGATCMTGGCTCAG-3′ and 1492R 5′-TACGGYTACCTTGTACGACTT-3′); moreover, F 5′-GCCCATATTTCCATTTCTCC-3′ and R 5′-GTGGTCGTTATGGAAATAAAGG-3′ were selected for amplification of the house-keeping gene rpoB. Thermocycling conditions were as follows: initial denaturation at 94 °C (2 min), followed by 30 cycles of denaturation at 94 °C (30 s), annealing at 55 °C (30 s) and extension at 72 °C (100 s), ultimately extending at 72 °C (2 min). The amplicons were checked by electrophoresis on a 1% agarose gel. MEGA5 was used to establish phylogenetic relationships among the obtained sequence and reference genes that were retrieved in NCBI GenBank through the neighbour-joining method.
Genomic DNA extraction and genome sequencing
Genomic DNA of B. velezensis LC1 was extracted according to the CTAB method based on the protocol described by Lin , and it was purified using the Wizard Genomic DNA Purification Kit (Vazyme Biotech Co., Ltd, Nanjing, China). Whole-genome sequencing, performed by Shanghai Majorbio Bio-pharm Technology Co., Ltd (Shanghai, China), was fulfilled through PacBio systems with an average genome coverage of 100 × for the raw data . HGAP 2.0 was used to filter and assemble reads to the scaffold.
Function annotation of B. velezensis LC1
Glimmer 3.02 was used to predict coding DNA sequences (CDSs). A BLAST search was then conducted for CDSs in some widely used databases: NCBI non-redundant (NR) database, Swiss-Prot, COGs, Gene Ontology (GO) and KEGGs [48, 49]. GO, an important bioinformatics tool, unified expressions of gene and genetic products in all species . KEGGs, a database resource to understand high-level functions, also could analyse metabolic pathways. Additionally, the CAZymes were identified, classified and annotated using CAZymes database (CAZyDB: http://www.cazy.org/).
Quantitative real‑time PCR
The primers used for qRT-PCR in this study were performed in Additional file 7: Table S6. PCR was performed under following conditions: 10 min initial denaturation at 95 °C, 45 cycles of 5 s denaturation at 95 °C, 50–65 °C anneal for 30 s, and 30 s extension at 55 °C, finally 10 s extension at 95 °C. All experiments were performed three times and analysed by 2− ΔΔCT Method. 16S rRNA was used as reference gene.
Bamboo shoots degradation by B. velezensis LC1
To obtain fermentable sugar from bamboo shoots, B. velezensis LC1 was used to degrade the bamboo shoot powder (BSP), which was prepared as described by Luo et al. . B. velezensis LC1 was cultured in liquid medium at a pH of 7.2, temperature of 37 °C and 200 rpm for 6 days. The culture medium comprised BSP 10 g/L, (NH4)2SO4 2 g/L, K2HPO4 1 g/L, KH2PO4 1 g/L, MgSO4 0.2 g/L, CaCl2 0.1 g/L, FeSO4·7H2O 0.05 g/L and MnSO4·H2O 0.02 g/L. The reaction mixture was incubated at 100 °C for 30 min to terminate the reaction and centrifuged at 13,000 rpm for 10 min, and the hydrolysate-containing supernatant and deposit were collected separately. The obtained deposit was dried and weighed to determine the cellulose levels using the Van Soest method . The hydrolysate-containing supernatant was used to determine the glucose and xylose contents according to the NREL methods . The BSP hydrolysate supernatant was sterilized and stored at − 20 °C for ethanol fermentation.
Saccharomyces cerevisiae, a glucose-fermenting strain, was pre-cultured in YPD at 30 °C for 24 h, and E. coli KO11, a xylose-fermenting strain, was cultured in LB at 37 °C for 24 h. S. cerevisiae (50 g/L) and E. coli KO11 (100 g/L) were prepared after centrifugation of pre-cultured cells. The initial cell concentrations were 0.33 g/L (S. cerevisiae) and 1.0 g/L (E. coli KO11) at the beginning of fermentation. 100 ml hydrolysate was used for ethanol fermentation in 250-ml serum bottles under anaerobic conditions. The fermentation was performed at 37 °C and 200 rpm for 96 h. After 48 h fermentation, ethanol production was monitored every 12 h. The ethanol concentration was determined via HPLC. All reactions were repeated three times.
Availability of data and materials
The sequence reads from this article have been deposited at the NCBI Sequence Read Archive under the accession PRJNA574012. The assembly data set supporting the results of this article has been deposited at GenBank under the accession CP044349. The version described in this paper is CP044349.
- B. velezensis :
The National Center for Biotechnology Information
Bamboo shoot powder
Yeast extract peptone dextrose medium
Sequence coding for amino acids in protein
Naik SN, Goud VV, Rout PK, Dalai AK. Production of first and second generation biofuels: a comprehensive review. Renew Sustain Energy Rev. 2010;14:578–97.
Littlewood J, Wang L, Turnbull C, Murphy RJ. Technoeconomic potential of bioethanol from bamboo in China. Biotechnol Biofuels. 2013;6:1.
Ragauskas AJ, Beckham GT, Biddy MJ, Chandra R, Chen F, Davis MF, et al. Lignin valorization: improving lignin processing in the biorefinery. Science. 2014;344:1246843.
Chen TY, Wen JL, Wang B, Wang HM, Liu CF, Sun RC. Assessment of integrated process based on autohydrolysis and robust delignification process for enzymatic saccharification of bamboo. Bioresour Technol. 2017;244:717–25.
Wi SG, Lee DS, Nguyen QA, Bae HJ. Evaluation of biomass quality in short-rotation bamboo (Phyllostachys pubescens) for bioenergy products. Biotechnol Biofuels. 2017;10:127.
Mansour AA, Da CA, Arnaud T, Lu-Chau TA, Fdz-Polanco M, Moreira MT, et al. Review of lignocellulolytic enzyme activity analyses and scale-down to microplate-based assays. Talanta. 2016;150:629–37.
Zhang HY, Jackson TA. Autochthonous bacterial fora indicated by PCR-DGGE of 16S rRNA gene fragments from the alimentary tract of Costelytra zealandica (Coleoptera: Scarabaeidae). J Appl Microbiol. 2008;105:1277–85.
de Gonzalo G, Colpa DI, Habib MH, Fraaije MW. Bacterial enzymes involved in lignin degradation. J Biotechnol. 2016;236:110–9.
Capolupo L, Faraco V. Green methods of lignocellulose pretreatment for biorefinery development. Appl Microbiol Biotechnol. 2016;100:9451–67.
Gong G, Lee SM, Woo HM, Park TH, Um Y. Influences of media compositions on characteristics of isolated bacteria exhibiting lignocellulolytic activities from various environmental sites. Appl Biochem Biotechnol. 2017;183:931–42.
Khelil O, Choubane S, Cheba BA. Polyphenols content of spent coffee grounds subjected to physico-chemical pretreatments influences lignocellulolytic enzymes production by Bacillus sp. R2. Bioresour Technol. 2016;211:769–73.
Gong G, Kim S, Lee SM, Woo HM, Park TH, Um Y. Complete genome sequence of Bacillus sp. 275, producing extracellular cellulolytic, xylanolytic and ligninolytic enzymes. J Biotechnol. 2017;254:59–62.
Dunlap CA, Kim SJ, Kwon SW, Rooney AP. Bacillus velezensis is not a later heterotypic synonym of Bacillus amyloliquefaciens; Bacillus methylotrophicus, Bacillus amyloliquefaciens sub sp. plantarum and ‘Bacillus oryzicola’ are later heterotypic synonyms of Bacillus velezensis based on phylogenomics. Int. J Syst Evol Microbiol. 2016;66:1212–7.
Kim SY, Song H, Sang MK, Weon HY, Song J. The complete genome sequence of Bacillus velezensis strain GH1-13 reveals agriculturally beneficial properties and a unique plasmid. J Biotechnol. 2017;259:221–7.
Liu G, Kong Y, Fan Y, Geng C, Peng D, Sun M. Whole-genome sequencing of Bacillus velezensis LS69, a strain with a broad inhibitory spectrum against pathogenic bacteria. J Biotechnol. 2017;249:20–4.
Niazi A, Manzoor S, Asari S, Bejai S, Meijer J, Bongcam-Rudloff E. Genome analysis of Bacillus amyloliquefaciens subsp. plantarum UCMB5113: a rhizobacterium that improves plant growth and stress management. PLoS ONE. 2014;9:e104651.
Luo CB, Li YQ, Chen Y, Fu C, Long WC, Xiao XM, et al. Bamboo lignocellulose degradation by gut symbiotic microbiota of the bamboo snout beetle Cyrtotrachelus buqueti. Biotechnol Biofuels. 2019;12:70.
Potprommanee L, Wang XQ, Han YJ, Nyobe D, Peng YP, Huang Q, et al. Characterization of a thermophilic cellulase from Geobacillus sp. HTA426, an efficient cellulase-producer on alkali pre-treated of lignocellulosic biomass. PLoS ONE. 2017;12:e0175004.
Yadav S, Dubey SK. Cellulose degradation potential of, Paenibacillus lautus, strain BHU3 and its whole genome sequence. Bioresour Technol. 2018;262:124–31.
Dutta N, Mukhopadhyay A, Dasgupta AK, Chakrabarti K. Improved production of reducing sugars from rice husk and rice straw using bacterial cellulase and xylanase activated with hydroxyapatite nanoparticles. Bioresour Technol. 2014;153:269–77.
Chen L, Gu W, Xu H, Yang G, Shan X, Chen G, et al. Complete genome sequence of Bacillus velezensis 157 isolated from Eucommia ulmoides with pathogenic bacteria inhibiting and lignocellulolytic enzymes production by SSF. 3 Biotech. 2018;8:114.
Zhang ZY, Raza MF, Zheng ZQ, Zhang XH, Dong XX, Zhang HY. Complete genome sequence of Bacillus velezensis ZY-1-1 reveals the genetic basis for its hemicellulosic/cellulosic substrate-inducible xylanase and cellulase activities. 3 Biotech. 2018;8:465.
Janeček Š, Gabriško M. Remarkable evolutionary relatedness among the enzymes and proteins from the α–amylase family. Cell Mol Life Sci. 2016;73:2707–25.
Chaudhary N, Kumari I, Sandhu P, Ahmed M, Akhter Y. Proteome scale census of major facilitator superfamily transporters in Trichoderma reesei using protein sequence and structure based classification enhanced ranking. Gene. 2016;585:166–76.
Madej MG, Sun L, Yan N, Kaback HR. Functional architecture of MFS d-glucose transporters. Proc Natl Acad Sci USA. 2014;111:E719–27.
Wang C, Dong D, Wang H, Müller K, Qin Y, Wang H, et al. Metagenomic analysis of microbial consortia enriched from compost: new insights into the role of Actinobacteria in lignocellulose decomposition. Biotechnol Biofuels. 2016;9:22.
Rashid N, Imanaka H, Fukui T, Atomi H, Imanaka T. Presence of a novel phosphopentomutase and a 2-deoxyribose 5-phosphate aldolase reveals a metabolic link between pentoses and central carbon metabolism in the hyperthermophilic archaeon Thermococcus kodakaraensis. J Bacteriol. 2004;186:4185–91.
Roth C, Weizenmann N, Bexten N, Saenger W, Zimmermann W, Maier T, et al. Amylose recognition and ring-size determination of amylomaltase. Sci Adv. 2017;3:e1601386.
Pandey S, Kushwa J, Tiwari R, Kumar R, Somvanshi VS, Nain L, et al. Cloning and expression of β-1, 4-endoglucanase gene from Bacillus subtilis isolated from soil long term irrigated with effluents of paper and pulp mill. Microbiol Res. 2014;169:693–8.
Graebin NG, Schoffer Jda N, Andrades D, Hertz PF, Ayub MA, Rodrigues RC. Immobilization of glycoside hydrolase families GH1, GH13, and GH70: state of the art and perspectives. Molecules. 2016;21:1074.
Bezzate S, Steinmetz M, Aymerich S. Cloning, sequencing, and disruption of a levanase gene of Bacillus polymyxa CF43. J Bacteriol. 1994;176:2177–83.
Romero AM, Mateo JJ, Maicas S. Characterization of an ethanol-tolerant 1,4-beta-xylosidase produced by Pichia membranifaciens. Lett Appl Microbiol. 2012;55:354–61.
Vanholme B, Jacob J, Cannoot B, Gheysen G, Haegeman A. Arabinogalactan endo-1,4-β-galactosidase: a putative plant cell wall-degrading enzyme of plant-parasitic nematodes. Nematology. 2009;11:739–47.
Onilude AA, Fadahunsi IF, Antia UE, Garuba EO, Inuwa M. Characterization of crude alkaline β-mannosidase produced by Bacillus sp. 3A isolated from degraded palm kernel cake. J Inorg Biochem. 2013;81:23–7.
Inoue Y, Yasutake N, Oshima Y, Yamamoto Y, Tomita T, Miyoshi S, et al. Cloning of the maltose phosphorylase gene from Bacillus sp. strain RK-1 and efficient production of the cloned gene and the trehalose phosphorylase gene from Bacillus stearothermophilus SK-1 in Bacillus subtilis. Biosci Biotechnol Biochem. 2001;66:2594–9.
Zhang J, Siika-Aho M, Tenkanen M, Viikari L. The role of acetyl xylan esterase in the solubilization of xylan and enzymatic hydrolysis of wheat straw and giant reed. Biotechnol Biofuels. 2011;4:60.
Mamo G, Hatti-Kaul R, Bo M. A thermostable alkaline active endo-β-1-4-xylanase from Bacillus halodurans S7: purification and characterization. Enzyme Microbial Technol. 2006;39:1492–8.
Zhao Z, Liu H, Wang C, Xu JR. Erratum to: comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genom. 2014;15:6.
Biely P. Microbial carbohydrate esterases deacetylating plant polysaccharides. Biotechnol Adv. 2012;30:1575–88.
See-Too WS, Chua KO, Lim YL, Chen JW, Convey P, Mohd Mohidin TB, et al. Complete genome sequence of Planococcus donghaensis JH1(T), a pectin-degrading bacterium. J Biotechnol. 2017;252:11–4.
van den Heuvel RH, Fraaije MW, Mattevi A, van Berkel WJ. Structure, function and redesign of vanillyl-alcohol oxidase. Int Congress Ser. 2002;1233:13–24.
Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol Biofuels. 2013;6:41.
Shimokawa T, Ishida M, Yoshida S, Nojiri M. Effects of growth stage on enzymatic saccharification and simultaneous saccharification and fermentation of bamboo shoots for bioethanol production. Bioresour Technol. 2009;100:6651–4.
Mondéjar LR, Zühlke D, Becher D, Riedel K, Baldrian P. Cellulose and hemicellulose decomposition by forest soil bacteria proceeds by the action of structurally variable enzymatic systems. Sci Rep. 2016;6:25279.
Teather RM, Wood PJ. Use of Congo red-polysaccharide interactions in enumeration and characterization of cellulolytic bacteria from bovine rumen. Appl Environ Microbiol. 1959;43:777–80.
Lin T, Zhu G, Zhang J, Xu X, Yu Q, Zheng Z, et al. Genomic analyses provide insights into the history of tomato breeding. Nat Genet. 2014;46:1220–8.
Westbrook CJ, Karl JA, Wiseman RW, Mate S, Koroleva G, Garcia K, et al. No assembly required: full-length MHC class I allele discovery by PacBio circular consensus sequencing. Hum Immunol. 2015;76:91–896.
Gao XY, Zhi XY, Li HW, Klenk HP, Li WJ. Comparative genomics of the bacterial genus Streptococcus illuminates evolutionary implications of species groups. PLoS ONE. 2014;9:e101229.
Sun C, Fu GY, Zhang CY, Hu J, Xu L, Wang RJ, et al. Isolation and complete genome sequence of Algibacter alginolytica sp. nov., a novel seaweed-degrading Bacteroidetes bacterium with diverse putative polysaccharide utilization Loci. Appl Environ Microbiol. 2016;82:2975–87.
Huntley RP, Sawford T, Mutowo-Meullenet P, Shypitsyna A, Bonilla C, Martin MJ, et al. The GOA database: gene ontology annotation updates for 2015. Nucleic Acids Res. 2015;43:D1057–63.
Van Soest PJ, Robertson JB, Lewis BA. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J Dairy Sci. 1991;74:3583–97.
Sluiter A, Hames B, Ruiz R, Scarlata C, Sluiter J, Templeton D, et al. Determination of structural carbohydrates and lignin in biomass. version 2012. Golden: National Renewable Energy Laboratory; 2012.
The data were analysed on the free online platform of Majorbio Cloud Platform (www.majorbio.com). We also thank other members of the laboratory for suggestions and discussion regarding this work and revision of the manuscript.
This work was supported by Sichuan science and technology program (2019YFG0139).
Ethics approval and consent to participate
Consent for publication
The authors have consented for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Cellulolytic activities of 4 isolates, PX9, PX10, PX11, and PX13, cultured on the CMC agar plate with congo red.
COG annotations of carbohydrate metabolism in genome of B. velezensis LC1.
GO annotations of carbohydrate metabolism in genome of B. velezensis LC1
KEGG Pathway classification of B. velezensis LC1.
Annotation of dioxygenase, monooxygenase, peroxidase, laccase and oxidoreductase genes identified in the genome of B. velezensis LC1.
CAZy annotation in the genome of B. velezensis LC1.
Primers used for RT-qPCR analysis.
About this article
Cite this article
Li, Y., Lei, L., Zheng, L. et al. Genome sequencing of gut symbiotic Bacillus velezensis LC1 for bioethanol production from bamboo shoots. Biotechnol Biofuels 13, 34 (2020). https://doi.org/10.1186/s13068-020-1671-9
- Bacillus velezensis LC1
- Carbohydrate-active enzyme
- Bamboo shoot