Skip to main content

Revealing the factors influencing a fermentative biohydrogen production process using industrial wastewater as fermentation substrate



Biohydrogen production through dark fermentation using organic waste as a substrate has gained increasing attention in recent years, mostly because of the economic advantages of coupling renewable, clean energy production with biological waste treatment. An ideal approach is the use of selected microbial inocula that are able to degrade complex organic substrates with simultaneous biohydrogen generation. Unfortunately, even with a specifically designed starting inoculum, there is still a number of parameters, mostly with regard to the fermentation conditions, that need to be improved in order to achieve a viable, large-scale, and technologically feasible solution. In this study, statistics-based factorial experimental design methods were applied to investigate the impact of various biological, physical, and chemical parameters, as well as the interactions between them on the biohydrogen production rates.


By developing and applying a central composite experimental design strategy, the effects of the independent variables on biohydrogen production were determined. The initial pH value was shown to have the largest effect on the biohydrogen production process. High-throughput sequencing-based metagenomic assessments of microbial communities revealed a clear shift towards a Clostridium sp.-dominated environment, as the responses of the variables investigated were maximized towards the highest H2-producing potential. Mass spectrometry analysis suggested that the microbial consortium largely followed hydrogen-generating metabolic pathways, with the simultaneous degradation of complex organic compounds, and thus also performed a biological treatment of the beer brewing industry wastewater used as a fermentation substrate.


Therefore, we have developed a complex optimization strategy for batch-mode biohydrogen production using a defined microbial consortium as the starting inoculum and beer brewery wastewater as the fermentation substrate. These results have the potential to bring us closer to an optimized, industrial-scale system which will serve the dual purpose of wastewater pre-treatment and concomitant biohydrogen production.


Our society today increasingly requires more energy to maintain overall ascending economic trends. Although the demand for energy is permanently growing, the reserves of our primary energy carriers will be depleted within a few decades [1]. In addition, our fossil fuel-based economy is dramatically accelerating the process of global warming with severe and permanent consequences for the environment [2]. Therefore, novel and safe energy carriers must be introduced. Hydrogen satisfies all the requirements for a clean and renewable fuel, producing only water as a by-product upon combustion or direct use in fuel-cell technology. Hydrogen has the highest energy content per unit weight of any known fuel (142 kJ/g or 61,000 Btu/lb) and can be transported for domestic/industrial consumption through conventional means [3],[4]. In addition to this, H2 gas is safer to handle than domestic natural gas, and can be used directly in internal combustion engines or in fuel cells to generate electricity [3]. Hydrogen use in fuel cells is inherently more efficient than the combustion currently required for the conversion of other potential fuels to mechanical energy [5]. However, most hydrogen is currently produced by conventional chemical or electrolytic methods, which require high amounts of energy and expensive technologies [6].

In the last few years, attention has shifted towards novel and less energy-intensive technologies for producing hydrogen [7]-[9]. Among the various hydrogen production processes, the biological methods (direct and indirect photolysis, photo-fermentation, and dark fermentation) appear to be the most promising [10]-[12]. In addition, certain methods of producing biological hydrogen such as dark fermentation can utilize various organic wastes as a substrate for fermentative hydrogen production, thus coupling organic waste treatment with renewable energy generation [13]-[16]. Recently, reducing the cost of wastewater treatment and finding ways to produce useful products from wastewater has been gaining importance with regard to environmental sustainability. One way to address both issues is to simultaneously generate bioenergy in the form of hydrogen by utilizing the organic matter present in wastewater [17]. In addition, certain types of wastewater generated by various industrial processes are considered ideal substrates because they contain high levels of easily degradable organic material [18].

The complete oxidation of glucose could yield a theoretical maximum of 12 moles of H2 per mole of glucose, but in this case no energy can be utilized to support the growth and metabolism of the hydrogen-producing organism [7],[19],[20]. Strictly anaerobic hydrogen-producing microbes are able to generate a maximum of 4 moles of H2 per mole of degraded glucose [21]-[23]. Fermentative H2 production using complex microbial communities therefore has the advantage of a high hydrogen production rate utilizing complex organic wastes as fermentation substrates with limited amounts of additional external energy input [24]. Bacteria and other microbes capable of hydrogen production exist widely in natural environments rich in organic nutrients such as soil, wastewater sludge, and compost [25]-[27]. Microbial populations sampled from these habitats can thus be used as cheap and highly efficient inocula for fermentative hydrogen production. In addition, dark hydrogen production processes using mixed microbial cultures as starting inocula are more efficient than those using pure cultures. The reason is that mixed cultures represent more simple systems to operate which are easier to control, and may be able to degrade a broader range of feedstock [28]. The use of mixed microbial cultures as starting inocula also allows the use of unsterile fermentation substrates, such as most types of wastewater. However, in a fermentative hydrogen production process using mixed cultures, the hydrogen produced by hydrogen-evolving bacteria can be consumed by hydrogen-consuming bacteria [29],[30]. Strategies for pretreatment of the parent inocula are therefore required to restrict or even terminate the methanogenic process to ensure that H2 remains the end product in the metabolic flow [31]-[33].

One of the major impediments in developing biohydrogen (bioH2) production processes for commercialization is low hydrogen production yield. The biological production rate of hydrogen and its molar yield, like most other bioprocesses, are dependent on several parameters including the activity rate of hydrogen-producing and -consuming bacteria, substrates, inorganic nutrients, and operational conditions of the bioreactor, among other factors [34]. Identifying these influencing factors and optimizing the fermentation conditions, particularly the nutritional and environmental parameters, is thus of primary importance in the bioprocess development. The most widely used screening and optimization strategy is the design of experiment (DOE) method, by which certain factors are selected and deliberately varied in a controlled manner, in order to study their effects, facilitate process comprehension, and even to improve performance [35]-[43]. The use of such DOE methods for process optimization in fermentative hydrogen production processes is critical due to the dynamic and complex nature of these systems.

In the present study, a central composite experimental design was used to investigate the influence of the process variables involved in batch-mode biohydrogen production, as well as to optimize their response. The experiments were performed using a defined mixed microbial consortium as the starting inoculum and wastewater obtained from a beer-brewing factory as the fermentation substrate. High-throughput metagenomic microbial community assessments as well as analysis of substrate degradation rates were performed during the experiments to understand the fermentation mechanisms involved. The results of this study have the potential to propel us closer to achieving and optimizing an industrial-scale system which will serve the dual purpose of providing wastewater pretreatment coupled with biohydrogen production.

Results and discussion

