A thermodynamically-based model for predicting microbial growth and community composition coupled to system geochemistry: Application to uranium bioreduction

Multiscale Modeling of Uranium Bioreduction in Porous Media by One-Dimensional Biofilms

We formulate a multiscale mathematical model that describes the bioreduction of uranium in porous media. On the mesoscale we describe the bioreduction of uranium [VI] to uranium [IV] using a multispecies one-dimensional biofilm model with suspended bacteria and thermodynamic growth inhibition. We upscale the mesoscopic (colony scale) model to the macroscale (reactor scale) and investigate the behavior of substrate utilization and production, attachment and detachment processes, and thermodynamic effects not usually considered in biofilm growth models. Simulation results of the reactor model indicate that thermodynamic inhibition quantitatively alters the dynamics of the model and neglecting thermodynamic effects may over- or underestimate chemical concentrations in the system. Furthermore, we numerically investigate uncertainties related to the specific choice of attachment and detachment rate coefficients and find that while increasing the attachment rate coefficient or decreasing the detachment rate coefficient leads to thicker biofilms, performance of the reactor remains largely unaffected.

Anaerobic respiration pathways and response to increased substrate availability of Arctic wetland soils

The availability of labile carbon (C) compounds in Arctic wetland soils is expected to increase due to thawing permafrost and increased fermentation as a result of decomposition of organic matter with warming. How microbial communities respond to this change will affect the balance of CO₂ and CH₄ emitted during anaerobic organic matter decomposition, and ultimately the net radiative forcing of greenhouse gas emissions from these soils. While soil water content limits aerobic respiration, the factors controlling methanogenesis and anaerobic respiration are poorly defined in suboxic Arctic soils. We conducted incubation experiments on two tundra soils from field sites on the Seward Peninsula, Alaska, with contrasting pH and geochemistry to determine the pathways of anaerobic microbial respiration and changes with increasing substrate availability upon warming. In incubation of soils from the circumneutral Teller site, the ratio of CO₂ to CH₄ dropped from 10 to <2 after 60 days, indicating rapid depletion of alternative terminal electron acceptors (TEAs). Addition of acetate stimulated production of CO₂ and CH₄ in a nearly 1 : 1 ratio, consistent with methanogenesis, and the composition of the microbial community shifted to favor clades capable of utilizing the added acetate such as the Fe(iii)-reducing Geobacter and the methanogenic archaea Methanosarcina. In contrast, both CO₂ and CH₄ production declined with acetate addition during incubation of soils from the more acidic Council site, and fermentative microorganisms increased in abundance despite the high availability of fermentation products. These results demonstrate that the degree to which increasing substrate availability stimulates greenhouse gas production in tundra wetlands will vary widely depending on soil pH and geochemistry.

Thermodynamic Inhibition in Chemostat Models : With an Application to Bioreduction of Uranium

We formulate a mathematical model of bacterial populations in a chemostat setting that also accounts for thermodynamic growth inhibition as a consequence of chemical reactions. Using only elementary mathematical and chemical arguments, we carry this out for two systems: a simple toy model with a single species, a single substrate, and a single reaction product, and a more involved model that describes bioreduction of uranium[VI] into uranium[IV]. We find that in contrast to most traditional chemostat models, as a consequence of thermodynamic inhibition the equilibria concentrations of nutrient substrates might depend on their inflow concentration and not only on reaction parameters and the reactor's dilution rate. Simulation results of the uranium degradation model indicate that thermodynamic growth inhibition quantitatively alters the solution of the model. This suggests that neglecting thermodynamic inhibition effects in systems where they play a role might lead to wrong model predictions and under- or over-estimate the efficacy of the process under investigation.

Engineering and kinetic aspects of bacterial uranium reduction for the remediation of uranium contaminated environments

