Deterministic influences exceed dispersal effects on hydrologically-connected microbiomes.

Subsurface groundwater-surface water mixing zones (hyporheic zones) have enhanced biogeochemical activity, but assembly processes governing subsurface microbiomes remain a critical uncertainty in understanding hyporheic biogeochemistry. To address this obstacle, we investigated (a) biogeographical patterns in attached and waterborne microbiomes across three hydrologically-connected, physicochemically-distinct zones (inland hyporheic, nearshore hyporheic and river); (b) assembly processes that generated these patterns; (c) groups of organisms that corresponded to deterministic changes in the environment; and (d) correlations between these groups and hyporheic metabolism. All microbiomes remained dissimilar through time, but consistent presence of similar taxa suggested dispersal and/or common selective pressures among zones. Further, we demonstrated a pronounced impact of deterministic assembly in all microbiomes as well as seasonal shifts from heterotrophic to autotrophic microorganisms associated with increases in groundwater discharge. The abundance of one statistical cluster of organisms increased with active biomass and respiration, revealing organisms that may strongly influence hyporheic biogeochemistry. Based on our results, we propose a conceptualization of hyporheic zone metabolism in which increased organic carbon concentrations during surface water intrusion support heterotrophy, which succumbs to autotrophy under groundwater discharge. These results provide new opportunities to enhance microbially-explicit ecosystem models describing hyporheic zone biogeochemistry and its influence over riverine ecosystem function.

The 'porin-cytochrome' model for microbe-to-mineral electron transfer.

Many species of bacteria can couple anaerobic growth to the respiratory reduction of insoluble minerals containing Fe(III) or Mn(III/IV). It has been suggested that in Shewanella species electrons cross the outer membrane to extracellular substrates via 'porin-cytochrome' electron transport modules. The molecular structure of an outer-membrane extracellular-facing deca-haem terminus for such a module has recently been resolved. It is debated how, once outside the cells, electrons are transferred from outer-membrane cytochromes to insoluble electron sinks. This may occur directly or by assemblies of cytochromes, perhaps functioning as 'nanowires', or via electron shuttles. Here we review recent work in this field and explore whether it allows for unification of the electron transport mechanisms supporting extracellular mineral respiration in Shewanella that may extend into other genera of Gram-negative bacteria.

Genome-scale modeling of light-driven reductant partitioning and carbon fluxes in diazotrophic unicellular cyanobacterium Cyanothece sp. ATCC 51142.

Genome-scale metabolic models have proven useful for answering fundamental questions about metabolic capabilities of a variety of microorganisms, as well as informing their metabolic engineering. However, only a few models are available for oxygenic photosynthetic microorganisms, particularly in cyanobacteria in which photosynthetic and respiratory electron transport chains (ETC) share components. We addressed the complexity of cyanobacterial ETC by developing a genome-scale model for the diazotrophic cyanobacterium, Cyanothece sp. ATCC 51142. The resulting metabolic reconstruction, iCce806, consists of 806 genes associated with 667 metabolic reactions and includes a detailed representation of the ETC and a biomass equation based on experimental measurements. Both computational and experimental approaches were used to investigate light-driven metabolism in Cyanothece sp. ATCC 51142, with a particular focus on reductant production and partitioning within the ETC. The simulation results suggest that growth and metabolic flux distributions are substantially impacted by the relative amounts of light going into the individual photosystems. When growth is limited by the flux through photosystem I, terminal respiratory oxidases are predicted to be an important mechanism for removing excess reductant. Similarly, under photosystem II flux limitation, excess electron carriers must be removed via cyclic electron transport. Furthermore, in silico calculations were in good quantitative agreement with the measured growth rates whereas predictions of reaction usage were qualitatively consistent with protein and mRNA expression data, which we used to further improve the resolution of intracellular flux values.

Biofilm shows spatially stratified metabolic responses to contaminant exposure.

