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.

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.

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.

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.

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.

Robustness of the division symmetry in Escherichia coli and functional consequences of symmetry breaking.

The morphological symmetry of the division process of Escherichia coli is well-known. Recent studies verified that, in optimal growth conditions, most divisions are symmetric, although there are exceptions. We investigate whether such morphological asymmetries in division introduce functional asymmetries between sister cells, and assess the robustness of the symmetry in division to mild chemical stresses and sub-optimal temperatures. First, we show that the difference in size between daughter cells at birth is positively correlated to the difference between the numbers of fluorescent protein complexes inherited from the parent cell. Next, we show that the degree of symmetry in division observed in optimal conditions is robust to mild acidic shift and to mild oxidative stress, but not to sub-optimal temperatures, in that the variance of the difference between the sizes of sister cells at birth is minimized at 37 °C. This increased variance affects the functionality of the cells in that, at sub-optimal temperatures, larger/smaller cells arising from asymmetric divisions exhibit faster/slower division times than the mean population division time, respectively. On the other hand, cells dividing faster do not do so at the cost of morphological symmetry in division. Finally we show that at suboptimal temperatures the mean distance between the nucleoids increases, explaining the increased variance in division. We conclude that the functionality of E. coli cells is not immune to morphological asymmetries at birth, and that the effectiveness of the mechanism responsible for ensuring the symmetry in division weakens at sub-optimal temperatures.

A library of reporters of the global regulators of gene expression in Escherichia coli.

The topology of the transcription factor network (TFN) of Escherichia coli is far from uniform, with 22 global regulator (GR) proteins controlling one-third of all genes. So far, their production rates cannot be tracked by comparable fluorescent proteins. We developed a library of fluorescent reporters for 16 GRs for this purpose. Each consists of a single-copy plasmid coding for green fluorescent protein (GFP) fused to the full-length copy of the native promoter. We tracked their activity in exponential and stationary growth, as well as under weak and strong stresses. We show that the reporters have high sensitivity and specificity to all stresses tested and detect single-cell variability in transcription rates. Given the influence of GRs on the TFN, we expect that the new library will contribute to dissecting global transcriptional stress-response programs of E. coli. Moreover, the library can be invaluable in bioindustrial applications that tune those programs to, instead of cell growth, favor productivity while reducing energy consumption.IMPORTANCECells contain thousands of genes. Many genes are involved in the control of cellular activities. Some activities require a few hundred genes to run largely synchronous transcriptional programs. To achieve this, cells have evolved global regulator (GR) proteins that can influence hundreds of genes simultaneously. We have engineered a library of Escherichia coli strains to track the levels over time of these, phenotypically critical, GRs. Each strain has a single-copy plasmid coding for a fast-maturing green fluorescent protein whose transcription is controlled by a copy of the natural GR promoter. By allowing the tracking of GR levels, with sensitivity and specificity, this library should become of wide use in scientific research on bacterial gene expression (from molecular to synthetic biology) and, later, be used in applications in therapeutics and bioindustries.

Hardware, Software, and Wetware Codesign Environment for Synthetic Biology.

Synthetic biology is the process of forward engineering living systems. These systems can be used to produce biobased materials, agriculture, medicine, and energy. One approach to designing these systems is to employ techniques from the design of embedded electronics. These techniques include abstraction, standards, modularity, automated design, and formal semantic models of computation. Together, these elements form the foundation of "biodesign automation," where software, robotics, and microfluidic devices combine to create exciting biological systems of the future. This paper describes a "hardware, software, wetware" codesign vision where software tools can be made to act as "genetic compilers" that transform high-level specifications into engineered "genetic circuits" (wetware). This is followed by a process where automation equipment, well-defined experimental workflows, and microfluidic devices are explicitly designed to house, execute, and test these circuits (hardware). These systems can be used as either massively parallel experimental platforms or distributed bioremediation and biosensing devices. Next, scheduling and control algorithms (software) manage these systems' actual execution and data analysis tasks. A distinguishing feature of this approach is how all three of these aspects (hardware, software, and wetware) may be derived from the same basic specification in parallel and generated to fulfill specific cost, performance, and structural requirements.

Genetic circuit design automation with Cello 2.0.

