Transcription Attenuation in Synthetic Promoters in Nonoverlapping Tandem Formation.

Closely spaced promoters are ubiquitous in prokaryotic and eukaryotic genomes. How their structure and dynamics relate remains unclear, particularly for tandem formations. To study their transcriptional interference, we engineered two pairs and one trio of synthetic promoters in nonoverlapping, tandem formation, in single-copy plasmids transformed into Escherichia coli cells. From in vivo measurements, we found that these promoters in tandem formation can have attenuated transcription rates. The attenuation strength can be widely fine-tuned by the promoters' positioning, natural regulatory mechanisms, and other factors, including the antibiotic rifampicin, which is known to hamper RNAP promoter escape. From this, and supported by in silico models, we concluded that the attenuation in these constructs emerges from premature terminations generated by collisions between RNAPs elongating from upstream promoters and RNAPs occupying downstream promoters. Moreover, we found that these collisions can cause one or both RNAPs to falloff. Finally, the broad spectrum of possible, externally regulated, attenuation strengths observed in our synthetic tandem promoters suggests that they could become useful as externally controllable regulators of future synthetic circuits.

Alteration of DNA supercoiling serves as a trigger of short-term cold shock repressed genes of E. coli.

Cold shock adaptability is a key survival skill of gut bacteria of warm-blooded animals. Escherichia coli cold shock responses are controlled by a complex multi-gene, timely-ordered transcriptional program. We investigated its underlying mechanisms. Having identified short-term, cold shock repressed genes, we show that their responsiveness is unrelated to their transcription factors or global regulators, while their single-cell protein numbers' variability increases after cold shock. We hypothesized that some cold shock repressed genes could be triggered by high propensity for transcription locking due to changes in DNA supercoiling (likely due to DNA relaxation caused by an overall reduction in negative supercoiling). Concomitantly, we found that nearly half of cold shock repressed genes are also highly responsive to gyrase inhibition (albeit most genes responsive to gyrase inhibition are not cold shock responsive). Further, their response strengths to cold shock and gyrase inhibition correlate. Meanwhile, under cold shock, nucleoid density increases, and gyrases and nucleoid become more colocalized. Moreover, the cellular energy decreases, which may hinder positive supercoils resolution. Overall, we conclude that sensitivity to diminished negative supercoiling is a core feature of E. coli's short-term, cold shock transcriptional program, and could be used to regulate the temperature sensitivity of synthetic circuits.

Analytical kinetic model of native tandem promoters in E. coli.

Closely spaced promoters in tandem formation are abundant in bacteria. We investigated the evolutionary conservation, biological functions, and the RNA and single-cell protein expression of genes regulated by tandem promoters in E. coli. We also studied the sequence (distance between transcription start sites 'dTSS', pause sequences, and distances from oriC) and potential influence of the input transcription factors of these promoters. From this, we propose an analytical model of gene expression based on measured expression dynamics, where RNAP-promoter occupancy times and dTSS are the key regulators of transcription interference due to TSS occlusion by RNAP at one of the promoters (when dTSS ≤ 35 bp) and RNAP occupancy of the downstream promoter (when dTSS > 35 bp). Occlusion and downstream promoter occupancy are modeled as linear functions of occupancy time, while the influence of dTSS is implemented by a continuous step function, fit to in vivo data on mean single-cell protein numbers of 30 natural genes controlled by tandem promoters. The best-fitting step is at 35 bp, matching the length of DNA occupied by RNAP in the open complex formation. This model accurately predicts the squared coefficient of variation and skewness of the natural single-cell protein numbers as a function of dTSS. Additional predictions suggest that promoters in tandem formation can cover a wide range of transcription dynamics within realistic intervals of parameter values. By accurately capturing the dynamics of these promoters, this model can be helpful to predict the dynamics of new promoters and contribute to the expansion of the repertoire of expression dynamics available to synthetic genetic constructs.

Efficiency and Robustness of Processes Driven by Nucleoid Exclusion in Escherichia coli.