Biofilms are core to a range of biological processes, including the bioremediation of environmental contaminants. Within a biofilm population, cells with diverse genotypes and phenotypes coexist, suggesting that distinct metabolic pathways may be expressed based on the local environmental conditions in a biofilm. However, metabolic responses to local environmental conditions in a metabolically active biofilm interacting with environmental contaminants have never been quantitatively elucidated. In this study, we monitored the spatiotemporal metabolic responses of metabolically active Shewanella oneidensis MR-1 biofilms to U(VI) (uranyl, UO(2)(2+)) and Cr(VI) (chromate, CrO(4) (2-)) using non-invasive nuclear magnetic resonance imaging (MRI) and spectroscopy (MRS) approaches to obtain insights into adaptation in biofilms during biofilm-contaminant interactions. While overall biomass distribution was not significantly altered upon exposure to U(VI) or Cr(VI), MRI and spatial mapping of the diffusion revealed localized changes in the water diffusion coefficients in the biofilms, suggesting significant contaminant-induced changes in structural or hydrodynamic properties during bioremediation. Finally, we quantitatively demonstrated that the metabolic responses of biofilms to contaminant exposure are spatially stratified, implying that adaptation in biofilms is custom-developed based on local microenvironments.

Fe(III) reduction and U(VI) immobilization by Paenibacillus sp. strain 300A, isolated from Hanford 300A subsurface sediments.

