Intracellular peroxynitrite perturbs redox balance, bioenergetics, and Fe-S cluster homeostasis in Mycobacterium tuberculosis.

The ability of Mycobacterium tuberculosis (Mtb) to tolerate nitric oxide (•NO) and superoxide (O2•-) produced by phagocytes contributes to its success as a human pathogen. Recombination of •NO and O2•- generates peroxynitrite (ONOO-), a potent oxidant produced inside activated macrophages causing lethality in diverse organisms. While the response of Mtb toward •NO and O2•- is well established, how Mtb responds to ONOO- remains unclear. Filling this knowledge gap is important to understand the persistence mechanisms of Mtb during infection. We synthesized a series of compounds that generate both •NO and O2•-, which should combine to produce ONOO-. From this library, we identified CJ067 that permeates Mtb to reliably enhance intracellular ONOO- levels. CJ067-exposed Mtb strains, including multidrug-resistant (MDR) and extensively drug-resistant (XDR) clinical isolates, exhibited dose-dependent, long-lasting oxidative stress and growth inhibition. In contrast, Mycobacterium smegmatis (Msm), a fast-growing, non-pathogenic mycobacterial species, maintained redox balance and growth in response to intracellular ONOO-. RNA-sequencing with Mtb revealed that CJ067 induces antioxidant machinery, sulphur metabolism, metal homeostasis, and a 4Fe-4S cluster repair pathway (suf operon). CJ067 impaired the activity of the 4Fe-4S cluster-containing TCA cycle enzyme, aconitase, and diminished bioenergetics of Mtb. Work with Mtb strains defective in SUF and IscS involved in Fe-S cluster biogenesis pathways showed that both systems cooperatively protect Mtb from intracellular ONOO- in vitro and inducible nitric oxide synthase (iNOS)-dependent growth inhibition during macrophage infection. Thus, Mtb is uniquely sensitive to intracellular ONOO- and targeting Fe-S cluster homeostasis is expected to promote iNOS-dependent host immunity against tuberculosis (TB).

Cysteine desulfurase (IscS)-mediated fine-tuning of bioenergetics and SUF expression prevents Mycobacterium tuberculosis hypervirulence.

Iron-sulfur (Fe-S) biogenesis requires multiprotein assembly systems, SUF and ISC, in most prokaryotes. M. tuberculosis (Mtb) encodes a complete SUF system, the depletion of which was bactericidal. The ISC operon is truncated to a single gene iscS (cysteine desulfurase), whose function remains uncertain. Here, we show that MtbΔiscS is bioenergetically deficient and hypersensitive to oxidative stress, antibiotics, and hypoxia. MtbΔiscS resisted killing by nitric oxide (NO). RNA sequencing indicates that IscS is important for expressing regulons of DosR and Fe-S-containing transcription factors, WhiB3 and SufR. Unlike wild-type Mtb, MtbΔiscS could not enter a stable persistent state, continued replicating in mice, and showed hypervirulence. The suf operon was overexpressed in MtbΔiscS during infection in a NO-dependent manner. Suppressing suf expression in MtbΔiscS either by CRISPR interference or upon infection in inducible NO-deficient mice arrests hypervirulence. Together, Mtb redesigned the ISC system to "fine-tune" the expression of SUF machinery for establishing persistence without causing detrimental disease in the host.

