Stochastic expression of the interferon-β gene

Single-cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10× Genomics Chromium single-cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human monocyte-derived dendritic cells infected with ZIKV at the single-cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN-dependent and -independent genes (the antiviral module). We modeled the ZIKV-specific antiviral state at the protein level, leveraging experimentally derived protein interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per-cell basis with experimental protein interaction data.

IMPORTANCE

Zika virus (ZIKV) remains a public health threat given its potential for re-emergence and the detrimental fetal outcomes associated with infection during pregnancy. Understanding the dynamics between ZIKV and its host is critical to understanding ZIKV pathogenesis. Through ZIKV-inclusive single-cell RNA sequencing (scRNA-seq), we demonstrate on the single-cell level the dynamic interplay between ZIKV and the host: the transcriptional program that restricts viral infection and ZIKV-mediated inhibition of that response. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool for gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

The population context is a driver of the heterogeneous response of epithelial cells to interferons

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 neighboring cells) of epithelial cells affects their response to interferons. We identified that cells at the edge of cellular colonies are more responsive than cells embedded within colonies. We determined that this spatial heterogeneity in interferon response resulted from the polarized basolateral interferon receptor distribution, making cells located in the center of cellular colonies less responsive to ectopic interferon stimulation. This was conserved across cell lines and primary cells originating from epithelial tissues. Importantly, cells embedded within cellular colonies were not protected from viral infection by apical interferon treatment, demonstrating that the population context-driven heterogeneous response to interferon influences the outcome of viral infection. Our data highlights that the behavior of isolated cells does not directly translate to their behavior in a population, placing the population context as one important factor influencing heterogeneity during interferon response in epithelial cells.

Probing T-cell activation in nanoliter tumor co-cultures using membrane displacement trap arrays

Immune responses against cancer are inherently stochastic, with small numbers of individual T cells within a larger ensemble of lymphocytes initiating the molecular cascades that lead to tumor cytotoxicity. A potential source of this intra-tumor variability is the differential ability of immune cells to respond to tumor cells. Classical microwell co-cultures of T cells and tumor cells are inadequate for reliably culturing and analyzing low cell numbers needed to probe this variability, and have failed in recapitulating the heterogeneous small domains observed in tumors. Here we leverage a membrane displacement trap array technology that overcomes limitations of conventional microwell plates for immunodynamic studies. The microfluidic platform supports on-demand formation of dense nanowell cultures under continuous perfusion reflecting the tumor microenvironment, with real-time monitoring of T cell proliferation and activation within each nanowell. The system enables selective ejection of cells for profiling by fluorescence activated cell sorting, allowing observed on-chip variability in immune response to be correlated with off-chip quantification of T cell activation. The technology offers new potential for probing the molecular origins of T cell heterogeneity and identifying specific cell phenotypes responsible for initiating and propagating immune cascades within tumors. Insight Box Variability in T cell activation plays a critical role in the immune response against cancer. New tools are needed to unravel the mechanisms that drive successful anti-tumor immune response, and to support the development of novel immunotherapies utilizing rare T cell phenotypes that promote effective immune surveillance. To this end, we present a microfluidic cell culture platform capable of probing differential T cell activation in an array of nanoliter-scale wells coupled with off-chip cell analysis, enabling a high resolution view of variable immune response within tumor / T cell co-cultures containing cell ensembles orders of magnitude smaller than conventional well plate studies.

Single cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10x Genomics Chromium single cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human moDCs infected with ZIKV at the single cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN dependent and independent genes (antiviral module). We modeled the ZIKV specific antiviral state at the protein level leveraging experimentally derived protein-interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per cell basis with experimental protein interaction data. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool to gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Early cellular mechanisms of type I interferon-driven susceptibility to tuberculosis

Mycobacterium tuberculosis (Mtb) causes 1.6 million deaths annually. Active tuberculosis correlates with a neutrophil-driven type I interferon (IFN) signature, but the cellular mechanisms underlying tuberculosis pathogenesis remain poorly understood. We found that interstitial macrophages (IMs) and plasmacytoid dendritic cells (pDCs) are dominant producers of type I IFN during Mtb infection in mice and non-human primates, and pDCs localize near human Mtb granulomas. Depletion of pDCs reduces Mtb burdens, implicating pDCs in tuberculosis pathogenesis. During IFN-driven disease, we observe abundant DNA-containing neutrophil extracellular traps (NETs) described to activate pDCs. Cell-type-specific disruption of the type I IFN receptor suggests that IFNs act on IMs to inhibit Mtb control. Single-cell RNA sequencing (scRNA-seq) indicates that type I IFN-responsive cells are defective in their response to IFNγ, a cytokine critical for Mtb control. We propose that pDC-derived type I IFNs act on IMs to permit bacterial replication, driving further neutrophil recruitment and active tuberculosis disease.

Antiviral responses are shaped by heterogeneity in viral replication dynamics

Antiviral signalling, which can be activated in host cells upon virus infection, restricts virus replication and communicates infection status to neighbouring cells. The antiviral response is heterogeneous, both quantitatively (efficiency of response activation) and qualitatively (transcribed antiviral gene set). To investigate the basis of this heterogeneity, we combined Virus Infection Real-time IMaging (VIRIM), a live-cell single-molecule imaging method, with real-time readouts of the dsRNA sensing pathway to analyse the response of human cells to encephalomyocarditis virus (EMCV) infection. We find that cell-to-cell heterogeneity in viral replication rates early in infection affect the efficiency of antiviral response activation, with lower replication rates leading to more antiviral response activation. Furthermore, we show that qualitatively distinct antiviral responses can be linked to the strength of the antiviral signalling pathway. Our analyses identify variation in early viral replication rates as an important parameter contributing to heterogeneity in antiviral response activation.

Antagonism between viral infection and innate immunity at the single-cell level

When infected with a virus, cells may secrete interferons (IFNs) that prompt nearby cells to prepare for upcoming infection. Reciprocally, viral proteins often interfere with IFN synthesis and IFN-induced signaling. We modeled the crosstalk between the propagating virus and the innate immune response using an agent-based stochastic approach. By analyzing immunofluorescence microscopy images we observed that the mutual antagonism between the respiratory syncytial virus (RSV) and infected A549 cells leads to dichotomous responses at the single-cell level and complex spatial patterns of cell signaling states. Our analysis indicates that RSV blocks innate responses at three levels: by inhibition of IRF3 activation, inhibition of IFN synthesis, and inhibition of STAT1/2 activation. In turn, proteins coded by IFN-stimulated (STAT1/2-activated) genes inhibit the synthesis of viral RNA and viral proteins. The striking consequence of these inhibitions is a lack of coincidence of viral proteins and IFN expression within single cells. The model enables investigation of the impact of immunostimulatory defective viral particles and signaling network perturbations that could potentially facilitate containment or clearance of the viral infection.

Comparison of antiviral responses in two bat species reveals conserved and divergent innate immune pathways

