Next Article in Journal
Self-Driving Car Location Estimation Based on a Particle-Aided Unscented Kalman Filter
Previous Article in Journal
Global Muscle Coactivation of the Sound Limb in Gait of People with Transfemoral and Transtibial Amputation
Previous Article in Special Issue
Portable Ultrasound-Based Device for Detecting Older Adults’ Sit-to-Stand Transitions in Unsupervised 30-Second Chair–Stand Tests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measurement and Modeling of Microbial Growth Using Timelapse Video

by
Konstantinos Delibasis
1,
Ifigenia Basanou
1 and
Alexandros-Apostolos A. Boulogeorgos
2,*
1
Department of Computer Science and Biomedical Informatics, University of Thessaly, 35100 Lamia, Greece
2
Department of Digital Systems, University of Piraeus, 18534 Piraeus, Greece
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(9), 2545; https://doi.org/10.3390/s20092545
Submission received: 29 February 2020 / Revised: 8 April 2020 / Accepted: 13 April 2020 / Published: 29 April 2020
(This article belongs to the Special Issue Sensor and Systems Evaluation for Telemedicine and eHealth)

Abstract

:
The development of timelapse videos for the investigation of growing microbial colonies has gained increasing interest due to its low cost and complexity implementation. In the present study, a simple experimental setup is proposed for periodic snapshot acquisition of a petri dish cultivating a fungus of the genus Candida SPP, thus creating a timelapse video. A computational algorithm, based on image processing techniques is proposed for estimating the microbial population and for extracting the experimental population curves, showing the time evolution of the population of microbes at any region of the dish. Likewise, a novel mathematical population evolution modeling approach is reported, which is based on the logistic function (LF). Parameter estimation of the aforementioned model is described and visually assessed, in comparison with the conventional and widely-used LF method. The effect of the image analysis parameterization is also highlighted. Our experiments take into account different area sizes, i.e., the number of pixels in the neighborhood, to generate population curves and calculate the model parameters. Our results reveal that, as the size of the area increases, the curve becomes smoother, the signal-to-noise-ratio increases and the estimation of model parameters becomes more accurate.

1. Introduction

During the last decades, image capturing and processing have evolved from scarce technology, used by professionals for niche applications, to a rapidly advancing research field with application in low-cost and complexity data acquisition of long duration physical and/or biological phenomena [1]. To this end, a great amount of research effort has been spent on designing image processing and machine learning-based approaches that not only quantify the population growth for image sequences, but also provide accurate predictions concerning its evolution.
Scanning the open technical literature, there are several published papers that report the use of images and logistic functions (LFs) for predicting a specific-cell type population growth [2,3,4,5,6,7,8,9,10]. In more detail, in [2], the use of eight consecutive images, which were acquired by a single reflex (SLR) camera, was presented, in conjunction with the probability of survival of hamster stem cells, which was modeled according to the LF approach. Additionally, in [3], the authors showed that the in-vitro cellular growth rate may be used as survival probability of the fetus in case of hamsters, while, in [4], a multi-camera approach for modeling larva population by means of LF was presented. Another example of timelapse application in population behavior is delivered in [5], where the insect Costelytra zealandica was infected with the bacteria Yersinia entomophaga and its behavior was observed and recorded by acquiring 1 image per 10 min , using a digital SLR camera. Fluorescence substances were used on the insects to facilitate analysis of the resulting image sequence. More recently, in [6], the construction of a special incubation chamber was described that uses two petri dishes to control humidity for a long duration. Images were acquired every 10 min for 40 days and used to record and study the morphology of bacterial colonies, such as streptomycin, bacillus subtillis E42, serratia marcescens, arthrobacter agilis, and nestekonia SP. The timelapse technique has also been applied to agriculture [7], embryo moprhokinetics [8] and microbial interaction [9,10].
From the data processing point of view, a number of LF-based functions have been employed for modeling and/or predicting the population evolution based on the birth and death rates [11,12]. In this direction, in [13,14], the bi-LF was introduced to model population evolution that show two overlapping or sequential phases of logistic growth. The two LFs were summed and fitted to rice crop evolution data to handle an initial growth increase followed by a decrease. The combined curve provided a good fit to the data under investigation. Of note, non-linear least squares are used to estimate the parameter values of the function. Meanwhile, in [15], the LF is augmented by using an additive constant to accommodate varying capacity. This method was applied to two actual cases of human population growth. Moreover, in [16], the bi-phasic LF was used for modeling the sudden changes in the evolution of research papers and patents. The combination of two LFs was treated in a systematic way, using the following categories: (i) Sequential (sum of two time translated curves), (ii) superposed simultaneous curves, (iii) converging (first logistic growth combined with a second logistic function, reaching the maximum at about the same time) and (iv) diverging (different growth rates and carrying capacities). In all aforementioned cases, the resulting curves were unable to accommodate data decreasing trends. Additionally, in [17], a number of different extensions of the LF were discussed including the generalized logistic regression, the Von Bertalanffy, and the Blumberg expressions. Again, these expressions were unable to model decreasing phases of the data. Finally, in [18], the combertzian population growth model, which is a double exponential model, with parameters that control the time shift, the maximum asymptotic population and the maximum growth rate, was applied to investigate the evolution of breast cancer cells.
To the best of the authors’ knowledge, the problems of modeling the decreasing phase of the data have not been adequately addressed in the technical literature. Motivated by this, the current contribution presents an appropriate growth approach that is able to model both the increasing and decreasing phases and verifies its efficiency by means of experiment-driven microbial population modeling. In particular, we propose an experimental setup for periodic snapshot recording of the evolution of a fungus on a petri dish. The selected fungus is the Candida species, which is the most common cause of fungal infections [19,20,21,22]. The acquired images are processed in the red–green–blue (RGB) color system to quantify the number of microbes in each pixel of the plate. A novel mathematical model of the time evolution of the fungus for individual pixel/regions of the plate is presented, which introduces an LF with second degree polynomial time exponent capable of modeling population curves that may eventually become decreasing. Finally, we apply parametric imaging to depict the spatial distribution of a quantity proportional to the growth rate of the microbial population. The proposed approach is expected to open a new road to automatic data acquisition, modeling, and prediction of different evolution processes. The results are further validated and the robustness of the proposed method assessed using two more publicly available timelapse videos with a different acquisition method and microbial colonies of different species.
The rest of this paper is organized as follows: Section 2 is focused on presenting the methodology that we follow. In more detail, the experimental setups accompanied by the developed theoretical framework are reported in this section. Next, in Section 3 experimental results that verify the proposed analysis and the accuracy of the proposed modeling function are presented. Finally, useful remarks and conclusions summarizing this work are provided in Section 4.
Notations: Unless otherwise stated, in what follows, lower and upper case bold letters denote vectors and matrices, respectively. Moreover, A B stand for the union of the sets A and B , while ⊕ and •, respectively, represent the morphological dilation and erosion. Additionally, ( · ) T and ( · ) 1 , respectively, stand for the transpose and inverse operator. Finally, exp · and ln · , respectively, denote the exponential and the natural logarithm.

