Live-cell imaging reveals single-cell and population-level infection strategies of Listeria monocytogenes in macrophages.

Pathogens have developed intricate strategies to overcome the host's innate immune responses. In this paper we use live-cell microscopy with a single bacterium resolution to follow in real time interactions between the food-borne pathogen L. monocytogenes and host macrophages, a key event controlling the infection in vivo. We demonstrate that infection results in heterogeneous outcomes, with only a subset of bacteria able to establish a replicative invasion of macrophages. The fate of individual bacteria in the same host cell was independent from the host cell and non-cooperative, being independent from co-infecting bacteria. A higher multiplicity of infection resulted in a reduced probability of replication of the overall bacterial population. By use of internalisation assays and conditional probabilities to mathematically describe the two-stage invasion process, we demonstrate that the higher MOI compromises the ability of macrophages to phagocytose bacteria. We found that the rate of phagocytosis is mediated via the secreted Listeriolysin toxin (LLO), while the probability of replication of intracellular bacteria remained constant. Using strains expressing fluorescent reporters to follow transcription of either the LLO-encoding hly or actA genes, we show that replicative bacteria exhibited higher PrfA regulon expression in comparison to those bacteria that did not replicate, however elevated PrfA expression per se was not sufficient to increase the probability of replication. Overall, this demonstrates a new role for the population-level, but not single cell, PrfA-mediated activity to regulate outcomes of host pathogen interactions.

Overexpression of IκB⍺ modulates NF-κB activation of inflammatory target gene expression.

Cells respond to inflammatory stimuli such as cytokines by activation of the nuclear factor-κB (NF-κB) signalling pathway, resulting in oscillatory translocation of the transcription factor p65 between nucleus and cytoplasm in some cell types. We investigate the relationship between p65 and inhibitor-κB⍺ (IκBα) protein levels and dynamic properties of the system, and how this interaction impacts on the expression of key inflammatory genes. Using bacterial artificial chromosomes, we developed new cell models of IκB⍺-eGFP protein overexpression in a pseudo-native genomic context. We find that cells with high levels of the negative regulator IκBα remain responsive to inflammatory stimuli and maintain dynamics for both p65 and IκBα. In contrast, canonical target gene expression is dramatically reduced by overexpression of IκBα, but can be partially rescued by overexpression of p65. Treatment with leptomycin B to promote nuclear accumulation of IκB⍺ also suppresses canonical target gene expression, suggesting a mechanism in which nuclear IκB⍺ accumulation prevents productive p65 interaction with promoter binding sites. This causes reduced target promoter binding and gene transcription, which we validate by chromatin immunoprecipitation and in primary cells. Overall, we show how inflammatory gene transcription is modulated by the expression levels of both IκB⍺ and p65. This results in an anti-inflammatory effect on transcription, demonstrating a broad mechanism to modulate the strength of inflammatory response.

Post-transcriptional regulatory feedback encodes JAK-STAT signal memory of interferon stimulation.

Immune cells fine tune their responses to infection and inflammatory cues. Here, using live-cell confocal microscopy and mathematical modelling, we investigate interferon-induced JAK-STAT signalling in innate immune macrophages. We demonstrate that transient exposure to IFN-γ stimulation induces a long-term desensitisation of STAT1 signalling and gene expression responses, revealing a dose- and time-dependent regulatory feedback that controls JAK-STAT responses upon re-exposure to stimulus. We show that IFN-α/ß1 elicit different level of desensitisation from IFN-γ, where cells refractory to IFN-α/ß1 are sensitive to IFN-γ, but not vice versa. We experimentally demonstrate that the underlying feedback mechanism involves regulation of STAT1 phosphorylation but is independent of new mRNA synthesis and cognate receptor expression. A new feedback model of the protein tyrosine phosphatase activity recapitulates experimental data and demonstrates JAK-STAT network's ability to decode relative changes of dose, timing, and type of temporal interferon stimulation. These findings reveal that STAT desensitisation renders cells with signalling memory of type I and II interferon stimulation, which in the future may improve administration of interferon therapy.

Quantitative Characterisation of Low Abundant Yeast Mitochondrial Proteins Reveals Compensation for Haplo-Insufficiency in Different Environments.

