Next Article in Journal
Antibacterial Effectiveness of Fecal Water and In Vitro Activity of a Multi-Strain Probiotic Formulation against Multi-Drug Resistant Microorganisms
Previous Article in Journal
The Uropathogenic Specific Protein Gene usp from Escherichia coli and Salmonella bongori is a Novel Member of the TyrR and H-NS Regulons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping

1
Software Engineering Institute, East China Normal University, 3663 Zhong Shan Road (N), Shanghai 200062, China
2
School of Life Science, East China Normal University, 3663 Zhong Shan Road (N), Shanghai 200062, China
3
Faculty of Life and Environmental Sciences, University of Tsukuba, 1-1-1 Tennoudai, Tsukuba, Ibaraki 305-8572, Japan
*
Authors to whom correspondence should be addressed.
Microorganisms 2020, 8(3), 331; https://doi.org/10.3390/microorganisms8030331
Submission received: 19 January 2020 / Revised: 18 February 2020 / Accepted: 25 February 2020 / Published: 26 February 2020
(This article belongs to the Section Systems Microbiology)

Abstract

:
Bacterial growth curves, representing population dynamics, are still poorly understood. The growth curves are commonly analyzed by model-based theoretical fitting, which is limited to typical S-shape fittings and does not elucidate the dynamics in their entirety. Thus, whether a certain growth condition results in any particular pattern of growth curve remains unclear. To address this question, up-to-date data mining techniques were applied to bacterial growth analysis for the first time. Dynamic time warping (DTW) and derivative DTW (DDTW) were used to compare the similarity among 1015 growth curves of 28 Escherichia coli strains growing in three different media. In the similarity evaluation, agglomerative hierarchical clustering, assessed with four statistic benchmarks, successfully categorized the growth curves into three clusters, roughly corresponding to the three media. Furthermore, a simple benchmark was newly proposed, providing a highly improved accuracy (~99%) in clustering the growth curves corresponding to the growth media. The biologically reasonable categorization of growth curves suggested that DTW and DDTW are applicable for bacterial growth analysis. The bottom-up clustering results indicate that the growth media determine some specific patterns of population dynamics, regardless of genomic variation, and thus have a higher priority of shaping the growth curves than the genomes do.

1. Introduction

Bacterial population dynamics are known to be affected by growth conditions [1,2,3,4,5,6]; however, whether and how population dynamics are linked to growth conditions remain unclear. Population dynamics are commonly recorded as growth curves, that is, temporal changes in bacterial population density under certain growth conditions. Growth curves have been studied to a large extent due to the simplicity and generality of their common features in a closed habitat as S-shaped curves [7] that are divided into four stages, i.e., the lag, exponential growth, stationary phase, and the death phase [2,8]. Such typical growth curves are often analyzed with various sigmoidal functions, e.g., Logistic, Gompertz, etc. [7,8,9,10], and their derivatives [11,12,13]. Nevertheless, model-based theoretical fitting was only applicable for growth curves with a regular S-shape and not for those with an irregular shape. The shape of a growth curve is unpredictable, even for defined strains growing under well-controlled conditions in the laboratory, because bacterial growth is attributed to not only the genotype but also the external factors, such as the composition of the medium. Our previous findings verified that bacterial growth was either coordinated with genome reduction [14,15] or determined by the chemical composition of the growth medium [16]. Thus, model-based growth analyses were unsuitable to create a linkage between the growth dynamics and the growth condition.
Analyses of the relationship between the growth curve and the growth environment are required for understanding bacterial dynamics in a quantitative manner. To date, no applicable methodology has been reported to classify growth curves in accordance with growth conditions. Here, we assumed that a direct comparison of the shapes of growth curves could provide us with some hints. Whether similar growth conditions resulted in similarly shaped growth curves was investigated. Considering that growth curves are presented in the form of time-series [17], that is, a series of data points indexed temporally, time series clustering might be able to meet the requirement of linking the growth curves to the growth conditions. In contrast to the model-based theoretical fitting, clustering analysis does not need to fit the growth curve to obtain the parameters but does need to divide the growth curves into the different groups according to the similarity of their shapes.
How can the similarity of growth curves be identified? In the present study, we developed a method that categorizes growth curves in response to growth media but not genetic information. Two state-of-the-art data-mining algorithms, dynamic time warping (DTW) [18,19] and derivative dynamic time warping (DDTW) [20,21,22], were first employed in bacterial growth studies. The similarity of the two growth curves was evaluated according to the distance between the two records of the two time series at the same time points. As the growth curves could be different in length and therefore the two records could not be compared, the combination of the time series algorithms of DTW and DDTW was developed. In addition, a novel evaluation criterion was proposed to achieve the proper combination of the number of clusters and the ratio of DTW/DDTW, which led to the best match of the growth curves to the growth media. The successful categorization of the growth curves in response to the growth media without considering the genomic variation strongly suggested that the contribution of growth media to the bacterial growth dynamics shared some similarity independent of genetic information.

