Arrested spread of vesicular stomatitis virus infections in vitro depends on interferon-mediated antiviral activity

Functional Analysis of Oligoadenylate Synthetase in the Emu (Dromaius novaehollandiae)

2'-5'-oligoadenylate synthetase (OAS) is one of the proteins that act as a defense mechanism against foreign RNA in cells. OAS has two functions: an antiviral effect against a wide range of virus species via the OAS/RNase L pathway with synthesized oligoadenylates and inhibition of viral replication specific to viruses of the genus Flavivirus, which is independent of enzymatic activity. Several birds have been reported to possess only one type of OAS family member, OASL, which has both enzymatic activity and inhibitory effects on flaviviral replication. However, the ostrich has two types of OASs, OAS1 and OASL, which show different functions-enzymatic and anti-flaviviral activities, respectively. In this study, emu OASs were cloned to investigate their sequence and function and elucidate the role of OASs in emus. The cloning results showed that emus had OAS1 and OASL, suggesting that emu OASs were more closely related to ostrich than to other birds. Functional investigations showed that emu OAS1 and OASL had enzymatic and anti-flaviviral activities, respectively, similar to those of the ostrich. Emus and ostriches are evolutionarily different from most birds and may be more closely related to mammalian OAS diversity.

Computational multigene interactions in virus growth and infection spread

Viruses persist in nature owing to their extreme genetic heterogeneity and large population sizes, which enable them to evade host immune defenses, escape antiviral drugs, and adapt to new hosts. The persistence of viruses is challenging to study because mutations affect multiple virus genes, interactions among genes in their impacts on virus growth are seldom known, and measures of viral fitness are yet to be standardized. To address these challenges, we employed a data-driven computational model of cell infection by a virus. The infection model accounted for the kinetics of viral gene expression, functional gene-gene interactions, genome replication, and allocation of host cellular resources to produce progeny of vesicular stomatitis virus, a prototype RNA virus. We used this model to computationally probe how interactions among genes carrying up to eleven deleterious mutations affect different measures of virus fitness: single-cycle growth yields and multicycle rates of infection spread. Individual mutations were implemented by perturbing biophysical parameters associated with individual gene functions of the wild-type model. Our analysis revealed synergistic epistasis among deleterious mutations in their effects on virus yield; so adverse effects of single deleterious mutations were amplified by interaction. For the same mutations, multicycle infection spread indicated weak or negligible epistasis, where single mutations act alone in their effects on infection spread. These results were robust to simulation in high- and low-host resource environments. Our work highlights how different types and magnitudes of epistasis can arise for genetically identical virus variants, depending on the fitness measure. More broadly, gene-gene interactions can differently affect how viruses grow and spread.

Patterns of virus growth across the diversity of life

Although viruses in their natural habitats add up to less than 10% of the biomass, they contribute more than 90% of the genome sequences [1]. These viral sequences or 'viromes' encode viruses that populate the Earth's oceans [2, 3] and terrestrial environments [4, 5], where their infections impact life across diverse ecological niches and scales [6, 7], including humans [8-10]. Most viruses have yet to be isolated and cultured [11-13], and surprisingly few efforts have explored what analysis of available data might reveal about their nature. Here, we compiled and analyzed seven decades of one-step growth and other data for viruses from six major families, including their infections of archaeal, bacterial and eukaryotic hosts [14-191]. We found that the use of host cell biomass for virus production was highest for archaea at 10%, followed by bacteria at 1% and eukarya at 0.01%, highlighting the degree to which viruses of archaea and bacteria exploit their host cells. For individual host cells, the yield of virus progeny spanned a relatively narrow range (10-1000 infectious particles per cell) compared with the million-fold difference in size between the smallest and largest cells. Furthermore, healthy and infected host cells were remarkably similar in the time they needed to multiply themselves or their virus progeny. Specifically, the doubling time of healthy cells and the delay time for virus release from infected cells were not only correlated (r = 0.71, p < 10-10, n = 101); they also spanned the same range from tens of minutes to about a week. These results have implications for better understanding the growth, spread and persistence of viruses in complex natural habitats that abound with diverse hosts, including humans and their associated microbes.

Pixel precise unsupervised detection of viral particle proliferation in cellular imaging data

Cellular and molecular imaging techniques and models have been developed to characterize single stages of viral proliferation after focal infection of cells in vitro. The fast and automatic classification of cell imaging data may prove helpful prior to any further comparison of representative experimental data to mathematical models of viral propagation in host cells. Here, we use computer generated images drawn from a reproduction of an imaging model from a previously published study of experimentally obtained cell imaging data representing progressive viral particle proliferation in host cell monolayers. Inspired by experimental time-based imaging data, here in this study viral particle increase in time is simulated by a one-by-one increase, across images, in black or gray single pixels representing dead or partially infected cells, and hypothetical remission by a one-by-one increase in white pixels coding for living cells in the original image model. The image simulations are submitted to unsupervised learning by a Self-Organizing Map (SOM) and the Quantization Error in the SOM output (SOM-QE) is used for automatic classification of the image simulations as a function of the represented extent of viral particle proliferation or cell recovery. Unsupervised classification by SOM-QE of 160 model images, each with more than three million pixels, is shown to provide a statistically reliable, pixel precise, and fast classification model that outperforms human computer-assisted image classification by RGB image mean computation. The automatic classification procedure proposed here provides a powerful approach to understand finely tuned mechanisms in the infection and proliferation of virus in cell lines in vitro or other cells.