The quantification of low abundant membrane-binding proteins such as transcriptional factors and chaperones has proven difficult, even with the most sophisticated analytical technologies. Here, we exploit and optimise the non-invasive Fluorescence Correlation Spectroscopy (FCS) for the quantitation of low abundance proteins, and as proof of principle, we choose two interacting proteins involved in the fission of mitochondria in yeast, Fis1p and Mdv1p. In Saccharomyces cerevisiae, the recruitment of Fis1p and Mdv1p to mitochondria is essential for the scission of the organelles and the retention of functional mitochondrial structures in the cell. We use FCS in single GFP-labelled live yeast cells to quantify the protein abundance in homozygote and heterozygote cells and to investigate the impact of the environments on protein copy number, bound/unbound protein state and mobility kinetics. Both proteins were observed to localise predominantly at mitochondrial structures, with the Mdv1p bound state increasing significantly in a strictly respiratory environment. Moreover, a compensatory mechanism that controls Fis1p abundance upon deletion of one allele was observed in Fis1p but not in Mdv1p, suggesting differential regulation of Fis1p and Mdv1p protein expression.

Quantification of protein abundance and interaction defines a mechanism for operation of the circadian clock.

The mammalian circadian clock exerts control of daily gene expression through cycles of DNA binding. Here, we develop a quantitative model of how a finite pool of BMAL1 protein can regulate thousands of target sites over daily time scales. We used quantitative imaging to track dynamic changes in endogenous labelled proteins across peripheral tissues and the SCN. We determine the contribution of multiple rhythmic processes coordinating BMAL1 DNA binding, including cycling molecular abundance, binding affinities, and repression. We find nuclear BMAL1 concentration determines corresponding CLOCK through heterodimerisation and define a DNA residence time of this complex. Repression of CLOCK:BMAL1 is achieved through rhythmic changes to BMAL1:CRY1 association and high-affinity interactions between PER2:CRY1 which mediates CLOCK:BMAL1 displacement from DNA. Finally, stochastic modelling reveals a dual role for PER:CRY complexes in which increasing concentrations of PER2:CRY1 promotes removal of BMAL1:CLOCK from genes consequently enhancing ability to move to new target sites.

Cryptochrome proteins regulate the circadian intracellular behavior and localization of PER2 in mouse suprachiasmatic nucleus neurons.

The ∼20,000 cells of the suprachiasmatic nucleus (SCN), the master circadian clock of the mammalian brain, coordinate subordinate cellular clocks across the organism, driving adaptive daily rhythms of physiology and behavior. The canonical model for SCN timekeeping pivots around transcriptional/translational feedback loops (TTFL) whereby PERIOD (PER) and CRYPTOCHROME (CRY) clock proteins associate and translocate to the nucleus to inhibit their own expression. The fundamental individual and interactive behaviors of PER and CRY in the SCN cellular environment and the mechanisms that regulate them are poorly understood. We therefore used confocal imaging to explore the behavior of endogenous PER2 in the SCN of PER2::Venus reporter mice, transduced with viral vectors expressing various forms of CRY1 and CRY2. In contrast to nuclear localization in wild-type SCN, in the absence of CRY proteins, PER2 was predominantly cytoplasmic and more mobile, as measured by fluorescence recovery after photobleaching. Virally expressed CRY1 or CRY2 relocalized PER2 to the nucleus, initiated SCN circadian rhythms, and determined their period. We used translational switching to control CRY1 cellular abundance and found that low levels of CRY1 resulted in minimal relocalization of PER2, but yet, remarkably, were sufficient to initiate and maintain circadian rhythmicity. Importantly, the C-terminal tail was necessary for CRY1 to localize PER2 to the nucleus and to initiate SCN rhythms. In CRY1-null SCN, CRY1Δtail opposed PER2 nuclear localization and correspondingly shortened SCN period. Through manipulation of CRY proteins, we have obtained insights into the spatiotemporal behaviors of PER and CRY sitting at the heart of the TTFL molecular mechanism.

Gene-Specific Linear Trends Constrain Transcriptional Variability of the Toll-like Receptor Signaling.

