Extended methods for spatial cell classification with DBSCAN-CellX.

Local cell densities and positioning within cellular monolayers and stratified epithelia have important implications for cell interactions and the functionality of various biological processes. To analyze the relationship between cell localization and tissue physiology, density-based clustering algorithms, such as DBSCAN, allow for a detailed characterization of the spatial distribution and positioning of individual cells. However, these methods rely on predefined parameters that influence the outcome of the analysis. With varying cell densities in cell cultures or tissues impacting cell sizes and, thus, cellular proximities, these parameters need to be carefully chosen. In addition, standard DBSCAN approaches generally come short in appropriately identifying individual cell positions. We therefore developed three extensions to the standard DBSCAN-algorithm that provide: (i) an automated parameter identification to reliably identify cell clusters, (ii) an improved identification of cluster edges; and (iii) an improved characterization of the relative positioning of cells within clusters. We apply our novel methods, which are provided as a user-friendly OpenSource-software package (DBSCAN-CellX), to cellular monolayers of different cell lines. Thereby, we show the importance of the developed extensions for the appropriate analysis of cell culture experiments to determine the relationship between cell localization and tissue physiology.

Adaptation to host cell environment during experimental evolution of Zika virus.

Zika virus (ZIKV) infection can cause important developmental and neurological defects in Humans. Type I/III interferon responses control ZIKV infection and pathological processes, yet the virus has evolved various mechanisms to defeat these host responses. Here, we established a pipeline to delineate at high-resolution the genetic evolution of ZIKV in a controlled host cell environment. We uncovered that serially passaged ZIKV acquired increased infectivity and simultaneously developed a resistance to TLR3-induced restriction. We built a mathematical model that suggests that the increased infectivity is due to a reduced time-lag between infection and viral replication. We found that this adaptation is cell-type specific, suggesting that different cell environments may drive viral evolution along different routes. Deep-sequencing of ZIKV populations pinpointed mutations whose increased frequencies temporally coincide with the acquisition of the adapted phenotype. We functionally validated S455L, a substitution in ZIKV envelope (E) protein, recapitulating the adapted phenotype. Its positioning on the E structure suggests a putative function in protein refolding/stability. Taken together, our results uncovered ZIKV adaptations to the cellular environment leading to accelerated replication onset coupled with resistance to TLR3-induced antiviral response. Our work provides insights into Zika virus adaptation to host cells and immune escape mechanisms.

Temporal control of the integrated stress response by a stochastic molecular switch.

Stress granules (SGs) are formed in the cytosol as an acute response to environmental cues and activation of the integrated stress response (ISR), a central signaling pathway controlling protein synthesis. Using chronic virus infection as stress model, we previously uncovered a unique temporal control of the ISR resulting in recurrent phases of SG assembly and disassembly. Here, we elucidate the molecular network generating this fluctuating stress response by integrating quantitative experiments with mathematical modeling and find that the ISR operates as a stochastic switch. Key elements controlling this switch are the cooperative activation of the stress-sensing kinase PKR, the ultrasensitive response of SG formation to the phosphorylation of the translation initiation factor eIF2α, and negative feedback via GADD34, a stress-induced subunit of protein phosphatase 1. We identify GADD34 messenger RNA levels as the molecular memory of the ISR that plays a central role in cell adaptation to acute and chronic stress.

Dengue virus is sensitive to inhibition prior to productive replication.

Uncovering vulnerable steps in the life cycle of viruses supports the rational design of antiviral treatments. However, information on viral replication dynamics obtained from traditional bulk assays with host cell populations is inherently limited as the data represent averages over a multitude of unsynchronized replication cycles. Here, we use time-lapse imaging of virus replication in thousands of single cells, combined with computational inference, to identify rate-limiting steps for dengue virus (DENV), a widespread human pathogen. Comparing wild-type DENV with a vaccine candidate mutant, we show that the viral spread in the mutant is greatly attenuated by delayed onset of productive replication, whereas wild-type and mutant virus have identical replication rates. Single-cell analysis done after applying the broad-spectrum antiviral drug, ribavirin, at clinically relevant concentrations revealed the same mechanism of attenuating viral spread. We conclude that the initial steps of infection, rather than the rate of established replication, are quantitatively limiting DENV spread.

Antiviral interferon response at single-cell resolution.

