Protection of bacteriophage-sensitive Escherichia coli by lysogens.

SignificanceSome viruses that infect bacteria, temperate bacteriophages, can confer immunity to infection by the same virus. Here we report λ-immune bacteria could protect λ-sensitive bacteria from killing by phage λ in mixed culture. The protection depended on the extent to which the immune bacteria were able to adsorb the phage. Reconciling modeling with experiment led to identifying a decline in protection as bacteria stopped growing. Adsorption of λ was compromised by inhibition of bacterial energy metabolism, explaining the loss of protection as bacterial growth ceased.

The dynamics of phage predation on a microcolony.

Phage predation is an important factor for controlling the bacterial biomass. At face value, dense microbial habitats are expected to be vulnerable to phage epidemics due to the abundance of fresh hosts immediately next to any infected bacteria. Despite this, the bacterial microcolony is a common habitat for bacteria in nature. Here, we experimentally quantify the fate of microcolonies of Escherichia coli exposed to virulent phage T4. It has been proposed that the outer bacterial layers of the colony will shield the inner layers from the phage invasion and thereby constrain the phage to the colony's surface. We develop a dynamical model that incorporates this shielding mechanism and fit the results with experimental measurements to extract important phage-bacteria interaction parameters. The analysis suggests that, while the shielding mechanism delays phage attack, T4 phage are able to diffuse so deep into the dense bacterial environment that colony-level survival of the bacterial community is challenged.

Host heterogeneity and epistasis explain punctuated evolution of SARS-CoV-2.

Identifying drivers of viral diversity is key to understanding the evolutionary as well as epidemiological dynamics of the COVID-19 pandemic. Using rich viral genomic data sets, we show that periods of steadily rising diversity have been punctuated by sudden, enormous increases followed by similarly abrupt collapses of diversity. We introduce a mechanistic model of saltational evolution with epistasis and demonstrate that these features parsimoniously account for the observed temporal dynamics of inter-genomic diversity. Our results provide support for recent proposals that saltational evolution may be a signature feature of SARS-CoV-2, allowing the pathogen to more readily evolve highly transmissible variants. These findings lend theoretical support to a heightened awareness of biological contexts where increased diversification may occur. They also underline the power of pathogen genomics and other surveillance streams in clarifying the phylodynamics of emerging and endemic infections. In public health terms, our results further underline the importance of equitable distribution of up-to-date vaccines.

Overdispersion in COVID-19 increases the effectiveness of limiting nonrepetitive contacts for transmission control.

Increasing evidence indicates that superspreading plays a dominant role in COVID-19 transmission. Recent estimates suggest that the dispersion parameter k for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is on the order of 0.1, which corresponds to about 10% of cases being the source of 80% of infections. To investigate how overdispersion might affect the outcome of various mitigation strategies, we developed an agent-based model with a social network that allows transmission through contact in three sectors: "close" (a small, unchanging group of mutual contacts as might be found in a household), "regular" (a larger, unchanging group as might be found in a workplace or school), and "random" (drawn from the entire model population and not repeated regularly). We assigned individual infectivity from a gamma distribution with dispersion parameter k We found that when k was low (i.e., greater heterogeneity, more superspreading events), reducing random sector contacts had a far greater impact on the epidemic trajectory than did reducing regular contacts; when k was high (i.e., less heterogeneity, no superspreading events), that difference disappeared. These results suggest that overdispersion of COVID-19 transmission gives the virus an Achilles' heel: Reducing contacts between people who do not regularly meet would substantially reduce the pandemic, while reducing repeated contacts in defined social groups would be less effective.

Establishment of heterochromatin in domain-size-dependent bursts.