Biosensor-integrated transposon mutagenesis reveals rv0158 as a coordinator of redox homeostasis in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (Mtb) is evolutionarily equipped to resist exogenous reactive oxygen species (ROS) but shows vulnerability to an increase in endogenous ROS (eROS). Since eROS is an unavoidable consequence of aerobic metabolism, understanding how Mtb manages eROS levels is essential yet needs to be characterized. By combining the Mrx1-roGFP2 redox biosensor with transposon mutagenesis, we identified 368 genes (redoxosome) responsible for maintaining homeostatic levels of eROS in Mtb. Integrating redoxosome with a global network of transcriptional regulators revealed a hypothetical protein (Rv0158) as a critical node managing eROS in Mtb. Disruption of rv0158 (rv0158 KO) impaired growth, redox balance, respiration, and metabolism of Mtb on glucose but not on fatty acids. Importantly, rv0158 KO exhibited enhanced growth on propionate, and the Rv0158 protein directly binds to methylmalonyl-CoA, a key intermediate in propionate catabolism. Metabolite profiling, ChIP-Seq, and gene-expression analyses indicate that Rv0158 manages metabolic neutralization of propionate toxicity by regulating the methylcitrate cycle. Disruption of rv0158 enhanced the sensitivity of Mtb to oxidative stress, nitric oxide, and anti-TB drugs. Lastly, rv0158 KO showed poor survival in macrophages and persistence defect in mice. Our results suggest that Rv0158 is a metabolic integrator for carbon metabolism and redox balance in Mtb.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models.

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, Emax, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with Emax but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and Emax values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Emergence of networks of shared restriction-modification systems in phage-bacteria ecosystems.

Restriction-modification (RM) systems are the most ubiquitous bacterial defence systems against bacteriophages. Using genome sequence data, we showed that RM systems are often shared among bacterial strains in a structured way. Examining the network of interconnections between bacterial strains within genera, we found that many strains share more RM systems than expected compared with a suitable null model. We also found that many genera have a larger than expected number of bacterial strains with unique RM systems. We used population dynamics models of closed and open phage-bacteria ecosystems to qualitatively understand the selection pressures that could lead to such network structures with enhanced overlap or uniqueness. In our models, we found that the phages impose a selection pressure that favours bacteria with greater number of RM systems, and higher overlap of RM systems with other strains, but in bacteria-dominated states, this is opposed by the increased cost-to-growth rate of these bacteria. Similar to what we observed in the genome data, we found that two distinct bacterial strategies emerge - strains either have a greater overlap than expected, or, at the other extreme, have unique RM systems. The former strategy appears to dominate when the repertoire of available RM systems is smaller but the average number of RM systems per strain is larger.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes.

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Global pleiotropic effects in adaptively evolved Escherichia coli lacking CRP reveal molecular mechanisms that define the growth physiology.

Evolution facilitates emergence of fitter phenotypes by efficient allocation of cellular resources in conjunction with beneficial mutations. However, system-wide pleiotropic effects that redress the perturbations to the apex node of the transcriptional regulatory networks remain unclear. Here, we elucidate that absence of global transcriptional regulator CRP in Escherichia coli results in alterations in key metabolic pathways under glucose respiratory conditions, favouring stress- or hedging-related functions over growth-enhancing functions. Further, we disentangle the growth-mediated effects from the CRP regulation-specific effects on these metabolic pathways. We quantitatively illustrate that the loss of CRP perturbs proteome efficiency, as evident from metabolic as well as ribosomal proteome fractions, that corroborated with intracellular metabolite profiles. To address how E. coli copes with such systemic defect, we evolved Δcrp mutant in the presence of glucose. Besides acquiring mutations in the promoter of glucose transporter ptsG, the evolved populations recovered the metabolic pathways to their pre-perturbed state coupled with metabolite re-adjustments, which altogether enabled increased growth. By contrast to Δcrp mutant, the evolved strains remodelled their proteome efficiency towards biomass synthesis, albeit at the expense of carbon efficiency. Overall, we comprehensively illustrate the genetic and metabolic basis of pleiotropic effects, fundamental for understanding the growth physiology.

Genome-scale reconstruction of Gcn4/ATF4 networks driving a growth program.

