Next Article in Journal
Effects of Extruded Linseed and Soybean Dietary Supplementation on Lactation Performance, First-Service Conception Rate, and Mastitis Incidence in Holstein Dairy Cows
Next Article in Special Issue
Live Feeds Used in the Larval Culture of Red Cusk Eel, Genypterus chilensis, Carry High Levels of Antimicrobial-Resistant Bacteria and Antibiotic-Resistance Genes (ARGs)
Previous Article in Journal
Genome-Wide Association Study Reveals Candidate Genes for Litter Size Traits in Pelibuey Sheep
Previous Article in Special Issue
Fatal Infection in a Wild Sandbar Shark (Carcharhinus plumbeus), Caused by Streptococcus agalactiae, Type Ia-ST7
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Insights into Modelling Bacterial Growth with Reference to the Fish Pathogen Flavobacterium psychrophilum

by
Christopher D. Powell
1,*,
Secundino López
2,3,* and
James France
1
1
Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
2
Departamento de Producción Animal, Universidad de León, E-24007 León, Spain
3
Instituto de Ganadería de Montaña, CSIC-Universidad de León, Finca Marzanas s/n, 24346 Grulleros, Spain
*
Authors to whom correspondence should be addressed.
Animals 2020, 10(3), 435; https://doi.org/10.3390/ani10030435
Submission received: 30 December 2019 / Revised: 26 February 2020 / Accepted: 28 February 2020 / Published: 5 March 2020
(This article belongs to the Special Issue Bacteria-Related Diseases in Fish Species)

Abstract

:

Simple Summary

Flavobacteriumpsychrophilum is a cold-water bacterium responsible for cold water disease and rainbow trout fry syndrome which has significant impacts on fish health and, by extension, negative economic impacts on aquaculture operations. Models can be applied to bacterial growth curves yielding parameter estimates describing rates of bacterial growth and the time it takes for a bacterium to reach its exponential phase of growth (lag time). These parameter estimates can be used to establish the relationship between microbial growth and environmental variables such as pH, temperature and effect of anti-microbial treatments. Two novel models are derived and their potential to describe bacterial growth assessed through their ability to mimic the growth of Flavobacterium psychrophilum on liquid media. Due to their mechanistic derivation, the proposed models result in flexible and robust growth functions that can be expressed as equations with biologically meaningful parameters. Based upon statistical measures of goodness-of-fit and cross-validation, the purposed models were able to describe satisfactorily the growth of Flavobacterium psychrophilum on various media. Furthermore, the proposed models also provide insight into underlying mechanisms that are driving microbial growth and how the current environment affects bacterial rate of growth.

Abstract

Two new models, based upon the principles promulgated by Baranyi and co-workers are presented and resulting growth functions evaluated based upon their ability to mimic bacterial growth of the fish pathogen Flavobacterium psychrophilum. These growth functions make use of a dampening function to suppress potential growth, represented by a logistic, and are derived from rate:state differential equations. Dampening effects are represented by a rectangular hyperbola or a simple exponential, incorporated into a logistic differential equation and solved analytically resulting in two newly derived growth equations, viz. logistic × hyperbola (log × hyp) and logistic × exponential (log × exp). These characteristics result in flexible and robust growth functions that can be expressed as equations with biologically meaningful parameters. The newly derived functions (log × hyp and log × exp), along with the Baranyi (BAR), simple logistic (LOG) and its modified form (MLOG) were evaluated based upon examination of residuals and measures of goodness-of-fit and cross-validation. Using these criteria, log × hyp, log × exp and BAR performed better than, or at least equally well as, LOG and MLOG. In contrast with log × exp and BAR, log × hyp can be easily manipulated mathematically allowing for simple algebraic expressions for time and microbial biomass at inflexion point, in addition to maximum and scaled maximum growth rates.

1. Introduction

With a move towards more exact methods in microbiology, the demand for less empirical models increases [1]. Although many equations are able to describe bacterial growth, to be considered a mathematical model and not just a convenient relationship for empirically fitting to data, an equation should have a sound physiological basis underlying the relationship [1,2,3]. The use of growth functions in the field of microbiology is widespread. In predictive food microbiology, these functions are used in so-called primary modelling to yield useful parameter estimates such as relative growth rate and lag time. These growth features are often then used for further model development (viz. secondary or tertiary modelling) to establish the relationship between microbial growth and environmental variables (e.g., pH, temperature). These estimates of growth parameters enable the prediction of the time required for pathogenic or foodborne microorganisms to reach a critical limit, resulting in clinical signs of disease or food spoilage [4]. Growth functions also find applications in the animal health field. With respect to the field of aquaculture, growth of the pathogenic bacteria Aeromonas hydrophila and Vibrio alginolyticus has been described using modified versions of the logistic and Gompertz equations [5].
Traditionally, microbial growth measurements have been made using the viable counts method. However, estimating microbial growth parameters from fitting growth functions to optical density measurements has become a more common practice, with the advantages of being rapid and inexpensive compared to the viable counts method [4]. Optical density measurements have been used extensively to describe bacterial growth in various fields of microbiology, in addition to allowing for accurate comparison of model-derived parameter estimates [4,6,7,8,9,10]. Through the correlation between the intensity of light absorbed/reflected and the number of cells in the culture medium, optical density measurements provide an estimate of microbial mass [11]. Regardless of the measure (viable count vs. optical density), microbial growth can be characterized by three distinct phases, viz. the lag, ‘exponential’ and stationary phases. The growth rate of the microbial population increases rapidly until a maximum is reached. The intercept between initial microbial biomass and the tangent to the growth curve at the time where growth rate reaches a maximum is generally referred to as the lag phase. Following the so-called exponential phase, growth rate decreases steadily towards zero reaching a stationary phase, resulting overall in a sigmoidal shape [12].
Microbial growth is commonly described using a number of sigmoidal functions, e.g., the logistic, Gompertz and the purpose-built Baranyi and Huang equations [2,13,14,15,16,17,18,19]. The logistic and the Gompertz equations were initially derived to describe population growth and human mortality, respectively [20,21]. Zwietering et al. [12] modified these functions, plotting the logarithm of the relative population size [ln(x/x0)] against time and attempting to ascribe increased biological meaning to some of the parameters. However, due to these modifications the resulting growth functions are not simple transformations of the original mechanistically-derived Gompertz and logistic equations and are better referred to as modified Gompertz and logistic equations [22]. As a result, mechanistic interpretation is no longer straightforward, and these modified models tend to give empirical descriptions of the sigmoidal pattern of microbial growth. Using a more detailed mechanistic approach, Baranyi and Roberts developed a six-parameter model [23], again plotting log-transformed data [ln(x)] against time. However, fitting this particular model can be challenging when using standard nonlinear regression programs [24,25]. Following re-parameterization and incorporation of a number of additional assumptions, a simplified four-parameter model was proposed which is the most widely applied of the Baranyi equations [2,7,15,24].
Flavobacterium psychrophilum is a gram-negative yellow-pigmented bacterium and the etiological agent of both cold-water disease (CWD) and rainbow trout fry syndrome (RTFS) [26,27]. F. psychrophilum can be found across a wide geographic range covering North and South America, Europe and parts of Asia [28]. Although found primarily in salmonid species, F. psychrophilum has been isolated in non-salmonids that inhabit water temperatures below 16 °C and affects both free-ranging and aquacultured fish [29]. Both CWD and RTFS significantly contribute to economic loss, with mortalities of fish affected by RTFS reaching up to 50% while CWD, a more chronic disease, results in 2–30% mortality in larger fish [28,30]. F. psychrophilum not only has ramifications in the aquaculture industry but also affects hatcheries used for stocking and rehabilitation purposes. To the authors’ knowledge, standard growth functions have yet to be applied to F. psychrophilum. The application of such functions can provide insight into the extent to which F. psychrophilum reacts to changes to its in vitro environment, including the effects of various antibiotics, through examination of the growth parameters in a given environment. This is especially relevant as antibiotic challenge tests are common practice in aquaculture, which is particularly important given the increasing resistance of F. psychrophilum to approved antibiotic drugs [31]. Given the limitations of standard functions modified to describe microbial growth and the challenges associated with the Baranyi equations, two novel models following an approach to derivation akin to that from which the Baranyi models were developed are presented and the resulting mathematical forms applied to data on F. psychrophilum.

2. Materials and Methods

2.1. Datasets

