- Open Access
Continuous production of biohythane from hydrothermal liquefied cornstalk biomass via two-stage high-rate anaerobic reactors
Biotechnology for Biofuels volume 9, Article number: 254 (2016)
The Erratum to this article has been published in Biotechnology for Biofuels 2017 10:33
Biohythane production via two-stage fermentation is a promising direction for sustainable energy recovery from lignocellulosic biomass. However, the utilization of lignocellulosic biomass suffers from specific natural recalcitrance. Hydrothermal liquefaction (HTL) is an emerging technology for the liquefaction of biomass, but there are still several challenges for the coupling of HTL and two-stage fermentation. One particular challenge is the limited efficiency of fermentation reactors at a high solid content of the treated feedstock. Another is the conversion of potential inhibitors during fermentation. Here, we report a novel strategy for the continuous production of biohythane from cornstalk through the integration of HTL and two-stage fermentation. Cornstalk was converted to solid and liquid via HTL, and the resulting liquid could be subsequently fed into the two-stage fermentation systems. The systems consisted of two typical high-rate reactors: an upflow anaerobic sludge blanket (UASB) and a packed bed reactor (PBR). The liquid could be efficiently converted into biohythane via the UASB and PBR with a high density of microbes at a high organic loading rate.
Biohydrogen production decreased from 2.34 L/L/day in UASB (1.01 L/L/day in PBR) to 0 L/L/day as the organic loading rate (OLR) of the HTL liquid products increased to 16 g/L/day. The methane production rate achieved a value of 2.53 (UASB) and 2.54 L/L/day (PBR), respectively. The energy and carbon recovery of the integrated HTL and biohythane fermentation system reached up to 79.0 and 67.7%, respectively. The fermentation inhibitors, i.e., 5-hydroxymethyl furfural (41.4–41.9% of the initial quantity detected) and furfural (74.7–85.0% of the initial quantity detected), were degraded during hydrogen fermentation. Compared with single-stage fermentation, the methane process during two-stage fermentation had a more efficient methane production rate, acetogenesis, and COD removal. The microbial distribution via Illumina MiSeq sequencing clarified that the biohydrogen process in the two-stage systems functioned not only for biohydrogen production, but also for the degradation of potential inhibitors. The higher distribution of the detoxification family Clostridiaceae, Bacillaceae, and Pseudomonadaceae was found in the biohydrogen process. In addition, a higher distribution of acetate-oxidizing bacteria (Spirochaetaceae) was observed in the biomethane process of the two-stage systems, revealing improved acetogenesis accompanied with an efficient conversion of acetate.
Biohythane production could be a promising process for the recovery of energy and degradation of organic compounds from hydrothermal liquefied biomass. The two-stage process not only contributed to the improved quality of the gas fuels but also strengthened the biotransformation process, which resulted from the function of detoxification during biohydrogen production and enhanced acetogenesis during biomethane production.
Hythane is regarded as a clean and efficient energy as it combines the advantages of both hydrogen and methane . Hythane production using biomass via two-stage anaerobic fermentation is respected as a promising direction . Compared with conventional methane fermentation, the biohythane production process improves energy recovery, reduces fermentation time, and leads to a better control of the microbial community due to the separation of the biohydrogen and biomethane processes [1–3].
Approximately 200 billion tons of lignocellulosic biomass is annually produced worldwide . The main components of lignocellulosic biomass are hemicellulose and cellulose, primarily consisting of C5 and C6 sugars, which could be used for the production of fuels and chemicals . Biohythane production using lignocellulosic biomass has been intensely investigated [3, 6, 7]. However, a long fermentation time and low gas production are observed when using lignocellulose biomass for fermentation [4, 6, 7]. This is mainly due to the natural recalcitrance of the lignocellulosic structure, which makes it difficult to directly and effectively use . Different kinds of pretreatments have been used to break down the structure of lignocellulosic biomass in order to make it less recalcitrant, including mechanical , chemical (alkaline, acid) [9, 10], biological , hydrothermal methods , or a combination of the preceding methods.
Among these, hydrothermal liquefaction (HTL) is a promising technology for the treatment and liquefaction of various biomass sources in which the water itself is an environmentally friendly solvent and reactant . The integration of HTL and anaerobic fermentation has been reported as a way that could enhance methane production from various lignocellulosic biomasses, including sunflower oil cake , sorghum forage , wheat straw [15, 16], sugar beet , rice straws , and sunflower stalks . The methane production was increased in a range of 6.5–222% in these studies [14–19]. However, there are still several challenges for successful integration. Firstly, the limited efficiency of fermentation reactors is a large bottleneck, as low HTL temperatures (100–200 °C) [14–19] in these studies resulted in a high solid content feedstock. Usually, batch reactors or conventional continuous stirred anaerobic reactors (CSTR) are used. However, these reactors are well known for their low efficiency, long retention time, and low organic loading rate. Another challenge is the inhibition of fermentation due to toxic organic compounds released during HTL. Furfural (0.08–13.32 g/L), 5-hydroxymethyl furfural (5-HMF) (0.032–4.3 g/L), and phenols (0.15–7.21 g/L) are produced during thermochemical treatment . A number of studies reported that these inhibitive compounds could be degraded during fermentation [21, 22]. However, contradictory results have also been mentioned, and biohydrogen production was suppressed by these inhibitors . Therefore, the degradation of these inhibitors during fermentation, especially two-stage biohythane systems, needs to be specifically investigated.
In this study, a novel strategy for the continuous production of biohythane from cornstalk through the integration of HTL and two-stage fermentation is proposed. HTL was conducted to get a high yield of liquid products from cornstalk. Our recent study  reported a high rate of liquefaction (up to 57.89%) and recovery of sugars and volatile fatty acids (VFAs) (up to 92.39% of aqueous products) from cornstalk via HTL. The liquid products from cornstalk after HTL were fed into the fermentation systems. The upflow anaerobic sludge blanket (UASB) and packed bed reactor (PBR) were used to build up the fermentation systems. A high density of microbes was developed in the bioreactors, further leading to an efficient performance of biohythane production [24–26]. By doing so, this study aims (1) to continuously produce biohythane from HTL liquid products via two-stage fermentation using UASB and PBR; (2) to investigate the conversion pathways of HTL liquid products, especially the inhibitors in the two-stage and single-stage processes; (3) to compare the recoveries of energy and carbon between two-stage and single-stage fermentations through batch and continuous operation; and (4) to study the structure of the microbial community during biohydrogen and biomethane production based on Illumina MiSeq sequencing.
Results and discussion
Hydrothermal liquefaction of cornstalk
The yield of aqueous phase reached up to 39.3 ± 1.8% of the dry mass of cornstalk. The carbon and nitrogen balance showed that 30% of the carbon and 58% of the nitrogen were distributed in the aqueous phase. This result indicated that the aqueous phase from the HTL process was one key stage to recover carbon and nutrients. The HTL liquid products mostly consisted of reducing sugars, VFAs, furfural, and 5-HMF (Table 1), which occupied 93.3% of all products (based on the COD). Compared with previous results , a low yield of reducing sugars was observed, probably due to the reactor scale used in this study. Compared with the 500 mL reactor in a previous study, the bigger reactor (1.8 L) used in this study had a slower heating and cooling rate. The xylose and glucose produced through the hydrolysis of the hemicellulose and cellulose were over decomposed into acids . The acids mainly consisted of acetic acid and lactic acid, which showed a similar distribution in the liquid products of hydrothermal liquefied beech wood . The inhibitors furfural and 5-HMF were produced in the treatment, and the concentration of 5-HMF was higher than furfural. This result suggested that cellulose was degraded as the 5-HMF formed from its degradation , whereas the furfural which was produced from hemicellulose was converted to acids. The produced gas mainly consisted of carbon dioxide (99.7%) and hydrogen (0.03%). The heating value (HV) of the solid residues (21.7 MJ/kg) was significantly improved compared with that of raw cornstalk (15.4 MJ/kg). The solid residue mainly consisted of lignin, as the degradation of cellulose and hemicellulose reached 100 and 50%, respectively, at this hydrothermal treatment severity . The solid residues could be used as solid fuels for combustion in a power plant .
Operational performance of continuous anaerobic reactors
Two-stage biohythane systems
The biohythane production setup consisted of two biohydrogen reactors (UASB-H, PBR-H) and biomethane reactors (UASB-M1, PBR-M1). The hydraulic retention time (HRT) for all reactors remained 12 h throughout all experiments unless elsewhere stated. A decrease of biohydrogen production in both UASB-H and PBR-H reactors was observed (Fig. 1). The hydrogen concentration in UASB-H was relatively stable (48.5 ± 8.5%) until the concentration of the HTL liquid products increased to 8 g COD/L (Additional file 1: Figure S1). However, the hydrogen concentration in PBRH increased from 24.8% in Phase 1 to 41.2% in Phase 4 (Additional file 1: Figure S1). This result was caused by the decrease of the initial pH which suppressed the hydrogen-consuming reactions in PBR-H . Biohydrogen production was undetectable when the concentration of the HTL liquid products reached 8 g COD/L (Phase 7). This result was probably caused by the low concentration of sugars in the HTL liquid products, which is the main resource for biohydrogen production during dark fermentation . In addition, the high concentrations of acetic and butyric acid in the substrate were considered to be inhibitive for biohydrogen production . What’s more, the initial concentration of 5-HMF and furfural in the feedstock increased as the loading of the HTL aqueous product increased. 5-HMF and furfural reached maximum values in Phase 7 as 142 and 15 mg/L, respectively. The presence of furan derivatives was reported to have a negative impact on biohydrogen production, which could lead to a metabolic shift from hydrogen-producing pathways (via acetate and butyrate) to non-hydrogen-producing pathways (via ethanol and lactate) . The change of potential pathway was supported by the lactate concentration in the biohydrogen reactors in Phase 7 (Fig. 1e, g), which reached values of 1080 mg/L in UASB and 1126 mg/L in PBR, respectively. The values were slightly increased, compared with the influent concentration (1024 mg/L). 5-HMF was partially degraded in both UASB-H (41.4% of the initial quantity detected) and PBR-H (41.9% of the initial quantity detected) in Phase 7, whereas most of the furfural was degraded in UASB-H (85.0% of the initial quantity detected) and PBR-H (74.7% of the initial quantity detected). One study reported a similar finding with an initial concentration below 1 g/L during biohydrogen fermentation . The result indicated that the biohydrogen process had the ability to degrade furfural and 5-HMF. COD removal in the biohydrogen reactors showed a similar change in the gas production trend. The COD removal rate in UASB-H decreased from 19.8% (Phase 1) to 2.2% (Phase 7). Similarly, this value also decreased from 21.9 to 2.8% in PBR-H.
The content of methane in UASB-M1 (73.9 ± 2.5%) and PBR-M1 (72.6 ± 3.6%) was relatively stable during operation (Additional file 1: Figure S1). The changes of pH before and after anaerobic digestion also supported a stable performance of methane production. Specifically, all the pH values of the effluents were above 7 (Fig. 1b), although the initial pH values of UASB-M1 and PBR-M1 were around 6, suggesting the acids in the HTL liquid products were converted to biogas (Fig. 1f, h). A slow decrease of the biomethane production rate was observed in the biohythane systems (Fig. 1). The methane production rates decreased to 2.53 (UASB-M1) and 2.54 (PBR-M1) L/L/day in Phase 7, respectively, corresponding to a decrease of the COD removal, at 67.2% in UASB-M1 and 68.6% in PBR-M1 (Phase 7), respectively. Most VFAs in the HTL liquid products were used for methane production, and the furfural and 5-HMF were undetectable in the effluent of the biomethane production reactors.
Single-stage biomethane systems
The single-stage systems for biomethane production (UASB-M2, PBR-M2) were set up using the same-scale reactors as the two-stage systems. The single-stage systems started with synthetic wastewater at Phase 1 and 2 with a HRT of 48 and 24 h, respectively. The total COD concentration of the influent was 8 g/L throughout all phases. The methane production in the single-stage systems showed a significant decrease when Phase 2 changed to Phase 3 (Fig. 2b), where the HTL liquid products were used instead of synthetic wastewater. By reducing the HRT from 24 to 12 h (Phase 4), methane production rates reached 2.27 in UASB-M2 and 2.07 L/L/day in PBR-M2. The COD removal decreased to 65.5 (UASB-M2) and 56.3% (PBR-M2), respectively. The methane content showed a stable performance in UASB-M2 (64.6 ± 4.9) and PBR-M2 (65.3 ± 4.1%) (Additional file 1: Figure S2). These values were lower than those in the methane reactors of the biohythane systems. Concentrations of furfural and 5-HMF were undetectable in the effluents, suggesting the complete degradation of these inhibitors in the single-stage system. The floating of granules was observed in the UASB-M2 in Phase 4 (Additional file 1: Figure S3). The floating granules accumulated around the gas–liquid–solid separator in the UASB, which probably resulted in the dysfunction of the separator and the wash out of the granules. The wash out of the granules caused by the presence of toxic compounds was also reported in a UASB used to treat the phenolic compounds . However, the detailed reason for the floating of granules awaits further investigation.
In comparison, a better performance of methane production was achieved in the two-stage process. One possible reason was the degradation of fermentation inhibitors during biohydrogen production. Moreover, the biomethane process in the biohythane systems enhanced the acetogenesis process as previously reported . The acetogenesis process is referred to as the rate-limiting procedure in the anaerobic fermentation of liquid products from HTL . This assumption was verified by the differences between the metabolic products in the effluent of the single-stage and two-stage processes. Compared with the two-stage process, the effluent of the single-stage process had a higher concentration of VFAs, mainly consisting of butyric acid, lactic acid, and acetic acid (Figs. 1f, h, 2c, d). Specifically, the concentration of lactic acid in the effluent of the PBR-M2 systems reached 773 mg/L, whereas it was undetectable in the two-stage system. The results of the batch experiments exhibited similar results. In addition, a shorter lag phase (Additional file 1: Figure S4, Table S1), higher methane production (Additional file 1: Figure S4, Table S1), and COD removal (Additional file 1: Figure S5) were observed in the two-stage batch fermentation.
Table 2 shows the performance of biohythane production through the integration of the two-stage process and HTL, compared with other studies. The biohydrogen production was limited in this study, which was also observed when using heat-treated sunflower stalks and Gelidium amansii as the feedstock [33, 36]. These studies revealed that the fermentation inhibitors produced from hydrothermal products, including 5-HMF and furfural, were supposed to change the hydrogen-producing pathway to the non-hydrogen-producing pathway. However, a hydrogen yield of 212 mL/g sugar and 109.6 and 288 mL/COD was achieved using the liquid products from pretreated switchgrass , Laminaria japonica , and wheat straw , respectively. This was probably due to the various feedstock and treatment conditions (i.e., temperature, retention time, chemicals, and reactors) which resulted in different inhibitor concentrations. The further decomposition of the produced sugars to inhibitors should be avoided. Previous studies for the hydrothermal pretreatment of lignocellulosic biomass were mostly conducted in batch reactors (Table 2), where a low heating and cooling rate may have resulted in the decomposition of produced sugars during the heating or cooling process. A continuous treatment may curtail the production of inhibitors, as the timely separation of sugars could effectively avoid their continued decomposition. Ji et al. reported a high yield of reducing sugar ratio (60.80%) and a low content of furfural in a continuous reactor . Hence, a better performance of biohydrogen production can be expected when glucose and xylose from lignocellulose are efficiently recovered under optimal HTL condition. The microbial community also plays an important role in the biogas production using HTL products. A high-rate reactor, which can retain a high density of microorganisms, seems to be more competitive. Kongjan et al. observed a higher hydrogen production rate in UASB and AF (anaerobic filter) reactor than conventional CSTR using the wheat straw hydrolysate from HTL treatment . As for the biomethane production, Table 2 shows the HRT (0.5 day) utilized in this study was much lower than previous reports (25–65 days), and a higher COD removal and methane yield were observed.
Figure 3 compares the carbon recovery (R carbon) and energy recovery (R energy) in different conversion processes. First, the integration of HTL and fermentation showed a higher carbon recovery (R carbon) and energy recovery (R energy) than direct fermentation of cornstalk. The integration of HTL and two-stage fermentation achieved a 79.0% R carbon and 67.7% R energy in continuous experiments. Batch experiments showed the potential of the HTL and two-stage fermentation, which could reach up to 84.4 and 79.0%, respectively. Second, R carbon and R energy in the two-stage process were improved in both the batch and continuous fermentation of cornstalk and HTL liquid products, compared to single-stage fermentation (Fig. 3). This result was also confirmed by a previous study . However, in the integration of the HTL and fermentation process, the solid residue from HTL contributed to the largest fraction of the total R carbon and R energy (Fig. 3). Note that the integration of HTL and two-stage fermentation did not exhibit a higher R carbon and R energy than direct fermentation if only the biogas production part is considered. The solid residues with a high heating value could be used for energy production via combustion in a power plant . In addition, the lignin-rich residue has also attracted increasing interest due to its potential to be utilized for value-added chemicals from the perspective of a biorefinery .
Although HTL is an energy intensive process, energy efficiency can be increased by recycling heat using a heat exchanger. One study reported that heat recovery could reach up to 90% . The recycled heat could be used to preheat the feedstock or maintain the temperature of anaerobic reactors . In addition, the dominant advantage of the current HTL and biohythane setup stems from its fast process (within 1 day) compared with direct fermentation of cornstalk (over 26 day) (Fig. 3). This process appears attractive for large-scale operation. The easily transportable HTL liquid products and the application of high-rate reactors would be helpful to improve the efficiency of biofuels production and reduce the capital investment and operating costs. However, detailed energy and economic analysis needs to be further evaluated.
Microbial community analysis
Scanning electron microscopy (SEM) images (Fig. 4) show the microbial morphology of granules and biofilms. As for the hydrogen reactors, the rod-shaped bacteria were dominant in both granules in UASB-H (Fig. 4a) and biofilms in PBR-H (Fig. 4b). However, a substantial amount of coccus-shaped bacteria were observed in the biofilms (Fig. 4b). As for the methane reactors, bamboo-like microorganisms were observed in the center of the methane-producing granules (Fig. 4c, e) and biofilms (Fig. 4d, f).
Illumina Miseq sequencing provided further analysis of the structure of the microbial community. Table 3 illustrates the differences in the microbial diversity. In the biohythane systems, the biohydrogen reactors (PBR-H, UASB-H) had a lower ACE, operational taxonomic units (OTUs), and Chao and Shannon indexes than the biomethane reactors. This result revealed the lower diversity of bacterial species in the biohydrogen process. Compared with PBRM2 and UASBM2, the lower ACE, OTUs, and Chao and Shannon indexes were observed in the PBRM1 and UASBM1, suggesting that the bacterial community of the methane reactors in the two-stage process had a lower diversity. However, the archaeal community showed a contrary result; the richness and diversity in the two-stage process were higher.
Further characterization of the microbial community is illustrated in Additional file 1: Figure S6. Significant differences in the microbial distribution between the biohydrogen and biomethane reactors in the two-stage process were observed (Additional file 1: Figure S6A). The biohydrogen reactors mainly consisted of the phylum Firmicutes, which occupied 99.5% in UASB-H and 94.8% in PBR-H, respectively. However, it occupied a much lower abundance in UASB-M1 (16.7%) and PBR-M1 (13.0%). The UASB-M2 and PBR-M2 had a higher abundance of Firmicutes, Proteobacteria, and Chloroflexi than UASB-M1 and PBR-M1. These bacteria were reported prevalent during the anaerobic degradation of aromatic organics, and were assumed relevant to the degradation of these inhibitors . This analysis suggested that the aromatic organics in the HTL liquid products had been degraded in UASBH and PBRH before being fed into UASB-M1 and PBR-M1. Table 4 illustrates the potential metabolic role of microbes depending on their representative species in order to further understand their microbial function. The bacterial community for hydrogen production was mainly from the family Clostridiaceae (Additional file 1: Figure S6B). The family Clostridiaceae is responsible for biohydrogen production. In addition, the family Clostridiaceae may play a very important role in the conversion process of HTL liquid products. Clostridiaceae was observed during the production of biohydrogen and acetone–butanol–ethanol at high concentrations of furfural and 5-HMF [21, 46]. Efficient conversion of cellobiose to hydrogen was observed using Clostridium sp. R1 in the presence of toxic phenolic compounds (0–1500 mg/L) . Another study proposed that the genus Clostridium was likely responsible for the conversion of phenol to benzoate, which was further degraded by acetogenesis bacteria . The family Bacillaceae and Pseudomonadaceae were found to be related to the degradation of aromatic compounds as well (Table 4). Compared with UASB-M1 and PBR-M1, UASB-H and PBR-H showed a higher Bacillaceae, Clostridiaceae, Leuconostocaceae, Planococcaceae, Pseudomonadaceae, and Streptococcaceae. Most of them are related to acidogenesis which is an important stage for the production of biohydrogen (Table 4). The higher abundance of the family Clostridiaceae, Bacillaceae, and Pseudomonadaceae revealed a detoxification function in the biohydrogen production.
In terms of the biomethane reactors, the higher abundance of potential detoxification families including Bacillaceae, Clostridiaceae, Geobacteraceae, Pseudomonadaceae, and Syntrophorhabdaceae were found in UASB-M2 and PBR-M2. Syntrophic bacteria also have an important role in the degradation of inhibitors. The genus Desulfovibrio was observed in all methane reactors and was responsible for the degradation of furfural and aromatic compounds [49, 50]. The family Syntrophorhabdaceae, which could syntrophically degrade aromatic compounds, was also found to exist in all methane reactors. However, these families were found to have a low distribution in the biohydrogen reactors (Additional file 1: Figure S6), probably due to the separation of hydrogen and methane production in the biohythane systems. The family Spirochaetaceae, which functions as an acetate-oxidizing agent , was much higher in UASB-M1 (27.9%) and PBR-M1 (10.8%) than UASB-M2 (3.3%) and PBR-M2 (3.7%). This family of bacteria enhanced the conversion of acetate to methane and strengthened acetogenesis in the methane reactors in the two-stage systems.
In addition, the archaeal community in the methane reactors (Additional file 1: Figure S6D) mainly belonged to the family Methanosaetaceae, which is an acetotrophic methanogen . Generally, UASB-M2 and PBR-M2 had a higher abundance of hydrogenotrophic methanogens. The hydrogenotrophic methanogens played an important role in the degradation of inhibitive compounds through oxidation by obligate syntrophs . The produced hydrogen from obligate syntrophs needs be consumed, which otherwise would lead to the inhibition of the end product. This result corresponds to the higher distribution of syntrophic detoxification bacteria Desulfovibrionaceae and Syntrophorhabdaceae in the UASB-M2 and PBR-M2.
This study demonstrated that biohythane production through two-stage fermentation was an attractive process for the recovery of energy and degradation of organic compounds from hydrothermal liquefied biomass. The energy and carbon recovery of the integrated HTL and continuous biohythane fermentation systems reached up to 79.0 and 67.7%, respectively. One critical challenge for biohythane production is the limited performance of hydrogen fermentation. Possible biological approaches to address this issue include the domestication of microbial consortium or upstream modification of metabolic pathways. Compared with the single-stage process, the two-stage process showed a more efficient gas production and COD removal. The two-stage process not only contributed to the improved quality of the gas fuels but also strengthened the biotransformation process due to the detoxifying function during biohydrogen production and stronger acetogenesis process during biomethane production. The higher detoxification bacteria in the biohydrogen process and the acetate-oxidizing bacteria in the biomethane process revealed by Illumina MiSeq sequencing supported the performance of gas biofuels production. The presented method might be a promising way to convert lignocellulosic biomass to biohythane.
Feedstock and HTL process
The cornstalk was collected from Golden Sun Farm (Beijing, China), and the content of cellulose, hemicellulose, and lignin in the cornstalk was 45.06 ± 0.70, 29.68 ± 0.31, and 5.65 ± 0.27%, respectively. HTL was conducted in a temperature controllable 1.8 L batch reactor (4578, Parr Instruments Co., Moline, IL, USA). In order to get an efficient degradation of cornstalk and recovery of reducing sugars and VFAs, the reaction temperature was set to 260 °C with a retention time of 0 min as previously described . The reactor was purged with nitrogen gas twice to ensure oxygen-free conditions and to maintain an initial pressure of 2.5 MPa. The liquid products were achieved by vacuum filtration of the mixture after HTL. The HTL liquid products were diluted to desired concentrations with tap water before being fed to fermentation reactors.
Biohythane and biomethane production
The UASB and PBR, which were made from transparent acrylic with a working volume of 2.5 L, were used to build up the fermentation systems. Carbon nanotubes (100 mg/L) were added into the UASB to accelerate granules formation, and polyethylene rings (1 cm diameter and 1 cm wide) were packed in the PBR. The reactors were maintained at 37 °C using a water jacket. The inoculum was obtained from the anaerobic reactor of the Xiaohongmen Municipal Wastewater Treatment Plant (Beijing, China). The inoculum of biohydrogen production was heat-pretreated before inoculation (100 °C, 15 min). The UASB and PBR were operated with synthetic wastewater before using the HTL liquid products. Regarding the synthetic wastewater, glucose was used as the carbon source, and NH4Cl served as the nitrogen source. The nutrients were added to the feedstock as previously described . NaHCO3 was added at 0.5 g/g COD for the hydrogen reactors and 1 g/g COD for the methane reactors. The pH of the substrate for biohydrogen production was controlled by adding 20 mL of 2 mol/L HCl per 1 L substrate, whereas it was not controlled for biomethane production. The reactors were operated over 200 days to enrich the microorganisms. Two series of biohythane systems were established by sequentially connecting the biohydrogen and biomethane reactors . One system consisted of two PBRs for biohydrogen and biomethane production in sequence, and the other system was composed of two UASB reactors. The concentration of HTL liquid products was increased from 0 to 8 g COD/L stepwise to replace the synthetic wastewater. Specifically, the concentration was 1 g COD/L at Phase 2, 2 g COD/L at Phase, 3 g COD/L at Phase 4, 4 g COD/L at Phase 5, 6 g COD/L at Phase 6, and 8 g COD/L at Phase 7, respectively.
The batch experiments were conducted using 250-mL glass flasks (200 mL working volume). The temperature was controlled at 37 °C by a water bath. The initial pH of biohydrogen production was controlled by adding HCl. Before biomethane production commenced, the pH was adjusted by adding 0.5 g NaHCO3. The gas was collected by gastight balloons, and the volume was measured by a syringe. All chemicals were of analytical grade and were purchased from Beijing Chemical Factory.
where P was the accumulative hydrogen or methane production (ml); P s was the hydrogen or methane production potential (mL); R m was the maximum hydrogen and methane production rate (mL/day); e is the exp(1) = 2.71828; λ was the lag time (day); and t was the incubation time (day).
where X a (%) and M a (g/mol) are the volume concentration and molar mass of the gases (hydrogen, methane, or carbon dioxide), respectively. V is the total volume of the produced gas. M cornstalk (kg) and M solid (kg) are the dry mass of cornstalk and the solid phase after HTL, respectively.
The energy yield of hydrogen (E H2, MJ/kg cornstalk) was calculated as
where Y H2 is the total hydrogen yield in HTL and fermentation (L/kg cornstalk), and HVH2 is the heating value of hydrogen (0.242 MJ/mol) .
The energy yield of methane (E CH4, MJ/kg cornstalk) was calculated as
where Y CH4 is the methane yield in HTL and fermentation (L/kg cornstalk), and HVCH4 is the heating value of methane (0.801 MJ/mol) .
The energy yield of the solid phase after HTL (E solid, MJ/kg cornstalk) was calculated as
where HVsolid is the heating value of solid phase (MJ/kg).
The energy recovery (R energy, %) was proposed to represent the ratio of produced energy to chemical energy of the feedstock. The produced energy included energy from the solid phase after HTL (E solid), biohydrogen (E H2), and biomethane (E CH4). The energy recovery can be calculated as
where HVcorstalk is the heating value of cornstalk (MJ/kg).
The heating values of the solid (HVsolid, MJ/kg) and corn straw (HVcornstalk, MJ/kg) were calculated according to the Dulong formula :
where C, H, and O are the mass percentages of carbon, hydrogen, and oxygen, respectively.
The carbon recovery (R carbon, %) was proposed to represent the ratio of carbon in the energy products (biomethane and solid phase) to the carbon in the feedstock. The carbon recovery can be calculated as
where C cornstalk, C solid, and C CH4 are the carbon content (%) in the cornstalk, solid phase, and biomethane, respectively.
Gas volume was monitored using gas meters at room temperature (25 ± 3 °C) and corrected under standard condition (273.15 K, 101.325 kPa). The gas content, including hydrogen, methane, and carbon dioxide, was determined by a gas chromatography (GC1490, Agilent Technologies, USA) equipped with a thermal conductivity detector and a stainless steel column packed with TDX-01. Nitrogen was used as the carrying gas at a flowrate of 50 mL/min. The temperature of the injector, column, and detector was 150, 120, and 150 °C, respectively. The acids, furfural, and 5-HMF were analyzed by high-performance liquid chromatography (10A, Shimadzu, Japan) equipped with an ultraviolet detector and a synergi 4u Hydro-RP (Phenomenex) column. 5 mmol/L H2SO4 was used as the mobile phase at a flowrate of 1 mL/min, and the oven temperature was 40° C. Reducing sugars were determined by the 3,5-dinitrosalicylic acid method as previously described . Element components of the cornstalk and solid residues were analyzed using a CHN analyzer (CE-440 Elemental Analyzer, Exeter Analytical, Inc. USA). The microbial morphology was observed by SEM (Quanta 200, FEI, USA) as previously described . The phylogenetic diversity of the microbial consortium was analyzed via Illumina MiSeq sequencing. Primers 515F (5′-barcode-GTGCCAGCMGCCGCGG-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′) for bacteria were used . Primers Arch344F (5′-ACGGGGYGCAGCAGGCGCGA-3′) and Arch915R (5′-GTGCTCCCCCGCCAATTCCT-3′) for archaea were used . The PCR process was conducted as previously described . Amplicons were extracted from 2% agarose gels, purified using the AxyPrep DNA gel extraction kit (Axygen Biosciences, USA), and quantified using QuantiFluor ST (Promega, USA). The purified amplicons were pooled in equimolar and paired-end sequenced on an Illumina MiSeq platform. The raw reads were deposited into the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database. The raw fastq files were demultiplexed and quality-filtered using Quantitative Insights into Microbial Ecology (QIIME). The phylogenetic affiliation of each 16S rRNA gene sequence was analyzed by a RDP Classifier (http://rdp.cme.msu.edu/) against the silva (SSU115) 16S rRNA database using a confidence threshold of 70%.
organic loading rate
upflow anaerobic sludge blanket
packed bed reactor
continuous stirred anaerobic reactors
scanning electron microscopy
volatile fatty acids
chemical oxygen demand
hydrogen production rate (L/L/day)
methane production rate (L/L/day)
hydraulic retention time
upflow anaerobic sludge blanket in two-stage system for hydrogen production
packed bed reactor in two-stage system for hydrogen production
- UASB-M1 :
upflow anaerobic sludge blanket in two-stage system for methane production
- PBR-M1 :
upflow anaerobic sludge blanket in two-stage system for methane production
- UASB-M2 :
upflow anaerobic sludge blanket in single-stage system for methane production
- PBR-M2 :
upflow anaerobic sludge blanket in single-stage system for methane production
National Center for Biotechnology
information sequence read archive
quantitative insights into microbial ecology
Liu Z, Zhang C, Lu Y, Wu X, Wang L, Wang L, et al. States and challenges for high-value biohythane production from waste biomass by dark fermentation technology. Bioresour Technol. 2013;135:292–303.
Lu F, Ji J, Shao L, He P. Bacterial bioaugmentation for improving methane and hydrogen production from microalgae. Biotechnol Biofuels. 2013;6(1):92.
Monlau F, Barakat A, Trably E, Dumas C, Steyer J, Carrère H. Lignocellulosic materials into biohydrogen and biomethane: impact of structural features and pretreatment. Crit Rev Env Sci Technol. 2013;43(3):260–322.
Monlau F, Sambusiti C, Barakat A, Guo XM, Latrille E, Trably E, et al. Predictive models of biohydrogen and biomethane production based on the compositional and structural features of lignocellulosic materials. Environ Sci Technol. 2012;46(21):12217–25.
Wettstein SG, Alonso DM, Gürbüz EI, Dumesic JA. A roadmap for conversion of lignocellulosic biomass to chemicals and fuels. Curr Opin Chem Eng. 2012;1(3):218–24.
Guo Y, Dai Y, Bai Y, Li Y, Fan Y, Hou H. Co-producing hydrogen and methane from higher-concentration of corn stalk by combining hydrogen fermentation and anaerobic digestion. Int J Hydrogen Energy. 2014;39(26):14204–11.
Kim M, Liu C, Noh J, Yang Y, Oh S, Shimizu K, et al. Hydrogen and methane production from untreated rice straw and raw sewage sludge under thermophilic anaerobic conditions. Int J Hydrogen Energy. 2013;38(21):8648–56.
Zheng Y, Zhao J, Xu F, Li Y. Pretreatment of lignocellulosic biomass for enhanced biogas production. Prog Energy Combust Sci. 2014;42(1):35–53.
Liu T, Williams DL, Pattathil S, Li M, Hahn MG, Hodge DB, et al. Coupling alkaline pre-extraction with alkaline-oxidative post-treatment of corn stover to enhance enzymatic hydrolysis and fermentability. Biotechnol Biofuels. 2014;7(1):48.
Shekiro J, Kuhn EM, Nagle N, Tucker M, Elander R, Schell D. Characterization of pilot-scale dilute acid pretreatment performance using deacetylated corn stover. Biotechnol Biofuels. 2014;7(1):23.
DeMartini JD, Wyman CE. Composition and hydrothermal pretreatment and enzymatic saccharification performance of grasses and legumes from a mixed-species prairie. Biotechnol Biofuels. 2011;4(1):52.
Li H, Chen X, Ren J, Deng H, Peng F, Sun R. Functional relationship of furfural yields and the hemicellulose-derived sugars in the hydrolysates from corncob by microwave-assisted hydrothermal pretreatment. Biotechnol Biofuels. 2015;8(1):1–12.
Jin F, Enomoto H. Rapid and highly selective conversion of biomass into value-added products in hydrothermal conditions: chemistry of acid/base-catalysed and oxidation reactions. Energy Environ Sci. 2011;4(2):382–97.
Fernández-Cegrí V, De la Rubia MA, Raposo F, Borja R. Effect of hydrothermal pretreatment of sunflower oil cake on biomethane potential focusing on fibre composition. Bioresour Technol. 2012;123:424–9.
Sambusiti C, Monlau F, Ficara E, Carr E, Re H, Malpei F. A comparison of different pre-treatments to increase methane production from two agricultural substrates. Appl Energy. 2013;104:62–70.
Chandra R, Takeuchi H, Hasegawa T, Kumar R. Improving biodegradability and biogas production of wheat straw substrates using sodium hydroxide and hydrothermal pretreatments. Energy. 2012;43(1):273–82.
Ziemi N, Ski K, Romanowska I, Kowalska-Wentel M, Cyran M. Effects of hydrothermal pretreatment of sugar beet pulp for methane production. Bioresour Technol. 2014;166:187–93.
Chandra R, Takeuchi H, Hasegawa T. Hydrothermal pretreatment of rice straw biomass: a potential and promising method for enhanced methane production. Appl Energy. 2012;94:129–40.
Hesami SM, Zilouei H, Karimi K, Asadinezhad A. Enhanced biogas production from sunflower stalks using hydrothermal and organosolv pretreatment. Ind Crop Prod. 2015;76:449–55.
Monlau F, Sambusiti C, Barakat A, Quéméneur M, Trably E, Steyer JP, et al. Do furanic and phenolic compounds of lignocellulosic and algae biomass hydrolyzate inhibit anaerobic mixed cultures?A comprehensive review. Biotechnol Adv. 2014;32(5):934–51.
Liu Z, Zhang C, Wang L, He J, Li B, Zhang Y, et al. Effects of furan derivatives on biohydrogen fermentation from wet steam-exploded cornstalk and its microbial community. Bioresour Technol. 2015;175:152–9.
Haroun BM, Nakhla G, Hafez H, Nasr FA. Impact of furfural on biohydrogen production from glucose and xylose in continuous-flow systems. Renew Energy. 2016;93:302–11.
Zhu Z, Liu Z, Zhang Y, Li B, Lu H, Duan N, et al. Recovery of reducing sugars and volatile fatty acids from cornstalk at different hydrothermal treatment severity. Bioresour Technol. 2016;199:220–7.
Si B, Liu Z, Zhang Y, Li J, Shen R, Zhu Z, et al. Towards biohythane production from biomass: influence of operational stage on anaerobic fermentation and microbial community. Int J Hydrogen Energy. 2016;41(7):4429–38.
Intanoo P, Rangsanvigit P, Malakul P, Chavadej S. Optimization of separate hydrogen and methane production from cassava wastewater using two-stage upflow anaerobic sludge blanket reactor (UASB) system under thermophilic operation. Bioresour Technol. 2014;173:256–65.
Park MJ, Jo JH, Park D, Lee DS, Park JM. Comprehensive study on a two-stage anaerobic digestion process for the sequential production of hydrogen and methane from cost-effective molasses. Int J Hydrogen Energy. 2010;35(12):6194–202.
Murakami K, Kasai K, Kato T, Sugawara K. Conversion of rice straw into valuable products by hydrothermal treatment and steam gasification. Fuel. 2012;93:37–43.
Yoshida K, Miyafuji H, Saka S. Methane production from organic acids obtained by supercritical water treatment of Japanese beech. Jpn Wood Res Soc. 2010;56(2):160–5.
Petersen MOS, Larsen J, Thomsen MH. Optimization of hydrothermal pretreatment of wheat straw for production of bioethanol at low water consumption without addition of chemicals. Biomass Bioenergy. 2009;33(5):834–40.
Shanmugam SR, Chaganti SR, Lalman JA, Heath DD. Statistical optimization of conditions for minimum H2 consumption in mixed anaerobic cultures: effect on homoacetogenesis and methanogenesis. Int J Hydrogen Energy. 2014;39(28):15433–45.
Lee H, Salerno MB, Rittmann BE. Thermodynamic evaluation on H2 production in glucose fermentation. Environ Sci Technol. 2008;42(7):2401–7.
Van Ginkel S, Logan BE. Inhibition of biohydrogen production by undissociated acetic and butyric acids. Environ Sci Technol. 2005;39(23):9351–6.
Monlau F, Aemig Q, Trably E, Hamelin JM, Steyer JP, Carrere H. Specific inhibition of biohydrogen-producing Clostridium sp. after dilute-acid pretreatment of sunflower stalks. Int J Hydrogen Energy. 2013;38(28):12273–82.
Subramanyam R, Mishra IM. Co-degradation of resorcinol and catechol in an UASB reactor. Bioresour Technol. 2008;99(10):4147–57.
Tommaso G, Chen W, Li P, Schideman L, Zhang Y. Chemical characterization and anaerobic biodegradability of hydrothermal liquefaction aqueous products from mixed-culture wastewater algae. Bioresour Technol. 2015;178:139–46.
Park JH, Yoon JJ, Park HD, Kim YJ, Lim DJ, Kim SH. Feasibility of biohydrogen production from Gelidium amansii. Int J Hydrogen Energy. 2011;36(21):13997–4003.
Veeravalli SS, Chaganti SR, Lalman JA, Heath DD. Optimizing hydrogen production from a switchgrass steam exploded liquor using a mixed anaerobic culture in an upflow anaerobic sludge blanket reactor. Int J Hydrogen Energy. 2014;39(7):3160–75.
Jung K, Kim D, Shin H. Fermentative hydrogen production from Laminaria japonica and optimization of thermal pretreatment conditions. Bioresour Technol. 2011;102(3):2745–50.
Kongjan P, Angelidaki I. Extreme thermophilic biohydrogen production from wheat straw hydrolysate using mixed culture fermentation: effect of reactor configuration. Bioresour Technol. 2010;101(20):7789–96.
Ji W, Shen Z, Wen Y. Hydrolysis of wheat straw by dilute sulfuric acid in a continuous mode. Chem Eng J. 2015;260:20–7.
Schievano A, Tenca A, Lonati S, Manzini E, Adani F. Can two-stage instead of one-stage anaerobic digestion really increase energy recovery from biomass? Appl Energy. 2014;124:335–42.
Azadi P, Inderwildi OR, Farnood R, King DA. Liquid fuels, hydrogen and chemicals from lignin: a critical review. Renew Sust Energy Rev. 2013;21(5):506–23.
Fortier MP, Roberts GW, Stagg-Williams SM, Sturm BSM. Life cycle assessment of bio-jet fuel from hydrothermal liquefaction of microalgae. Appl Energy. 2014;122:73–82.
Zhang G, Li C, Ma D, Zhang Z, Xu G. Anaerobic digestion of antibiotic residue in combination with hydrothermal pretreatment for biogas. Bioresour Technol. 2015;192:257–65.
Liang D, Fang H, Zhang T. Microbial characterization and quantification of an anaerobic sludge degrading dimethyl phthalate. J Appl Microbiol. 2009;106(1):296–305.
Zhang Y, Han B, Ezeji TC. Biotransformation of furfural and 5-hydroxymethyl furfural (HMF) by Clostridium acetobutylicum ATCC 824 during butanol fermentation. New Biotechnol. 2012;29(3):345–51.
Ho K, Chen Y, Lee D. Biohydrogen production from cellobiose in phenol and cresol-containing medium using Clostridium sp. R1. Int J Hydrogen Energy. 2010;35(19):10239–44.
Fang H, Liu Y, Ke SZ, Zhang T. Anaerobic degradation of phenol in wastewater at ambient temperature. Water Sci Technol. 2004;49(1):95–102.
Boopathy R, Daniels L. Isolation and characterization of a furfural degrading sulfate-reducing bacterium from an anaerobic digester. Curr Microbiol. 1991;23(6):327–32.
Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F. The prokaryotes: Deltaproteobacteria and Epsilonproteobacteria. 4th ed. Berlin: Springer; 2014.
Lee S, Park J, Kang H, Lee YH, Lee TJ, Park H. Distribution and abundance of Spirochaetes in full-scale anaerobic digesters. Bioresour Technol. 2013;145:25–32.
Liu Y, Whitman WB. Metabolic, phylogenetic, and ecological diversity of the methanogenic archaea. Ann NY Acad Sci. 2008;1125:171–89.
Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F. The prokaryotes: other major lineages of bacteria and the archaea. 4th ed. Berlin: Springer; 2014.
Si B, Liu Z, Zhang Y, Li J, Xing X, Li B, et al. Effect of reaction mode on biohydrogen production and its microbial diversity. Int J Hydrogen Energy. 2015;40(8):3191–200.
Zhou Y, Schideman L, Yu G, Zhang Y. A synergistic combination of algal wastewater treatment and hydrothermal biofuel production maximized by nutrient and carbon recycling. Energy Environ Sci. 2013;6(12):3765–79.
Liu Z, Lv F, Zheng H, Zhang C, Wei F, Xing X. Enhanced hydrogen production in a UASB reactor by retaining microbial consortium onto carbon nanotubes (CNTs). Int J Hydrogen Energy. 2012;37(14):10619–26.
Zheng X, Su Y, Li X, Xiao N, Wang D, Chen Y. Pyrosequencing reveals the key microorganisms involved in sludge alkaline fermentation for efficient short-chain fatty acids production. Environ Sci Technol. 2013;47(9):4262–8.
Zhou Y, Schideman L, Zheng M, Martin-Ryals A, Li P, Tommaso G, et al. Anaerobic digestion of post-hydrothermal liquefaction wastewater for improved energy efficiency of hydrothermal bioenergy processes. Water Sci Technol. 2015;72(12):2139–47.
Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F. The prokaryotes: Alphaproteobacteria and Betaproteobacteria. 4th ed. Berlin Heidelberg: Springer; 2014.
Yamada T, Sekiguchi Y, Hanada S, Imachi H, Ohashi A, Harada H, et al. Anaerolinea thermolimosa sp. nov., Levilinea saccharolytica gen. nov., sp nov and Leptolinea tardivitalis gen. nov., so. nov., novel filamentous anaerobes, and description of the new classes anaerolineae classis nov and Caldilineae classis nov. in the bacterial phylum Chloroflexi. Int J Syst Evol Microbiol. 2006;56(6):1331–40.
Rosenkranz F, Cabrol L, Carballa M, Donoso-Bravo A, Cruz L, Ruiz-Filippi G, et al. Relationship between phenol degradation efficiency and microbial community structure in an anaerobic SBR. Water Res. 2013;47(17):6739–49.
De VP, Garrity GM, Jones D, Krieg NR, Ludwig W, Rainey FA, et al. Bergey’s manual® of systematic bacteriology: the Firmicutes. 2nd ed. New York: Springer; 2009.
Pineiro SA, Stine OC, Chauhan A, Steyert SR, Smith R, Williams HN. Global survey of diversity among environmental saltwater Bacteriovoracaceae. Environ Microbiol. 2007;9(10):2441–50.
Gregoire P, Bohli M, Cayol J, Joseph M, Guasco S, Dubourg K, et al. Caldilinea tarbellica sp. nov., a filamentous, thermophilic, anaerobic bacterium isolated from a deep hot aquifer in the Aquitaine Basin. Int J Syst Evol Microbio. 2011;61(6):1436–41.
Morotomi M, Nagai F, Watanabe Y. Description of Christensenella minuta gen. nov., sp nov., isolated from human faeces, which forms a distinct branch in the order Clostridiales, and proposal of Christensenellaceae fam. nov. Int J Syst Evol Microbiol. 2012;62(1):144–9.
Cheng HH, Whang LM, Lin CA, Liu IC, Wu CW. Metabolic flux network analysis of fermentative hydrogen production: using Clostridium tyrobutyricum as an example. Bioresour Technol. 2013;141:233–9.
Lin P, Whang L, Wu Y, Ren W, Hsiao C, Li S, et al. Biological hydrogen production of the genus Clostridium: metabolic study and mathematical model simulation. Int J Hydrogen Energy. 2007;32(12):1728–35.
Zhao C, Gao Z, Qin Q, Li F, Ruan L. Desulfobaculum xiamenensis gen. nov., sp nov., a member of the family Desulfovibrionaceae isolated from marine mangrove sediment. Int J Syst Evol Microbiol. 2012;62(7):1570–5.
Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F. The prokaryotes: Firmicutes and Tenericutes. 4th ed. Berlin: Springer; 2014.
Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F. The prokaryotes: Actinobacteria. 4th ed. Berlin: Springer; 2014.
Altenschmidt U, Georg G. Anaerobic degradation of toluene in denitrifying pseudomonas sp.: indication for toluene methylhydroxylation and benzoyl-CoA as central aromatic intermediate. Arch Microbiol. 1991;156(2):152–8.
Honda T, Fujita T, Tonouchi A. Aminivibrio pyruvatiphilus gen. nov., sp nov., an anaerobic, amino-acid-degrading bacterium from soil of a Japanese rice field. Int J Syst Evol Microbio. 2013;63(10):3679–986.
Mavromatis K, Stackebrandt E, Held B, Lapidus A, Nolan M, Lucas S, et al. Complete genome sequence of the moderate thermophile Anaerobaculum mobile type strain (NGAT). Stand Genomic Sci. 2013;8(1):47–57.
Moertelmaier C, Li C, Winter J, Gallert C. Fatty acid metabolism and population dynamics in a wet biowaste digester during re-start after revision. Bioresour Technol. 2014;166:479–84.
Ren J, Yuan X, Li J, Ma X, Zhao Y, Zhu W, et al. Performance and microbial community dynamics in a two-phase anaerobic co-digestion system using cassava dregs and pig manure. Bioresour Technol. 2014;155:342–51.
Sousa DZ, Pereira MA, Smidt H, Stams AJM, Alves MM. Molecular assessment of complex microbial communities degrading long chain fatty acids in methanogenic bioreactors. FEMS Microbiol Ecol. 2007;60(2):252–65.
Djao ODN, Zhang X, Lucas S, Lapidus A, Del Rio TG, Nolan M, et al. Complete genome sequence of Syntrophothermus lipocalidus type strain (TGB-C1T). Stand Genomic Sci. 2010;3(3):267–75.
Hao J, Wang H. Volatile fatty acids productions by mesophilic and thermophilic sludge fermentation: biological responses to fermentation temperature. Bioresour Technol. 2015;175:367–73.
ZDL and BCS designed the experiment. BCS and JML operated the UASB and PBR, and did the analysis of this part. BCS, ZBZ, and JWL operated the HTL experiments and did the analysis of this part. BCS and RXS did the analysis of the microbial community. BCS wrote the manuscript and analyzed the data. ZDL, YHZ, CZ, and XHX contributed to the data analysis and revised the manuscript. All authors read and approved the final manuscript.
The authors are extremely grateful to the financial support by the National Natural Science Foundation of China (5151101167, U1562107), the National Key Technology Support Program of China (2014BAD02B03), and Beijing Youth Top-notch Talent Support Project (2015000026833ZK10). The authors thank Jimson Watson (University of Illinois at Urbana-Champaign) for improving the manuscript.
The authors declare that they have no competing interests.
Availability of supporting data
All data generated or analyzed during this study are included in this published article and its supplementary information files (Additional file 1).
Consent for publication
All authors have approved the manuscript to be published.
This work was financially supported by the National Natural Science Foundation of China (5151101167, U1562107), the National Key Technology Support Program of China (2014BAD02B03), and Beijing Youth Top-notch Talent Support Project (2015000026833ZK10).
An erratum to this article is available at http://dx.doi.org/10.1186/s13068-017-0716-1.