Relationship between Microorganisms Inhabiting Alkaline Siliceous Hot Spring Mat Communities and Overflowing Water.

The compositions of Octopus Spring and Mushroom Spring (Yellowstone National Park, Wyoming, USA) microbial mats have been thoroughly studied, but the compositions of the effluent waters that flow above the mats have not. In this study, cells in the mats and overflowing waters of both springs were investigated at multiple sites where Synechococcus spp. are the dominant cyanobacteria (ca. 72°C to ca. 50°C), and on several dates. In addition to microscopic analyses of stained and autofluorescent cells, 16S rRNA gene sequencing was used to characterize the major taxa present and a protein-encoding gene (psaA) was sequenced and analyzed by ecotype simulation to predict species of Synechococcus The mats of both springs were similar in terms of the downstream distribution of predominant taxa detected previously. However, waters above these mats were predominated by taxa that reside in upstream mats or communities above the upper-temperature limit of the mat. A disturbance/recolonization study was performed at a site normally predominated by Synechococcus species adapted to low temperatures. After removing indigenous Synechococcus cells, Synechococcus species adapted to higher temperatures, which were predominant in the water overflowing this site, colonized the newly forming mat. Differences in recolonization under reduced and UV-screened irradiance suggested that, in addition to physical transport, environmental conditions likely select for species that are better adapted to these different conditions and can influence mat recovery. A transport model was developed and used to predict that, in Mushroom Spring, erosion predominates in the narrower and deeper upstream effluents and deposition predominates over erosion in wider and shallower downstream effluents.IMPORTANCE In flowing aquatic systems, cell erosion and deposition are important to the dispersal of cells from one location to another. Very little is known about microbial dispersal and the physical processes that underlie it. This study demonstrates its importance to colonization of downstream surfaces and especially to the recolonization and functioning of disturbed sites. Ecological systems in flowing environments are often, roughly speaking, pseudosteady, in that nutrients enter the system and by-products leave at relatively steady rates. Over time, material inputs and outputs must balance. Measurements of input fluxes (e.g., growth rates and proxies, such as photosynthesis rates) are frequent. However, erosion and deposition of cells are seldom measured and ecological significance is sometimes neglected. The importance of these parameters is immediately evident in any attempt to construct a model of long-time community behavior, as spatial ecological structure is significantly impacted and can be dominated by migration of organisms, even in small numbers.

Direct Microscopic Observation of Human Neutrophil-Staphylococcus aureus Interaction In Vitro Suggests a Potential Mechanism for Initiation of Biofilm Infection on an Implanted Medical Device.

The ability of human neutrophils to clear newly attached Staphylococcus aureus bacteria from a serum-coated glass surface was examined in vitro using time-lapse confocal scanning laser microscopy. Quantitative image analysis was used to measure the temporal change in bacterial biomass, neutrophil motility, and fraction of the surface area policed by neutrophils. In control experiments in which the surface was inoculated with bacteria but no neutrophils were added, prolific bacterial growth was observed. Neutrophils were able to control bacterial growth but only consistently when the neutrophil/bacterium number ratio exceeded approximately 1. When preattached bacteria were given a head start and allowed to grow for 3 h prior to neutrophil addition, neutrophils were unable to maintain control of the nascent biofilm. In these head-start experiments, aggregates of bacterial biofilm with areas of 50 µm2 or larger formed, and the growth of such aggregates continued even when multiple neutrophils attacked a cluster. These results suggest a model for the initiation of a biofilm infection in which a delay in neutrophil recruitment to an abiotic surface allows surface-attached bacteria time to grow and form aggregates that become protected from neutrophil clearance. Results from a computational model of the neutrophil-biofilm surface contest supported this conceptual model and highlighted the stochastic nature of the interaction. Additionally, we observed that both neutrophil motility and clearance of bacteria were impaired when oxygen tension was reduced to 0% or 2% O2.