Experimental studies of the innate immune response of mammalian cells to viruses reveal pervasive heterogeneity at the level of single cells. Interferons are induced only in a fraction of virus-infected cells; subsequently a fraction of cells exposed to interferons upregulate interferon-stimulated genes. Nevertheless, quantitative experiments and linked mathematical models show that the interferon response can be effective in curbing viral spread through two distinct mechanisms. First, paracrine interferon signals from scattered source cells can protect many uninfected cells, and the self-amplification of interferon production might serve to calibrate response amplitude to strength of viral infection. Second, models of the tug-of-war between viral replication and the innate interferon response imply a pivotal role of interferon action on already infected cells in curbing viral spread, through effectively lowering virus replication rate. This finding is in line with the observation that several pathogenic viruses selectively abrogate interferon action on infected cells. Thus, interferons may delay viral spread in acute infections by acting as sentinels, warning uninfected cells of imminent danger, or as negative feedback regulators of virus replication in infected cells. The timing of the interferon response relative to the onset of viral replication is critical for its effectiveness in curbing viral spread.

Systems Level Analysis of the Yeast Osmo-Stat.

Adaptation is an important property of living organisms enabling them to cope with environmental stress and maintaining homeostasis. Adaptation is mediated by signaling pathways responding to different stimuli. Those signaling pathways might communicate in order to orchestrate the cellular response to multiple simultaneous stimuli, a phenomenon called crosstalk. Here, we investigate possible mechanisms of crosstalk between the High Osmolarity Glycerol (HOG) and the Cell Wall Integrity (CWI) pathways in yeast, which mediate adaptation to hyper- and hypo-osmotic challenges, respectively. We combine ensemble modeling with experimental investigations to test in quantitative terms different hypotheses about the crosstalk of the HOG and the CWI pathways. Our analyses indicate that for the conditions studied i) the CWI pathway activation employs an adaptive mechanism with a variable volume-dependent threshold, in contrast to the HOG pathway, whose activation relies on a fixed volume-dependent threshold, ii) there is no or little direct crosstalk between the HOG and CWI pathways, and iii) its mainly the HOG alone mediating adaptation of cellular osmotic pressure for both hyper- as well as hypo-osmotic stress. Thus, by iteratively combining mathematical modeling with experimentation we achieved a better understanding of regulatory mechanisms of yeast osmo-homeostasis and formulated new hypotheses about osmo-sensing.

To senesce or not to senesce: how primary human fibroblasts decide their cell fate after DNA damage.

Excessive DNA damage can induce an irreversible cell cycle arrest, called senescence, which is generally perceived as an important tumour-suppressor mechanism. However, it is unclear how cells decide whether to senesce or not after DNA damage. By combining experimental data with a parameterized mathematical model we elucidate this cell fate decision at the G1-S transition. Our model provides a quantitative and conceptually new understanding of how human fibroblasts decide whether DNA damage is beyond repair and senesce. Model and data imply that the G1-S transition is regulated by a bistable hysteresis switch with respect to Cdk2 activity, which in turn is controlled by the Cdk2/p21 ratio rather than cyclin abundance. We experimentally confirm the resulting predictions that to induce senescence i) in healthy cells both high initial and elevated background DNA damage are necessary and sufficient, and ii) in already damaged cells much lower additional DNA damage is sufficient. Our study provides a mechanistic explanation of a) how noise in protein abundances allows cells to overcome the G1-S arrest even with substantial DNA damage, potentially leading to neoplasia, and b) how accumulating DNA damage with age increasingly sensitizes cells for senescence.

Mathematical modelling of arsenic transport, distribution and detoxification processes in yeast.

Arsenic has a dual role as causative and curative agent of human disease. Therefore, there is considerable interest in elucidating arsenic toxicity and detoxification mechanisms. By an ensemble modelling approach, we identified a best parsimonious mathematical model which recapitulates and predicts intracellular arsenic dynamics for different conditions and mutants, thereby providing novel insights into arsenic toxicity and detoxification mechanisms in yeast, which could partly be confirmed experimentally by dedicated experiments. Specifically, our analyses suggest that: (i) arsenic is mainly protein-bound during short-term (acute) exposure, whereas glutathione-conjugated arsenic dominates during long-term (chronic) exposure, (ii) arsenic is not stably retained, but can leave the vacuole via an export mechanism, and (iii) Fps1 is controlled by Hog1-dependent and Hog1-independent mechanisms during arsenite stress. Our results challenge glutathione depletion as a key mechanism for arsenic toxicity and instead suggest that (iv) increased glutathione biosynthesis protects the proteome against the damaging effects of arsenic and that (v) widespread protein inactivation contributes to the toxicity of this metalloid. Our work in yeast may prove useful to elucidate similar mechanisms in higher eukaryotes and have implications for the use of arsenic in medical therapy.