Growth and starvation are considered opposite ends of a spectrum. To sustain growth, cells use coordinated gene expression programs and manage biomolecule supply in order to match the demands of metabolism and translation. Global growth programs complement increased ribosomal biogenesis with sufficient carbon metabolism, amino acid and nucleotide biosynthesis. How these resources are collectively managed is a fundamental question. The role of the Gcn4/ATF4 transcription factor has been best studied in contexts where cells encounter amino acid starvation. However, high Gcn4 activity has been observed in contexts of rapid cell proliferation, and the roles of Gcn4 in such growth contexts are unclear. Here, using a methionine-induced growth program in yeast, we show that Gcn4/ATF4 is the fulcrum that maintains metabolic supply in order to sustain translation outputs. By integrating matched transcriptome and ChIP-Seq analysis, we decipher genome-wide direct and indirect roles for Gcn4 in this growth program. Genes that enable metabolic precursor biosynthesis indispensably require Gcn4; contrastingly ribosomal genes are partly repressed by Gcn4. Gcn4 directly binds promoter-regions and transcribes a subset of metabolic genes, particularly driving lysine and arginine biosynthesis. Gcn4 also globally represses lysine and arginine enriched transcripts, which include genes encoding the translation machinery. The Gcn4 dependent lysine and arginine supply thereby maintains the synthesis of the translation machinery. This is required to maintain translation capacity. Gcn4 consequently enables metabolic-precursor supply to bolster protein synthesis, and drive a growth program. Thus, we illustrate how growth and starvation outcomes are both controlled using the same Gcn4 transcriptional outputs that function in distinct contexts.

Whole-Genome Sequencing of Sphingobium sp. Strain RSMS, a Highly Efficient Tributyl Phosphate-Degrading Bacterium.

Sphingobium sp. strain RSMS was described earlier as an efficient degrader of tributyl phosphate, an organic pollutant. This report describes the generation and annotation of the genome sequence of Sphingobium sp. strain RSMS, which will facilitate future studies to identify genetic elements responsible for the degradation of tributyl phosphate.

Evolutionary and Comparative Analysis of Bacterial Nonhomologous End Joining Repair.

DNA double-strand breaks (DSBs) are a threat to genome stability. In all domains of life, DSBs are faithfully fixed via homologous recombination. Recombination requires the presence of an uncut copy of duplex DNA which is used as a template for repair. Alternatively, in the absence of a template, cells utilize error-prone nonhomologous end joining (NHEJ). Although ubiquitously found in eukaryotes, NHEJ is not universally present in bacteria. It is unclear as to why many prokaryotes lack this pathway. Toward understanding what could have led to the current distribution of bacterial NHEJ, we carried out comparative genomics and phylogenetic analysis across ∼6,000 genomes. Our results show that this pathway is sporadically distributed across the phylogeny. Ancestral reconstruction further suggests that NHEJ was absent in the eubacterial ancestor and can be acquired via specific routes. Integrating NHEJ occurrence data for archaea, we also find evidence for extensive horizontal exchange of NHEJ genes between the two kingdoms as well as across bacterial clades. The pattern of occurrence in bacteria is consistent with correlated evolution of NHEJ with key genome characteristics of genome size and growth rate; NHEJ presence is associated with large genome sizes and/or slow growth rates, with the former being the dominant correlate. Given the central role these traits play in determining the ability to carry out recombination, it is possible that the evolutionary history of bacterial NHEJ may have been shaped by requirement for efficient DSB repair.

A Mutant RNA Polymerase Activates the General Stress Response, Enabling Escherichia coli Adaptation to Late Prolonged Stationary Phase.

Escherichia coli populations undergo repeated replacement of parental genotypes with fitter variants deep in stationary phase. We isolated one such variant, which emerged after 3 weeks of maintaining an E. coli K-12 population in stationary phase. This variant displayed a small colony phenotype and slow growth and was able to outcompete its ancestor over a narrow time window in stationary phase. The variant also shows tolerance to beta-lactam antibiotics, though not previously exposed to the antibiotic. We show that an RpoC(A494V) mutation confers the slow growth and small colony phenotype on this variant. The ability of this mutation to confer a growth advantage in stationary phase depends on the availability of the stationary-phase sigma factor σS The RpoC(A494V) mutation upregulates the σS regulon. As shown over 20 years ago, early in prolonged stationary phase, σS attenuation, but not complete loss of activity, confers a fitness advantage. Our study shows that later mutations enhance σS activity, either by mutating the gene for σS directly or via mutations such as RpoC(A494V). The balance between the activities of the housekeeping major sigma factor and σS sets up a trade-off between growth and stress tolerance, which is tuned repeatedly during prolonged stationary phase.IMPORTANCE An important general mechanism of a bacterium's adaptation to its environment involves adjusting the balance between growing fast and tolerating stresses. One paradigm where this plays out is in prolonged stationary phase: early studies showed that attenuation, but not complete elimination, of the general stress response enables early adaptation of the bacterium E. coli to the conditions established about 10 days into stationary phase. We show here that this balance is not static and that it is tilted back in favor of the general stress response about 2 weeks later. This can be established by direct mutations in the master regulator of the general stress response or by mutations in the core RNA polymerase enzyme itself. These conditions can support the development of antibiotic tolerance although the bacterium is not exposed to the antibiotic. Further exploration of the growth-stress balance over the course of stationary phase will necessarily require a deeper understanding of the events in the extracellular milieu.

Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria.

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

