_{1}

^{*}

The paper reviews previous publications and reports some comments about a semi empirical model of the growth and decay process of a planktonic microbial culture. After summarizing and reshaping some fundamental mathematical expressions, the paper highlights the reasons for the choice of a suitable time origin that makes the parameters of the model self-consistent. Besides the potential applications to predictive microbiology studies and to effects of bactericidal drugs, the model allows a suitable proxy of the fitness of the microbial culture, which can be of interest for the studies on the evolution across some thousand generations of a Long Term Evolution Experiment.

Every microbial culture behaves like a factory where substrates become living organisms [

In such perspective, the microbial growth looks like an auto-catalyzed process. However, this is not what really happens, since the growth rate rises up to a maximum and then declines toward zero, in spite of the largely increased number of the cells.

Such a behavior is reflected by the so-called growth curve, log (N)-vs-t, that describes the trend of the population density, N, as function of the time, t.

Since the first decades of the last century, many authors, including the Nobel laureate Jaques Monod [

Recent investigations addressed the attention to the underlying biochemical mechanism and to the molecular peculiarities that characterize the growth of various microbial species in different mediums [

Finally, the awareness of possible genomic mutations that can either enhance or depress the fitness of a given microbial culture, led to a vast specialized literature on the evolution of the microbes and the related effects, like the emergence of resistance to antibiotic drugs, or the improved fitness after thousand generations in Long Term Evolution Experiments (LTEE) [

These approaches could mislead to the conclusion that the early models are just naïve descriptions of the microbial growth. However, in spite of the wide variability of biochemical cycles within the cells and chemical composition of the surrounding medium, all the microbial cultures show the same kind of “macroscopic” behavior, namely, the growth curve, which therefore must obey to laws that prevail biological and biochemical differences and still deserves a reliable interpretation [

The present paper is a review of previous works [

The present paper summarizes the fundamental aspects of the model, the details being reported elsewhere [_{o} cells. No cell dies during the growth progress. The basic assumption of the model is a variable generation time, τ, that allows reproduction of the sigmoid trend of the growth curve as the best fit of the experimental log (N) data that significantly (confidence > 95%) differ from the log (N_{o}) starting population. Such empirical basis therefore reflects a duplication process in action and does not directly concern any earlier adjustment of the cells.

This selection of the experimental data differentiates the present model from all the previous others, which propose best fits that gather data that clearly reflect a no-growth condition together with data that indicate the increase of the cell population. This apparent contradiction comes from the common opinion that the microbes must first “self-adjust” to the medium at the beginning of any plate count experiment [

The present model takes into account such a gap between no-duplication and duplication phases through a wider vision of the issue, stating that the time origin of the whole process (thus including the self-adjustment phase) does not necessarily coincide with the start of the experiment. This statement rises the problem of defining an ideal time origin of the growth. It was noticed that the condition N = 1 is the lowest requirement to obey the constraint that the system has to host a trigger, without which no duplication is possible. For this reason, the “ideal” time origin, θ = 0, for the growth process should comply with such requirement.

In order to single out the “ideal” time origin using the available experimental data, one has to address the attention to the best evidence of the duplication process, namely, the fastest duplication rate. This corresponds to the largest specific growth rate, N ˙ / N , namely the flex point of the growth curve. This condition, dubbed “balanced growth” [

This choice creates a self-consistency between the time scale and the progress of the cell duplication. Now the experimental data to treat must be those that actually reflect the increase of N, neglecting those that do not significantly differ from the starting value, N_{o}. These data enter the best-fit routine based on the equation for the cell duplication,

N = N o 2 ( θ − θ o ) / τ ( θ ) (1)

where τ is the variable generation time and θ_{o} is the time lapse that precedes the onset of the duplication process. The experimental evidence indicates that the duplication rate, 1/τ, is null for θ = θ_{o} and for θ → ∞.

A simple function that obeys this condition [

τ = α ϑ − θ o + ϑ − θ o β (2)

Equation (2) allows rewriting Equation (1) in the log_{2} scale,

log 2 ( N N o ) = β ( ϑ − θ o ) 2 α β + ( ϑ − θ o ) 2 (3)

Equation (3) states that the largest extent of the cell duplication (for θ → ∞) is N_{max} = N_{o}2^{β}, which identifies β as the number of duplication steps experienced by each of the N_{o} duplication lines. A straightforward algebra [

( θ ∗ − θ o ) = ( α β / 3 ) 1 / 2 , when log 2 ( N N o ) ∗ = β 4 . The value of μ is

μ = 3 3 8 β α (4)

This means that the straight line tangent to the growth curve in the plot log_{2} (N)-vs-θ corresponds to the equation

log 2 ( N ) = ( log e ( 2 ) 3 3 8 β α ) θ = log e ( 2 ) μ θ (5)

This straight line goes through log 2 ( N ) = [ log 2 ( N o ) − β / 8 ] for θ = θ_{o}, while crosses the log_{2} (N_{o}) and the log_{2} (N_{max}) levels at θ (0) and θ_{end}, respectively (

For θ < 0, any cellular activity does not directly aim at the duplication. This is another main difference from all the other models reported in the literature.

It is important to remind that the above parameters are properties of the whole system (cells + medium). This means that microbial cultures prepared pouring a given microbial population in different mediums will show different values of the parameters. All the above relationships hold for the growth curves of every prokaryotic microbial culture so far checked by the author, including those quite far from planktonic conditions, once the time scale of the experiment is reported to their own time origin θ = 0.

When so, one can use the reduced variables θ R = ( θ − θ o ) / ( θ ∗ − θ o ) and ξ = log 2 ( N / N o ) / β to gather all the growth curves in a single ξ-vs-θ_{R} master plot [

ξ = θ R 2 3 + θ R 2 (6)

Such master plot reflects the collective self-consistent behavior underlying the shape of the growth curve of every microbial culture, in spite of the physical, chemical, biochemical and biological peculiarities of the system.

Another important relationship easily achievable through a straightforward algebra is:

μ θ o = log 2 ( N o ) − β / 8 log e ( 2 ) (7)

that interconnects three fundamental parameters of the model, namely, μ, β and θ_{o}. In particular, Equation (7) states that, for given N_{o} and β, the shorter the

latency phase that precedes the duplication onset, the faster the specific duplication rate. This makes sense as long as θ_{o} indicates the promptness of the cells to duplicate in the surrounding medium. Finally, since (μθ_{o}) > 0, β < 8log_{2} (N_{o}).

Finally, it can be of some interest to notice that the formal treatment used to describe the duplication process that is typical of prokaryotic microbes holds also for eukaryotic yeasts and molds that show a sigmoid growth curve with a flex point. What one needs to do is just to put

log e ( N ) = log e ( N o ) + [ ( θ − θ o ) / t ] log e ( n ) (8)

with n > 1, and replace β and log_{2} with νβ, where ν = log_{e} (n), and log_{e}, respectively, in the above equations.

The experimental evidence shows that most microbial cultures undergo a decline once they have attained a maximum level of the population density. Between ascending and descending trends of N, a steady intermediate phase or a broad maximum can occur. Literature reports either kind of trends. Since the relationship between growth and decay could be of some relevance in food predictive microbiology as well as in the antibiotic pharmacology, an overall description seemed of some interest.

The above growth model finds a natural extension to the decay phase through the assumption of a basic principle. In the absence of any external adverse agent, death primarily hits the oldest cells [

The N_{o} starting cell of the model planktonic culture have the same age, but, after few duplication steps, the population hosts cells of different age: the newest born represent 50% of the whole population (no cell death occurs during the growth process) and soon overwhelm the starting N_{o} ancestors.

For this reason, if one accepts the above principle for the cell death, the first to dye will be just the N_{o} ancestors, namely a negligible fraction of the population. Small fractions will follow, until death hits the younger generations that correspond to major fractions of N_{max}.

Once reported in a log scale, the start of the decay is almost undetectable and remains so for some while, giving the impression of a culture in a steady condition if the average life span of the cells is large compared to the duplication time, τ. Conversely, if the life span is comparable or lower than τ, the observed trend of N shows a broad maximum between ascending and descending branches. This approach justifies the “cascade” effect of the decay trend that shows an increasingly downward slope [

Formally, one can describe such a behavior through a pseudo exponential function of the time, t, elapsed after the attainment of the N_{max} threshold,

N surv = N max exp ( − t 2 d ) (9)

where d is a constant.

To be reconciled with the ascending trend of N, Equation (9) requires an adjustment of the time scale to align it with the selection of the time origin θ = 0 described above.

The real start of the decay, θ_{s}, remains rather uncertain and therefore becomes a further parameter to determine through a best-fit routine based on the expression for the surviving population,

N surv = N θ = θ s 2 θ − θ s τ d (10)

where τ_{d} is the decay pace. Just like the generation time, τ, τ_{d} depends on the conditions of the medium that likely worsen for increasing (θ − θ_{s}). Putting τ d = d / ( θ − θ s ) , one finally obtains the expression

N surv = N θ = θ s 2 ( θ − θ s ) 2 d (11)

Equation (11) accounts for the above picture of the decay process through a ranking of the cell ages that parallels, in the reverse direction, the growth process [

Finally, if one identifies N ( θ − θ s ) = N max , the number of viable cells throughout the decay phase is

log 2 ( N surv ) = log 2 ( N o ) + β − ( θ − θ s ) 2 d (12)

To define an expression for the overall (growth + decay) profile of N, one can use an expression like

log 2 ( N ) = log 2 ( N o ) + β ( θ − θ o ) 2 α β + ( θ − θ o ) 2 − ( θ − θ o ) 2 δ (13)

where δ ≠ d is an adjustable parameter. If δ ≫ ( θ s − θ o ) 2 , Equation (13) leads to practically the same value as Equation (3), for θ = θ_{s}, namely, log_{2} (N_{max}). Otherwise, the curve bends downward after going through a broad maximum below the level [log_{2} (N_{o}) + β].

When facing adverse conditions, some microorganisms are able to modify their own phenotype or physiological behavior [

In other cases, the changes are very subtle, requiring specific approaches to perceive their effects. These approaches deal with observations of the evolution of a given microbial culture in the course of several thousand generations, as in the Lenski’s LTEE [

As long as β is in relation with μ and θ_{o} (Equation (7)), which are both representative of the efficiency of the cell in a given surrounding medium, β can be an effective proxy of the fitness. The increase of fitness after some thousand generations in a LTEE [_{o} level in the same medium) implies increase of μ and decrease of θ_{o} (

In spite of its semi-empirical nature and the assumed idealized conditions, the growth and decay model summarized in the present paper seems rather flexible and adequate to describe the behavior of real microbial cultures. It suggests a vision of the evolution of a given culture, in a given medium, that relates the duplication activity with the preceding no-growth phase, thanks to a suitable choice of a virtual time origin. The phenomenological description of the growth curve, through the assumption of a variable generation time, reveals some important interconnections between the parameters of the model (all determined with best-fit treatments of experimental data). These suggest an overall vision of the growth process as the result of a cooperative behavior of the microbes since the no-growth phase that precedes the duplication onset. Ranking the population fractions according to the respective age allows a naïve, but reasonable, interpretation of the steady plateau, or broad maximum, between rising and declining branches of the population density. The satisfactory check of the model once applied to a number of experimental data leaves me confident about its reliability.

The proposed approach and related model are easy to use thanks to the simple mathematical treatment of the experimental data and have a number of possible applications to microbial spoilage of food, pharmaceutical investigations about the efficacy of bactericidal and bacteriostatic drugs, as well as to studies on the evolution of microbial organisms in chemostat cultures [

The author did not receive funds to sustain the work and the publication of the paper.

The author declares no conflicts of interest regarding the publication of this paper.

Schiraldi, A. (2024) Microbial Growth and Decay: A Commented Review of the Model. Advances in Microbiology, 14, 1-10. https://doi.org/10.4236/aim.2024.141001