Bats host a range of disease-causing viruses without displaying clinical symptoms. The mechanisms behind this are a continuous source of interest. Here, we studied the antiviral response in the Egyptian fruit bat and Kuhl's pipistrelle, representing two subordinal clades. We profiled the antiviral response in fibroblasts using RNA sequencing and compared bat with primate and rodent responses. Both bats upregulate similar genes; however, a subset of these genes is transcriptionally divergent between them. These divergent genes also evolve rapidly in sequence, have specific promoter architectures, and are associated with programs underlying tolerance and resistance. Finally, we characterized antiviral genes that expanded in bats, with duplicates diverging in sequence and expression. Our study reveals a largely conserved antiviral program across bats and points to a set of genes that rapidly evolve through multiple mechanisms. These can contribute to bat adaptation to viral infection and provide directions to understanding the mechanisms behind it.

Control of enhancer and promoter activation in the type I interferon response by the histone demethylase Kdm4d/JMJD2d

Introduction

Transcriptional activation depends on the interplay of chromatin modifiers to establish a permissive epigenetic landscape. While histone 3 lysine 9 (H3K9) methylation has long been associated with gene repression, there is limited evidence to support a role for H3K9 demethylases in gene activation.

Methods

We leveraged knockdown and overexpression of JMJD2d / Kdm4d in mouse embryonic fibroblasts, coupled with extensive epigenomic analysesm to decipher the role of histone 3 lysine 9 demethylases in the innate immune response.

Results

Here we describe the H3K9 demethylase Kdm4d/JMJD2d as a positive regulator of type I interferon responses. In mouse embryonic fibroblasts (MEFs), depletion of JMJD2d attenuates the transcriptional response, conferring increased viral susceptibility, while overexpression of the demethylase results in more robust IFN activation. We find that the underlying mechanism of JMJD2d in type I interferon responses consists of an effect both on the transcription of enhancer RNAs (eRNAs) and on dynamic H3K9me2 at associated promoters. In support of these findings, we establish that JMJD2d is associated with enhancer regions throughout the genome prior to stimulation but is redistributed to inducible promoters in conjunction with transcriptional activation.

Discussion

Taken together, our data reveal JMJD2d as a chromatin modifier that connects enhancer transcription with promoter demethylation to modulate transcriptional responses.

Transiently heritable fates and quorum sensing drive early IFN-I response dynamics

Type I interferon (IFN-I)-mediated antiviral responses are central to host defense against viral infections. Crucial is the tight and well-orchestrated control of cellular decision-making leading to the production of IFN-Is. Innovative single-cell approaches revealed that the initiation of IFN-I production is limited to only fractions of 1-3% of the total population, both found in vitro, in vivo, and across cell types, which were thought to be stochastically regulated. To challenge this dogma, we addressed the influence of various stochastic and deterministic host-intrinsic factors on dictating early IFN-I responses, using a murine fibroblast reporter model. Epigenetic drugs influenced the percentage of responding cells. Next, with the classical Luria-Delbrück fluctuation test, we provided evidence for transient heritability driving responder fates, which was verified with mathematical modeling. Finally, while studying varying cell densities, we substantiated an important role for cell density in dictating responsiveness, similar to the phenomenon of quorum sensing. Together, this systems immunology approach opens up new avenues to progress the fundamental understanding on cellular decision-making during early IFN-I responses, which can be translated to other (immune) signaling systems.

Caspase-Mediated Regulation and Cellular Heterogeneity of the cGAS/STING Pathway in Kaposi's Sarcoma-Associated Herpesvirus Infection

As a result of the ongoing virus-host arms race, viruses have evolved numerous immune subversion strategies, many of which are aimed at suppressing the production of type I interferons (IFNs). Apoptotic caspases have recently emerged as important regulators of type I IFN signaling both in noninfectious contexts and during viral infection. Despite being widely considered antiviral factors since they can trigger cell death, several apoptotic caspases promote viral replication by suppressing innate immune response. Indeed, we previously discovered that the AIDS-associated oncogenic gammaherpesvirus Kaposi's sarcoma-associated herpesvirus (KSHV) exploits caspase activity to suppress the antiviral type I IFN response and promote viral replication. However, the mechanism of this novel viral immune evasion strategy is poorly understood, particularly with regard to how caspases antagonize IFN signaling during KSHV infection. Here, we show that caspase activity inhibits the DNA sensor cGAS during KSHV lytic replication to block type I IFN induction. Furthermore, we used single-cell RNA sequencing to reveal that the potent antiviral state conferred by caspase inhibition is mediated by an exceptionally small percentage of IFN-β-producing cells, thus uncovering further complexity of IFN regulation during viral infection. Collectively, these results provide insight into multiple levels of cellular type I IFN regulation that viruses co-opt for immune evasion. Unraveling these mechanisms can inform targeted therapeutic strategies for viral infections and reveal cellular mechanisms of regulating interferon signaling in the context of cancer and chronic inflammatory diseases. IMPORTANCE Type I interferons are key factors that dictate the outcome of infectious and inflammatory diseases. Thus, intricate cellular regulatory mechanisms are in place to control IFN responses. While viruses encode their own immune-regulatory proteins, they can also usurp existing cellular interferon regulatory functions. We found that caspase activity during lytic infection with the AIDS-associated oncogenic gammaherpesvirus Kaposi's sarcoma-associated herpesvirus inhibits the DNA sensor cGAS to block the antiviral type I IFN response. Moreover, single-cell RNA sequencing analyses unexpectedly revealed that an exceptionally small subset of infected cells (<5%) produce IFN, yet this is sufficient to confer a potent antiviral state. These findings reveal new aspects of type I IFN regulation and highlight caspases as a druggable target to modulate cGAS activity.

Revising immune cell coordination: Origins and importance of single-cell variation

Moving from the optimalization of single-cell technologies to the interpretation of the multi-complex single-cell data, the field of immunoengineering is granted with numerous important insights into the coordination of immune cell activation and how to modulate it for therapeutic purposes. However, insights come with additional follow-up questions that challenge our perception on how immune responses are generated and fine-tuned to fight a wide array of pathogens in ever-changing and often unpredictable microenvironments. Are immune responses really either being tightly regulated by molecular determinants, or highly flexible attributed to stochasticity? What exactly makes up the basic rules by which single cells cooperate to establish tissue-level immunity? Taking the type I IFN system and its newest insights as a main example throughout this review, we revise the basic concepts of (single) immune cell coordination, redefine the concepts of noise, stochasticity and determinism, and highlight the importance of single-cell variation in immunology and beyond.

Library-based analysis reveals segment and length dependent characteristics of defective influenza genomes

Found in a diverse set of viral populations, defective interfering particles are parasitic variants that are unable to replicate on their own yet rise to relatively high frequencies. Their presence is associated with a loss of population fitness, both through the depletion of key cellular resources and the stimulation of innate immunity. For influenza A virus, these particles contain large internal deletions in the genomic segments which encode components of the heterotrimeric polymerase. Using a library-based approach, we comprehensively profile the growth and replication of defective influenza species, demonstrating that they possess an advantage during genome replication, and that exclusion during population expansion reshapes population composition in a manner consistent with their final, observed, distribution in natural populations. We find that an innate immune response is not linked to the size of a deletion; however, replication of defective segments can enhance their immunostimulatory properties. Overall, our results address several key questions in defective influenza A virus biology, and the methods we have developed to answer those questions may be broadly applied to other defective viruses.

Single-cell secretion analysis reveals a dual role for IL-10 in restraining and resolving the TLR4-induced inflammatory response