Eight datasets pertaining to F. psychrophilum grown on eight different liquid media, taken from two studies, were used for model evaluation. In Study 1, the first four datasets encompass growth of Strain SH3-81 of F. psychrophilum isolated from the kidney of juvenile coho salmon exhibiting classical bacterial cold-water disease (Sandy Hatchery, Oregon, USA). F. psychrophilum was plated on (1) tryptone yeast extract salts (TYES), (2) Shieh, (3) modified Cytophaga and (4) Cy7 broths [29,32,33,34]. The bacterium was inoculated in each of the four broth media and incubated for 24 h at 17 °C on a shaker prior to the start of the growth study. In duplicate, 1 mL of each broth culture was added to 50 mL of the corresponding broth medium using Nephelo culture flasks (Belco Glass, Inc., Vineland, NJ). Using a gyratory shaker (624 Environmental incubator shaker, New Brunswick Scientific Co., Inc., Edison, NJ), flasks were incubated at 17 °C and 1.136 g. Optical density was measured at 525 nm wavelength and recorded every 3–4 h for 72 h with datasets comprising 16 to 20 data points. See Holt [29] for further details concerning bacterial growth conditions.
In Study 2, four datasets encompassing growth of F. psychrophilum strain Trouw 224 were isolated from diseased rainbow trout in Spain. Bacteria were pre-cultured in 25 mL of four liquid media, namely (5) tryptone–yeast extract–salts broth supplemented with glucose (FLPB), (6) TYES broth, (7) Cytophaga broth supplemented with carbohydrates and skimmed milk (CBCM), and (8) modified Anacker and Ordal broth (MAOB). These starting cultures were then inoculated (2 mL) into 200 mL of the respective liquid medium and incubated with 100 rpm shaking at 18 °C. Optical density was measured at 580 nm for 40 h with measurements being taken every 5 h. See Cepeda et al. [35] for further details concerning bacterial growth conditions.
The entirety of the observed growth values (Datasets 1–8) used in this study can be found in the Supplementary Materials Table S1.

2.2. Mathematical Considerations

2.2.1. Potential Growth

Following Baranyi et al. [23], we assume the potential rate of bacterial growth follows a logistic pattern. We formalize our assumption thus:
d x d t = μ x ,   0 t <
μ = μ ( 1 x x m ) ,   0 < x x m ,   x ( 0 ) = x 0 ,   μ > 0
where x denotes microbial biomass at time t, x 0 is microbial biomass at time zero, x m is maximal microbial biomass as t   , μ is a constant representing relative growth rate and μ relates rate of growth ( d x / d t ) to a specific microbial biomass ( x ). Substituting for μ’ in Equation (1) using (2) and integrating analytically by the method of partial fractions gives f(t), the potential growth logistic function:
x = x 0 x m x 0 + ( x m x 0 ) e μ t = f ( t )
This equation can be re-parameterized as follows:
x = p 1 1 + p 2 e p 3 t
p 1 = x m ,   p 2 = x m x 0 1 ,     p 3 = μ
Substituting for μ’ in Equation (1) and differentiating yields:
1 μ d 2 x d t 2 = ( 1 2 x x m ) d x d t
Equating Equation (6) to zero gives the point of inflexion (t*, x*) for f(t) at:
x = x * = x m 2
This occurs when:
t = t * = 1 μ ln ( x m x 0 x 0 )

2.2.2. Actual Growth

We determine the actual rate of bacterial growth by dampening the potential rate, i.e.,
d x d t = α μ x , 0 t <
μ = μ ( 1 x x m ) , 0 < x x m , x ( 0 ) = x 0 , μ > 0
α α ( t ) , 0 < α 1 , α ( 0 ) = 1 ,     α ( ) 0
The dampening factor α, which gives the instantaneous ratio of actual to potential growth rate, proposed herein is a monotonically decreasing function of time that mimics, for example, the effect of a build-up of an inhibitor on growth rate.
The lag phase may be identified as follows. The slope (s) of the actual growth curve g(t) at its point of inflexion (t*, x*) is obtained from Equations (9) to (11) as the product:
s = α ( t * ) × μ ( x * ) × x *
s / x * = α ( t * ) × μ ( x * )
Note that this product (Equation (12)) represents maximum growth rate. The equation of the tangent line, y(t), touching the steepest part of the curve (i.e., at the point of inflexion) is therefore:
y x * = s ( t t * )
The length of the lag phase, T, is determined by the point of intersection between this tangent line and the abscissa x = x0, with Equation (14) yielding:
T = t * ( x * x 0 ) s

Rectangular Hyperbola

The dampening effect α can be represented by a rectangular hyperbola:
α = λ λ + t , λ > 0
where λ is the point whereby the dampening effect is half maximal. Equation (16) is akin to the expression from enzyme kinetics for relative velocity in the presence of a noncompetitive inhibitor (see Segel [36]), assuming the instantaneous inhibitor concentration is directly proportional to time, and substituting for t in Equation (16). The actual growth function g(t) is obtained by substituting for μ’ and α in Equation (9) using (10) and (16), then integrating analytically using partial fractions to give:
x = x 0 x m x 0 + ( x m x 0 ) ( λ λ + t ) λ μ g ( t )
Equation (17) can be expressed in the re-parameterized form:
x = p 1 1 + p 2 ( 1 + p 3 t ) p 4
p 1 = x m ,   p 2 = x m x 0 1 ,   p 3 = 1 λ ,     p 4 = λ μ
Differentiating Equation (9) having substituted for μ’ and α yields:
d 2 x d t 2 = μ ( λ λ + t ) ( 1 1 λ μ 2 x x m ) d x d t
Therefore, a point of inflexion would exist at:
x = x * = ( 1 1 λ μ ) x m 2
This would occur when:
t = t * = λ { [ ( λ μ 1 ) ( x m x 0 ) ( λ μ + 1 ) x 0 ] 1 / ( λ μ ) 1 }

Simple Exponential

The dampening effect can also be represented by a simple exponential function:
α = e λ t , λ > 0
where λ is a rate parameter. Equation (23) states the ratio of actual to potential growth rate declines in a negative exponential manner over the growth period, possibly reflecting the deleterious build-up of some inhibitor. Substituting for μ’ and α in Equation (9) using (10) and (23), then integrating analytically using partial fractions gives g(t):
x = x 0 x m x 0 + ( x m x 0 ) exp [ μ λ ( 1 e λ t ) ]   g ( t )
Re-parameterizing Equation (24) gives:
x = p 1 1 + p 2 exp [ p 3 ( 1 e p 4 t ) ]
p 1 = x m ,     p 2 = x m x 0 1 ,       p 3 = μ λ ,     p 4 = λ    
Differentiating Equation (9) having substituted for μ’ and α yields:
d 2 x d t 2 = [ μ e λ t ( 1 2 x x m ) λ ] d x d t = 0
Therefore, the point of inflexion (t*, x*) satisfies the pair of equations:
x * = x m 2 ( 1 λ μ e λ t * )
t * = 1 λ ln { 1 + λ μ ln [ x 0 ( x m x * ) x * ( x m x 0 ) ] }
This pair of equations has to be solved using numerical methods [37].
The newly derived functions, herein referred to as logistic × hyperbola (log × hyp) and logistic × exponential (log × exp), were investigated along with popular growth functions (Table 1). The simple logistic was modified according to Zwietering et al. [12], where xm and T are as described above. In the Baranyi model, parameter µmax is defined as the theoretical maximum relative rate, with µmax also representing maximum relative rate in the modified logistic.

2.3. Model Fitting

For each of the eight datasets, microbial growth data were expressed on an absolute (untransformed) basis (optical density, OD) in addition to being expressed on a logarithmic basis [ln(OD/OD0)]. The BAR and MLOG models were derived and modified, respectively, to describe logarithmic transformed microbial growth data and therefore these models were fitted to the logarithmic transformed data. In contrast, the LOG, log × hyp and log × exp models were derived to describe absolute growth data and thus fitted to the untransformed data. All models were fitted by means of nonlinear regression using the Newton algorithm in the NLIN procedure of SAS [38]. A range of starting values for each parameter was determined through visual inspection of the growth curve. For parameter values that were difficult to estimate using visual inspection, xm and x0 were fixed in the first attempt at model fitting, providing an estimated starting value for the remaining parameters for future iterations whereby all parameters of the model were estimated. Using the range of starting values, PROC NLIN forms a grid and evaluates the model at each point on the grid. The values on the grid that yield the smallest objective function are used as initial parameter estimates for the first iteration of the fitting processes [38]. The range of starting values was the same between models that contained common parameters.

2.4. Statistical Analysis