•

Automatic classification of cell imaging data for in vitro viral propagation in host cells.

•

Economical in terms of computation times for rapid change/no change detection in large image data prior to human decision making.

•

Simulating viral proliferation/cell recovery by progressive and selective single pixel changes in contrast polarity and/or intensity.

Identification of Natural Molecular Determinants of Ross River Virus Type I Interferon Modulation

Ross River virus (RRV) belongs to the genus Alphavirus and is prevalent in Australia. RRV infection can cause arthritic symptoms in patients and may include rash, fever, arthralgia, and myalgia. Type I interferons (IFN) are the primary antiviral cytokines and trigger activation of the host innate immune system to suppress the replication of invading viruses. Alphaviruses are able to subvert the type I IFN system, but the mechanisms used are ill defined. In this study, seven RRV field strains were analyzed for induction of and sensitivity to type I IFN. The sensitivities of these strains to human IFN-β varied significantly and were highest for the RRV 2548 strain. Compared to prototype laboratory strain RRV-T48, RRV 2548 also induced higher type I IFN levels both in vitro and in vivo and caused milder disease. To identify the determinants involved in type I IFN modulation, the region encoding the nonstructural proteins (nsPs) of RRV 2548 was sequenced, and 42 amino acid differences from RRV-T48 were identified. Using fragment swapping and site-directed mutagenesis, we discovered that substitutions E402A and R522Q in nsP1 as well as Q619R in nsP2 were responsible for increased sensitivity of RRV 2548 to type I IFN. In contrast, substitutions A31T, N219T, S580L, and Q619R in nsP2 led to induction of higher levels of type I IFN. With exception of E402A, all these variations are common for naturally occurring RRV strains. However, they are different from all known determinants of type I IFN modulation reported previously in nsPs of alphaviruses.IMPORTANCE By identifying natural Ross River virus (RRV) amino acid determinants for type I interferon (IFN) modulation, this study gives further insight into the mechanism of type I IFN modulation by alphaviruses. Here, the crucial role of type I IFN in the early stages of RRV disease pathogenesis is further demonstrated. This study also provides a comparison of the roles of different parts of the RRV nonstructural region in type I IFN modulation, highlighting the importance of nonstructural protein 1 (nsP1) and nsP2 in this process. Three substitutions in nsP1 and nsP2 were found to be independently associated with enhanced type I IFN sensitivity, and four independent substitutions in nsP2 were important in elevated type I IFN induction. Such evidence has clear implications for RRV immunobiology, persistence, and pathology. The identification of viral proteins that modulate type I IFN may also have importance for the pathogenesis of other alphaviruses.

The Curious Case of the Nidovirus Exoribonuclease: Its Role in RNA Synthesis and Replication Fidelity

Among RNA viruses, the order Nidovirales stands out for including viruses with the largest RNA genomes currently known. Nidoviruses employ a complex RNA-synthesizing machinery comprising a variety of non-structural proteins (nsps). One of the postulated drivers of the expansion of nidovirus genomes is the presence of a proofreading 3′-to-5′ exoribonuclease (ExoN) belonging to the DEDDh family. ExoN may enhance the fidelity of RNA synthesis by correcting nucleotide incorporation errors made by the RNA-dependent RNA polymerase. Here, we review our current understanding of ExoN evolution, structure, and function. Most experimental data are derived from studies of the ExoN domain of coronaviruses (CoVs), which were triggered by the bioinformatics-based identification of ExoN in the genome of severe acute respiratory syndrome coronavirus (SARS-CoV) and its relatives in 2003. Although convincing data supporting the proofreading hypothesis have been obtained, from biochemical assays and studies with CoV mutants lacking ExoN functionality, the features of ExoN from most other nidovirus families remain to be characterized. Remarkably, viable ExoN knockout mutants were obtained only for two CoVs, mouse hepatitis virus (MHV) and SARS-CoV, whose RNA synthesis and replication kinetics were mildly affected by the lack of ExoN function. In several other CoV species, ExoN inactivation was not tolerated, and knockout mutants could not be rescued when launched using a reverse genetics system. This suggests that ExoN is also critical for primary viral RNA synthesis, a property that poorly matches the profile of an enzyme that would merely boost long-term replication fidelity. In CoVs, ExoN resides in a bifunctional replicase subunit (nsp14) whose C-terminal part has (N7-guanine)-methyltransferase activity. The crystal structure of SARS-CoV nsp14 has shed light on the interplay between these two domains, and on nsp14’s interactions with nsp10, a co-factor that strongly enhances ExoN activity in vitro assays. Further elucidation of the structure-function relationships of ExoN and its interactions with other (viral and/or host) members of the CoV replication machinery will be key to understanding the enzyme’s role in viral RNA synthesis and pathogenesis, and may contribute to the design of new approaches to combat emerging nidoviruses.

HCV and flaviviruses hijack cellular mechanisms for nuclear STAT2 degradation: Up-regulation of PDLIM2 suppresses the innate immune response