Productivity and equilibrium in simple biofilm models.

Biofilms are dense, sessile collections of microorganisms with complicated internal structures. However, in many applications internal details are less important, rather basic, averaged information such as overall community productivity are of most interest. This paper studies averaged community functions in the context of one dimensional, single species, single limiting substrate biofilm models. In particular, using a derived formula for flux of substrate into the biofilm as a function of biofilm height and substrate loading, overall community production can be calculated and system equilibria can be characterized. Consequences for equilibria dependence on a number of mechanisms for balancing growth are considered.

Complex and unexpected outcomes of antibiotic therapy against a polymicrobial infection.

Antibiotics are our primary approach to treating complex infections, yet we have a poor understanding of how these drugs affect microbial communities. To better understand antimicrobial effects on host-associated microbial communities we treated cultured sputum microbiomes from people with cystic fibrosis (pwCF, n = 24) with 11 different antibiotics, supported by theoretical and mathematical modeling-based predictions in a mucus-plugged bronchiole microcosm. Treatment outcomes we identified in vitro that were predicted in silico were: 1) community death, 2) community resistance, 3) pathogen killing, and 4) fermenter killing. However, two outcomes that were not predicted when antibiotics were applied were 5) community profile shifts with little change in total bacterial load (TBL), and 6) increases in TBL. The latter outcome was observed in 17.8% of samples with a TBL increase of greater than 20% and 6.8% of samples with an increase greater than 40%, demonstrating significant increases in community carrying capacity in the presence of an antibiotic. An iteration of the mathematical model showed that TBL increase was due to antibiotic-mediated release of pH-dependent inhibition of pathogens by anaerobe fermentation. These dynamics were verified in vitro when killing of fermenters resulted in a higher community carrying capacity compared to a no antibiotic control. Metagenomic sequencing of sputum samples during antibiotic therapy revealed similar dynamics in clinical samples. This study shows that the complex microbial ecology dictates the outcomes of antibiotic therapy against a polymicrobial infection.

An exclusion principle and the importance of mobility for a class of biofilm models.

Much of the earth's microbial biomass resides in sessile, spatially structured communities such as biofilms and microbial mats, systems consisting of large numbers of single-celled organisms living within self-secreted matrices made of polymers and other molecules. As a result of their spatial structure, these communities differ in important ways from well-mixed (and well-studied) microbial systems such as those present in chemostats. Here we consider a widely used class of 1D biofilm models in the context of a description of their basic ecology. It will be shown via an exclusion principle resulting from competition for space that these models lead to restrictions on ecological structure. Mathematically, this result follows from a classification of steady-state solutions based on a 0-stability condition: 0-stable solutions are in some sense determined by competitive balance at the biofilm base, whereas solutions that are not 0-stable, while less dependent on conditions at the biofilm base, are unstable at the base. As a result of the exclusion principle, it is argued that some form of downward mobility, against the favorable substrate gradient direction, is needed at least in models and possibly in actuality.

Measurements of accumulation and displacement at the single cell cluster level in Pseudomonas aeruginosa biofilms.

Quantitative descriptions of biofilm growth and dynamics at the individual cell level are largely missing from the literature. To fill this gap, research was done to describe growth, accumulation and displacement patterns in developing Pseudomonas aeruginosa biofilms. A parent strain of PAO1 was labelled with either a cyan or yellow fluorescent protein. These were then grown in a flow cell biofilm together so that pockets of dividing cells could be identified and their accumulation and displacement tracked. This analysis revealed a pattern of exponential accumulation for all clusters followed by a stationary accumulation phase. A background 'carpet' layer of cells uniformly colonizing the surface exhibited zero net accumulation of bio-volume. The individual clusters were found to have a mean accumulation rate of 0.34 h(-1) with a range of 0.28-0.41 h(-1). Cluster accumulation rates were negatively correlated with cluster size; larger clusters accumulated volume at a slower rate (P < 0.001). Pockets of cells on the inside of clusters initially accumulated at a comparable rate to the cluster within which they resided, but later invariably exhibited zero to slightly negative accumulation despite continued exponential (positive) accumulation of the cluster. Expanding clusters were able to displace neighbouring cells from the surface, and larger clusters displaced smaller clusters. This work provides a more detailed quantitative experimental observation of biofilm behaviour than has been described previously.