Models were evaluated for goodness-of-fit along with analysis of their residuals using a range of statistical tests. Mean square prediction error (MSPE) is a measure of the expected mean squared difference between predicted and observed values. MSPE was calculated as the sum of the squared difference between predicted and observed values divided by the number of observations [39]. Agreement between model predictions and observed values was further evaluated using the concordance correlation coefficient (CCC), a single statistic which contains both accuracy and precision indicators [40]. CCC values were calculated according to Lin [40] with values of 1 representing perfect agreement, −1 perfect disagreement, and 0 no agreement between predicted and observed values [40]. Models’ goodness-of-fit for individual datasets were compared to one another using the Akaike information criterion (AIC). AIC is a test for model selection which accounts for goodness-of-fit while penalizing for over-fitting [41]. AIC was calculated for each model for a given dataset according to Akaike [41] with the model resulting in the smallest AIC being considered the best [41]. The accuracy factor (AF) averages the minimum ‘distance’ between each predicted data point and the line of equivalence. Thus, the accuracy factor is a measure of the average deviation of a model’s predictions and is used as a simple measure of the level of confidence in these predictions [42]. AF was calculated according to Baranyi et al. [43]. Calculated AF will always be greater than or equal to one, with larger values resulting from poorer fits.
The ability of each model to predict growth without systematically over- or under-estimating microbial mass was determined using the runs and Durbin–Watson (DW) tests. The number of runs test examines a sequence of residuals for unusual groupings of positive or negative residuals and tests against the null hypothesis that arrangement of signs is random. A run is a sequence of residuals with the same sign, positive or negative. Using the runs test, the probability of too few runs (indicating clustering of residuals with the same sign leading to systematic bias) or too many runs (indicating negative serial correlation) was determined [44]. Unlike the runs test, which ignores the actual size of the residual, the DW test examines dependencies in the error terms by testing for correlations between a residual and the residuals immediately before and after it in the sequence. The DW statistic (d), and upper (du) and lower (dl) critical values were calculated according to Draper and Smith [44]. When d is less than the lower critical value dl, evidence of positive autocorrelation occurs, when the d value is greater than the upper critical value du, evidence of negative autocorrelation occurs. The test is inconclusive when d falls between the upper and lower critical values. Positive serial autocorrelation in the residuals suggests the model has the tendency to systemically over- or under-estimate projected values [44]. The quantitative bias factor (BF) was calculated according to Ross [42]. Perfect agreement between predictions and observations will result in a BF of one. Values of BF greater than one occur when a model’s predictions are on average greater than observed values, while a value of less than unity occurs when a models’ predictions are on average less than observed values. When calculating the AF and BF statistics, both observed and predicted values underwent a ’x + 1’ transformation. This transformation allowed for the inclusion of data points whereby x = 0 when calculating these statistics. Likewise, data points that approached zero, and due to the logarithmic nature of the calculation result in nonsensical AF and BF values, were able to be included in the calculations following ‘x + 1’ transformation.

2.5. Model Validation

As a means of cross-validation of proposed models, the statistic predicted residual error sum of squares (PRESS) was calculated. PRESS is based upon the leave-one-out or jackknife technique which evaluates the fit of a model to a sample of observations which in their entirety were not used for model development [45,46,47]. The basis of the statistic is the summing of the squared residuals resulting from systematically omitting individual observations, one by one, from a given dataset, while simultaneously refitting the proposed model to these reduced datasets. For example, if a dataset contained n data points, n datasets would be created from this one dataset, all differing by exactly one data point, with the model being fitted to these n datasets. After each individual fitting, the value of the y-variable is predicted at the omitted data point. The PRESS statistic was calculated according to Allen [45] providing an unbiased estimate of the model’s future prediction performance [47]. When fitting multiple models to the same dataset, the lowest PRESS value indicates the superior model. Additionally, CCC was used as a measure of goodness-of-fit between model predictions when fitted to datasets encompassing omitted data points and the corresponding observed values.

3. Results

The suitability of the two new models derived (viz. log × hyp and log × exp) was assessed through their ability to mimic the growth of Flavobacterium psychrophilum cultivated on eight liquid media based upon absolute OD growth data. In addition, the simple logistic was fitted to these OD data, and its modified form plus the four-parameter Baranyi were fitted to the logarithmic transformation [ln(OD/OD0)] of the datasets as a means of affirmation. All five models examined in this study were evaluated based upon fitting behavior, examination of residuals, measures of goodness-of-fit and cross-validation.

3.1. Paramater Estimates

In general, all five growth functions fitted the data without major problems. Through visual inspection of all datasets, a range of initial parameter values was obtained for use with the grid search technique in the statistical software SAS. The resulting final values of lag time (T), maximum growth rate (s) and scaled maximum growth rate (s/x*) for the two new derivations in addition to the simple logistic (along with µmax and T for BAR and MLOG) are presented in Table 2. Unlike MLOG and BAR, lag times for log × hyp and log × exp are not explicitly represented as a parameter in the equation and must be calculated using Equation (15). Likewise, maximum growth rate for log × hyp and log × exp is not explicitly represented as an equation parameter and therefore was calculated using Equation (12), while scaled maximum growth rate was calculated by dividing maximum growth rate by microbial biomass at the point of inflexion (Equation (13)). Estimates of s, s/x* and T were influenced by dataset. In Study 1, using OD data, the LOG, log × hyp and log × exp estimates of s, s/x* and T were in close agreement with one another. These three models determined the longest lag time and slowest maximum growth rate to occur when F. psychrophilum was cultured on the Cy7 medium (Table 2—Dataset 4). All three models were in agreement that fastest maximum growth rate and highest scaled maximum growth rate occurred when grown on the modified Cytophaga medium (Table 2—Dataset 3), with the shortest lag time and lowest scaled maximum growth occurring with the Shieh medium (Table 2—Dataset 2). When applied to logarithmic transformed OD data in Study 1, close agreement was found between the µmax and T of the BAR and MLOG. In agreement with the log × hyp, log × exp and LOG, the models BAR and MLOG determined the longest lag time when F. psychrophilum was cultured on Cy7 (Table 2—Dataset 4). Both BAR and MLOG were in agreement that the shortest lag time occurred with the TYES medium (Table 2—Dataset 1). Highest relative growth rate with BAR and MLOG occurred on Shieh (Table 2—Dataset 2), with lowest relative growth rate occurring on Cy7 (Table 2—Dataset 4). It is important to note that µmax of the Baranyi represents a maximum relative rate (units of per unit time, h−1), whereas the s of log × hyp and log × exp is a maximum rate (units of microbial biomass or its surrogate per unit time) and s/x* is a scaled maximum rate (units of per time, h−1).
In Study 2, LOG, log × hyp and log × exp were in agreement that the lowest maximum growth rate and shortest lag time occurred with the MAOB medium (Table 2—Dataset 8), with both log × hyp and log × exp determining lowest scaled maximum growth rate to occur on this medium. These three models determined maximum growth rate (absolute or scaled) was highest on the TYESB medium (Table 2—Dataset 6) and lag time longest on the CBCM medium (Table 2—Dataset 7). With both MLOG and BAR maximum relative growth rate occurred on CBCM (Table 2—Dataset 7), though the two values (0.847 and 1.66, respectively) were significantly different from one another. Also, with both models the lowest maximum relative growth rate occurred on the FLPB medium (Table 2—Dataset 5). Both models were in agreement that longest lag time occurred on MAOB (Table 2—Dataset 8) and shortest on FLPB (Table 2—Dataset 5).
Refer to Tables S2 and S3 in the Supplementary Materials for all of parameter estimates generated using log × hyp and log × exp when fitted to the eight liquid media datasets.

3.2. Growth Prediction

Following parameter estimation, growth of F. psychrophilum was predicted as a function of time using the five growth equations, viz. Equations (3), (17) and (24) and those for MLOG and BAR displayed in Table 1. An example of observed and predicted absorbance values given by each growth function is presented in Figure 1.
As seen in this figure, all five models appear to describe the observed growth data well with apparently little distinction between models. This figure also illustrates the effect that a logarithmic transformation has on the shape of the growth curve. The dampening effect of the log × hyp and log × exp functions, along with predicted growth, is shown in Figure 2.
In this figure, potential growth rate and potential growth, represented by Equations (1) and (3) respectively, are presented for log × hyp and log × exp. Using either a rectangular hyperbola (Equation (16)) or a simple exponential (Equation (23)), potential growth rate is dampened resulting in actual growth rate (Equation (9)) for log × hyp and log × exp. Integrating Equation (9), the actual rate of growth, with the associated dampening functions (Equations (16) and (23)) results in Equations (17) and (24), representing OD as a function of time for log × hyp and log × exp, respectively. Refer to Figures S1–S6 in the Supplementary Materials for the figures corresponding to Figure 1 for the remaining liquid media. Comparison of growth and associated growth rates of F. psychrophilum between the two liquid media that resulted in the highest and lowest maximum growth rates (s) in Studies 1 and 2 obtained using Equation (12) for log × hyp and log × exp is presented in Figure 3. The interaction between growth rate and resulting growth can clearly be observed, in addition to the time at which s occurs coinciding with the point of inflexion of the respective growth curve.

3.3. Model Evaluation