Host encounters with viruses lead to an innate immune response that must be rapid and broadly targeted but also tightly regulated to avoid the detrimental effects of unregulated interferon expression. Viral stimulation of host negative regulatory mechanisms is an alternate method of suppressing the host innate immune response. We examined three key mediators of the innate immune response: NF-KB, STAT1 and STAT2 during HCV infection in order to investigate the paradoxical induction of an innate immune response by HCV despite a multitude of mechanisms combating the host response. During infection, we find that all three are repressed only in HCV infected cells but not in uninfected bystander cells, both in vivo in chimeric mouse livers and in cultured Huh7.5 cells after IFNα treatment. We show here that HCV and Flaviviruses suppress the innate immune response by upregulation of PDLIM2, independent of the host interferon response. We show PDLIM2 is an E3 ubiquitin ligase that also acts to stimulate nuclear degradation of STAT2. Interferon dependent relocalization of STAT1/2 to the nucleus leads to PDLIM2 ubiquitination of STAT2 but not STAT1 and the proteasome-dependent degradation of STAT2, predominantly within the nucleus. CRISPR/Cas9 knockout of PDLIM2 results in increased levels of STAT2 following IFNα treatment, retention of STAT2 within the nucleus of HCV infected cells after IFNα stimulation, increased interferon response, and increased resistance to infection by several flaviviruses, indicating that PDLIM2 is a global regulator of the interferon response.

The response of cells to an invading pathogen must be swift and well controlled because of the detrimental effects of chronic inflammation. However, viruses often hijack host control mechanisms. HCV and flaviviruses are known to suppress the innate immune response in cells by a variety of mechanisms. This study clarifies and expands a specific cellular mechanism for global control of the antiviral response after the induction of interferon expression. It shows how several viruses hijack this control mechanism to suppress the innate interferon response.

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.

Spatial-Temporal Patterns of Viral Amplification and Interference Initiated by a Single Infected Cell

When viruses infect their host cells, they can make defective virus-like particles along with intact virus. Cells coinfected with virus and defective particles often exhibit interference with virus growth caused by the competition for resources by defective genomes. Recent reports of the coexistence and cotransmission of such defective interfering particles (DIPs) in vivo, across epidemiological length and time scales, suggest a role in viral pathogenesis, but it is not known how DIPs impact infection spread, even under controlled culture conditions. Using fluorescence microscopy, we quantified coinfections of vesicular stomatitis virus (VSV) expressing a fluorescent reporter protein and its DIPs on BHK-21 host cell monolayers. We found that viral gene expression was more delayed, infections spread more slowly, and patterns of spread became more “patchy” with higher DIP inputs to the initial cell. To examine how infection spread might depend on the behavior of the initial coinfected cell, we built a computational model, adapting a cellular automaton (CA) approach to incorporate kinetic data on virus growth for the first time. Specifically, changes in observed patterns of infection spread could be directly linked to previous high-throughput single-cell measures of virus-DIP coinfection. The CA model also provided testable hypotheses on the spatial-temporal distribution of the DIPs, which remain governed by their predator-prey interaction. More generally, this work offers a data-driven computational modeling approach for better understanding of how single infected cells impact the multiround spread of virus infections across cell populations.

IMPORTANCE Defective interfering particles (DIPs) compete with intact virus, depleting host cell resources that are essential for virus growth and infection spread. However, it is not known how such competition, strong or weak, ultimately affects the way in which infections spread and cause disease. In this study, we address this unmet need by developing an integrated experimental-computational approach, which sheds new light on how infections spread. We anticipate that our approach will also be useful in the development of DIPs as therapeutic agents to manage the spread of viral infections.

Rapid induction and persistence of paracrine-induced cellular antiviral states arrest viral infection spread in A549 cells

The virus/host interaction is a complex interplay between pro- and anti-viral factors that ultimately determines the spread or halt of virus infections in tissues. This interplay develops over multiple rounds of infection. The purpose of this study was to determine how cellular-level processes combine to impact the spatial spread of infection. We measured the kinetics of virus replication (VSV), antiviral paracrine signal upregulation and secretion, spatial spread of virus and paracrine antiviral signaling, and inhibition of virus production in antiviral-exposed A549 human lung epithelial cells. We found that initially infected cells released antiviral signals 4-to-7h following production of virus. However, the subsequent rapid dissemination of signal and fast induction of a robust and persistent antiviral state ultimately led to a suppression of infection spread. This work shows how cellular responses to infection and activation of antiviral responses can integrate to ultimately control infection spread across host cell populations.

High-Throughput Single-Cell Kinetics of Virus Infections in the Presence of Defective Interfering Particles

Defective interfering particles (DIPs) are virus mutants that lack essential genes for growth. In coinfections with helper virus, the diversion of viral proteins to the replication and packaging of DIP genomes can interfere with virus production. Mounting cases of DIPs and DIP-like genomes in clinical and natural isolates, as well as growing interest in DIP-based therapies, underscore a need to better elucidate how DIPs work. DIP activity is primarily measured by its inhibition of virus infection yield, an endpoint that masks the dynamic and potentially diverse individual cell behaviors. Using vesicular stomatitis virus (VSV) as a model, we coinfected BHK cells with VSV DIPs and recombinant helper virus carrying a gene encoding a red fluorescent protein (RFP) whose expression correlates with the timing and level of virus release. For single cells within a monolayer, 10 DIPs per cell suppressed the reporter expression in only 1.2% of the cells. In most cells, it slowed and reduced viral gene expression, manifested as a shift in mean latent time from 4 to 6 h and reduced virus yields by 10-fold. For single cells isolated in microwells, DIP effects were more pronounced, reducing virus yields by 100-fold and extending latent times to 12 h, including individual instances above 20 h. Together, these results suggest that direct or indirect cell-cell interactions prevent most coinfected cells from being completely suppressed by DIPs. Finally, a gamma distribution model captures well how the infection kinetics quantitatively depends on the DIP dose. Such models will be useful for advancing a predictive biology of DIP-associated virus growth and infection spread.

