2  Experiment analysis

2.1 Introduction

The previous chapter (see Chapter 1) detailed the steps necessary to extract data from a set of microfluidic images through image analysis techniques and fluorescence microscopy. Each step was instrumental in creating a dataset that was easy to explore and ask questions. With the help of computational biology, systems biology, and data analysis techniques, we could process these files to help us in the search to find the role of filamentation in cell survival.

Computational biology and systems biology contributed to the development of this analysis. In principle, computational biology originated after the origin of computer science with the British mathematician and logistician Alan Turing (regularly known as the father of computing) (TURING 1950). Over time, systems biology emerged as an area that synergistically combines models and experimental data to understand biological processes (Bruggeman et al. 2007). Thus, giving a step towards creating models that, in general, are phenomenological but sometimes serve to discover new ideas about the process under study. Without the computer’s power, modern ideas and aspects of studying biological sciences would otherwise be unthinkable.

Here, we divide the experimental analysis into two main parts: 1) at the cell level or measurements at specific points in time and 2) at the population level and time series. The first level allowed us to identify the individual contribution of each variable understudy to determine cell survival. The second level allowed us to understand how the population behaves according to the passage of time in the face of exposure to a harmful agent (in this case, beta-lactam antibiotics). Together, both visions of the same study phenomenon allowed us to extract the main ideas for postulating a mathematical model that seeks to show how filamentation is a factor for cell survival in stressful environments (see Chapter 3).

2.2 General preprocessing of data

The raw data processing consisted mainly of creating two levels of observation for the cells of both chromosomal strains and multicopy plasmids. The first level is at a cell granularity, point properties. The second level consists of the cells over time, thus observing properties at the population level. We did this because it would allow us to understand what factors affect filamentation and why.

We normalized the fluorescence values of DsRed and GFP for both experiments based on the values observed before antibiotic exposure. It allowed us to have a basis to work with and compare expressions between cells. In the case of DsRed environment drug concentration, we also applied a logarithmic transformation to observe subtle changes in fluorescence intensity that would allow us to detect cell death.

Ultimately, we decided to classify cells into four fundamental groups based on whether the cell filamented and survived (see Figure 2.1). We define a filamented cell as a cell with more than two standard deviations from the mean concerning the lengths observed before introducing antibiotics into the system. On the other hand, although there are multiple ways to define death from single-cell observations (Trevors 2012; Kroemer et al. 2008), we considered a cell dead or missing when we stopped having information about it, either because of fluorescence in the red channel was above a given threshold (resulting from an increase in cell membrane permeability and the introduction of fluorescent dye into the cell) or because it left the field of observation. Therefore, we defined a surviving cell as a cell observed before and after antibiotic exposure that did not surpass the DsRed death threshold.

Figure 2.1: Cell classification and its distribution across experiments. We define a filamented cell as a cell whose length exceeded two standard deviations from the mean at any time during the experiment. A surviving cell is a cell we observed before and after exposure to the antibiotic and did not surpass the DsRed death threshold. Accordingly, we removed from the analysis those cells that died before or were born after antibiotic exposure. Therefore, we delimited the effect caused by antibiotic exposure.

2.3 Results

2.3.1 Cell length and the amount of GFP are crucial in determining cell survival

We evaluated the DsRed, GFP, and length values for each cell at different time points: initial, filamentation, and end. This preprocessing allowed us to observe and quantify each cell at critical times in the experiment and eliminate noise or signals outside the scope of this investigation.

We define the initial time as the first time we observed the cell in the experiment. Filamentation time equals when a cell reaches the filamentation threshold (see Figure 2.4) for the first time. We defined the end time as the time of the last observation of the cell. We decided to bound the end time for surviving cells to one frame (10 min) after the end of antibiotic exposure so that the observed signal would reflect the final stress responses.