Methylation of histone H3K9 is a hallmark of epigenetic silencing in eukaryotes. Nucleosome modifications often rely on positive feedback where enzymes are recruited by modified nucleosomes. A combination of local and global feedbacks has been proposed to account for some dynamic properties of heterochromatin, but the range at which the global feedbacks operate and the exact mode of heterochromatin propagation are not known. We investigated these questions in fission yeast. Guided by mathematical modeling, we incrementally increased the size of the mating-type region and profiled heterochromatin establishment over time. We observed exponential decays in the proportion of cells with active reporters, with rates that decreased with domain size. Establishment periods varied from a few generations in wild type to >200 generations in the longest region examined, and highly correlated silencing of two reporters located outside the nucleation center was observed. On a chromatin level, this indicates that individual regions are silenced in sudden bursts. Mathematical modeling accounts for these bursts if heterochromatic nucleosomes facilitate a deacetylation or methylation reaction at long range, in a distance-independent manner. A likely effector of three-dimensional interactions is the evolutionarily conserved Swi6HP1 H3K9me reader, indicating the bursting behavior might be a general mode of heterochromatin propagation.

Epigenetic switching with asymmetric bridging interactions.

Gene expression states are often stably sustained in cis despite massively disruptive events like DNA replication. This is achieved by on-going enzymatic activity that maintains parts of the DNA in either heterochromatic (packed) or euchromatic (free) states, each of which is stabilized by both positive feedback and bridging interactions between individual nucleosomes. In contrast to condensed matter, however, the dynamics is not only governed by equilibrium binding interactions but is also mediated by enzymes that recognize and act on specific amino acid tails of the nucleosomes. The mechanical result is that some nucleosomes can bind to one another and form tightly packed polymer configurations, whereas others remain unbound and form free, noncompact polymer configurations. Here, we study the consequences of such an asymmetric interaction pattern on the dynamics of epigenetic switching. We develop a 3D polymer model and show that traits associated with epigenetic switching, such as bistability and epigenetic memory, are permitted by such a model. We find, however, that the experimentally observed burst-like nature of some epigenetic switches is difficult to reproduce by this biologically motivated interaction. Instead, the behavior seen in experiments can be explained by introducing partial confinement, which particularly affects the euchromatic regions of the chromosome.

From kill the winner to eliminate the winner in open phage-bacteria systems.

Phages and bacteria manage to coexist and sustain ecosystems with a high diversity of strains, despite limited resources and heavy predation. This diversity can be explained by the "kill the winner" model where virulent phages predominantly prey on fast-growing bacteria and thereby suppress the competitive exclusion of slower-growing bacteria. Here we computationally investigate the robustness of these systems against invasions, where new phages or bacteria may interact with more than one of the resident strains. The resulting interaction networks were found to self-organize into a network with strongly interacting specialized predator-prey pairs, resembling that of the "kill the winner" model. Furthermore, the "kill the winner" dynamics is enforced with the occasional elimination of even the fastest-growing bacteria strains due to a phage infecting the fast and slow growers. The frequency of slower-growing strains was increased with the introduction of even a few non-diagonal interactions. Hence, phages capable of infecting multiple hosts play significant roles both in the evolution of the ecosystem by eliminating the winner and in supporting diversity by allowing slow growers to coexist with faster growers.

COVID-19 Superspreading Suggests Mitigation by Social Network Modulation.

Although COVID-19 has caused severe suffering globally, the efficacy of nonpharmaceutical interventions has been greater than typical models have predicted. Meanwhile, evidence is mounting that the pandemic is characterized by superspreading. Capturing this phenomenon theoretically requires modeling at the scale of individuals. Using a mathematical model, we show that superspreading drastically enhances mitigations which reduce the overall personal contact number and that social clustering increases this effect.

Deconstructing the principles of ductal network formation in the pancreas.

The mammalian pancreas is a branched organ that does not exhibit stereotypic branching patterns, similarly to most other glands. Inside branches, it contains a network of ducts that undergo a transition from unconnected microlumen to a mesh of interconnected ducts and finally to a treelike structure. This ductal remodeling is poorly understood, both on a microscopic and macroscopic level. In this article, we quantify the network properties at different developmental stages. We find that the pancreatic network exhibits stereotypic traits at each stage and that the network properties change with time toward the most economical and optimized delivery of exocrine products into the duodenum. Using in silico modeling, we show how steps of pancreatic network development can be deconstructed into two simple rules likely to be conserved for many other glands. The early stage of the network is explained by noisy, redundant duct connection as new microlumens form. The later transition is attributed to pruning of the network based on the flux of fluid running through the pancreatic network into the duodenum.

A growing microcolony can survive and support persistent propagation of virulent phages.