IMPORTANCE

During the last century, basic studies in virology have focused on developing a molecular mechanistic understanding of how infectious viruses reproduce in their living host cells. However, over the last 10 years, the advent of deep sequencing and other powerful technologies has revealed in natural and patient infections that viruses do not act alone. Instead, viruses are often accompanied by defective virus-like particles that carry large deletions in their genomes and fail to replicate on their own. Coinfections of viable and defective viruses behave in unpredictable ways, but they often interfere with normal virus growth, potentially enabling infections to evade host immune surveillance. In the current study, controlled levels of defective viruses are coinfected with viable viruses that have been engineered to express a fluorescent reporter protein during infection. Unique profiles of reporter expression acquired from thousands of coinfected cells reveal how interference acts at multiple stages of infection.

Genome rearrangement affects RNA virus adaptability on prostate cancer cells

Gene order is often highly conserved within taxonomic groups, such that organisms with rearranged genomes tend to be less fit than wild type gene orders, and suggesting natural selection favors genome architectures that maximize fitness. But it is unclear whether rearranged genomes hinder adaptability: capacity to evolutionarily improve in a new environment. Negative-sense non-segmented RNA viruses (order Mononegavirales) have specific genome architecture: 3' UTR - core protein genes - envelope protein genes - RNA-dependent RNA-polymerase gene - 5' UTR. To test how genome architecture affects RNA virus evolution, we examined vesicular stomatitis virus (VSV) variants with the nucleocapsid (N) gene moved sequentially downstream in the genome. Because RNA polymerase stuttering in VSV replication causes greater mRNA production in upstream genes, N gene translocation toward the 5' end leads to stepwise decreases in N transcription, viral replication and progeny production, and also impacts the activation of type 1 interferon mediated antiviral responses. We evolved VSV gene-order variants in two prostate cancer cell lines: LNCap cells deficient in innate immune response to viral infection, and PC-3 cells that mount an IFN stimulated anti-viral response to infection. We observed that gene order affects phenotypic adaptability (reproductive growth; viral suppression of immune function), especially on PC-3 cells that strongly select against virus infection. Overall, populations derived from the least-fit ancestor (most-altered N position architecture) adapted fastest, consistent with theory predicting populations with low initial fitness should improve faster in evolutionary time. Also, we observed correlated responses to selection, where viruses improved across both hosts, rather than suffer fitness trade-offs on unselected hosts. Whole genomics revealed multiple mutations in evolved variants, some of which were conserved across selective environments for a given gene order.

Phospholipase D facilitates efficient entry of influenza virus, allowing escape from innate immune inhibition

Lipid metabolism plays a fundamental role during influenza virus replication, although key regulators of lipid-dependent trafficking and virus production remain inadequately defined. This report demonstrates that infection by influenza virus stimulates phospholipase D (PLD) activity and that PLD co-localizes with influenza during infection. Both chemical inhibition and RNA interference of PLD delayed viral entry and reduced viral titers in vitro. Although there may be contributions by both major isoenzymes, the effects on viral infectivity appear to be more dependent on the PLD2 isoenzyme. In vivo, PLD2 inhibition reduced virus titer and correlated with significant increases in transcription of innate antiviral effectors. The reduction in viral titer downstream of PLD2 inhibition was dependent on Rig-I (retinoic acid-inducible gene-1), IRF3, and MxA (myxovirus resistance gene A) but not IRF7. Inhibition of PLD2 accelerated the accumulation of MxA in foci as early as 30 min postinfection. Together these data suggest that PLD facilitates the rapid endocytosis of influenza virus, permitting viral escape from innate immune detection and effectors that are capable of limiting lethal infection.

Visualizing infection spread: dual-color fluorescent reporting of virus-host interactions

Although the molecular mechanisms by which host cells defend themselves against viral infection have been studied in great depth, and countermeasures viruses employ to suppress such defensive responses have been widely documented, relatively little attention has been devoted toward elucidating how such interactions between virus and host are resolved over multiple rounds of infection. Here, we describe the design, synthesis, and validation of a dual-color fluorescent reporter system to study how viral infections spread through a host cell monolayer and how the cellular innate immune system mounts an antiviral response. We employed recombinant, red fluorescent protein expressing mutants of a prototypical RNA virus, vesicular stomatitis virus to enable identification and tracking of infected cells. Further, we generated stable reporter cells that use green fluorescent protein to report on the expression of IFIT2, an interferon stimulated gene involved in the interference of viral protein translation, and a marker of antiviral defense activation. The presence of the fluorescent protein reporters had minimal effects on the normal behavior of the cells or viruses. Moreover, expression of the virus and cell reporters correlated with the kinetics of viral replication and activation of an anti-viral response, respectively. This two-color system enabled us to track and quantify in live cells how viral replication and activation of host defensive responses play out over multiple rounds of infection. Initial study of propagating infections demonstrated that antiviral activation over multiple rounds was critical for slowing and ultimately halting the spread of infection.

A quantitative infection assay for human type I, II, and III interferon antiviral activities

Background

Upon virus infection, cells secrete a diverse group of antiviral molecules that signal proximal cells to enter into an antiviral state, slowing or preventing viral spread. These paracrine signaling molecules can work synergistically, so measurement of any one antiviral molecule does not reflect the total antiviral activity of the system.

Results

We have developed an antiviral assay based on replication inhibition of an engineered fluorescent vesicular stomatitis virus reporter strain on A549 human lung epithelial cells. Our assay provides a quantitative functional readout of human type I, II, and III interferon activities, and it provides better sensitivity, intra-, and inter-assay reproducibility than the traditional crystal violet based assay. Further, it eliminates cell fixation, rinsing, and staining steps, and is inexpensive to implement.

Conclusions

A dsRed2-strain of vesicular stomatitis virus that is sensitive to type I, II, and III interferons was used to develop a convenient and sensitive assay for interferon antiviral activity. We demonstrate use of the assay to quantify the kinetics of paracrine antiviral signaling from human prostate cancer (PC3) cells in response to viral infection. The assay is applicable to high-throughput screening for anti-viral compounds as well as basic studies of cellular antiviral signaling.

Increased early RNA replication by chimeric West Nile virus W956IC leads to IPS-1-mediated activation of NF-κB and insufficient virus-mediated counteraction of the resulting canonical type I interferon signaling

Although infections with "natural" West Nile virus (WNV) and the chimeric W956IC WNV infectious clone virus produce comparable peak virus yields in type I interferon (IFN) response-deficient BHK cells, W956IC infection produces higher levels of "unprotected" viral RNA at early times after infection. Analysis of infections with these two viruses in IFN-competent cells showed that W956IC activated NF-κB, induced higher levels of IFN-β, and produced lower virus yields than WNV strain Eg101. IPS-1 was required for both increased induction of IFN-β and decreased yields of W956IC. In Eg101-infected cells, phospho-STAT1/STAT2 nuclear translocation was blocked at all times analyzed, while some phospho-STAT1/STAT2 nuclear translocation was still detected at 8 h after infection in W956IC-infected mouse embryonic fibroblasts (MEFs), and early viral protein levels were lower in these cells. A set of additional chimeras was made by replacing various W956IC gene regions with the Eg101 equivalents. As reported previously, for three of these chimeras, the low early RNA phenotype of Eg101 was restored in BHK cells. Analysis of infections with two of these chimeric viruses in MEFs detected lower early viral RNA levels, higher early viral protein levels, lower early IFN-β levels, and higher virus yields similar to those seen after Eg101 infection. The data suggest that replicase protein interactions directly or indirectly regulate genome switching between replication and translation at early times in favor of translation to minimize NF-κB activation and IFN induction by decreasing the amount of unprotected viral RNA, to produce sufficient viral protein to block canonical type I IFN signaling, and to efficiently remodel cell membranes for exponential genome amplification.

Towards multiscale modeling of influenza infection

Aided by recent advances in computational power, algorithms, and higher fidelity data, increasingly detailed theoretical models of infection with influenza A virus are being developed. We review single scale models as they describe influenza infection from intracellular to global scales, and, in particular, we consider those models that capture details specific to influenza and can be used to link different scales. We discuss the few multiscale models of influenza infection that have been developed in this emerging field. In addition to discussing modeling approaches, we also survey biological data on influenza infection and transmission that is relevant for constructing influenza infection models. We envision that, in the future, multiscale models that capitalize on technical advances in experimental biology and high performance computing could be used to describe the large spatial scale epidemiology of influenza infection, evolution of the virus, and transmission between hosts more accurately.

Kinetics of virus production from single cells

The production of virus by infected cells is an essential process for the spread and persistence of viral diseases, the effectiveness of live-viral vaccines, and the manufacture of viruses for diverse applications. Yet despite its importance, methods to precisely measure virus production from cells are lacking. Most methods test infected-cell populations, masking how individual cells behave. Here we measured the kinetics of virus production from single cells. We combined simple steps of liquid-phase infection, serial dilution, centrifugation, and harvesting, without specialized equipment, to track the production of virus particles from BHK cells infected with vesicular stomatitis virus. Remarkably, cell-to-cell differences in latent times to virus release were within a factor of two, while production rates and virus yields spanned over 300-fold, highlighting an extreme diversity in virus production for cells from the same population. These findings have fundamental and technological implications for health and disease.

West nile virus infections suppress early viral RNA synthesis and avoid inducing the cell stress granule response

West Nile virus (WNV) recently became endemic in the United States and is a significant cause of human morbidity and mortality. Natural WNV strain infections do not induce stress granules (SGs), while W956IC (a lineage 2/1 chimeric WNV infectious clone) virus infections produce high levels of early viral RNA and efficiently induce SGs through protein kinase R (PKR) activation. Additional WNV chimeric viruses made by replacing one or more W956IC genes with the lineage 1 Eg101 equivalent in the W956IC backbone were analyzed. The Eg-NS4b+5, Eg-NS1+3+4a, and Eg-NS1+4b+5 chimeras produced low levels of viral RNA at early times of infection and inefficiently induced SGs, suggesting the possibility that interactions between viral nonstructural proteins and/or between viral nonstructural proteins and cell proteins are involved in suppressing early viral RNA synthesis and membrane remodeling during natural WNV strain infections. Detection of exposed viral double-stranded RNA (dsRNA) in W956IC-infected cells suggested that the enhanced early viral RNA synthesis surpassed the available virus-induced membrane protection and allowed viral dsRNA to activate PKR.

Viral epidemics in a cell culture: novel high resolution data and their interpretation by a percolation theory based model

Because of its relevance to everyday life, the spreading of viral infections has been of central interest in a variety of scientific communities involved in fighting, preventing and theoretically interpreting epidemic processes. Recent large scale observations have resulted in major discoveries concerning the overall features of the spreading process in systems with highly mobile susceptible units, but virtually no data are available about observations of infection spreading for a very large number of immobile units. Here we present the first detailed quantitative documentation of percolation-type viral epidemics in a highly reproducible in vitro system consisting of tens of thousands of virtually motionless cells. We use a confluent astroglial monolayer in a Petri dish and induce productive infection in a limited number of cells with a genetically modified herpesvirus strain. This approach allows extreme high resolution tracking of the spatio-temporal development of the epidemic. We show that a simple model is capable of reproducing the basic features of our observations, i.e., the observed behaviour is likely to be applicable to many different kinds of systems. Statistical physics inspired approaches to our data, such as fractal dimension of the infected clusters as well as their size distribution, seem to fit into a percolation theory based interpretation. We suggest that our observations may be used to model epidemics in more complex systems, which are difficult to study in isolation.

Virus infection speeds: theory versus experiment

In order to explain the speed of Vesicular Stomatitis Virus (VSV) infections, we develop a simple model that improves previous approaches to the propagation of virus infections. For VSV infections, we find that the delay time elapsed between the adsorption of a viral particle into a cell and the release of its progeny has a very important effect. Moreover, this delay time makes the adsorption rate essentially irrelevant in order to predict VSV infection speeds. Numerical simulations are in agreement with the analytical results. Our model satisfactorily explains the experimentally measured speeds of VSV infections.

Dynamics of virus spread in the presence of fluid flow

The dynamics of viral infection spread, whether in laboratory cultures or in naturally infected hosts, reflects a coupling of biological and physical processes that remain to be fully elucidated. Biological processes include the kinetics of virus growth in infected cells while physical processes include transport of virus progeny from infected cells, where they are produced, to susceptible cells, where they initiate new infections. Mechanistic models of infection spread have been widely developed for systems where virus growth is coupled with transport of virus particles by diffusion, but they have yet to be developed for systems where viruses move under the influence of fluid flows. Recent experimental observations of flow-enhanced infection spread in laboratory cultures motivate here the development of initial continuum and discrete virus-particle models of infection spread. The magnitude of a dimensionless group, the Damköhler number, shows how parameters that characterize particle adsorption to cells, strain rates that reflect flow profiles, and diffusivities of virus particles combine to influence the spatial pattern of infection spread.

Infection on a chip: a microscale platform for simple and sensitive cell-based virus assays

The plaque assay has long served as the "gold standard" to measure virus infectivity and test antiviral drugs, but the assay is labor-intensive, lacks sensitivity, uses excessive reagents, and is hard to automate. Recent modification of the assay to exploit flow-enhanced virus spread with quantitative imaging has increased its sensitivity. Here we performed flow-enhanced infection assays in microscale channels, employing passive fluid pumping to inoculate cell monolayers with virus and drive infection spread. Our test of an antiviral drug (5-fluorouracil) against vesicular stomatitis virus infections of BHK cell monolayers yielded a two-fold improvement in sensitivity, relative to the standard assay based on plaque counting. The reduction in scale, simplified fluid handling, image-based quantification, and higher assay sensitivity will enable infection measurements for high-throughput drug screening, sero-conversion testing, and patient-specific diagnosis of viral infections.

Agent-based modeling of host-pathogen systems: The successes and challenges

Agent-based models have been employed to describe numerous processes in immunology. Simulations based on these types of models have been used to enhance our understanding of immunology and disease pathology. We review various agent-based models relevant to host–pathogen systems and discuss their contributions to our understanding of biological processes. We then point out some limitations and challenges of agent-based models and encourage efforts towards reproducibility and model validation.

Multiple-hit inhibition of infection by defective interfering particles

Defective interfering particles (DIPs) are virus-like particles that arise during virus growth, fail to grow in the absence of virus, and replicate at the expense of virus during co-infections. The inhibitory effects of DIPs on virus growth are well established, but little is known about how DIPs influence their own growth. Here vesicular stomatitis virus (VSV) and its DIPs were used to co-infect BHK cells, and the effect of DIP dose on virus and DIP production was measured using a yield-reduction assay. The resulting dose-response data were used to fit and evaluate mathematical models that employed different assumptions. Our analysis supports a multiple-hit process where DIPs inhibit or promote virus and DIP production, depending on dose. Specifically, three regimes of co-infection were apparent: (i) low DIP - where both virus and DIPs are amplified, (ii) medium DIP - where amplification of both virus and DIPs is inhibited, and (iii) high DIP - with limited recovery of virus production and further inhibition of DIP growth. In addition, serial-passage infections enabled us to estimate the frequency of de novo DIP generation during virus amplification. Our combined experiments and models provide a means to understand better how DIPs quantitatively impact the growth of viruses and the spread of their infections.

Enhanced adaptation of vesicular stomatitis virus in cells infected with vaccinia virus

