- Open Access
Genomic comparison of Clostridium species with the potential of utilizing red algal biomass for biobutanol production
Biotechnology for Biofuels volume 11, Article number: 42 (2018)
Sustainable biofuels, which are widely considered as an attractive alternative to fossil fuels, can be generated by utilizing various biomass from the environment. Marine biomass, such as red algal biomass, is regarded as one potential renewable substrate source for biofuels conversion due to its abundance of fermentable sugars (e.g., galactose). Previous studies focused on the enhancement of biofuels production from different Clostridium species; however, there has been limited investigation into their metabolic pathways, especially on the conversion of biofuels from galactose, via whole genomic comparison and evolutionary analysis.
Two galactose-utilizing Clostridial strains were examined and identified as Clostridium acetobutylicum strain WA and C. beijerinckii strain WB. Via the genomic sequencing of both strains, the comparison of the whole genome together with the relevant protein prediction of 33 other Clostridium species was established to reveal a clear genome profile based upon various genomic features. Among them, five representative strains, including C. beijerinckii NCIMB14988, C. diolis DSM 15410, C. pasteurianum BC1, strain WA and WB, were further discussed to demonstrate the main differences among their respective metabolic pathways, especially in their carbohydrate metabolism. The metabolic pathways involved in the generation of biofuels and other potential products (e.g., riboflavin) were also reconstructed based on the utilization of marine biomass. Finally, a batch fermentation process was performed to verify the fermentative products from strains WA and WB using 60 g/L of galactose, which is the main hydrolysate from algal biomass. It was observed that strain WA and WB could produce up to 16.98 and 12.47 g/L of biobutanol, together with 21,560 and 10,140 mL/L biohydrogen, respectively.
The determination of the production of various biofuels by both strains WA and WB and their genomic comparisons with other typical Clostridium species on the analysis of various metabolic pathways was presented. Through the identification of their metabolic pathways, which are involved in the conversion of galactose into various potential products, such as biobutanol, the obtained results extend the current insight into the potential capability of utilizing marine red algal biomass and provide a systematic investigation into the relationship between this genus and the generation of sustainable bioenergy.
Increasing concerns about greenhouse gas-mediated climate change and the current high energy demands are driving the development of renewable and sustainable sources that can replace non-renewable fossil fuels [1,2,3]. As one of the most promising renewable fuels, biofuels can be generated via the microbial fermentation process using various biomass from the environment [4, 5]. Biobutanol, bioethanol and biohydrogen are crucial biofuels that are regarded as major optional substitutes for fossil fuels [6,7,8], and the development of this sustainable and renewable biomass should facilitate and advance of biofuels production.
Marine biomass, which has gradually attracted attentions, is considered to be one of the potential sources for biofuels conversion due to the high amount of carbohydrates [9,10,11]. Seaweeds represent an abundant and renewable biomass with fast-growing characteristics that are beneficial for the production of third generation biofuels [12, 13], which are often referred as marine macro-algae including red, brown and green algae . As one of predominant sources of marine biomass, red algal biomass is comprised of agar and cellulose that can be hydrolyzed into various simple sugars (e.g., glucose or galactose) for biofuels fermentation [15,16,17] and it has been utilized as a sustainable and environmentally friendly feedstock for biohydrogen production over the past few decades . Furthermore, there is no lignin found within most red algal biomass, which further reduces the cost of the pretreatment to generate utilizable substrates for microorganisms to produce biofuels .
Microorganisms from the Clostridium genus are diverse and include a large group of anaerobic, endospore-forming bacteria, which possesses hundreds of species , and the majority are recognized as the most notable native cellular factories due to their vast range of substrates utilization and metabolic diversity for the generation of various bio-products [21, 22]. Therefore, one of the most feasible options is the application of Clostridial strains to produce biobutanol by utilizing algal biomass through the acetone–butanol–ethanol (ABE) fermentation process . C. pasteurianum was first reported to be capable of converting algal biomass into limited butanol (0.13 g/L) with the presence of 4% of glycerol . In recent studies, Ellis et al.  adopted C. saccharoperbutylacetonicum to utilize wastewater algal biomass to produce butanol with the addition of xylanases and cellulases, whereas green seaweed was also used for butanol production by strains C. acetobutylicum and C. beijerinckii with the co-metabolism of glucose and xylose . Potts et al.  further processed 15.2 g/L of reducing sugar recovered from macro-algae to enhance butanol production up to 4 g/L using both strains of C. beijerinckii and C. saccharoperbutylacetonicum. Moreover, other solventogenic Clostridial strains, such as C. tetanomorphum ATCC 49273, C. aurantibutyricum NCIMB 10659 and C. beijerinckii NCIMB 8052, were also reported to obtain biobutanol from another marine macroalgae feedstock, namely, Ceylon moss .
In addition to biobutanol, other bio-products can also be obtained via Clostridial fermentation from algal biomass. Park et al.  employed anaerobic sewage sludge microflora for biohydrogen production using macro-algae biomass (Laminaria japonica), and wastewater algal biomass was used to ferment bioethanol via the cellulolytic strain C. phytofermentans DSM 1183 . Clostridium species was also reported to co-produce butanol with riboflavin (vitamin B2), a yellow water-soluble vitamin used as an important cofactor in cells, which also provides an economically practicable way to further exploit the process using algal biomass . In addition, the production of butyric acid using red algae Gelidium amansii as the carbon source was also presented . Sivagurunathan et al.  applied the combined inoculation strategy to improve biohydrogen production from galactose, which is the main hydrolysate from algal biomass, and Sund et al.  evaluated the different roles of C. acetobutylicum in the galactose utilization pathway. Therefore, the Clostridial strains, especially C. acetobutylicum and C. beijerinckii, possess the potential to generate value-added bio-products by using galactose, algal hydrolysate or even algal biomass as substrates.
Although there are large amounts of natural microbial isolates with various metabolic pathways involved in the utilization of biomass and the conversion of biofuels , the lack of well-developed genetic tools and the complicated physiological characteristics from various microbial strains resulted in a limited understanding and development of certain microbial groups , and the comprehensive comparison of the respective strains and their specific capabilities are still lacking . With further investigation via genetic and genomic analysis and the recent efforts for the metabolically engineered Clostridial strains, their innate capabilities, especially the possible potential metabolites and the utilization of recalcitrant substrates, can be demonstrated [22, 34]. In addition, the phylogeny function based methodology that was well known for the study of genomic libraries , special functional enzymes  and ecosystem analysis  can also be established to investigate the relationship between their phenotypes and genotypes. Therefore, an approach dealing with the phylogenetic tree based on whole genomic sequences and a functional comparison on genomic scale could also be applied to analyze the metabolic pathway involved in the generation of biofuels or bio-products by Clostridium species.
In this study, two newly isolated galactose-utilizing Clostridium strains were identified as C. acetobutylicum strain WA and C. beijerinckii strain WB via the whole genomic sequencing. In addition to make comparisons of the genome profiles based upon the genomic features of the other 33 Clostridial strains, three representative strains, including C. beijerinckii NCIMB 14988, C. diolis DSM 15410 and C. pasteurianum BC1, were selected to reveal the critical differences among their respective metabolic potential in utilizing algal biomass for various biofuels and/or biochemicals production by comparison with strains WA and WB, which was further experimentally verified via the ABE fermentation process. In total, this work not only presents the metabolic pathway of the bioconversion of galactose to biobutanol by Clostridium sp. strain WA and WB, but also builds up a comprehensive investigation on the metabolic potential of other industrial bio-products using Clostridial strains and algal biomass through whole genomic comparison and evolutional analysis.
Results and discussion
Genomic features of strains WA and WB
Through genomic sequencing and annotation, the metabolic pathways of biofuels/biochemicals production by the galactose-utilizing strains, Clostridium sp. WA and WB, were analyzed and compared to achieve better insights. It was observed from Table 1 that the genome of strain WA is comprised of a circular chromosome and a mega-plasmid, consisting of a genomic size of 4.07 Mbp with a G+C content of 30.8%. However, the genome of strain WB displayed a much larger genome of 5.78 Mbp with a similar G+C content of 29.7% and there was no plasmid detected during the assembly of strain WB. The final annotation of strain WA and WB resulted in 3878 and 5085 coding sequences (CDSs), respectively. A total of 104 RNA sequences, including 31 rRNA genes (5S, 16S and 23S) and 72 tRNA genes, were found in strain WA; however, strain WB appeared to have less rRNA or tRNA genes detected within its genome.
Whole genome-based phylogenetic analysis of strains WA and WB
By a composition-heterogeneous model in the P4 software package, a randomized axelerated maximum likelihood (RAxML) phylogenetic tree  was constructed based on the whole genomes from strains WA and WB (Additional file 1: Tables S1 and S2) together with 33 other available Clostridial strains according to the concatenated alignment of 129 bacterial single copy marker genes with a total of 10793 amino acid sites, and Bacillus licheniformis ATCC 14580 was set as an outgroup. It was observed that these strains were phylogenetically placed into five clades (Fig. 1). The genomes from strains C. acetobutylicum WA, C. arbusti SL206, C. akagii DSM 12554, C. cellulovorans 743B and C. pasteurianum BC1 formed a robustly monophyletic group with the same evolutionary clade (Clade 1), whereas strains C. beijerinckii NCIMB 14988, C. diolis DSM 15410 and C. beijerinckii WB were clustered together into another single subclade (Clade 5). These two newly sequenced strains (strains WA and WB) reveal a relatively far evolutionary relationship/distance, which indicates the possible distinction in their metabolism related to biomass utilization and the conversion of bio-products.
Clostridium strains with the potential of utilizing various biomass (e.g., corn cobs, cassava and rice bran) and the production of biofuels (e.g., butanol) were mainly classified within Clade 1 and Clade 5 such as C. cellulovorans 743B and C. saccharoperbutylacetonicum N1-4 [39,40,41]. Genomes from the same clustered group usually appear to have similar metabolic functions, which demonstrates a comprehensive way to understand certain isolated strain based on the available reference genomes from other strains. Therefore, according to the phylogenetic location and genomic characteristics, the reference genomes from three other Clostridial strains, including C. beijerinckii NCIMB14988, C. diolis DSM 15410 and C. pasteurianum BC1, together with the genomes of strains WA and WB, were selected to determine the differences from their metabolic pathways and provide directions for future fermentation. Via the classification from annotated genes using the COG database (Fig. 2), the distributions of the functional proteins were observed to be the same trend for the affiliation of the selected strains (Fig. 1). The result also suggests that the phylogenetic analysis based on single-copy gene families could be utilized to cluster the individual microbial strains from the same genera into the specific clade using the interaction of the functional genes, which provides more insights rather than using only 16S rRNA gene sequences .
Determination of the genomic characteristics from diverse Clostridial strains
To better understand the basic divergence among reference strains, the genomic features of 35 total Clostridial strains were analyzed and assigned to a relatively clear genome classification (Fig. 3), and it was observed that different species, even from the same genus of Clostridium, possess the significant divergences such as ~ 2.5–7.0 Mb range for the genome size, ~ 2000–6000 range for the genes/proteins and ~ 50–100 range for the tRNA. This possibly led to the distinction on their respective metabolism. The genome of strain WA was observed to be within the median level from those analyzed genomes, whereas strain WB showed a relatively larger genome with more genes/CDSs but less RNAs from its own genome. In addition to the chromosomal DNA, 8 out of 35 strains were determined to have their own separate circular plasmid. The mega-plasmid of WA (pWA), which encoded a total of 178 proteins with the sol operon, is similar to the reported mega-plasmid (pSOL1) from C. acetobutylicum ATCC 824 that contained one vital gene, namely, aldehyde/alcohol dehydrogenase (aad) involved in biobutanol generation . C. saccharoperbutylacetonicum N1-4(HMT), which is another biosolvents-producing strain, carried a similar size of mega-plasmid (Csp_135p, 0.136 Mbp) that is apparently not related to the formation of solvents. However, those smaller plasmids (50–750 kbp) from C. aceticum DSM 1496, C. kluyveri DSM 555, C. pasteurianum BC1, C. tetani E88 and C. tyrobutyricum KCTC 5387 were found without a known role in the Clostridial physiological process .
Comparison of the genome-wide metabolic pathway from representative Clostridial strains
The result obtained from the alignment using the COG database indicated that the distinction of gene counts among the five strains involved the pathways in carbohydrate transport and metabolism (Fig. 2). Both strains WA and BC1 have more genes participating in various metabolic processes; however, their genes referring to the carbohydrate transport and metabolism were found to be much less than that from the other three strains. Clostridium species exhibit a broad substrate range but there have been limited studies of the mechanisms involved in regulation of uptake and metabolism of fermentable carbohydrates . The presence of numerous phosphotransferase systems (PTS) is reported to be significantly related to the uptake of sugars . For example, C. acetobutylicum utilized phosphotransferase system (PTS) transporters for the uptake of disaccharides and hexoses, whereas pentoses were primarily taken up by ATP-binding cassette (ABC) transporters . However, apart from common PTS systems (e.g., PTS-Glc-crr, glucose-specific PTS) involved in both strains WA and WB, more PTS systems were exclusively present in the metabolism of strain WB, such as PTS-Dga-dgaB/dgaC/dgaD (d-glucosaminate specific PTS), PTS-Gam-agaC (galactosamine-specific PTS) and PTS-Ula-ulaA/sgaT (ascorbate-specific PTS), which contributed to the higher number of genes that participated in carbohydrate transport and metabolism in strain WB.
When processing genomic comparison and metabolic reconstruction, analysis of the metabolic capabilities of different Clostridial strains is necessary to consider the relationship between microorganisms and substrates to understand their requirements of carbon catch and energy delivery . The entire genomic distinction among the above five selected Clostridial strains together with two plasmids were compared, and the locations of those obvious syntenic blocks are both highlighted within a circle co-assembly map (Fig. 4). It is worth mentioning that most regions of strain DSM 15410 can match to partial genomic regions of strain WB and NCIMB 14988. However, few similar regions were found from strains WA and BC1 when comparing with the other three strains and even when supplementing their plasmids into the whole genomes. Similarly, few of the same regions could be detected between the galactose-utilizing strain WA and WB (Additional file 2: Figure S1), and most of the genes were involved in maintaining the basic metabolic process such as glycolysis (galactose consumption), TCA cycle and butanoate metabolism (butanol generation). As judged from the genome repertoire (Table 2), five Clostridial strains were able to utilize various saccharides, such as mannose, fructose, glucose and cellulose, to produce various bio-products, and genes from the plasmids in both strains WA and BC1 were also found to participate in the metabolic pathways involved in biofuels generation, especially those from strain WA. However, there are still many differences observed with each other in a metabolism series, such as carbohydrate and energy metabolism.
As determined from Fig. 2, the overall trend of gene distribution in the five strains was similar. Although the number of genes is not directly proportional to the total number of metabolic pathways, there are obvious differences found from the pathway of carbohydrate transport and metabolism. For example, the quantity of genes in strain WB, DSM 15410 or NCIMB 14988 is twice more as much compared to that from strain WA or BC1, which is probably caused by their evolutionary relationships (Fig. 1). Microbial strains from the Clostridium genus have a wide and efficient utilization of a variety of carbon sources such as glucose, galactose, xylan and other polysaccharides. However, as indicated in Table 2, there are many significant differences among these five strains, especially in strains WA and WB. It was observed that strain WA is able to participate in the metabolic pathways involved not only in lactose, trehalose utilization and pyruvate conversion into butanoate, but also in the sulfoquinovosyl diacylglycerol biosynthesis, superoxide radicals degradation and thioredoxin pathway, whereas the number of metabolic pathways in strain WB was observed to be much higher than that in other strains. In addition, strains WA and BC1 had specific pathways involved in the biosynthesis of cellulose, fatty acids and 2,3-dihydroxybenzoate, whereas the metabolic pathways of gluconeogenesis and CDP-diacylglycerol biosynthesis were exclusively observed in strain WB, DSM 15410 and NCIMB 14988. These findings indicate the strong correlation between functional evolutionary relationships of five strains and their metabolic pathways.
In addition, five Clostridial strains can participate in a series of energy metabolism including biofuels production, nitrogen fixation, iron-sulfur cluster biosynthesis, etc. (Table 2). The energy metabolism involved in the sulfoquinovosyl diacylglycerol biosynthesis and the thioredoxin pathway were only discovered in strain WA; however, no such pathways were observed in strain WB, similar to that found in strain DSM 15410 or NCIMB 14988. Moreover, strain WA and BC1 have PTS-Dgl-gamP (d-glucosamine specific PTS), whereas PTS-Mal-malX (maltose specific PTS) is only involved in strain WB, DSM 15410 and NCIMB 14988. There is no doubt that all of the microorganisms, including Clostridium species, should possess vital metabolic pathways in the membrane transport, signal transduction and signaling molecules and interaction to adapt and respond to the culture environments; thus, the comparative genomic analysis of strains WA and WB with the other Clostridial species offers a better understanding of the substrate utilization process and bio-products generation .
Reconstruction of the biofuels-related metabolic pathway via Clostridial genomes
So far, studies on the mechanism and regulation of sugar uptake and transport still remain limited; therefore, the newly isolated strains WA and WB were used to evaluate the feasibility of biofuels (e.g., biobutanol and biohydrogen) generation from galactose, the main component from red algal biomass [49, 50]. Via the reconstruction of the biofuels/biochemicals-related metabolic pathways of strains WA and WB with the references of three other representative Clostridial strains (Fig. 5), the obvious differences between the genomes of strains WA and WB are demonstrated in Table 5. It was observed that the genes of galM (aldose-1-epimerase), galK (galactokinase) and galT (galactose-1-phosphate uridylyltransferase) within the Leloir (LP) pathway for the conversion of galactose into α-d-glucose-1-P were crucial in both strains WA and WB . Moreover, the number of genes participating in the carbohydrate metabolism from strain WB was also found to be higher than that from WA (Table 5), such as phosphoglucomutase (pgm) (Fig. 5).
In addition, C. pasteurianum BC1 was also reported to produce various biofuels (e.g., biobutanol and biohydrogen) [51, 52], and C. diolis DSM 15410 was identified as a 1,3-propanediol (1,3-PDO) producer ; however, C. beijerinckii NCIMB 14988 was used as a common strain for the evolutionary analysis without available bio-products shown . Therefore, the metabolic pathways involved in the production of butanol, hydrogen and 1,3-PDO, together with another potential product (riboflavin, VB2), were reconstructed to better elaborate the possible generation of bio-products by strain WA and WB (Fig. 5). Same as strains NCIMB 14988, DSM 15410 and BC1, both strains WA and WB also possess all the crucial genes to complete the biosynthesis of riboflavin (Fig. 5, Table 2), which indicated the potential of these strains to produce riboflavin directly from galactose when an optimal cultivating condition is provided. Zhao et al.  reported that C. acetobutylicum ATCC 824 could generate riboflavin as a by-product during its ABE fermentation process via the supplementation of sodium acetate, and they also mentioned that the synthetic rate of GTP (precursor of riboflavin) could lead to the generation of riboflavin. Through the over-expression of riboflavin biosynthesis-related operon genes (ribGBAH) from strain ATCC 824, 20 mg/L of riboflavin could be determined extracellularly from E. coli . Therefore, with the identification of the riboflavin metabolic pathway from strains WA and WB, it should be emphasized that the occurrence of the complete RBP (riboflavin biosynthetic pathway) system in both strains WA and WB would improve their possibility of producing riboflavin by regulating various external factors, such as culture pH, precursors or transporters [55, 57]. In total, the findings from the reconstruction of the powerful and complicated metabolic systems in Clostridial strains will facilitate the exploration of possible carbohydrate utilization and value-added products generation via the biochemical and/or molecular methodology .
Determination of biofuels production by galactose-utilizing strains WA and WB
To verify the potential of biofuels production from galactose, a batch fermentation process was conducted using strains WA and WB in the defined culture medium  supplemented with 60 g/L of galactose, respectively. As shown in Fig. 6, strain WA produced 0.95 g/L butanol, 0.11 g/L ethanol and 190 mL hydrogen at the acidogenic stage with the decrement of pH value from 6.5 to 4.0 and relatively fast bacterial growth (OD600 nm = 3.5) during the first 24 h (Fig. 6a). Followed by a second solventogenic stage, strain WA can produce up to 16.98 g/L butanol, 0.88 g/L ethanol and 1077.67 mL hydrogen by consuming almost all of the substrates after 120 h of fermentation (Fig. 6b). Strain WB generated a higher amount of butanol (1.72 g/L) with less acidic intermediates in the first 24 h (Fig. 6c), and finally produced 12.47 g/L of butanol with a very limited amount of ethanol (0.47 g/L) and 506.67 mL hydrogen after 120 h of fermentation (Fig. 6d). The yields of biohydrogen (353.21 and 165.99 mL/g) and biobutanol (0.283 and 0.207 g/g) from WA and WB, respectively, appeared to be much higher than that from other reported butanol-producing bacterial strains (Tables 3 and 4), which indicated a high possibility for adopting these two strains for bioconversion of biofuels from those galactose-rich substrates such as biomass from red alga. However, there were some obvious differences observed between these two strains, such as (i) the rate of galactose utilization by strain WA was much higher rather than that by strain WB (90.93% versus 66.20%), which led to a high production of butanol and hydrogen (Fig. 6a, c); (ii) a neglectable amount of ethanol was detected during the fermentation by strain WB, which was highly related to the lack of the pyruvate decarboxylase (pdc) gene within the genome of strain WB (Fig. 5, Table 5); and (iii) more butanol was formed from strain WA rather than that from strain WB during the fermentation, probably due to the lower amount of butanol synthesis-related gene (bdhB) in the genome of strain WB (Table 5). The butanol dehydrogenases A/B (bdhA/B), which are the critical enzymes for biobutanol synthesis, are both NADH dependent. However, the generation of butyrate during the ABE fermentation could consume ATP and inhibit the synthesis of NADH, which would further affect butanol yield [59, 60], and the pH adjustment during the fermentation could inhibit butyrate kinase (buk) to effectively reduce this influence . With the occurrence of both bdhA and bdhB genes in strain WA, a significantly high amount of butanol was observed from the fermentation of strain WA with the pH adjustment (Fig. 6).
In this study, the genomes of two newly isolated Clostridial strains, WA and WB, with efficient galactose utilization were analyzed to illustrate their capability of synthesizing biofuels and/or biochemicals production, which might be improved via a future gene-engineering process without changing the native functional features. The genomic comparisons with other typical Clostridium species involved in various metabolic pathways also facilitate extension of the current limited understanding of their potential capability of using marine algal biomass as a sustainable substrate. In addition, CRISPR/Cas9-based genome editing using the plasmid of pNICKclos was recently developed to achieve an editing efficiency up to 100% in both C. acetobutylicum ATCC 824 and C. beijerinckii NCIMB 8052, two Clostridial strains with a complete genomic interpretation . Therefore, with a more in-depth elaboration, strains WA and WB will demonstrate their demands of genome editing and regulation control, which also accelerates the progress in further metabolic engineering of the solventogenic Clostridia, regardless of the utilization of sustainable substrates or the generation of potentially industrial products.
Via the comprehensive study of phylogeny and the genomic comparisons for two galactose-utilizing Clostridium strains identified to be C. acetobutylicum strain WA and C. beijerinckii strain WB, we provide a useful approach to highlight significant differences in biofuels-related metabolism. Furthermore, the results also demonstrate that potential products, such as riboflavin, were identified in the Clostridium metabolic pathway with marine biomass. Finally, the present work further extends our current understanding of Clostridia and provides a systematic investigation into the relationship between this genus and the generation of sustainable bioenergy.
Bacterial strains and cultivation conditions
The bacterial strains WA and WB with the capability of utilizing galactose as the sole carbon source were both isolated from mangrove sediments via enrichment using the reinforced clostridial medium (RCM). The cultivation medium of these two strains was prepared by using the defined culture medium amended with 60 g/L of galactose as described by Wu et al. .
Genome sequencing and re-annotation of strains WA and WB
The genomic DNA of strains WA and WB was extracted using the E.Z.N.A.® Bacterial DNA Kit (Omega Bio-Tek, USA) according to manufacturers’ instructions, and applied to whole shotgun sequencing using the Illumina paired-end sequencing technology at the Beijing Genomics Institute (BGI, China). The obtained reads were assembled into contigs in different scales by using SOAPdenovo (V1.05). The re-annotation of whole genomes, including the functional genes and RNA prediction, was performed using the prokaryotic genome annotation system pipeline program (V1.11) , and the identification and classification of the encoded functional proteins was determined based on the Clusters of Orthologous Groups (COG) database.
Whole genome-based phylogenetic analysis of Clostridial strains
Genomes from strains WA and WB, together with those from 33 other bacterial strains from the Clostridium genus and an outgroup strain (Bacillus licheniformis ATCC 14580) affiliated to the same Clostridiaceae family were analyzed to establish their phylogenetic relationship. By referring to the method from Sun et al. , a Clostridial phylogenetic tree was finally constructed with the LG substitution matrix and the gamma model using the RAxML tool (V8.0) based on a concatenation of 10793 amino acid sites over 129 single-copy gene families shared by 36 available genomes  and finally demonstrated via MEGA6 .
Comparative genomic analysis on the genome-wide metabolic pathway
To verify the differential genome-wide metabolic pathways among the Clostridial strains, C. acetobutylicum strain WA and C. beijerinckii strain WB together with their phylogenetically close strains, C. beijerinckii NCIMB 14988, C. pasteurianum BC1 and C. diolis DSM 15410, were selected as the representative strains for genomic comparison. Circos  and Mauve  software was used to compare the assembly differences. The MetaPathways software (V2.0) was further adopted to re-construct the genome-wide metabolic pathways with the following parameters: (i) ORFs detection by Prodigal with minimal length of 60 amino acids and (ii) functional annotation via BLAST with an e value of 10−5 and a Blast-score ratio of 0.4 [69, 70] using the protein databases of KEGG, CAZY, COG, MetaCyc and RefSeq. In addition, the metabolic pathways involved in algal biomass utilization and biochemicals/biofuels synthesis were further validated through the databases of KEGG and TCDB, and reconstructed using Adobe Illustrator CS6 software.
Determination of fermentative products by strain WA and WB
To determine both strains WA and WB and their potential for fermentative products generation using galactose, the activated strains WA and WB were inoculated into 50 mL of galactose-supplemented culture medium and cultivated at 37 °C in a shaker with a rotary speed of 150 rpm for 96 h, respectively. The pH of the culture medium was manually maintained at 5.0–5.5 using a 3 M NaOH solution during the entire fermentation process. Starting from 0 h of fermentation, the production of biohydrogen and biosolvents (i.e., acetone, butanol, ethanol, acetic acid and butyric acid), the bacterial growth as well as the substrate consumption were recorded every 12 or 24 h. Hydrogen was collected by gas sampling bags and determined using a GC-2010Plus gas chromatograph (GC) equipped with a thermal conductivity detector (TCD) and a Supelco 80/100, Porapak-N column (Shimadzu, Japan) as described by the method of Wu et al. . The oven temperature was maintained at 110 °C for 5 min, and argon was used as the carrier gas with a column flow rate of 1 mL/min. Standard gaseous mixtures, which consists of nitrogen (60%), carbon dioxide (15%), carbon monoxide (15%), hydrogen (5%) and methane (5%), were used for quantification. The concentration of the biosolvents was also measured by a GC-2010Plus GC equipped with a flame ionization detector (FID) and a DB-WAXetr column (30 m × 0.25 mm × 0.25 μm ID) based on the method of Xin et al.  with minor modifications. The oven temperature was initially held at 60 °C for 2 min, increased at 15 °C/min to 230 °C, and was then held for 1 min. Helium was used as the carrier gas with a column flow rate of 1 mL/min and a mixed biosolvents standard curve was established for quantification. The bacterial biomass was assessed using a UVmini-1240 UV–visible spectrophotometer (Shimadzu, Japan) at a wavelength of 600 nm, whereas the concentration of galactose was measured using the 3,5-dinitrosalicylic acid (DNS) method .
Munir R, Levin DB. Enzyme systems of anaerobes for biomass conversion. Adv Biochem Eng Biotechnol. 2016;156:113–38.
Gaida SM, Liedtke A, Jentges AH, Engels B, Jennewein S. Metabolic engineering of Clostridium cellulolyticum for the production of n-butanol from crystalline cellulose. Microb Cell Fact. 2016;15:6.
Zhao X, Hwang BG, Gao Y. A fuzzy synthetic evaluation approach for risk assessment: a case of Singapore’s green projects. J Clean Prod. 2016;115:203–13.
Lan EI, Liao JC. Microbial synthesis of n-butanol, isobutanol, and other higher alcohols from diverse resources. Bioresour Technol. 2013;135:339–49.
Shanmugama S, Hari A, Ulaganathan P, Yang F, Krishnaswamy S, Wu YR. Potential of biohydrogen generation using the delignified lignocellulosic biomass by a newly identified thermostable laccase from Trichoderma asperellum strain BPLMBT1. Int J Hydrogen Energy. 2018. https://doi.org/10.1016/j.ijhydene.2018.01.016.
Xue C, Zhao XQ, Liu CG, Chen LJ, Bai FW. Prospective and development of butanol as an advanced biofuel. Biotechnol Adv. 2013;31:1575–84.
Demirbas A. Progress and recent trends in biofuels. Prog Energy Combust Sci. 2007;33:1–18.
Hallenbeck PC, Ghosh D. Advances in fermentative biohydrogen production: the way forward? Trends Biotechnol. 2009;27:287–97.
Hannon M, Gimpel J, Tran M, Rasala B, Mayfield S. Biofuels from algae: challenges and potential. Biofuels. 2010;1:763–84.
Menetrez MY. An overview of algae biofuel production and potential environmental impact. Environ Sci Technol. 2012;46:7073–85.
Wei N, Quarterman J, Jin YS. Marine macroalgae: an untapped resource for producing fuels and chemicals. Trends Biotechnol. 2013;31:70–7.
Raheem A, Wan Azlina WAKG, Taufiq Yap YH, Danquah MK, Harun R. Thermochemical conversion of microalgal biomass for biofuel production. Renew Sustain Energy Rev. 2015;49:990–9.
Hou X, From N, Angelidaki I, Huijgen WJJ, Bjerre AB. Butanol fermentation of the brown seaweed Laminaria digitata by Clostridium beijerinckii DSM-6422. Bioresour Technol. 2017;238:16–21.
Park JI, Lee J, Sim SJ, Lee JH. Production of hydrogen from marine macro-algae biomass using anaerobic sewage sludge microflora. Biotechnol Bioprocess Eng. 2009;14:307–15.
Jol CN, Neiss TG, Penninkhof B, Rudolph B, De Ruiter GA. A novel high-performance anion-exchange chromatographic method for the analysis of carrageenans and agars containing 3,6-anhydrogalactose. Anal Chem. 1999;268:213–22.
Park JH, Yoon JJ, Park HD, Kim YJ, Lim DJ, Kim SH. Feasibility of biohydrogen production from Gelidium amansii. Int J Hydrog Energy. 2011;36:13997–4003.
Dixon PS. Biology of the Rhodophyta. Edinburgh: Oliver & Boyd Edinburgh; 1973.
Sambusiti C, Bellucci M, Zabaniotou A, Beneduce L, Monlau F. Algae as promising feedstocks for fermentative biohydrogen production according to a biorefinery approach: a comprehensive review. Renew Sustain Energy Rev. 2015;44:20–36.
Park JH, Hong JY, Jang HC, Oh SG, Kim SH, Yoon JJ, Kim YJ. Use of Gelidium amansii as a promising resource for bioethanol: a practical approach for continuous dilute-acid hydrolysis and fermentation. Bioresour Technol. 2012;108:83–8.
Bruce T, Leite FG, Miranda M, Thompson CC, Pereira N, Faber M, Thompson FL. Insights from genome of Clostridium butyricum INCQS635 reveal mechanisms to convert complex sugars for biofuel production. Arch Microbiol. 2016;198:115–27.
Poehlein A, Solano JDM, Flitsch SK, Krabben P, Winzer K, Reid SJ, Jones DT, Green E, Minton NP, Daniel R, Durre P. Microbial solvent formation revisited by comparative genome analysis. Biotechnol Biofuels. 2017;10:58.
Tracy BP, Jones SW, Fast AG, Indurthi DC, Papoutsakis ET. Clostridia: the importance of their exceptional substrate and metabolite diversity for biofuel and biorefinery applications. Curr Opin Biotechnol. 2012;23:364–81.
Jang YS, Malaviya A, Cho C, Lee J, Lee SY. Butanol production from renewable biomass by clostridia. Bioresour Technol. 2012;123:653–63.
Nakas J, Schaedle M, Parkinson C, Coonley C, Tanenbaum S. System development for linked-fermentation production of solvents from algal biomass. Appl Environ Microbiol. 1983;46:1017–23.
Ellis JT, Hengge NN, Sims RC, Miller CD. Acetone, butanol, and ethanol production from wastewater algae. Bioresour Technol. 2012;111:491–5.
van der Wal H, Sperber BL, Houweling-Tan B, Bakker RR, Brandenburg W, López-Contreras AM. Production of acetone, butanol, and ethanol from biomass of the green seaweed Ulva lactuca. Bioresour Technol. 2013;128:431–7.
Potts T, Du J, Paul M, May P, Beitle R, Hestekin J. The production of butanol from Jamaica bay macro algae. Environ Prog Sustain. 2012;31:29–36.
Kim Y, Kim D, Kim T, Shin M-K, Kim YJ, Yoon J-J, Chang IS. Use of red algae, Ceylon moss (Gelidium amansii), hydrolyzate for clostridial fermentation. Biomass Bioenergy. 2013;56:38–42.
Fathima AA, Sanitha M, Kumar T, Iyappan S, Ramya M. Direct utilization of waste water algal biomass for ethanol production by cellulolytic Clostridium phytofermentans DSM1183. Bioresour Technol. 2016;202:253–6.
Qureshi N, Blaschek H. ABE production from corn: a recent economic evaluation. J Ind Microbiol Biotechnol. 2001;27:292–7.
Lee KM, Choi O, Kim KY, Woo HM, Kim Y, Han SO, Sang BI, Um Y. Butyric acid production from red algae by a newly isolated Clostridium sp. S1. Biotechnol Lett. 2015;37:1837–44.
Sivagurunathan P, Anburajan P, Kumar G, Arivalagan P, Bakonyi P, Kim SH. Improvement of hydrogen fermentation of galactose by combined inoculation strategy. J Biosci Bioeng. 2017;123:353–7.
Sund CJ, Servinsky MD, Gerlach ES. Differing roles for Clostridium acetobutylicum’s galactose utilization pathways. Adv Microbiol. 2013;03:490–7.
Alper H, Stephanopoulos G. Engineering for biofuels: exploiting innate microbial capacity or importing biosynthetic potential? Nat Rev Microbiol. 2009;7:715–23.
Leveau JH, Gerards S, De Boer W, Van Veen JA. Phylogeny–function analysis of (meta) genomic libraries: screening for expression of ribosomal RNA genes by large-insert library fluorescent in situ hybridization (LIL-FISH). Environ Microbiol. 2004;6:990–8.
Trachsel J, Bayles DO, Looft T, Levine UY, Allen HK. Function and phylogeny of bacterial butyryl coenzyme A: acetate transferases and their diversity in the proximal colon of swine. Appl Environ Microbiol. 2016;82:6788–98.
Preheim SP, Olesen SW, Spencer SJ, Materna A, Varadharajan C, Blackburn M, Friedman J, Rodríguez J, Hemond H, Alm EJ. Surveys, simulation and single-cell assays relate function and phylogeny in a lake ecosystem. Nat Microbiol. 2016;1:16130.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
Thang VH, Kanda K, Kobayashi G. Production of acetone–butanol–ethanol (ABE) in direct fermentation of Cassava by Clostridium saccharoperbutylacetonicum N1-4. Appl Biochem Biotechnol. 2010;161:157–70.
Wen Z, Wu M, Lin Y, Yang L, Lin J, Cen P. Artificial symbiosis for acetone-butanol-ethanol (ABE) fermentation from alkali extracted deshelled corn cobs by co-culture of Clostridium beijerinckii and Clostridium cellulovorans. Microb Cell Fact. 2014;13:92.
Al-Shorgani NKN, Kalil MS, Yusoff WMW. Biobutanol production from rice bran and de-oiled rice bran by Clostridium saccharoperbutylacetonicum N1-4. Bioprocess Biosyst Eng. 2012;35:817–26.
Simon M, Scheuner C, Meier-Kolthoff JP, Brinkhoff T, Wagner-Dobler I, Ulbrich M, Klenk HP, Schomburg D, Petersen J, Goker M. Phylogenomics of Rhodobacteraceae reveals evolutionary adaptation to marine and non-marine habitats. ISME J. 2017;11:1483–99.
Cornillot E, Nair RV, Papoutsakis ET, Soucaille P. The genes for butanol and acetone formation in Clostridium acetobutylicum ATCC 824 reside on a large plasmid whose loss leads to degeneration of the strain. J Bacteriol. 1997;179:5442–7.
Poehlein A, Cebulla M, Ilg MM, Bengelsdorf FR, Schiel-Bengelsdorf B, Whited G, Andreesen JR, Gottschalk G, Daniel R, Dürre P. The complete genome sequence of Clostridium aceticum: a missing link between Rnf-and cytochrome-containing autotrophic acetogens. MBio. 2015;6:e01168-01115.
Mitchell WJ. Carbohydrate uptake and utilization by Clostridium beijerinckii NCIMB 8052. Anaerobe. 1996;2:379–84.
Mitchell WJ. Sugar uptake by the solventogenic clostridia. World J Microbiol Biotechnol. 2016;32:32.
Servinsky MD, Kiel JT, Dupuy NF, Sund CJ. Transcriptional analysis of differential carbohydrate utilization by Clostridium acetobutylicum. Microbiology. 2010;156:3478–91.
Pyne ME, Sokolenko S, Liu X, Srirangan K, Bruder MR, Aucoin MG, Moo-Young M, Chung DA, Chou CP. Disruption of the reductive 1,3-propanediol pathway triggers production of 1,2-propanediol for sustained glycerol fermentation by Clostridium pasteurianum. Appl Environ Microbiol. 2016;82:5375–88.
Ehsaan M, Kuit W, Zhang Y, Cartman ST, Heap JT, Winzer K, Minton NP. Mutant generation by allelic exchange and genome resequencing of the biobutanol organism Clostridium acetobutylicum ATCC 824. Biotechnol Biofuels. 2016;9:4.
Xia A, Jacob A, Herrmann C, Murphy JD. Fermentative bio-hydrogen production from galactose. Energy. 2016;96:346–54.
Gao W, Francis AJ. Fermentation and hydrogen metabolism affect uranium reduction by clostridia. ISRN Biotechnol. 2013;2013:657160.
Gu Y, Li J, Zhang L, Chen J, Niu L, Yang Y, Yang S, Jiang W. Improvement of xylose utilization in Clostridium acetobutylicum via expression of the talA gene encoding transaldolase from Escherichia coli. J Biotechnol. 2009;143:284–7.
Otte B, Grunwaldt E, Mahmoud O, Jennewein S. Genome shuffling in Clostridium diolis DSM 15410 for improved 1,3-propanediol production. Appl Environ Microbiol. 2009;75:7610–6.
Sedlar K, Kolek J, Provaznik I, Patakova P. Reclassification of non-type strain Clostridium pasteurianum NRRL B-598 as Clostridium beijerinckii NRRL B-598. J Biotechnol. 2017;244:1–3.
Zhao X, Kasbi M, Chen J, Peres S, Jolicoeur M. A dynamic metabolic flux analysis of ABE (acetone-butanol-ethanol) fermentation by Clostridium acetobutylicum ATCC 824, with riboflavin as a by-product. Biotechnol Bioeng. 2017;114:2907–19.
Cai X, Bennett GN. Improving the Clostridium acetobutylicum butanol fermentation by engineering the strain for co-production of riboflavin. J Ind Microbiol Biotechnol. 2011;38:1013–25.
Gutierrez-Preciado A, Torres AG, Merino E, Bonomi HR, Goldbaum FA, Garcia-Angulo VA. Extensive identification of bacterial riboflavin transporters and their distribution across bacterial species. PLoS ONE. 2015;10:e0126124.
Wu YR, Mao A, Sun C, Shanmugam S, Li J, Zhong M, Hu Z. Catalytic hydrolysis of starch for biohydrogen production by using a newly identified amylase from a marine bacterium Catenovulum sp. X3. Int J Biol Macromol. 2017;104:716–23.
Herrmann G, Jayamani E, Mai G, Buckel W. Energy Conservation via electron-transferring flavoprotein in anaerobic bacteria. J Bacteriol. 2007;190:784–91.
Tashiro Y, Shinto H, Hayashi M, Baba S, Kobayashi G, Sonomoto K. Novel high-efficient butanol production from butyrate by non-growing Clostridium saccharoperbutylacetonicum N1-4 (ATCC 13564) with methyl viologen. J Biosci Bioeng. 2007;104:238–40.
Bahl H, Andersch W, Braun K, Gottschalk G. Effect of pH and butyrate concentration on the production of acetone and butanol by Clostridium acetobutylicum grown in continuous culture. Eur J Appl microbiol Biotechnol. 1982;14:17–20.
Li Q, Chen J, Minton NP, Zhang Y, Wen Z, Liu J, Yang H, Zeng Z, Ren X, Yang J. CRISPR-based genome editing and expression control systems in Clostridium acetobutylicum and Clostridium beijerinckii. Biotechnol J. 2016;11:961–72.
Wu YR, Zhang M, Zhong M, Hu Z. Synergistic enzymatic saccharification and fermentation of agar for biohydrogen production. Bioresour Technol. 2017;241:369–73.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
Sun Y, Zhang Y, Hollibaugh JT, Luo H. Ecotype diversification of an abundant Roseobacter lineage. Environ Microbiol. 2017;19(4):1625–38.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
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:1639–45.
Darling AE, Mau B, Perna NT. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE. 2010;5:e11147.
Konwar KM, Hanson NW, Pagé AP, Hallam SJ. MetaPathways: a modular pipeline for constructing pathway/genome databases from environmental sequence information. BMC Bioinform. 2013;14:202.
Hanson NW, Konwar KM, Hawley AK, Altman T, Karp PD, Hallam SJ. Metabolic pathways for the whole community. BMC Genom. 2014;15:619.
Xin F, Chen T, Jiang Y, Dong W, Zhang W, Zhang M, Wu H, Ma J, Jiang M. Strategies for improved isopropanol–butanol production by a Clostridium strain from glucose and hemicellulose through consolidated bioprocessing. Biotechnol Biofuels. 2017;10:118.
Miller GL. Use of dinitrosalicylic acid reagent for determination of reducing sugar. Anal Chem. 1959;31:426–8.
Lo YC, Chen WM, Hung CH, Chen SD, Chang JS. Dark H2 fermentation from sucrose and xylose using H2-producing indigenous bacteria: feasibility and kinetic studies. Water Res. 2008;42:827–42.
Trchounian K, Müller N, Schink B, Trchounian A. Glycerol and mixture of carbon sources conversion to hydrogen by Clostridium beijerinckii DSM791 and effects of various heavy metals on hydrogenase activity. Int J Hydrog Energy. 2017;42:7875–82.
Moreno R, Fierro J, Fernández C, Cuetos M, Gómez X. Biohydrogen production from lactose: influence of substrate and nitrogen concentration. Environ Technol. 2015;36:2401–9.
Wu H, Wang C, Chen P, He A-Y, Xing F-X, Kong X-P, Jiang M. Effects of pH and ferrous iron on the coproduction of butanol and hydrogen by Clostridium beijerinckii IB4. Int J Hydrog Energy. 2017;42:6547–55.
Chen W, Tseng Z, Lee K, Chang J. Fermentative hydrogen production with CGS5 isolated from anaerobic sewage sludge. Int J Hydrog Energy. 2005;30:1063–70.
Sarma S, Anand A, Dubey VK, Moholkar VS. Metabolic flux network analysis of hydrogen production from crude glycerol by Clostridium pasteurianum. Bioresour Technol. 2017;242:169–77.
Zhang J, Yu L, Lin M, Yan Q, Yang ST. n-Butanol production from sucrose and sugarcane juice by engineered Clostridium tyrobutyricum overexpressing sucrose catabolism genes and adhE2. Bioresour Technol. 2017;233:51–7.
Xin F, Wang C, Dong W, Zhang W, Wu H, Ma J, Jiang M. Comprehensive investigations of biobutanol production by a non-acetone and 1,3-propanediol generating Clostridium strain from glycerol and polysaccharides. Biotechnol Biofuels. 2016;9:220.
Malaviya A, Jang YS, Lee SY. Continuous butanol production with reduced byproducts formation from glycerol by a hyper producing mutant of Clostridium pasteurianum. Appl Microbiol Biotechnol. 2012;93:1485–94.
CS, SZ, FX, SS and YRW participated in the design of the work, organized and interpreted the data, and drafted the manuscript. CS and YRW identified two isolated strains. SZ performed the genomic analysis and comparison. CS performed the anaerobic fermentation and identified the bio-products. FX, SS and YRW participated in the data analysis and the manuscript revision. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The draft genomes of Clostridial strains WA and WB in this study are publicly deposited at DDBJ/ENA/GenBank under the Accession Numbers of PIQQ00000000 and PIQP00000000, respectively.
Consent for publication
All authors consent to publication.
Ethics approval and consent to participate
This work was financially supported by the Start-Up Funding of Shantou University (No. NTF15007), the “Sail Plan” Program for the Introduction of Outstanding Talents of Guangdong Province of China (No. 14600601) and the Major University Research Foundation of Guangdong Province of China (No. 2015KQNCX041).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.