When we compared the distributions of DsRed, GFP, and length for both experiments, we observed the changes in their role in cell survival. In Figure 2.2, we show that indistinctly and, as expected, surviving cells managed to eliminate the antibiotic by the end time. In contrast, dead cells presented higher levels of antibiotics (measured by proxy through the mean DsRed intensity of the cell).

Figure 2.2: DsRed temporal distribution. To evaluate the incident effect of the antibiotic marked by DsRed on cells by class, we show its values at three key moments: start, filamentation (SOS), and end. The upper asterisks represent the significance value when comparing a group X to the filamented and surviving cell reference. Asterisks in a line indicate whether or not there is a significant difference in the survival of non-filamented cells. Dots represent the mean of each group. The line bars represents the distribution of the data. Although, at the initial time, we observe multiple significant differences, this is likely due to the intrinsic noise of the system since, as expected, the values are close to zero. We observed a difference between the surviving and non-filamented cells for the chromosomal strain for the SOS time, but the same did not occur for the plasmid strain. The final amount of DsRed makes a clear difference between survival and death.

On the other hand, GFP observations in Figure 2.3 showed us that filamented cells had low fluorescent intensities (low plasmid copy-number) at the beginning of the experiment. In comparison, the chromosomal strain did not exhibit noticeable changes in GFP levels. For the final observation times, GFP measurements indicated that among the cells that did not filament, the ones that survived exhibited a reduced GFP expression concerning cells killed by the antibiotic. Meanwhile, for the filamented cells, whether surviving or dead, their GFP measurements indicated no difference at the beginning or the end of the experiment, suggesting the presence of other determinants of cell survival.

Figure 2.3: GFP temporal distribution. To evaluate the incident effect of the GFP on cells by class, we used the same notation as in Figure 2.2. The chromosomal strain exhibits variability in GFP at different time points, mainly due to experimental noise resulting from low fluorescent intensity values. As expected, filamented cells had a lower initial GFP in the plasmid strain. At the time of filamentation, there appear to be differences in fluorescence between surviving and dead cells. However, in the end time, we observed that the surviving non-filamented cells have lower GFP values than the non-filamented dead cells and alive filamented cells.

Cell length was one of the factors that GFP expression levels could not explain for cell survival. In Figure 2.4, we show that the conclusions regarding filamentation were applicable for both chromosomal and plasmid strains. For the initial times, filamented and survived cells were shorter in length than those that died but longer than not filamented cells of both classes, while non-filamented cells did not differ. We observed no length differences between cells at filamentation time. Thus, survival could depend on other factors, such as growth rate. In the final time, the results were well-defined. Surviving cells had a greater length than their non-surviving pair (i.e., dead filamented and non-filamented cells). However, for filamented cells, surviving cells generally represent a distribution of higher final length values but are not as extensive as their dead counterpart. Which we could explain as a length limit to which cells can grow without dying. Nevertheless, we had no information to evaluate such a hypothesis.

Figure 2.4: Length temporal distribution. To evaluate the incident effect of length on cells by class, we use the same notation as in Figure 2.2. The observations for both strains, chromosomal or plasmid, are the same. In the beginning, the surviving filamented cells already have a difference in length from the rest of the classes. At the time of filamentation, there is no difference to help determine whether the cell will survive or not. Finally, in the final time, it seems that the surviving filamented cells have a greater length than the rest of the groups. However, this length is moderate compared to the excess length shown by non-surviving filamented cells. On the other hand, we highlighted the growth of the surviving non-filamented cells. Therefore, although they did not reach a length for us to classify as filamented, the cells did resort to filamentation.

Once we observed the effects of GFP expression levels and lengths in determining whether a cell lives or dies, we projected the cells onto the plane. We painted them with their class status (See Figure 2.1) to determine whether these two variables contained the necessary information to cluster the data correctly. In Figure 2.5, we show the initial GFP and length values projection. While, with some work, we could contextually place the results in Figure 2.3 and Figure 2.4, the initial values did not appear to determine the classes. Therefore, we explored the final versus initial values differences in Figure 2.6. With this new representation of the cells in the plane, we contextualized the statistical results presented in Figure 2.3 and Figure 2.4. Besides, it showed us that differences in length (i.e., filamentation) and reductions in GFP expression are essential in determining cell survival. Though the clustering of cell status is not entirely separated, other variables affect the experimental results in cell survival.