Batch experiments were conducted using mixed microbial consortia able to degrade complex organic substrates like wastewater in addition to simultaneous biohydrogen production. The involvement of biological, chemical, and physical factors influencing the fermentative biohydrogen production process as well as the interactions of these parameters were assessed by experiments designed by statistics-based methods and metagenomic monitoring of this complex ecosystem.

Effect of various influencing factors and the interactions between them on the biohydrogen production process

To assess the influence of various factors on the biohydrogen production process, as well as the level of interaction between these factors, a cybernetic representation of the dark biohydrogen fermentation process was developed, with an input-output structure (Figure 1). The input data are defined as influencing factors (IF) while the system output data are defined as objective functions (OF). By developing and applying this experimental modeling approach, connecting relationships (most frequently with a predefined polynomial shape) between the IF and the OF (in this case, H2 production rates) can be identified for an experimental domain previously defined as being of interest. The data obtained by applying this experimental design strategy can be used to generate predictions regarding the behavior of the system under investigation. They can also provide deep insight into the degree of influence of the investigated variables on the objective functions. In addition, by applying advanced statistical approaches, the optimum region of a system for a specific response can be identified (in this case, biohydrogen production rate).

Figure 1

Cybernetic representation of the investigated dark biohydrogen fermentation process with an input-output structure. The anaerobic fermentation of the wastewater is considered a black-box system influenced by different variables, with consequences on the objective function (in this case, bioH2 production).

Because of the steep curve generated by the response surface approach during the investigation of an optimal area within a specific system, using linear methods to adequately identify the stationary point becomes obsolete. Therefore, depending on the number of investigated variables, further additional experimental points are required. One of the most recommended approaches to address this issue is the central composite experimental design method. By applying this strategy, one can use the information obtained through a first-order mathematical model, completing these insights by consecutively extending the experimental approach to explain a second-order mathematical model. Thus, a central composite design (CCD) consists of a standard first-order design with a 2kor 2k-p orthogonal factorial point, where k is the number of selected factors and p the number of interactions replaced with influencing factors. 2k axial points are displayed in a "star pattern" at distance α from the center of the design with n0 center points. Based on preliminary screening experiments, three IF were selected; fermentation temperature, starting pH value, and degradable substrate availability (glucose addition; Table 1). Each IF was defined between two physical values, low and high, coded as -1 and +1 respectively, with a center value coded as 0.

Table 1 Physical and coded values of the variables used in the central composite design

The relationship between the physical and coded values is given by the following equation:

x jcod = x jphys - x j 0 phys I jphys ,j=1,2,3


xjcod is the coded value of the j factor

xjphys is the physical value of the j factor

xj0phys is the central level of the j factor

Ijphys is the variation interval for the j factor; double this value gives the variation range D (Equation 2)

D= x jphys high - x jphys low =2 I jphys

The factorial portion of the CCD is a complete 23 factorial with eight runs, which contains all the possible combinations within the defined levels of the investigated variables (runs 1-8; Table 2). The additional experimental runs (9-14) represent supplementary axial points displayed in a "star pattern" around the center of the design, at distance α of 1.287 from the center. The design is completed with n0 = two observations at the experimental center (runs 15 and 16). Because of the different nature of the experimental design approaches (runs 1-8 and runs 9-14 together with runs 15 and 16), we will separately discuss the results obtained for these two main experimental groups. Biohydrogen production was monitored every 24 h for the duration of all experiments.

Table 2 Central composite experimental design matrix of the three investigated variables, with the total measured H 2 production for each of the experimental runs

Even though the same starting microbial inocula and the same wastewater were used, notable differences in biohydrogen production rates were observed for the factorial portion of the CCD (runs 1-8; Figure 2A). This is the first indication that the selected variables manifest, beyond doubt, a clear influence on the OF in the investigated system. The highest biohydrogen production rates observed during the factorial portion of the CCD were measured during experimental run 4 (a total production of 25.2 mL h2 was measured at the end of the experiment), followed closely by experimental run 3 (a total production of 22.3 mL h2 was measured at the end of the experiment). The other six experimental runs generated considerably smaller amounts of H2, ranging from 14.6 mL (experimental run 7) to 8.2 mL (experimental run 5) H2 at the end of the experiment (Figure 2A). Strong similarities in the H2 production rates were observed between similar experimental conditions (Table 2). This was particularly true for most of the cases, where only the glucose addition variable (X3) differed.

Figure 2

Biohydrogen evolution measured during the full factorial (A) and additional orthogonal central composite (B) multifactorial experiments. Clear differences can be observed between the different experimental lines, suggesting a strong influence of the investigated variables on the OF.

Notable differences in volumetric biohydrogen production were also observed for the supplementary star and central points (experimental runs 9-16; Figure 2B). The highest biohydrogen production rates were obtained for experimental runs 12 and 14 (a total H2 production of 27.1 mL and 25.9 mL was measured at the end of the experiments respectively), while the lowest biohydrogen production rate was obtained for experimental run 11 (a total H2 production of 2.9 mL was measured at the end of the experiment). Regarding the fermentation conditions of the highest hydrogen-producing experimental run (12), the fermentation temperature (X1) as well as the glucose addition (X3) values were situated in the center of the experimental model, while the initial pH (X2) was fixed at a value of 6.74, representing a distance of 1.287 from the experimental center (Table 2). These preliminary findings suggest that the initial value of the environmental pH has a strong influence on the biohydrogen production rate.

The analysis of variance (ANOVA) test performed on the central composite experimental design showed that both linear and quadratic effects of the temperature and pH, as well as their first-order interactions, are statistically significant parameters with regard to the OF (Table 3). However, the main influence of initial glucose addition (X3), its quadratic effect, and first-order interactions on biohydrogen production are regarded as statistically insignificant. The lack of fit has a P-value slightly higher than α = 0.05, which makes it marginally significant (meaning that it is significant at the 0.10 α - level).

Table 3 ANOVA test performed on the second-order polynomial model used to discriminate between the significant linear (L) and quadratic (Q) effects of the investigated variables, as well as their interactions, for the investigated system and response (biohydrogen production)

The estimated effects of the IF on the OF together with their order of magnitude were revealed through additional data analysis (Figure 3A). The largest effect on biohydrogen production for the investigated system was caused by a switch in initial pH value from lower to higher values. Fermentation temperature also had a large effect on the OF, followed by the interaction between the fermentation temperature (X1) and initial pH value (X2). The interaction plot for X1X2 indicates a considerable increase in biohydrogen production at the point when initial pH values move from low to high levels, and temperature is low (Figure 3B). This demonstrates that understanding the interactions among the various process parameters in complex systems is crucial for developing and optimizing an efficient biohydrogen production process.