Bacteria form colonies and secrete extracellular polymeric substances that surround the individual cells. These spatial structures are often associated with collaboration and quorum sensing between the bacteria. Here we investigate the mutual protection provided by spherical growth of a monoclonal colony during exposure to phages that proliferate on its surface. As a proof of concept we exposed growing colonies of Escherichia coli to a virulent mutant of phage P1. When the colony consists of less than [Formula: see text]50,000 members it is eliminated, while larger initial colonies allow long-term survival of both phage-resistant mutants and, importantly, colonies of mostly phage-sensitive members. A mathematical model predicts that colonies formed solely by phage-sensitive bacteria can survive because the growth of bacteria throughout the colony exceeds the killing of bacteria on the surface and pinpoints how the critical colony size depends on key parameters in the phage infection cycle.

On Phage Adsorption to Bacterial Chains.

Bacteria often arrange themselves in various spatial configurations, which changes how they interact with their surroundings. In this work, we investigate how the structure of the bacterial arrangements influences the adsorption of bacteriophages. We quantify how the adsorption rate scales with the number of bacteria in the arrangement and show that the adsorption rates for microcolonies (increasing with exponent ∼1/3) and bacterial chains (increasing with exponent ∼0.5-0.8) are substantially lower than for well-mixed bacteria (increasing with exponent 1). We further show that, after infection, the spatially clustered arrangements reduce the effective burst size by more than 50% and cause substantial superinfections in a very short time interval after phage lysis.

Optimal Transport Flows for Distributed Production Networks.

Network flows often exhibit a hierarchical treelike structure that can be attributed to the minimization of dissipation. The common feature of such systems is a single source and multiple sinks (or vice versa). In contrast, here we study networks with only a single source and sink. These systems can arise from secondary purposes of the networks, such as blood sugar regulation through insulin production. Minimization of dissipation in these systems leads to vascular shunting, a single vessel connecting the inlet and outlet. We show instead how optimizing the transport time yields network topologies that match those observed in the insulin-producing pancreatic islets. These are patterns of periphery-to-center and center-to-periphery flows. The obtained flow networks are broadly independent of how the flow velocity depends on the flow flux, but continuous and discontinuous phase transitions appear at extreme flux dependencies. Lastly, we show how constraints on flows can lead to buckling of the branches of the network, a feature that is also observed in pancreatic islets.

Four simple rules that are sufficient to generate the mammalian blastocyst.

Early mammalian development is both highly regulative and self-organizing. It involves the interplay of cell position, predetermined gene regulatory networks, and environmental interactions to generate the physical arrangement of the blastocyst with precise timing. However, this process occurs in the absence of maternal information and in the presence of transcriptional stochasticity. How does the preimplantation embryo ensure robust, reproducible development in this context? It utilizes a versatile toolbox that includes complex intracellular networks coupled to cell-cell communication, segregation by differential adhesion, and apoptosis. Here, we ask whether a minimal set of developmental rules based on this toolbox is sufficient for successful blastocyst development, and to what extent these rules can explain mutant and experimental phenotypes. We implemented experimentally reported mechanisms for polarity, cell-cell signaling, adhesion, and apoptosis as a set of developmental rules in an agent-based in silico model of physically interacting cells. We find that this model quantitatively reproduces specific mutant phenotypes and provides an explanation for the emergence of heterogeneity without requiring any initial transcriptional variation. It also suggests that a fixed time point for the cells' competence of fibroblast growth factor (FGF)/extracellular signal-regulated kinase (ERK) sets an embryonic clock that enables certain scaling phenomena, a concept that we evaluate quantitatively by manipulating embryos in vitro. Based on these observations, we conclude that the minimal set of rules enables the embryo to experiment with stochastic gene expression and could provide the robustness necessary for the evolutionary diversification of the preimplantation gene regulatory network.

The uneasy coexistence of predators and pathogens.