Single-cell gene expression is inherently variable, but how this variability is controlled in response to stimulation remains unclear. Here, we use single-cell RNA-seq and single-molecule mRNA counting (smFISH) to study inducible gene expression in the immune toll-like receptor system. We show that mRNA counts of tumor necrosis factor α conform to a standard stochastic switch model, while transcription of interleukin-1ß involves an additional regulatory step resulting in increased heterogeneity. Despite different modes of regulation, systematic analysis of single-cell data for a range of genes demonstrates that the variability in transcript count is linearly constrained by the mean response over a range of conditions. Mathematical modeling of smFISH counts and experimental perturbation of chromatin state demonstrates that linear constraints emerge through modulation of transcriptional bursting along with gene-specific relationships. Overall, our analyses demonstrate that the variability of the inducible single-cell mRNA response is constrained by transcriptional bursting.

Quantitative live imaging of Venus::BMAL1 in a mouse model reveals complex dynamics of the master circadian clock regulator.

Evolutionarily conserved circadian clocks generate 24-hour rhythms in physiology and behaviour that adapt organisms to their daily and seasonal environments. In mammals, the suprachiasmatic nucleus (SCN) of the hypothalamus is the principal co-ordinator of the cell-autonomous clocks distributed across all major tissues. The importance of robust daily rhythms is highlighted by experimental and epidemiological associations between circadian disruption and human diseases. BMAL1 (a bHLH-PAS domain-containing transcription factor) is the master positive regulator within the transcriptional-translational feedback loops (TTFLs) that cell-autonomously define circadian time. It drives transcription of the negative regulators Period and Cryptochrome alongside numerous clock output genes, and thereby powers circadian time-keeping. Because deletion of Bmal1 alone is sufficient to eliminate circadian rhythms in cells and the whole animal it has been widely used as a model for molecular disruption of circadian rhythms, revealing essential, tissue-specific roles of BMAL1 in, for example, the brain, liver and the musculoskeletal system. Moreover, BMAL1 has clock-independent functions that influence ageing and protein translation. Despite the essential role of BMAL1 in circadian time-keeping, direct measures of its intra-cellular behaviour are still lacking. To fill this knowledge-gap, we used CRISPR Cas9 to generate a mouse expressing a knock-in fluorescent fusion of endogenous BMAL1 protein (Venus::BMAL1) for quantitative live imaging in physiological settings. The Bmal1Venus mouse model enabled us to visualise and quantify the daily behaviour of this core clock factor in central (SCN) and peripheral clocks, with single-cell resolution that revealed its circadian expression, anti-phasic to negative regulators, nuclear-cytoplasmic mobility and molecular abundance.

Heat shock response regulates stimulus-specificity and sensitivity of the pro-inflammatory NF-κB signalling.

BACKGROUND: Ability to adapt to temperature changes trough the Heat Shock Response (HSR) pathways is one of the most fundamental and clinically relevant cellular response systems. Heat Shock (HS) affects the signalling and gene expression responses of the Nuclear Factor κB (NF-κB) transcription factor, a critical regulator of proliferation and inflammation, however, our quantitative understanding of how cells sense and adapt to temperature changes is limited. METHODS: We used live-cell time-lapse microscopy and mathematical modelling to understand the signalling of the NF-κB system in the human MCF7 breast adenocarcinoma cells in response to pro-inflammatory Interleukin 1ß (IL1ß) and Tumour Necrosis Factor α (TNFα) cytokines, following exposure to a 37-43 °C range of physiological and clinical temperatures. RESULTS: We show that exposure to 43 °C 1 h HS inhibits the immediate NF-κB signalling response to TNFα and IL1ß stimulation although uptake of cytokines is not impaired. Within 4 h after HS treatment IL1ß-induced NF-κB responses return to normal levels, but the recovery of the TNFα-induced responses is still affected. Using siRNA knock-down of Heat Shock Factor 1 (HSF1) we show that this stimulus-specificity is conferred via the Inhibitory κB kinase (IKK) signalosome where HSF1-dependent feedback regulates TNFα, but not IL1ß-mediated IKK recovery post HS. Furthermore, we demonstrate that through the temperature-dependent denaturation and recovery of IKK, TNFα and IL1ß-mediated signalling exhibit different temperature sensitivity and adaptation to repeated HS when exposed to a 37-43 °C temperature range. Specifically, IL1ß-mediated NF-κB responses are more robust to temperature changes in comparison to those induced by TNFα treatment. CONCLUSIONS: We demonstrate that the kinetics of the NF-κB system following temperature stress is cytokine specific and exhibit differential adaptation to temperature changes. We propose that this differential temperature sensitivity is mediated via the IKK signalosome, which acts as a bona fide temperature sensor trough the HSR cross-talk. This novel quantitative understanding of NF-κB and HSR interactions is fundamentally important for the potential optimization of therapeutic hyperthermia protocols. Video Abstract.