Laboratory Evolution Experiments Help Identify a Predominant Region of Constitutive Stable DNA Replication Initiation.

The bacterium Escherichia coli can initiate replication in the absence of the replication initiator protein DnaA and/or the canonical origin of replication oriC in a ΔrnhA background. This phenomenon, which can be primed by R-loops, is called constitutive stable DNA replication (cSDR). Whether DNA replication during cSDR initiates in a stochastic manner through the length of the chromosome or at specific sites and how E. coli can find adaptations to loss of fitness caused by cSDR remain inadequately answered. We use laboratory evolution experiments of ΔrnhA-ΔdnaA strains followed by deep sequencing to show that DNA replication preferentially initiates within a broad region located ∼0.4 to 0.7 Mb clockwise of oriC. This region includes many bisulfite-sensitive sites, which have been previously defined as R-loop-forming regions, and includes a site containing sequence motifs that favor R-loop formation. Initiation from this region would result in head-on replication-transcription conflicts at rRNA loci. Inversions of these rRNA loci, which can partly resolve these conflicts, help the bacterium suppress the fitness defects of cSDR. These inversions partially restore the gene expression changes brought about by cSDR. The inversion, however, increases the possibility of conflicts at essential mRNA genes, which would utilize only a minuscule fraction of RNA polymerase molecules, most of which transcribe rRNA genes. Whether subsequent adaptive strategies would attempt to resolve these conflicts remains an open question.IMPORTANCE The bacterium E. coli can replicate its DNA even in the absence of the molecules that are required for canonical replication initiation. This often requires the formation of RNA-DNA hybrid structures and is referred to as constitutive stable DNA replication (cSDR). Where on the chromosome does cSDR initiate? We answer this question using laboratory evolution experiments and genomics and show that selection favors cSDR initiation predominantly at a region ∼0.6 Mb clockwise of oriC. Initiation from this site will result in more head-on collisions of DNA polymerase with RNA polymerase operating on rRNA loci. The bacterium adapts to this problem by inverting a region of the genome including several rRNA loci such that head-on collisions between the two polymerases are minimized. Understanding such evolutionary strategies in the context of cSDR can provide insights into the potential causes of resistance to antibiotics that target initiation of DNA replication.

Early fate of exogenous promoters in E. coli.

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Targeting redox heterogeneity to counteract drug tolerance in replicating Mycobacterium tuberculosis.

The capacity of Mycobacterium tuberculosis (Mtb) to tolerate multiple antibiotics represents a major problem in tuberculosis (TB) management. Heterogeneity in Mtb populations is one of the factors that drives antibiotic tolerance during infection. However, the mechanisms underpinning this variation in bacterial population remain poorly understood. Here, we show that phagosomal acidification alters the redox physiology of Mtb to generate a population of replicating bacteria that display drug tolerance during infection. RNA sequencing of this redox-altered population revealed the involvement of iron-sulfur (Fe-S) cluster biogenesis, hydrogen sulfide (H2S) gas, and drug efflux pumps in antibiotic tolerance. The fraction of the pH- and redox-dependent tolerant population increased when Mtb infected macrophages with actively replicating HIV-1, suggesting that redox heterogeneity could contribute to high rates of TB therapy failure during HIV-TB coinfection. Pharmacological inhibition of phagosomal acidification by the antimalarial drug chloroquine (CQ) eradicated drug-tolerant Mtb, ameliorated lung pathology, and reduced postchemotherapeutic relapse in in vivo models. The pharmacological profile of CQ (C max and AUClast) exhibited no major drug-drug interaction when coadministered with first line anti-TB drugs in mice. Our data establish a link between phagosomal pH, redox metabolism, and drug tolerance in replicating Mtb and suggest repositioning of CQ to shorten TB therapy and achieve a relapse-free cure.