Figure 3

The effects and the interactions of the investigated variables on the biohydrogen production process. Pareto histogram depicting the estimated linear (L) and quadratic (Q) effects of each of the analyzed independent variables (in decreasing order of magnitude) on the biohydrogen production rate (A). Analysis of the effect of the interactions between the X1 and X2 variables on biohydrogen production rate (B).

By analyzing the evolution of the main effects throughout the experimental process (120 h), an explicit shift in the direction and intensity of the effect of these variables on the biohydrogen production rate was detected (Figure 4). The main effects of initial pH (X2) and temperature (X1) tended to increase over time, while the effects of other variables, like glucose addition (X3), tended to maintain a relatively constant value. These findings emphasize even more the dynamic characteristics of such a complex environment. Thus, the initial environmental conditions as well as the fermentation process settings greatly influenced the development of microbiological processes, and thereby the metabolic end products, in the investigated system. Because of the differences in the direction and intensity of the effects that investigated variables exerted on the objective functions, multiple optimal situations may be identified depending on the point selected in the fermentation development process. Consequently, application of these findings to an industrial-scale continuous biohydrogen fermentation system is crucial for establishing an optimum process development strategy with regard to influent feeding of the system among other factors. These findings particularly support the necessity for a complete exploration of the effects manifested over time on the biohydrogen production process by the different influencing factors.

Figure 4

The evolution of the main effects of the investigated variables and their interactions on the biohydrogen production rate. The three investigated variables as well as the statistically significant interactions between the operating temperature and the initial pH are represented.

These crucial insights show that in complex biotechnological systems, both the effect of a particular variable at a specific time on the OF and the global effect of multiple influencing factors together with their different interactions must be considered. In addition, the evolution of the main effects of the IF throughout the biohydrogen fermentation process must be considered. Using these insights, an optimized industrial-scale biohydrogen production system can be successfully designed and operated in a feasible economic context.

Identification of the optimum conditions for biohydrogen production in the investigated system

Once the estimated effects on biohydrogen production rates and their magnitudes were established for each of the analyzed independent variables as well as for their interactions, a strategy to identify the optimum conditions of the investigated system to maximize hydrogen production rates was developed. To understand the specific hydrogen production conditions and to confirm the validity of the statistical experimental strategies applied, a response surface and contour plot methodology was developed. Prediction of H2 production at any tested parameter within the range of the applied experimental design is achieved by employing a second-order polynomial regression equation obtained from experimental data (Equation 3 and Figure 5). With respect to the coded factor levels, the second-order model used to fit the experimental data is:

Y=20.95-2.83 X 1 +6.22 X 2 +1.35 X 3 -2.42 X 1 X 2 -0.82 X 1 X 3 -0.11 X 2 X 3 -3.97 X 1 2 -3.35 X 2 2 -0.12 X 3 2
Figure 5

Modeling the biohydrogen production using beer-brewing wastewater as the fermentation substrate. Observed versus predicted values of the biohydrogen production process using a complex microbial consortium as the starting inoculum and beer-brewing wastewater as the fermentation substrate.

with 60% of the variation in biohydrogen production explained by the model (R2adj = 0.60).

When considering the observed versus predicted values of biohydrogen production, it can be assumed that the statistical model developed is reasonably accurate. The response surface and contour plot analysis was developed for an initial glucose concentration of 10 g/L (Figure 6). The stationary point was found to be outside the investigated experimental domain. The canonical analysis of the response surface indicates a possible rising ridge, since two of the canonical coefficients are close to zero (-1 = -0.002995; -2 = -0.099786; -3 = -4.65529). In this type of ridge system, inferences about the true surface of the stationary point cannot be drawn because it is outside the region where the model has been fitted. The contour plots for three levels of X3 (initial glucose concentration; low, medium, and high) display an optimum around X1 = 27°C, X2 = 6.7, which also indicates that the influence of X3 on biohydrogen production in the investigated experimental domain is limited. Therefore, if the initial glucose addition is pooled out from the model (according to Equation 4), considering that it was found statistically insignificant in the investigated domain, the stationary point is found at X1 = 26.7°C and X2 = 6.66, with a predicted optimal value of biohydrogen production of 25.57 mL with a ±95% confidence interval.

Y=-343.22+9.05 X 1 +74.42 X 2 -0.47 X 1 X 2 -0.11 X 1 2 -4.64 X 2 2
Figure 6

Prediction of the optimum area for the highest biohydrogen production yields. Response surface (A) and contour plot (B) analysis of H2 production as a function of temperature and pH, with a constant glucose addition value of 10 g/L.

where the coefficients are calculated with uncoded (natural) factor values.

The stationary point was verified by carrying out experiments in triplicate. The mean biohydrogen production obtained was 28 mL, and this value was within the ±95% confidence interval for the predicted maximum. Therefore, the new empirical model, with a determination coefficient of R2adj = 0.62 and from which all the terms containing the X3 independent variable have been pooled out, was able to predict the behavior of the system within the experimental domain (Figure 7).

Figure 7

Confirmation of the predicted stationary point for the investigated system. The mean of the maximum biohydrogen production yield of 28 mL measured at the end of the confirmation experiments is within the ±95% confidence interval for the predicted optimum area. *The microbial population used as a starting inoculum during the wastewater degradation experiments was subjected to heat pretreatment prior to the inoculation.

Wastewater degradation

During the dark fermentative biohydrogen production process, organic matter is converted from complex long-chain molecules to simple compounds. By using wastewater as a fermentative organic substrate, at least partial biodegradation of this waste can therefore be achieved. Mass spectrometry (MS) was used to evaluate the wastewater biodegradation efficiency. Several organic compounds including lactose, glucose, acetic acid, propionic acid, and furfurol, among others, were identified and monitored during the wastewater degradation experiments. The results allowed direct comparison of the experimental runs by revealing the different degradation rates and metabolic pathways followed by the microbial communities involved.

The analysis of three different experimental situations, one with low hydrogen production (experimental run 11), one with medium hydrogen production (experimental run 13), and one with high hydrogen production (experimental run 14), revealed significant differences in the wastewater composition at the end of the biohydrogen fermentation process (Figure 8). As expected, comparison of low and high hydrogen-producing experimental runs revealed a general decrease in concentration of most of the measured components in the fermented wastewater with increasing hydrogen production. This tendency was especially marked for lactose, glucose, capric acid, maltose, lactic acid, furfurol, and caproic acid (Figure 8). These observations suggest that most of the macromolecules in the wastewater are metabolized, with biodegradation efficiency clearly correlating with the amount of hydrogen generated. Some components (glucose, maltose, and furfurol) were fully consumed during fermentation, suggesting that these compounds represent easily accessible energy sources which drive the hydrogen-producing fermentation pathways. At the same time, accumulation of propionic acid and galactose was observed concomitant with the increasing hydrogen yield (Figure 8).