Disease and predation are both highly important in ecology, and pathogens with multiple host species have turned out to be common. Nonetheless, the interplay between multi-host epidemics and predation has received relatively little attention. Here, we analyse a model of a predator-prey system with disease in both prey and predator populations and determine reasonable parameter values using allometric mass scaling relations. Our analysis focuses on the possibility of extinction events rather than the linear stability of the model equations, and we derive approximate relations for the parameter values at which we expect these events to occur. We find that if the predator is a specialist, epidemics frequently drive the predator species to extinction. If the predator has an additional, immune prey species, predators will usually survive. Coexistence of predator and disease is impossible in the single-prey model. We conclude that for the prey species, carrying a pathogen can be an effective weapon against predators, and that being a generalist is a major advantage for a predator in the event of an epidemic affecting the prey or both species.

Efficient chromosomal-scale DNA looping in Escherichia coli using multiple DNA-looping elements.

Genes are frequently regulated by interactions between proteins that bind to the DNA near the gene and proteins that bind to DNA sites located far away, with the intervening DNA looped out. But it is not understood how efficient looping can occur when the sites are very far apart. We develop a simple theoretical framework that relates looping efficiency to the energetic cost and benefit of looping, allowing prediction of the efficiency of single or multiple nested loops at different distances. Measurements of absolute loop efficiencies for Lac repressor and λ CI using gene expression reporters in Escherichia coli cells show that, as predicted by the model, long-range DNA looping between a pair of sites can be strongly enhanced by the use of nested DNA loops or by the use of additional protein-binding sequences. A combination of these approaches was able to generate efficient DNA looping at a 200 kb distance.

DNA methylation in human epigenomes depends on local topology of CpG sites.

In vertebrates, methylation of cytosine at CpG sequences is implicated in stable and heritable patterns of gene expression. The classical model for inheritance, in which individual CpG sites are independent, provides no explanation for the observed non-random patterns of methylation. We first investigate the exact topology of CpG clustering in the human genome associated to CpG islands. Then, by pooling genomic CpG clusters on the basis of short distances between CpGs within and long distances outside clusters, we show a strong dependence of methylation on the number and density of CpG organization. CpG clusters with fewer, or less densely spaced, CpGs are predominantly hyper-methylated, while larger clusters are predominantly hypo-methylated. Intermediate clusters, however, are either hyper- or hypo-methylated but are rarely found in intermediate methylation states. We develop a model for spatially-dependent collaboration between CpGs, where methylated CpGs recruit methylation enzymes that can act on CpGs over an extended local region, while unmethylated CpGs recruit demethylation enzymes that act more strongly on nearby CpGs. This model can reproduce the effects of CpG clustering on methylation and produces stable and heritable alternative methylation states of CpG clusters, thus providing a coherent model for methylation inheritance and methylation patterning.

Models of life: epigenetics, diversity and cycles.

This review emphasizes aspects of biology that can be understood through repeated applications of simple causal rules. The selected topics include perspectives on gene regulation, phage lambda development, epigenetics, microbial ecology, as well as model approaches to diversity and to punctuated equilibrium in evolution. Two outstanding features are repeatedly described. One is the minimal number of rules to sustain specific states of complex systems for a long time. The other is the collapse of such states and the subsequent dynamical cycle of situations that restitute the system to a potentially new metastable state.

Food Web Assembly Rules for Generalized Lotka-Volterra Equations.

In food webs, many interacting species coexist despite the restrictions imposed by the competitive exclusion principle and apparent competition. For the generalized Lotka-Volterra equations, sustainable coexistence necessitates nonzero determinant of the interaction matrix. Here we show that this requirement is equivalent to demanding that each species be part of a non-overlapping pairing, which substantially constrains the food web structure. We demonstrate that a stable food web can always be obtained if a non-overlapping pairing exists. If it does not, the matrix rank can be used to quantify the lack of niches, corresponding to unpaired species. For the species richness at each trophic level, we derive the food web assembly rules, which specify sustainable combinations. In neighboring levels, these rules allow the higher level to avert competitive exclusion at the lower, thereby incorporating apparent competition. In agreement with data, the assembly rules predict high species numbers at intermediate levels and thinning at the top and bottom. Using comprehensive food web data, we demonstrate how omnivores or parasites with hosts at multiple trophic levels can loosen the constraints and help obtain coexistence in food webs. Hence, omnivory may be the glue that keeps communities intact even under extinction or ecological release of species.