Redox-dependent condensation of the mycobacterial nucleoid by WhiB4.

Oxidative stress response in bacteria is mediated through coordination between the regulators of oxidant-remediation systems (e.g. OxyR, SoxR) and nucleoid condensation (e.g. Dps, Fis). However, these genetic factors are either absent or rendered non-functional in the human pathogen Mycobacterium tuberculosis (Mtb). Therefore, how Mtb organizes genome architecture and regulates gene expression to counterbalance oxidative imbalance is unknown. Here, we report that an intracellular redox-sensor, WhiB4, dynamically links genome condensation and oxidative stress response in Mtb. Disruption of WhiB4 affects the expression of genes involved in maintaining redox homeostasis, central metabolism, and respiration under oxidative stress. Notably, disulfide-linked oligomerization of WhiB4 in response to oxidative stress activates the protein's ability to condense DNA. Further, overexpression of WhiB4 led to hypercondensation of nucleoids, redox imbalance and increased susceptibility to oxidative stress, whereas WhiB4 disruption reversed this effect. In accordance with the findings in vitro, ChIP-Seq data demonstrated non-specific binding of WhiB4 to GC-rich regions of the Mtb genome. Lastly, data indicate that WhiB4 deletion affected the expression of ~ 30% of genes preferentially bound by the protein, suggesting both direct and indirect effects on gene expression. We propose that WhiB4 structurally couples Mtb's response to oxidative stress with genome organization and transcription.

Regulation of Global Transcription in Escherichia coli by Rsd and 6S RNA.

In Escherichia coli, the sigma factor σ70 directs RNA polymerase to transcribe growth-related genes, while σ38 directs transcription of stress response genes during stationary phase. Two molecules hypothesized to regulate RNA polymerase are the protein Rsd, which binds to σ70, and the non-coding 6S RNA which binds to the RNA polymerase-σ70 holoenzyme. Despite multiple studies, the functions of Rsd and 6S RNA remain controversial. Here we use RNA-Seq in five phases of growth to elucidate their function on a genome-wide scale. We show that Rsd and 6S RNA facilitate σ38 activity throughout bacterial growth, while 6S RNA also regulates widely different genes depending upon growth phase. We discover novel interactions between 6S RNA and Rsd and show widespread expression changes in a strain lacking both regulators. Finally, we present a mathematical model of transcription which highlights the crosstalk between Rsd and 6S RNA as a crucial factor in controlling sigma factor competition and global gene expression.

Using stochastic cell division and death to probe minimal units of cellular replication.

The invariant cell initiation mass measured in bacterial growth experiments has been interpreted as a minimal unit of cellular replication. Here we argue that the existence of such minimal units induces a coupling between the rates of stochastic cell division and death. To probe this coupling we tracked live and dead cells in Escherichia coli populations treated with a ribosome-targeting antibiotic. We find that the growth exponent from macroscopic cell growth or decay measurements can be represented as the difference of microscopic first-order cell division and death rates. The boundary between cell growth and decay, at which the number of live cells remains constant over time, occurs at the minimal inhibitory concentration (MIC) of the antibiotic. This state appears macroscopically static but is microscopically dynamic: division and death rates exactly cancel at MIC but each is remarkably high, reaching 60% of the antibiotic-free division rate. A stochastic model of cells as collections of minimal replicating units we term 'widgets' reproduces both steady-state and transient features of our experiments. Sub-cellular fluctuations of widget numbers stochastically drive each new daughter cell to one of two alternate fates, division or death. First-order division or death rates emerge as eigenvalues of a stationary Markov process, and can be expressed in terms of the widget's molecular properties. High division and death rates at MIC arise due to low mean and high relative fluctuations of widget number. Isolating cells at the threshold of irreversible death might allow molecular characterization of this minimal replication unit.