Infections involving different viruses (multiple infections) are common in nature and can take place between different strains of the same virus or between different virus species, including DNA and RNA viruses. The influence of multiple infections on viral evolution has been previously studied using different populations of the same virus. Here, we took a step forward by studying the evolution of an RNA virus (vesicular stomatitis virus, VSV) in the presence of a resident DNA virus (vaccinia virus, VV). Cell cultures were infected with a constant amount of VV, and VSV was added at four different post-VV-inoculation times and four different population sizes. The results showed that the presence of VV accelerates the adaptation of VSV to a cellular environment, especially at high population sizes. The effect of VV on VSV evolution was stronger when cells were incubated for longer times with VV prior to the addition of VSV. Our results suggest that cooperation between the two viruses rather than competition might be responsible for the enhanced rate of adaptation of VSV. Further studies are needed to discern whether infections involving different viruses could have an increased ability to escape antiviral strategies.

Susceptibility to viral infection is enhanced by stable expression of 3A or 3AB proteins from foot-and-mouth disease virus

The foot-and-mouth disease virus (FMDV) 3A protein is involved in virulence and host range. A distinguishing feature of FMDV 3B among picornaviruses is that three non-identical copies are encoded in the viral RNA and required for optimal replication in cell culture. Here, we have studied the involvement of the 3AB region on viral infection using constitutive and transient expression systems. BHK-21 stably transformed clones expressed low levels of FMDV 3A or 3A(B) proteins in the cell cytoplasm. Transformed cells stably expressing these proteins did not exhibit inner cellular rearrangements detectable by electron microscope analysis. Upon FMDV infection, clones expressing either 3A alone or 3A(B) proteins showed a significant increase in the percentage of infected cells, the number of plaque forming units and the virus yield. The 3A-enhancing effect was specific for FMDV as no increase in viral multiplication was observed in transformed clones infected with another picornavirus, encephalomyocarditis virus, or the negative-strand RNA virus vesicular stomatitis virus. A potential role of 3A protein in viral RNA translation was discarded by the lack of effect on FMDV IRES-dependent translation. Increased viral susceptibility was not caused by a released factor; neither the supernatant of transformed clones nor the addition of purified 3A protein to the infection medium was responsible for this effect. Unlike stable expression, high levels of 3A or 3A(B) protein transient expression led to unspecific inhibition of viral infection. Therefore, the effect observed on viral yield, which inversely correlated with the intracellular levels of 3A protein, suggests a transacting role operating on the FMDV multiplication cycle.

Image-guided modeling of virus growth and spread

Although many tools of cellular and molecular biology have been used to characterize single intracellular cycles of virus growth, few culture methods exist to study the dynamics of spatially spreading viruses over multiple generations. We have previously developed a method that addresses this need by tracking the spread of focal infections using immunocytochemical labeling and digital imaging. Here, we build reaction-diffusion models to account for spatio-temporal patterns formed by the spreading viral infection front as well as data from a single cycle of virus growth (one-step growth). Systems with and without the interferon-mediated antiviral response of the host cells are considered. Dynamic images of the spreading infections guide iterative model refinement steps that lead to reproduction of all of the salient features contained in the images, not just the velocity of the infection front. The optimal fits provide estimates for key parameters such as virus-host binding and the production rate of interferon. For the examined data, highly-lumped infection models that ignore the one-step growth dynamics provide a comparable fit to models that more accurately account for these dynamics, highlighting the fact that increased model complexity does not necessarily translate to improved fit. This work demonstrates how model building can facilitate the interpretation of experiments by highlighting contributions from both biological and methodological factors.

Superior infectivity for mosquito vectors contributes to competitive displacement among strains of dengue virus

BACKGROUND

Competitive displacement of a weakly virulent pathogen strain by a more virulent strain is one route to disease emergence. However the mechanisms by which pathogens compete for access to hosts are poorly understood. Among vector-borne pathogens, variation in the ability to infect vectors may effect displacement. The current study focused on competitive displacement in dengue virus serotype 3 (DENV3), a mosquito-borne pathogen of humans. In Sri Lanka in the 1980's, a native DENV3 strain associated with relatively mild dengue disease was displaced by an invasive DENV3 strain associated with the most severe disease manifestations, dengue hemorrhagic fever/dengue shock syndrome (DHF/DSS), resulting in an outbreak of DHF/DSS. Here we tested the hypothesis that differences between the invasive and native strain in their infectivity for Aedes aegypti mosquitoes, the primary vector of DENV, contributed to the competitive success of the invasive strain

RESULTS

To be transmitted by a mosquito, DENV must infect and replicate in the midgut, disseminate into the hemocoel, infect the salivary glands, and be released into the saliva. The ability of the native and invasive DENV3 strains to complete the first three steps of this process in Aedes aegypti mosquitoes was measured in vivo. The invasive strain infected a similar proportion of mosquitoes as the native strain but replicated to significantly higher titers in the midgut and disseminated with significantly greater efficiency than the native strain. In contrast, the native and invasive strain showed no significant difference in replication in cultured mosquito, monkey or human cells.

CONCLUSION

The invasive DENV3 strain infects and disseminates in Ae. aegypti more efficiently than the displaced native DENV3 strain, suggesting that the invasive strain is transmitted more efficiently. Replication in cultured cells did not adequately characterize the known phenotypic differences between native and invasive DENV3 strains. Infection dynamics within the vector may have a significant impact on the spread and replacement of dengue virus lineages.