Following Toll-like receptor 4 (TLR4) stimulation of macrophages, negative feedback mediated by the anti-inflammatory cytokine interleukin-10 (IL-10) limits the inflammatory response. However, extensive cell-to-cell variability in TLR4-stimulated cytokine secretion raises questions about how negative feedback is robustly implemented. To explore this, we characterize the TLR4-stimulated secretion program in primary murine macrophages using a single-cell microwell assay that enables evaluation of functional autocrine IL-10 signaling. High-dimensional analysis of single-cell data reveals three tiers of TLR4-induced proinflammatory activation based on levels of cytokine secretion. Surprisingly, while IL-10 inhibits TLR4-induced activation in the highest tier, it also contributes to the TLR4-induced activation threshold by regulating which cells transition from non-secreting to secreting states. This role for IL-10 in restraining TLR4 inflammatory activation is largely mediated by intermediate interferon (IFN)-β signaling, while TNF likely mediates response resolution by IL-10. Thus, cell-to-cell variability in cytokine regulatory motifs provides a means to tailor the TLR4-induced inflammatory response.

Virus Infection Variability by Single-Cell Profiling

Cell-to-cell variability of infection has long been known, yet it has remained one of the least understood phenomena in infection research. It impacts on disease onset and development, yet only recently underlying mechanisms have been studied in clonal cell cultures by single-virion immunofluorescence microscopy and flow cytometry. In this review, we showcase how single-cell RNA sequencing (scRNA-seq), single-molecule RNA-fluorescence in situ hybridization (FISH), and copper(I)-catalyzed azide-alkyne cycloaddition (click) with alkynyl-tagged viral genomes dissect infection variability in human and mouse cells. We show how the combined use of scRNA-FISH and click-chemistry reveals highly variable onsets of adenoviral gene expression, and how single live cell plaques reveal lytic and nonlytic adenovirus transmissions. The review highlights how scRNA-seq profiling and scRNA-FISH of coxsackie, influenza, dengue, zika, and herpes simplex virus infections uncover transcriptional variability, and how the host interferon response tunes influenza and sendai virus infections. We introduce the concept of "cell state" in infection variability, and conclude with advances by single-cell simultaneous measurements of chromatin accessibility and mRNA counts at high-throughput. Such technology will further dissect the sequence of events in virus infection and pathology, and better characterize the genetic and genomic stability of viruses, cell autonomous innate immune responses, and mechanisms of tissue injury.

Organoid modeling of Zika and herpes simplex virus 1 infections reveals virus-specific responses leading to microcephaly

Viral infection in early pregnancy is a major cause of microcephaly. However, how distinct viruses impair human brain development remains poorly understood. Here we use human brain organoids to study the mechanisms underlying microcephaly caused by Zika virus (ZIKV) and herpes simplex virus (HSV-1). We find that both viruses efficiently replicate in brain organoids and attenuate their growth by causing cell death. However, transcriptional profiling reveals that ZIKV and HSV-1 elicit distinct cellular responses and that HSV-1 uniquely impairs neuroepithelial identity. Furthermore, we demonstrate that, although both viruses fail to potently induce the type I interferon system, the organoid defects caused by their infection can be rescued by distinct type I interferons. These phenotypes are not seen in 2D cultures, highlighting the superiority of brain organoids in modeling viral infections. These results uncover virus-specific mechanisms and complex cellular immune defenses associated with virus-induced microcephaly.

Decoding the dynamics of multilayered stochastic antiviral IFN-I responses

Type I Interferon (IFN-I) responses were first recognized for their role in antiviral immunity, but it is now widely appreciated that IFN-Is have many immunomodulatory functions, influencing antitumor responses, autoimmune manifestations, and antimicrobial defenses. Given these pivotal roles, it may be surprising that multilayered stochastic events create highly heterogeneous, but tightly regulated, all-or-nothing cellular decisions. Recently, mathematical models have provided crucial insights into the stochastic nature of antiviral IFN-I responses, which we critically evaluate in this review. In this context, we emphasize the need for innovative single-cell technologies combined with mathematical models to further reveal, understand, and predict the complexity of the IFN-I system in physiological and pathological conditions that may be relevant to a plethora of diseases.

Agent-based modeling reveals benefits of heterogeneous and stochastic cell populations during cGAS-mediated IFNβ production

MOTIVATION

The cGAS pathway is a component of the innate immune system responsible for the detection of pathogenic DNA and upregulation of interferon beta (IFNβ). Experimental evidence shows that IFNβ signaling occurs in highly heterogeneous cells and is stochastic in nature; however, the benefits of these attributes remain unclear. To investigate how stochasticity and heterogeneity affect IFNβ production, an agent-based model is developed to simulate both DNA transfection and viral infection.

RESULTS

We show that heterogeneity can enhance IFNβ responses during infection. Furthermore, by varying the degree of IFNβ stochasticity, we find that only a percentage of cells (20-30%) need to respond during infection. Going beyond this range provides no additional protection against cell death or reduction of viral load. Overall, these simulations suggest that heterogeneity and stochasticity are important for moderating immune potency while minimizing cell death during infection.

AVAILABILITY AND IMPLEMENTATION

Model repository is available at: https://github.com/ImmuSystems-Lab/AgentBasedModel-cGASPathway.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Causes and Consequences of Spatial Organization of the Genome in Regulation of Gene Expression

Regulation of gene expression in time, space and quantity is orchestrated by the functional interplay of cis-acting elements and trans-acting factors. Our current view postulates that transcription factors recognize enhancer DNA and read the transcriptional regulatory code by cooperative DNA binding to specific DNA motifs, thus instructing the recruitment of transcriptional regulatory complexes forming a plethora of higher-ordered multi-protein-DNA and protein-protein complexes. Here, we reviewed the formation of multi-dimensional chromatin assemblies implicated in gene expression with emphasis on the regulatory role of enhancer hubs as coordinators of stochastic gene expression. Enhancer hubs contain many interacting regulatory elements and represent a remarkably dynamic and heterogeneous network of multivalent interactions. A functional consequence of such complex interaction networks could be that individual enhancers function synergistically to ensure coordination, tight control and robustness in regulation of expression of spatially connected genes. In this review, we discuss fundamental paradigms of such inter- and intra- chromosomal associations both in the context of immune-related genes and beyond.

Phenotypical Diversification of Early IFNα-Producing Human Plasmacytoid Dendritic Cells Using Droplet-Based Microfluidics

Plasmacytoid dendritic cells (pDCs) are a rare type of highly versatile immune cells that besides their specialized function of massive type I interferon (IFN-I) production are able to exert cytotoxic effector functions. However, diversification upon toll like receptor (TLR)-induced activation leads to highly heterogeneous responses that have not been fully characterized yet. Using droplet-based microfluidics, we showed that upon TLR7/8 and TLR9-induced single-cell activation only 1-3% secretes IFNα, and only small fractions upregulate cytotoxicity markers. Interestingly, this 1-3% of early IFN-producing pDCs, also known as first responders, express high levels of programmed death-ligand 1 (PD-L1) and TNF-related apoptosis-inducing ligand (TRAIL), which makes these hybrid cells similar to earlier described IFN-I producing killer pDCs (IKpDCs). IFN-I priming increases the numbers of IFNα producing cells up to 40%, but does not significantly upregulate the cytotoxicity markers. Besides, these so-called second responders do not show a cytotoxic phenotype as potent as observed for the first responders. Overall, our results indicate that the first responders are the key drivers orchestrating population wide IFN-I responses and possess high cytotoxic potential.