Population Dynamics of Phage and Bacteria in Spatially Structured Habitats Using Phage λ and Escherichia coli.

UNLABELLED: Bacteria living in physically structured habitats are exposed heterogeneously to both resources and different types of phages. While there have been numerous experimental approaches to examine spatially distributed bacteria exposed to phages, there is little theory to guide the design of these experiments, interpret their results, or expand the inferences drawn to a broader ecological and evolutionary context. Plaque formation provides a window into understanding phage-bacterium interactions in physically structured populations, including surfaces, semisolids, and biofilms. We develop models to address the plaque dynamics for a temperate phage and its virulent mutants. The models are compared with phage λ-Escherichia coli system to quantify their applicability. We found that temperate phages gave an increasing number of gradually smaller colonies as the distance increased from the plaque center. For low-lysogen frequency this resulted in plaques with most of the visible colonies at an intermediate distance between the center and periphery. Using spot inoculation, where phages in excess of bacteria were inoculated in a circular area, we measured the frequency and spatial distribution of lysogens. The spot morphology of cII-negative (cII(-)) and cIII(-) mutants of phage λ displays concentric rings of high-density lysogenic colonies. The simplest of these ring morphologies was reproduced by including multiplicity of infection (MOI) sensitivity in lysis-lysogeny decisions, but its failure to explain the occasional observation of multiple rings in cIII(-) mutant phages highlights unknown features of this phage. Our findings demonstrated advantages of temperate phages over virulent phages in exploiting limited resources in spatially distributed microbial populations. IMPORTANCE: Phages are the most abundant organisms on earth, and yet little is known about how phages and bacterial hosts are influencing each other in density and evolution. Phages can be either virulent or temperate, a difference that is highlighted when a spatially structured bacterial population is infected. Phage λ is a temperate phage, with a capacity for dormancy that can be modified by single gene knockouts. The stochastic bias in the lysis-lysogeny decision's probability is reflected in plaque morphologies on bacterial lawns. We present a model for plaque morphology of both virulent and temperate phages, taking into account the underlying survival of bacterial microcolonies. It reproduces known plaque morphologies and speaks to advantages of temperate phages in a spatially structured environment.

Rapid Curtailing of the Stringent Response by Toxin-Antitoxin Module-Encoded mRNases.

UNLABELLED: Escherichia coli regulates its metabolism to adapt to changes in the environment, in particular to stressful downshifts in nutrient quality. Such shifts elicit the so-called stringent response, coordinated by the alarmone guanosine tetra- and pentaphosphate [(p)ppGpp]. On sudden amino acid (aa) starvation, RelA [(p)ppGpp synthetase I] activity is stimulated by binding of uncharged tRNAs to a vacant ribosomal site; the (p)ppGpp level increases dramatically and peaks within the time scale of a few minutes. The decrease of the (p)ppGpp level after the peak is mediated by the decreased production of mRNA by (p)ppGpp-associated transcriptional regulation, which reduces the vacant ribosomal A site and thus constitutes negative feedback to the RelA-dependent (p)ppGpp synthesis. Here we showed that on sudden isoleucine starvation, this peak was higher in an E. coli strain that lacks the 10 known mRNase-encoding toxin-antitoxin (TA) modules present in the wild-type (wt) strain. This observation suggested that toxins are part of the negative-feedback mechanism to control the (p)ppGpp level during the early stringent response. We built a ribosome trafficking model to evaluate the fold increase in RelA activity just after the onset of aa starvation. Combining this with a feedback model between the (p)ppGpp level and the mRNA level, we obtained reasonable fits to the experimental data for both strains. The analysis revealed that toxins are activated rapidly, within a minute after the onset of starvation, reducing the mRNA half-life by ∼30%. IMPORTANCE: The early stringent response elicited by amino acid starvation is controlled by a sharp increase of the cellular (p)ppGpp level. Toxin-antitoxin module-encoded mRNases are activated by (p)ppGpp through enhanced degradation of antitoxins. The present work shows that this activation happens over a very short time scale and that the activated mRNases negatively affect the (p)ppGpp level. The proposed mathematical model of (p)ppGpp regulation through the mRNA level highlights the importance of several feedback loops in early (p)ppGpp regulation.