Cells interact with their environment, communicate among themselves, track time and make decisions through functions controlled by natural regulatory genetic circuits consisting of interacting biological components. Synthetic programmable circuits used in therapeutics and other applications can be automatically designed by computer-aided tools. The Cello software designs the DNA sequences for programmable circuits based on a high-level software description and a library of characterized DNA parts representing Boolean logic gates. This process allows for design specification reuse, modular DNA part library curation and formalized circuit transformations based on experimental data. This protocol describes Cello 2.0, a freely available cross-platform software written in Java. Cello 2.0 enables flexible descriptions of the logic gates' structure and their mathematical models representing dynamic behavior, new formal rules for describing the placement of gates in a genome, a new graphical user interface, support for Verilog 2005 syntax and a connection to the SynBioHub parts repository software environment. Collectively, these features expand Cello's capabilities beyond Escherichia coli plasmids to new organisms and broader genetic contexts, including the genome. Designing circuits with Cello 2.0 produces an abstract Boolean network from a Verilog file, assigns biological parts to each node in the Boolean network, constructs a DNA sequence and generates highly structured and annotated sequence representations suitable for downstream processing and fabrication, respectively. The result is a sequence implementing the specified Boolean function in the organism and predictions of circuit performance. Depending on the size of the design space and users' expertise, jobs may take minutes or hours to complete.

Buildout and integration of an automated high-throughput CLIA laboratory for SARS-CoV-2 testing on a large urban campus.

In 2019, the first cases of SARS-CoV-2 were detected in Wuhan, China, and by early 2020 the first cases were identified in the United States. SARS-CoV-2 infections increased in the US causing many states to implement stay-at-home orders and additional safety precautions to mitigate potential outbreaks. As policies changed throughout the pandemic and restrictions lifted, there was an increase in demand for COVID-19 testing which was costly, difficult to obtain, or had long turn-around times. Some academic institutions, including Boston University (BU), created an on-campus COVID-19 screening protocol as part of a plan for the safe return of students, faculty, and staff to campus with the option for in-person classes. At BU, we put together an automated high-throughput clinical testing laboratory with the capacity to run 45,000 individual tests weekly by Fall of 2020, with a purpose-built clinical testing laboratory, a multiplexed reverse transcription PCR (RT-qPCR) test, robotic instrumentation, and trained staff. There were many challenges including supply chain issues for personal protective equipment and testing materials in addition to equipment that were in high demand. The BU Clinical Testing Laboratory (CTL) was operational at the start of Fall 2020 and performed over 1 million SARS-CoV-2 PCR tests during the 2020-2021 academic year.

Transcription closed and open complex formation coordinate expression of genes with a shared promoter region.

Many genes are spaced closely, allowing coordination without explicit control through shared regulatory elements and molecular interactions. We study the dynamics of a stochastic model of a gene-pair in a head-to-head configuration, sharing promoter elements, which accounts for the rate-limiting steps in transcription initiation. We find that only in specific regions of the parameter space of the rate-limiting steps is orderly coexpression exhibited, suggesting that successful cooperation between closely spaced genes requires the coevolution of compatible rate-limiting step configuration. The model predictions are validated using in vivo single-cell, single-RNA measurements of the dynamics of pairs of genes sharing promoter elements. Our results suggest that, in E. coli, the kinetics of the rate-limiting steps in active transcription can play a central role in shaping the dynamics of gene-pairs sharing promoter elements.

Estimating RNA numbers in single cells by RNA fluorescent tagging and flow cytometry.

Estimating the statistics of single-cell RNA numbers has become a key source of information on gene expression dynamics. One of the most informative methods of in vivo single-RNA detection is MS2d-GFP tagging. So far, it requires microscopy and laborious semi-manual image analysis, which hampers the amount of collectable data. To overcome this limitation, we present a new methodology for quantifying the mean, standard deviation, and skewness of single-cell distributions of RNA numbers, from flow cytometry data on cells expressing RNA tagged with MS2d-GFP. The quantification method, based on scaling flow-cytometry data from microscopy single-cell data on integer-valued RNA numbers, is shown to readily produce precise, big data on in vivo single-cell distributions of RNA numbers and, thus, can assist in studies of transcription dynamics.

Chromosome and plasmid-borne PLacO3O1 promoters differ in sensitivity to critically low temperatures.

Temperature shifts trigger genome-wide changes in Escherichia coli's gene expression. We studied if chromosome integration impacts on a gene's sensitivity to these shifts, by comparing the single-RNA production kinetics of a PLacO3O1 promoter, when chromosomally-integrated and when single-copy plasmid-borne. At suboptimal temperatures their induction range, fold change, and response to decreasing temperatures are similar. At critically low temperatures, the chromosome-integrated promoter becomes weaker and noisier. Dissection of its initiation kinetics reveals longer lasting states preceding open complex formation, suggesting enhanced supercoiling buildup. Measurements with Gyrase and Topoisomerase I inhibitors suggest hindrance to escape supercoiling buildup at low temperatures. Consistently, similar phenomena occur in energy-depleted cells by DNP at 30 °C. Transient, critically-low temperatures have no long-term consequences, as raising temperature quickly restores transcription rates. We conclude that the chromosomally-integrated PLacO3O1 has higher sensitivity to low temperatures, due to longer-lasting super-coiled states. A lesser active, chromosome-integrated native lac is shown to be insensitive to Gyrase overexpression, even at critically low temperatures, indicating that the rate of escaping positive supercoiling buildup is temperature and transcription rate dependent. A genome-wide analysis supports this, since cold-shock genes exhibit atypical supercoiling-sensitivities. This phenomenon might partially explain the temperature-sensitivity of some transcriptional programs of E. coli.