Detecting single cell interferon-beta production using a fluorescent reporter telomerase-immortalized human fibroblast cell line

Recent data suggest that cells respond to infection by upregulating the antiviral cytokine interferon-beta (IFN-ß) in a fraction of infected cells. Approaches are thus needed to study these responses on a single-cell level rather than bulk population. Here, we describe a protocol to analyze the IFN-ß response of individual cells using flow cytometry and immunofluorescence microscopy. We show the heterogeneous IFN-ß response to inactivated Sendai virus and human cytomegalovirus, but this protocol can be adapted to other viruses.

For complete details on the use and execution of this protocol, please refer to Hare et al. (2020).

•

Single-cell assays are needed to measure IFN production in virus-infected cells

•

Immortalized THF cells with IFN-b reporter are used to measure IFN production

•

Individual cells are analyzed via flow cytometry and fluorescence microscopy

•

Simple single-cell assays can uncover data obscured in bulk population measurements

Mutational pressure by host APOBEC3s more strongly affects genes expressed early in the lytic phase of herpes simplex virus-1 (HSV-1) and human polyomavirus (HPyV) infection

Herpes-Simplex Virus 1 (HSV-1) infects most humans when they are young, sometimes with fatal consequences. Gene expression occurs in a temporal order upon lytic HSV-1 infection: immediate early (IE) genes are expressed, then early (E) genes, followed by late (L) genes. During this infection cycle, the HSV-1 genome has the potential for exposure to APOBEC3 (A3) proteins, a family of cytidine deaminases that cause C>U mutations on single-stranded DNA (ssDNA), often resulting in a C>T transition. We developed a computational model for the mutational pressure of A3 on the lytic cycle of HSV-1 to determine which viral kinetic gene class is most vulnerable to A3 mutations. Using in silico stochastic methods, we simulated the infectious cycle under varying intensities of A3 mutational pressure. We found that the IE and E genes are more vulnerable to A3 than L genes. We validated this model by analyzing the A3 evolutionary footprints in 25 HSV-1 isolates. We find that IE and E genes have evolved to underrepresent A3 hotspot motifs more so than L genes, consistent with greater selection pressure on IE and E genes. We extend this model to two-step infections, such as those of polyomavirus, and find that the same pattern holds for over 25 human Polyomavirus (HPyVs) genomes. Genes expressed earlier during infection are more vulnerable to mutations than those expressed later.

Cell-to-cell variation in gene expression and the aging process

There is tremendous variation in biological traits, and much of it is not accounted for by variation in DNA sequence, including human diseases and lifespan. Emerging evidence points to differences in the execution of the genetic program as a key source of variation, be it stochastic variation or programmed variation. Here we discuss variation in gene expression as an intrinsic property and how it could contribute to variation in traits, including the rate of aging. The review is divided into sections describing the historical context and evidence to date for nongenetic variation, the different approaches that may be used to detect nongenetic variation, and recent findings showing that the amount of variation in gene expression can be both genetically programmed and epigenetically controlled. Finally, we present evidence that changes in cell-to-cell variation in gene expression emerge as part of the aging process and may be linked to disease vulnerability as a function of age. These emerging concepts are likely to be important across the spectrum of biomedical research and may well underpin what we understand as biological aging.

Importance of Type I and III Interferons at Respiratory and Intestinal Barrier Surfaces

Interferons (IFNs) constitute the first line of defense against microbial infections particularly against viruses. They provide antiviral properties to cells by inducing the expression of hundreds of genes known as interferon-stimulated genes (ISGs). The two most important IFNs that can be produced by virtually all cells in the body during intrinsic innate immune response belong to two distinct families: the type I and type III IFNs. The type I IFN receptor is ubiquitously expressed whereas the type III IFN receptor's expression is limited to epithelial cells and a subset of immune cells. While originally considered to be redundant, type III IFNs have now been shown to play a unique role in protecting mucosal surfaces against pathogen challenges. The mucosal specific functions of type III IFN do not solely rely on the restricted epithelial expression of its receptor but also on the distinct means by which type III IFN mediates its anti-pathogen functions compared to the type I IFN. In this review we first provide a general overview on IFNs and present the similarities and differences in the signal transduction pathways leading to the expression of either type I or type III IFNs. By highlighting the current state-of-knowledge of the two archetypical mucosal surfaces (e.g. the respiratory and intestinal epitheliums), we present the differences in the signaling cascades used by type I and type III IFNs to uniquely induce the expression of ISGs. We then discuss in detail the role of each IFN in controlling pathogen infections in intestinal and respiratory epithelial cells. Finally, we provide our perspective on novel concepts in the field of IFN (stochasticity, response heterogeneity, cellular polarization/differentiation and tissue microenvironment) that we believe have implications in driving the differences between type I and III IFNs and could explain the preferences for type III IFNs at mucosal surfaces.

Virus-Intrinsic Differences and Heterogeneous IRF3 Activation Influence IFN-Independent Antiviral Protection

Type 1 interferon (IFN) plays a critical role in early antiviral defense and priming of adaptive immunity by signaling upregulation of host antiviral IFN-stimulated genes (ISGs). Certain stimuli trigger strong activation of IFN regulatory factor 3 (IRF3) and direct upregulation of ISGs in addition to IFN. It remains unclear why some stimuli are stronger activators of IRF3 and how this leads to IFN-independent antiviral protection. We found that UV-inactivated human cytomegalovirus (HCMV) particles triggered an IFN-independent ISG signature that was absent in cells infected with UV-inactivated Sendai virus particles. HCMV particles triggered mostly uniform activation of IRF3 and low-level IFN-β production within the population while SeV particles triggered a small fraction of cells producing abundant IFN-β. These findings suggest that population-level activation of IRF3 and antiviral protection emerges from a diversity of responses occurring simultaneously in single cells. Moreover, this occurs in the absence of virus replication.

•

The antiviral response to virus particles requires low levels of interferon

•

Cells respond differently to HCMV or SeV particles

•

Heterogeneous IRF3 activation influences the response to virus

Poly I:C stimulation in-vitro as a marker for an antiviral response in different cell types generated from Buffalo (Bubalus bubalis)

The innate immune system is activated upon virus invasion of a host cell by recognizing viral component, such as dsRNA through specific receptors, resulting in the production of type- I IFNs, which confer an antiviral state within the invaded as well as surrounding cells. In the present study, fibroblast, monocyte and macrophage cells derived from water Buffalo (Bubalus bubalis) were exposed to a synthetic dsRNA analogue, poly I:C to mimic viral invasion in each cell type. Recognition of poly I:C through cytosolic helicase receptors RIG-I and MDA5 molecule lead to the activation of the RLR pathway, subsequently activating the MAVS-IRF3/7 cascade and the production of antiviral effector molecule like IFNβ and ISGs. Within the different cell types, we identified variability in RLR receptor and IFNβ expression after poly I:C administration. Fibroblasts responded quickly and strongly with IFNβ production, followed by macrophages and monocytes. Despite absolute expression variability among different cell types the expression trend of RLRs pathway genes were similar. Length of poly I:C molecule also influence IFNβ expression in response of RLR pathway. Short (LMW) poly I:C induce stronger IFN-β expression in myeloid (macrophage and monocyte) cells. In contrast long (HMW) poly I:C preferably elicit higher IFNβ expression in non-myeloid (fibroblast) cell. Therefore, MDA5 and RIG-1 plays an indispensable role in eliciting antiviral response in non- immune (fibroblast) host cell. Thus, stimulation of RLR pathway with suitable and potentially cell-type specific agonist molecules successfully elicit antiviral state in the host animal, with fibroblasts conferring a stronger antiviral state compared with the monocytes and macrophages.