The internal spatial organization of prokaryotic organisms, including Escherichia coli, is essential for the proper functioning of processes such as cell division. One source of this organization in E. coli is the nucleoid, which causes the exclusion of macromolecules - e.g. protein aggregates and the chemotaxis network - from midcell. Similarly, following DNA replication, the nucleoid(s) assist in placing the Z-ring at midcell. These processes need to be efficient in optimal conditions and robust to suboptimal conditions. After reviewing recent findings on these topics, we make use of past data to study the efficiency of the spatial constraining of Z-rings, chemotaxis networks, and protein aggregates, as a function of the nucleoid(s) morphology. Also, we compare the robustness of these processes to nonoptimal temperatures. We show that Z-rings, Tsr clusters, and protein aggregates have temperature-dependent spatial distributions along the major cell axis that are consistent with the nucleoid(s) morphology and the volume-exclusion phenomenon. Surprisingly, the consequences of the changes in nucleoid size with temperature are most visible in the kurtosis of these spatial distributions, in that it has a statistically significant linear correlation with the mean nucleoid length and, in the case of Z-rings, with the distance between nucleoids prior to cell division. Interestingly, we also find a negative, statistically significant linear correlation between the efficiency of these processes at the optimal condition and their robustness to suboptimal conditions, suggesting a trade-off between these traits.

SCIP: a single-cell image processor toolbox.

Summary: Each cell is a phenotypically unique individual that is influenced by internal and external processes, operating in parallel. To characterize the dynamics of cellular processes one needs to observe many individual cells from multiple points of view and over time, so as to identify commonalities and variability. With this aim, we engineered a software, 'SCIP', to analyze multi-modal, multi-process, time-lapse microscopy morphological and functional images. SCIP is capable of automatic and/or manually corrected segmentation of cells and lineages, automatic alignment of different microscopy channels, as well as detect, count and characterize fluorescent spots (such as RNA tagged by MS2-GFP), nucleoids, Z rings, Min system, inclusion bodies, undefined structures, etc. The results can be exported into *mat files and all results can be jointly analyzed, to allow studying not only each feature and process individually, but also find potential relationships. While we exemplify its use on Escherichia coli, many of its functionalities are expected to be of use in analyzing other prokaryotes and eukaryotic cells as well. We expect SCIP to facilitate the finding of relationships between cellular processes, from small-scale (e.g. gene expression) to large-scale (e.g. cell division), in single cells and cell lineages. Availability and implementation: http://www.ca3-uninova.org/project_scip. Supplementary information: Supplementary data are available at Bioinformatics online.

Reversible Disruption of Neuronal Mitochondria by Ischemic and Traumatic Injury Revealed by Quantitative Two-Photon Imaging in the Neocortex of Anesthetized Mice.

Mitochondria play a variety of functional roles in cortical neurons, from metabolic support and neuroprotection to the release of cytokines that trigger apoptosis. In dendrites, mitochondrial structure is closely linked to their function, and fragmentation (fission) of the normally elongated mitochondria indicates loss of their function under pathological conditions, such as stroke and brain trauma. Using in vivo two-photon microscopy in mouse brain, we quantified mitochondrial fragmentation in a full spectrum of cortical injuries, ranging from severe to mild. Severe global ischemic injury was induced by bilateral common carotid artery occlusion, whereas severe focal stroke injury was induced by Rose Bengal photosensitization. The moderate and mild traumatic injury was inflicted by focal laser lesion and by mild photo-damage, respectively. Dendritic and mitochondrial structural changes were tracked longitudinally using transgenic mice expressing fluorescent proteins localized either in cytosol or in mitochondrial matrix. In response to severe injury, mitochondrial fragmentation developed in parallel with dendritic damage signified by dendritic beading. Reconstruction from serial section electron microscopy confirmed mitochondrial fragmentation. Unlike dendritic beading, fragmentation spread beyond the injury core in focal stroke and focal laser lesion models. In moderate and mild injury, mitochondrial fragmentation was reversible with full recovery of structural integrity after 1-2 weeks. The transient fragmentation observed in the mild photo-damage model was associated with changes in dendritic spine density without any signs of dendritic damage. Our findings indicate that alterations in neuronal mitochondria structure are very sensitive to the tissue damage and can be reversible in ischemic and traumatic injuries. SIGNIFICANCE STATEMENT: During ischemic stroke or brain trauma, mitochondria can either protect neurons by supplying ATP and adsorbing excessive Ca2+, or kill neurons by releasing proapoptotic factors. Mitochondrial function is tightly linked to their morphology: healthy mitochondria are thin and long; dysfunctional mitochondria are thick (swollen) and short (fragmented). To date, fragmentation of mitochondria was studied either in dissociated cultured neurons or in brain slices, but not in the intact living brain. Using real-time in vivo two-photon microscopy, we quantified mitochondrial fragmentation during acute pathological conditions that mimic severe, moderate, and mild brain injury. We demonstrated that alterations in neuronal mitochondria structural integrity can be reversible in traumatic and ischemic injuries, highlighting mitochondria as a potential target for therapeutic interventions.