Figure 2.5: Experiment’s initial values. By positioning a cell in space based on its initial length and GFP values, we can see that class separation occurs, but not as a strong signal. Therefore, we concluded that although the initial state influences the result, this is not everything. For this, we have the example of the length changes throughout the experiment caused by filamentation. In this graph, the GFP scale is at log10 to help us observe those minor differences between the experiments.

Figure 2.6: Experiment’s initial values differences. By comparing the metric differences of the last observation and the first observation of a cell, we can separate mainly the surviving filamented cells from those that did not do it in both experiments (blue dots). Meanwhile, cells with plasmids form a small accumulation of surviving cells that did not produce filament (green dots). However, this has made a breakthrough in understanding what is affecting cell survival. There are still variables that we can include to understand this phenomenon better.

2.3.2 Number of divisions and cell age do not appear to play a clear role in determining cell survival

In Section 2.3.1, we explored the effect on cell survival through GFP variability and cell length.However, Figure 2.5 and Figure 2.6 showed us the possibility of other factors relevant to the phenomenon under study. As some papers in the literature suggest, some of these other factors may be cell division and chronological age (i.e., how much time has passed since the last cell division at the time of exposure to a toxic agent) (Moger-Reischer and Lennon 2019; Roostalu et al. 2008; Heinrich, Leslie, and Jonas 2015). Therefore, we chose to observe these two metrics in experiments at a purely qualitative level, i.e., without the inclusion of, e.g., metrics of membrane or cell cycle properties (Joseleau-Petit, Vinella, and D’Ari 1999).

Although we expected to see a small contribution, either by the number of divisions or cell age, in Figure 2.7 and Figure 2.8, we could not observe a precise effect of these variables on cell survival. Although they could have an explanation or biological significance, we decided to omit as relevant in the characterization of our cells, since the signal was not clear. However, we derived from this analysis a slightly simpler variable that tells us whether a cell underwent a cell division event or not. So it gives us a more generalized picture of the contribution of division to cell survival (see Figure 2.14).

Figure 2.7: Cell’s number of divisions. Both chromosomal and plasmid cells exhibited a wider distribution of divisions for the surviving cells against non-surviving cells. However, we did not appreciate a significant change between the chromosome filamented cells. Therefore, the number of cell divisions’ contribution to filamentation remains uncertain.

Figure 2.8: Time elapsed since the last division at the beginning of the experiment. The mean time of the last division before starting the experiment indicates that it did not influence the final result for chromosomal cells. There is a slight difference between the filamented-not survived cells and the rest for cells with plasmids. However, the signal does not appear to be strong on the survival role. Therefore, we conclude that we have no evidence to support that the time of the last division at the beginning of the experiment influences the final classification results.

2.3.3 Time to reach filamentation matters in determining cell survival

In Figure 2.2, Figure 2.3, and Figure 2.4, we showed how, at the time of filamentation, DsRed and GFP levels appeared indifferent to the cells. Therefore, we hypothesized that a possible variable determining cell survival could be its time to activate its anti-stress response system that causes filamentation. Furthermore, we also guided our hypothesis by previous reports showing us how the gene expression level can induce filamentation with tight temporal coordination.

While, for our analyses, we did not measure the concentration of antibiotic that triggers filamentation per se, we indirectly quantified its effect by using the time it took for a cell to reach a length at which it is already considered a filamentating cell. Furthermore, to recognize that the observed effect was a product of the experiment, we decided to keep only filamented cells just once antibiotic exposure began.

Figure 2.9 shows how filamentation times are narrower for chromosomal cells than for plasmid-bearing cells. Then, we hypothesize that the effect could come from the heterogeneity in the plasmid copy number in the population. Also, interestingly, we observed that, for both experiments, cells that survived had longer filamentation times than those that died. These differences in response times suggest the following: 1) if the cell grows too fast, it will reach a limit and start to accumulate antibiotics constantly, and 2) if the cell grows too fast, likely, the cost of maintaining an ample length for prolonged periods of exposure will become counterproductive.

Figure 2.9: Time to filamentation filtered. We only keep cells that filamented during the antibiotic exposure to quantify their time to filamentation and its effect on survival. In this way, we normalize the start times for the calculation of the filamentation time. For both strains, the filamentation time had a more significant delay in the surviving cells.

In Figure Figure 2.10, we decided to project the results of Figure Figure 2.9 in a space similar to the one described in Figure Figure 2.5). Thus, we separated our data into cells that survived and cells that did not and painted them when it took them to reach their filamented state. We realized that, by adding this temporal component to the initial variables of length and GFP, we could separate surviving cells from dead cells to a greater degree. However, it may still not be enough, and there are still many other variables that play a crucial role in understanding the ecology of stress and how some cells will be survivors or not.

Figure 2.10: Experiment initial values with time to filamentation. As in Figure 2.5, including the time it will take for cells to filament allows us to better understand the phenomenon of survival. Cells that filamented and survived generally have a much higher delay than their non-filamented peers for both strains (see Figure 2.9).

2.3.4 Increasing the system’s complexity and analyzing it in an unsupervised way allows a correct classification of cell states

In the experiments, we observed the importance of GFP filamentation and variability for cell survival. Similarly, we realized that other variables must be affecting the final results. Filamentation and GFP variability alone did not fully recapitulate the expected behavior of the data. That is, the target variables did not capture the system’s heterogeneity.

The inability to reproduce cell classification led us to question two things: 1) the possibility that our sorting was wrong beforehand, and 2) we did not have enough variables to capture the study phenomenon. We decided to take the unsupervised learning way to answer these subjects because it allows us to project our data without prior knowledge.

We opted for the path of dimensionality reduction techniques where each variable or feature is equivalent to one dimension. The essence of dimensionality reduction is that it is not feasible to analyze each dimension with many dimensions. Furthermore, dimensionality reduction helps us counteract several problems, such as reducing the complexity of a model, reducing the possibility of overfitting a model, removing all correlated variables, and visualizing our data in a two- or three-dimensional space for better appreciation. Improved visualization and identification of essential variables are the main reasons to guide and complement our research with this technique.

2.3.4.1 Principal Component Analysis (PCA) emphasizes the importance of cell length and its GFP in cell survival

The first dimensionality reduction technique we decided to use was Principal Component Analysis (PCA) (Pearson 1901; Hotelling 1936). Scientist mainly uses PCA to create predictive models or in Exploratory Data Analysis (EDA). In our case, we only use it as an EDA.

For chromosomal and plasmid strain, in Figure 2.11 and Figure 2.12, we show the projection of the first two principal components (PCs), respectively. Figure 2.11 separates the manually annotated classes, surviving cells separated from non-surviving cells. However, for Figure 2.12, the class separation was a bit rougher but allowed us to separate the surviving filament cells from the dead ones.

Figure 2.11: Principal Component Analysis of chromosomal strain. When integrating the information of different variables in a dimensionality reduction analysis, we observed a clear separation between the surviving cells and those that did not. The contributions that determined this phenomenon come mainly from the last amount of DsRed, GFP, and cell length (see Figure 2.13). Although it seems obvious, it effectively confirms that the temporal classification that we carry out makes sense. Longer length represents a greater uptake of antibiotics, but in a much larger volume, so the net effect is an internal reduction of antibiotics (see Figure 3.1).

Figure 2.12: Principal Component Analysis of plasmid strain. By integrating the information from different variables in a dimensionality reduction analysis, we observed a clear separation between the filamented and non-filamented cells. Said class separation is given by component 2 (Y-axis), which is determined primarily by the initial and final lengths of the cells (see Figure 2.14). Furthermore, the classification also allows us to separate those filamented cells that died from those that survived. Therefore, despite the increase in the system’s complexity, length plays a role in determining survival.

For their part, in Figure 2.13 and Figure 2.14, we show the total contribution of each variable per PC for the chromosomal and plasmid strain, respectively. Finding that filamentation plays a crucial role in determining cell survival. For example, for PC2, we appreciated how the variable end DsRed directed the dots to the positive side, while the variable end and start length directed the dots to the opposing side. Therefore, we can support that filamentation has a role in moving cells away from having higher amounts of DsRed.

Figure 2.13: Variables contribution of Principal Component Analysis of chromosomal strain. In Figure 2.11, we see that the classes we created manually reflected what we observed when performing a reduction of dimensions analysis. Here we show the individual contribution of each variable for the first two components. The variables that most affected components 1 and 2 (X-axis and Y-axis, respectively) are the final measurements of DsRed, GFP, length, and the initial amount of GFP. Given that they are chromosomal strains, we should note that this variability could be produced by intrinsic experimental noise that we could not remove. With that in mind, having the DsRed and the final length highlights the inherent role of cells by having increased their size.

Figure 2.14: Variables contribution of Principal Component Analysis of plasmid strain. In Figure 2.12, we saw that we could separate the filamented cells from the non-filamented ones. The reduction analysis also shows a slight difference between surviving and dead cells within the small group of filamented cells. Here we offer the individual contribution of each variable for the first two components. For the first component (x-axis in Figure 2.11), the initial and final GFP measurements mainly received the component’s variability. We expected this component’s importance since it is a chromosomal strain, so we hope its inherent variation will be inherited. On the other hand, the second component (Y-axis in Figure 2.11) was determined by the length of the cell. Factors that, in the chromosomal strain (see Figure 2.13), determined with the help of DsRed the separation between surviving and dead cells.

2.3.4.2 Uniform Manifold Approximation and Projection (UMAP) correctly represents the local structure of cell states

Staying with only a one-dimensionality reduction technique was not an option, so we used the UMAP technique (McInnes, Healy, and Melville 2018). We mainly decided to use UMAP for clustering purposes and see if the annotated clusters corresponded to the manually annotated ones. UMAP has certain advantages for these purposes, e.g., it preserves the global structure across the whole space, so the distances between clusters matter.

In Figure 2.15 and Figure 2.16, we show how, using the same variables used in the “PCA” section, UMAP accomplished clustering the four proposed classes correctly. Interestingly, in Figure 2.15, UMAP formed three general groups and four for Figure 2.16. However, in general, UMAP clustered the surviving cells from those that did not survive. On investigating why this separation occurred, we found that the large groups coalesced into one another if we eliminated the division variable. So, in a way, the division also has a role in determining survival, but it is not essential or at least not over-represented in our data.

Figure 2.15: UMAP coordinates of chromosome strain. We represented the cells in a low-dimensional space. This new projection allowed it to group the cells that survived and those that did not. Therefore, as in PCA Figure 2.11, this technique supports the manual classification that we carry out.

Figure 2.16: UMAP coordinates of plasmid strain. As in Figure 2.15, the representation in a low-dimensional space helped classify the cells into four groups, two survivors and two non-survivors. The variable division marks the separation of classes. The division variable indicates whether a cell divided during its lifetime or not. Together, the UMAP represents the manually assigned classes.

2.3.5 Population dynamics reveal how filamentation contributes cell survival

From the full tracking dataset, we evaluated how the different cell states behaved over time—for example, understanding how the cells absorbed antibiotics or how they elongated in time. In contrast to the dataset generated in the Section 2.3.1, we did not truncate the results 10 minutes after the antibiotic exposure. In this way, we could observe cell behavior before and after the presence of the toxic agent.

In Figure 2.17, we observed a small fraction of filamentous cells without exposure to the toxic agent in both cell strains. However, after antibiotic exposure at minute 60, we observed increases in the proportion of filamented cells. It is interesting to note how filamented cells grew after antibiotic exposure for the chromosomal strain. We speculate that this post-antibiotic growth exists because, once the SOS system that triggers filamentation is activated, the system continues to grow until it reaches a limit regardless of whether the damaging agent is still present (Justice et al. 2008; Mückl et al. 2018). Moreover, we observed how the cells start to divide again after some time because the proportion of non-filament cells starts to grow while the filament cells start to divide. We observed the same effects for the plasmid strain. However, the number of filament cells expected was much lower by experimental design.

Figure 2.17: Population status over time. We calculate how many cells of each type existed for each time point: non-filamented and filamented living cells (green and orange areas, respectively) and dead cells (red area; we considered dead cells as those that existed at one time and then stopped tracking). The gray vertical lines represent each experiment’s start and end of antibiotic exposure. The experiment was finalized with the resolution of the cells when they returned to their non-filamented state. The effect of filamentation and its spread after exposure to the antibiotic is evident for the chromosomal strain. For its part, for the plasmid strain, it is observed how the filamented cells begin to appear slowly. Their proportion is as expected, given that the population had a wide distribution of GFP that allowed them to combat exposure to the antibiotic.

In Figure 2.18, we showed that once antibiotics exposure began, those cells that died had a much faster increase in DsRed than those that did manage to live, regardless of whether they were filamented. On the other hand, surviving cells maintained their relatively stable DsRed levels. We noted that length was critical for the surviving cells for the chromosomal strain by turning to the GFP and length variables for a temporal explanation. Even cells categorized as non-filamented reached the filamentation threshold minutes after antibiotic exposure. However, the distinction between live or dead filamented cells was not as evident as expected. As for cells with plasmids, the effect on GFP for surviving cells was maintained for filamented cells and decreased for non-filamented cells. For the filament cells that died, we showed that they had, on average, a much longer initial length than the surviving cells. We also consider it necessary to understand which variables affect cell survival.

Figure 2.18: Population measurements over time. The colored lines symbolize the average value of each metric at each instant of time, while its surrounding shaded area represents the 95% confidence interval. The vertical lines represent the start and end of antibiotic exposure. The horizontal line in the length metric symbolizes the threshold to consider a cell filament. Regarding the GFP metric, the behavior is relatively stable for the chromosomal strain. We observed a faster increase of DsRed for the non-surviving populations in both experiments. In contrast, for the plasmid strain, a decline in GFP is observed for the population that did not survive. For the length metric, it is interesting to note how the chromosome cells that did not filament continued to grow past the filamentation threshold once the exposure to the antibiotic in the chromosomal strain had ended. On the other hand, the filamented and dead cells seem to have a greater length from the beginning for the plasmid strain.

2.3.6 Heterogeneity in plasmid copy-number allows various forms of survival in addition to filamentation

We are confident that filamentation has a fundamental role in determining cell survival, as we have shown so far. However, for plasmid cells, we have a component of our complete interest; heterogeneity. Each cell can possess a different plasmid copy number; thus, each could show a different behavior under stress (San Millan et al. 2016). For instance, heterogeneity can produce resistant cells that do not suffer damage, susceptible cells, and cells that form filaments to mitigate environmental stress.

To study the effect of variability in plasmid copy number on the survival probability of the population, we decided to group cells by the proportion of initial GFP with respect to the population maximum. We defined 100% of the population as the number of total cells at the onset of antibiotic exposure. Figure 2.19 shows how the cells with the highest amount of GFP remained unchanged once antibiotic exposure began, while the rest of the cells started to decrease their percentage of surviving cells. However, the decrease was not linear. On the contrary, we observed a bi-modal distribution in the reduction of live cells. An average GFP point provided higher survival than a point below or above the average (except for cells very close to the population maximum).