2. Materials and Methods

2.1. Bacterial Growth Data

The growth data set was acquired from our previous study [14]. In brief, a total of 28 E. coli strains with genomes of various lengths were cultured in three different media (a complete medium (LB), a minimal medium (M63), and a minimal medium supplemented with 20 amino acids (MAA)), as previously described [14]. The growth curves were acquired by measuring the optical turbidity at an absorbance of 600 nm (OD600) at 30-min or 1-h intervals for 24 to 48 h. A total of 11–30 biological repeats were analyzed for each strain in each medium, which resulted in a total of 1015 growth curves (Figure 1A), that is, 336, 333, and 346 growth curves in the LB, MAA, and M63 media, respectively.

2.2. Growth Data Processing

The temporal changes in bacterial growth recorded by means of optical turbidity were acquired at an interval of 30 min for growth in the LB medium and of 1 h for growth in the M63 and MAA media. As the similarity analysis required the identical temporal interval of the time series for comparison (described in the following section), data thinning for an interval of 1 h was performed with the growth records in LB. In addition, the bacteria replicated by binary fission, so the temporal changes in growth data were on a logarithmic scale. To reflect this biological feature properly, the raw growth data on a linear scale were converted into a logarithmic scale for similarity analysis (Figure 1A, bottom; Figure 1B, step 1).

2.3. Similarity Comparison of the Growth Curves

Two data-mining algorithms, dynamic time warping (DTW), and derivative dynamic time warping (DDTW), were used to measure the similarity between two growth curves (Figure 1B, step 2.1). Briefly, DTW [18] calculated the distance of two growth curves (Equation (1)) and gave an optimal match between them.
DTW ( X ,   Y ) =   Γ ( n ,   m )
Here Γ is the Cost Matrix calculated according to the two time series, the growth curves of X and Y; n and m represent the lengths of X and Y, respectively. DDTW [22] used the first derivative data (X’, Y’) instead of the original time series (X, Y) for calculation, as follows:
X ( i ) = x i + 1 x i ,   ( i = 1 ,   2 ,   3 ,   ,   n 1 ) ,
DDTW ( X ,   Y ) = DTW ( X ,   Y ) ,
According to the calculation algorithms, DTW was largely contributed by the growth lag time and the growth maximum; DDTW was largely attributed to the growth rate (Figure S1). To achieve a reasonable and generally unbiased comparison of the bacterial growth dynamics, we introduced a parameter, α, to combine DTW and DDTW (Equation (4)), as previously proposed [23].
Distance   ( X ,   Y ) = ( 1 α )   DTW ( X ,   Y ) + α   DDTW ( X ,   Y ) ,   α [ 0 ,   1 ] ,
To adjust the weight of DTW and DDTW, the changing interval of α was set to 0.01, because the calculated values (distances) of DTW and DDTW were differed in order. The analyses were performed with Python. The script was implemented originally.

2.4. Hierarchical Clustering of the Growth Curves

The agglomerative algorithm, with the Ward’s linkage and the given number of clusters (N), was used for hierarchical clustering [24] (Figure 1B, step 2.2). The criterion of Ward’s minimum variance method, which calculated the distance between the newly formed cluster and the other clusters, was applied in the hierarchical cluster analysis. The analyses were performed with Python with the package scipy.cluster.hierarchy.

2.5. Evaluation of the Goodness of Clustering

A total of four methods were used for the evaluation. To determine the optimal value of parameter α and cluster number N, we tested one internal benchmark, i.e., the silhouette coefficient (SC), and three external benchmarks, i.e., the adjusted Rand index (ARI), adjusted mutual information (AMI), and the V-measure (V) (Figure 1B, step 2.3). Silhouette (S(i)) interpreted and validated the consistency of the data within clusters [25], where a(i) and b(i) were the intra-cluster and nearest-cluster distances, respectively (Equation (5)).
S ( i ) = { 1 a ( i ) b ( i ) , a ( i ) < b ( i ) 0 , a ( i ) = b ( i ) b ( i ) a ( i ) 1 , a ( i ) > b ( i ) ,
The mean silhouette coefficient (SC) of all calculates of S(i) (Equation (6)) was used for the evaluation of the hierarchical clustering.
SC =   mean ( 0 i = n S ( i ) ) ,
The adjusted Rand index (ARI) was the corrected-for-chance version of the Rand index [26,27], a measure of the similarity between two clusters, where ai, bi, and nij represent the total number of growth curves assigned in each category (medium), the number of growth curves categorized in each cluster (C), and the number of growth curves correctly clustered in corresponding to the assigned category, respectively (Equation (7)).
ARI = i , j ( n i j 2 ) [ i ( a i 2 ) j ( b j 2 ) ] / ( n 2 ) 1 2 [ i ( a i 2 ) + j ( b j 2 ) ] [ i ( a i 2 ) j ( b j 2 ) ] / ( n 2 ) ,
The adjusted mutual information (AMI), similar to ARI, estimated the degree of coincidence between two data distributions [26,27], where C, G, and MI represented the category, the cluster and the mutual information, respectively (Equation (8)).
AMI ( C ,   G ) = M I E [ M I ] max ( H ( C ) ,   H ( G ) ) E [ M I ] ,
The V-measure (V) was the harmonic mean of homogeneity (h) and completeness (c) (Equation (9)). It evaluated whether each cluster contained only a single type of growth curve and whether all growth curves of a given type were assigned to the same cluster [28].
V = 2 h × c h + c ,
The larger values of SC, ARI, AMI, and V indicated that the clusters (the growth curves) were more consistent with the real situation (the growth media). The analyses were performed with Python, with the packages of sklearn.metrics.silhouette_score and sklearn.metrics.

3. Results

3.1. Features of the Growth Curves and Analytical Approaches

Population dynamics of genome-differentiated E. coli strains growing in three different media were acquired from a previous study [14,15]. Since the growth curve of a bacterial population is usually expressed indirectly by measuring the optical turbidity [29,30], the most common method for growth assays, the growth data used as the time series data in the following analyses were the temporal changes in optical density (OD) detected at a wavelength of 600 nm. The temporal changes in the OD600 of E. coli cells growing in LB, M63, and MAA media were recorded as growth curves (Figure 1A). Due to nutritional variation, the growth curves were largely differentiated in response to the growth media, e.g., a rapid increase in OD600 and a short lag time in the rich LB medium (Figure 1A, gray), whereas a gradual increase in OD600 and a long lag time were observed in the poor M63 medium (Figure 1A, blue). These features were thought to be beneficial for clustering the growth curves in correspondence to the growth media, regardless of the analytical algorithms.
On the other hand, a significant variation in the length and shape of the growth curves would largely disturb the clustering accuracy, depending on the analytical approaches. As the different length of growth curves was either attributed to the population dynamics itself or simply because of the recording time of the assay experiments, common normalization of the growth curves would lead to a misrepresented analytical result. To avoid prospective bias in comparing the similarity of the length-differentiated growth curves, we attempted to apply the data-mining algorithms of DTW/DDTW for the first time (Figure 1B). A similarity comparison of the growth curves in the logarithmic scale (Figure 1B, step 1) was performed with the combination of DTW and DDTW (Figure 1B, step 2.1), followed by clustering (Figure 1B, step 2.2). The weight ratios between DTW and DDTW (α) and the number of clusters (N) were determined according to four different statistical benchmarks (Figure 1B, step 2.3). Following the biological assessment of the clustering results, a simple benchmark was newly proposed for improved clustering (Figure 1B, step 3).

3.2. Similarity Evaluation and Hierarchical Clustering

A similarity comparison among growth curves was performed with both DTW and DDTW. Using the combination of the two algorithms of DTW and DDTW, but not any individual, one could reduce the biases/errors in evaluating vertically and/or horizontally differentiated growth curves (Figure S1), which were caused by the experimental operations but not the bacterial growth itself. Hierarchical clustering was performed with the bottom-up approach to avoid bias. The clustering results were determined by both the weight ratio α of DTW/DDTW and the number of clusters N. To determine the optimal combination of α and N, four statistical benchmarks (ARI, AMI, V, and SC) were used for the evaluation. The calculations/simulations of α and N were performed with ranges of 0–1 at a 0.01 interval and 2–28 in units, respectively. The resultant values in the matrix showed the changes in α and N in combination (Figure 2, gray lines). The index for the goodness of clustering was that the larger values represented a higher accuracy. Accordingly, the number of clusters N (Figure 2, red lines) and the corresponding values of α (Figure 2, broken lines) were determined by four benchmarks.
Intriguingly, three statistical benchmarks resulted in a common number of clusters, n = 3, which perfectly agreed with the number of media variations. In particular, ARI and AMI resulted in identical values of α and N as the optimal combination (Figure 2A,B, insets). The number of clusters was equivalent to the variation of the growth media (n = 3) but not to the variation of the genomes (n = 28), strongly indicating that environmental differentiation played a more significant role in bacterial population dynamics than genomic variation.
To verify whether the clustering results were biologically reasonable, how the 1015 growth curves were categorized in the three clusters was validated. The growth curves acquired in the LB and M63 growth media were mostly categorized in a single cluster (C1 or C2), and those acquired in the MAA medium were almost equally divided between the two clusters, according to the benchmark of the SC (Figure 3I). Comparatively, the growth curves acquired in the M63 growth medium were separated into two clusters (C2 and C3), and those acquired in the LB and MAA media were perfectly categorized in the different clusters of C1 and C3, respectively, as evaluated by V-measure (Figure 3II). The differentiated clustering results indicated that the benchmarks of SC and V were good at the varied features (shapes) of the growth curves, although both did not correspond well to the growth media. On the other hand, the growth curves were well categorized into three individual clusters and evaluated by AMI and ARI (Figure 3III,IV). Approximately 90% of the growth curves were clustered according to the growth media, indicating a biologically reasonable categorization. The meaningful clustering revealed that the first attempt to apply DTW/DDTW to bacterial growth analyses was successful.

3.3. A novel Algorithm for Improved Clustering

To improve the accuracy of clustering the growth curves in correspondence to the media, a novel algorithm was proposed as follows: the growth curves categorized in the cluster of the majority were considered as the truth and the rest as errors. For instance, according to the benchmark of AMI and ARI, 300 and 46 growth curves in M63 were categorized in C2 and C3 (Figure 3), which were considered the majority and the errors, respectively. The novel benchmark was that fewer errors represented better clustering. Note that the number of clusters was unnecessarily equal to the media variation, and the growth curves from the same medium could be clustered into multiple clusters.
The similarity evaluation with DTW/DDTW and the hierarchical clustering were performed as described above, resulting in a combination matrix of α and N (Figure 4A, gray). The errors that occurred in all combinations were calculated, and the one that showed the fewest errors was determined to be the best one. Accordingly, the best choice with the fewest errors should be the combination of α = 0.66 and n = 11 (Figure 4A, yellow, arrow). Taking into account the fact that the increase in the number of clusters led to the decrease in the number of errors, whereas the decrease in the rate of errors was gradual (Figure 4B), an alternative choice was made (Figure 4A, red). As the errors decreased most significantly when N turned from 5 to 6, as well as to avoid overclustering, the optimal combination was chosen with 11 errors (Figure 4B, horizontal arrow), instead of the best combination of only three errors (Figure 4B, vertical arrow). Consequently, the combinations of α and N were decided (Figure 4C). As the value of α remained constant and was the highest when N was between 6 and 10, the optimal number of clusters was determined to be six (Figure 4C).
Subsequently, whether the six clusters, with the combination of n = 6 and α = 0.78, achieved a reasonable categorization of the growth curves was examined. The clustering results show that these 11 errors occurred mainly in the MAA medium (Figure 4D). Thus, the newly proposed benchmark, which made only 11 errors, largely improved the accuracy in clustering with biological meaningfulness. Additionally, in the case of the lowest errors, the present clusters were subdivided between the MAA and M63 media (Figure S2). This case was thought to be over clustered or over segmented because the present six clusters were already equivalent to an accuracy of ~99%, corresponding to the growth media.
Intriguingly, the growth curves were categorized almost equally in two clusters, regardless of the media variation (Figure 4D). As the growth curves were acquired from two different initial cell concentrations in all three media, the two clusters might reveal the contribution of the initial population to the shape of the growth curves. Whether there is any significant relationship between the initial concentrations (low or high) and the six clusters was further investigated. Contrary to the assumption, the results show that all the clusters (C1–C6) contained the growth curves of both high and low concentrations (Figure 4E). The results indicate that the bacterial growth curves were not simply decided by any individual parameter (e.g., initial concentration); they were the consequence of the unpredictable and complex interactions among the genomic and environmental factors. Another technical reason was that the changes in initial concentrations usually caused the vertical and/or horizonal shift of growth curves, which was ignored when evaluated with the DTW/DDTW algorithms, as described in the Discussion section (Figure S1).

4. Discussion

In the present study, data-mining algorithms were applied for the first time to investigate bacterial growth curves, which previously were analyzed largely by model-based theoretical fitting [7,8,9,10]. Since the growth curves were not always in a typical S-shape, computational approaches other than theoretical fitting were required to address the questions of whether and how the growth curves were linked to the growth conditions. The growth rate during the exponential phase and the growth maximum at the stationary phase were known to be the important features of the growth curves. Nevertheless, other undiscovered features were also responsible for the shape of the growth curves, e.g., the entire length of the growth curve and the decline in the death phase. The similarity comparison rather than mathematic fitting was used to consider the features of growth as much as possible. Therefore, the similarity comparison (i.e., distance measurement) of the time series, a key part of time series clustering, was employed here. The choice of methodology could greatly influence the clustering results. The commonly used methods include the Euclidean distance [31], Pearson correlation [32], Hausdorff distance [33], DTW [18,19], and DDTW [20,21,22]. As the Euclidean distance and Hausdorff distance were only available for the time series of identical lengths [17], DTW and DDTW were suitable for the analysis, considering the fact that the growth curves were often varied in time scale (length).
The reason for using the combination of DTW and DDTW was that the two algorithms had different weak points. DTW aligned the original X-axis misalignment by dynamically extending and shortening the time series, and DDTW was its first derivative. If the two growth curves were horizontally shifted, both DTW and DDTW produced a high similarity (Figure S1A). However, when they were vertically shifted, even with a similar shape, DTW and DDTW produced conflicting evaluations of low and high similarity, respectively (Figure S1C). If the two growth curves were varied in length with a significant decline (Figure S1B) or showed completely different shapes (Figure S1D), distinctive conclusions would be drawn by DTW and DDTW. Combining DTW with DDTW by introducing the weight ratio α was reasonable to offset some technical errors in the growth curves. As the clustering results were biologically reasonable, the DTW/DDTW method was appropriate for the bacterial growth analyses.
The most valuable point of the present study was that the agglomerative hierarchical clustering relied on the similarity evaluation with DTW/DDTW, resulting in the optimal categorization of 1015 growth curves of 28 different genomes into clusters corresponding to the growth media. Although hierarchical clustering [24] is one of the commonly used clustering algorithms, it was surprising that the final number of clusters was equivalent to the media variation because the clustering was performed statistically without predetermining the number of clusters. In addition, clustering was used to divide the time series data into groups based on similarity, without advanced knowledge of the definitions of the groups [34], and has been applied in the fields of climate, environment, economy, finance, medicine, and biology [35,36,37,38,39]. Most of the applications in biology were related to genes [37,40,41,42,43]. The present study provided the first application for bacterial growth analysis. The success of clustering with biological relevance not only verified the usefulness of the clustering approaches in clarifying bacterial population dynamics but also indicated that the growth media played a deterministic role in the growth dynamics.

5. Conclusions

Temporal changes in bacterial growth analyzed with DTW/DDTW were reported for the first time. More than one thousand growth curves were clustered in a biologically reasonable manner. This result not only demonstrates the advantageousness of up-to-date data-mining algorithms in dynamic studies on living organisms but it also indicates a predictable linkage between bacterial growth patterns and growth conditions. Although the genome was considered the decision-making factor for growth, bottom-up clustering, ignoring genomic differentiation, successfully categorized the growth curves corresponding to the growth media, which strongly suggests that the growth media had a higher priority of shaping the growth curves than the genomes did.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/8/3/331/s1, Figure S1: Comparison of DTW and DDTW; Figure S2: Clusters of the growth curves.

Author Contributions

Analyses, Y.-Y.C.; conceptualization, B.-W.Y.; methodology, Y.-Y.C. and B.-W.Y.; resources, B.-W.Y.; software, Y.-Y.C.; validation, Y.-Y.C., T.Y. and B.-W.Y.; writing, Y.-Y.C. and B.-W.Y.; visualization, Y.-Y.C. and B.-W.Y.; supervision, T.Y. and B.-W.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the JSPS KAKENHI, Grant-in-Aid for Scientific Research (B), grant number 19H03215 (to B.-W.Y.).

Acknowledgments

The authors thank the lab members for technical support and donations in materials and tools used for research.

Conflicts of Interest

The authors declare no conflict of interest. The funder 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. Swinnen, I. Predictive modelling of the microbial lag phase: A review. Int. J. Food Microbiol. 2004, 94, 137–159. [Google Scholar] [CrossRef]
  2. Maier, R.M.; Pepper, I.L. Bacterial Growth. In Environmental Microbiology, 3rd ed.; Academic Press: London, UK, 2015; pp. 37–56. [Google Scholar]
  3. Gibson, A.M.; Bratchell, N.; Roberts, T.A. Predicting microbial growth: Growth responses of salmonellae in a laboratory medium as affected by pH, sodium chloride and storage temperature. Int. J. Food Microbiol. 1988, 6, 155–178. [Google Scholar] [CrossRef]
  4. Baranyi, J.; Roberts, T.A. Mathematics of predictive food microbiology. Int. J. Food Microbiol. 1995, 26, 199–218. [Google Scholar] [CrossRef] [Green Version]
  5. Ponciano, J.M.; Vandecasteele, F.P.J.; Hess, T.F.; Forney, L.J.; Crawford, R.L.; Joyce, P. Use of Stochastic Models To Assess the Effect of Environmental Factors on Microbial Growth. Appl. Environ. Microbiol. 2005, 71, 2355–2364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. An, H.B.; Li, F.; Zhao, M.L.; Liu, Y.J. Optimized Spectral Indices Based Estimation of Forage Grass Biomass. Guang Pu Xue Yu Guang Pu Fen Xi 2015, 35, 3155–3160. [Google Scholar] [PubMed]
  7. 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] [PubMed]
  8. Peleg, M.; Corradini, M.G. Microbial growth curves: What the models tell us and what they cannot. Crit. Rev. Food Sci. Nutr. 2011, 51, 917–945. [Google Scholar] [CrossRef] [PubMed]
  9. 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]
  10. Merks, R.M.H.; Tjørve, K.M.C.; Tjørve, E. The use of Gompertz models in growth analyses, and new Gompertz-model approach: An addition to the Unified-Richards family. PLoS ONE 2017, 12, e0178691. [Google Scholar]
  11. Yates, G.T.; Smotzer, T. On the lag phase and initial decline of microbial growth curves. J. Theor. Biol. 2007, 244, 511–517. [Google Scholar] [CrossRef]
  12. Tsuchiya, K.; Cao, Y.-Y.; Kurokawa, M.; Ashino, K.; Yomo, T.; Ying, B.-W. A decay effect of the growth rate associated with genome reduction in Escherichia coli. BMC Microbiol. 2018, 18, 101. [Google Scholar] [CrossRef] [PubMed]
  13. Fujikawa, H.; Kai, A.; Morozumi, S. A new logistic model for Escherichia coli growth at constant and dynamic temperatures. Food Microbiol. 2004, 21, 501–509. [Google Scholar] [CrossRef]
  14. Kurokawa, M.; Seno, S.; Matsuda, H.; Ying, B.-W. Correlation between genome reduction and bacterial growth. DNA Res. 2016, 23, 517–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kurokawa, M.; Ying, B.W.J.M. Experimental Challenges for Reduced Genomes: The Cell Model Escherichia coli. Microorganisms 2020, 8, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Ashino, K.; Sugano, K.; Amagasa, T.; Ying, B.W. Predicting the decision making chemicals used for bacterial growth. Sci. Rep. 2019, 9, 7251. [Google Scholar] [CrossRef]
  17. Aghabozorgi, S.; Shirkhorshidi, A.S.; Wah, T.Y. Time-series clustering—A decade review. Inf. Syst. 2015, 53, 16–38. [Google Scholar] [CrossRef]
  18. Senin, P.J.I.; Computer Science Department University of Hawaii at Manoa Honolulu, U. Dynamic time warping algorithm review. Int. J. Res. Electr. Comput. Eng. 2008, 855, 40. [Google Scholar]
  19. Sakoe, H.; Chiba, S. Dynamic programming algorithm optimization for spoken word recognition. Inf. Comput. Sci. Dep. Univ. Hawaii Manoa Honol. 1978, 26, 43–49. [Google Scholar] [CrossRef] [Green Version]
  20. Górecki, T.; Łuczak, M. First and Second Derivatives in Time Series Classification Using DTW. Commun. Stat. Simul. Comput. 2014, 43, 2081–2092. [Google Scholar]
  21. Górecki, T.; Łuczak, M. Using derivatives in time series classification. Data Min. Knowl. Discov. 2012, 26, 310–331. [Google Scholar] [CrossRef] [Green Version]
  22. Keogh, E.J.; Pazzani, M.J. Derivative dynamic time warping. In Proceedings of the 2001 SIAM international conference on data mining, Chicago, IL, USA, 5–7 April 2001; pp. 1–11. [Google Scholar]
  23. Łuczak, M. Hierarchical clustering of time series data with parametric derivative dynamic time warping. Expert Syst. Appl. 2016, 62, 116–130. [Google Scholar] [CrossRef]
  24. Kaufman, L.; Rousseeuw, P.J. Finding Groups in Data: An Introduction to Cluster Analysis, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2009; pp. 1–162. [Google Scholar]
  25. Rousseeuw, P.J. Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. J. Comput. Appl. Math. 1987, 20, 53–65. [Google Scholar] [CrossRef] [Green Version]
  26. Vinh, N.X.; Epps, J.; Bailey, J. Information theoretic measures for clusterings comparison: Variants, properties, normalization and correction for chance. J. Mach. Learn. Res. 2010, 11, 2837–2854. [Google Scholar]
  27. Vinh, N.X.; Epps, J.; Bailey, J. Information Theoretic Measures for Clusterings Comparison: Is a Correction for Chance Necessary? In Proceedings of the 26th Annual International Conference on Machine Learning, Montreal, QC, Canada, 14–18 June 2009; pp. 1073–1080. [Google Scholar]
  28. Rosenberg, A.; Hirschberg, J. V-Measure: A Conditional Entropy-Based External Cluster Evaluation Measure. In Proceedings of the 2007 Joint Conference on Empirical Methods in Natural Language Processing and Computational Natural Language Learning (EMNLP-CoNLL), Prague, Czech Republic, 28–30 June 2007; pp. 410–420. [Google Scholar]
  29. Myers, J.A.; Curtis, B.S.; Curtis, W.R. Improving accuracy of cell and chromophore concentration measurements using optical density. BMC Biophys. 2013, 6, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Mytilinaios, I.; Salih, M.; Schofield, H.K.; Lambert, R.J.W. Growth curve prediction from optical density data. Int. J. Food Microbiol. 2012, 154, 169–176. [Google Scholar] [CrossRef] [Green Version]
  31. Deza, M.M.; Deza, E. Encyclopedia of Distances. In Encyclopedia of Distances, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 1–583. [Google Scholar]
  32. Lee, R.J.; Nicewander, W.A. Thirteen Ways to Look at the Correlation Coefficient. Am. Stat. 1988, 42, 59–66. [Google Scholar]
  33. Huttenlocher, D.P.; Klanderman, G.A.; Rucklidge, W.J. Comparing images using the Hausdorff distance. IEEE Trans. Pattern Anal. Mach. Intell. 1993, 15, 850–863. [Google Scholar] [CrossRef] [Green Version]
  34. Rai, P.; Singh, S. A survey of clustering techniques. Int. J. Comput. Appl. 2010, 7, 1–5. [Google Scholar] [CrossRef]
  35. Aghabozorgi, S.; Teh, Y.W. Stock market co-movement assessment using a three-phase clustering method. Expert Syst. Appl. 2014, 41, 1301–1314. [Google Scholar] [CrossRef]
  36. Elangasinghe, M.A.; Singhal, N.; Dirks, K.N.; Salmond, J.A.; Samarasinghe, S. Complex time series analysis of PM10 and PM2.5 for a coastal site using artificial neural network modelling and k-means clustering. Atmos. Environ. 2014, 94, 106–116. [Google Scholar] [CrossRef]
  37. Ernst, J.; Nau, G.J.; Bar-Joseph, Z. Clustering short time series gene expression data. Bioinformatics 2005, 21, i159–i168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Steinbach, M.; Tan, P.-N.; Kumar, V.; Klooster, S.; Potter, C. Discovery of climate indices using clustering. Proc. Proc. Ninth Acm Sigkdd Int. Conf. Knowl. Discov. Data Min. 2003, 3, 446–455. [Google Scholar]
  39. Wismüller, A.; Lange, O.; Dersch, D.R.; Leinsinger, G.L.; Hahn, K.; Pütz, B.; Auer, D. Cluster analysis of biomedical image time-series. Int. J. Comput. Vis. 2002, 46, 103–128. [Google Scholar] [CrossRef]
  40. Bishopric, N.H.; Pyatnitskiy, M.; Mazo, I.; Shkrob, M.; Schwartz, E.; Kotelnikova, E. Clustering Gene Expression Regulators: New Approach to Disease Subtyping. PLoS ONE 2014, 9, e84955. [Google Scholar]
  41. Möller-Levet, C.S.; Klawonn, F.; Cho, K.H.; Wolkenhauer, O. Fuzzy clustering of short time-series and unevenly distributed sampling points. In International Symposium on Intelligent Data Analysis; Springer: Berlin/Heidelberg, Germany, 2003; pp. 330–340. [Google Scholar]
  42. Subhani, N.; Rueda, L.; Ngom, A.; Burden, C.J. Multiple gene expression profile alignment for microarray time-series data clustering. Bioinformatics 2010, 26, 2281–2288. [Google Scholar] [CrossRef] [Green Version]
  43. Fujita, A.; Severino, P.; Kojima, K.; Sato, J.R.; Patriota, A.G.; Miyano, S. Functional clustering of time series gene expression data by Granger causality. BMC Syst. Biol. 2012, 6, 137. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Overview of the growth data and analytical approach. (A) Bacterial growth curves used in the present study. A total of 28 E. coli strains were grown in the LB (gray), MAA (yellow) and M63 media (blue). The top and bottom panels represent the growth curves shown in the linear and logarithmic scales, respectively. (B) Scheme of the analytical approach used in the present study. Three steps (Steps 1–3), including data processing, growth analyses, and a novel proposal were used.