The clock gene Bmal1 inhibits macrophage motility, phagocytosis, and impairs defense against pneumonia.

The circadian clock regulates many aspects of immunity. Bacterial infections are affected by time of day, but the mechanisms involved remain undefined. Here we show that loss of the core clock protein BMAL1 in macrophages confers protection against pneumococcal pneumonia. Infected mice show both reduced weight loss and lower bacterial burden in circulating blood. In vivo studies of macrophage phagocytosis reveal increased bacterial ingestion following Bmal1 deletion, which was also seen in vitro. BMAL1-/- macrophages exhibited marked differences in actin cytoskeletal organization, a phosphoproteome enriched for cytoskeletal changes, with reduced phosphocofilin and increased active RhoA. Further analysis of the BMAL1-/- macrophages identified altered cell morphology and increased motility. Mechanistically, BMAL1 regulated a network of cell movement genes, 148 of which were within 100 kb of high-confidence BMAL1 binding sites. Links to RhoA function were identified, with 29 genes impacting RhoA expression or activation. RhoA inhibition restored the phagocytic phenotype to that seen in control macrophages. In summary, we identify a surprising gain of antibacterial function due to loss of BMAL1 in macrophages, associated with a RhoA-dependent cytoskeletal change, an increase in cell motility, and gain of phagocytic function.

The circadian clock protein REVERBα inhibits pulmonary fibrosis development.

Pulmonary inflammatory responses lie under circadian control; however, the importance of circadian mechanisms in the underlying fibrotic phenotype is not understood. Here, we identify a striking change to these mechanisms resulting in a gain of amplitude and lack of synchrony within pulmonary fibrotic tissue. These changes result from an infiltration of mesenchymal cells, an important cell type in the pathogenesis of pulmonary fibrosis. Mutation of the core clock protein REVERBα in these cells exacerbated the development of bleomycin-induced fibrosis, whereas mutation of REVERBα in club or myeloid cells had no effect on the bleomycin phenotype. Knockdown of REVERBα revealed regulation of the little-understood transcription factor TBPL1. Both REVERBα and TBPL1 altered integrinß1 focal-adhesion formation, resulting in increased myofibroblast activation. The translational importance of our findings was established through analysis of 2 human cohorts. In the UK Biobank, circadian strain markers (sleep length, chronotype, and shift work) are associated with pulmonary fibrosis, making them risk factors. In a separate cohort, REVERBα expression was increased in human idiopathic pulmonary fibrosis (IPF) lung tissue. Pharmacological targeting of REVERBα inhibited myofibroblast activation in IPF fibroblasts and collagen secretion in organotypic cultures from IPF patients, thus suggesting that targeting of REVERBα could be a viable therapeutic approach.

Quantitative analysis of competitive cytokine signaling predicts tissue thresholds for the propagation of macrophage activation.

Toll-like receptor (TLR) signaling regulates macrophage activation and effector cytokine propagation in the constrained environment of a tissue. In macrophage populations, TLR4 stimulates the dose-dependent transcription of nuclear factor κB (NF-κB) target genes. However, using single-RNA counting, we found that individual cells exhibited a wide range (three orders of magnitude) of expression of the gene encoding the proinflammatory cytokine tumor necrosis factor-α (TNF-α). The TLR4-induced TNFA transcriptional response correlated with the extent of NF-κB signaling in the cells and their size. We compared the rates of TNF-α production and uptake in macrophages and mouse embryonic fibroblasts and generated a mathematical model to explore the heterogeneity in the response of macrophages to TLR4 stimulation and the propagation of the TNF-α signal in the tissue. The model predicts that the local propagation of the TLR4-dependent TNF-α response and cellular NF-κB signaling are limited to small distances of a few cell diameters between neighboring tissue-resident macrophages. In our predictive model, TNF-α propagation was constrained by competitive uptake of TNF-α from the environment, rather than by heterogeneous production of the cytokine. We propose that the highly constrained architecture of tissues enables effective localized propagation of inflammatory cues while avoiding out-of-context responses at longer distances.