Figure 2.19: Population survivals binned by initial GFP over time. We categorized the cells’ GFP into ranges of proportions 0.05 concerning the maximum amount of GFP in the population. 100% cells per bin of GFP was taken as the number of cells one frame before the start of exposure to the antibiotic (minute 50). Therefore, dark to light colors represent a generation of new cells, and light to dark colors the death of cells. The gray vertical bars represent the start and end of antibiotic exposure. Bar’s size and color on the right represent the percentage of the living cells 10 minutes after the end of the experiment. As shown in Figure 2.20, we showed that the surviving cells appear to follow something similar to a bimodal distribution. More cells survive with a moderate amount of GFP or with an amount close to the maximum of the population.

Therefore, what we observed was a bimodal distribution for GFP-dependent cell survival. In order to show this effect more clearly, in Figure 2.20, we plotted the survival probability for each GFP bin without normalizing for the population maximum. This new plot allowed us to observe how the bimodal survival distribution occurs for cells that did not grow as filaments, whereas cells that filament increase their survival probability gradually as they have more initial GFP (see also Figure 2.3).-distribution).

Figure 2.20: Plasmid initial GFP survival probability. We calculated the survival probability after comparing the population distributions of GFP with those of the cells that managed to survive. To assess survival by GFP, we only used plasmid cells. A bell forms with an upturned tail for non-filamented cells (green dots). On the other hand, for the filamented cells (orange dots), a continuous increase in survival is shown just when it seems that the probability of the non-filamented cells has decreased. In global, much GFP has higher resistance, but an average GFP value without filamentation also increases the probability of survival.

As in Figure 2.20, in Figure 2.21, we show the survival probability given an initial length. We observe that survival is higher for cells that did not grow as filaments if the initial length was less than the average. In contrast, for filamented cells, the survival probability increased as cell length was longer at the beginning of the experiment (see also Figure 2.4). However, it is noteworthy that the probability of survival had a limit in which a higher initial length meant a lower probability of survival (see red dotted lines in Figure 2.21).

Figure 2.21: Plasmid initial length survival probability. We calculated the survival probability after comparing the population distributions of length with those of the cells that managed to survive. For non-filamented cells (blue dots), the survival probability is higher for those cells with small initial lengths, while It seems to decrease with a more extensive initial size. For their part, for filamented cells (red dots), the probability of survival increases according to their length but then declines when the cells are too long at first (see red dotted line). Therefore, generally, a small and moderate length or an initial length already filamented from the beginning increases the chances of survival.

2.4 Discussion

Here, we evaluated different variables that could determine cell survival upon exposure to toxic agents by studying two experimental populations of E. coli, one strain with a resistance gene on the chromosome and the other on multicopy plasmids. We identified two variables that are predominantly responsible for cell survival: cell length and GFP amount related to the cell’s inherent resistance to the toxic agent and heterogeneity in response times.

On the other hand, as other studies have already mentioned (Heinrich, Leslie, and Jonas 2015; Wang and Levin 2009), we examined cell activity and youth in a minimalistic way. While the distribution of the number of divisions exemplifies a broader and more uniform range for the surviving cells, the cells that died tended to have fewer divisions. However, for the study of cellular youth at the time of exposure to the toxic agent, the results did not show a clear pattern of behavior for cell fate determination. Therefore, it would be interesting to study cellular youth at a higher level of complexity in future studies to understand its contribution to cell survival.

Interestingly, when we used temporal measurements of cell length, GFP, DsRed, and if a cell divided, we could recapitulate, for the most part, the fates of cellular states (see Section 2.3.1 and Section 2.3.4). Thus, increasing the system’s complexity led to better clustering of cell states, but not how these factors interact biologically in determining cell survival. Therefore, we decided to postulate a mathematical model that helps us understand the critical components of cell survival.