Figure 1. Overview of the growth data and analytical approach. (A) Bacterial growth curves used in the present study. A total of 28 E. coli strains were grown in the LB (gray), MAA (yellow) and M63 media (blue). The top and bottom panels represent the growth curves shown in the linear and logarithmic scales, respectively. (B) Scheme of the analytical approach used in the present study. Three steps (Steps 1–3), including data processing, growth analyses, and a novel proposal were used.
Microorganisms 08 00331 g001
Figure 2. Similarity comparison with dynamic time warping (DTW) and derivative DTW (DDTW). The similarity of the growth curves was evaluated with DTW/DDTW, and the clustering was evaluated by four different statistical benchmarks, the adjusted Rand index (ARI) (A), adjusted mutual information (AMI) (B), the V-measure (V) (C), and the silhouette coefficient (SC) (D) Each line represents the changes in the statistical benchmark along with the changes in the weight ratio (α) of DTW/DDTW under a certain number (N) of clusters. A total of 28 lines in gray gradation are shown for all four benchmarks, representing n = 2–28. The red lines represent the optimal N determined with the four statistical benchmarks. The best combinations of α and N decided by the statistical benchmarks are indicated with the broken lines and the insets.
Figure 2. Similarity comparison with dynamic time warping (DTW) and derivative DTW (DDTW). The similarity of the growth curves was evaluated with DTW/DDTW, and the clustering was evaluated by four different statistical benchmarks, the adjusted Rand index (ARI) (A), adjusted mutual information (AMI) (B), the V-measure (V) (C), and the silhouette coefficient (SC) (D) Each line represents the changes in the statistical benchmark along with the changes in the weight ratio (α) of DTW/DDTW under a certain number (N) of clusters. A total of 28 lines in gray gradation are shown for all four benchmarks, representing n = 2–28. The red lines represent the optimal N determined with the four statistical benchmarks. The best combinations of α and N decided by the statistical benchmarks are indicated with the broken lines and the insets.
Microorganisms 08 00331 g002
Figure 3. Clusters of the growth curves. The four benchmarks, SC, V, AMI and ARI are indicated in IIV, respectively, according to the best combinations of α and N indicated in Figure 2. The determined clusters are shown as C1, C2, and C3. The growth curves categorized in different clusters are counted and correspond to the LB, MAA, and M63 growth media. The resultant numbers of the growth curves are indicated with both the values and the color bars.
Figure 3. Clusters of the growth curves. The four benchmarks, SC, V, AMI and ARI are indicated in IIV, respectively, according to the best combinations of α and N indicated in Figure 2. The determined clusters are shown as C1, C2, and C3. The growth curves categorized in different clusters are counted and correspond to the LB, MAA, and M63 growth media. The resultant numbers of the growth curves are indicated with both the values and the color bars.
Microorganisms 08 00331 g003
Figure 4. New benchmark for better clustering. (A). A similarity comparison with DTW/DDTW. The similarity of the growth curves was evaluated with DTW/DDTW, and the clustering was evaluated by the new method. Each line represents the changes in the benchmark along with the changes in the weight ratio (α) of DTW/DDTW under a certain number (N) of clusters. A total of 28 lines in gray gradation are shown, representing n = 2–28. The red and yellow lines represent the optimal N values of 6 and 11, respectively. The best combination of α and N determined by the new benchmark is indicated with the arrow. (B) The relationship between the number of errors and the number of clusters. The decreased number of errors along with the increased number of clusters when α = 0.78 is shown. The horizontal and vertical arrows indicate the positions of n = 6 and 11, respectively. (C) The relationship between the weight ratio of DTW/DDTW and the number of clusters. The changes in α are shown with the changes in N when the number of errors is 11. (D) The categorization of the growth curves. The determined combination of α and N and the number of errors are indicated. The growth curves categorized in different clusters (C1–C6) are counted, corresponding to the LB, MAA, and M63 growth media. The resultant numbers of growth curves are indicated with both the values and the color bars. (E) A breakdown of the six clusters upon the initial cell concentration. The growth curves initiated from low and high concentrations are indicated in light and dark, respectively.
Figure 4. New benchmark for better clustering. (A). A similarity comparison with DTW/DDTW. The similarity of the growth curves was evaluated with DTW/DDTW, and the clustering was evaluated by the new method. Each line represents the changes in the benchmark along with the changes in the weight ratio (α) of DTW/DDTW under a certain number (N) of clusters. A total of 28 lines in gray gradation are shown, representing n = 2–28. The red and yellow lines represent the optimal N values of 6 and 11, respectively. The best combination of α and N determined by the new benchmark is indicated with the arrow. (B) The relationship between the number of errors and the number of clusters. The decreased number of errors along with the increased number of clusters when α = 0.78 is shown. The horizontal and vertical arrows indicate the positions of n = 6 and 11, respectively. (C) The relationship between the weight ratio of DTW/DDTW and the number of clusters. The changes in α are shown with the changes in N when the number of errors is 11. (D) The categorization of the growth curves. The determined combination of α and N and the number of errors are indicated. The growth curves categorized in different clusters (C1–C6) are counted, corresponding to the LB, MAA, and M63 growth media. The resultant numbers of growth curves are indicated with both the values and the color bars. (E) A breakdown of the six clusters upon the initial cell concentration. The growth curves initiated from low and high concentrations are indicated in light and dark, respectively.
Microorganisms 08 00331 g004

Share and Cite

MDPI and ACS Style

Cao, Y.-Y.; Yomo, T.; Ying, B.-W. Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping. Microorganisms 2020, 8, 331. https://doi.org/10.3390/microorganisms8030331

AMA Style

Cao Y-Y, Yomo T, Ying B-W. Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping. Microorganisms. 2020; 8(3):331. https://doi.org/10.3390/microorganisms8030331

Chicago/Turabian Style

Cao, Yang-Yang, Tetsuya Yomo, and Bei-Wen Ying. 2020. "Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping" Microorganisms 8, no. 3: 331. https://doi.org/10.3390/microorganisms8030331

APA Style

Cao, Y. -Y., Yomo, T., & Ying, B. -W. (2020). Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping. Microorganisms, 8(3), 331. https://doi.org/10.3390/microorganisms8030331

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