Figure 8

Wastewater chemical composition at the end of the fermentation period. Comparison between experimental runs which generated low (experimental run 11), medium (experimental run 13), and high (experimental run 14) hydrogen production rates.

The clear differences between experimental runs observed during the biohydrogen-producing wastewater degradation experiments indicate a strong correlation between the consumption of the organic constituents and hydrogen production rates. This suggests that the technological process developed here has the potential to generate biohydrogen using wastewater as an organic substrate.

Assessment of microbial community composition

The same mixed microbial consortium was used as a starting inoculum in all experimental runs and stages. This consortium was designed according to our previous research on the ability of various complex consortia to degrade complex organic substrates and simultaneously produce biohydrogen [44]. An organic nutrient-rich habitat (generated during the biological denitrification step at a communal wastewater treatment plant) was sampled in order to isolate the required microbial inoculum. To further increase the biohydrogen-generating abilities of the population, heat pretreatment was applied. As a result, the majority of the spore-forming hydrogen-producing bacteria survived the treatment, while the hydrogen-consuming methanogens were largely eliminated.

To better understand the hydrogen-producing dark fermentation processes, a sequencing-based metagenomic analysis was carried out on samples selected from different experimental runs. This approach allows for a detailed taxonomic and functional characterization of the selected ecosystem as a function of time, which is essential when planning strategies for dark fermentation-based biohydrogen production using wastewater as a fermentative substrate.

To fully understand the experimental model, the complete factorial portion of the CCD was subjected to metagenomic analysis. Additional samples from the center as well as axial experimental points were analyzed metagenomically, with the microbial composition of samples with the lowest and highest H2-producing capabilities being investigated. This approach provided valuable insight into the transitions and structural rearrangements occurring in the microbial communities under different fermentation conditions. The influence of these rearrangements on the OF (biohydrogen production rates) was also elucidated.

The metagenomes extracted from the full factorial experimental design approach (experimental runs 1-8) revealed similarities in the microbial populations collected from the different experimental setups (Figure 9). Two major clusters formed at the genus level. Cluster A contained experimental runs 2, 3, 4, 6, and 7, while cluster B contained experimental runs 1, 5, and 8. In addition, both clusters could be further divided into smaller units. Most of the cluster A microbial communities developed under similar experimental conditions with regard to the initial pH value (X2) and glucose addition (X3). A notable difference between these clustered groups was the fermentation temperature (X1). It seems that in the temperature range investigated (25 to 37°C) the microbial composition in each experimental run was similar at the same levels of X2 and X3, regardless of X1 levels. Interestingly, the microbial communities located at the closest Bray-Curtis distance (experimental runs 5 and 8), which are isolated from the experimental runs in cluster B, developed under the same temperature conditions (37°C) and different levels of X2 and X3. However, even in this cluster, the microbial community isolated from experimental run 1 developed at the same levels of X2 and X3 but at levels of X1 that differed from the community in experimental run 5. These findings suggest a strong correlation between the microenvironmental conditions inside the bioreactor and the different developmental pathways followed by the microbial community.

Figure 9

Partial heatmap calculated for samples from the full factorial experimental design approach (experimental runs 1-8). The heatmap was redrawn using normalized values clustered by genus using a Bray-Curtis distance metric (divided into clusters A and B for practical reasons).

The cluster formed by samples from experimental runs 3, 4, and 7 displayed the highest H2 production potential (22.24 mL, 25.23 mL, and 14.67 mL, respectively). Certain key microorganisms identified in this cluster, such as Aeromonas spp., Clostridium spp., Thermoanaerobacter spp., and Alkaliphilus spp., were generally more abundant in experimental runs 3 and 4, while Bacillus spp. and Lactobacillus spp. were poorly represented in these runs (Figure 9).

Samples from the supplementary points (star and central points) in the experimental matrix (experimental runs 9-16) were also analyzed using metagenomics, with samples producing the lowest and highest amounts of H2 being compared (Figure 2B and Figure 10). Strong correlations were observed between the structural changes within these microbial communities and their biohydrogen-producing potential. The microbial community with the lowest H2 production capacity was dominated by Lactobacillus spp. and Tetragenococcus spp. (Figure 10A). The microbial consortium which produced moderate quantities of H2 was dominated by Citrobacter spp. and Aeromonas spp. (Figure 10B). The microbial population which produced the highest amounts of H2 was clearly dominated by Clostridium spp. (Figure 10C). A clear shift from a Lactobacillus spp.-dominated microbial population towards a Clostridium spp.-dominated population concomitant with increasing H2 production rates was thus observed. Interestingly, the only factor that varied between these three experimental runs was the initial pH value (X2; Table 2). Therefore, changing only one IF in the investigated system can result in significant changes in the OF (the biohydrogen production rate).

Figure 10

Microbial population shifts during the biohydrogen production experiments. Microbial community composition of samples from experimental runs 11 (A), 16 (B), and 12 (C), part of the central composite fractional factorial experimental design approach.

The results obtained from metagenomic investigations have led to a better understanding of the impact of each physicochemical parameter and the interactions between them, all of which are essential for the fermentative biohydrogen production process. To maximize the biohydrogen production potential of a microbial consortium during wastewater degradation, special emphasis has to be made not only on the careful selection and pretreatment of these populations but also on the optimal fermentation conditions for each system.


A central composite experimental design was used to underline the involvement of the various biological, chemical, and physical factors influencing the fermentative biohydrogen production process. The experiments were performed using a mixed microbial consortium as the starting inoculum and beer-brewing wastewater as the fermentation substrate.

It was demonstrated that the selected variables have a clear influence on the objective function in the investigated system. Both linear and quadratic effects of the fermentation temperature and initial pH, as well as their first-order interactions, were statistically significant with regard to the biohydrogen production rates for the system considered. The largest effect was caused by a change in initial pH value from its lower to its higher level. The fermentation temperature also had a strong effect on the objective function, followed by the interaction between the fermentation temperature (X1) and initial pH value (X2). Analysis of the evolution of these main effects over time during the experiment revealed a strong shift in the direction and intensity of the effect of these variables on the biohydrogen production rate. Because of this, further optimal situations can be identified depending on the point in the fermentation process being investigated. These crucial insights show that this complex biotechnological system is governed by the combined effect of several influencing factors as well as by their interactions with each other.

We successfully used response surface and contour plot methodology to understand the specific hydrogen production conditions of the system under study, and to confirm the validity of the statistical experimental strategies applied. In addition, confirmation experiments yielded a mean biohydrogen production of 28 mL, a value situated within the ±95% confidence interval of the predicted maximum. The new empirical model was thus able to predict the system behavior within the experimental domain.

The mass spectrometry analyses performed on the degraded wastewater during each biohydrogen production experimental run revealed significant differences in the biochemical composition between runs, confirming a strong correlation between the consumption of most of the organic constituents and a high level of hydrogen generation. This indicates the potential of the developed technological process to generate biohydrogen by using wastewater as an organic substrate.

To better understand the microbiological factors driving the dark fermentative biohydrogen production processes in the system studied, high-throughput next generation metagenomic analyses were carried out on samples taken from different experimental runs. These analyses revealed strong correlations between the micro-environmental conditions inside the bioreactors and the developmental pathways taken by the microbial communities, even though the same microbial consortium was used as a starting inoculum in all experimental runs. By analyzing the metagenomes of the lowest and highest H2-producing samples, population shifts in hydrogen production potential within the microbial consortia triggered by the different fermentation conditions can be traced. These results have led to a better understanding of the impact which a mixed microbial consortium has on the fermentative biohydrogen production process, and at the same time, the influence of the factors involved on the development of these populations over time.

Together the data generated by the present study reveal a strong interconnection between the investigated variables, as well as a permanent and shifting influence on the biohydrogen production process under the investigated conditions. Only by understanding these phenomena can an optimized industrial-scale biohydrogen production system be successfully designed and operated in a feasible economic context. This will bring us one step closer to a clean, fossil fuel-free future and a developed H2-based economy.


Seed inocula

Biological samples were collected from the denitrification step at a municipal sewage wastewater treatment plant. This ecosystem has a high microbial biodiversity consisting of naturally formed populations of microflora suitable for biodegradation of complex organic substrates. Once collected, the samples were stored at 4°C until inoculation.

Prior to inoculating the bioreactors, microbial samples were subjected to a heat pretreatment process at 70°C for one hour as described in our previous work [44]. The purpose of the pre-treatment strategy is to enrich with spore-forming hydrogen-producing microbes as well as to reduce the abundance of hydrogen-consuming microorganisms.

Bioreactor design and operation

The experimental setup was conducted using wastewater generated by the beer-brewing industry as the fermentation substrate in different experimental combinations as detailed in Table 2. The composition of the wastewater was: total COD - 6558 mg/L, soluble COD - 5066 mg/L, total N - 75.8 mg/L, and total P - 58.4 mg/L. The batch-mode experiments were conducted in 100-mL serum vials with 50 mL of wastewater and 10 mL pretreated sediment samples as a starting inoculum. A granular biofilm support material (sterilized river-bed rocks) was introduced into the bioreactors to increase the contact surface area of the microbial populations. The bottles were capped with rubber septum stoppers and aluminum rings. Incubation was performed at varying temperature conditions for a period of 120 h (Table 2). During the fermentation period, response parameters including biogas production and composition, metabolite concentration, substrate degradation, and total microbial community composition were monitored as described below. All batch experiments were performed in triplicate.

Statistical experimental design methods

A central composite experimental design approach was used to assess the IF and their effect on the OF (biohydrogen production rate). In the early stage of the experiments, a statistically based full factorial experimental design approach was utilized to determine the influence of the factors under study (fermentation temperature, starting pH value, and glucose addition) on the response (H2 production; Table 1). Use of this strategy helped to avoid overlap of different effects and interactions among these variables. The screened factors were chosen based on our previous work and were tested at low, medium, and high levels, coded as -1, 0, and +1 (Table 1) [44]. The factorial portion of the design is a complete 23 factorial with eight runs, which contains all the possible combinations within the defined levels of the investigated variables (runs 1-8; Table 2). In addition, a second statistically based factorial experimental design method was applied in order to obtain a complete understanding of the investigated process. These experimental runs (9- 14) represent additional axial points displayed in a "star pattern" around the center of the design, at α distance of 1.287 from the center, a value that ensures its orthogonality. The design also contains two observations at the experimental center (runs 15 and 16; Table 2). This approach provided N0 - 1 degrees of freedom in estimating the experimental error, and at the same time established the estimation precision of the OF (biohydrogen production rate) around the central experimental point. The biohydrogen production rate was monitored every 24 h for the duration of the experiments. Experimental design approaches were developed and analyzed using the Statistica 8 software suite (StatSoft Inc., USA).

Analytical methods

The quantity and composition of the biogas produced was directly measured by gas chromatography in 24-h intervals using an Agilent Technologies 7890A GC system equipped with a thermal conductivity detector and argon as a carrier gas. The temperatures of the injector, detector, and column were kept at 30°C, 200°C, and 230°C, respectively. An HP PLOTQ column (15 m × 530 mm × 40 mm) was used. Since a concentration gradient of H2 gas can form in the headspace, gas samples (0.5 mL) were taken out after mixing of the headspace gas by sparging several times with a gas-tight syringe.

Microbial degradation of wastewater was monitored at the end of the experimental runs by mass spectrometry conducted on a High Capacity Ion Trap Ultra mass spectrometer (HCT Ultra, PTM discovery; Bruker Daltonics, Bremen, Germany). All mass spectra were acquired in the mass range 100 to 3000 m/z, with a scan speed of 2.1 scans/s. Tandem mass spectrometry was carried out by collision-induced dissociation (CID) using He as the collision gas. For MS/MS sequencing, precursor ions were selected within an isolation width of 2 μm. The fully automated process was performed using a NanoMate 400 robot incorporating ESI chip technology (Advion BioSciences, Ithaca, USA) couplet on a High Capacity Ion Trap Ultra mass spectrometer (HCT Ultra, PTM discovery; Bruker Daltonics, Bremen, Germany). The robot was controlled and manipulated by ChipSoft software operated under the Windows system. The position of the electrospray chip was adjusted to the sampling cone potential to give rise to an optimal transfer of the ionic species into the mass spectrometer. To avoid contamination, a glass-coated microtiter plate was used in all experiments. Five microliter aliquots of working sample solutions were loaded onto a 96-well plate. The robot was programmed to aspirate the whole sample volume into the pipette tip followed by 2 μL of air, and then to deliver the sample to the inlet side of the microchip. Each nozzle had an internal diameter of 2.5 μm, and under the given conditions delivered a flow rate of approximately 200 nL/min. The nano-ESI process was initiated by applying voltages of 1.5 to 1.8 kV and a head pressure of 0.5 to 0.7 PSI. After spray initialization, infusion parameters (ESI voltage in the pipette tip, voltage, and desolvation gas flow) were optimized.