Single-Cell Expression Variability Implies Cell Function

As single-cell RNA sequencing (scRNA-seq) data becomes widely available, cell-to-cell variability in gene expression, or single-cell expression variability (scEV), has been increasingly appreciated. However, it remains unclear whether this variability is functionally important and, if so, what are its implications for multi-cellular organisms. Here, we analyzed multiple scRNA-seq data sets from lymphoblastoid cell lines (LCLs), lung airway epithelial cells (LAECs), and dermal fibroblasts (DFs) and, for each cell type, selected a group of homogenous cells with highly similar expression profiles. We estimated the scEV levels for genes after correcting the mean-variance dependency in that data and identified 465, 466, and 364 highly variable genes (HVGs) in LCLs, LAECs, and DFs, respectively. Functions of these HVGs were found to be enriched with those biological processes precisely relevant to the corresponding cell type’s function, from which the scRNA-seq data used to identify HVGs were generated—e.g., cytokine signaling pathways were enriched in HVGs identified in LCLs, collagen formation in LAECs, and keratinization in DFs. We repeated the same analysis with scRNA-seq data from induced pluripotent stem cells (iPSCs) and identified only 79 HVGs with no statistically significant enriched functions; the overall scEV in iPSCs was of negligible magnitude. Our results support the “variation is function” hypothesis, arguing that scEV is required for cell type-specific, higher-level system function. Thus, quantifying and characterizing scEV are of importance for our understating of normal and pathological cellular processes.

Order by chance: origins and benefits of stochasticity in immune cell fate control

To protect against diverse challenges, the immune system must continuously generate an arsenal of specialized cell types, each of which can mount a myriad of effector responses upon detection of potential threats. To do so, it must generate multiple differentiated cell populations with defined sizes and proportions, often from rare starting precursor cells. Here, we discuss the emerging view that inherently probabilistic mechanisms, involving rare, rate-limiting regulatory events in single cells, control fate decisions and population sizes and fractions during immune development and function. We first review growing evidence that key fate control points are gated by stochastic signaling and gene regulatory events that occur infrequently over decision-making timescales, such that initially homogeneous cells can adopt variable outcomes in response to uniform signals. We next discuss how such stochastic control can provide functional capabilities that are harder to achieve with deterministic control strategies, and may be central to robust immune system function.

Regulation of interferon stimulated gene expression levels at homeostasis

Interferon stimulated genes (ISGs), a collection of genes important in the early innate immune response, are upregulated in response to stimulation by extracellular type I interferons. The regulation of ISGs has been extensively studied in cells exposed to significant interferon stimulation, but less is known about ISG regulation in homeostatic regimes in which extracellular interferon levels are low. Using a collection of pre-existing, publicly available microarray datasets, we investigated ISG regulation at homeostasis in CD4, pulmonary epithelial, fibroblast and macrophage cells. We used a linear regression model to predict ISG expression levels from regulator expression levels. Our results suggest significant regulation of ISG expression at homeostasis, both through the ISGF3 molecule and through IRF7 and IRF8 associated pathways. We find that roughly 50% of ISGs have expression levels significantly correlated with ISGF3 expression levels at homeostasis, supporting previous results suggesting that homeostatic IFN levels have broad functional consequences. We find that ISG expression levels varied in their correlation with ISGF3, with epithelial and macrophage cells showing more correlation than CD4 and fibroblast cells. Our analysis provides a novel approach for decomposing and quantifying ISG regulation.

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.

Porcine Interferon Complex and Co-Evolution with Increasing Viral Pressure after Domestication

Consisting of nearly 60 functional genes, porcine interferon (IFN)-complex represents an evolutionary surge of IFN evolution in domestic ungulate species. To compare with humans and mice, each of these species contains about 20 IFN functional genes, which are better characterized using the conventional IFN-α/β subtypes as examples. Porcine IFN-complex thus represents an optimal model for studying IFN evolution that resulted from increasing viral pressure during domestication and industrialization. We hypothesize and justify that porcine IFN-complex may extend its functionality in antiviral and immunomodulatory activity due to its superior molecular diversity. Furthermore, these unconventional IFNs could even confer some functional and signaling novelty beyond that of the well-studied IFN-α/β subtypes. Investigations into porcine IFN-complex will further our understanding of IFN biology and promote IFN-based therapeutic designs to confront swine viral diseases.

HSV-1 single-cell analysis reveals the activation of anti-viral and developmental programs in distinct sub-populations

Viral infection is usually studied at the population level by averaging over millions of cells. However, infection at the single-cell level is highly heterogeneous, with most infected cells giving rise to no or few viral progeny while some cells produce thousands. Analysis of Herpes Simplex virus 1 (HSV-1) infection by population-averaged measurements has taught us a lot about the course of viral infection, but has also produced contradictory results, such as the concurrent activation and inhibition of type I interferon signaling during infection. Here, we combine live-cell imaging and single-cell RNA sequencing to characterize viral and host transcriptional heterogeneity during HSV-1 infection of primary human cells. We find extreme variability in the level of viral gene expression among individually infected cells and show that these cells cluster into transcriptionally distinct sub-populations. We find that anti-viral signaling is initiated in a rare group of abortively infected cells, while highly infected cells undergo cellular reprogramming to an embryonic-like transcriptional state. This reprogramming involves the recruitment of β-catenin to the host nucleus and viral replication compartments, and is required for late viral gene expression and progeny production. These findings uncover the transcriptional differences in cells with variable infection outcomes and shed new light on the manipulation of host pathways by HSV-1.

West Nile Virus-Inclusive Single-Cell RNA Sequencing Reveals Heterogeneity in the Type I Interferon Response within Single Cells

West Nile virus (WNV) is a neurotropic mosquito-borne flavivirus of global importance. Neuroinvasive WNV infection results in encephalitis and can lead to prolonged neurological impairment or death. Type I interferon (IFN-I) is crucial for promoting antiviral defenses through the induction of antiviral effectors, which function to restrict viral replication and spread. However, our understanding of the antiviral response to WNV infection is mostly derived from analysis of bulk cell populations. It is becoming increasingly apparent that substantial heterogeneity in cellular processes exists among individual cells, even within a seemingly homogenous cell population. Here, we present WNV-inclusive single-cell RNA sequencing (scRNA-seq), an approach to examine the transcriptional variation and viral RNA burden across single cells. We observed that only a few cells within the bulk population displayed robust transcription of IFN-β mRNA, and this did not appear to depend on viral RNA abundance within the same cell. Furthermore, we observed considerable transcriptional heterogeneity in the IFN-I response, with genes displaying high unimodal and bimodal expression patterns. Broadly, IFN-stimulated genes negatively correlated with viral RNA abundance, corresponding with a precipitous decline in expression in cells with high viral RNA levels. Altogether, we demonstrated the feasibility and utility of WNV-inclusive scRNA-seq as a high-throughput technique for single-cell transcriptomics and WNV RNA detection. This approach can be implemented in other models to provide insights into the cellular features of protective immunity and identify novel therapeutic targets.IMPORTANCE West Nile virus (WNV) is a clinically relevant pathogen responsible for recurrent epidemics of neuroinvasive disease. Type I interferon is essential for promoting an antiviral response against WNV infection; however, it is unclear how heterogeneity in the antiviral response at the single-cell level impacts viral control. Specifically, conventional approaches lack the ability to distinguish differences across cells with varying viral abundance. The significance of our research is to demonstrate a new technique for studying WNV infection at the single-cell level. We discovered extensive variation in antiviral gene expression and viral abundance across cells. This protocol can be applied to primary cells or in vivo models to better understand the underlying cellular heterogeneity following WNV infection for the development of targeted therapeutic strategies.

