Anti-CRISPR Phages Cooperate to Overcome CRISPR-Cas Immunity.

Some phages encode anti-CRISPR (acr) genes, which antagonize bacterial CRISPR-Cas immune systems by binding components of its machinery, but it is less clear how deployment of these acr genes impacts phage replication and epidemiology. Here, we demonstrate that bacteria with CRISPR-Cas resistance are still partially immune to Acr-encoding phage. As a consequence, Acr-phages often need to cooperate in order to overcome CRISPR resistance, with a first phage blocking the host CRISPR-Cas immune system to allow a second Acr-phage to successfully replicate. This cooperation leads to epidemiological tipping points in which the initial density of Acr-phage tips the balance from phage extinction to a phage epidemic. Furthermore, both higher levels of CRISPR-Cas immunity and weaker Acr activities shift the tipping points toward higher initial phage densities. Collectively, these data help elucidate how interactions between phage-encoded immune suppressors and the CRISPR systems they target shape bacteria-phage population dynamics.

Six distinct NFκB signaling codons convey discrete information to distinguish stimuli and enable appropriate macrophage responses.

Macrophages initiate inflammatory responses via the transcription factor NFκB. The temporal pattern of NFκB activity determines which genes are expressed and thus, the type of response that ensues. Here, we examined how information about the stimulus is encoded in the dynamics of NFκB activity. We generated an mVenus-RelA reporter mouse line to enable high-throughput live-cell analysis of primary macrophages responding to host- and pathogen-derived stimuli. An information-theoretic workflow identified six dynamical features-termed signaling codons-that convey stimulus information to the nucleus. In particular, oscillatory trajectories were a hallmark of responses to cytokine but not pathogen-derived stimuli. Single-cell imaging and RNA sequencing of macrophages from a mouse model of Sjögren's syndrome revealed inappropriate responses to stimuli, suggestive of confusion of two NFκB signaling codons. Thus, the dynamics of NFκB signaling classify immune threats through six signaling codons, and signal confusion based on defective codon deployment may underlie the etiology of some inflammatory diseases.

Differential Expression of the Transcription Factor GATA3 Specifies Lineage and Functions of Innate Lymphoid Cells.

Lymphoid tissue inducer (LTi) cells are regarded as a subset of innate lymphoid cells (ILCs). However, these cells are not derived from the ILC common progenitor, which generates other ILC subsets and is defined by the expression of the transcription factor PLZF. Here, we examined transcription factor(s) determining the fate of LTi progenitors versus non-LTi ILC progenitors. Conditional deletion of Gata3 resulted in the loss of PLZF+ non-LTi progenitors but not the LTi progenitors that expressed the transcription factor RORγt. Consistently, PLZF+ non-LTi progenitors expressed high amounts of GATA3, whereas GATA3 expression was low in RORγt+ LTi progenitors. The generation of both progenitors required the transcriptional regulator Id2, which defines the common helper-like innate lymphoid progenitor (ChILP), but not cytokine signaling. Nevertheless, low GATA3 expression was necessary for the generation of functionally mature LTi cells. Thus, differential expression of GATA3 determines the fates and functions of distinct ILC progenitors.

Resolving Cell Fate Decisions during Somatic Cell Reprogramming by Single-Cell RNA-Seq.

Somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs), which is a highly heterogeneous process. Here we report the cell fate continuum during somatic cell reprogramming at single-cell resolution. We first develop SOT to analyze cell fate continuum from Oct4/Sox2/Klf4- or OSK-mediated reprogramming and show that cells bifurcate into two categories, reprogramming potential (RP) or non-reprogramming (NR). We further show that Klf4 contributes to Cd34+/Fxyd5+/Psca+ keratinocyte-like NR fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the RP cells. We analyze more than 150,000 single cells from both OSK and chemical reprograming and identify additional NR/RP bifurcation points. Our work reveals a generic bifurcation model for cell fate decisions during somatic cell reprogramming that may be applicable to other systems and inspire further improvements for reprogramming.

A dynamical systems treatment of transcriptomic trajectories in hematopoiesis.

Inspired by Waddington's illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein exogenous signaling dynamics couple to the enormously complex signaling and transcriptional machinery of a cell to elicit qualitative transitions in its collective state. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the variations of a myriad of genes. Here, we present a mathematical formalism for rigorously evaluating, from a dynamical systems perspective, whether scRNA-seq trajectories display statistical signatures consistent with bifurcations and, as a case study, pinpoint regions of multistability along the neutrophil branch of hematopoeitic differentiation. Additionally, we leverage the geometric features of linear instability to identify the low-dimensional phase plane in gene expression space within which the multistability unfolds, highlighting novel genetic players that are crucial for neutrophil differentiation. Broadly, we show that a dynamical systems treatment of scRNA-seq data provides mechanistic insights into the high-dimensional processes of cellular differentiation, taking a step toward systematic construction of mathematical models for transcriptomic dynamics.

Diverse branching forms regulated by a core auxin transport mechanism in plants.

Diverse branching forms have evolved multiple times across the tree of life to facilitate resource acquisition and exchange with the environment. In the vascular plant group, the ancestral pattern of branching involves dichotomy of a parent shoot apex to form two new daughter apices. The molecular basis of axillary branching in Arabidopsis is well understood, but few regulators of dichotomous branching are known. Through analyses of dichotomous branching in the lycophyte, Selaginella kraussiana, we identify PIN-mediated auxin transport as an ancestral branch regulator of vascular plants. We show that short-range auxin transport out of the apices promotes dichotomy and that branch dominance is globally coordinated by long-range auxin transport. Uniquely in Selaginella, angle meristems initiate at each dichotomy, and these can develop into rhizophores or branching angle shoots. We show that long-range auxin transport and a transitory drop in PIN expression are involved in angle shoot development. We conclude that PIN-mediated auxin transport is an ancestral mechanism for vascular plant branching that was independently recruited into Selaginella angle shoot development and seed plant axillary branching during evolution.