Integration of Kinase and Calcium Signaling at the Level of Chromatin Underlies Inducible Gene Activation in T Cells.

TCR signaling pathways cooperate to activate the inducible transcription factors NF-κB, NFAT, and AP-1. In this study, using the calcium ionophore ionomycin and/or PMA on Jurkat T cells, we show that the gene expression program associated with activation of TCR signaling is closely related to specific chromatin landscapes. We find that calcium and kinase signaling cooperate to induce chromatin remodeling at ∼2100 chromatin regions, which demonstrate enriched binding motifs for inducible factors and correlate with target gene expression. We found that these regions typically function as inducible enhancers. Many of these elements contain composite NFAT/AP-1 sites, which typically support cooperative binding, thus further reinforcing the need for cooperation between calcium and kinase signaling in the activation of genes in T cells. In contrast, treatment with PMA or ionomycin alone induces chromatin remodeling at far fewer regions (∼600 and ∼350, respectively), which mostly represent a subset of those induced by costimulation. This suggests that the integration of TCR signaling largely occurs at the level of chromatin, which we propose plays a crucial role in regulating T cell activation.

Stochasticity in the miR-9/Hes1 oscillatory network can account for clonal heterogeneity in the timing of differentiation.

Recent studies suggest that cells make stochastic choices with respect to differentiation or division. However, the molecular mechanism underlying such stochasticity is unknown. We previously proposed that the timing of vertebrate neuronal differentiation is regulated by molecular oscillations of a transcriptional repressor, HES1, tuned by a post-transcriptional repressor, miR-9. Here, we computationally model the effects of intrinsic noise on the Hes1/miR-9 oscillator as a consequence of low molecular numbers of interacting species, determined experimentally. We report that increased stochasticity spreads the timing of differentiation in a population, such that initially equivalent cells differentiate over a period of time. Surprisingly, inherent stochasticity also increases the robustness of the progenitor state and lessens the impact of unequal, random distribution of molecules at cell division on the temporal spread of differentiation at the population level. This advantageous use of biological noise contrasts with the view that noise needs to be counteracted.

Fenamate NSAIDs inhibit the NLRP3 inflammasome and protect against Alzheimer's disease in rodent models.

Non-steroidal anti-inflammatory drugs (NSAIDs) inhibit cyclooxygenase-1 (COX-1) and COX-2 enzymes. The NLRP3 inflammasome is a multi-protein complex responsible for the processing of the proinflammatory cytokine interleukin-1ß and is implicated in many inflammatory diseases. Here we show that several clinically approved and widely used NSAIDs of the fenamate class are effective and selective inhibitors of the NLRP3 inflammasome via inhibition of the volume-regulated anion channel in macrophages, independently of COX enzymes. Flufenamic acid and mefenamic acid are efficacious in NLRP3-dependent rodent models of inflammation in air pouch and peritoneum. We also show therapeutic effects of fenamates using a model of amyloid beta induced memory loss and a transgenic mouse model of Alzheimer's disease. These data suggest that fenamate NSAIDs could be repurposed as NLRP3 inflammasome inhibitors and Alzheimer's disease therapeutics.

Signal transduction controls heterogeneous NF-κB dynamics and target gene expression through cytokine-specific refractory states.