A facultative iron-reducing [Fe(III)-reducing] Paenibacillus sp. strain was isolated from Hanford 300A subsurface sediment biofilms that was capable of reducing soluble Fe(III) complexes [Fe(III)-nitrilotriacetic acid and Fe(III)-citrate] but unable to reduce poorly crystalline ferrihydrite (Fh). However, Paenibacillus sp. 300A was capable of reducing Fh in the presence of low concentrations (2 µM) of either of the electron transfer mediators (ETMs) flavin mononucleotide (FMN) or anthraquinone-2,6-disulfonate (AQDS). Maximum initial Fh reduction rates were observed at catalytic concentrations (<10 µM) of either FMN or AQDS. Higher FMN concentrations inhibited Fh reduction, while increased AQDS concentrations did not. We also found that Paenibacillus sp. 300A could reduce Fh in the presence of natural ETMs from Hanford 300A subsurface sediments. In the absence of ETMs, Paenibacillus sp. 300A was capable of immobilizing U(VI) through both reduction and adsorption. The relative contributions of adsorption and microbial reduction to U(VI) removal from the aqueous phase were ∼7:3 in PIPES [piperazine-N,N'-bis(2-ethanesulfonic acid)] and ∼1:4 in bicarbonate buffer. Our study demonstrated that Paenibacillus sp. 300A catalyzes Fe(III) reduction and U(VI) immobilization and that these reactions benefit from externally added or naturally existing ETMs in 300A subsurface sediments.

Redox reactions of reduced flavin mononucleotide (FMN), riboflavin (RBF), and anthraquinone-2,6-disulfonate (AQDS) with ferrihydrite and lepidocrocite.

Flavins are secreted by the dissimilatory iron-reducing bacterium Shewanella and can function as endogenous electron transfer mediators. To assess the potential importance of flavins in Fe(III) bioreduction, we investigated the redox reaction kinetics of reduced flavin mononucleotide (i.e., FMNH(2)) and reduced riboflavin (i.e., RBFH(2)) with ferrihydrite and lepidocrocite. The organic reductants rapidly reduced and dissolved ferrihydrite and lepidocrocite in the pH range 4-8. The rate constant k for 2-line ferrihydrite reductive dissolution by FMNH(2) was 87.5 ± 3.5 M(-1)·s(-1) at pH 7.0 in batch reactors, and k was similar for RBFH(2). For lepidocrocite, k was 500 ± 61 M(-1)·s(-1) for FMNH(2) and 236 ± 22 M(-1)·s(-1) for RBFH(2). The surface area normalized initial reaction rates (r(a)) were between 0.08 and 77 µmol·m(-2)·s(-1) for various conditions in stopped-flow experiments. Initial rates (r(o)) were first-order with respect to iron(III) oxide concentration, and r(a) increased with decreasing pH. Poorly crystalline 2-line ferrihydrite yielded the highest r(a), followed by more crystalline 6-line ferrihydrite and crystalline lepidocrocite. Compared to a previous whole-cell study with Shewanella oneidensis strain MR-1, our findings suggest that the reduction of electron transfer mediators by the Mtr (i.e., metal-reducing) pathway coupled to lactate oxidation is rate limiting, rather than heterogeneous electron transfer to the iron(III) oxide.

Pore-scale characterization of biogeochemical controls on iron and uranium speciation under flow conditions.

Etched silicon microfluidic pore network models (micromodels) with controlled chemical and redox gradients, mineralogy, and microbiology under continuous flow conditions are used for the incremental development of complex microenvironments that simulate subsurface conditions. We demonstrate the colonization of micromodel pore spaces by an anaerobic Fe(III)-reducing bacterial species (Geobacter sulfurreducens) and the enzymatic reduction of a bioavailable Fe(III) phase within this environment. Using both X-ray microprobe and X-ray absorption spectroscopy, we investigate the combined effects of the precipitated Fe(III) phases and the microbial population on uranium biogeochemistry under flow conditions. Precipitated Fe(III) phases within the micromodel were most effectively reduced in the presence of an electron shuttle (AQDS), and Fe(II) ions adsorbed onto the precipitated mineral surface without inducing any structural change. In the absence of Fe(III), U(VI) was effectively reduced by the microbial population to insoluble U(IV), which was precipitated in discrete regions associated with biomass. In the presence of Fe(III) phases, however, both U(IV) and U(VI) could be detected associated with biomass, suggesting reoxidation of U(IV) by localized Fe(III) phases. These results demonstrate the importance of the spatial localization of biomass and redox active metals, and illustrate the key effects of pore-scale processes on contaminant fate and reactive transport.

Genomic reconstruction of Shewanella oneidensis MR-1 metabolism reveals a previously uncharacterized machinery for lactate utilization.

The ability to use lactate as a sole source of carbon and energy is one of the key metabolic signatures of Shewanellae, a diverse group of dissimilatory metal-reducing bacteria commonly found in aquatic and sedimentary environments. Nonetheless, homology searches failed to recognize orthologs of previously described bacterial d- or l-lactate oxidizing enzymes (Escherichia coli genes dld and lldD) in any of the 13 analyzed genomes of Shewanella spp. By using comparative genomic techniques, we identified a conserved chromosomal gene cluster in Shewanella oneidensis MR-1 (locus tag: SO_1522-SO_1518) containing lactate permease and candidate genes for both d- and l-lactate dehydrogenase enzymes. The predicted d-LDH gene (dld-II, SO_1521) is a distant homolog of FAD-dependent lactate dehydrogenase from yeast, whereas the predicted l-LDH is encoded by 3 genes with previously unknown functions (lldEGF, SO_1520-SO_1518). Through a combination of genetic and biochemical techniques, we experimentally confirmed the predicted physiological role of these novel genes in S. oneidensis MR-1 and carried out successful functional validation studies in Escherichia coli and Bacillus subtilis. We conclusively showed that dld-II and lldEFG encode fully functional d-and l-LDH enzymes, which catalyze the oxidation of the respective lactate stereoisomers to pyruvate. Notably, the S. oneidensis MR-1 LldEFG enzyme is a previously uncharacterized example of a multisubunit lactate oxidase. Comparative analysis of >400 bacterial species revealed the presence of LldEFG and Dld-II in a broad range of diverse species accentuating the potential importance of these previously unknown proteins in microbial metabolism.

Characterization of an electron conduit between bacteria and the extracellular environment.

A number of species of Gram-negative bacteria can use insoluble minerals of Fe(III) and Mn(IV) as extracellular respiratory electron acceptors. In some species of Shewanella, deca-heme electron transfer proteins lie at the extracellular face of the outer membrane (OM), where they can interact with insoluble substrates. To reduce extracellular substrates, these redox proteins must be charged by the inner membrane/periplasmic electron transfer system. Here, we present a spectro-potentiometric characterization of a trans-OM icosa-heme complex, MtrCAB, and demonstrate its capacity to move electrons across a lipid bilayer after incorporation into proteoliposomes. We also show that a stable MtrAB subcomplex can assemble in the absence of MtrC; an MtrBC subcomplex is not assembled in the absence of MtrA; and MtrA is only associated to the membrane in cells when MtrB is present. We propose a model for the modular organization of the MtrCAB complex in which MtrC is an extracellular element that mediates electron transfer to extracellular substrates and MtrB is a trans-OM spanning beta-barrel protein that serves as a sheath, within which MtrA and MtrC exchange electrons. We have identified the MtrAB module in a range of bacterial phyla, suggesting that it is widely used in electron exchange with the extracellular environment.

Extracellular polymeric substances from Shewanella sp. HRCR-1 biofilms: characterization by infrared spectroscopy and proteomics.

The composition of extracellular polymeric substances (EPS) from Shewanella sp. HRCR-1 biofilms was investigated using infrared spectroscopy and proteomics to provide insight into potential ecophysiological functions and redox activity of the EPS. Both bound and loosely associated EPS were extracted from Shewanella sp. HRCR-1 biofilms prepared using a hollow-fibre membrane biofilm reactor. Fourier transform infrared spectra revealed the presence of proteins, polysaccharides, nucleic acids, membrane lipids and fatty acids in the EPS fractions. Using a global proteomic approach, a total of 58 extracellular and outer membrane proteins were identified in the EPS. These included homologues of multiple Shewanella oneidensis MR-1 proteins that potentially contribute to key physiological biofilm processes, such as biofilm-promoting protein BpfA, surface-associated serine protease, nucleotidases (CpdB and UshA), an extracellular lipase, and oligopeptidases (PtrB and a M13 family oligopeptidase lipoprotein). In addition, 20 redox proteins were found in extracted EPS. Among the detected redox proteins were the homologues of two S. oneidensis MR-1 c-type cytochromes, MtrC and OmcA, which have been implicated in extracellular electron transfer. Given their detection in the EPS of Shewanella sp. HRCR-1 biofilms, c-type cytochromes may contribute to the possible redox activity of the biofilm matrix and play important roles in extracellular electron transfer reactions.

Pyruvate and lactate metabolism by Shewanella oneidensis MR-1 under fermentation, oxygen limitation, and fumarate respiration conditions.

Shewanella oneidensis MR-1 is a facultative anaerobe that derives energy by coupling organic matter oxidation to the reduction of a wide range of electron acceptors. Here, we quantitatively assessed the lactate and pyruvate metabolism of MR-1 under three distinct conditions: electron acceptor-limited growth on lactate with O(2), lactate with fumarate, and pyruvate fermentation. The latter does not support growth but provides energy for cell survival. Using physiological and genetic approaches combined with flux balance analysis, we showed that the proportion of ATP produced by substrate-level phosphorylation varied from 33% to 72.5% of that needed for growth depending on the electron acceptor nature and availability. While being indispensable for growth, the respiration of fumarate does not contribute significantly to ATP generation and likely serves to remove formate, a product of pyruvate formate-lyase-catalyzed pyruvate disproportionation. Under both tested respiratory conditions, S. oneidensis MR-1 carried out incomplete substrate oxidation, whereby the tricarboxylic acid (TCA) cycle did not contribute significantly. Pyruvate dehydrogenase was not involved in lactate metabolism under conditions of O(2) limitation but was required for anaerobic growth, likely by supplying reducing equivalents for biosynthesis. The results suggest that pyruvate fermentation by S. oneidensis MR-1 cells represents a combination of substrate-level phosphorylation and respiration, where pyruvate serves as an electron donor and an electron acceptor. Pyruvate reduction to lactate at the expense of formate oxidation is catalyzed by a recently described new type of oxidative NAD(P)H-independent d-lactate dehydrogenase (Dld-II). The results further indicate that pyruvate reduction coupled to formate oxidation may be accompanied by the generation of proton motive force.

Constraint-based model of Shewanella oneidensis MR-1 metabolism: a tool for data analysis and hypothesis generation.

Shewanellae are gram-negative facultatively anaerobic metal-reducing bacteria commonly found in chemically (i.e., redox) stratified environments. Occupying such niches requires the ability to rapidly acclimate to changes in electron donor/acceptor type and availability; hence, the ability to compete and thrive in such environments must ultimately be reflected in the organization and utilization of electron transfer networks, as well as central and peripheral carbon metabolism. To understand how Shewanella oneidensis MR-1 utilizes its resources, the metabolic network was reconstructed. The resulting network consists of 774 reactions, 783 genes, and 634 unique metabolites and contains biosynthesis pathways for all cell constituents. Using constraint-based modeling, we investigated aerobic growth of S. oneidensis MR-1 on numerous carbon sources. To achieve this, we (i) used experimental data to formulate a biomass equation and estimate cellular ATP requirements, (ii) developed an approach to identify cycles (such as futile cycles and circulations), (iii) classified how reaction usage affects cellular growth, (iv) predicted cellular biomass yields on different carbon sources and compared model predictions to experimental measurements, and (v) used experimental results to refine metabolic fluxes for growth on lactate. The results revealed that aerobic lactate-grown cells of S. oneidensis MR-1 used less efficient enzymes to couple electron transport to proton motive force generation, and possibly operated at least one futile cycle involving malic enzymes. Several examples are provided whereby model predictions were validated by experimental data, in particular the role of serine hydroxymethyltransferase and glycine cleavage system in the metabolism of one-carbon units, and growth on different sources of carbon and energy. This work illustrates how integration of computational and experimental efforts facilitates the understanding of microbial metabolism at a systems level.

Contribution of extracellular polymeric substances from Shewanella sp. HRCR-1 biofilms to U(VI) immobilization.

The goal of this study was to quantify the contribution of extracellular polymeric substances (EPS) to U(VI) immobilization by Shewanella sp. HRCR-1. Through comparison of U(VI) immobilization using cells with bound EPS (bEPS) and cells with minimal EPS, we show that (i) bEPS from Shewanella sp. HRCR-1 biofilms contribute significantly to U(VI) immobilization, especially at low initial U(VI) concentrations, through both sorption and reduction; (ii) bEPS can be considered a functional extension of the cells for U(VI) immobilization and they likely play more important roles at lower initial U(VI) concentrations; and (iii) the U(VI) reduction efficiency is dependent upon the initial U(VI) concentration and decreases at lower concentrations. To quantify the relative contributions of sorption and reduction to U(VI) immobilization by EPS fractions, we isolated loosely associated EPS (laEPS) and bEPS from Shewanella sp. HRCR-1 biofilms grown in a hollow fiber membrane biofilm reactor and tested their reactivity with U(VI). We found that, when reduced, the isolated cell-free EPS fractions could reduce U(VI). Polysaccharides in the EPS likely contributed to U(VI) sorption and dominated the reactivity of laEPS, while redox active components (e.g., outer membrane c-type cytochromes), especially in bEPS, possibly facilitated U(VI) reduction.

Conserved synteny at the protein family level reveals genes underlying Shewanella species' cold tolerance and predicts their novel phenotypes.

Bacteria of the genus Shewanella can thrive in different environments and demonstrate significant variability in their metabolic and ecophysiological capabilities including cold and salt tolerance. Genomic characteristics underlying this variability across species are largely unknown. In this study, we address the problem by a comparison of the physiological, metabolic, and genomic characteristics of 19 sequenced Shewanella species. We have employed two novel approaches based on association of a phenotypic trait with the number of the trait-specific protein families (Pfam domains) and on the conservation of synteny (order in the genome) of the trait-related genes. Our first approach is top-down and involves experimental evaluation and quantification of the species' cold tolerance followed by identification of the correlated Pfam domains and genes with a conserved synteny. The second, a bottom-up approach, predicts novel phenotypes of the species by calculating profiles of each Pfam domain among their genomes and following pair-wise correlation of the profiles and their network clustering. Using the first approach, we find a link between cold and salt tolerance of the species and the presence in the genome of a Na(+)/H(+) antiporter gene cluster. Other cold-tolerance-related genes include peptidases, chemotaxis sensory transducer proteins, a cysteine exporter, and helicases. Using the bottom-up approach, we found several novel phenotypes in the newly sequenced Shewanella species, including degradation of aromatic compounds by an aerobic hybrid pathway in Shewanella woodyi, degradation of ethanolamine by Shewanella benthica, and propanediol degradation by Shewanella putrefaciens CN32 and Shewanella sp. W3-18-1.

In situ effective diffusion coefficient profiles in live biofilms using pulsed-field gradient nuclear magnetic resonance.

Diffusive mass transfer in biofilms is characterized by the effective diffusion coefficient. It is well documented that the effective diffusion coefficient can vary by location in a biofilm. The current literature is dominated by effective diffusion coefficient measurements for distinct cell clusters and stratified biofilms showing this spatial variation. Regardless of whether distinct cell clusters or surface-averaging methods are used, position-dependent measurements of the effective diffusion coefficient are currently: (1) invasive to the biofilm, (2) performed under unnatural conditions, (3) lethal to cells, and/or (4) spatially restricted to only certain regions of the biofilm. Invasive measurements can lead to inaccurate results and prohibit further (time-dependent) measurements which are important for the mathematical modeling of biofilms. In this study our goals were to: (1) measure the effective diffusion coefficient for water in live biofilms, (2) monitor how the effective diffusion coefficient changes over time under growth conditions, and (3) correlate the effective diffusion coefficient with depth in the biofilm. We measured in situ two-dimensional effective diffusion coefficient maps within Shewanella oneidensis MR-1 biofilms using pulsed-field gradient nuclear magnetic resonance methods, and used them to calculate surface-averaged relative effective diffusion coefficient (D(rs)) profiles. We found that (1) D(rs) decreased from the top of the biofilm to the bottom, (2) D(rs) profiles differed for biofilms of different ages, (3) D(rs) profiles changed over time and generally decreased with time, (4) all the biofilms showed very similar D(rs) profiles near the top of the biofilm, and (5) the D(rs) profile near the bottom of the biofilm was different for each biofilm. Practically, our results demonstrate that advanced biofilm models should use a variable effective diffusivity which changes with time and location in the biofilm.

Prediction of Neighbor-Dependent Microbial Interactions From Limited Population Data.

Modulation of interspecies interactions by the presence of neighbor species is a key ecological factor that governs dynamics and function of microbial communities, yet the development of theoretical frameworks explicit for understanding context-dependent interactions are still nascent. In a recent study, we proposed a novel rule-based inference method termed the Minimal Interspecies Interaction Adjustment (MIIA) that predicts the reorganization of interaction networks in response to the addition of new species such that the modulation in interaction coefficients caused by additional members is minimal. While the theoretical basis of MIIA was established through the previous work by assuming the full availability of species abundance data in axenic, binary, and complex communities, its extension to actual microbial ecology can be highly constrained in cases that species have not been cultured axenically (e.g., due to their inability to grow in the absence of specific partnerships) because binary interaction coefficients - basic parameters required for implementing the MIIA - are inestimable without axenic and binary population data. Thus, here we present an alternative formulation based on the following two central ideas. First, in the case where only data from axenic cultures are unavailable, we remove axenic populations from governing equations through appropriate scaling. This allows us to predict neighbor-dependent interactions in a relative sense (i.e., fractional change of interactions between with versus without neighbors). Second, in the case where both axenic and binary populations are missing, we parameterize binary interaction coefficients to determine their values through a sensitivity analysis. Through the case study of two microbial communities with distinct characteristics and complexity (i.e., a three-member community where all members can grow independently, and a four-member community that contains member species whose growth is dependent on other species), we demonstrated that despite data limitation, the proposed new formulation was able to successfully predict interspecies interactions that are consistent with experimentally derived results. Therefore, this technical advancement enhances our ability to predict context-dependent interspecies interactions in a broad range of microbial systems without being limited to specific growth conditions as a pre-requisite.

Minimal Interspecies Interaction Adjustment (MIIA): Inference of Neighbor-Dependent Interactions in Microbial Communities.

An intriguing aspect in microbial communities is that pairwise interactions can be influenced by neighboring species. This creates context dependencies for microbial interactions that are based on the functional composition of the community. Context dependent interactions are ecologically important and clearly present in nature, yet firmly established theoretical methods are lacking from many modern computational investigations. Here, we propose a novel network inference method that enables predictions for interspecies interactions affected by shifts in community composition and species populations. Our approach first identifies interspecies interactions in binary communities, which is subsequently used as a basis to infer modulation in more complex multi-species communities based on the assumption that microbes minimize adjustments of pairwise interactions in response to neighbor species. We termed this rule-based inference minimal interspecies interaction adjustment (MIIA). Our critical assessment of MIIA has produced reliable predictions of shifting interspecies interactions that are dependent on the functional role of neighbor organisms. We also show how MIIA has been applied to a microbial community composed of competing soil bacteria to elucidate a new finding that - in many cases - adding fewer competitors could impose more significant impact on binary interactions. The ability to predict membership-dependent community behavior is expected to help deepen our understanding of how microbiomes are organized in nature and how they may be designed and/or controlled in the future.

Utilization of DNA as a sole source of phosphorus, carbon, and energy by Shewanella spp.: ecological and physiological implications for dissimilatory metal reduction.

The solubility of orthophosphate (PO4(3-)) in iron-rich sediments can be exceedingly low, limiting the bioavailability of this essential nutrient to microbial populations that catalyze critical biogeochemical reactions. Here we demonstrate that dissolved extracellular DNA can serve as a sole source of phosphorus, as well as carbon and energy, for metal-reducing bacteria of the genus Shewanella. Shewanella oneidensis MR-1, Shewanella putrefaciens CN32, and Shewanella sp. strain W3-18-1 all grew with DNA but displayed different growth rates. W3-18-1 exhibited the highest growth rate with DNA. While strain W3-18-1 displayed Ca2+-independent DNA utilization, both CN32 and MR-1 required millimolar concentrations of Ca2+ for growth with DNA. For S. oneidensis MR-1, the utilization of DNA as a sole source of phosphorus is linked to the activities of extracellular phosphatase(s) and a Ca2+-dependent nuclease(s), which are regulated by phosphorus availability. Mass spectrometry analysis of the extracellular proteome of MR-1 identified one putative endonuclease (SO1844), a predicted UshA (bifunctional UDP-sugar hydrolase/5' nucleotidase), a predicted PhoX (calcium-activated alkaline phosphatase), and a predicted CpdB (bifunctional 2',3' cyclic nucleotide 2' phosphodiesterase/3' nucleotidase), all of which could play important roles in the extracellular degradation of DNA under phosphorus-limiting conditions. Overall, the results of this study suggest that the ability to use exogenous DNA as the sole source of phosphorus is widespread among the shewanellae, and perhaps among all prokaryotes, and may be especially important for nutrient cycling in metal-reducing environments.

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process.

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community's traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators.

Indirect Interspecies Regulation: Transcriptional and Physiological Responses of a Cyanobacterium to Heterotrophic Partnership.

The mechanisms by which microbes interact in communities remain poorly understood. Here, we interrogated specific interactions between photoautotrophic and heterotrophic members of a model consortium to infer mechanisms that mediate metabolic coupling and acclimation to partnership. This binary consortium was composed of a cyanobacterium, Thermosynechococcus elongatus BP-1, which supported growth of an obligate aerobic heterotroph, Meiothermus ruber strain A, by providing organic carbon, O2, and reduced nitrogen. Species-resolved transcriptomic analyses were used in combination with growth and photosynthesis kinetics to infer interactions and the environmental context under which they occur. We found that the efficiency of biomass production and resistance to stress induced by high levels of dissolved O2 increased, beyond axenic performance, as a result of heterotrophic partnership. Coordinated transcriptional responses transcending both species were observed and used to infer specific interactions resulting from the synthesis and exchange of resources. The cyanobacterium responded to heterotrophic partnership by altering expression of core genes involved with photosynthesis, carbon uptake/fixation, vitamin synthesis, and scavenging of reactive oxygen species (ROS). IMPORTANCE This study elucidates how a cyanobacterial primary producer acclimates to heterotrophic partnership by modulating the expression levels of key metabolic genes. Heterotrophic bacteria can indirectly regulate the physiology of the photoautotrophic primary producers, resulting in physiological changes identified here, such as increased intracellular ROS. Some of the interactions inferred from this model system represent putative principles of metabolic coupling in phototrophic-heterotrophic partnerships.