Dissecting the stochastic transcription initiation process in live Escherichia coli.

We investigate the hypothesis that, in Escherichia coli, while the concentration of RNA polymerases differs in different growth conditions, the fraction of RNA polymerases free for transcription remains approximately constant within a certain range of these conditions. After establishing this, we apply a standard model-fitting procedure to fully characterize the in vivo kinetics of the rate-limiting steps in transcription initiation of the Plac/ara-1 promoter from distributions of intervals between transcription events in cells with different RNA polymerase concentrations. We find that, under full induction, the closed complex lasts ∼788 s while subsequent steps last ∼193 s, on average. We then establish that the closed complex formation usually occurs multiple times prior to each successful initiation event. Furthermore, the promoter intermittently switches to an inactive state that, on average, lasts ∼87 s. This is shown to arise from the intermittent repression of the promoter by LacI. The methods employed here should be of use to resolve the rate-limiting steps governing the in vivo dynamics of initiation of prokaryotic promoters, similar to established steady-state assays to resolve the in vitro dynamics.

Kinetics of the cellular intake of a gene expression inducer at high concentrations.

From in vivo single-event measurements of the transient and steady-state transcription activity of a single-copy lac-ara-1 promoter in Escherichia coli, we characterize the intake kinetics of its inducer (IPTG) from the media. We show that the empirical data are well-fit by a model of intake assuming a bilayer membrane, with the passage through the second layer being rate-limiting, coupled to a stochastic, sub-Poissonian, multi-step transcription process. Using this model, we show that for a wide range of extracellular inducer levels (up to 1.25 mM) the intake process is diffusive-like, suggesting unsaturated membrane permeability. Inducer molecules travel from the periplasm to the cytoplasm in, on average, 31.7 minutes, strongly affecting cells' response time. The novel methodology followed here should aid the study of cellular intake mechanisms at the single-event level.

Single-cell kinetics of a repressilator when implemented in a single-copy plasmid.

Synthetic genetic clocks, such as the Elowitz-Leibler repressilator, will be key regulatory components of future synthetic circuits. We constructed a single-copy repressilator (SCR) by implementing the original repressilator circuit on a single-copy F-plasmid. After verifying its functionality, we studied its behaviour as a function of temperature and compared it with that of the original low-copy-number repressilator (LCR). Namely, we compared the period of oscillations, functionality (the fraction of cells exhibiting oscillations) and robustness to internal fluctuations (the fraction of expected oscillations that would occur). We found that, under optimal temperature conditions, the dynamics of the two systems differs significantly, although qualitatively they respond similarly to temperature changes. Exception to this is in the functionality, in which the SCR is higher at lower temperatures but lower at higher temperatures. Next, by adding IPTG to the medium at low and high concentrations during microscopy sessions, we showed that the functionality of the SCR is more robust to external perturbations, which indicates that the oscillatory behaviour of the LCR can be disrupted by affecting only a few of the copies in a cell. We conclude that the SCR, the first functional, synthetic, single-copy, ring-type genetic clock, is more robust to lower temperatures and to external perturbations than the original LCR. The SCR will be of use in future synthetic circuits, since it complements the array of tasks that the LCR can perform.

Effects of temperature on the dynamics of the LacI-TetR-CI repressilator.

We studied the behaviour of the repressilator at 28 °C, 30 °C, 33 °C, and 37 °C. From the fluorescence in each cell over time, we determined the period of oscillations, the functionality (fraction of cells exhibiting oscillations) and the robustness (fraction of expected oscillations that occur) of this circuit. We show that the oscillatory dynamics differs with temperature. Functionality is maximized at 30 °C. Robustness decreases beyond 30 °C, as most cells exhibit 'failed' oscillations. These failures cause the distribution of periods to become bimodal, with an 'apparent period' that is minimal at 30 °C, while the true period decreases with increasing temperature. Based on previous studies, we hypothesized that the failures are due to a loss of functionality of one protein of the repressilator, CI. To test this, we studied the kinetics of a genetic switch, formed by the proteins CI and Cro, whose expression is controlled by PRM and PR, respectively. By probing the activity of PRM by in vivo detection of MS2-GFP tagged RNA, we find that, beyond 30 °C, the production of the CI-coding RNA changes from sub-Poissonian to super-Poissonian. Given this, we suggest that the decrease in efficiency of CI as a repressor with temperature hinders the robustness of the repressilator beyond 30 °C. We conclude that the repressilator is sensitive but not robust to temperature. Replacing CI for a less temperature-dependent protein should enhance robustness.