Hypoxia arising from concerted oxygen consumption by neutrophils and microorganisms in biofilms.

Infections associated with microbial biofilms are often found to involve hypoxic or anoxic conditions within the biofilm or its vicinity. To shed light on the phenomenon of local oxygen depletion, mathematical reaction-diffusion models were derived that integrated the two principal oxygen sinks, microbial respiration and neutrophil consumption. Three simple one-dimensional problems were analyzed approximating biofilm near an air interface as in a dermal wound or mucus layer, biofilm on an implanted medical device, or biofilm aggregates dispersed in mucus or tissue. In all three geometries considered, hypoxia at the biofilm-neutrophil interface or within the biofilm was predicted for a subset of plausible parameter values. The finding that oxygen concentration at the biofilm-neutrophil juncture can be diminished to hypoxic levels is biologically relevant because oxygen depletion will reduce neutrophil killing ability. The finding that hypoxia can readily establish in the interior of the biofilm is biologically relevant because this change will alter microbial metabolism and persistence.

Niche partitioning of a pathogenic microbiome driven by chemical gradients.

Environmental microbial communities are stratified by chemical gradients that shape the structure and function of these systems. Similar chemical gradients exist in the human body, but how they influence these microbial systems is more poorly understood. Understanding these effects can be particularly important for dysbiotic shifts in microbiome structure that are often associated with disease. We show that pH and oxygen strongly partition the microbial community from a diseased human lung into two mutually exclusive communities of pathogens and anaerobes. Antimicrobial treatment disrupted this chemical partitioning, causing complex death, survival, and resistance outcomes that were highly dependent on the individual microorganism and on community stratification. These effects were mathematically modeled, enabling a predictive understanding of this complex polymicrobial system. Harnessing the power of these chemical gradients could be a drug-free method of shaping microbial communities in the human body from undesirable dysbiotic states.

Challenges in microbial ecology: building predictive understanding of community function and dynamics.

The importance of microbial communities (MCs) cannot be overstated. MCs underpin the biogeochemical cycles of the earth's soil, oceans and the atmosphere, and perform ecosystem functions that impact plants, animals and humans. Yet our ability to predict and manage the function of these highly complex, dynamically changing communities is limited. Building predictive models that link MC composition to function is a key emerging challenge in microbial ecology. Here, we argue that addressing this challenge requires close coordination of experimental data collection and method development with mathematical model building. We discuss specific examples where model-experiment integration has already resulted in important insights into MC function and structure. We also highlight key research questions that still demand better integration of experiments and models. We argue that such integration is needed to achieve significant progress in our understanding of MC dynamics and function, and we make specific practical suggestions as to how this could be achieved.

Estimation of a biofilm-specific reaction rate: kinetics of bacterial urea hydrolysis in a biofilm.

BACKGROUND/OBJECTIVES: Biofilms and specifically urea-hydrolysing biofilms are of interest to the medical community (for example, urinary tract infections), scientists and engineers (for example, microbially induced carbonate precipitation). To appropriately model these systems, biofilm-specific reaction rates are required. A simple method for determining biofilm-specific reaction rates is described and applied to a urea-hydrolysing biofilm. METHODS: Biofilms were grown in small silicon tubes and influent and effluent urea concentrations were determined. Immediately after sampling, the tubes were thin sectioned to estimate the biofilm thickness profile along the length of the tube. Urea concentration and biofilm thickness data were used to construct an inverse model for the estimation of the urea hydrolysis rate. RESULTS/CONCLUSIONS: It was found that urea hydrolysis in Escherichia coli MJK2 biofilms is well approximated by first-order kinetics between urea concentrations of 0.003 and 0.221 mol/l (0.186 and 13.3 g/l). The first-order rate coefficient (k1) was estimated to be 23.2±6.2 h-1. It was also determined that advection dominated the experimental system rather than diffusion, and that urea hydrolysis within the biofilms was not limited by diffusive transport. Beyond the specific urea-hydrolysing biofilm discussed in this work, the method has the potential for wide application in cases where biofilm-specific rates must be determined.