Models were evaluated for goodness-of-fit using four criteria: Akaike information criterion, mean square prediction error, concordance correlation coefficient and accuracy factor. The goodness-of-fit values for the five models were compared over studies in addition to being averaged over the eight datasets.
The averaged AIC, MSPE, CCC and AF values resulting from fitting the five models to the four datasets of Study 1 and the four datasets of Study 2, in addition to the totaled average over the eight datasets, are given in Table 3. The BAR resulted in a smaller averaged AIC (−55.7 ± 5.6) value in comparison with MLOG (−51.6 ± 3.7), when comparing across all eight datasets. In Study 1, BAR (−67.5 ± 5.5) resulted in a smaller averaged AIC value compared with MLOG (−58.8 ± 3.1), while very similar AIC values were observed between BAR (−44.0 ± 2.9) and MLOG (−44.4 ± 3.6) in Study 2. Comparing AIC values resulting from models fitted to the original, rather than the logarithmic transformed data, viz. LOG, log × hyp and log × exp, the LOG (–80.8 ± 15.1) gave the smallest AIC value averaged across all eight datasets, followed by log × exp (–78.3 ± 14.0) and log × hyp (−77.0 ± 14.0). The LOG (−117.1 ± 11.5) resulted in the lowest averaged AIC value in Study 1, followed by log × exp (−111.6 ± 11.3) and log × hyp (−110.0 ± 11.5). In contrast, log × exp (−45.0 ± 2.9) resulted in the lowest AIC value in Study 2 followed by LOG (−44.6 ± 3.9) and log × hyp (–43.9 ± 2.7). Based upon MSPE, the BAR (0.009 ± 0.003) outperformed MLOG (0.013 ± 0.004) in both studies. LOG (0.001 ± 0.0002), log × hyp (0.001 ± 0.0003) and log × exp (0.001 ± 0.0003) performed equally well in Study 1, while in Study 2 both log × exp (0.003 ± 0.001) and log × hyp (0.004 ± 0.001) resulted in lower averaged MSPE compared with LOG (0.005 ± 0.002). Both log × exp (0.002 ± 0.001) and log × hyp (0.002 ± 0.001) resulted in lower averaged MSPE values compared to LOG (0.003 ± 0.001) when averaged over all eight datasets.
Similar to the results found for MSPE, the LOG, log × hyp and log × exp performed equally well (0.998 ± 0.001) on the basis of CCC when fitting the datasets of Study 1. When fitting the datasets of Study 2, both log × hyp (0.988 ± 0.005) and log × exp (0.989 ± 0.005) gave higher CCC values compared to LOG (0.987 ± 0.004), with higher averaged CCC values for both log × hyp (0.993 ± 0.003) and log × exp (0.993 ± 0.003) compared to LOG (0.992 ± 0.003) over the eight datasets. The BAR resulted in higher CCC values in both Study 1 and 2 compared to MLOG with averaged values of 0.998 ± 0.001 and 0.997 ± 0.001, respectively. On the basis of AF, BAR (1.024 ± 0.005) performed better than MLOG (1.039 ± 0.005) in Study 1 in addition to Study 2, BAR (1.016 ± 0.003) MLOG (1.022 ± 0.007). When averaged over the eight studies, BAR (1.020 ± 0.003) resulted in a lower AF compared to MLOG (1.031 ± 0.005). In Study 1, LOG (1.013 ± 0.002) performed slightly better than log × exp (1.015 ± 0.002) and log × hyp (1.016 ± 0.003) on the basis of AF index. In Study 2, log × exp (1.027 ± 0.002) resulted in the smallest AF, followed by log × hyp (1.030 ± 0.003) and LOG (1.033 ± 0.005). When averaged over all eight datasets, log × exp (1.021 ± 0.002) resulted in the smallest AF statistic, followed by log × hyp (1.023 ± 0.003) and LOG (1.023 ± 0.005).
Residuals were examined through clustering of same signed residuals and through serial correlation, with results of the runs test and the DW test for serial correlation presented in Table 4. Distribution of runs of signs was examined according to Draper and Smith [44]. None of the five models showed a significant number of too many runs of signs (data not shown). Of the eight datasets, runs of signs were deemed to be random in all eight datasets for both BAR and MLOG. In contrast, too few runs were observed with LOG, log × hyp and log × exp in 2, 3 and 2 of the datasets, respectively. Runs of signs were found to be random with the remaining datasets. Serial correlation was investigated using the DW statistic and the number of times a specific model gave positive or no serial correlation when fitted to a given dataset is shown in Table 4. Both MLOG and BAR displayed no serial correlation when fitted. The remaining three models showed no serial correlation for seven of the eight datasets, with these models exhibiting positive serial correlation for one of these datasets. Interestingly, evidence of serial correlation in the residuals for LOG, log × hyp and log × exp all occurred with the same dataset, viz. Dataset 4. Comparing MLOG to BAR on the basis of BF index, the two models had a tendency to overestimate observations in Study 1 and in Study 2, with BAR resulting in a BF closer to 1 in both studies. In Study 1, LOG, log × hyp and log × exp had a tendency to slightly underestimate observed values with all three models resulting in BF below 1, 0.996 ± 0.002, 0.993 ± 0.001 and 0.994 ± 0.000, respectively. In Study 2, LOG slightly overestimated observed values (1.001 ± 0.003) while both log × hyp (0.994 ± 0.002) and log × exp (0.994 ± 0.003) underestimated them. When averaged over all eight datasets, LOG resulted in a value just below 1, (0.999 ± 0.002) while log × hyp (0.994 ± 0.001) and log × exp (0.991 ± 0.001) also slightly underestimated observed values on average.

3.4. Model Validation

Using the leave-one-out/jackknife technique, the models were successively fitted to sub-datasets of Dataset 1 and 2 from Study 1 and Dataset 5 and 6 from Study 2. PRESS and CCC statistics were calculated from model predictions and observed values for these sub-datasets (Table 5). When examining the PRESS statistic, the lowest value indicates the best future prediction performance while CCC presents a clear indication of prediction precision and accuracy. The log × exp resulted in the lowest PRESS value in three of the four datasets with LOG resulting in the lowest value in the other. Examining CCC, log × exp resulted in the best CCC value in three of the four datasets with LOG giving the highest in the other. Both log × exp (0.9814, 0.9808) followed by log × hyp (0.9782, 0.9807) resulted in higher CCC values when fitted to the sub-datasets from Dataset 5 and 6 compared to LOG (0.9726, 0.9806).

4. Discussion