Temperature-dependence of the single-cell variability in the kinetics of transcription activation in Escherichia coli.

From in vivo single-cell, single-RNA measurements of the activation times and subsequent steady-state active transcription kinetics of a single-copy Lac-ara-1 promoter in Escherichia coli, we characterize the intake kinetics of the inducer (IPTG) from the media, following temperature shifts. For this, for temperature shifts of various degrees, we obtain the distributions of transcription activation times as well as the distributions of intervals between consecutive RNA productions following activation in individual cells. We then propose a novel methodology that makes use of deconvolution techniques to extract the mean and the variability of the distribution of intake times. We find that cells, following shifts to low temperatures, have higher intake times, although, counter-intuitively, the cell-to-cell variability of these times is lower. We validate the results using a new methodology for direct estimation of mean intake times from measurements of activation times at various inducer concentrations. The results confirm that E. coli's inducer intake times from the environment are significantly higher following a shift to a sub-optimal temperature. Finally, we provide evidence that this is likely due to the emergence of additional rate-limiting steps in the intake process at low temperatures, explaining the reduced cell-to-cell variability in intake times.

The precision of the symmetry in Z-ring placement in Escherichia coli is hampered at critical temperatures.

Cell division in Escherichia coli is morphologically symmetric due to, among other things, the ability of these cells to place the Z-ring at midcell. Studies have reported that, at sub-optimal temperatures, this symmetry decreases at the single-cell level, but the causes remain unclear. Using fluorescence microscopy, we observe FtsZ-GFP and DAPI-stained nucleoids to assess the robustness of the symmetry of Z-ring formation and positioning in individual cells under sub-optimal and critical temperatures. We find the Z-ring formation and positioning to be robust at sub-optimal temperatures, as the Z-ring's mean width, density and displacement from midcell maintain similar levels of correlation to one another as at optimal temperatures. However, at critical temperatures, the Z-ring displacement from midcell is greatly increased. We present evidence showing that this is due to enhanced distance between the replicated nucleoids and, thus, reduced Z-ring density, which explains the weaker precision in setting a morphologically symmetric division site. This also occurs in rich media and is cumulative, i.e. combining richer media and critically high temperatures enhances the asymmetries in division, which is evidence that the causes are biophysical. To further support this, we show that the effects are reversible, i.e. shifting cells from optimal to critical, and then to optimal again, reduces and then enhances the symmetry in Z-ring positioning, respectively, as the width and density of the Z-ring return to normal values. Overall, our findings show that the Z-ring positioning in E. coli is a robust biophysical process under sub-optimal temperatures, and that critical temperatures cause significant asymmetries in division.

Effects of σ factor competition are promoter initiation kinetics dependent.

In Escherichia coli, the expression of a σ factor is expected to indirectly down-regulate the expression of genes recognized by another σ factor, due to σ factor competition for a limited pool of RNA polymerase core enzymes. Evidence suggests that the sensitivity of genes to indirect down-regulation differs widely. We studied the variability in this sensitivity in promoters primarily recognized by RNAP holoenzymes carrying σ(70). From qPCR and live single-cell, single-RNA measurements of the transcription kinetics of several σ(70)-dependent promoters in various conditions and from the analysis of σ factors population-dependent models of transcription initiation, we find that, the smaller is the time-scale of the closed complex formation relative to the open complex formation, the weaker is a promoter's responsiveness to changes in σ(38) numbers. We conclude that, in E. coli, a promoter's responsiveness to indirect regulation by σ factor competition is determined by the sequence-dependent kinetics of the rate limiting steps of transcription initiation.

Increased cytoplasm viscosity hampers aggregate polar segregation in Escherichia coli.

In Escherichia coli, under optimal conditions, protein aggregates associated with cellular aging are excluded from midcell by the nucleoid. We study the functionality of this process under sub-optimal temperatures from population and time lapse images of individual cells and aggregates and nucleoids within. We show that, as temperature decreases, aggregates become homogeneously distributed and uncorrelated with nucleoid size and location. We present evidence that this is due to increased cytoplasm viscosity, which weakens the anisotropy in aggregate displacements at the nucleoid borders that is responsible for their preference for polar localisation. Next, we show that in plasmolysed cells, which have increased cytoplasm viscosity, aggregates are also not preferentially located at the poles. Finally, we show that the inability of cells with increased viscosity to exclude aggregates from midcell results in enhanced aggregate concentration in between the nucleoids in cells close to dividing. This weakens the asymmetries in aggregate numbers between sister cells of subsequent generations required for rejuvenating cell lineages. We conclude that the process of exclusion of protein aggregates from midcell is not immune to stress conditions affecting the cytoplasm viscosity. The findings contribute to our understanding of E. coli's internal organisation and functioning, and its fragility to stressful conditions.

Characterizing rate limiting steps in transcription from RNA production times in live cells.

MOTIVATION: Single-molecule measurements of live Escherichia coli transcription dynamics suggest that this process ranges from sub- to super-Poissonian, depending on the conditions and on the promoter. For its accurate quantification, we propose a model that accommodates all these settings, and statistical methods to estimate the model parameters and to select the relevant components. RESULTS: The new methodology has improved accuracy and avoids overestimating the transcription rate due to finite measurement time, by exploiting unobserved data and by accounting for the effects of discrete sampling. First, we use Monte Carlo simulations of models based on measurements to show that the methods are reliable and offer substantial improvements over previous methods. Next, we apply the methods on measurements of transcription intervals of different promoters in live E. coli, and show that they produce significantly different results, both in low- and high-noise settings, and that, in the latter case, they even lead to qualitatively different results. Finally, we demonstrate that the methods can be generalized for other similar purposes, such as for estimating gene activation kinetics. In this case, the new methods allow quantifying the inducer uptake dynamics as opposed to just comparing them between cases, which was not previously possible. We expect this new methodology to be a valuable tool for functional analysis of cellular processes using single-molecule or single-event microscopy measurements in live cells. AVAILABILITY AND IMPLEMENTATION: Source code is available under Mozilla Public License at http://www.cs.tut.fi/%7Ehakkin22/censored/ CONTACT: andre.ribeiro@tut.fi or andre.sanchesribeiro@tut.fi SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Temperature-Dependent Model of Multi-step Transcription Initiation in Escherichia coli Based on Live Single-Cell Measurements.

Transcription kinetics is limited by its initiation steps, which differ between promoters and with intra- and extracellular conditions. Regulation of these steps allows tuning both the rate and stochasticity of RNA production. We used time-lapse, single-RNA microscopy measurements in live Escherichia coli to study how the rate-limiting steps in initiation of the Plac/ara-1 promoter change with temperature and induction scheme. For this, we compared detailed stochastic models fit to the empirical data in maximum likelihood sense using statistical methods. Using this analysis, we found that temperature affects the rate limiting steps unequally, as nonlinear changes in the closed complex formation suffice to explain the differences in transcription dynamics between conditions. Meanwhile, a similar analysis of the PtetA promoter revealed that it has a different rate limiting step configuration, with temperature regulating different steps. Finally, we used the derived models to explore a possible cause for why the identified steps are preferred as the main cause for behavior modifications with temperature: we find that transcription dynamics is either insensitive or responds reciprocally to changes in the other steps. Our results suggests that different promoters employ different rate limiting step patterns that control not only their rate and variability, but also their sensitivity to environmental changes.

Polar Localization of the Serine Chemoreceptor of Escherichia coli Is Nucleoid Exclusion-Dependent.

We studied whether nucleoid exclusion contributes to the segregation and retention of Tsr chemoreceptor clusters at the cell poles. Using live time-lapse, single-cell microscopy measurements, we show that the single-cell spatial distributions of Tsr clusters have heterogeneities and asymmetries that are consistent with nucleoid exclusion and cannot be explained by the diffusion-and-capture mechanism supported by Tol-Pal complexes at the poles. Also, in cells subjected to ampicillin, which enhances relative nucleoid lengths, Tsr clusters locate relatively closer to the cell extremities, whereas in anucleated cells (deletion mutants for mukB), the Tsr clusters are closer to midcell. In addition, we find that the fraction of Tsr clusters at the poles is smaller in deletion mutants for Tol-Pal than in wild-type cells, although it is still larger than would be expected by chance. Also in deletion mutants, the distribution of Tsr clusters differs widely between cells with relatively small and large nucleoids, in a manner consistent with nucleoid exclusion from midcell. This comparison further showed that diffusion-and-capture by Tol-Pal complexes and nucleoid exclusion from the midcell have complementary effects. Subsequently, we subjected deletion mutants to suboptimal temperatures that are known to enhance cytoplasm viscosity, which hampers nucleoid exclusion effects. As the temperature was lowered, the fraction of clusters at the poles decreased linearly. Finally, a stochastic model including nucleoid exclusion at midcell and diffusion-and-capture due to Tol-Pal at the poles is shown to exhibit a cluster dynamics that is consistent with the empirical data. We conclude that nucleoid exclusion also contributes to the preference of Tsr clusters for polar localization.