The molecular dimension of microbial species: 2. Synechococcus strains representative of putative ecotypes inhabiting different depths in the Mushroom Spring microbial mat exhibit different adaptive and acclimative responses to light.

Closely related strains of thermophilic Synechococcus were cultivated from the microbial mats found in the effluent channels of Mushroom Spring, Yellowstone National Park (YNP). These strains have identical or nearly identical 16S rRNA sequences but are representative of separate, predicted putative ecotype (PE) populations, which were identified by using the more highly resolving psaA locus and which predominate at different vertical positions within the 1-mm-thick upper-green layer of the mat. Pyrosequencing confirmed that each strain contained a single, predominant psaA genotype. Strains differed in growth rate as a function of irradiance. A strain with a psaA genotype corresponding to a predicted PE that predominates near the mat surface grew fastest at high irradiances, whereas strains with psaA genotypes representative of predominant subsurface populations grew faster at low irradiance and exhibited greater sensitivity to abrupt shifts to high light. The high-light-adapted and low-light-adapted strains also exhibited differences in pigment content and the composition of the photosynthetic apparatus (photosystem ratio) when grown under different light intensities. Cells representative of the different strains had similar morphologies under low-light conditions, but under high-light conditions, cells of low-light-adapted strains became elongated and formed short chains of cells. Collectively, the results presented here are consistent with the hypothesis that closely related, but distinct, ecological species of Synechococcus occupy different light niches in the Mushroom Spring microbial mat and acclimate differently to changing light environments.

Diel metabolomics analysis of a hot spring chlorophototrophic microbial mat leads to new hypotheses of community member metabolisms.

Dynamic environmental factors such as light, nutrients, salt, and temperature continuously affect chlorophototrophic microbial mats, requiring adaptive and acclimative responses to stabilize composition and function. Quantitative metabolomics analysis can provide insights into metabolite dynamics for understanding community response to such changing environmental conditions. In this study, we quantified volatile organic acids, polar metabolites (amino acids, glycolytic and citric acid cycle intermediates, nucleobases, nucleosides, and sugars), wax esters, and polyhydroxyalkanoates, resulting in the identification of 104 metabolites and related molecules in thermal chlorophototrophic microbial mat cores collected over a diel cycle in Mushroom Spring, Yellowstone National Park. A limited number of predominant taxa inhabit this community and their functional potentials have been previously identified through metagenomic and metatranscriptomic analyses and in situ metabolisms, and metabolic interactions among these taxa have been hypothesized. Our metabolomics results confirmed the diel cycling of photorespiration (e.g., glycolate) and fermentation (e.g., acetate, propionate, and lactate) products, the carbon storage polymers polyhydroxyalkanoates, and dissolved gasses (e.g., H2 and CO2) in the waters overlying the mat, which were hypothesized to occur in major mat chlorophototrophic community members. In addition, we have formulated the following new hypotheses: (1) the morning hours are a time of biosynthesis of amino acids, DNA, and RNA; (2) photo-inhibited cells may also produce lactate via fermentation as an alternate metabolism; (3) glycolate and lactate are exchanged among Synechococcus and Roseiflexus spp.; and (4) fluctuations in many metabolite pools (e.g., wax esters) at different times of day result from species found at different depths within the mat responding to temporal differences in their niches.