Simple mechanistic growth models are derived from rate:state principles represented as a differential equation relating growth rate to organism or population size. Equations expressed in this form allow straightforward biological interpretation of the parameters, and in many cases, can be solved analytically. Following this methodology, classical growth equations have been derived including the logistic, Gompertz, and Richards [20,21,48]. In an effort to apply them to bacterial growth curves, Zwietering et al. [12] modified these functions in such a way as to ascribe increased biological meaning to some of the parameters, resulting in the modified logistic, Gompertz and Richards. The modified logistic has been shown to describe microbial growth with more reliable kinetic parameters estimates compared to the modified Gompertz [6,49]. As such the modified logistic, viz. MLOG, was chosen as a basis for model comparison along with BAR in this study.
The Baranyi family of equations [2,15,23], in addition to the model of Huang [18,19], represent two of very few attempts to derive growth functions mechanistically to describe bacterial growth. The Baranyi models are based upon the concept of potential and actual growth, with actual growth being less than potential as the bacteria adapt from the pre- to post-inoculation environment. Potential growth is represented by the first-order differential equation d x / d t = μ ( x ) x whereby x denotes cell concentration and µ(x) denotes relative growth rate [23]. The differential equation adopted in the Baranyi models is akin to the logistic. However, due to food microbiologists’ preference for using the logarithm of cell concentration, the logarithmic transformation of microbial cell concentration was related to incubation time [2]. This transformation departs from the mechanistic derivation of the logistic and hinders the model’s application to absolute, non-transformed growth data. As actual growth is less than potential, an adjustment function is incorporated into the differential equation representing growth potential. The resulting differential equation for actual bacterial growth therefore takes the form d x / d t = μ max α ( t ) u ( x ) x with α ( t ) representing the adjustment function and u ( x ) the so-called inhibition function. The authors direct readers to Baranyi et al. [23] and Baranyi and Roberts [2] for derivation of this function and all other mathematical considerations. The differential equation representing actual growth can be solved analytically resulting in a six-parameter equation [2]. However, due to fitting difficulties associated with the large number of parameters, and sensitivity to number of data points, certain assumptions were implemented to reduce the original Baranyi to a four-parameter equation [2,24,25]. The resulting four-parameter model BAR is the most widely applied of the Baranyi family of models and one of the models investigated herein [7]. Like the Baranyi family of models, the Huang model is also mechanistically derived to describe bacterial growth. While the Baranyi takes into account the pre- and post-inoculation state of the cells on growth, Huang [18] aimed at developing a model that described bacterial growth on the basis of current growth environment. However, the discontinuous nature of the model in describing the three phases of growth, viz. lag, exponential and stationary, requires the use of a transitional function. In addition to complicating the mathematics in manipulating the analytical solution to the growth function, the transitional function introduces a non-biological parameter into the model.
Two new models, based upon the principles promulgated by Baranyi and co-workers, aimed at describing bacterial growth are presented in this paper. These growth functions, making use of the concepts of potential and actual growth, utilize a mathematical function to dampen potential logistic growth. Dampening effects are represented by either a rectangular hyperbola or a simple exponential. These dampening functions are incorporated into the basic differential equation which is then solved analytically resulting in two growth equations (viz. log × hyp and log × exp). The use of these seemingly empirical dampening functions results in growth equations that can be expressed as relatively simple analytical expressions. The Baranyi makes use of an adjustment function to take into account the effects of pre- and post-inoculation environments on growth. Like the log × hyp and log × exp models, potential growth in the Baranyi is represented by a first-order differential equation. Actual growth is given by incorporating an adjustment function into this first-order differential equation. The adjustment function of Baranyi and Roberts [2] is based upon the principle that growth during the lag phase is inhibited by a ‘bottle-neck’ intracellular substance following Michaelis-Menten kinetics [2]. The adjustment function takes the form: α ( t ) = q 0 q 0 + e v t , a two-parameter function with a rate parameter (v) and a constant ( q 0 ) representing the physiological state of the inoculum at time zero. However, as the adjustment function gives no indication of what the critical factor causing the ‘bottle-neck’ could be, the adjustment function of the Baranyi, like dampening functions of log × hyp and log × exp, is essentially a quasi-mechanistic expression with biological meaning being prescribed to it [24]. In contrast to the two novel models proposed (viz. log × hyp and log × exp), the adjustment function and ordinary differential equation of the Baranyi are integrated separately resulting in A(T), the explicit solution to the adjustment function, and y(t) the explicit solution to the ordinary differential equation and are derived specifically to describe logarithmic transformed microbial concentration data. These two equations result in the six-parameter model of Baranyi and Roberts [2]. Although more biologically thorough, the six-parameter Baranyi model is generally difficult to fit, leading to the previously mentioned re-parameterization and subsequent fixing of two shape parameters, resulting in the more widely used and robust four-parameter Baranyi model BAR. See Table 1 for explicit solutions to BAR. Regardless, the solution to BAR is still quite convoluted in comparison to the explicit solutions to log × hyp and log × exp, complicating the mathematics involved when manipulating this equation.
In contrast to BAR and MLOG, T is not explicitly represented in our LOG, log × hyp and log × exp models. Instead it has to be calculated as the intercept between the tangent to the steepest part of the growth equation and initial microbial biomass. When grown on the TYES liquid medium (Dataset 1), LOG, log × hyp and log × exp determined lag times to be 18.8, 18.7 and 18.3 h respectively, with all three models in agreement with maximum growth rate (units of microbial biomass per h) of 0.122 achieved on this medium. In terms of scaled maximum growth rate (h−1), LOG, log × hyp and log × exp were in close agreement with values of 0.17, 0.18 and 0.17 when grown on the TYES liquid medium. In comparison, when MLOG and BAR were fitted to the logarithmic transformation of this same dataset, lag times were determined to be 4.6 and 3.7 h, respectively. Maximum relative growth rates (h−1) of MLOG and BAR were determined to be 0.314 and 0.292, respectively. Due to the logarithmic transformation effecting the overall shape of a growth curve, it is clear that lag times are not comparable between models applied to absolute OD data and logarithmic transformed data, viz. LOG, log × hyp and log × exp versus MLOG and BAR, respectively. However, BAR and MLOG estimate lag times that are very much shorter than lag times when using the absolute untransformed OD data. Lag time estimates of MLOG and BAR are much shorter than the time at which OD starts to rise over the baseline OD0 using the actual untransformed data; therefore these two models appear to underestimate the real lag time. This can be seen through comparison of lag time values in Table 2 and visualization of Figure 1, in addition to the corresponding figures (Figures S1–S6 in the Supplementary Materials). The extent of underestimation of lag time for BAR and MLOG is more apparent in Study 1 compared to Study 2, whereby Study 1 is characterized by lower x0 values and a more prolonged lag phase.
The models MLOG and BAR represent the maximum growth parameter as a relative rate (h−1), while LOG, log × hyp and log × exp represent maximum growth in an absolute sense (units of microbial biomass per h) in addition to being able to represent maximum growth on a relative basis through scaled maximum growth rate (h−1). Very little information exists concerning parameters estimates resulting from models fitted to bacterial growth data in the field of aquaculture, with even less being available concerning F. psychrophilum. Stenholm et al. [50] reported maximum growth rates between 0.062 and 0.098 h−1 for 26 strains of F. psychrophilum. It is unclear however if these values were generated using a model, with the authors describing growth rates as being determined from the exponential increase in OD at 525 nm. Growth rates of two pathogenic bacteria Aeromonas hydrophila and Vibrio alginolyticus, obtained by fitting modified forms of the logistic and Gompertz, have been published [5]. However, the authors omitted the units when describing parameter estimates. Clearly, many factors must be standardized before meaningful comparisons of model-derived growth parameter estimates can be made. These factors include but are not limited to: method by which bacterial growth is measured, transformation applied to the dependent variable, and actual definition of growth rate [7].
The shortcomings of measuring microbial growth using the OD technique are worthy of note. The ability to derive lag and rate parameters from OD experiments is limited by the relatively high (~106–107 CFU/mL detection thresholds of OD measurement devices [4]. Therefore, this technique is best suited in circumstances whereby high cell densities are reached. In the field of predictive food microbiology spoilage resulting from bacterial load is a concern at high levels, making the absorbance technique practical in this field of study [4]. Likewise, high levels of the pathogenic bacteria used in this current study, F. psychrophilum, are commonly incubated to elicit typical clinical signs of rainbow trout fry syndrome. In a challenge test, Aoki et al. [51] inoculated ~106 CFU/mL of F. psychrophilum on a modified Cytophaga liquid medium, the same media as in Dataset 3 of this study. Bacteria reached the exponential phase of growth, at CFU/mL of ~107, and subsequently experimental fish were challenged with this load resulting in clinical signs of disease. Therefore, it is it feasible that the relatively high level of F. psychrophilum used in studies permits the use of the absorbance technique to determine lag and rate parameter estimates when studying this bacterium.
On the basis of goodness-of-fit and cross-validation, all the models fitted the data well. Compared to LOG and MLOG, the models log × hyp, log × exp and BAR are characterized by a flexible point of inflexion. In contrast, LOG and MLOG have fixed points of inflexion. For a logistic function, the fixed point of inflexion results in a symmetric sigmoidal curve. The major drawback of a fixed point of inflexion is that growth parameters are estimated at this maximum rate which always occurs at 0.5 of xm for the logistic [14]. The remaining equations display flexible points of inflexion that result in a more robust and versatile growth curve. The point of inflexion of these growth functions can all be calculated by equating the second derivative of the function to zero and solving for x* (the microbial biomass at the point of inflexion). However, log × hyp is the only equation out of log × hyp, log × exp and BAR whereby t* (time at which the point of inflexion occurs) and x* can be solved for analytically, with log × exp requiring numerical methods to solve for t* and the second derivative of the BAR unable to be solved in such a way that it is explicit in t, the independent variable [7]. Therefore, without the ability to solve for t* for the Baranyi, it is not possible to calculate maximum rate of actual growth using d x / d t = μ max α ( t * ) u ( x * ) x * . The ability of a growth function to be easily manipulated mathematically has important ramifications for their application and is an example of the advantage of log × hyp over the other functions examined.
The dampening functions of log × hyp and log × exp and the adjustment function of BAR both delay and suppress the exponential phase of growth. This is can be seen through visualization of Figure 2, where potential growth, represented by a logistic, is suppressed by a dampening function resulting in actual growth represented by log × hyp and log × exp. In addition, the resulting actual growth rate of log × hyp and log × exp is presented following suppression of the potential growth rate of the logistic. Using the parameter estimates from fitting log × hyp and log × exp to a given dataset, actual growth rate (Equation (9)) can be calculated and plotted as a function of time. The ability of log × hyp and log × exp to be manipulated mathematically allows for comparison of how growing conditions influence maximum rate of growth (s), absolute growth and rate of growth as displayed in Figure 3. The dampening functions of log × hyp and log × exp, and the adjustment function of BAR both suppress potential logistic growth but do so in different ways. The adjustment function of the Baranyi monotonically increases from a positive fraction to one, and represents the gradual adaptation of the bacteria to attain the maximum specific growth rate, as affected by their physiological state before and after inoculation. In contrast, the dampening functions of log × hyp and log × exp monotonically decrease from one to zero and represent the effect of culture conditions on growth rate. While the dampening functions of log × hyp and log × exp, viz. a rectangular hyperbola and a simple exponential, are seemingly empirical in nature some mechanistic meaning can be ascribed to them. For example, they may represent the effect of the build-up of an inhibitor on growth rate, or the effect of changes in growing conditions (e.g., acidity, temperature, water activity, oxygen levels) on growth rate over time. Indeed, a rectangular hyperbola is the standard representation in enzyme kinetics for relative velocity in the presence of a noncompetitive inhibitor [36]. The ability to describe these effects on bacterial growth rate add a level of novelty to this family of models. Without the adjustment function in BAR, and the dampening function of the log × models, these equations revert to a standard logistic form.

5. Conclusions

The objective of our research was to assess the ability of two novel functions, whose derivations are based upon principles promulgated by Baranyi and co-workers, to describe bacterial growth of the fish pathogen F. psychrophilum. Models were evaluated on the basis of statistical measures of goodness-of-fit, cross-validation and examination of residuals. Due to the nature of their derivation, log × hyp, log × exp and LOG were fitted to bacterial growth data in terms of measured OD over time, while BAR and MLOG were fitted to the logarithmic transformation of these data. Not only do log × hyp, log × exp and BAR perform better, or equally well, compared to their logistic counterparts (LOG and MLOG) in their ability to describe the time evolution in microbial mass, they also provide insight into the mechanisms that are driving bacterial growth through the concepts of potential and actual growth. The adjustment function of the Baranyi represents the suppression of growth that occurs when a bacterium adjusts from the pre- to post-inoculation environment, with the Baranyi model successfully describing bacterial growth under changing temperatures. Similarly, the two new models presented herein, log × hyp and log × exp, contain a dampening function, which represents the depression in the potential logistic growth, and theoretically encompasses a bacterium’s growth response to changing growing conditions. The two novel models have successfully described bacterial growth under constant external environmental conditions. The two elementary dampening functions represented in these equations, a rectangular hyperbola and a single exponential function, result in simple analytical equations in comparison with the existing Baranyi family of models in addition to describing rates of growth in both absolute and relative terms. From LOG and log × hyp functions it is possible to derive simple algebraic expressions to calculate inflexion point, maximum growth rate and scaled maximum growth rate. Based on fitting and predictive performance, the two newly derived models can be considered suitable alternatives to classical functions to describe bacterial growth and estimate key growth parameters.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/10/3/435/s1, Table S1: Observed growth data of Flavobacterium psychrophilum on eight liquid mediums, (1) TYES, (2) Shieh, (3) modified Cytophaga, (4) Cy7, (5) FLPB, (6) TYESB, (7) CBCM and (8) MAOB; Table S2: Resulting parameter estimates of log × hyp fit to Flavobacterium psychrophilum grown on eight liquid mediums: (1) TYES, (2) Shieh, (3) modified Cytophaga, (4) Cy7, (5) FLPB, (6) TYESB, (7) CBCM and (8) MAOB; Table S3: Resulting parameter estimates of log × exp fit to Flavobacterium psychrophilum grown on eight liquid mediums: (1) TYES, (2) Shieh, (3) modified Cytophaga, (4) Cy7, (5) FLPB, (6) TYESB, (7) CBCM and (8) MAOB; Figure S1: Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log ×exp (b) MLOG and BAR grown on a Sheih liquid medium; Figure S2: Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp (b) MLOG and BAR grown on a modified Cytophaga liquid medium; Figure S3: Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp (b) MLOG and BAR grown on a Cy7 liquid medium; Figure S4: Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp (b) MLOG and BAR grown on a TYESB liquid medium; Figure S5: Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp (b) MLOG and BAR grown on a CBCM liquid medium; Figure S6: Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp (b) MLOG and BAR grown on a MAOB liquid medium.

Author Contributions

Conceptualization, J.F.; methodology, J.F. and S.L.; software, C.D.P.; validation, J.F., S.L. and C.D.P; formal analysis, C.D.P; investigation, S.L. and C.D.P.; resources, C.D.P.; data curation, C.D.P.; writing—original draft preparation, C.D.P.; writing—review and editing, J.F., S.L. and C.D.P.; visualization, C.D.P.; supervision, J.F. and S.L.; project administration, J.F.; funding acquisition, J.F. All authors have read and agree to the published version of the manuscript.

Funding

This research was funded by The Canada Research Chairs program, grant number 045867 (Natural Sciences and Engineering Research Council of Canada, Ottawa).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Whiting, R.C. Notes on reparameterization of bacterial growth curves—A reply to J. Baranyi and W. E. Garthright. Food Microbiol. 1992, 9, 173–174. [Google Scholar] [CrossRef]
  2. Baranyi, J.; Roberts, T.A. A dynamic approach to predicting bacterial growth in food. Int. J. Food Microbiol. 1994, 23, 277–294. [Google Scholar] [CrossRef]
  3. López, S.; Prieto, M.; Dijkstra, J.; Dhanoa, M.S.; France, J. Statistical evaluation of mathematical models for microbial growth. Int. J. Food Microbiol. 2004, 96, 289–300. [Google Scholar] [CrossRef]
  4. Dalgaard, P.; Koutsoumanis, K. Comparison of maximum specific growth rates and lag times estimated from absorbance and viable count data by different mathematical models. J. Microbiol. Methods 2001, 43, 183–196. [Google Scholar] [CrossRef]
  5. Vine, N.; Leukes, W.; Kaiser, H. In vitro growth characteristics of five candidate aquaculture probiotics and two fish pathogens grown in fish intestinal mucus. FEMS Microbiol. Lett. 2004, 231, 145–152. [Google Scholar] [CrossRef] [Green Version]
  6. Dalgaard, P.; Ross, T.; Kamperman, L.; Neumeyer, K.; McMeekin, T.A. Estimation of bacterial growth rates from turbidimetric and viable count data. Int. J. Food Microbiol. 1994, 23, 391–404. [Google Scholar] [CrossRef]
  7. Perni, S.; Andrew, P.W.; Shama, G. Estimating the maximum growth rate from microbial growth curves: Definition is everything. Food Microbiol. 2005, 22, 491–495. [Google Scholar] [CrossRef] [Green Version]
  8. Campana, R.; Favi, G.; Baffone, W.; Lucarini, S. Marine alkaloid 2,2-bis(6-bromo-3-indolyl) ethylamine and its synthetic derivatives inhibit microbial biofilms formation and disaggregate developed biofilms. Microorganisms 2019, 7, 28. [Google Scholar] [CrossRef] [Green Version]
  9. Farber, R.; Rosenberg, A.; Rozenfeld, S.; Benet, G.; Cahan, R. Bioremediation of artificial diesel-contaminated soil using bacterial consortium immobilized to plasma-pretreated wood waste. Microorganisms 2019, 7, 497. [Google Scholar] [CrossRef] [Green Version]
  10. Deeb, J.G.; Smith, J.; Belvin, B.R.; Grzech-Leśniak, K.; Lewis, J. Er:YAG laser irradiation reduces microbial viability when used in combination with irrigation with sodium hypochlorite, chlorhexidine, and hydrogen peroxide. Microorganisms 2019, 7, 612. [Google Scholar] [CrossRef] [Green Version]
  11. Mauerhofer, L.M.; Pappenreiter, P.; Paulik, C.; Seifert, A.H.; Bernacchi, S.; Rittmann, S.K.M.R. Methods for quantification of growth and productivity in anaerobic microbiology and biotechnology. Folia Microbiol. (Praha). 2019, 321–360. [Google Scholar] [CrossRef] [Green Version]
  12. Zwietering, M.H.; Jongenburger, I.; Rombouts, F.M.; van ’T Riet, K. Modeling of the bacterial growth curve. Appl. Environ. Microbiol. 1990, 56, 1875–1881. [Google Scholar] [CrossRef] [Green Version]
  13. Buchanan, R.L.; Phillips, J.G. Response surface model for predicting the effects of temperature pH, sodium chloride content, sodium nitrite concentration and atmosphere on the growth of Listeria monocytogenes. J. Food Prot. 1990, 53, 370–376. [Google Scholar] [CrossRef]
  14. Garthright, W.E. Refinements in the prediction of microbial growth curves. Food Microbiol. 1991, 8, 239–248. [Google Scholar] [CrossRef]
  15. Baranyi, J.; Robinson, T.P.; Kaloti, A.; Mackey, B.M. Predicting growth of Brochothrix thermosphacta at changing temperature. Int. J. Food Microbiol. 1995, 27, 61–75. [Google Scholar] [CrossRef]
  16. Peleg, M. A model of microbial growth and decay in a closed habitat based on combined Fermi’s and the logistic equations. J. Sci. Food Agric. 1996, 71, 225–230. [Google Scholar] [CrossRef]
  17. Fujikawa, H.; Kai, A.; Morozumi, S. A new logistic model for bacterial growth. J. Food Hyg. Soc. Japan 2003, 44, 155–160. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, L. Growth kinetics of Listeria monocytogenes in broth and beef frankfurters - Determination of lag phase duration and exponential growth rate under isothermal conditions. J. Food Sci. 2008, 73, 235–242. [Google Scholar] [CrossRef]
  19. Huang, L. A new mechanistic growth model for simultaneous determination of lag phase duration and exponential growth rate and a new Bělehdrádek-type model for evaluating the effect of temperature on growth rate. Food Microbiol. 2011, 28, 770–776. [Google Scholar] [CrossRef]
  20. Gompertz, B. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. Philos. Trans. R. Soc. London 1825, 115, 513–583. [Google Scholar]
  21. Verhulst, P.F. Notice sur la loi que la population suit dans son accroissement. Corresp. Mathématique Phys. 1838, 10, 113–121. [Google Scholar]
  22. McMeekin, T.A. Predictive Microbiology: Theory and Application; Innovation in Microbiology Series 5; Research Studies Press Ltd.: Taunton, Somerset, UK, 1993; ISBN 0863801323. [Google Scholar]
  23. Baranyi, J.; Roberts, T.A.; McClure, P. A non-autonomous differential equation to model bacterial growth. Food Microbiol. 1993, 10, 43–59. [Google Scholar] [CrossRef]
  24. Baranyi, J. Simple is good as long as it is enough. Food Microbiol. 1997, 14, 189–192. [Google Scholar] [CrossRef]
  25. Buchanan, R.L.; Whiting, R.C.; Damert, W.C. When is simple good enough: A comparison of the Gompertz, Baranyi, and three-phase linear models for fitting bacterial growth curves. Food Microbiol. 1997, 14, 313–326. [Google Scholar] [CrossRef]
  26. Lorenzen, E.; Dalsgaard, I.; Bernardet, J. Characterization of isolates of Flavobacterium psychrophilum associated with coldwater disease or rainbow trout fry syndrome: Phenotypic and genomic studies. Dis. Aquat. Organ. 1997, 31, 197–208. [Google Scholar] [CrossRef]
  27. Lafrentz, B.R.; Lapatra, S.E.; Jones, G.R.; Cain, K.D. Passive immunization of rainbow trout, Oncorhynchus mykiss (Walbaum), against Flavobacterium psychrophilum, the causative agent of bacterial coldwater disease and rainbow trout fry syndrome. J. Fish Dis. 2003, 26, 377–384. [Google Scholar] [CrossRef]
  28. Starliper, C.E. Bacterial coldwater disease of fishes caused by Flavobacterium psychrophilum. J. Adv. Res. 2011, 2, 97–108. [Google Scholar] [CrossRef] [Green Version]
  29. Holt, R.A. Cytophaga psychrophila, the Causative Agent of Bacterial Cold-Water Disease in Salmonid Fish. Ph.D. Thesis, Oregon State University, Corvallis, OR, USA, 1987. [Google Scholar]
  30. Wiens, G.D.; LaPatra, S.E.; Welch, T.J.; Evenhuis, J.P.; Rexroad III, C.E.; Leeds, T.D. On-farm performance of rainbow trout (Oncorhynchus mykiss) selectively bred for resistance to bacterial cold water disease: Effect of rearing environment on survival phenotype. Aquaculture 2013, 388, 128–136. [Google Scholar] [CrossRef]
  31. Bruun, M.S.; Schmidt, A.; Madsen, L.; Dalsgaard, I. Antimicrobial resistance patterns in Danish isolates of Flavobacterium psychrophilum. Aquaculture 2000, 187, 201–212. [Google Scholar] [CrossRef]
  32. Shieh, H.S. Studies on the nutrition of a fish pathogen, Flexibacter columnaris. Microbios Lett. 1980, 13, 129–133. [Google Scholar]
  33. Wakabayashi, H. Characteristics of myxobacteria associated with some freshwater fish diseases in Japan. Bull. Japanese Soc. Sci. Fish 1974, 40, 751–757. [Google Scholar] [CrossRef]
  34. Reichenbach, H.; Dworkin, M. Introduction to the gliding bacteria. In The Prokaryotes; Starr, M.P., Stolp, H., Trüper, H.G., Balows, A., Schlegel, H.G., Eds.; Springer: Berlin, 1981; pp. 315–327. [Google Scholar]
  35. Cepeda, C.; García-Márquez, S.; Santos, Y. Improved growth of Flavobacterium psychrophilum using a new culture medium. Aquaculture 2004, 238, 75–82. [Google Scholar] [CrossRef]
  36. Segel, I.H. Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady-State Enzyme Systems; Wiley Classics Library; Wiley: New York, NY, USA, 1993; ISBN 0471303097. [Google Scholar]
  37. Miller, G. Numerical Analysis for Engineers and Scientists; Cambridge University Press: Cambridge, UK, 2014; ISBN 9781107021082. [Google Scholar]
  38. SAS Institute Inc. SAS/STAT® 9.2 User’s Guide, 2nd ed.; SAS Institute Inc.: Cary, NC, USA, 2009; ISBN 978-1607644392. [Google Scholar]
  39. Bibby, J.; Toutenburg, H. Prediction and Improved Estimation in Linear Models; John Wiley and Sons: Chichester, UK, 1977; ISBN 0-471-01656-X. [Google Scholar]
  40. Lin, L.I.-K. A concordance correlation coefficient to evaluate reproducibility. Biometrics 1989, 45, 255–268. [Google Scholar] [CrossRef]
  41. Akaike, H. A new look at the statistical model identification. IEEE Trans. Automat. Contr. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  42. Ross, T. Indices for performance evaluation of predictive models in food microbiology. J. Appl. Bacteriol. 1996, 81, 501–508. [Google Scholar] [CrossRef]
  43. Baranyi, J.; Pin, C.; Ross, T. Validating and comparing predictive models. Int. J. Food Microbiol. 1999, 48, 159–166. [Google Scholar] [CrossRef]
  44. Draper, N.R.; Smith, H. Applied Regression Analysis, 3rd ed.; Wiley Series in Probability and Statistics; John Wiley & Sons: New York, NY, YSA, 1998; ISBN 9780471170822. [Google Scholar]
  45. Allen, D.M. The relationship between variable selection and data augmentation and a method for prediction. Technometrics 1974, 16, 125–127. [Google Scholar] [CrossRef]
  46. Tarpey, T. A note on the prediction sum of squares statistic for restricted least squares. Am. Stat. 2000, 54, 116–118. [Google Scholar] [CrossRef]
  47. Holiday, D.B.; Ballard, J.E.; McKeown, B.C. PRESS-related statistics: Regression tools for cross-validation and case diagnostics. Med. Sci. Sports Exerc. 1995, 27, 612–620. [Google Scholar] [CrossRef]
  48. Richards, F.J. A flexible growth function for empirical use. J. Exp. Bot. 1959, 10, 290–301. [Google Scholar] [CrossRef]
  49. Dalgaard, P. Modelling of microbial activity and prediction of shelf life for packed fresh fish. Int. J. Food Microbiol. 1995, 305–317. [Google Scholar] [CrossRef]
  50. Stenholm, A.R.; Dalsgaard, I.; Middelboe, M. Isolation and characterization of bacteriophages infecting the fish pathogen Flavobacterium psychrophilum. Appl. Environ. Microbiol. 2008, 74, 4070–4078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Aoki, M.; Kondo, M.; Kawai, K.; Oshima, S.I. Experimental bath infection with Flavobacterium psychrophilum, inducing typical signs of rainbow trout Oncorhynchus mykiss fry syndrome. Dis. Aquat. Organ. 2005, 67, 73–79. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp and (b) MLOG and BAR, grown on a TYES liquid medium, and (c) LOG, log × hyp, log × exp and (d) MLOG and BAR, grown on a FLBP liquid medium. In panels (a) and (c), models are fitted to untransformed optical density measurements (OD) while (b) and (d) models are fitted to [ln(OD/OD0)].
Figure 1. Optical density growth data and model predictions for Flavobacterium psychrophilum. Growth predictions using (a) LOG, log × hyp, log × exp and (b) MLOG and BAR, grown on a TYES liquid medium, and (c) LOG, log × hyp, log × exp and (d) MLOG and BAR, grown on a FLBP liquid medium. In panels (a) and (c), models are fitted to untransformed optical density measurements (OD) while (b) and (d) models are fitted to [ln(OD/OD0)].
Animals 10 00435 g001
Figure 2. The effect of dampening potential growth rate, resulting actual growth, and actual rate of growth (AGR) of Flavobacterium psychrophilum. Dampening effect is represented by: a rectangular hyperbola (a), a simple exponential (b) when grown on a TYES liquid medium and a rectangular hyperbola (c), a simple exponential (d) when grown on FLBP liquid medium.
Figure 2. The effect of dampening potential growth rate, resulting actual growth, and actual rate of growth (AGR) of Flavobacterium psychrophilum. Dampening effect is represented by: a rectangular hyperbola (a), a simple exponential (b) when grown on a TYES liquid medium and a rectangular hyperbola (c), a simple exponential (d) when grown on FLBP liquid medium.
Animals 10 00435 g002
Figure 3. Comparison of actual growth rate (AGR), predicted growth and maximum actual growth rate (s) from liquid media resulting in highest and lowest maximum growth rates (s) from Study 1 and Study 2 when fitting log × hyp and log × exp. Panels (a) and (b) are taken from Study 1, and panels (c) and (d) from Study 2.
Figure 3. Comparison of actual growth rate (AGR), predicted growth and maximum actual growth rate (s) from liquid media resulting in highest and lowest maximum growth rates (s) from Study 1 and Study 2 when fitting log × hyp and log × exp. Panels (a) and (b) are taken from Study 1, and panels (c) and (d) from Study 2.
Animals 10 00435 g003
Table 1. Evaluated growth functions.
Table 1. Evaluated growth functions.
Functional Form
Simple logistic
(LOG)
x = x 0 x m x 0 + ( x m x 0 ) e μ t
Modified logistic (MLOG) x = x m { 1 + exp [ 4 μ max x m ( T t ) + 2 ] }
Baranyi four parameter
(BAR)
x = x 0 + μ max A ( t ) ln [ 1 + e μ max A ( t ) 1 e ( x m x 0 ) ]
A ( t ) = t + 1 μ max ln [ e μ max t + e h 0 e μ max t h 0 ]
where h 0 = µmaxT = ln[1 + 1 q 0 ]
log × hyp x = x 0 x m x 0 + ( x m x 0 ) ( λ λ + t ) λ μ
log × exp x = x 0 x m x 0 + ( x m x 0 ) e x p [ μ λ ( 1 e λ t ) ]
Table 2. Comparison of maximum growth rate (s, units of microbial mass per h), scaled maximum growth rate (s/x*, h−1), lag time (T, h) and maximum relative growth rate (µmax, h−1), generated by the standard logistic (LOG), logistic × hyperbola (log × hyp) and logistic × exponential (log × exp) fitted to optical density (OD) growth data and the modified logistic (MLOG) and four-parameter Baranyi (BAR) fitted to logarithmic transformed growth data [ln(OD/OD0)].
Table 2. Comparison of maximum growth rate (s, units of microbial mass per h), scaled maximum growth rate (s/x*, h−1), lag time (T, h) and maximum relative growth rate (µmax, h−1), generated by the standard logistic (LOG), logistic × hyperbola (log × hyp) and logistic × exponential (log × exp) fitted to optical density (OD) growth data and the modified logistic (MLOG) and four-parameter Baranyi (BAR) fitted to logarithmic transformed growth data [ln(OD/OD0)].
LOGlog × hyplog × expMLOGBAR
Datasetss/x*Tss/x*Tss/x*TµmaxTµmaxT
10.1220.16818.80.1220.18218.70.1220.17018.30.3144.60.2923.7
20.0960.13816.80.0950.15116.60.0970.16317.10.3625.60.3485.3
30.1220.19119.80.1240.20719.80.1230.19719.50.3025.50.2915.1
40.0520.16524.50.0530.18124.70.0520.16824.20.2197.80.2217.7
50.1110.1033.10.1240.1741.70.1370.1611.60.2851.70.3341.8
60.1330.1613.70.1450.2253.50.1480.2373.60.3092.90.5064.4
70.1200.2064.10.1190.2174.00.1070.2794.80.8477.61.6568.3
80.0310.1301.40.0310.1461.20.0300.1300.70.4577.71.4528.7
Table 3. Goodness-of- fit: Akaike information criterion (AIC), mean square prediction error (MSPE) concordance correlation coefficient (CCC) and accuracy factor (AF) values obtained when fitting models to microbial growth data from Study 1 and 2.
Table 3. Goodness-of- fit: Akaike information criterion (AIC), mean square prediction error (MSPE) concordance correlation coefficient (CCC) and accuracy factor (AF) values obtained when fitting models to microbial growth data from Study 1 and 2.
Measure of Goodness-of-FitLOGlog × hyplog × expMLOGBAR
AIC
 Study 1-Average
 (±SE)
−117.1
(11.5)
−110.0
(11.5)
−111.6
(11.3)
−58.8
(3.1)
−67.5
(5.5)
 SStudy 2-Average
 (±SE)
−44.6
(3.9)
−43.9
(2.7)
−45.0
(2.9)
−44.4
(3.6)
−44.0
(2.9)
 SOverall-Average
 (±SE)
−80.8
(15.1)
−77.0
(14.0)
−78.3
(14.0)
−51.6
(3.7)
−55.7
(5.6)
MSPE
 SStudy 1-Average
 (±SE)
0.001
(0.0002)
0.001
(0.0003)
0.001
(0.0002)
0.024
(0.003)
0.014
(0.004)
 SStudy 2-Average
 (±SE)
0.005
(0.002)
0.004
(0.001)
0.003
(0.001)
0.005
(0.002)
0.004
(0.001)
 SOverall-Average
 (±SE)
0.003
(0.001)
0.002
(0.001)
0.002
(0.001)
0.013
(0.004)
0.009
(0.003)
CCC
 SStudy 1-Average
 (±SE)
0.998
(0.001)
0.998
(0.001)
0.998
(0.001)
0.997
(0.001)
0.998
(0.001)
 SStudy 2-Average
 (±SE)
0.987
(0.004)
0.988
(0.005)
0.989
(0.005)
0.997
(0.001)
0.998
(0.001)
 SOverall-Average
 (±SE)
0.992
(0.003)
0.993
(0.003)
0.993
(0.003)
0.997
(0.001)
0.998
(0.001)
AF
 SStudy 1-Average
 (±SE)
1.013
(0.002)
1.016
(0.003)
1.015
(0.002)
1.039
(0.005)
1.024
(0.005)
 SStudy 2-Average
 (±SE)
1.033
(0.005)
1.030
(0.003)
1.027
(0.002)
1.022
(0.007)
1.016
(0.003)
 SOverall-Average
 (±SE)
1.023
(0.005)
1.023
(0.003)
1.021
(0.002)
1.031
(0.005)
1.020
(0.003)
Table 4. Examination of residuals using the runs test, Durbin–Watson (DW) and averaged bias factor (BF) statistics obtained when fitting models to growth data from Study 1 and 2.
Table 4. Examination of residuals using the runs test, Durbin–Watson (DW) and averaged bias factor (BF) statistics obtained when fitting models to growth data from Study 1 and 2.
Test for Examination of ResidualsLOGlog × hyplog × expMLOGBAR
Runs test
 Runs were random65688
 Too few runs23200
No. curves exhibiting serial correlation determined by DW statistic (α = 0.01)
 No serial correlation77788
 Positive Correlation11101
BF
 SStudy 1-Average
 (±SE)
0.996
(0.002)
0.993
(0.001)
0.994
(0.000)
1.017
(0.001)
1.003
(0.001)
 SStudy 2-Average
 (±SE)
1.001
(0.003)
0.994
(0.002)
0.994
(0.003)
1.006
(0.003)
1.001
(0.000)
 SOverall-Average
 (±SE)
0.999
(0.002)
0.994
(0.001)
0.994
(0.001)
1.011
(0.003)
1.002
(0.001)
Table 5. Predicted residual error sum of squares (PRESS) and concordance correlation coefficient (CCC) resulting from cross-validation of Datasets 1-2 from Study 1 and Datasets 5-6 from Study 2.
Table 5. Predicted residual error sum of squares (PRESS) and concordance correlation coefficient (CCC) resulting from cross-validation of Datasets 1-2 from Study 1 and Datasets 5-6 from Study 2.
Cross-Validation TestModel
LOGlog × hyplog × exp
PRESS
 Dataset 10.07170.07180.0599 *
 Dataset 20.0110 *0.01250.0134
 Dataset 50.26400.21350.1858 *
 Dataset 60.11840.11530.1149 *
CCC
 Dataset 10.99450.99450.9953 *
 Dataset 20.9990 *0.99890.9988
 Dataset 50.97260.97820.9814 *
 Dataset 60.98060.98070.9808 *
(*) represents the model which performed best, on the basis of PRESS and CCC, for a given dataset.

Share and Cite

MDPI and ACS Style

Powell, C.D.; López, S.; France, J. New Insights into Modelling Bacterial Growth with Reference to the Fish Pathogen Flavobacterium psychrophilum. Animals 2020, 10, 435. https://doi.org/10.3390/ani10030435

AMA Style

Powell CD, López S, France J. New Insights into Modelling Bacterial Growth with Reference to the Fish Pathogen Flavobacterium psychrophilum. Animals. 2020; 10(3):435. https://doi.org/10.3390/ani10030435

Chicago/Turabian Style

Powell, Christopher D., Secundino López, and James France. 2020. "New Insights into Modelling Bacterial Growth with Reference to the Fish Pathogen Flavobacterium psychrophilum" Animals 10, no. 3: 435. https://doi.org/10.3390/ani10030435

APA Style

Powell, C. D., López, S., & France, J. (2020). New Insights into Modelling Bacterial Growth with Reference to the Fish Pathogen Flavobacterium psychrophilum. Animals, 10(3), 435. https://doi.org/10.3390/ani10030435

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop