Kenneth J. Locey
Melany C. Fisk
J. T. Lennon- 1Department of Biology, Indiana University, Bloomington, IN, USA
- 2Department of Biology, Miami University, Oxford, OH, USA
Microbial dormancy leads to the emergence of seed banks in environmental, engineered, and host-associated ecosystems. These seed banks act as reservoirs of diversity that allow microbes to persist under adverse conditions, including extreme limitation of resources. While microbial seed banks may be influenced by macroscale factors, such as the supply of resources, the importance of microscale encounters between organisms and resource particles is often overlooked. We hypothesized that dimensions of spatial, trophic, and resource complexity determine rates of encounter, which in turn, drive the abundance, productivity, and size of seed banks. We tested this using >10,000 stochastic individual based models (IBMs) that simulated energetic, physiological, and ecological processes across combinations of resource, spatial, and trophic complexity. These IBMs allowed realistic dynamics and the emergence of seed banks from ecological selection on random variation in species traits. Macroscale factors like the supply and concentration of resources had little effect on resource encounter rates. In contrast, encounter rates were strongly influenced by interactions between dispersal mode and spatial structure, and also by the recalcitrance of resources. In turn, encounter rates drove abundance, productivity, and seed bank dynamics. Time series revealed that energetically costly traits can lead to large seed banks and that recalcitrant resources can lead to greater stability through the formation of seed banks and the slow consumption of resources. Our findings suggest that microbial seed banks emerge from microscale dimensions of ecological complexity and their influence on resource limitation and energetic costs.
Introduction
Most microorganisms live in environments where they experience energy limitation, resource limitation, or both (Hoehler and Jørgensen, 2013; Moore et al., 2013). Microorganisms have evolved an expansive repertoire of traits that allow them to persist under conditions of resource scarcity (Follows and Dutkiewicz, 2011; Lever et al., 2015; Litchman et al., 2015). One strategy that is important for microorganisms that experience resource limitation is dormancy, i.e., a reversible state of reduced metabolic activity (Jones and Lennon, 2010; Aanderud et al., 2016). Dormant microorganisms make up a seed bank, which contributes to the maintenance of diversity (Lennon and Jones, 2011; Aanderud et al., 2015) and the functioning of ecosystems (Wang et al., 2015). Transitions into and out of dormancy are often driven by the availability of energy and nutrients (Lennon and Jones, 2011); yet seed banks still accumulate in otherwise resource-rich environments. For example, >90 % of microbial biomass in soils can be dormant even though organic matter content in these habitats can be quite high (Alvarez et al., 1998; Lennon and Jones, 2011; Blagodatskaya and Kuzyakov, 2013). Therefore, seed-bank dynamics may be influenced by factors other than macroscale properties such as the bulk concentration or supply rate of resources.
In an idealized system with few trophic interactions and where labile substrates are homogenously distributed, encounter rates between individual microorganisms and resources are governed by physical processes, such as turbulence and diffusion (Dusenbery, 2009; Rusconi and Stocker, 2015). However, these idealized well-mixed conditions are rarely met in nature. Instead, microorganisms live in complex habitats where aggregated particles of many resource types can vary in size, energetic yield, and spatial distribution (e.g., Hernández and Hobbie, 2010; Macalady et al., 2013). Such complexities modify the rate at which microorganisms encounter consumable resource particles (Kiørboe et al., 2002; Andersen et al., 2016; Großkopf and Soyer, 2016). Because it is challenging to integrate this fine-scale complexity into empirical studies, microorganisms are often investigated at spatial scales that exceed the scales of their individual interactions (Fierer and Lennon, 2011; Vos et al., 2013). For this reason, microorganisms may be nutrient- or energy-limited even though macroscale measurements would suggest that their habitat is replete with resources (Don et al., 2013; Allison et al., 2014). This phenomenon has led to the hypothesis that there is an advantage to maintaining large but inactive populations (i.e., seed banks) that are able to maximize the probability of encountering resources that vary in time or space (Vaqué et al., 1989).
The development of microbial seed banks may be influenced in part by microbe-resource encounter rates. These encounter rates are likely driven by interacting dimensions of ecological complexity. For example, microorganisms have highly aggregated spatial distributions in physically structured habitats (Raynaud and Nunan, 2014) but also in seemingly well-mixed systems (Azam and Malfatti, 2007). Such patterns may reflect the non-random distribution of resources and the capacity of microorganisms to disperse (Mitchell and Kogure, 2006; Smriga et al., 2016). Encounter rates may also be affected by the resource pool, which often comprises diverse substrates with complex molecular structures (Muscarella et al., 2014; Logue et al., 2016). Some of these resources may only be accessible to specialized taxa that produce extracellular enzymes (Lennon, 2007), which require energy that could otherwise be used for maintenance and growth (Traving et al., 2015). Last, encounters between microorganisms and resources may be influenced by trophic interactions, such as competition, predation, and parasitism (e.g., Hibbing et al., 2010). Resource-limited microorganisms also engage in the consumption of dead cells (i.e., scavenging) and the division of labor that arises from the sharing of metabolites with neighboring cells (i.e., cross-feeding) (Rozen et al., 2009; Pande et al., 2015). Together, varying degrees of spatial, trophic, and resource complexity should influence encounter rates in ways that drive growth, abundance, and activity of microbial communities, in part, through the emergence of seed banks.
Studying complex interactions at the microscale is a profound challenge for microbial ecology (Cordero and Datta, 2016). Individual-based models (IBMs) provide a way to address this challenge by exploring how individual-level interactions give rise to higher-order phenomena at the scale of populations, communities, and ecosystems (Grimm et al., 2005; Hellweger et al., 2016). IBMs can include degrees of realism and complexity that are unattainable with other types of modeling, and which lead to the emergence of unexpected but insightful phenomena (DeAngelis and Gross, 1992; Grimm, 1999; Grimm et al., 2005; Rosindell et al., 2015). In this study, we developed an IBM framework that explicitly simulates the physiology, life history, energetics, and metabolic activity of microorganisms while exploring dimensions of resource, spatial, and trophic complexity. Using randomized model parameters and random combinations of species traits, and then iterating stochastic life history processes over thousands of generations of ecological selection, we expected microbial seed banks and dynamics of growth and abundance to emerge from individual-level interactions. We show how microscale aspects of resource, spatial, and trophic complexity can drive encounter rates which, in turn, influence aspects of microbial community structure, including abundance, productivity, and seed bank dynamics. Finally, our models provide evidence of the functions that seed banks perform in energy- and nutrient-limited environments.
Methods
Overview of Individual-Based Modeling
We constructed an automated source-code that built and ran stochastic individual-based models (IBMs) starting with random values of resource supply and species traits as well as random combinations of spatial, resource, and trophic complexity. These IBMs simulated physical encounters between organisms and resource particles of realistic size within spatially explicit three-dimensional environments. The purpose of this computationally demanding approach was to examine how spatial, resource, and trophic complexity affect cell-resource encounter rates and how encounter rates influence the abundance, productivity, and activity of microbial communities. In the following sections we describe how the models were parameterized with random species-, resource-, and individual-level parameters, and how the models simulated (1) realistic microscales of space and time, (2) organisms, species, and resource particles, (3) life history processes and their energetic costs as well as encounter-limited growth, and (4) levels of spatial complexity, resource complexity, and trophic complexity.
Randomized Model Parameterization
The parameter values of each IBM were randomly chosen to fulfill several requirements (Table 1). These parameters included upper limits on energetic constraints (e.g., maintenance energy) and species vital rates (e.g., dispersal, growth), and also included resource-related variables (e.g., inflow rate, inflowing resource diversity), as well as random combinations of spatial, trophic, and resource complexity (Table 1). Once assembled, each IBM was populated with 100 individuals whose species identities, vital rates, resource use, and maintenance energies were drawn at random (Table 1).
Table 1. Parameter values for individual based models (IBMs).
Our extensive randomization accomplished several objectives. First, the IBMs were able to generate diverse trait combinations and combinations of ecological complexity across thousands of models. Second, the IBMs were able to simulate realistic but variable capacities for dormancy and realistic relationships between rates of reproduction and rates of dispersal. Third, the ranges of model parameters naturally limited organism abundance and the total number of resource particles to less than 10,000, which was necessary for practical computational overhead.
In addition to allowing traits, trade-offs, and life-history strategies to emerge as a result of ecological selection operating over thousands of generations, our randomized parameterization and stochastic simulation addressed a critical challenge to ecological modeling. Specifically, when traits, trade-offs, or life-history strategies are enforced within a model, they cannot provide evidence of the forces that produce them. However, the emergence of traits, trade-offs, and life-history strategies from initially random conditions within our models would provide evidence about the importance of microscale factors in driving rates of encounter between individual organisms and resources and, in turn, the importance of microbe-resource encounter rates in driving abundance, productivity, and the emergence of seed banks.
Spatial and Temporal Scale
We chose spatial and temporal scales that were relevant to individual microorganisms, but which also established a space-time equivalency between rates of reproduction and dispersal. For example, one of the most well-studied microorganisms (Escherichia coli) is believed to disperse up to 36 microns (μm) per second (Milo and Phillips, 2015). Likewise, under ideal conditions, E. coli is able to double every 20 min (Wang et al., 2010). Using 1200 seconds (20 min) as a minimum doubling time (i.e., nothing doubling faster than E. coli) and 36 μm per second as a maximum specific dispersal rate, our models obeyed a space-time equivalency such that each time step allowed, at most, one doubling per individual for the fastest reproducing species or 36 μm of dispersal for the most quickly dispersing species. We then constrained each time step to 0.12 s, translating to 1200 simulated real-world seconds and, hence, a three-dimensional spatial extent of 43,200 μm in each direction for a total volume of 8·104 cm3. For context, densities of bacteria can reach several billion per cm3 in resource-rich environments, such as soils, sediments, and gut contents (Whitman et al., 1998). In our modeling, we aimed to simulate at most a few tens of thousands of individual organisms and resource particles, primarily due to computational limits but also to simulate sparsely populated environments that were highly limited in energy-yielding resources.
Simulating Individuals, Species, and Resources
Our modeling simulated different types of resource particles that could be consumed by heterotrophic microorganisms belonging to different species. We used data objects (e.g., lists and dictionaries) that are native to most computing languages. Individuals were distinguished by collections of elements within dictionaries, which are data objects that use key-value pairs to assign values (e.g., growth rate) to a “key” (e.g., individual ID). For example, a dictionary of information for individual “1” could be as follows: {“quota”: 0.5, “x”: 50, “y”: 100}. Here, “quota” refers to the individual's cell quota, i.e., amount of endogenous resource. This individual would have half of its maximum cell quota (i.e., 0.5), and would be located at the x-y coordinates of 50–100. Each individual and its dictionary of information was then held in a community-level dictionary (INDs):
Species-specific information including specific rates of growth and dispersal, values of cellular maintenance, probability of randomly reactivating, resource specificity, etc., were held in a community-level species dictionary (SPs):
In simulating an ecological process (e.g., growth, dispersal, cell maintenance) a model would choose an individual (e.g., 1, 2, …) at random and then use the information for the individual in INDs and its species to determine how much the individual grows when it encounters a resource, how fast it disperses, how much cell quota is lost to maintenance, etc. Like individuals, resource particles were distinguished by collections of elements within dictionaries. When resource particles entered the system, they were assigned a resource identity (e.g., a, b, c) at random, a value between 1000 and 10,000 representing the size of the resource particle, as well as a set of three-dimensional spatial coordinates. Each particle was also given a unique identity (ID) for referencing the resource dictionary.
More information on the simulation of individual organisms, resource particles, and all other aspects our modeling can be found in a highly detailed version of the standard IBM protocol of Grimm et al. (2006), which is available as a supplemental file (see Supplemental File 1) and in a public GitHub repository (https://github.com/LennonLab/Micro-Encounter).
Simulating Life History Processes
Our IBMs used uniform random sampling to choose which individuals at a particular moment would undergo specific life history processes, such as growth, reproduction, death, cellular maintenance, dispersal, and transitions into an out of dormancy. Whether a randomly chosen individual underwent a specific life history process was determined, in part, by the individual's cell quota, species-specific trait values, and distance to resource particles. Modeling in this probabilistic way simulated the partly probabilistic and partly deterministic nature of environmental filtering and individual-level interactions. The order in which processes occurred was allowed to vary at random from time step to time step, which prevented artifacts from arising that might be driven by the ordering of simulated life history processes. Specifically, we implemented life history process as follows:
Immigration
Individuals entered the system at any point in the three-dimensional environment. Species identities of inflowing propagules were chosen at random from a uniform distribution. Individuals were assigned a unique ID and a randomly chosen cell quota between 0.5 and 1.0, which translated to starting cell diameters of 1 to 1.25 μm.
Consumption, Growth, and Cell Size
Consumption increased an individual's cell quota and decreased the size of the resource that an individual had encountered. Growth was simulated as an increase in cell quota that corresponded to a proportional increase in biomass. That is, we assumed that cell quota (Q) and cell volume (V) were proportional and hence, that Q was related to a cell's radius through a geometric relationship: Q ∝ V = 4/3 π r3. Estimating cell size was important as the ability to consume a resource particle relied on whether a cell and a resource particle were in physical contact. The energetic cost of growth was proportional to the product of growth rate and individual cell quota, such that larger individuals required more energy to increase by a fraction of their body mass.
Physiological Maintenance
Individual cell quotas diminished according to a randomly chosen species-specific maintenance cost, i.e., a numeric value subtracted from the cell quota. This cost was the implied energetic cost of maintaining basal metabolism and cellular maintenance. When active individuals transitioned to a dormant state, the cell maintenance cost decreased by a species-specific factor that was likewise chosen at random. This cell maintenance cost was reversed when individuals resuscitated to a metabolically active state.
Reproduction
Reproduction was clonal and occurred via binary fission and without mutation. Individuals reproduced with a probability directly determined by the ratio of cell quota to maximum cell quota (i.e., 1000). Hence, the probability of reproducing when randomly sampled was between 0.0 and 1.0. The daughter cell was given a unique individual ID along with the species ID and geographical coordinates of its mother cell.
Death
An individual sampled at random died if its cell quota was less than the species-specific maintenance energy, whether in an active or dormant state. Depending on the model, the remains of dead individuals could be scavenged and consumed by other cells. Otherwise, dead individuals were effectively lost to the system, as in most community ecology models.
Simulating Resource Breakdown and Consumption
Resource particles entered the system with diameters ranging between 5000 and 20,000 μm. When encountered, individuals would consume the resource according to species specific rates of consumption. The remaining portion of the resource particle would then be randomly broken in two parts. Below, we describe how, depending on the model, the specific details of resource breakdown and consumption varied with regard to whether the model included scavenging, cross-feeding, or resources of varying complexity or recalcitrance, i.e., the difficulty of being broken down.
Simulating Ecological Complexity
In addition to explicitly simulating the spatial environment and individual-level changes in organisms and resource particles, we constructed our modeling framework to allow random combinations of various levels of spatial, resource, and trophic complexity. Each IBM was parameterized at random with one of 72 combinations of complexity (4 trophic, 3 resource, 3 dispersal, 2 spatial resource distribution) (Figures 1, 2). We then explored how these dimensions of complexity affected encounter rates, along with attributes, such as total abundance (i.e., total number of individual in a community; N), the abundance of active individuals, production of individuals per time step (i.e., productivity), and the relative size of the microbial seed bank (i.e., percent dormancy).
Figure 1. Conceptual diagram depicting conditions of spatial complexity simulated in individual based models (IBMs). Columns: The spatial distribution of resources and organisms was modeled in two ways. (Left) Locations of resource particles and individual organisms changed randomly between time steps, resulting in a well-mixed model of uncorrelated changes in position. (Right) Aggregated resource particles remained largely fixed in position and only made slight movement under Brownian motion, resulting in a highly structured and spatially heterogeneous environment. Rows: Whether in a well-mixed or structured environment, individuals could change position under three modes of dispersal. (Top) Individuals could disperse passively (no active dispersal), relying entirely on either the well-mixed conditions or Brownian motion imposed by their environment. (Center) Individuals could actively disperse along a trajectory and change direction upon hitting a barrier or edge, in a “run and tumble” fashion. (Bottom) Most energetically expensive, individuals could actively disperse in a directed manner after sensing the presence of consumable resources, i.e., chemotaxis. In the modeling, each condition of resource distribution was compatible with all conditions of individual dispersal, resulting in 6 possible combinations.
Figure 2. Conceptual diagram depicting aspects of resource complexity and trophic complexity simulated in individual based models (IBMs). Resource complexity (top row) was simulated in three ways. (Left) All resources supplied from outside the system to an individual microbe (solid red circle) were of the same type (open red circles), resulting in monoculture conditions. (Center) All resources supplied from outside the system to microbial taxa (solid red circles, solid blue triangles, solid green squares) were of three different types (open blue triangles, open green squares, open red circles) supplied in equal proportions and resulting in a polyculture. (Right) Resources required energetically costly efforts to breakdown and use, much like an enzymatic “lock-and-key” in the breakdown of complex macromolecules. Models were defined by either monoculture or polyculture resource conditions, and also had either “lock-and-key” resource complexity or simple resources that were immediately assimilated without an energetic cost. Trophic complexity (bottom row) was also simulated in three ways. (Left) A consumer-resource condition where all species belonged to the same trophic level and where individuals only consumed externally supplied resources. (Center) A scavenging condition where the remains of dead individuals were consumed by any other individual in the system that was in close enough physical proximity. (Right) A cross-feeding condition where individuals produced metabolites or by-products that were then consumed by individuals belonging to another taxon. In the modeling, cross-feeding could emerge as a one-way relationship (pictured) or as a two-way relationship.
Spatial Complexity
We simulated levels of spatial complexity in terms of the spatial structure of the environment and modes of individual dispersal (Figures 1, 2). The first level of spatial complexity in the environment was a well-mixed model in which the locations of individual organisms and resource molecules changed at each time step in an uncorrelated way. Hence, every organism and resource particle had the same chance of moving to any location within the environment at each time step in the model. While this is unrealistic in highly structured systems of relatively large size, the microscale spatial dynamics of more fluid systems can be highly complex, rapidly mixed, and unpredictable (Rusconi et al., 2014). This model created a well-mixed environment and allowed for passive dispersal. In the second level of spatial complexity, resource molecules only changed their location via Brownian motion, resulting in a system that lacked mixing and where, through the breakdown of large resource particles, smaller resource particles became aggregated. We refer to this second level as “spatially structured.”
The first mode of individual dispersal excluded active dispersal and was entirely passive. While passive dispersal has the benefit of not requiring an energetic investment, we expected this mode of dispersal to only have a benefit under well-mixed conditions. The second mode of dispersal allowed active movement, but did not allow chemoreception. This “run and tumble” strategy allowed individuals to move along a straight line in a randomly chosen direction and, once hitting the edge of the environment, turned and moved in a new, randomly chosen direction. By chance alone, individuals would be able to encounter resource particles. The third mode of dispersal simulated chemotaxis, i.e., the ability to sense the location of resources and actively move toward them. Because chemotaxis involves sense perception, we encoded it as a more energetically expensive mode of dispersal. In this way, one could envision a potential trade-off between the greater energetic cost of chemotaxis and the efficiency associated with more rapid encounter of resources.
Resource Complexity
Our simulations captured important features of resource complexity that we hypothesized would affect encounter rates. First, we considered resource diversity, which refers to the different types of resources that are supplied to a system. In the “monoculture” simulation, only one type of resource was supplied. Second, we considered a “polyculture” scenario, where three different types of resources were supplied and each could only be used by specialist consumers (Figure 1). Third, we considered that some resource particles are more difficult to break down, or recalcitrant, than others. To simulate recalcitrance, we imposed a “lock and key” scenario that required individuals to invest time and energy to break down complex molecules for consumption. Importantly, a given model could include this “lock and key” feature while being a monoculture or polyculture.
Trophic Complexity
We simulated three aspects of trophic complexity that we hypothesized would influence resource encounter (Figure 1). First, we simulated a simple heterotrophic “consumer-resource” interaction where all individuals were solely consumers of inflowing resources. The second level of trophic complexity allowed for “scavenging” which is specified in the model as the consumption of resources from dead individuals (e.g., Rozen et al., 2009). Our third level of trophic complexity simulated a situation wherein individuals produced metabolites as a result of breaking down and consuming resources. These metabolites could potentially be consumed by other species that, in turn, produced metabolites that also be consumed. This situation was meant to simulate conditions that are characteristic of cross-feeding or syntrophy (Pande et al., 2015). Additionally, because the metabolites that individuals of particular species produced were chosen at random at the start of the model, it was possible for mutualistic cross-feeding to occur, i.e., where two species produced metabolites that the other could consume.
Modeling Workflow
Each model was run to a state of mean reversion, i.e., a point where the total abundance of the microbial community fluctuated around a given value, as the model iterated over thousands of generations of randomly ordered life history processes. This burn-in period was then discarded and each model was run for 2000 additional generations. We recorded information for each tenth generation past the burn-in period (Table 1).
Results and Discussion
Dormancy allows microorganisms to persist in low-resource environments, yet seed banks emerge under other conditions as well (Lennon and Jones, 2011; Blagodatskaya and Kuzyakov, 2013). While it is generally assumed that nutrient and energy limitation can lead to seed banks, the interrelated variables and dimensions of ecological complexity that influence the transition of individuals into dormancy are poorly understood and are difficult to study in natural environments. Our results from >10,000 individual-based models (IBMs) suggest that microscale factors modify microbe-resource encounters in ways that influence fundamental properties of microbial communities. For example, dispersal, spatial structure, and resource complexity had a greater influence on encounter rates than macroscale properties, such as the supply rate and bulk concentration of resources in the environment. In turn, these microscale drivers served as important controls on microbial community attributes, such as abundance, productivity, and seed banks.
Realistic Model Behavior
Not even the most sophisticated models can accurately simulate all factors affecting growth, abundance, and activity in complex ecological communities. Nevertheless, our IBMs reproduced several realistic features of microbial systems. First, resource particles were broken down to mean particle sizes of 3.5 to 412 μm. These ranges approximate the sizes of marine snow aggregates (Alldredge, 1998), macro- and micro-aggregates of soil (Puget et al., 2000), and suspended and sedimented detritus in marsh ecosystems (Marsh and Odum, 1979). Second, mean cell diameters ranged from 0.2 to 2.0 μm, which approximates the range of sizes from the extremely small marine bacteria SAR11 (~0.2 μm) (Rappé et al., 2002) to the relatively large diameter (0.4–0.8 μm) and length of E. coli (1.0–3.5 μm) (Trueba and Woldringh, 1980). Third, we found that the ecological model that has been shown to best explain patterns of commonness and rarity among microbial communities from aquatic, terrestrial, host-related systems often explained >80% of variation among species in our IBMs (Figure S1; also see Shoemaker et al., 2016). Fourth, as discussed in the following sections, we observed intuitive relationships of growth, abundance, and dormancy. These relationship reinforce reasonable model behavior and provide insights into how microscale complexity drives encounter and how encounter rates drive growth, abundance, and the emergence of microbial seed banks.
Microscale vs. Macroscale Drivers of Resource Encounter
The consumption of resources needed for growth and reproduction are regulated by factors that influence encounter rates. While our IBMs reflected this basic assumption, we found that encounter rates were differentially influenced by macroscale and microscale properties of the simulated environment. For example, there was no relationship between encounter rate and the resource supply rate (Figure S1) and though encounter rate increased with resource concentration, the relationship was noisy and complicated (Figure S2). We also found that the diversity of inflowing resources, which ranged from 1 to 10 different types, had little-to-no influence on encounter rates (Figure 3). In contrast, several microscale properties had substantial influence on encounter rates, including dispersal mode, the recalcitrance of individual resources, and to a lesser degree, trophic complexity (i.e., scavenging, cross-feeding) (Figure 3). We discuss the nature of these microscale factors on microbe-resource encounter rates in the following subsections.
Figure 3. Box plots revealing the influence resource complexity (left), trophic complexity (right), spatial distributions (top and bottom), and dispersal conditions (different colors) on frequencies of encounter between organisms and resource particles. Top row: In well-mixed environments, “lock-and-key” chemical complexity (in combination with either monoculture or polyculture conditions) had a strong influence on encounter rate, resulting in a substantial decrease. Monoculture and polyculture box plots did not include “lock-and-key” resource complexity. In the top row, (i.e., well-mixed environments) chemotaxis and “run and tumble” do not substantially change encounter rates. Bottom row: In structured environments, active dispersal overrode the influence of chemical complexity, with chemotaxis providing the greatest increase in encounter.
Active Dispersal
Encounter rates were most affected by interactions between dispersal mode and the spatial structure of the environment (Figure 3). In well-mixed environments, neither chemotaxis nor “run and tumble” produced greater encounters than the energetically free strategy of passive dispersal. However, in structured environments, chemotaxis produced 10–100 times more encounters while passive dispersal frequently resulted in no encounters. Chemotaxis is an important microbial trait that has strong effects on resource encounter rates (e.g., Datta et al., 2016; Smriga et al., 2016). The energetic cost of this directed movement must have been offset by the energy saved in not encountering resources at random, as with our simulated “run and tumble” mode of dispersal. It is generally assumed that chemotaxis is a trait that is favored in productive habitats and selected against in oligotrophic environments where searching for sparse chemically complex resources could be energetically wasteful (Ottemann and Miller, 1997). This view is supported by genomic evidence suggesting that bacteria from productive environments tend to harbor more motility genes than bacteria from oligotrophic environments (Giovannoni et al., 2005; Lauro et al., 2009). Though our models did not support these observations, the energetic cost of chemotaxis, as modeled, may not have been great enough to substantially tax individual cell quotas. More likely, however, individuals that became energetically depleted through chemotaxis simply went dormant.
Spatial Distribution
Habitat heterogeneity is thought to limit the probability of decomposition and to influence the accumulation of resources in spatially structured environments like soil (Schmidt et al., 2011; Dungait et al., 2012; Lehmann and Kleber, 2015). Our simulations support this view. Except in the case of chemotaxis, we found that encounter rates were lower in spatially structured environments, which resulted in the aggregation and persistence of resources in the environment (Figure 3). This pattern is most clearly seen in the comparisons between mixed and structured environments, where the “run and tumble” strategy or passive dispersal resulted in greatly reduced encounter rates (Figure 3). Our results highlight the benefits of chemotaxis in overcoming the challenges of resource acquisition in spatially structured environments. On the other hand, in well-mixed environments, chemotaxis did not result in greater encounters compared to the energetically less expensive “run and tumble” mode of dispersal and the energetically free passive dispersal strategy. These model-based findings suggest the importance of considering energetic trade-offs related to modes of dispersal and the spatial distribution of resources.
Resource Complexity
The recalcitrance of resources reduced encounter rates in both mixed and structured environments (Figure 3). In order for microorganisms to consume complex resources, they needed to invest energy in cleaving resource particles into consumable parts, slowing encounter with those breakdown products. This finding confirms how structurally complex or recalcitrant resources are thought to influence the growth and activity of microorganisms within habitats, such as soil (Schimel and Weintraub, 2003; Allison et al., 2011). Along with the spatial distribution of resources, chemical complexity slows the loss of energy from the resource pool by limiting encounter. This “slow release” effect has been hypothesized to affect ecosystem dynamics (e.g., Wetzel, 1999) and its interaction with spatial structure is consistent with emerging views on the controls of organic matter persistence in soil environments (Lehmann and Kleber, 2015). Furthermore, time series of our simulated data suggest that the slow turnover of recalcitrant substrates contributes to a more stable dynamic of resource encounter (Figure 4). Whereas the inclusion of complex resources (i.e., “lock-and-key” models) had little-to-no influence on average encounter in well-mixed models, the presence of recalcitrant resources did impart greater stability in encounter rates over time, regardless of dispersal mode or trophic complexity (Figure 4).
Figure 4. Time series from structured environments reveal a stabilizing influence of “lock-and-key” resource chemical complexity on microbe-resource encounters. Both “monoculture” and “polyculture” hulls (75% confidence intervals) pertain to simple (i.e., labile) resources, and have higher variability in encounter rates than the “lock-and-key” models. This was the case across all modes of dispersal and all modeled forms of trophic complexity (i.e., simple consumer-resource model, cross-feeding, scavenging).
Trophic Complexity
Our simulations revealed that trophic complexity had only modest effects on encounter rates (Figure 3). Both scavenging and cross-feeding led to slight increases in resource encounter in both well-mixed and structured environments across all modes of dispersal. However, the small effect of trophic complexity on resource encounter should be cautiously interpreted. First, we only considered a few types of trophic interactions and we encoded them in fairly specific forms. In addition to scavenging and cross-feeding, microbial communities engage in a plethora of trophic interactions, which could affect encounter and consumer-resource interaction strength. Second, aspects of our modeling may have dampened the effects of scavenging and cross-feeding on encounter rates. For example, because our models included dormancy which reduces mortality, seed banks may have reduced the importance of scavenging by decreasing available necromass. Finally, it is possible that scavenging and cross-feeding affected resource encounter for a small cross-section of models but that the signal was overshadowed by the larger effects of spatial and resource complexity. We envision two ways that future studies could more closely examine the roles of trophic complexity in driving microbe-resource encounters. First, studies could explore the influence of trophic levels or ask whether trophic interactions that lower the abundances of primary consumers (e.g., predator-prey, parasitism) drive resource encounters. Second, the possibility that scavenging could increase the availability of resources to organisms warrants greater attention to the energetic value derived from necromass.
Influence of Encounter on Microbial Dynamics
Aspects of abundance, productivity, and activity were strongly related to encounter rates (Figure 5). We found that microbial seed banks emerged from multiple mechanisms related to resource-mediated energetic limitations. Finally, both microbial seed banks and a slow release effect of recalcitrant resources may stabilize microbial systems against large fluctuations in growth and resource availability.
Figure 5. Heat maps revealing a negative relationship of encounter to the relative size of the seed bank (% dormancy) and positive relationships of encounter to the production of new individuals (productivity), and to the abundance of the total and active communities. In each plot, blue heat maps represent the results of models that included “lock-and-key” chemical complexity (complex molecules) and excluded active dispersal. Red heat maps represent the results of models that excluded chemical complexity (simple molecules) but included chemotaxis. Greater heat, i.e., areas within heat maps that have lighter colors, corresponds to a greater number of model results.
Abundance, Productivity, and Emergence of Seed Banks
Greater resource encounters led to greater productivity and greater total abundance, and smaller seed banks (Figure 5). While these relationships can be expected since resource encounters are required to fuel growth, we also observed that modes of dispersal influenced the general forms of these relationships. For example, in spatially structured models, chemotaxis not only led to greater encounters but also some of the largest seed banks. In contrast, the “run and tumble” strategy produced the overall smallest seed banks despite encounter rates that were often very low (Figure 5). We attribute the stark difference between these two strategies to their energetic cost, where chemotaxis was twice as energetically costly as “run and tumble”. Resource complexity and trophic complexity had no effect on relationships of abundance and productivity to rates of encounter regardless of whether the environment was structured or well-mixed.
Consistent with energy- and nutrient limitation as a control on microbial activity, we observed that seed banks were largest in simulations that included recalcitrant resources and passive dispersal (Figure 6). In contrast, the relative size of seed banks was not related to the macroscale properties of total resource supply and bulk resource concentration (Figure 6). Finally, we often found that passively dispersing individuals did not encounter resource particles (Figure 7). On closer inspection, this appeared to be due to the small size and sparse spatial distribution of many resource particles, i.e., “crumbs,” which has been suggested as an alternative to resource recalcitrance as a reason for why large amounts of total resources seem to be resistant to microbial consumption (Arrieta et al., 2015). In soil systems, low concentrations of organic matter have been found to limit its use by microbes (Don et al., 2013), presumably because dilute resources limit the encounter with organisms. Consequently, it may be that natural systems in which we find large seed banks are likely to be those characterized by recalcitrant or chemically complex resources, such as plant structural materials, as well as sparse distributions of relatively small resource particles.
Figure 6. Heat maps reveal that the relative size of the seed bank (% Dormant) was strongly related to the microscale property of encounter but was independent of the macroscale properties of total resources and resource supply. In each plot, blue heat maps represent the results of models that included chemical complexity (recalcitrant; complex molecules) and excluded active dispersal. Red heat maps are results of models that excluded chemical complexity (labile; simple molecules) but included chemotaxis. Greater heat, i.e., areas within heat maps that have lighter colors, corresponds to a greater number of model results.
Figure 7. Time series of structured environments reveal the influence of dispersal mode (i.e., chemotaxis, “run and tumble,” no active dispersal) on seed bank size (% Dormant), total abundance (N), total resources (R), and the number of encounters. Chemotaxis generally produced the largest seed banks while “run and tumble” produced the smallest seed banks, even though chemotaxis resulted in the greatest encounters; a potential consequence of the energetic cost of chemotaxis. Under passive dispersal, the community was largely dormant, small in abundance, and experienced low encounters even as resource particles accumulated.
Influence of Seed Banks on Temporal Dynamics
Time series from our IBMs revealed important influences of microbial seed banks across 2000 generations. Active dispersal and a strong capacity for dormancy helped maintain large and stable communities over time (Figures 7, 8). However, while chemotaxis resulted in the highest rates of encounter and greatest abundances, it also produced substantially larger seed banks than models that incorporated the less energetically expensive “run and tumble” strategy (Figure 7). This result suggests that while chemotaxis may be effective for locating resources, its energetic demands can cause organisms to enter dormancy as a consequence of exhausting their cell quotas in pursuit of resources. We found that communities went extinct less often and fluctuated less wildly when there was a high capacity for dormancy, i.e., low maintenance energy and low probability of randomly resuscitating (Figure 8). Under a weak capacity for dormancy, communities of active and passive dispersers often went extinct even when the bulk concentration of resources was high and relatively stable (Figure 8). These findings not only confirm our support the importance of seed banks to microbial communities in a variety of systems, but also indicate that extinction can occur under relatively high resource concentrations, i.e., if individuals fail to encounter and consume resources.
Figure 8. Time series results from models with a strong dormancy response (dormancy decreases maintenance energy by a factor of 100, probability of random resuscitation = 0.001) and a weak dormancy response (dormancy decreases maintenance energy by a factor of 10, probability of random resuscitation = 0.1). Top row: 95% confidence hulls reveal that a strong dormancy response almost always prevents extinction while a weak dormancy response leads to greater activity but not necessarily greater total abundance (N). Bottom row: A random selection of two models. (Left) A strong dormancy response leads to a more stable N, despite high variability in or absence of total resources (R). (Right) A weak dormancy response can lead to extinction, even in the presence of increasing resources.
Conclusion
Our study demonstrates that microbial seed banks can emerge from microscale factors that influence the accessibility of resources imposed by the spatial structure of the environment, energetic trade-offs among modes of dispersal, the structural complexity of resources, and to a lesser degree, cross-feeding and scavenging. Our study also illustrates that macroscale factors like the bulk supply and concentration of resources do not always influence encounter rates nor the emergence of microbial seed banks. These findings help explain the presence of large microbial seed banks throughout environmental, engineered, and host-associated ecosystems of high resource concentrations and supply. Likewise, our findings show that while resource limitation may be a primary driver of microbial seed banks, these limitations can drive the emergence of seed banks via multiple mechanisms. These mechanisms include the slow and consistent release of energy from resource that are not highly labile, the sparse spatial distribution of dilute resource particles, the energetic costs of foraging, and the spatial mixing of the environment. To better understand the role of microscale mechanisms in microbial systems and the importance of individual-level interactions and the role of the local environment, we suggest that ecologists pursue microscale work through the intersection of modeling and empirical studies. These efforts may require single-cell metabolic analyses of microscale samples and continued refinement of individual-based models for microbial ecology.
Author Contributions
KL, MF, and JTL designed the study. KL created models and performed analyses. KL, MF, and JTL wrote the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We acknowledge constructive feedback on this research from members of the Lennon Lab at Indiana University. This work was supported by National Science Foundation Dimensions of Biodiversity Grant 1442246 and US Army Research Office Grant W911NF-14-1-0411.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.02040/full#supplementary-material
References
Aanderud, Z. T., Jones, S. E., Fierer, N., and Lennon, J. T. (2015). Resuscitation of the rare boisphere contributes to pulses of ecosystem activity. Front. Microbiol. 6:24. doi: 10.3389/fmicb.2015.00024
Aanderud, Z. T., Vert, J. C., Lennon, J. T., Magnusson, T. W., Breakwell, D. P., and Harker, A. R. (2016). Bacterial dormancy is more prevalent in freshwater than hypersaline lakes. Front. Microbiol. 7:853. doi: 10.3389/fmicb.2016.00853
Alldredge, A. (1998). The carbon, nitrogen and mass content of marine snow as a function of aggregate size. Deep Sea Res. I Oceanogr. Res. Pap. 45, 529–541. doi: 10.1016/S0967-0637(97)00048-4
Allison, S. D., Chacon, S. S., and German, D. P. (2014). Substrate concentration constraints on microbial decomposition. Soil Biol. Biochem. 79, 43–49. doi: 10.1016/j.soilbio.2014.08.021
Allison, S. D., Weintraub, M. N., Gartner, T. B., and Waldrop, M. P. (2011). “Evolutionary-economic principles as regulators of soil enzyme production and ecosystem function,” in Soil Enzymology, eds G. C. Shukla and A. Varma (Berlin: Springer-Verlag), 229–243.
Alvarez, C. R., Alvarez, R., Grigera, S., and Lavado, R. S. (1998). Associations between organic matter fractions and the active soil microbial biomass. Soil Biol. Biochem. 30, 767–773. doi: 10.1016/S0038-0717(97)00168-5
Andersen, K. H., Berge, T., Gonçalves, R. J., Hartvig, M., Heuschele, J., Hylander, S., et al. (2016). Characteristic sizes of life in the oceans, from bacteria to whales. Ann. Rev. Mar. Sci. 8, 217–241. doi: 10.1146/annurev-marine-122414-034144
Arrieta, J. M., Mayol, E., Hansman, R. L., Herndl, G. J., Dittmar, T., and Duarte, C. M. (2015). Dilution limits dissolved organic carbon utilization in the deep ocean. Science 348, 331–333. doi: 10.1126/science.1258955
Azam, F., and Malfatti, F. (2007). Microbial structuring of marine ecosystems. Nat. Rev. Microbiol. 5, 782–791. doi: 10.1038/nrmicro1747
Blagodatskaya, E., and Kuzyakov, Y. (2013). Active microorganisms in soil: critical review of estimation criteria and approaches. Soil Biol. Biochem. 67, 192–211. doi: 10.1016/j.soilbio.2013.08.024
Cordero, O. X., and Datta, M. S. (2016). Microbial interactions and community assembly at microscales. Curr. Opin. Microbiol. 31, 227–234. doi: 10.1016/j.mib.2016.03.015
Datta, M. S., Sliwersdka, E., Gore, J., Polz, M. F., and Cordero, O. X. (2016). Microbial interactions lead to rapid micro-scale successions on model marine particles. Nat. Commun. 7:11965. doi: 10.1038/ncomms11965
DeAngelis, D. L., and Gross, L. J. (1992). Individual-Based Models And Approaches in Ecology: Populations, Communities and Ecosystems. New York, NY: Chapman & Hall.
Don, A., Rodenbeck, C., and Gleixner, G. (2013). Unexpected control of soil carbon turnover by soil carbon concentration. Environ. Chem. Lett. 11, 407–413. doi: 10.1007/s10311-013-0433-3
Dungait, J. A. J., Hopkins, D. W., Gregory, A. S., and Whitmore, A. P. (2012). Soil organic matter turnover is governed by accessibility not recalcitrance. Glob. Chang. Biol. 18, 1781–1796. doi: 10.1111/j.1365-2486.2012.02665.x
Dusenbery, D. B. (2009). Living at the Micro Scale: the Unexpected Physics of Being Small. Cambridge, MA: Harvard University Press.
Fierer, N., and Lennon, J. T. (2011). The generation and maintenance of diversity in microbial communities. Am. J. Bot. 98, 439–448. doi: 10.3732/ajb.1000498
Follows, M. J., and Dutkiewicz, S. (2011). Modeling diverse communities of marine microbes. Ann. Rev. Mar. Sci. 3, 427–451. doi: 10.1146/annurev-marine-120709-142848
Giovannoni, S. J., Tripp, H. J., Givan, S., Podar, M., Vergin, K. L., Baptista, D., et al. (2005). Genome streamlining in a cosmopolitan oceanic bacterium. Science 309, 1242–1245. doi: 10.1126/science.1114057
Grimm, V. (1999). Ten years of individual-based modelling in ecology: what have we learned and what could we learn in the future? Ecol. Model. 115, 129–148. doi: 10.1016/S0304-3800(98)00188-4
Grimm, V., Revilla, E., Berger, U., Jeltsch, F., Mooij, W. M., Railsback, S. F., et al. (2005). Pattern-oriented modeling of agent-based complex systems: lessons from ecology. Science 310, 987–991. doi: 10.1126/science.1116681
Grimm, V., Berger, U., Bastiansen, F., Eliassen, S., Ginot, V., Giske, J., et al. (2006). A standard protocol for describing individual-based and agent-based models. Ecol. Model. 198, 115–126. doi: 10.1016/j.ecolmodel.2006.04.023
Großkopf, T., and Soyer, O. S. (2016). Microbial diversity arising from thermodyamic constraints. ISME J. 10, 2725–2733. doi: 10.1038/ismej.2016.49
Hellweger, F. L., Clegg, R. J., Clark, J. R., Plugge, C. M., and Kreft, J. U. (2016). Advancing microbial sciences by individual-based modelling. Nat. Rev. Microbiol. 14, 461–471. doi: 10.1038/nrmicro.2016.62
Hernández, D. L., and Hobbie, S. E. (2010). The effects of substrate composition, quantity, and diversity on microbial activity. Plant Soil 335, 397–411. doi: 10.1007/s11104-010-0428-9
Hibbing, M. E., Fuqua, C., Parsek, M. R., and Peterson, S. B. (2010). Bacterial competition: surviving and thriving in the microbial jungle. Nat. Rev. Microbiol. 8, 15–25. doi: 10.1038/nrmicro2259
Hoehler, T. M., and Jørgensen, B. B. (2013). Microbial life under extreme energy limitation. Nat. Rev. Microbiol. 11, 83–94. doi: 10.1038/nrmicro2939
Jones, S. E., and Lennon, J. T. (2010). Dormancy contributes to the maintenance of microbial diversity. Proc. Natl. Acad. Sci. U.S.A. 107, 5881–5886. doi: 10.1073/pnas.0912765107
Kiørboe, T., Grossart, H. P., Ploug, H., and Tang, K. (2002). Mechanisms and rates of bacterial colonization of sinking aggregates. Appl. Environ. Microbiol. 68, 3996–4006. doi: 10.1128/AEM.68.8.3996-4006.2002
Lauro, F. M., Mcdougald, D., Thomas, T., Williams, T. J., Egan, S., Rice, S., et al. (2009). The genomic basis of trophic strategy in marine bacteria. Proc. Natl. Acad. Sci. U.S.A. 106, 15527–15533. doi: 10.1073/pnas.0903507106
Lehmann, J., and Kleber, M. (2015). The contentious nature of soil organic matter. Nature 538, 60–68. doi: 10.1038/nature16069
Lennon, J. T. (2007). Diversity and metabolism of marine bacteria cultivated on dissolved DNA. Appl. Environ. Microbiol. 73, 2799–2805. doi: 10.1128/AEM.02674-06
Lennon, J. T., and Jones, S. E. (2011). Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat. Rev. Microbiol. 9, 119–130. doi: 10.1038/nrmicro2504
Lever, M. A., Rogers, K. L., Lloyd, K. G., Overmann, J., Schink, B., Thauer, R. K., et al. (2015). Life under extreme energy limitation: a synthesis of laboratory- and field-based investigations. FEMS Microbiol. Rev. 39, 688–728. doi: 10.1093/femsre/fuv020
Litchman, E., Edwards, K. F., and Klausmeier, C. A. (2015). Microbial resource utilization traits and trade-offs: implications for community structure, functioning, and biogeochemical impacts at present and in the future. Front. Microbiol. 6:254. doi: 10.3389/fmicb.2015.00254
Logue, J. B., Stedmon, C. A., Kellerman, A. M., Nielsen, N. J., Andersson, A. F., Laudon, H., et al. (2016). Experimental insights into the importance of aquatic bacterial community composition to the degradation of dissolved organic matter. ISME J. 10, 533–545. doi: 10.1038/ismej.2015.131
Macalady, J. L., Hamilton, T. L., Grettenberger, C. L., Jones, D. S., Tsao, L. E., and Burgos, W. D. (2013). Energy, ecology and the distribution of microbial life. Philos. Trans. R. Soc. Lond. B Biol. Sci. 368:20120383. doi: 10.1098/rstb.2012.0383
Marsh, D. H., and Odum, W. E. (1979). Effect of suspension and sedimentation on the amount of microbial colonization of salt marsh microdetritus. Estuaries 2, 184–188. doi: 10.2307/1351733
Mitchell, J. G., and Kogure, K. (2006). Bacterial motility: links to the environment and a driving force for microbial physics. FEMS Microbiol. Ecol. 55, 3–16. doi: 10.1111/j.1574-6941.2005.00003.x
Muscarella, M. E., Bird, K. C., Larsen, M. L., Placella, S. A., and Lennon, J. T. (2014). Phosphorus resource heterogeneity in microbial food webs. Aquat. Microb. Ecol. 73, 259–272. doi: 10.3354/ame01722
Moore, C. M., Mills, M. M., Arrigo, K. R., Berman-Frank, I., Bopp, L., Boyd, P. W., et al (2013). Processes and patterns of oceanic nutrient limitation. Nat. Geosci. 6, 701–710. doi: 10.1038/ngeo1765
Ottemann, K. M., and Miller, J. F. (1997). Roles for motility in bacterial-host interactions. Mol. Microbiol. 24, 1109–1107. doi: 10.1046/j.1365-2958.1997.4281787.x
Pande, S., Shitut, S., Freund, L., Westermann, M., Bertels, F., Colesie, C., et al. (2015). Metabolic cross-feeding via intercellular nanotubes among bacteria. Nat. Commun. 6:6238. doi: 10.1038/ncomms7238
Puget, P., Chenu, C., and Balesdent, J. (2000). Dynamics of soil organic matter associated with particle-size fractions of water-stable aggregates. Eur. J. Soil Sci. 51, 595–605. doi: 10.1111/j.1365-2389.2000.00353.x
Rappé, M. S., Connon, S. A., Vergin, K. L., and Giovannoni, S. J. (2002). Cultivation of the ubiquitous SAR11 marine bacterioplankton clade. Nature 418, 630–633. doi: 10.1038/nature00917
Raynaud, X., and Nunan, N. (2014). Spatial ecology of bacteria at the microscale in soil. PLPS ONE 9:e87217. doi: 10.1371/journal.pone.0087217
Rosindell, J., Harmon, L. J., and Etienne, R. S. (2015). Unifying ecology and macroevolution with individual-based theory. Ecol. Lett. 18, 472–482. doi: 10.1111/ele.12430
Rozen, D. E., Philippe, N., Arjan de Visser, J., Lenski, R. E., and Schneider, D. (2009). Death and cannibalism in a seasonal environment facilitate bacterial coexistence. Ecol. Lett. 12, 34–44. doi: 10.1111/j.1461-0248.2008.01257.x
Rusconi, R., Garren, M., and Stocker, R. (2014). Microfluidics expanding the frontiers of microbial ecology. Ann. Rev. Biophys. 6, 65–91. doi: 10.1146/annurev-biophys-051013-022916
Rusconi, R., and Stocker, R. (2015). Microbes in flow. Curr. Opin. Microbiol. 25, 1–8. doi: 10.1016/j.mib.2015.03.003
Schimel, J. P., and Weintraub, M. N. (2003). The implications of exoenzyme activity on microbial carbon and nitrogen limitation in soil: a theoretical model. Soil Biol. Biochem. 35, 549–563. doi: 10.1016/s0038-0717(03)00015-4
Schmidt, M. W., Torn, M. S., Abiven, S., Dittmar, T., Guggenberger, G., Janssens, I. A., et al. (2011). Persistence of soil organic matter as an ecosystem property. Nature 478, 49–56. doi: 10.1038/nature10386
Shoemaker, W. R., Locey, K. J., and Lennon, J. T. (2016). A macroecological theory of microbial biodiversity. Peer J. Preprints 4:e1450ve14. doi: 10.7287/peerj.preprints.1450v4
Smriga, S., Fernandez, V. I., Mitchell, J. G., and Stocker, R. (2016). Chemotaxis toward phytoplankton drives organic matter partitioning among marine bacteria. Proc. Natl. Acad. Sci. U.S.A. 113, 1576–1581. doi: 10.1073/pnas.1512307113
Traving, S. J., Thygesen, U. H., Riemann, L., and Stedmon, C. A. (2015). A model of extracellular enzymes in free-living microbes: which strategy pays off? Appl. Environ. Microbiol. 81, 7385–7393. doi: 10.1128/AEM.02070-15
Trueba, F. J., and Woldringh, C. L. (1980). Changes in cell diameter during the division cycle of Escherichia coli. J. Bacteriol. 142, 869–878.
Vaqué, D., Duarte, C. M., and Marrase, C. (1989). Phytoplankton colonization by bacteria: encounter probability as a limiting factor. Mar. Ecol. Prog. Ser. 54, 137–140. doi: 10.3354/meps054137
Vos, M., Wolf, A. B., Jennings, S. J., and Kowalchuk, G. A. (2013). Micro-scale determinants of bacterial diversity in soil. FEMS Microbiol. Rev. 37, 936–954. doi: 10.1111/1574-6976.12023
Wang, P., Robert, L., Pelletier, J., Dang, W. L., Taddei, F., Wright, A., et al. (2010). Robust growth of Escherichia coli. Curr. Biol. 20, 1099–1103. doi: 10.1016/j.cub.2010.04.045
Wang, G. S., Jagadamma, S., Mayes, M. A., Schadt, C. W., Steinweg, J. M., Gu, L., et al. (2015). Microbial dormancy improves development and experimental validation of ecosystem model. ISME J. 9, 226–237. doi: 10.1038/ismej.2014.120
Wetzel, R. G. (1999). Biodiversity and shifting energetic stability within freshwater ecosystems. Arch. Hydrobiol. Spec. Issues Advanc. Limnol. 54, 19–32.
Keywords: individual based models, microbial diversity, dormany, encounter rate, scaling, energy limitation, deep biosphere, seed bank
Citation: Locey KJ, Fisk MC and Lennon JT (2017) Microscale Insight into Microbial Seed Banks. Front. Microbiol. 7:2040. doi: 10.3389/fmicb.2016.02040
Received: 27 July 2016; Accepted: 05 December 2016;
Published: 09 January 2017.
Edited by:
Mark Alexander Lever, ETH Zurich, SwitzerlandReviewed by:
Jens Kallmeyer, Helmholtz Zentrum Potsdam-GFZ, GermanyJoseph Vallino, Marine Biological Laboratory, USA
Copyright © 2017 Locey, Fisk and Lennon. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: J. T. Lennon, bGVubm9uakBpbmRpYW5hLmVkdQ==