Type I interferon response impairs differentiation potential of pluripotent stem cells

Unlike all differentiated cells, pluripotent stem cells do not elicit a productive antiviral response when infected by a pathogen. This observation seems at odds with the importance of pluripotent stem cells given their absolute requirement for the development of life. Here we investigate why this antiviral response is not utilized in these unique cells. We find that the factors required to maintain pluripotency are incompatible with those involved in eliciting the canonical interferon-based response to virus infection.

Upon virus infection, pluripotent stem cells neither induce nor respond to canonical type I interferons (IFN-I). To better understand this biology, we characterized induced pluripotent stem cells (iPSCs) as well as their differentiated parental or rederived counterparts. We confirmed that only iPSCs failed to respond to viral RNA, IFN-I, or viral infection. This lack of response could be phenocopied in fibroblasts with the expression of a reprogramming factor which repressed the capacity to induce canonical antiviral pathways. To ascertain the consequences of restoring the antiviral response in the context of pluripotency, we engineered a system to engage these defenses in iPSCs. Inducible expression of a recombinant virus-activated transcription factor resulted in the successful reconstitution of antiviral defenses through the direct up-regulation of IFN-I–stimulated genes. Induction of the antiviral signature in iPSCs, even for a short duration, resulted in the dysregulation of genes associated with all three germ layers despite maintaining pluripotency markers. Trilineage differentiation of these same cells showed that engagement of the antiviral defenses compromised ectoderm and endoderm formation and dysregulated the development of mesodermal sublineages. In all, these data suggest that the temporal induction of the antiviral response primes iPSCs away from pluripotency and induces numerous aberrant gene products upon differentiation. Together these results suggest that the IFN-I system and pluripotency may be incompatible with each other and thus explain why stem cells do not utilize the canonical antiviral system.

Gene expression variability across cells and species shapes innate immunity

As the first line of defence against pathogens, cells mount an innate immune response, which varies widely from cell to cell. The response must be potent but carefully controlled to avoid self-damage. How these constraints have shaped the evolution of innate immunity remains poorly understood. Here we characterize the innate immune response's transcriptional divergence between species and variability in expression among cells. Using bulk and single-cell transcriptomics in fibroblasts and mononuclear phagocytes from different species, challenged with immune stimuli, we map the architecture of the innate immune response. Transcriptionally diverging genes, including those that encode cytokines and chemokines, vary across cells and have distinct promoter structures. Conversely, genes that are involved in the regulation of this response, such as those that encode transcription factors and kinases, are conserved between species and display low cell-to-cell variability in expression. We suggest that this expression pattern, which is observed across species and conditions, has evolved as a mechanism for fine-tuned regulation to achieve an effective but balanced response.

Investigating Functional Roles for Positive Feedback and Cellular Heterogeneity in the Type I Interferon Response to Viral Infection

Secretion of type I interferons (IFN) by infected cells mediates protection against many viruses, but prolonged or excessive type I IFN secretion can lead to immune pathology. A proper type I IFN response must therefore maintain a balance between protection and excessive IFN secretion. It has been widely noted that the type I IFN response is driven by positive feedback and is heterogeneous, with only a fraction of infected cells upregulating IFN expression even in clonal cell lines, but the functional roles of feedback and heterogeneity in balancing protection and excessive IFN secretion are not clear. To investigate the functional roles for feedback and heterogeneity, we constructed a mathematical model coupling IFN and viral dynamics that extends existing mathematical models by accounting for feedback and heterogeneity. We fit our model to five existing datasets, reflecting different experimental systems. Fitting across datasets allowed us to compare the IFN response across the systems and suggested different signatures of feedback and heterogeneity in the different systems. Further, through numerical experiments, we generated hypotheses of functional roles for IFN feedback and heterogeneity consistent with our mathematical model. We hypothesize an inherent tradeoff in the IFN response: a positive feedback loop prevents excessive IFN secretion, but also makes the IFN response vulnerable to viral antagonism. We hypothesize that cellular heterogeneity of the IFN response functions to protect the feedback loop from viral antagonism. Verification of our hypotheses will require further experimental studies. Our work provides a basis for analyzing the type I IFN response across systems.

Single-cell analysis reveals that stochasticity and paracrine signaling control interferon-alpha production by plasmacytoid dendritic cells

Type I interferon (IFN) is a key driver of immunity to infections and cancer. Plasmacytoid dendritic cells (pDCs) are uniquely equipped to produce large quantities of type I IFN but the mechanisms that control this process are poorly understood. Here we report on a droplet-based microfluidic platform to investigate type I IFN production in human pDCs at the single-cell level. We show that type I IFN but not TNFα production is limited to a small subpopulation of individually stimulated pDCs and controlled by stochastic gene regulation. Combining single-cell cytokine analysis with single-cell RNA-seq profiling reveals no evidence for a pre-existing subset of type I IFN-producing pDCs. By modulating the droplet microenvironment, we demonstrate that vigorous pDC population responses are driven by a type I IFN amplification loop. Our study highlights the significance of stochastic gene regulation and suggests strategies to dissect the characteristics of immune responses at the single-cell level.

Plasmacytoid dendritic cells (pDC) are a pivotal component of the immune system. Here, the authors utilize single-cell microfluidics to interrogate the human pDC compartment and reveal a subset of type I IFN secreting pDCs that is regulated by stochastic gene expression and amplified by microenvironmental cues.

The Transcription Factor ThPOK Orchestrates Stochastic Interchromosomal Interactions Required for IFNB1 Virus-Inducible Gene Expression

Virus infection induces stochastic activation of the interferon-β gene. Three previously identified Alu-like DNA elements called NRCs (NF-κB reception centers) function by capturing and delivering NF-κB to the IFNB1 enhancer via stochastic interchromosomal interactions. We show that the transcription factor ThPOK binds cooperatively with NF-κB to NRCs and mediates their physical proximity with the IFNB1 gene via its ability to oligomerize when bound to DNA. ThPOK knockdown significantly decreased the frequency of interchromosomal interactions, NF-κB DNA binding to the IFNB1 enhancer, and virus-induced IFNB1 gene activation. We also demonstrate that cooperative DNA binding between ThPOK and NF-κB on the same face of the double DNA helix is required for interchromosomal interactions and distinguishes NRCs from various other Alu elements bearing κB sites. These studies show how DNA binding cooperativity of stereospecifically aligned transcription factors provides the necessary ultrasensitivity for the all-or-none stochastic cell responses to virus infection.

Cell fate in antiviral response arises in the crosstalk of IRF, NF-κB and JAK/STAT pathways