Robustness of the Process of Nucleoid Exclusion of Protein Aggregates in Escherichia coli.

UNLABELLED: Escherichia coli segregates protein aggregates to the poles by nucleoid exclusion. Combined with cell divisions, this generates heterogeneous aggregate distributions in subsequent cell generations. We studied the robustness of this process with differing medium richness and antibiotics stress, which affect nucleoid size, using multimodal, time-lapse microscopy of live cells expressing both a fluorescently tagged chaperone (IbpA), which identifies in vivo the location of aggregates, and HupA-mCherry, a fluorescent variant of a nucleoid-associated protein. We find that the relative sizes of the nucleoid's major and minor axes change widely, in a positively correlated fashion, with medium richness and antibiotic stress. The aggregate's distribution along the major cell axis also changes between conditions and in agreement with the nucleoid exclusion phenomenon. Consequently, the fraction of aggregates at the midcell region prior to cell division differs between conditions, which will affect the degree of asymmetries in the partitioning of aggregates between cells of future generations. Finally, from the location of the peak of anisotropy in the aggregate displacement distribution, the nucleoid relative size, and the spatiotemporal aggregate distribution, we find that the exclusion of detectable aggregates from midcell is most pronounced in cells with mid-sized nucleoids, which are most common under optimal conditions. We conclude that the aggregate management mechanisms of E. coli are significantly robust but are not immune to stresses due to the tangible effect that these have on nucleoid size. IMPORTANCE: Escherichia coli segregates protein aggregates to the poles by nucleoid exclusion. From live single-cell microscopy studies of the robustness of this process to various stresses known to affect nucleoid size, we find that nucleoid size and aggregate preferential locations change concordantly between conditions. Also, the degree of influence of the nucleoid on aggregate positioning differs between conditions, causing aggregate numbers at midcell to differ in cell division events, which will affect the degree of asymmetries in the partitioning of aggregates between cells of future generations. Finally, we find that aggregate segregation to the cell poles is most pronounced in cells with mid-sized nucleoids. We conclude that the energy-free process of the midcell exclusion of aggregates partially loses effectiveness under stressful conditions.

Estimation of GFP-tagged RNA numbers from temporal fluorescence intensity data.

MOTIVATION: MS2-GFP-tagging of RNA is currently the only method to measure intervals between consecutive transcription events in live cells. For this, new transcripts must be accurately detected from intensity time traces. RESULTS: We present a novel method for automatically estimating RNA numbers and production intervals from temporal data of cell fluorescence intensities that reduces uncertainty by exploiting temporal information. We also derive a robust variant, more resistant to outliers caused e.g. by RNAs moving out of focus. Using Monte Carlo simulations, we show that the quantification of RNA numbers and production intervals is generally improved compared with previous methods. Finally, we analyze data from live Escherichia coli and show statistically significant differences to previous methods. The new methods can be used to quantify numbers and production intervals of any fluorescent probes, which are present in low copy numbers, are brighter than the cell background and degrade slowly. AVAILABILITY: Source code is available under Mozilla Public License at http://www.cs.tut.fi/%7ehakkin22/jumpdet/.

Estimation of fluorescence-tagged RNA numbers from spot intensities.

MOTIVATION: Present research on gene expression using live cell imaging and fluorescent proteins or tagged RNA requires accurate automated methods of quantification of these molecules from the images. Here, we propose a novel automated method for classifying pixel intensities of fluorescent spots to RNA numbers. RESULTS: The method relies on a new model of intensity distributions of tagged RNAs, for which we estimated parameter values in maximum likelihood sense from measurement data, and constructed a maximum a posteriori classifier to estimate RNA numbers in fluorescent RNA spots. We applied the method to estimate the number of tagged RNAs in individual live Escherichia coli cells containing a gene coding for an RNA with MS2-GFP binding sites. We tested the method using two constructs, coding for either 96 or 48 binding sites, and obtained similar distributions of RNA numbers, showing that the method is adaptive. We further show that the results agree with a method that uses time series data and with quantitative polymerase chain reaction measurements. Lastly, using simulated data, we show that the method is accurate in realistic parameter ranges. This method should, in general, be applicable to live single-cell measurements of low-copy number fluorescence-tagged molecules. AVAILABILITY AND IMPLEMENTATION: MATLAB extensions written in C for parameter estimation and finding decision boundaries are available under Mozilla public license at http://www.cs.tut.fi/%7ehakkin22/estrna/ CONTACT: andre.ribeiro@tut.fi.