2. Methodology

This section aims at reporting the methodology that is used for deriving and modeling the experimental data with the aid of LF with second degree (LFSD) exponent. In particular, in Section 2.1, the experimental setup is described, while the approach that we use to estimate the microbe population evolution and generate its curves is reported in Section 2.2. Finally, Section 2.3 presents the proposed time-shifted LFSD exponent approach, while the parametric imaging is reported in Section 2.4.

2.1. Experimental Setup

The experiment took place at the Department of Computer Science and Biomedical Informatics. The experimental setup consisted of a Logitec 920 camera, used to acquire images of the petri dish that was placed inside the incubator. As illustrated in Figure 1, an artificial light source of a white, 6 Watt light emitting diode (LED) was used for illumination. Images were captured on an inexpensive laptop using the free version of video-velocity software (http://www.candylabs.com/videovelocity). The Candida SPP was used as a fungus. The evolution of Candida SPP took place on a Biomeriex cos petri dish. The red color of the plate is convenient for the post processing of the acquired images. The camera was placed directly above the petri dish, at an approximate distance of 15 cm . Camera focus was set manually and motion detection was disactivated. Image resolution was set to 360 × 640 pixels. The period of frame acquisition was set to 15 s and the fungus was allowed to grow for 50 h , 24 min and 36 s , resulting in a total of 12 , 099 jpeg color images. The first 3800 of them contained no evidence of fungus evolution; thus, they were excluded from the study to facilitate further processing, leaving 8300 to participate in further processing. Indicative examples of the microbial evolution are presented in Figure 2.
In order to extent the validation of the proposed algorithm, as well as test its robustness, two more time lapse movies of microbial/fungi growth in petri dish were used, namely petri dish 2 and 3, respectively. In petri dish 2, a number of different cultures were recorded on a blood sheep of similar color as the one that we used in our experiment. The time lapse movie was obtained from Youtube (https://www.youtube.com/watch?v=XLmk0zYIFjE) and consists of 880 frames of 720 lines by 1280 columns. For petri dish 3, another timelapse movie was obtained as a publicly available file, described in [23]. The movie was acquired using a scanner, with a period of one frame every 5 min . A total number of 2250 frames consisting of 720 lines by 1280 columns were extracted and processed. The frames were cropped inside the dish.

2.2. Image-Based Measurement of Microbe Population and Generation of Experimental Population Curves

As can be observed in Figure 2, the sheep blood agar plate 1 and 2 has a distinct red color. The fungal colonies color is different than the plate’s; thus, we make the fundamental assumption that the number of microbes of each imaged pixel is proportional to difference between the two colors. Formally, the number of microbes at each pixel ( i , j ) at time t, can be estimated by the difference of the red channel value R i , j t from the corresponding value at the initial frame R i , j 0 . The total number of microbes are computed for all pixels in the current frame t.
The pixels that lie outside the plate should be excluded. In this direction, a binary mask, M i , j , is defined with size equal to the one of the frame, which ideally should be 1 in all pixels inside the plate and 0 otherwise. In other words, the pixels inside the plate can be detected by using the following heuristic condition:
M i , j = 1 R i , j 0 2 G R i , j 0 2 B i , j 0 ,
where G i , j 0 and B i , j 0 are, respectively, the values of the green and blue channel of pixel ( i , j ) in the initial frame. The binary mask in Equation (1) is further processed using the following morphological operation:
M ˜ = M B 1 B 2 ,
where M is a matrix with elements M i , j , while B 1 and B 2 are, respectively, square structuring elements of size 11 × 11 and 31 × 31 .
The total number of microbes in the current frame t can be obtained as
N t = i , j Δ i , j t M ˜ i , j ,
where M ˜ i , j is the i , j element of the matrix M ˜ , while
Δ i , j t = R i , j t R i , j 0 , R i , j t R i , j 0 T 0 , otherwise ,
with T being a threshold, which is applied in order to prevent insignificant color changes from contributing to the population calculation. The threshold is experimentally determined and is set to 10. The total estimated number of microbes in the whole plate for each frame is illustrated in Figure 3.
In the case of the third petri dish, which is radically different than the red-colored circular petri dishes 1 and 2, the following minor modifications were made to the image-based measurement algorithm. Since the whole frame contains microbial colonies there is no need to define the binary mask as in Equations (1)–(3) can be simplified as M i j = 1 one for all i , j . Finally, since the original color of the plate is black, Equation (4) should be modified, in order to take into account all chromatic channels, as
Δ i j = max r i j , g i j , b i j ,
where
r i j = | R i j t R i j 0 | , for | R i j t R i j 0 | T 0 , otherwise ,
g i j = | G i j t G i j 0 | , for | G i j t G i j 0 | T 0 , otherwise ,
and
b i j = | B i j t B i j 0 | , for | B i j t B i j 0 | T 0 , otherwise .
From Equation (3), it becomes evident that the population of microbes in a region of interest (ROI), A, of the plate can be evaluated as
N t i , j A Δ i , j .
The ROI is a square defined round its central pixel of size 2 r + 1 × 2 r + 1 . An indicative example is presented in Figure 4a, where the result of applying Expression (9) to each pixel of the last frame for r = 8 pixel is provided. Figure 4b depicts a typical example of an experimental population curve (blue) accompanied by the smoothed population one (black). The experimental curves show statistical fluctuations, due to inaccuracies of the estimation method of the population, which are considered as noise. In order to increase the accuracy of population modeling function, smoothing of the extracted curve is necessary. In this direction, we employ a 27-element median filter.
The instantaneous signal-to-noise ratio (SNR) is defined as
γ ( t ) = s ( t ) n ( t ) ,
where s ( t ) and n ( t ) are, respectively, the signal and noise amplitudes. Note that the signal amplitude is measured as the mean value of the curve, while the noise amplitude is calculated as the standard deviation of the curve at the same time interval. Let us consider t p points before and after the current frame, t, then s ( t ) and n ( t ) can be, respectively, obtained as
s ( t ) = 1 2 t p + 1 τ = t t p t + t p N τ ,
n ( t ) = 1 2 t p + 1 τ = t t p t + t p N τ s ( τ ) .

2.3. The Proposed Time-Shifted LFSD Exponent

At the onset of the LF model, the population grows exponentially, but after some time, the rate of increase reduces until it reaches zero. Thus, the LF approximates asymptotically a fixed value that represents a theoretical maximum limit to the population size under given conditions, such as the resource availability decrease, toxic products concentration, etc. Let p LF ( t ) denote the percentage of the population maximum at time t, which can be expressed as
p LF ( t ) = 1 1 + exp b 1 t + b 0 ,
where b 0 and b 1 are LF-specific parameters.
As the colony evolves further, it is possible for its actual population curve to exhibit a downward trend and become decreasing, due to lack of nutrients, or the accumulation of toxic metabolic products. In order to model this type of population trend, we modify the conventional LF by introducing a second degree polynomial as exponent. The LFSD can be written as
p LFSD ( t ) = 1 1 + exp b 2 t 2 + b 1 t + b 0 .
In contrary to the conventional LF, the LFSD cannot be easily adjusted to fit time-translated experimental data. That would require the exponent of LFSD to become b 2 t t s 2 + b 1 t t s + b 0 , introducing an extra parameter t s , which cannot be determined by solving a linear system of equations. Instead, we alleviate this problem by measuring experimentally the time t s of the first appearance of non-zero experimental data. Thus, given the values of the population N i , for i = 1 , , n , where i represents different time values, LFSD can be rewritten as
p LFSD ( t + t s ) = 1 1 + exp b 2 t 2 + b 1 t + b 0 ,
where b 0 , b 1 , and b 2 are LFSD-specific parameters. In order to use Equation (15) for modeling the experimental population curve of a colony on the plate, we have to normalize the population curve in the range 0 , 1 . In this direction, we set
p t = N t N max + e ,
where N max is the maximum value of the population curve and e is a small positive number.
From Equation (15), we can obtain a linear form that facilitates the calculation of parameter b 0 , b 1 , and b 2 , as
y i = ln p LFSD ( t i + t s ) 1 p LFSD ( t i + t s ) ,
or equivalently
y i = b 2 t i 2 + b 1 t i + b 0 .
Given the values of ( t i , y i ) with i = 1 , , N , the parameters b 0 , b 1 , and b 2 can be obtained as the solution of
T b = y ,
where b = b 1 , b 2 , b 3 T , y = y 1 , , y N , and
T = t 1 2 t 1 1 t N 2 t N 1 .
From Equation (19), b can be calculated as
b = T T T 1 T T y .
Theorem 1.
The maximum value of the population curve can be achieved at
t max = b 1 2 b 2 + t s .
Proof. 
The growth rate can be evaluated as
ρ = d y i d t i ,
or equivalently
ρ = d y i d p L F S D ( t i + t s ) d p L F S D ( t i + t s ) d t i ,
which can be analytically written as
ρ = C 2 b 2 t i + b 1 exp b 2 t i 2 + b 1 t i + b 0 ,
where
C = 1 p L F S D ( t i + t s ) 1 p L F S D ( t i + t s ) .
In order to evaluate t max , we set ρ = 0 and we get Equation (22). □
Remark 1.
Note that if the population curve is monotonically increasing during the time period of image acquisition, then Equation (22) returns a value outside the time frame. This value should be discarded.
An example of determining the model parameters for the experimental population curve of the plate that exhibits a late decreasing phase is shown in Figure 5. Modeling the experimental data using LF and LFSD is performed, respectively shown in pink and green curve. From this figure, it becomes evident that the LFSD achieves more accurate modeling compared to the LF, because the former is able to model the last phase, when the population declines, while, in the deceasing phase, the latter increases monotonically with its slope approaching zero.

2.4. The Parametric Image

It is possible to apply the experimental population curve extraction and its parameter estimation to each pixel of the frame, so that the values of the model parameters can be calculated separately for each non-zero pixel and displayed as a parametric image. In this study, parametric images are constructed for the parameters b 2 and b 1 of the LFSD, b 1 of the LF, t max as calculated by Equation (22) from the LFSD and t 0 . Therefore, the parametric image shows regions of the plate with a different population growth rate. Figure 6 shows the value of the coefficient b 2 for the LFSD model for ROI with r = 4 .

3. Results

This section is devoted to presenting experimental results that evaluate the effectiveness and accuracy of the proposed approach. In this direction, in Figure 7, an indicative example is depicted, where, in the final acquired frame, two regions of the plate, which contain several colonies, are selected. The extracted population curves are also presented in Figure 7d,e. It is evident that the curve of the smaller blue area, which is depicted in Figure 7b,d, shows a decreasing trend. Finally, as illustrated in Figure 7c,e, the evolution curve of the yellow area has a monotonically increasing trend.
Next, the influence of the selected area size is explored. First, a pixel is selected and the experimental population curve is calculated for different sizes of the square region centered at the selected pixel. As depicted in Figure 8, for each experimental curve, the modeling is performed using both the LF and LFSD approaches. Moreover, the LFSD approach is observed to be more accurate compared to the LF. Interestingly, it becomes obvios that as the value of r decreases, the population curve becomes more noisy.
For each r parameter, the coefficients of the exponential polynomial in the case of LF and LFSD are, respectively, reported in Table 1 and Table 2. Notice that b 1 in the case of LF and b 2 , b 1 of LFSD are population growth rate metrics. The consistency of the curve extraction and proposed modeling is shown in Figure 9, where t max is calculated using Equation (22). It can be observed that t max does not vary substantially even for drastically r changes.
To highlight the superiority of LFSD against the conventional LF, in Figure 10a, we present a population evolution, which after a specific frame is decreasing. Figure 10b illustrates the estimated population against the LF-based model. From this figure, it becomes evident that the conventional LF is incapable of modeling this behavior. This highlights the importance of LFSD for modeling such trends. The number of pixels with negative b 1 in the conventional LF in the whole petri dish is also reported in Table 3. We notice that the number of pixels having negative exponent b 1 decreases as the radius of the ROI area increases (Figure 11). However, the proposed LFSD does not suffer from this shortcoming.
The extracted population curves exhibit noticeable fluctuation, which appears to increase as the local mean value of the curve increases. Figure 9 plots the instantaneous signal-to-noise ratio (SNR) for the population curve shown in Figure 8. The parameter t p is set to 80. The logarithmic definition of SNR was used. It becomes apparent that as the population increases with time, both the amplitude of signal and noise increases, however the signal amplitude increases faster than the noise (curve’s local standard deviation); thus, the SNR overall increases. Furthermore, as expected, the noise decreases as the size of the ROI used for the calculation increases.
The time t s of the first appearance of non-zero experimental data for each pixel is shown as a parametric image in Figure 12.
Similarly, t max is shown in Figure 13a and also superimposed on the petri dish in Figure 13b, using shades of green.
Finally, Figure 14 shows parametric images of b 2 of LF, with their color scale, for two different radii r = 4 and 8 pixels, respectively. Visual inspection of the parametric images b 1 and considering the influence of the b 1 factor on the population curve, suggests that areas with a high parameter values exhibit a low population density.
In Figure 15, four locations are shown on different colonies of petri dish 2 in the final frame. The population curves that were extracted and modeled using the LF and the proposed LFSD, are shown in Figure 16. Again, it can be observed that the LFSD achieves more accurate population modeling in the cases of populations reaching their maximum during the acquisition, as well as in the case of population’s continuous grow.
The parametric image of petri dish 2 that depicts the time in which the population is maximized, is shown in Figure 17. Blue color represents pixels with continuously increasing population that did not exhibit maximum value during the acquisition, yellow color marks pixels whose t max cannot be determined because the curve is too noisy to be modeled (or pixels outside the dish) and finally shades of green color (from black to light green) encode the determined value of t m a x for the remaining pixels, in the range [ 0 , 900 ] frame. Note that even the yellow pixels inside the petri dish appear in areas where no visible colony growth could be detected
In Figure 18, four locations are shown on different colonies of petri dish 3. The extracted population curves and the extracted LF and LFSD-based models are depicted in Figure 19. It can be observed that the LFSD achieves more accurate population modeling.

4. Conclusions

In this paper, a low-complexity approach for petri dish image acquisition, automatic measure and mathematically modeling fungus evolution was presented. The images were acquired using an inexpensive web camera. An image-based heuristic method for microbial number estimation was described and applied for the generation of experimental population curves. A novel modeling approach based on the LFSD was presented for the population evolution and compared with the conventional LF. Our results highlighted that the the proposed approach outperforms the commonly-used LF one, in terms of accuracy, since it can also capture decreasing trends of the population evolution. The concept of parametric imaging was studied in order to further improve the visualization of the LFSD model parameters. Finally, the different origins of the three petri dishes, the variety of image acquisition process, image quality and resolution and the different species of growing microorganisms, show the robustness of the proposed population measurement and modeling method.
Future work may include comparison between the growth measures identified in this work, such as the rate of growth and the time of maximum population, between different cultures of the same or different micro-organisms. Such comparisons also require identical image acquisition conditions and equipment and possibly color histogram normalization of the images, adapted to the background color of each individual dish.

Author Contributions

Conceptualization, K.D. and I.B.; methodology, K.D. and I.B.; software, K.D. and I.B.; validation, K.D., I.B. and A.-A.A.B.; formal analysis, A.-A.A.B.; data curation, I.B.; writing–original draft preparation, I.B., K.D.; writing–review and editing, A.-A.A.B.; visualization, A.-A.A.B.; supervision, K.D. All authors have read and agree to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sargent, J.R.; Bartley, C.; Dille, P.S.; Keller, J.M.; Nourbakhsh, I.; LeGrand, R. Timelapse GigaPan: Capturing, Sharing, and Exploring Timelapse Gigapixel Imagery. In Proceedings of the Fine International Conference on Gigapixel Imaging for Science, Pittsburgh, PA, USA, 11–13 November 2010; pp. 1–12. [Google Scholar]
  2. Gonzales, D.S.; Pinheiro, J.C.; Bavister, B.D. Prediction of the developmental potential of hamster embryos in vitro by precise timing of the third cell cycle. Reproduction 1995, 105, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. McKiernan, S.H.; Bavister, B.D. Environmental variables influencing in vitro development of hamster 2-cell embryos to the blastocyst stage. Biol. Reprod. 1990, 43, 404–413. [Google Scholar] [CrossRef] [PubMed]
  4. Schenk, D.; Bacher, S. Functional response of a generalist insect predator to one of its prey species in the field. J. Anim. Ecol. 2002, 71, 524–531. [Google Scholar] [CrossRef]
  5. Hurst, M.R.H.; van Koten, C.; Jackson, T.A. Pathology of Yersinia entomophaga MH96 towards Costelytra zealandica (Coleoptera; Scarabaeidae) larvae. J. Invertebr. Pathol. 2014, 115, 102–107. [Google Scholar] [CrossRef]
  6. Cobo, M.P.; Libro, S.; Marechal, N.; D’Entremont, D.; Cobo, D.P.; Berkmen, M. Visualizing bacterial colony morphologies using timelapse imaging chamber MOCHA. J. Bacteriol. 2018, 200, e00413-17. [Google Scholar]
  7. Dresbøll, D.B.; Thorup-Kristensen, K.; McKenzie, B.M.; Dupuy, L.X.; Bengough, A.G. Timelapse scanning reveals spatial variation in tomato (Solanum lycopersicum L.) root elongation rates during partial waterlogging. Plant Soil 2013, 369, 467–477. [Google Scholar] [CrossRef] [Green Version]
  8. Desai, N.; Ploskonka, S.; Goodman, L.R.; Austin, C.; Goldberg, J.; Falcone, T. Analysis of embryo morphokinetics, multinucleation and cleavage anomalies using continuous timelapse monitoring in blastocyst transfer cycles. Reprod. Biol. Endocrinol. 2014, 12, 54. [Google Scholar] [CrossRef] [PubMed]
  9. Shapiro, J.A.; Hsu, C. Escherichia coli K-12 cell-cell interactions seen by timelapse video. J. Bacteriol. 1989, 171, 5963–5974. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Delpretti, S.; Luisier, F.; Ramani, S.; Blu, T.; Unser, M. Multiframe sure-let denoising of timelapse fluorescence microscopy images. In Proceedings of the 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro, Paris, France, 14–17 May 2008. [Google Scholar] [CrossRef] [Green Version]
  11. Zhao, L.; Chen, Y.; Schaffner, D.W. Comparison of Logistic Regression and Linear Regression in Modeling Percentage Data. Appl. Environ. Microbiol. 2001, 67, 2129–2135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Sartor, M.A.; Leikauf, G.D.; Medvedovic, M. LRpath: A logistic regression approach for identifying enriched biological groups in gene expression data. Bioinformatics 2008, 25, 211–217. [Google Scholar] [CrossRef] [PubMed]
  13. Meyer, P. Bi-logistic growth. Technol. Forecas. Soc. Chang. 1994, 47, 89–102. [Google Scholar] [CrossRef]
  14. Sheehy, J.E. Bi-Phasic Growth Patterns in Rice. Ann. Bot. 2004, 94, 811–817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Meyer, P.S.; Ausubel, J.H. Carrying Capacity: A Model with Logistically Varying Limits. Technol. Forecast. Soc. Chang. 1999, 61, 209–214. [Google Scholar] [CrossRef]
  16. Wong, C.Y.; Goh, K.L. Catch-up models of science and technology: A theorization of the Asian experience from bi-logistic growth trajectories. Technol. Forecast. Soc. Chang. 2015, 95, 312–327. [Google Scholar] [CrossRef]
  17. Tsoularis, A.; Wallace, J. Analysis of logistic growth models. Math. Biosci. 2002, 179, 21–55. [Google Scholar] [CrossRef] [Green Version]
  18. Norton, L. A Gompertzian model of human breast cancer growth. Cancer Res. 1988, 48, 7067–7071. [Google Scholar] [PubMed]
  19. Achkar, J.M.; Fries, B.C. Candida Infections of the Genitourinary Tract. Clin. Microbiol. Rev. 2010, 23, 253–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Bauters, T.G.; Dhont, M.A.; Temmerman, M.I.; Nelis, H.J. Prevalence of vulvovaginal candidiasis and susceptibility to fluconazole in women. Am. J. Obstet. Gynecol. 2002, 187, 569–574. [Google Scholar] [CrossRef] [PubMed]
  21. Morrell, M.; Fraser, V.J.; Kollef, M.H. Delaying the Empiric Treatment of Candida Bloodstream Infection until Positive Blood Culture Results Are Obtained: A Potential Risk Factor for Hospital Mortality. Antimicrob. Agents Chemother. 2005, 49, 3640–3645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ahmad, A.; Khan, A.U. Prevalence of Candida species and potential risk factors for vulvovaginal candidiasis in Aligarh, India. Eur. J. Obstetr. Gynecol. Reprod. Biol. 2009, 144, 68–71. [Google Scholar] [CrossRef] [PubMed]
  23. Van Eck, W.; Lamers, M.H. Biological Content Generation: Evolving Game Terrains Through Living Organisms. In Evolutionary and Biologically Inspired Music, Sound, Art and Design; Springer International Publishing: Berlin, Germany, 2015; pp. 224–235. [Google Scholar] [CrossRef]
Figure 1. Experimental setup.
Figure 1. Experimental setup.
Sensors 20 02545 g001
Figure 2. Two petri dish 1 frames in different microbial evolution states.
Figure 2. Two petri dish 1 frames in different microbial evolution states.
Sensors 20 02545 g002
Figure 3. The evolution of total number of microbes in the whole plate, according to Equation (3).
Figure 3. The evolution of total number of microbes in the whole plate, according to Equation (3).
Sensors 20 02545 g003
Figure 4. (a) The estimated number of microbes according to Equation (9) for each pixel, encoded in shades of gray. (b) A curve example and its smoothing effect at a random point of plate by using median filter.
Figure 4. (a) The estimated number of microbes according to Equation (9) for each pixel, encoded in shades of gray. (b) A curve example and its smoothing effect at a random point of plate by using median filter.
Sensors 20 02545 g004
Figure 5. The experimental curve of the fungi on the whole plate is illustrated in blue and the modeling of the logistic function (LF) with linear and square exponent with pink and green, respectively.
Figure 5. The experimental curve of the fungi on the whole plate is illustrated in blue and the modeling of the logistic function (LF) with linear and square exponent with pink and green, respectively.
Sensors 20 02545 g005
Figure 6. The parametric image resulting from the determination of the linear exponent b 2 of the LF with second degree (LFSD), for each non-zero pixel of the image in the final frame, using the region of interest (ROI) with r = 4 pixels.
Figure 6. The parametric image resulting from the determination of the linear exponent b 2 of the LF with second degree (LFSD), for each non-zero pixel of the image in the final frame, using the region of interest (ROI) with r = 4 pixels.
Sensors 20 02545 g006
Figure 7. (a) The last snapshot of the petri dish, along with two selected areas containing many colonies, depicted in yellow and blue; (b,c) the selected areas of the plate magnified; (d,e) the corresponding experimental population curves of the blue and yellow areas, respectively.
Figure 7. (a) The last snapshot of the petri dish, along with two selected areas containing many colonies, depicted in yellow and blue; (b,c) the selected areas of the plate magnified; (d,e) the corresponding experimental population curves of the blue and yellow areas, respectively.
Sensors 20 02545 g007
Figure 8. The modeling of the generated experimental curves with LF and LFSD for a different size area r, around the central point. (a) r = 15, (b) r = 14, (c) r = 13, (d) r = 12, (e) r = 11, (f) r = 10, (g) r = 9, (h) r = 8, (i) r = 7.
Figure 8. The modeling of the generated experimental curves with LF and LFSD for a different size area r, around the central point. (a) r = 15, (b) r = 14, (c) r = 13, (d) r = 12, (e) r = 11, (f) r = 10, (g) r = 9, (h) r = 8, (i) r = 7.
Sensors 20 02545 g008
Figure 9. The instantaneous signal-to-noise ratio (SNR) for the population curve in Figure 8, calculated for ROI size r = 5 and r = 12 .
Figure 9. The instantaneous signal-to-noise ratio (SNR) for the population curve in Figure 8, calculated for ROI size r = 5 and r = 12 .
Sensors 20 02545 g009
Figure 10. (a) An example of a population curve with negative b1 of the LF. (b) The logarithmic quotient and its classic LF modeling.
Figure 10. (a) An example of a population curve with negative b1 of the LF. (b) The logarithmic quotient and its classic LF modeling.
Sensors 20 02545 g010
Figure 11. The value of b 1 coefficient of the LF as a function of the number of pixels in ROI.
Figure 11. The value of b 1 coefficient of the LF as a function of the number of pixels in ROI.
Sensors 20 02545 g011
Figure 12. Parametric image of the time of first microbial appearance t s .
Figure 12. Parametric image of the time of first microbial appearance t s .
Sensors 20 02545 g012
Figure 13. The pixels with observed zero growth rate during the experiment, shown in green, superimposed on the final frame of the petri dish.
Figure 13. The pixels with observed zero growth rate during the experiment, shown in green, superimposed on the final frame of the petri dish.
Sensors 20 02545 g013
Figure 14. The effect of r to the determination of b2 (a,c) and b1 (b,d) for the LFSD population model are shown for r = 4, 8.
Figure 14. The effect of r to the determination of b2 (a,c) and b1 (b,d) for the LFSD population model are shown for r = 4, 8.
Sensors 20 02545 g014
Figure 15. The final frame of petri dish 2, with four locations marked on different colonies, from which population curves are extracted and modeled.
Figure 15. The final frame of petri dish 2, with four locations marked on different colonies, from which population curves are extracted and modeled.
Sensors 20 02545 g015
Figure 16. The modeling of the generated experimental curves with LF and LFSD for area of size r = 5, around the central point. The curves are extracted from locations marked as (a) A, (b) B, (c) C, and (d) D on the final frame of the petri dish 2.
Figure 16. The modeling of the generated experimental curves with LF and LFSD for area of size r = 5, around the central point. The curves are extracted from locations marked as (a) A, (b) B, (c) C, and (d) D on the final frame of the petri dish 2.
Sensors 20 02545 g016
Figure 17. The final frame of petri dish 2, with four locations marked on different colonies, from which population curves are extracted and modeled.
Figure 17. The final frame of petri dish 2, with four locations marked on different colonies, from which population curves are extracted and modeled.
Sensors 20 02545 g017
Figure 18. The final frame of petri dish 3, with four locations marked on different colonies.
Figure 18. The final frame of petri dish 3, with four locations marked on different colonies.
Sensors 20 02545 g018
Figure 19. The modeling of the generated experimental curves with LF and LFSD for area of size r = 5, around the central point. The curves are extracted from locations marked as (a) A, (b) B, (c) C, and (d) D on the final frame of the petri dish 2.
Figure 19. The modeling of the generated experimental curves with LF and LFSD for area of size r = 5, around the central point. The curves are extracted from locations marked as (a) A, (b) B, (c) C, and (d) D on the final frame of the petri dish 2.
Sensors 20 02545 g019
Table 1. The values of b 1 of LF for different values of r, for the selected pixel of the dish.
Table 1. The values of b 1 of LF for different values of r, for the selected pixel of the dish.
Size of Are (Pixel)Number of Pixels in the Area b 1 ( × 10 3 ) b 0
7225 8.39 1.918
8289 8.44 1.962
9361 8.46 1.811
10441 8.61 2.040
11529 8.51 1.976
12625 8.68 2.004
13729 8.20 1.828
14841 8.68 2.263
15961 8.72 2.386
Table 2. The values of b 2 , b 1 , b 0 of LFSD for different values of r, for the selected pixel of the dish.
Table 2. The values of b 2 , b 1 , b 0 of LFSD for different values of r, for the selected pixel of the dish.
Size of Are (Pixel)Number of Pixels in the Area b 2 ( × 10 7 ) b 1 ( × 10 3 ) b 0
7225 1.06 1.59 2.934
8289 0.78 1.35 2.647
9361 1.01 1.57 2.786
10441 1.26 1.78 3.316
11529 1.25 1.77 3.246
12625 1.23 1.78 3.273
13729 1.31 1.79 3.182
14841 1.60 2.11 4.052
15961 1.78 2.27 4.431
Table 3. The effect of r to the determination of b 1 for the conventional LF population model.
Table 3. The effect of r to the determination of b 1 for the conventional LF population model.
rNumber of Pixels with b 1 < 0
2376
4229
885

Share and Cite

MDPI and ACS Style

Delibasis, K.; Basanou, I.; Boulogeorgos, A.-A.A. Measurement and Modeling of Microbial Growth Using Timelapse Video. Sensors 2020, 20, 2545. https://doi.org/10.3390/s20092545

AMA Style

Delibasis K, Basanou I, Boulogeorgos A-AA. Measurement and Modeling of Microbial Growth Using Timelapse Video. Sensors. 2020; 20(9):2545. https://doi.org/10.3390/s20092545

Chicago/Turabian Style

Delibasis, Konstantinos, Ifigenia Basanou, and Alexandros-Apostolos A. Boulogeorgos. 2020. "Measurement and Modeling of Microbial Growth Using Timelapse Video" Sensors 20, no. 9: 2545. https://doi.org/10.3390/s20092545

APA Style

Delibasis, K., Basanou, I., & Boulogeorgos, A. -A. A. (2020). Measurement and Modeling of Microbial Growth Using Timelapse Video. Sensors, 20(9), 2545. https://doi.org/10.3390/s20092545

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