Biological reduction of soluble uranium from U(VI) to insoluble U(IV) coupled to the oxidation of an electron donor (hydrogen or organic compounds) is a potentially cost-efficient way to reduce the U concentrations in contaminated waters to below regulatory limits. A variety of microorganisms originating from both U contaminated and non-contaminated environments have demonstrated U(VI) reduction capacity under anaerobic conditions. Bioreduction of U(VI) is considered especially promising for in situ remediation, where the activity of indigenous microorganisms is stimulated by supplying a suitable electron donor to the subsurface to contain U contamination to a specific location in a sparingly soluble form. Less studied microbial biofilm-based bioreactors and bioelectrochemical systems have also shown potential for efficient U(VI) reduction to remove U from contaminated water streams. This review compares the advantages and challenges of U(VI)-reducing in situ remediation processes, bioreactors and bioelectrochemical systems. In addition, the current knowledge of U(VI) bioreduction mechanisms and factors affecting U(VI) reduction kinetics (e.g. pH, temperature, and the chemical composition of the contaminated water) are discussed, as both of these aspects are important in designing efficient remediation processes.

Targeted quantification of functional enzyme dynamics in environmental samples for microbially mediated biogeochemical processes

Microbial enzymes catalytically drive biogeochemical processes in environments. The dynamic linkage between functional enzymes and biogeochemical species transformation has, however, rarely been investigated for decades because of the challenges to directly quantify enzymes in environmental samples. The diversity of microorganisms, the low amount of available biomass and the complexity of chemical composition in environmental samples represent the main challenges. To address the diversity challenge, we first identify several signature peptides that are conserved in the targeted enzymes with the same functionality across many phylogenetically diverse microorganisms using metagenome-based protein sequence data. Quantification of the signature peptides then allows estimation of the targeted enzyme abundance. To achieve analyses of the requisite sensitivity for complex environmental samples with low available biomass, we adapted a recently developed ultrasensitive targeted quantification technology, termed high-pressure high-resolution separations with intelligent selection and multiplexing (PRISM) by improving peptide separation efficiency and method detection sensitivity. Nitrate reduction dynamics catalyzed by dissimilatory and assimilatory enzymes in a hyporheic zone sediment was used as an example to demonstrate the application of the enzyme quantification approach. Together with the measurements of biogeochemical species, the approach enables investigating the dynamic linkage between functional enzymes and biogeochemical processes.

An improved cellular automaton method to model multispecies biofilms

Biomass-spreading rules used in previous cellular automaton methods to simulate multispecies biofilm introduced extensive mixing between different biomass species or resulted in spatially discontinuous biomass concentration and distribution; this caused results based on the cellular automaton methods to deviate from experimental results and those from the more computationally intensive continuous method. To overcome the problems, we propose new biomass-spreading rules in this work: Excess biomass spreads by pushing a line of grid cells that are on the shortest path from the source grid cell to the destination grid cell, and the fractions of different biomass species in the grid cells on the path change due to the spreading. To evaluate the new rules, three two-dimensional simulation examples are used to compare the biomass distribution computed using the continuous method and three cellular automaton methods, one based on the new rules and the other two based on rules presented in two previous studies. The relationship between the biomass species is syntrophic in one example and competitive in the other two examples. Simulation results generated using the cellular automaton method based on the new rules agree much better with the continuous method than do results using the other two cellular automaton methods. The new biomass-spreading rules are no more complex to implement than the existing rules.

U(VI) bioreduction with emulsified vegetable oil as the electron donor--microcosm tests and model development

We conducted microcosm tests and biogeochemical modeling to study U(VI) reduction in contaminated sediments amended with emulsified vegetable oil (EVO). Indigenous microorganisms in the sediments degraded EVO and stimulated Fe(III), U(VI), and sulfate reduction, and methanogenesis. Acetate concentration peaked in 100-120 days in the EVO microcosms versus 10-20 days in the oleate microcosms, suggesting that triglyceride hydrolysis was a rate-limiting step in EVO degradation and subsequent reactions. Acetate persisted 50 days longer in oleate- and EVO- than in ethanol-amended microcosms, indicating that acetate-utilizing methanogenesis was slower in the oleate and EVO than ethanol microcosms. We developed a comprehensive biogeochemical model to couple EVO hydrolysis, production, and oxidation of long-chain fatty acids (LCFA), glycerol, acetate, and hydrogen, reduction of Fe(III), U(VI) and sulfate, and methanogenesis with growth and decay of multiple functional microbial groups. By estimating EVO, LCFA, and glycerol degradation rate coefficients, and introducing a 100 day lag time for acetoclastic methanogenesis for oleate and EVO microcosms, the model approximately matched observed sulfate, U(VI), and acetate concentrations. Our results confirmed that EVO could stimulate U(VI) bioreduction in sediments and the slow EVO hydrolysis and acetate-utilizing methanogens growth could contribute to longer term bioreduction than simple substrates (e.g., ethanol, acetate, etc.) in the subsurface.

Uranium and other heavy metal resistance and accumulation in bacteria isolated from uranium mine wastes

Ten bacterial strains isolated from uranium mine wastes were characterized in terms of their uranium and other metal resistance and accumulation. 16S rRNA gene sequence analysis identified the strains as members of genera Bacillus, Serratia, and Arthrobacter. Strains were able to utilize various carbon sources, particularly aromatic hydrocarbons, grow at broad pH and temperature ranges and produce non specific acid phosphatase relevant for metal phosphate precipitation in contaminated environment. The isolates exhibited high uranium and other heavy metals (Ni, Co, Cu and Cd) resistance and accumulation capacities. Particularly, Arthrobacter sp. J001 and Bacillus sp. J003 were superior in terms of U resistance at low pH (pH 4.0) along with metals and actinides (U and Th) removal with maximum cell loading of 1088 μmol U, 1293 μmol Th, 425 μmol Cu, 305 μmol Cd, 377 μmol Zn, 250 μmol Ni g(-1) cell dry wt. Genes encoding P(1B)-type ATPases (Cu-CPx and Zn-CPx) and ABC transporters (nik) as catalytic tools for maintaining cellular metal homeostasis were detected within several Bacillus spp., with possible incidence of horizontal gene transfer for the later gene showing phylogenetic lineage to α Proteobacteria members. The study provides evidence on intrinsic abilities of indigenous bacteria from U-mine suitable for survival and cleaning up of contaminated mine sites.

Microbial physiology-based model of ethanol metabolism in subsurface sediments

A biogeochemical reaction model was developed based on microbial physiology to simulate ethanol metabolism and its influence on the chemistry of anoxic subsurface environments. The model accounts for potential microbial metabolisms that degrade ethanol, including those that oxidize ethanol directly or syntrophically by reducing different electron acceptors. Out of the potential metabolisms, those that are active in the environment can be inferred by fitting the model to experimental observations. This approach was applied to a batch sediment slurry experiment that examined ethanol metabolism in uranium-contaminated aquifer sediments from Area 2 at the U.S. Department of Energy Field Research Center in Oak Ridge, TN. According to the simulation results, complete ethanol oxidation by denitrification, incomplete ethanol oxidation by ferric iron reduction, ethanol fermentation to acetate and H(2), hydrogenotrophic sulfate reduction, and acetoclastic methanogenesis: all contributed significantly to the degradation of ethanol in the aquifer sediments. The assemblage of the active metabolisms provides a frame work to explore how ethanol amendment impacts the chemistry of the environment, including the occurrence and levels of uranium. The results can also be applied to explore how diverse microbial metabolisms impact the progress and efficacy of bioremediation strategies.

Thermodynamics of microbial growth coupled to metabolism of glucose, ethanol, short-chain organic acids, and hydrogen

A literature compilation demonstrated a linear relationship between microbial growth yield and the free energy of aerobic and anaerobic (respiratory and/or fermentative) metabolism of glucose, ethanol, formate, acetate, lactate, propionate, butyrate, and H(2). This relationship provides a means to estimate growth yields for modeling microbial redox metabolism in soil and sedimentary environments.