Cells respond dynamically to pulsatile cytokine stimulation. Here we report that single, or well-spaced pulses of TNFα (>100 min apart) give a high probability of NF-κB activation. However, fewer cells respond to shorter pulse intervals (<100 min) suggesting a heterogeneous refractory state. This refractory state is established in the signal transduction network downstream of TNFR and upstream of IKK, and depends on the level of the NF-κB system negative feedback protein A20. If a second pulse within the refractory phase is IL-1ß instead of TNFα, all of the cells respond. This suggests a mechanism by which two cytokines can synergistically activate an inflammatory response. Gene expression analyses show strong correlation between the cellular dynamic response and NF-κB-dependent target gene activation. These data suggest that refractory states in the NF-κB system constitute an inherent design motif of the inflammatory response and we suggest that this may avoid harmful homogenous cellular activation.

Investigating IL-1ß Secretion Using Real-Time Single-Cell Imaging.

The pro-inflammatory cytokine interleukin (IL)-1ß is an important mediator of the inflammatory response. In order to perform its role in the inflammatory cascade, IL-1ß must be secreted from the cell, yet it lacks a signal peptide that is required for conventional secretion, and the exact mechanism of release remains undefined. Conventional biochemical methods have limited the investigation into the processes involved in IL-1ß secretion to population dynamics, yet heterogeneity between cells has been observed at a single-cell level. Here, greater sensitivity is achieved with the use of a newly developed vector that codes for a fluorescently labelled version of IL-1ß. Combining this with real-time single-cell confocal microscopy using the methods described here, we have developed an effective protocol for investigating the mechanisms of IL-1ß secretion and the testing of the hypothesis that IL-1ß secretion requires membrane permeabilisation.

Modeling the dynamics of hypoxia inducible factor-1α (HIF-1α) within single cells and 3D cell culture systems.

HIF (hypoxia inducible factor) is an oxygen-regulated transcription factor that mediates the intracellular response to hypoxia in human cells. There is increasing evidence that cell signaling pathways encode temporal information, and thus cell fate may be determined by the dynamics of protein levels. We have developed a mathematical model to describe the transient dynamics of the HIF-1α protein measured in single cells subjected to hypoxic shock. The essential characteristics of these data are modeled with a system of differential equations describing the feedback inhibition between HIF-1α and prolyl hydroxylases (PHD) oxygen sensors. Heterogeneity in the single-cell data is accounted through parameter variation in the model. We previously identified the PHD2 isoform as the main PHD sensor responsible for controlling the HIF-1α transient response, and make here testable predictions regarding HIF-1α dynamics subject to repetitive hypoxic pulses. The model is further developed to describe the dynamics of HIF-1α in cells cultured as 3D spheroids, with oxygen dynamics parameterized using experimental measurements of oxygen within spheroids. We show that the dynamics of HIF-1α and transcriptional targets of HIF-1α display a non-monotone response to the oxygen dynamics. Specifically we demonstrate that the dynamic transient behavior of HIF-1α results in differential dynamics in transcriptional targets.

Tight control of hypoxia-inducible factor-α transient dynamics is essential for cell survival in hypoxia.

Intracellular signaling involving hypoxia-inducible factor (HIF) controls the adaptive responses to hypoxia. There is a growing body of evidence demonstrating that intracellular signals encode temporal information. Thus, the dynamics of protein levels, as well as protein quantity and/or localization, impacts on cell fate. We hypothesized that such temporal encoding has a role in HIF signaling and cell fate decisions triggered by hypoxic conditions. Using live cell imaging in a controlled oxygen environment, we observed transient 3-h pulses of HIF-1α and -2α expression under continuous hypoxia. We postulated that the well described prolyl hydroxylase (PHD) oxygen sensors and HIF negative feedback regulators could be the origin of the pulsatile HIF dynamics. We used iterative mathematical modeling and experimental analysis to scrutinize which parameter of the PHD feedback could control HIF timing and we probed for the functional redundancy between the three main PHD proteins. We identified PHD2 as the main PHD responsible for HIF peak duration. We then demonstrated that this has important consequences, because the transient nature of the HIF pulse prevents cell death by avoiding transcription of p53-dependent pro-apoptotic genes. We have further shown the importance of considering HIF dynamics for coupling mathematical models by using a described HIF-p53 mathematical model. Our results indicate that the tight control of HIF transient dynamics has important functional consequences on the cross-talk with key signaling pathways controlling cell survival, which is likely to impact on HIF targeting strategies for hypoxia-associated diseases such as tumor progression and ischemia.