A quantitative comet assay: imaging and analysis of virus plaques formed with a liquid overlay

Although the plaque assay defines a "gold-standard" for measuring virus infectivity, its reliance on plaque counting limits its sensitivity. When the assay is performed with a liquid overlay, instead of agar overlay, spontaneous flows can promote a uni-directional spread of infection, creating elongated regions of cytopathology that resemble comets. As a model system comet and plaque cultures of vesicular stomatitis virus (VSV) on baby hamster kidney (BHK-21) cells were compared. Host-cell monolayers were infected with VSV particles, incubated 15 h in the presence of liquid or agar overlays and stained. VSV formed significantly larger comets than plaques, consistent with a mechanism of flow-enhanced spread. When antiviral drug (5-fluorouracil) was incorporated into the liquid overlay, comet sizes were reduced in a dose-dependent manner. Images of infected monolayers, acquired using a simple digital scanner, enabled a quantification of the inhibitory effect of the drug on infectivity. The resulting measure of drug susceptibility was found to be 18-fold more sensitive than the IC(50) measure attained by the traditional plaque-reduction assay. This quantitative comet assay has the potential to similarly enhance the sensitivity of infection measures for other plaque-forming viruses.

Cell-type-specific growth restriction of vesicular stomatitis virus polR mutants is linked to defective viral polymerase function

Vesicular stomatitis virus polR mutants synthesize defective RNA replication products in vitro and display growth restriction in some cultured cells (J. L. Chuang, R. L. Jackson, and J. Perrault, Virology 229:57-67, 1997). We show here that a recombinant virus carrying the polR N protein mutation (R179H) yielded approximately 100-fold- and approximately 40-fold-lower amounts of infectious virus than the wild type in mouse L-929 and rat 3Y1 cells, respectively, but only approximately 3-fold less in hamster BHK cells. Virus genome accumulation was inhibited 6- to 10-fold in restricting cells, but transcription was not affected. No defect in encapsidation of replication products was detected, but virus protein accumulation was reduced two- to threefold in both restricting and nonrestricting cells. polR virus particles released from the latter were 5- to 10-fold less infectious than the wild type but showed no difference in protein composition. Phosphorylation of the alpha subunit of eukaryotic translation initiation factor 2 (eIF-2alpha) was enhanced approximately 3-fold in polR versus wild-type virus-infected L-929 cells, but neither inhibition of host gene transcription nor inhibition of double-stranded RNA (dsRNA)-activated protein kinase showed significant effects on restriction. Conditioned medium studies revealed no evidence for secretion of antiviral factors from restricting cells. We conclude that the block in polR growth is due to the combined effect of reduced genome replication and lower infectivity of released virus particles and may be due to overproduction of dsRNA. An accompanying paper (D. Ostertag, T. M. Hoblitzell-Ostertag, and J. Perrault, J. Virol. 81:503-513, 2007) provides compelling evidence for the role of dsRNA in this unique restriction phenomenon.

Spatial patterns of protein expression in focal infections of human cytomegalovirus

Human cytomegalovirus (HCMV) is a medically significant human pathogen that infects a wide range of cell and tissue types. During infection, HCMV activates a variety of signal transduction pathways that induce profound changes in cellular processes and dramatically affect cellular gene expression patterns. To better define how these virus-host interactions affect the local microenvironment and influence the spatial and temporal spread of HCMV, we initiated HCMV focal infections on normal human dermal fibroblast monolayers and monitored viral gene expression patterns and infection spread over 45 days. To establish baseline temporal measurements of HCMV infection and spread in cell monolayers, we characterized the influence of three experimental variables on viral gene expression: cell plating density, the presence of serum, and neutralization of cellular antiviral responses with an antibody against interferon-beta. We found that high cell plating density or the inclusion of serum correlated with enhanced HCMV infection spread. Dramatic differences in the expression pattern of the viral immediate early 2 (IE2) gene were observed under these conditions as compared to low plating density or the absence of serum. In the latter case round, uniform foci were observed with a clear wave of IE2 expression visible in advance of a late stage viral protein, envelope glycoprotein B. By contrast, larger irregular foci with arms of IE2 expression were observed in the presence of serum. Addition of the antibody had little effect on the rate of spread, which is consistent with the knowledge that HCMV represses antiviral responses during infection. This experimental system provides a useful means to visualize and quantify complex virus-host interactions.

Probing the effects of the well-mixed assumption on viral infection dynamics

Viral kinetics have been extensively studied in the past through the use of spatially well-mixed ordinary differential equations describing the time evolution of the diseased state. However, emerging spatial structures such as localized populations of dead cells might adversely affect the spread of infection, similar to the manner in which a counter-fire can stop a forest fire from spreading. In a previous publication [Beauchemin, C., Samuel, J., Tuszynski, J., 2005. A simple cellular automaton model for influenza A viral infections. J. Theor. Biol. 232(2), 223-234], a simple two-dimensional cellular automaton model was introduced and shown to be accurate enough to model an uncomplicated infection with influenza A. Here, this model is used to investigate the effects of relaxing the well-mixed assumption. Particularly, the effects of the initial distribution of infected cells, the regeneration rule for dead epithelial cells, and the proliferation rule for immune cells are explored and shown to have an important impact on the development and outcome of the viral infection in our model.