The innate immune system processes pathogen-induced signals into cell fate decisions. How information is turned to decision remains unknown. By combining stochastic mathematical modelling and experimentation, we demonstrate that feedback interactions between the IRF3, NF-κB and STAT pathways lead to switch-like responses to a viral analogue, poly(I:C), in contrast to pulse-like responses to bacterial LPS. Poly(I:C) activates both IRF3 and NF-κB, a requirement for induction of IFNβ expression. Autocrine IFNβ initiates a JAK/STAT-mediated positive-feedback stabilising nuclear IRF3 and NF-κB in first responder cells. Paracrine IFNβ, in turn, sensitises second responder cells through a JAK/STAT-mediated positive feedforward pathway that upregulates the positive-feedback components: RIG-I, PKR and OAS1A. In these sensitised cells, the 'live-or-die' decision phase following poly(I:C) exposure is shorter-they rapidly produce antiviral responses and commit to apoptosis. The interlinked positive feedback and feedforward signalling is key for coordinating cell fate decisions in cellular populations restricting pathogen spread.

Single-cell analysis of early antiviral gene expression reveals a determinant of stochastic IFNB1 expression

RIG-I-like receptors (RLRs) are cytoplasmic sensors of viral RNA that trigger the signaling cascade that leads to type I interferon (IFN) production. Transcriptional induction of RLRs by IFN is believed to play the role of positive feedback to further amplify viral sensing. We found that RLRs and several other IFN-stimulated genes (ISGs) are induced early in viral infection independent of IFN. Expression of these early ISGs requires IRF3/IRF7 and is highly correlated amongst them. Simultaneous detection of mRNA of IFNB1, viral replicase, and ISGs revealed distinct populations of IFNB1 expressing and non-expressing cells which are highly correlated with the levels of early ISGs but are uncorrelated with IFN-dependent ISGs and viral gene expression. Individual expression of RLRs made IFNB1 expression more robust and earlier, suggesting a causal relation between levels of RLR and induction of IFN.

Integrated modeling and analysis of intracellular and intercellular mechanisms in shaping the interferon response to viral infection

The interferons (IFNs) responses to viral infection are heterogeneous, while the underlying mechanisms for variability among cells are still not clear. In this study, we developed a hybrid model to systematically identify the sources of IFN induction heterogeneity. The experiment-integrated simulation demonstrated that the viral dose/type, the diversity in transcriptional factors activation and the intercellular paracrine signaling could strikingly shape the heterogeneity of IFN expression. We further determined that the IFNβ and IFNλ1 induced diverse dynamics of IFN-stimulated genes (ISGs) production. Collectively, our findings revealed the intracellular and intercellular mechanisms contributing to cell-to-cell variation in IFN induction, and further demonstrated the significant effects of IFN heterogeneity on antagonizing viruses.

Replication defective viral genomes exploit a cellular pro-survival mechanism to establish paramyxovirus persistence

Replication defective viral genomes (DVGs) generated during virus replication are the primary triggers of antiviral immunity in many RNA virus infections. However, DVGs can also facilitate viral persistence. Why and how these two opposing functions of DVGs are achieved remain unknown. Here we report that during Sendai and respiratory syncytial virus infections DVGs selectively protect a subpopulation of cells from death, thereby promoting the establishment of persistent infections. We find that during Sendai virus infection this phenotype results from DVGs stimulating a mitochondrial antiviral-signaling (MAVS)-mediated TNF response that drives apoptosis of highly infected cells while extending the survival of cells enriched in DVGs. The pro-survival effect of TNF depends on the activity of the TNFR2/TRAF1 pathway that is regulated by MAVS signaling. These results identify TNF as a pivotal factor in determining cell fate during a viral infection and delineate a MAVS/TNFR2-mediated mechanism that drives the persistence of otherwise acute viruses.

Replication defective viral genomes (DVGs) can facilitate persistence of paramyxoviruses, but the underlying mechanisms are unclear. Using FISH, Xu et al. here analyze the cellular response to DVGs on a single cell level and show that a MAVS-mediated TNF response specifically extends survival of cells enriched in DVGs.

Quantitative profiling of innate immune activation by viral infection in single cells

Cells infected by viruses can exhibit diverse patterns of viral and cellular gene expression. The patterns arise in part from the stochastic or noisy reaction kinetics associated with the small number of genomes, enzymes, and other molecules that typically initiate virus replication and activate cellular anti-viral defenses. It is not known what features, if any, of the early viral or cellular gene expression correlate with later processes of viral replication or cell survival. Here we used two fluorescent reporters to visualize innate immune activation of human prostate cancer (PC3) cells against infection by vesicular stomatitis virus. The cells were engineered to express green-fluorescent protein under control of the promoter for IFIT2, an interferon-sensitive component of the anti-viral response, while red-fluorescent protein was expressed as a byproduct of virus infection. To isolate and quantitatively analyze single-cells, we used a unique microwell array device and open-source image processing software. Kinetic analysis of viral and cellular reporter profiles from hundreds of cells revealed novel relationships between gene expression and the outcome of infection. Specifically, the relative timing rather than the magnitude of the viral gene expression and innate immune activation correlated with the infection outcome. Earlier viral or anti-viral gene expression favored or hindered virus growth, respectively. Further, analysis of kinetic parameters estimated from these data suggests a trade-off between robust antiviral signaling and cell death, as indicated by a higher rate of detectable cell lysis in infected cells with a detectable immune response. In short, cells that activate an immune response lyse at a higher rate. More broadly, we demonstrate how the intrinsic heterogeneity of individual cell behaviors can be exploited to discover features of viral and host gene expression that correlate with single-cell outcomes, which will ultimately impact whether or not infections spread.

A new efficient approach to fit stochastic models on the basis of high-throughput experimental data using a model of IRF7 gene expression as case study

Background

Mathematical models are used to gain an integrative understanding of biochemical processes and networks. Commonly the models are based on deterministic ordinary differential equations. When molecular counts are low, stochastic formalisms like Monte Carlo simulations are more appropriate and well established. However, compared to the wealth of computational methods used to fit and analyze deterministic models, there is only little available to quantify the exactness of the fit of stochastic models compared to experimental data or to analyze different aspects of the modeling results.

Results

Here, we developed a method to fit stochastic simulations to experimental high-throughput data, meaning data that exhibits distributions. The method uses a comparison of the probability density functions that are computed based on Monte Carlo simulations and the experimental data. Multiple parameter values are iteratively evaluated using optimization routines. The method improves its performance by selecting parameters values after comparing the similitude between the deterministic stability of the system and the modes in the experimental data distribution. As a case study we fitted a model of the IRF7 gene expression circuit to time-course experimental data obtained by flow cytometry. IRF7 shows bimodal dynamics upon IFN stimulation. This dynamics occurs due to the switching between active and basal states of the IRF7 promoter. However, the exact molecular mechanisms responsible for the bimodality of IRF7 is not fully understood.

Conclusions

Our results allow us to conclude that the activation of the IRF7 promoter by the combination of IRF7 and ISGF3 is sufficient to explain the observed bimodal dynamics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0406-4) contains supplementary material, which is available to authorized users.

Protein synthesis inhibition and GADD34 control IFN-β heterogeneous expression in response to dsRNA