Values for ESI capillary, cone potential, and desolvation gas (nitrogen) were optimized to achieve an efficient ionization and to produce the optimum transfer of ions during MS. Measurement parameters were: capillary voltage of 1 kV, contraelectrode voltage (cone voltage) of 60 V, acquisition time 2 min, scan speed 2.1 scans/s, and a mass range of 100 to 3000 m/z. The NanoMate HCT MS system was tuned to operate in the positive ion mode. This technique was chosen because glucide ionization is highly efficient in this mode. The source block maintained at a constant temperature of 80°C provided optimal desolvation of the generated droplets without the need for desolvation gas. To prevent any cross-contamination or carry-over, the pipette tip was ejected and replaced after every sample infusion and MS analysis. All mass spectra were processed using Data Analysis 3.4 software (Bruker Daltonik, Bremen, Germany). The mass spectra were calibrated using sodium iodide. Accurate determination of the average mass was 20 ppm. Samples were dissolved in methanol at a concentration of approximately 5 pmol/μL. At the acquisition time of 2 min, the required volume of sample was approximately 2 pmol, a value indicating a highly sensitive analysis.

Total DNA extraction from samples

DNA from the complex samples was extracted and purified according to described methods with some modifications [45]. Samples (0.5 g) were extracted with 1.3 mL extraction buffer (100 mM Tris-Cl pH 8.0, 100 mM EDTA pH 8.0, 1.5 m NaCl, 100 mM sodium phosphate pH 8.0, 1% CTAB). After thorough mixing, 7 μL of proteinase K (20.2 mg/mL) was added. After incubation for 45 min, 160 μL 20% SDS was added and mixed by inversion several times with further incubation at 60°C for 1 h with intermittent shaking every 15 min. Samples were centrifuged at 13,000 RPM for 5 min, and the supernatant was transferred into new Eppendorf tubes. The remaining soil pellets were treated three times with 400 μL extraction buffer and 60 μL SDS (20%) and kept at 60°C for 15 min with intermittent shaking every 5 min. Supernatants collected from all four extractions were mixed with an equal quantity of chloroform and isoamyl alcohol (25:24:1). The aqueous layer was separated and precipitated with 0.7 vol isopropanol. After centrifugation at 13,000 RPM for 15 min, the brown pellets were washed with 70% ethanol, dried at room temperature, and dissolved in TE (10 mM Tris-Cl, 1 mM EDTA, pH 8.0).

Metagenomic characterization of microbial communities

The total DNA from selected samples was prepared for high-throughput next generation sequencing analysis performed on the Ion Torrent PGM platform (Life Technologies). An average of 291.322 sequencing reads were generated for each sample, with a mean read length of 161 nucleotides. Bioinformatic analyses (taxonomic profiling and assessment of metabolic potential) were conducted using the public MG-RAST software package, which is a modified version of RAST (Rapid Annotations based on Subsystem Technology) [46]. The sequence data were compared to M5NR using a maximum e-value of 1 × 10-5, a minimum identity of 95%, and a minimum alignment length of 15, measured in amino acids for proteins and base pairs for RNA databases.

Authors' contributions

IZB participated in the conception, design, experimental work, and data collection and analysis, and also drafted the manuscript. MI participated in the design of the experiment approaches and mathematical modeling. VDG carried out the microbiological sampling and pretreatment procedures. IM conceived the design and operation of the bioreactors. BP elaborated and carried out the total DNA extractions as well as community composition assessments by the 16S rRNA method and metagenomic approaches. AN designed and performed the wastewater analysis. EK participated in the critical discussion of the results. TB contributed to the evaluation of the analytics data. GM performed the metagenomics experiments, designed the study, and participated in the critical discussion of the results. All authors read and approved the final manuscript.





central composite design


chemical oxygen demand


design of experiments


influencing factors


objective functions