General theory for integrated analysis of growth, gene, and protein expression in biofilms.

A theory for analysis and prediction of spatial and temporal patterns of gene and protein expression within microbial biofilms is derived. The theory integrates phenomena of solute reaction and diffusion, microbial growth, mRNA or protein synthesis, biomass advection, and gene transcript or protein turnover. Case studies illustrate the capacity of the theory to simulate heterogeneous spatial patterns and predict microbial activities in biofilms that are qualitatively different from those of planktonic cells. Specific scenarios analyzed include an inducible GFP or fluorescent protein reporter, a denitrification gene repressed by oxygen, an acid stress response gene, and a quorum sensing circuit. It is shown that the patterns of activity revealed by inducible stable fluorescent proteins or reporter unstable proteins overestimate the region of activity. This is due to advective spreading and finite protein turnover rates. In the cases of a gene induced by either limitation for a metabolic substrate or accumulation of a metabolic product, maximal expression is predicted in an internal stratum of the biofilm. A quorum sensing system that includes an oxygen-responsive negative regulator exhibits behavior that is distinct from any stage of a batch planktonic culture. Though here the analyses have been limited to simultaneous interactions of up to two substrates and two genes, the framework applies to arbitrarily large networks of genes and metabolites. Extension of reaction-diffusion modeling in biofilms to the analysis of individual genes and gene networks is an important advance that dovetails with the growing toolkit of molecular and genetic experimental techniques.

A multidimensional multispecies continuum model for heterogeneous biofilm development.

We propose a multidimensional continuum model for heterogeneous growth of biofilm systems with multiple species and multiple substrates. The new model provides a deterministic framework for the study of the interactions between several spe1cies and their effects on biofilm heterogeneity. It consists of a system of partial differential equations derived on the basis of conservation laws and reaction kinetics. The derivation and key assumptions are presented. The assumptions used are a combination of those used in the established one dimensional model, due to Wanner and Gujer, and for the viscous fluid model, of Dockery and Klapper. The work of Wanner and Gujer in particular has been extensively used through the years, and thus this new model is an extension to several spatial dimensions of an already proven working model. The model equations are solved using numerical techniques, for purposes of simulation and verification. The new model is applied to two different biofilm systems in several spatial dimensions, one of which is equivalent to a system originally studied by Wanner and Gujer. Dimensionless formulations for these two systems are given, and numerical simulation results with varying initial conditions are presented.

Interaction between arthropod filiform hairs in a fluid environment.

Many arthropods use filiform hairs as mechanoreceptors to detect air motion. In common house crickets (Acheta domestica) the hairs cover two antenna-like appendages called cerci at the rear of the abdomen. The biomechanical stimulus-response properties of individual filiform hairs have been investigated and modeled extensively in several earlier studies. However, only a few previous studies have considered viscosity-mediated coupling between pairs of hairs, and only in particular configurations. Here, we present a model capable of calculating hair-to-hair coupling in arbitrary configurations. We simulate the coupled motion of a small group of mechanosensory hairs on a cylindrical section of cercus. We have found that the coupling effects are non-negligible, and likely constrain the operational characteristics of the cercal sensory array.

Adaptive responses to antimicrobial agents in biofilms.

Bacterial biofilms demonstrate adaptive resistance in response to antimicrobial stress more effectively than corresponding planktonic populations. We propose here that, in biofilms, reaction-diffusion limited penetration may result in only low levels of antimicrobial exposure to deeper regions of the biofilm. Sheltered cells are then able to enter an adapted resistant state if the local time scale for adaptation is faster than that for disinfection. This mechanism is not available to a planktonic population. A mathematical model is presented to illustrate. Results indicate that, for a sufficiently thick biofilm, cells in the biofilm implement adaptive responses more effectively than do freely suspended cells. Effective disinfection requires applied biocide concentration that increases quadratically or exponentially with biofilm thickness.