Bifurcation instructed design of multistate machines.

We propose a design paradigm for multistate machines where transitions from one state to another are organized by bifurcations of multiple equilibria of the energy landscape describing the collective interactions of the machine components. This design paradigm is attractive since, near bifurcations, small variations in a few control parameters can result in large changes to the system's state providing an emergent lever mechanism. Further, the topological configuration of transitions between states near such bifurcations ensures robust operation, making the machine less sensitive to fabrication errors and noise. To design such machines, we develop and implement a new efficient algorithm that searches for interactions between the machine components that give rise to energy landscapes with these bifurcation structures. We demonstrate a proof of concept for this approach by designing magnetoelastic machines whose motions are primarily guided by their magnetic energy landscapes and show that by operating near bifurcations we can achieve multiple transition pathways between states. This proof of concept demonstration illustrates the power of this approach, which could be especially useful for soft robotics and at the microscale where typical macroscale designs are difficult to implement.

Functionally redundant formate dehydrogenases enable formate-dependent growth in Methanococcus maripaludis.

Methanogens are essential for the complete remineralization of organic matter in anoxic environments. Most cultured methanogens are hydrogenotrophic, using H2 as an electron donor to reduce CO2 to CH4, but in the absence of H2 many can also use formate. Formate dehydrogenase (Fdh) is essential for formate oxidation, where it transfers electrons for the reduction of coenzyme F420 or to a flavin-based electron bifurcating reaction catalyzed by heterodisulfide reductase (Hdr), the terminal reaction of methanogenesis. Furthermore, methanogens that use formate encode at least two isoforms of Fdh in their genomes, but how these different isoforms participate in methanogenesis is unknown. Using Methanococcus maripaludis, we undertook a biochemical characterization of both Fdh isoforms involved in methanogenesis. Both Fdh1 and Fdh2 interacted with Hdr to catalyze the flavin-based electron bifurcating reaction, and both reduced F420 at similar rates. F420 reduction preceded flavin-based electron bifurcation activity for both enzymes. In a Δfdh1 mutant background, a suppressor mutation was required for Fdh2 activity. Genome sequencing revealed that this mutation resulted in the loss of a specific molybdopterin transferase (moeA), allowing for Fdh2-dependent growth, and the metal content of the proteins suggested that isoforms are dependent on either molybdenum or tungsten for activity. These data suggest that both isoforms of Fdh are functionally redundant, but their activities in vivo may be limited by gene regulation or metal availability under different growth conditions. Together these results expand our understanding of formate oxidation and the role of Fdh in methanogenesis.

The rapid-reaction kinetics of an electron-bifurcating flavoprotein, the crotonyl-CoA-dependent NADH:ferredoxin oxidoreductase EtfAB:bcd.

We have investigated the kinetic behavior of the electron-bifurcating crotonyl-CoA-dependent NADH: ferredoxin oxidoreductase EtfAB:bcd from Megasphaera elsdenii. The overall behavior of the complex in both the reductive and the oxidative half-reactions is consistent with that previously determined for the individual EtfAB and bcd components. This includes an uncrossing of the half-potentials of the bifurcating flavin of the EtfAB component in the course of ferredoxin-reducing catalysis, ionization of the bcd flavin semiquinone and the appearance of a charge transfer complex upon binding of the high potential acceptor crotonyl-CoA. The observed rapid-reaction rates of ferredoxin reduction are independent of [NADH], [crotonyl-CoA], and [ferredoxin], with an observed rate of ∼0.2 s-1, consistent with the observed steady-state kinetics. In enzyme-monitored turnover experiments, an approach to steady-state where the complex's flavins become reduced but no ferredoxin is generated is followed by a steady-state phase characterized by extensive ferredoxin reduction but little change in overall levels of flavin reduction. The approach to the steady-state phase can be eliminated by prior reduction of the complex, in which case there is no lag in the onset of ferredoxin reduction; this is consistent with the et FAD needing to be reduced to the level of the (anionic) semiquinone for bifurcation and concomitant ferredoxin reduction to occur. Single-turnover experiments support this conclusion, with the accumulation of the anionic semiquinone of the et FAD apparently required to prime the system for subsequent bifurcation and ferredoxin reduction.

Collective cell migration is spatiotemporally regulated during mammary epithelial bifurcation.

Branched epithelial networks are generated through an iterative process of elongation and bifurcation. We sought to understand bifurcation of the mammary epithelium. To visualize this process, we utilized three-dimensional (3D) organotypic culture and time-lapse confocal microscopy. We tracked cell migration during bifurcation and observed local reductions in cell speed at the nascent bifurcation cleft. This effect was proximity dependent, as individual cells approaching the cleft reduced speed, whereas cells exiting the cleft increased speed. As the cells slow down, they orient both migration and protrusions towards the nascent cleft, while cells in the adjacent branches orient towards the elongating tips. We next tested the hypothesis that TGF-ß signaling controls mammary branching by regulating cell migration. We first validated that addition of TGF-ß1 (TGFB1) protein increased cleft number, whereas inhibition of TGF-ß signaling reduced cleft number. Then, consistent with our hypothesis, we observed that pharmacological inhibition of TGF-ß1 signaling acutely decreased epithelial migration speed. Our data suggest a model for mammary epithelial bifurcation in which TGF-ß signaling regulates cell migration to determine the local sites of bifurcation and the global pattern of the tubular network.