Rapid Annotations based on Subsystem Technology


  1. 1.

    Matsakas L, Kekos D, Loizidou M, Christakopoulos P: Utilization of household food waste for the production of ethanol at high dry material content. Biotechnol Biofuels. 2014, 7: 4-10.1186/1754-6834-7-4.

    Article  Google Scholar 

  2. 2.

    Dziugan P, Balcerek M, Pielech-Przybylska K, Patelski P: Evaluation of the fermentation of high gravity thick sugar beet juice worts for efficient bioethanol production. Biotechnol Biofuels. 2013, 6: 158-10.1186/1754-6834-6-158.

    Article  Google Scholar 

  3. 3.

    Das D: Advances in biohydrogen production processes: an approach towards commercialization. Int J Hydrog Energy. 2009, 34: 7349-7357. 10.1016/j.ijhydene.2008.12.013.

    Article  Google Scholar 

  4. 4.

    Sevilla M, Mokaya R: Energy storage applications of activated carbons: supercapacitors and hydrogen storage. Energy Environ Sci. 2014, 7: 1250-1280. 10.1039/c3ee43525c.

    Article  Google Scholar 

  5. 5.

    Levin DB, Chahine R: Challenges for renewable hydrogen production from biomass. Int J Hydrog Energy. 2010, 35: 4962-4969. 10.1016/j.ijhydene.2009.08.067.

    Article  Google Scholar 

  6. 6.

    De Vrije T, Budde MA, Lips SJ, Bakker RR, Mars AE, Claassen PA: Hydrogen production from carrot pulp by the extreme thermophiles Caldicellulosiruptor saccharolyticus and Thermotoga neapolitana. Int J Hydrog Energy. 2010, 35: 13206-13213. 10.1016/j.ijhydene.2010.09.014.

    Article  Google Scholar 

  7. 7.

    Ghosh D, Bisaillon A, Hallenbeck P: Increasing the metabolic capacity of Escherichia coli for hydrogen production through heterologous expression of the Ralstonia eutropha SH operon. Biotechnol Biofuels. 2013, 6: 122-10.1186/1754-6834-6-122.

    Article  Google Scholar 

  8. 8.

    Sarma SJ, Brar SK, Le Bihan Y, Buelna G, Soccol CR: Mitigation of the inhibitory effect of soap by magnesium salt treatment of crude glycerol - a novel approach for enhanced biohydrogen production from the biodiesel industry waste. Bioresour Technol. 2014, 151: 49-53. 10.1016/j.biortech.2013.10.042.

    Article  Google Scholar 

  9. 9.

    Li W-W, Yu H-Q, He Z: Towards sustainable wastewater treatment by using microbial fuel cells-centered technologies. Energy Environ Sci. 2014, 7: 911-924. 10.1039/c3ee43106a.

    Article  Google Scholar 

  10. 10.

    Motte J-C, Trably E, Escudie R, Hamelin J, Steyer J-P, Bernet N, Delgenes J-P, Dumas C: Total solids content: a key parameter of metabolic pathways in dry anaerobic digestion. Biotechnol Biofuels. 2013, 6: 164-10.1186/1754-6834-6-164.

    Article  Google Scholar 

  11. 11.

    Roy S, Kumar K, Ghosh S, Das D: Thermophilic biohydrogen production using pre-treated algal biomass as substrate. Biomass Bioenergy. 2014, 61: 157-166. 10.1016/j.biombioe.2013.12.006.

    Article  Google Scholar 

  12. 12.

    Lin X, Xia Y, Yan Q, Shen W, Zhao M: Acid tolerance response (ATR) of microbial communities during the enhanced biohydrogen process via cascade acid stress. Bioresour Technol. 2014, 155: 98-103. 10.1016/j.biortech.2013.12.097.

    Article  Google Scholar 

  13. 13.

    Wicher E, Seifert K, Zagrodnik R, Pietrzyk B, Laniecki M: Hydrogen gas production from distillery wastewater by dark fermentation. Int J Hydrog Energy. 2013, 38: 7767-7773. 10.1016/j.ijhydene.2013.04.008.

    Article  Google Scholar 

  14. 14.

    Fernández C, Carracedo B, Martínez EJ, G mez X, Morán A: Application of a packed bed reactor for the production of hydrogen from cheese whey permeate: effect of organic loading rate. J Environ Sci Health A. 2014, 49: 210-217. 10.1080/10934529.2013.838885.

    Article  Google Scholar 

  15. 15.

    Mahmoud N, Tawfik A, Ookawara S, Suzuki M: Biological hydrogen production from starch wastewater using a novel up-flow anaerobic staged reactor. BioResources. 2013, 8: 4951-4968.

    Google Scholar 

  16. 16.

    Deng L, Huang B, Zhao HY, Du JY, Gao L: Study on modified PVA-H3BO3 immobilization microorganism method for hydrogen production from wastewater. Adv Mater Res. 2013, 634: 280-285. 10.4028/

    Article  Google Scholar 

  17. 17.

    Sekoai PT, Gueguim Kana EB: Semi-pilot scale production of hydrogen from Organic Fraction of Solid Municipal Waste and electricity generation from process effluents. Biomass Bioenergy. 2014, 60: 156-163. 10.1016/j.biombioe.2013.11.008.

    Article  Google Scholar 

  18. 18.

    Qin Z, Qin Q, Yang Y: Continuous biohydrogen production with CSTR reactor under high organic loading rate condition. Adv Mater Res. 2014, 864: 225-228.

    Google Scholar 

  19. 19.

    McKinlay JB, Oda Y, R hl M, Posto AL, Sauer U, Harwood CS: Non-growing Rhodopseudomonas palustris increases the hydrogen gas yield from acetate by shifting from the glyoxylate shunt to the tricarboxylic acid cycle. J Biol Chem. 2014, 289: 1960-1970. 10.1074/jbc.M113.527515.

    Article  Google Scholar 

  20. 20.

    Cai G, Jin B, Monis P, Saint C: A genetic and metabolic approach to redirection of biochemical pathways of Clostridium butyricum for enhancing hydrogen production. Biotechnol Bioeng. 2013, 110: 338-342. 10.1002/bit.24596.

    Article  Google Scholar 

  21. 21.

    Cha M, Chung D, Elkins J, Guss A, Westpheling J: Metabolic engineering of Caldicellulosiruptor bescii yields increased hydrogen production from lignocellulosic biomass. Biotechnol Biofuels. 2013, 6: 85-10.1186/1754-6834-6-85.

    Article  Google Scholar 

  22. 22.

    Maru BT, Bielen AAM, Constanti M, Medina F, Kengen SWM: Glycerol fermentation to hydrogen by Thermotoga maritima: proposed pathway and bioenergetic considerations. Int J Hydrog Energy. 2013, 38: 5563-5572. 10.1016/j.ijhydene.2013.02.130.

    Article  Google Scholar 

  23. 23.

    Kumaraswamy GK, Guerra T, Qian X, Zhang S, Bryant DA, Dismukes GC: Reprogramming the glycolytic pathway for increased hydrogen production in cyanobacteria: metabolic engineering of NAD + -dependent GAPDH. Energy Environ Sci. 2013, 6: 3722-3731. 10.1039/c3ee42206b.

    Article  Google Scholar 

  24. 24.

    Varrone C, Van Nostrand JD, Liu W, Zhou B, Wang Z, Liu F, He Z, Wu L, Zhou J, Wang A: Metagenomic-based analysis of biofilm communities for electrohydrogenesis: from wastewater to hydrogen. Int J Hydrog Energy. 2014, 39: 4222-4233. 10.1016/j.ijhydene.2014.01.001.

    Article  Google Scholar 

  25. 25.

    Abreu A, Karakashev D, Angelidaki I, Sousa D, Alves M: Biohydrogen production from arabinose and glucose using extreme thermophilic anaerobic mixed cultures. Biotechnol Biofuels. 2012, 5: 6-10.1186/1754-6834-5-6.

    Article  Google Scholar 

  26. 26.

    Rafrafi Y, Trably E, Hamelin J, Latrille E, Meynial-Salles I, Benomar S, Giudici-Orticoni M-T, Steyer J-P: Sub-dominant bacteria as keystone species in microbial communities producing bio-hydrogen. Int J Hydrog Energy. 2013, 38: 4975-4985. 10.1016/j.ijhydene.2013.02.008.

    Article  Google Scholar 

  27. 27.

    Ning Y-Y, Wang S-F, Jin D-W, Harada H, Shi X-Y: Formation of hydrogen-producing granules and microbial community analysis in a UASB reactor. Renew Energy. 2013, 53: 12-17. 10.1016/j.renene.2012.10.051.

    Article  Google Scholar 

  28. 28.

    Gadhe A, Sonawane SS, Varma MN: Kinetic analysis of biohydrogen production from complex dairy wastewater under optimized condition. Int J Hydrog Energy. 2014, 39: 1306-1314. 10.1016/j.ijhydene.2013.11.022.

    Article  Google Scholar 

  29. 29.

    Masset J, Calusinska M, Hamilton C, Hiligsmann S, Joris B, Wilmotte A, Thonart P: Fermentative hydrogen production from glucose and starch using pure strains and artificial co-cultures of Clostridium spp. Biotechnol Biofuels. 2012, 5: 35-10.1186/1754-6834-5-35.

    Article  Google Scholar 

  30. 30.

    Gadow SI, Jiang H, Hojo T, Li Y-Y: Cellulosic hydrogen production and microbial community characterization in hyper-thermophilic continuous bioreactor. Int J Hydrog Energy. 2013, 38: 7259-7267. 10.1016/j.ijhydene.2013.04.041.

    Article  Google Scholar 

  31. 31.

    Park J-H, Lee S-H, Yoon J-J, Kim S-H, Park H-D: Predominance of cluster Clostridium in hydrogen fermentation of galactose seeded with various heat-treated anaerobic sludges. Bioresour Technol. 2014, 157: 98-106. 10.1016/j.biortech.2014.01.081.

    Article  Google Scholar 

  32. 32.

    Penteado ED, Lazaro CZ, Sakamoto IK, Zaiat M: Influence of seed sludge and pretreatment method on hydrogen production in packed-bed anaerobic reactors. Int J Hydrog Energy. 2013, 38: 6137-6145. 10.1016/j.ijhydene.2013.01.067.

    Article  Google Scholar 

  33. 33.

    El-Bery H, Tawfik A, Kumari S, Bux F: Effect of thermal pre-treatment on inoculum sludge to enhance bio-hydrogen production from alkali hydrolysed rice straw in a mesophilic anaerobic baffled reactor. Environ Technol. 2013, 34: 1965-1972. 10.1080/09593330.2013.824013.

    Article  Google Scholar 

  34. 34.

    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 Hydrog Energy. 2014, 39: 3160-3175. 10.1016/j.ijhydene.2013.12.057.

    Article  Google Scholar 

  35. 35.

    Guo L, Zhao J, She Z, Lu M, Zong Y: Statistical key factors optimization of conditions for hydrogen production from S-TE (Solubilization by thermophilic enzyme) waste sludge. Bioresour Technol. 2013, 137: 51-56. 10.1016/j.biortech.2013.03.102.

    Article  Google Scholar 

  36. 36.

    Gadhe A, Sonawane SS, Varma MN: Optimization of conditions for hydrogen production from complex dairy wastewater by anaerobic sludge using desirability function approach. Int J Hydrog Energy. 2013, 38: 6607-6617. 10.1016/j.ijhydene.2013.03.078.

    Article  Google Scholar 

  37. 37.

    Whiteman J, Kana EG: Comparative assessment of the artificial neural network and response surface modelling efficiencies for biohydrogen production on sugar cane molasses. BioEnergy Res. 2013, 7: 295-305. 10.1007/s12155-013-9375-7.

    Article  Google Scholar 

  38. 38.

    Xiao Y, Zhang X, Zhu M, Tan W: Effect of the culture media optimization, pH and temperature on the biohydrogen production and the hydrogenase activities by Klebsiella pneumoniae ECU-15. Bioresour Technol. 2013, 137: 9-17. 10.1016/j.biortech.2013.03.109.

    Article  Google Scholar 

  39. 39.

    Wu X, Lin H, Zhu J: Optimization of continuous hydrogen production from co-fermenting molasses with liquid swine manure in an anaerobic sequencing batch reactor. Bioresour Technol. 2013, 136: 351-359. 10.1016/j.biortech.2013.02.109.

    Article  Google Scholar 

  40. 40.

    Mullai P, Yogeswari M, Sridevi K: Optimisation and enhancement of biohydrogen production using nickel nanoparticles - a novel approach. Bioresour Technol. 2013, 141: 212-219. 10.1016/j.biortech.2013.03.082.

    Article  Google Scholar 

  41. 41.

    Pakshirajan K, Mal J: Biohydrogen production using native carbon monoxide converting anaerobic microbial consortium predominantly Petrobacter sp. Int J Hydrog Energy. 2013, 38: 16020-16028. 10.1016/j.ijhydene.2013.09.129.

    Article  Google Scholar 

  42. 42.

    Boonsayompoo O, Reungsang A: Thermophilic biohydrogen production from the enzymatic hydrolysate of cellulose fraction of sweet sorghum bagasse by Thermoanaerobacterium thermosaccharolyticum KKU19: optimization of media composition. Int J Hydrog Energy. 2013, 38: 15777-15786. 10.1016/j.ijhydene.2013.04.129.

    Article  Google Scholar 

  43. 43.

    Jollet V, Gissane C, Schlaf M: Optimization of the neutralization of Red Mud by pyrolysis bio-oil using a design of experiments approach. Energy Environ Sci. 2014, 7: 1125-1133. 10.1039/c3ee43585g.

    Article  Google Scholar 

  44. 44.

    Boboescu IZ, Gherman VD, Mirel I, Pap B, Tengölics R, Rákhely G, Kovács KL, Kondorosi É, Maróti G: Simultaneous biohydrogen production and wastewater treatment based on the selective enrichment of the fermentation ecosystem. Int J Hydrog Energy. 2014, 39: 1502-1510. 10.1016/j.ijhydene.2013.08.139.

    Article  Google Scholar 

  45. 45.

    Sharma PK, Capalash N, Kaur J: An improved method for single step purification of metagenomic DNA. Mol Biotechnol. 2007, 36: 61-63. 10.1007/s12033-007-0015-3.

    Article  Google Scholar 

  46. 46.

    Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A, Wilkening J, Edwards RA: The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008, 9: 386-10.1186/1471-2105-9-386.

    Article  Google Scholar 

Download references


This work was supported by the following international (EU) and domestic funding bodies: "SYMBIOTICS" ERC AdG EU Grant, "BIOSIM" PN-II-PT-PCCA-2011-3.1-1129 European Fund (Romania), and by PIAC_13-1-2013-0145 supported by the Hungarian Government and financed by the Research and Technology Innovation Fund.

Author information



Corresponding author

Correspondence to Gergely Maróti.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Boboescu, I.Z., Ilie, M., Gherman, V.D. et al. Revealing the factors influencing a fermentative biohydrogen production process using industrial wastewater as fermentation substrate. Biotechnol Biofuels 7, 139 (2014).

Download citation


  • Biohydrogen
  • Central composite experimental design
  • Microbial inocula
  • Metagenomics
  • Industrial wastewater