In silico analysis of division times of Escherichia coli populations as a function of the partitioning scheme of non-functional proteins.

Recent evidence suggests that cells employ functionally asymmetric partitioning schemes in division to cope with aging. We explore various schemes in silico, with a stochastic model of Escherichia coli that includes gene expression, non-functional proteins generation, aggregation and polar retention, and molecule partitioning in division. The model is implemented in SGNS2, which allows stochastic, multi-delayed reactions within hierarchical, transient, interlinked compartments. After setting parameter values of non-functional proteins' generation and effects that reproduce realistic intracellular and population dynamics, we investigate how the spatial organization of non-functional proteins affects mean division times of cell populations in lineages and, thus, mean cell numbers over time. We find that division times decrease for increasingly asymmetric partitioning. Also, increasing the clustering of non-functional proteins decreases division times. Increasing the bias in polar segregation further decreases division times, particularly if the bias favors the older pole and aggregates' polar retention is robust. Finally, we show that the non-energy consuming retention of inherited non-functional proteins at the older pole via nucleoid occlusion is a source of functional asymmetries and, thus, is advantageous. Our results suggest that the mechanisms of intracellular organization of non-functional proteins, including clustering and polar retention, affect the vitality of E. coli populations.

In vivo single-molecule kinetics of activation and subsequent activity of the arabinose promoter.

Using a single-RNA detection technique in live Escherichia coli cells, we measure, for each cell, the waiting time for the production of the first RNA under the control of PBAD promoter after induction by arabinose, and subsequent intervals between transcription events. We find that the kinetics of the arabinose intake system affect mean and diversity in RNA numbers, long after induction. We observed the same effect on Plac/ara-1 promoter, which is inducible by arabinose or by IPTG. Importantly, the distribution of waiting times of Plac/ara-1 is indistinguishable from that of PBAD, if and only if induced by arabinose alone. Finally, RNA production under the control of PBAD is found to be a sub-Poissonian process. We conclude that inducer-dependent waiting times affect mean and cell-to-cell diversity in RNA numbers long after induction, suggesting that intake mechanisms have non-negligible effects on the phenotypic diversity of cell populations in natural, fluctuating environments.

In vivo kinetics of segregation and polar retention of MS2-GFP-RNA complexes in Escherichia coli.

The cytoplasm of Escherichia coli is a crowded, heterogeneous environment. From single cell live imaging, we investigated the spatial kinetics and heterogeneities of synthetic RNA-protein complexes. First, although their known tendency to accumulate at the cell poles does not appear to introduce asymmetries between older and newer cell poles within a cell lifetime, these emerge with cell divisions. This suggests strong polar retention of the complexes, which we verified in their history of positions and mean escape time from the poles. Next, we show that the polar retention relies on anisotropies in the displacement distribution in the region between midcell and poles, whereas the speed is homogeneous along the major cell axis. Afterward, we establish that these regions are at the border of the nucleoid and shift outward with cell growth, due to the nucleoid's replication. Overall, the spatiotemporal kinetics of the complexes, which is robust to suboptimal temperatures, suggests that nucleoid occlusion is a source of dynamic heterogeneities of macromolecules in E. coli that ultimately generate phenotypic differences between sister cells.

CellAging: a tool to study segregation and partitioning in division in cell lineages of Escherichia coli.

MOTIVATION: Cell division in Escherichia coli is morphologically symmetric. However, as unwanted protein aggregates are segregated to the cell poles and, after divisions, accumulate at older poles, generate asymmetries in sister cells' vitality. Novel single-molecule detection techniques allow observing aging-related processes in vivo, over multiple generations, informing on the underlying mechanisms. RESULTS: CellAging is a tool to automatically extract information on polar segregation and partitioning in division of aggregates in E.coli, and on cellular vitality. From time-lapse, parallel brightfield and fluorescence microscopy images, it performs cell segmentation, alignment of brightfield and fluorescence images, lineage construction and pole age determination, and it computes aging-related features. We exemplify its use by analyzing spatial distributions of fluorescent protein aggregates from images of cells across generations. AVAILABILITY: CellAging, instructions and an example are available at http://www.cs.tut.fi/%7esanchesr/cellaging/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.