In innate immune responses, induction of type-I interferons (IFNs) prevents virus spreading while viral replication is delayed by protein synthesis inhibition. We asked how cells perform these apparently contradictory activities. Using single fibroblast monitoring by flow cytometry and mathematical modeling, we demonstrate that type-I IFN production is linked to cell's ability to enter dsRNA-activated PKR-dependent translational arrest and then overcome this inhibition by decreasing eIF2α phosphorylation through phosphatase 1c cofactor GADD34 (Ppp1r15a) expression. GADD34 expression, shown here to be dependent on the IRF3 transcription factor, is responsible for a biochemical cycle permitting pulse of IFN synthesis to occur in cells undergoing protein synthesis inhibition. Translation arrest is further demonstrated to be key for anti-viral response by acting synergistically with MAVS activation to amplify TBK1 signaling and IFN-β mRNA transcription, while GADD34-dependent protein synthesis recovery contributes to the heterogeneous expression of IFN observed in dsRNA-activated cells.

Gene Expression Correlates with the Number of Herpes Viral Genomes Initiating Infection in Single Cells

Viral gene expression varies significantly among genetically identical cells. The sources of these variations are not well understood and have been suggested to involve both deterministic host differences and stochastic viral host interactions. For herpesviruses, only a limited number of incoming viral genomes initiate expression and replication in each infected cell. To elucidate the effect of this limited number of productively infecting genomes on viral gene expression in single cells, we constructed a set of fluorescence-expressing genetically tagged herpes recombinants. The number of different barcodes originating from a single cell is a good representative of the number of incoming viral genomes replicating (NOIVGR) in that cell. We identified a positive correlation between the NOIVGR and viral gene expression, as measured by the fluorescent protein expressed from the viral genome. This correlation was identified in three distinct cell-types, although the average NOIVGR per cell differed among these cell-types. Among clonal single cells, high housekeeping gene expression levels are not supportive of high viral gene expression, suggesting specific host determinants effecting viral infection. We developed a model to predict NOIVGR from cellular parameters, which supports the notion that viral gene expression is tightly linked to the NOIVGR in single-cells. Our results support the hypothesis that the stochastic nature of viral infection and host cell determinants contribute together to the variability observed among infected cells.

Single cell variation is of major interest in understanding key biological processes, like cancer, development and host pathogen interaction. During viral infection, these cell to cell variations can change the outcome of the whole organism infection. We suggested that differences in the number of parental viral genomes that initiate the replication process alter the outcome of infection among single cells. In this work we present a method based on genetically barcoded herpesvirus recombinants to identify the number of viral genomes initiating replication in individual cells. Our results indicate that viral gene expression is tightly linked to the number of viral genomes replicating per cell. Remarkably, we found that high cellular gene expression was an indicator for a lower viral gene expression in a given cell. We suggest that variations among single cells result from preexisting differences among cells, as well as from random viral host interactions.

Impact of Viral Infections on Hematopoiesis: From Beneficial to Detrimental Effects on Bone Marrow Output

The ability of the bone marrow (BM) to generate copious amounts of blood cells required on a daily basis depends on a highly orchestrated process of proliferation and differentiation of hematopoietic stem and progenitor cells (HSPCs). This process can be rapidly adapted under stress conditions, such as infections, to meet the specific cellular needs of the immune response and the ensuing physiological changes. This requires a tight regulation in order to prevent either hematopoietic failure or transformation. Although adaptation to bacterial infections or systemic inflammation has been studied and reviewed in depth, specific alterations of hematopoiesis to viral infections have received less attention so far. Viruses constantly pose a significant health risk and demand an adequate, balanced response from our immune system, which also affects the BM. In fact, both the virus itself and the ensuing immune response can have a tremendous impact on the hematopoietic process. On one hand, this can be beneficial: it helps to boost the cellular response of the body to resolve the viral infection. But on the other hand, when the virus and the resulting antiviral response persist, the inflammatory feedback to the hematopoietic system will become chronic, which can be detrimental for a balanced BM output. Chronic viral infections frequently have clinical manifestations at the level of blood cell formation, and we summarize which viruses can lead to BM pathologies, like aplastic anemia, pancytopenia, hemophagocytic lymphohistiocytosis, lymphoproliferative disorders, and malignancies. Regarding the underlying mechanisms, we address specific effects of acute and chronic viral infections on blood cell production. As such, we distinguish four different levels in which this can occur: (1) direct viral infection of HSPCs, (2) viral recognition by HSPCs, (3) indirect effects on HSPCs by inflammatory mediators, and (4) the role of the BM microenvironment on hematopoiesis upon virus infection. In conclusion, this review provides a comprehensive overview on how viral infections can affect the formation of new blood cells, aiming to advance our understanding of the underlying cellular and molecular mechanisms to improve the treatment of BM failure in patients.

Multilayered regulations of RIG-I in the anti-viral signaling pathway

RIG-I is a cytosolic receptor recognizing virus-specific RNA structures and initiates an antiviral signaling that induces the production of interferons and proinflammatory cytokines. Because inappropriate RIG-I signaling affects either viral clearance or immune toxicity, multiple regulations of RIG-I have been investigated since its discovery as the viral RNA detector. In this review, we describe the recent progress in research on the regulation of RIG-I activity or abundance. Specifically, we focus on the mechanism that modulates RIG-I-dependent antiviral response through post-translational modifications of or protein-protein interactions with RIG-I.

Type I Interferon-Mediated Induction of Antiviral Genes and Proteins Fails to Protect Cells from the Cytopathic Effects of Sendai Virus Infection

Sendai virus (SeV), a murine paramyxovirus, has been used to study the induction of type I interferon (IFN) subtypes in robust quantities. Few studies have measured whether the IFN that SeV induces actually fulfills its intended purpose of interfering with virus-mediated effects in the cells in which it is produced. We determined the effects of IFN on SeV-mediated cytopathic effects (CPE) and the ability of IFN to protect against virus infection. SeV-induced biologically active IFN resulted in Jak/STAT activation and the production of a number of interferon-stimulated genes (ISGs). However, these responses did not inhibit SeV replication or CPE. This observation was not due to SeV effects on canonical IFN signaling. Furthermore, pretreating cells with type I IFN and establishing an antiviral state before infection did not mediate SeV effects. Therefore, the induction of canonical IFN signaling pathways and ISGs does not always confer protection against the IFN-inducing virus. Because type I IFNs are approved to treat various infections, our findings suggest that typical markers of IFN activity may not be indicative of a protective antiviral response and should not be used alone to determine whether an antiviral state against a particular virus is achieved.

End Sequence Analysis Toolkit (ESAT) expands the extractable information from single-cell RNA-seq data

RNA-seq protocols that focus on transcript termini are well suited for applications in which template quantity is limiting. Here we show that, when applied to end-sequencing data, analytical methods designed for global RNA-seq produce computational artifacts. To remedy this, we created the End Sequence Analysis Toolkit (ESAT). As a test, we first compared end-sequencing and bulk RNA-seq using RNA from dendritic cells stimulated with lipopolysaccharide (LPS). As predicted by the telescripting model for transcriptional bursts, ESAT detected an LPS-stimulated shift to shorter 3′-isoforms that was not evident by conventional computational methods. Then, droplet-based microfluidics was used to generate 1000 cDNA libraries, each from an individual pancreatic islet cell. ESAT identified nine distinct cell types, three distinct β-cell types, and a complex interplay between hormone secretion and vascularization. ESAT, then, offers a much-needed and generally applicable computational pipeline for either bulk or single-cell RNA end-sequencing.