Search

Isogenic cells respond in a heterogeneous manner to interferon. Using a micropatterning approach combined with high-content imaging and spatial analyses, we characterized how the population context (position of a cell with respect to the neighboring cells) of human intestinal epithelial cells affects single cell response to interferons. We identified that cells at the edge of a cellular colony are significantly more responsive than cells embedded within this colony. We determined that this spatial heterogeneity in IFN response was the result of the polarized basolateral distribution of the IFN receptors making cells located in the center of a cellular colony not responsive to ectopic IFN stimulation. We could demonstrate that this population context driven cell-to-cell variability influences the outcome of viral infection as cells embedded in a cellular colony are not protected by interferons and therefore more susceptible to infection. Our data highlights that the behavior of individual isolated cells does not directly translate to their behavior in a population, placing the population context as a key driver of cell-to-cell heterogeneity in IFN response.

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 mRNA levels as the molecular memory of the ISR that plays a central role in cell adaptation to acute and chronic stress.

Zika virus (ZIKV) infection can cause developmental and neurological defects and represents a threat for human health. 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, defined as the probability for one virus to initiate infection, 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 quasi-species pinpointed mutations whose increased frequencies temporally coincide with the acquisition of the adapted phenotype. We functionally validated a point-mutation in ZIKV envelope (E) protein recapitulating the adapted phenotype. Its positioning on the E structure suggests a putative function in protein refolding/stability. Altogether, our results uncovered ZIKV adaptations to the cell environment leading to an accelerated replication onset coupled with resistance to TLR3-induced antiviral response. Our work provides insights into viral escape mechanisms and interactions with host cell and can serve as a framework to study other viruses.

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.

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.

It is currently believed that type I and III interferons (IFNs) have redundant functions. However, the preferential distribution of type III IFN receptor on epithelial cells suggests functional differences at epithelial surfaces. Here, using human intestinal epithelial cells we could show that although both type I and type III IFNs confer an antiviral state to the cells, they do so with distinct kinetics. Type I IFN signaling is characterized by an acute strong induction of interferon stimulated genes (ISGs) and confers fast antiviral protection. On the contrary, the slow acting type III IFN mediated antiviral protection is characterized by a weaker induction of ISGs in a delayed manner compared to type I IFN. Moreover, while transcript profiling revealed that both IFNs induced a similar set of ISGs, their temporal expression strictly depended on the IFNs, thereby leading to unique antiviral environments. Using a combination of data-driven mathematical modeling and experimental validation, we addressed the molecular reason for this differential kinetic of ISG expression. We could demonstrate that these kinetic differences are intrinsic to each signaling pathway and not due to different expression levels of the corresponding IFN receptors. We report that type III IFN is specifically tailored to act in specific cell types not only due to the restriction of its receptor but also by providing target cells with a distinct antiviral environment compared to type I IFN. We propose that this specific environment is key at surfaces that are often challenged with the extracellular environment., Author summary The human intestinal tract plays two important roles in the body: first it is responsible for nutrient absorption and second it is the primary barrier which protects the human body from the outside environment. This complex tissue is constantly exposed to commensal bacteria and is often exposed to both bacterial and viral pathogens. To protect itself, the gut produces, among others, secreted agents called interferons which help to fight against pathogen attacks. There are several varieties (type I, II, and III) of interferons and our work aims at understanding how type I and III interferon act to protect human intestinal epithelial cells (hIECs) during viral infection. In this study, we confirmed that both interferons can protect hIECs against viral infection but with different kinetics. We determined that type I confer an antiviral state to hIECs faster than type III interferons. We uncovered that these differences were intrinsic to each pathway and not the result of differential abundance of the respective interferon receptors. The results of this study suggest that type III interferon may provide a different antiviral environment to the epithelium target cells which is likely critical for maintaining gut homeostasis. Our findings will also help us to design therapies to aid in controlling and eliminating viral infections of the gut.

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.

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.

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.

Mathematical modelling has been instrumental to understand kinetics of radiation-induced DNA damage repair and associated secondary cancer risk. The widely accepted two-lesion kinetic (TLK) model assumes two kinds of double strand breaks, simple and complex ones, with different repair rates. Recently, persistent DNA damage associated with telomeres was reported as a new kind of DNA damage. We therefore extended existing versions of the TLK model by new categories of DNA damage and re-evaluated those models using extensive data. We subjected different versions of the TLK model to a rigorous model discrimination approach. This enabled us to robustly select a best approximating parsimonious model that can both recapitulate and predict transient and persistent DNA damage after ionizing radiation. Models and data argue for i) nonlinear dose-damage relationships and ii) negligible saturation of repair kinetics even for high doses. Additionally, we show that simulated radiation-induced persistent telomere-associated DNA damage foci (TAF) can be used to predict excess relative risk (ERR) of developing secondary leukemia after fractionated radiotherapy. We suggest that TAF may serve as an additional measure to predict cancer risk after radiotherapy using high dose rates. This may improve predicting risk-dose dependency of ionizing radiation especially for long-term therapies.

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.