Cell-to-Cell Variation in Defective Virus Expression and Effects on Host Responses during Influenza Virus Infection

Influenza A defective viral genome production is altered by metabolites, metabolic signaling molecules, and cyanobacteria extracts

RNA virus infections are composed of a diverse mix of viral genomes that arise from low fidelity in replication within cells. The interactions between "defective" and full-length viral genomes have been shown to shape pathogenesis, leading to intense research into employing these to develop novel antivirals. In particular, Influenza A defective viral genomes (DVGs) have been associated with milder clinical outcomes. Yet, the full potential of DVGs as broad-spectrum antivirals remains untapped due to the unknown mechanisms of their de novo production. Much of the research into the factors affecting defective viral genome production has focused on the virus, while the role of the host has been neglected. We recently showed that altering host cell metabolism away from pro-growth pathways using alpelisib increased the production of Influenza A defective viral genomes. To uncover other drugs that could induce infections to create more DVGs, we subjected active influenza infections of the two circulating human subtypes (A/H1N1 & A/H3N2) to a screen of metabolites, metabolic signaling molecules, and cyanobacteria-derived biologics, after which we quantified the defective viral genomes (specifically deletion-containing viral genomes, DelVGs) and total viral genomes using third generation long-read sequencing. Here we show that metabolites and signaling molecules of host cell central carbon metabolism can significantly alter DelVG production early in Influenza A infection. Adenosine, emerged as a potent inducer of defective viral genomes, significantly amplifying DelVG production across both subtypes. Insulin had similar effects, albeit subtype-specific, predominantly enhancing polymerase segment DVGs in TX12 infections. Tricarboxylic Acid (TCA) cycle inhibitors 4-octyl itaconate and UK5099, along with the purine analog favipiravir, increased total viral genome production across subtypes. Cyanobacterial extracts primarily affected DVG and total viral genome production in TX12, with a specific, almost complete shutdown of influenza antigenic segments. These results underscore the influence of host metabolic pathways on DVG production and suggest new avenues for antiviral intervention, including PI3K-AKT and Ras-MAPK signaling pathways, TCA cycle metabolism, purine-pyrimidine metabolism, polymerase inhibition, and cyanotherapeutic approaches. More broadly, our findings suggest that the social interactions observed between defective and full-length viral genomes, depend not only on the viral actors, but can be altered by the stage provided by the host. Our study advances our fundamental understanding of DVG production mechanisms and highlights the potential of targeting host metabolism to develop broad-spectrum influenza therapeutics.

Influenza A defective viral genomes and non-infectious particles are increased by host PI3K inhibition via anti-cancer drug alpelisib

RNA viruses produce abundant defective viral genomes during replication, setting the stage for interactions between viral genomes that alter the course of pathogenesis. Harnessing these interactions to develop antivirals has become a recent goal of intense research focus. Despite decades of research, the mechanisms that regulate the production and interactions of Influenza A defective viral genomes are still unclear. The role of the host is essentially unexplored; specifically, it remains unknown whether host metabolism can influence the formation of defective viral genomes and the particles that house them. To address this question, we manipulated host cell anabolic signaling activity and monitored the production of defective viral genomes and particles by A/H1N1 and A/H3N2 strains, using a combination of single-cell immunofluorescence quantification, third-generation long-read sequencing, and the cluster-forming assay, a method we developed to titer defective and fully-infectious particles simultaneously. Here we show that alpelisib (Piqray), a highly selective inhibitor of mammalian Class 1a phosphoinositide-3 kinase (PI3K) receptors, significantly changed the proportion of defective particles and viral genomes (specifically deletion-containing viral genomes) in a strain-specific manner, under conditions that minimize multiple cycles of replication. Alpelisib pre-treatment of cells led to an increase in defective particles in the A/H3N2 strain, while the A/H1N1 strain showed a decrease in total viral particles. In the same infections, we found that defective viral genomes of polymerase and antigenic segments increased in the A/H1N1 strain, while the total particles decreased suggesting defective interference. We also found that the average deletion size in polymerase complex viral genomes increased in both the A/H3N2 and A/H1N1 strains. The A/H1N1 strain, additionally showed a dose-dependent increase in total number of defective viral genomes. In sum, we provide evidence that host cell metabolism can increase the production of defective viral genomes and particles at an early stage of infection, shifting the makeup of the infection and potential interactions among virions. Given that Influenza A defective viral genomes can inhibit pathogenesis, our study presents a new line of investigation into metabolic states associated with less severe flu infection and the potential induction of these states with metabolic drugs.

Stromal Rigidity Stress Accelerates Pancreatic Intraepithelial Neoplasia Progression and Chromosomal Instability via Nuclear Protein Tyrosine Kinase 2 Localization

Because the mechanotransduction by stromal stiffness stimulates the rupture and repair of the nuclear envelope in pancreatic progenitor cells, accumulated genomic aberrations are under selection in the tumor microenvironment. Analysis of cell growth, micronuclei, and phosphorylated Ser-139 residue of the histone variant H2AX (γH2AX) foci linked to mechanotransduction pressure in vivo during serial orthotopic passages of mouse KrasLSL-G12D/+;Trp53flox/flox;Pdx1-Cre (KPC) cancer cells in the tumor and in migrating through the size-restricted 3-μm micropores. To search for pancreatic cancer cell-of-origin, analysis of single-cell data sets revealed that the extracellular matrix shaped an alternate route of acinar-ductal transdifferentiation of acinar cells into topoisomerase II α (TOP2A)-overexpressing cancer cells and derived subclusters with copy number amplifications in MYC-PTK2 (protein tyrosine kinase 2) locus and PIK3CA. High-PTK2 expression is associated with 171 differentially methylated CpG loci, 319 differentially expressed genes, and poor overall survival in The Cancer Genome Atlas-Pancreatic Adenocarcinoma cohort. Abolished RGD-integrin signaling by disintegrin KG blocked the PTK2 phosphorylation, increased cancer apoptosis, decreased vav guanine nucleotide exchange factor 1 (VAV1) expression, and prolonged overall survival in the KPC mice. Reduction of α-smooth muscle actin deposition in the CD248 knockout KPC mice remodeled the tissue stroma and down-regulated TOP2A expression in the epithelium. In summary, stromal stiffness induced the onset of cancer cells-of-origin by ectopic TOP2A expression, and the genomic amplification of MYC-PTK2 locus via alternative transdifferentiation of pancreatic progenitor cells is the vulnerability useful for disintegrin KG treatment.

Structure and dynamics of enterovirus genotype networks

Like all biological populations, viral populations exist as networks of genotypes connected through mutation. Mapping the topology of these networks and quantifying population dynamics across them is crucial to understanding how populations adapt to changes in their selective environment. The influence of mutational networks is especially profound in viral populations that rapidly explore their mutational neighborhoods via high mutation rates. Using a single-cell sequencing method, scRNA-seq-enabled acquisition of mRNA and consensus haplotypes linking individual genotypes and host transcriptomes (SEARCHLIGHT), we captured and assembled viral haplotypes from hundreds of individual infected cells, revealing the complexity of viral population structures. We obtained these genotypes in parallel with host cell transcriptome information, enabling us to link host cell transcriptional phenotypes to the genetic structures underlying virus adaptation. Our examination of these structures reveals the common evolutionary dynamics of enterovirus populations and illustrates how viral populations reach through mutational "tunnels" to span evolutionary landscapes and maintain connection with multiple adaptive genotypes simultaneously.

Selection of nonstandard viral genomes during the evolution of RNA viruses: A virus survival strategy or a pesky inconvenience?

RNA viruses are some of the most successful biological entities due their ability to adapt and evolve. Despite their small genome and parasitic nature, RNA viruses have evolved many mechanisms to ensure their survival and maintenance in the host population. We propose that one of these mechanisms of survival is the generation of nonstandard viral genomes (nsVGs) that accumulate during viral replication. NsVGs are often considered to be accidental defective byproducts of the RNA virus replication, but their ubiquity and the plethora of roles they have during infection indicate that they are an integral part of the virus life cycle. Here we review the different types of nsVGs and discuss how their multiple roles during infection could be beneficial for RNA viruses to be maintained in nature. By shifting our perspectives on what makes a virus successful, we posit that nsVG generation is a conserved phenomenon that arose during RNA virus evolution as an essential component of a healthy virus community.

Exposure of negative-sense viral RNA in the cytoplasm initiates innate immunity to West Nile virus

For many RNA viruses, immunity is triggered when RIG-I-like receptors (RLRs) detect viral RNA. However, only a minority of infected cells undergo innate immune activation. By examining these "first responder" cells during West Nile virus infection, we found that specific accumulation of anti- genomic negative-sense viral RNA (-vRNA) underlies innate immune activation and that RIG-I preferentially interacts with -vRNA. However, flaviviruses sequester -vRNA into membrane-bound replication compartments away from cytosolic sensors. We found that single-stranded -vRNA accumulates outside of replication compartments in "first responder" cells, rendering it accessible to RLRs. Exposure of this -vRNA occurs at late timepoints of infection, is linked to viral assembly, and depends on the expression of viral structural proteins. These findings reveal that while most infected cells replicate high levels of vRNA, release of -vRNA from replication compartments during assembly occurs at low frequency and is critical for initiation of innate immunity during flavivirus infection.

Defective viral genomes: advances in understanding their generation, function, and impact on infection outcomes

Defective viral genomes (DVGs) are truncated derivatives of their parental viral genomes generated during an aberrant round of viral genomic replication. Distinct classes of DVGs have been identified in most families of both positive- and negative-sense RNA viruses. Importantly, DVGs have been detected in clinical samples from virally infected individuals and an emerging body of association studies implicates DVGs in shaping the severity of disease caused by viral infections in humans. Consequently, there is growing interest in understanding the molecular mechanisms of de novo DVG generation, how DVGs interact with the innate immune system, and harnessing DVGs as novel therapeutics and vaccine adjuvants to attenuate viral pathogenesis. This minireview focuses on single-stranded RNA viruses (excluding retroviridae), and summarizes the current knowledge of DVG generation, the functions and diversity of DVG species, the roles DVGs play in influencing disease progression, and their application as antivirals and vaccine adjuvants.

The reproduction number and its probability distribution for stochastic viral dynamics

We consider stochastic models of individual infected cells. The reproduction number, R, is understood as a random variable representing the number of new cells infected by one initial infected cell in an otherwise susceptible (target cell) population. Variability in R results partly from heterogeneity in the viral burst size (the number of viral progeny generated from an infected cell during its lifetime), which depends on the distribution of cellular lifetimes and on the mechanism of virion release. We analyse viral dynamics models with an eclipse phase: the period of time after a cell is infected but before it is capable of releasing virions. The duration of the eclipse, or the subsequent infectious, phase is non-exponential, but composed of stages. We derive the probability distribution of the reproduction number for these viral dynamics models, and show it is a negative binomial distribution in the case of constant viral release from infectious cells, and under the assumption of an excess of target cells. In a deterministic model, the ultimate in-host establishment or extinction of the viral infection depends entirely on whether the mean reproduction number is greater than, or less than, one, respectively. Here, the probability of extinction is determined by the probability distribution of R, not simply its mean value. In particular, we show that in some cases the probability of infection is not an increasing function of the mean reproduction number.

The use of single-cell RNA-seq to study heterogeneity at varying levels of virus-host interactions

The outcome of viral infection depends on the diversity of the infecting viral population and the heterogeneity of the cell population that is infected. Until almost a decade ago, the study of these dynamic processes during viral infection was challenging and limited to certain targeted measurements. Presently, with the use of single-cell sequencing technology, the complex interface defined by the interactions of cells with infecting virus can now be studied across the breadth of the transcriptome in thousands of individual cells simultaneously. In this review, we will describe the use of single-cell RNA sequencing (scRNA-seq) to study the heterogeneity of viral infections, ranging from individual virions to the immune response between infected individuals. In addition, we highlight certain key experimental limitations and methodological decisions that are critical to analyzing scRNA-seq data at each scale.

Influenza virus transcription and progeny production are poorly correlated in single cells

The ultimate success of a viral infection at the cellular level is determined by the number of progeny virions produced. However, most single-cell studies of infection quantify the expression of viral transcripts and proteins, rather than the amount of progeny virions released from infected cells. Here, we overcome this limitation by simultaneously measuring transcription and progeny production from single influenza virus-infected cells by embedding nucleotide barcodes in the viral genome. We find that viral transcription and progeny production are poorly correlated in single cells. The cells that transcribe the most viral mRNA do not produce the most viral progeny and often represent aberrant infections that fail to express the influenza NS gene. However, only some of the discrepancy between transcription and progeny production can be explained by viral gene absence or mutations: there is also a wide range of progeny production among cells infected by complete unmutated virions. Overall, our results show that viral transcription is a relatively poor predictor of an infected cell's contribution to the progeny population.

Broad-Spectrum Antiviral Activity of Influenza A Defective Interfering Particles against Respiratory Syncytial, Yellow Fever, and Zika Virus Replication In Vitro

New broadly acting and readily available antiviral agents are needed to combat existing and emerging viruses. Defective interfering particles (DIPs) of influenza A virus (IAV) are regarded as promising options for the prevention and treatment of IAV infections. Interestingly, IAV DIPs also inhibit unrelated viral infections by stimulating antiviral innate immunity. Here, we tested the ability of IAV DIPs to suppress respiratory syncytial, yellow fever and Zika virus infections in vitro. In human lung (A549) cells, IAV DIP co-infection inhibited the replication and spread of all three viruses. In contrast, we observed no antiviral activity in Vero cells, which are deficient in the production of interferon (IFN), demonstrating its importance for the antiviral effect. Further, in A549 cells, we observed an enhanced type-I and type-III IFN response upon co-infection that appears to explain the antiviral potential of IAV DIPs. Finally, a lack of antiviral activity in the presence of the Janus kinase 1/2 (JAK1/2) inhibitor ruxolitinib was detected. This revealed a dependency of the antiviral activity on the JAK/signal transducers and activators of transcription (STAT) signaling pathway. Overall, this study supports the notion that IAV DIPs may be used as broad-spectrum antivirals to treat infections with a variety of IFN-sensitive viruses, particularly respiratory viruses.

Influenza Defective Interfering Virus Promotes Multiciliated Cell Differentiation and Reduces the Inflammatory Response in Mice

Influenza defective interfering (DI) viruses have long been considered promising antiviral candidates because of their ability to interfere with replication-competent viruses and induce antiviral immunity. However, the mechanisms underlying DI-mediated antiviral immunity have not been extensively explored. Here, we demonstrated the interferon (IFN)-independent protection conferred by the influenza DI virus against homologous virus infection in mice deficient in type I and III IFN signaling. We identified unique host signatures responding to DI coinfection by integrating transcriptional and posttranscriptional regulatory data. DI-treated mice exhibited reduced viral transcription, less intense inflammatory and innate immune responses, and primed multiciliated cell differentiation in their lungs at an early stage of infection, even in the absence of type I or III IFNs. This increased multiciliogenesis could also be detected at the protein level via the immunofluorescence staining of lung tissue from DI-treated mice. Overall, our study provides mechanistic insight into the protection mediated by DIs, implying a unifying theme involving inflammation and multiciliogenesis in maintaining respiratory homeostasis and revealing their IFN-independent antiviral activity. IMPORTANCE During replication, the influenza virus generates genetically defective viruses. These are found in natural infections as part of the virus population within the infected host. Some versions of these defective viruses are thought to have protective effects through their interference with replication-competent viruses and induction of antiviral immunity. To better determine the mechanisms underlying the protective effects of these defective interfering (DI) viruses, we tested a DI that we previously identified in vitro with mice. Mice that were infected with a mix of wild-type influenza and DI viruses had less intense inflammatory and innate immune responses than did mice that were infected with the wild-type virus only, even when type I or III interferons, which are cytokines that play a prominent role in defending the respiratory epithelial barrier, were absent. More interestingly, the DI-infected mice had primed multiciliated cell differentiation in their lungs, indicating the potential promotion of epithelial repair by DIs.

Quantifying the effect of defective viral genomes in respiratory syncytial virus infections

Defective viral genomes (DVGs) are viral genomes that contain only a partial viral RNA and so cannot replicate within cells on their own. If a cell containing DVGs is subsequently infected with a complete viral genome, the DVG can then use the missing proteins expressed by the full genome in order to replicate itself. Since the cell is producing defective genomes, it has less resources to produce fully functional virions and thus release of complete virions is often suppressed. Here, we use data from challenge studies of respiratory syncytial virus (RSV) in healthy adults to quantify the effect of DVGs. We use a mathematical model to fit the data, finding that late onset of DVGs and prolonged DVG detection are associated with lower infection rates and higher clearance rates. This result could have implications for the use of DVGs as a therapeutic.

Multiscale network analysis identifies potential receptors for SARS-CoV-2 and reveals their tissue-specific and age-dependent expression

The coronavirus disease 2019 (COVID-19) pandemic has affected tens of millions of individuals and caused hundreds of thousands of deaths worldwide. Here, we present a comprehensive, multiscale network analysis of the transcriptional response to the virus. In particular, we focused on key regulators, cell receptors, and host processes that were hijacked by the virus for its advantage. ACE2-controlled processes involved CD300e (a TYROBP receptor) as a key regulator and the activation of IL-2 pro-inflammatory cytokine signaling. We further investigated the age dependency of such receptors in different tissues. In summary, this study provides novel insights into the gene regulatory organization during the SARS-CoV-2 infection and the tissue-specific, age-dependent expression of the cell receptors involved in COVID-19.

Maximal interferon induction by influenza lacking NS1 is infrequent owing to requirements for replication and export

Influenza A virus exhibits high rates of replicative failure due to a variety of genetic defects. Most influenza virions cannot, when acting as individual particles, complete the entire viral life cycle. Nevertheless influenza is incredibly successful in the suppression of innate immune detection and the production of interferons, remaining undetected in >99% of cells in tissue-culture models of infection. Notably, the same variation that leads to replication failure can, by chance, inactivate the major innate immune antagonist in influenza A virus, NS1. What explains the observed rarity of interferon production in spite of the frequent loss of this, critical, antagonist? By studying how genetic and phenotypic variation in a viral population lacking NS1 correlates with interferon production, we have built a model of the "worst-case" failure from an improved understanding of the steps at which NS1 acts in the viral life cycle to prevent the triggering of an innate immune response. In doing so, we find that NS1 prevents the detection of de novo innate immune ligands, defective viral genomes, and viral export from the nucleus, although only generation of de novo ligands appears absolutely required for enhanced detection of virus in the absence of NS1. Due to this, the highest frequency of interferon production we observe (97% of infected cells) requires a high level of replication in the presence of defective viral genomes with NS1 bearing an inactivating mutation that does not impact its partner encoded on the same segment, NEP. This is incredibly unlikely to occur given the standard variation found within a viral population, and would generally require direct, artificial, intervention to achieve at an appreciable rate. Thus from our study, we procure at least a partial explanation for the seeming contradiction between high rates of replicative failure and the rarity of the interferon response to influenza infection.

Tissue factor links inflammation, thrombosis, and senescence in COVID-19

COVID-19 is a highly contagious respiratory infection caused by the SARS-CoV-2 virus. The infected lung epithelial cells secrete a group of chemokines and cytokines, which triggers harmful cytokine storms and hyper-thrombotic responses. Recent studies have proposed that viral-induced senescence is responsible for cytokine release and inflammation in COVID-19 patients. However, it is unknown whether cellular senescence is commonly triggered after viral infection and how inflammation and thrombosis, hyper-activated in these patients, are functionally connected. To address these questions, we conducted a bioinformatics-based meta-analysis using single-cell and bulk RNA sequencing datasets obtained from human patient studies, animal models, and cell lines infected with SARS-CoV-2 and other respiratory viruses. We found that the senescence phenotype is robustly upregulated in most SARS-CoV-2-infected patients, especially in the infected lung epithelial cells. Notably, the upregulation of Tissue factor (F3), a key initiator of the extrinsic blood coagulation pathway, occurs concurrently with the upregulation of the senescence-associated secretory phenotype (SASP) factors. Furthermore, F3 levels are positively correlated with the senescence and hyper-coagulation gene signatures in COVID-19 patients. Together, these data demonstrate the prevalence of senescence in respiratory viral infection and suggest F3 as a critical link between inflammation, thrombosis, and senescence in these disease states.

Levels of Influenza A Virus Defective Viral Genomes Determine Pathogenesis in the BALB/c Mouse Model

Defective viral genomes (DVGs), which are generated by the viral polymerase in error during RNA replication, can trigger innate immunity and are implicated in altering the clinical outcome of infection. Here, we investigated the impact of DVGs on innate immunity and pathogenicity in a BALB/c mouse model of influenza virus infection. We generated stocks of influenza viruses containing the internal genes of an H5N1 virus that contained different levels of DVGs (indicated by different genome-to-PFU ratios). In lung epithelial cells, the high-DVG stock was immunostimulatory at early time points postinfection. DVGs were amplified during virus replication in myeloid immune cells and triggered proinflammatory cytokine production. In the mouse model, infection with the different virus stocks produced divergent outcomes. The high-DVG stock induced an early type I interferon (IFN) response that limited viral replication in the lungs, resulting in minimal weight loss. In contrast, the virus stock with low levels of DVGs replicated to high titers and amplified DVGs over time, resulting in elevated levels of proinflammatory cytokines accompanied by rapid weight loss and increased morbidity and mortality. Our results suggest that the timing and levels of immunostimulatory DVGs generated during infection contribute to H5N1 pathogenesis. IMPORTANCE Mammalian infections with highly pathogenic avian influenza viruses (HPAIVs) cause severe disease associated with excessive proinflammatory cytokine production. Aberrant replication products, such as defective viral genomes (DVGs), can stimulate the antiviral response, and cytokine induction is associated with their emergence in vivo. We show that stocks of a recombinant virus containing HPAIV internal genes that differ in their amounts of DVGs have vastly diverse outcomes in a mouse model. The high-DVG stock resulted in extremely mild disease due to suppression of viral replication. Conversely, the stock that contained low DVGs but rapidly accumulated DVGs over the course of infection led to severe disease. Therefore, the timing of DVG amplification and proinflammatory cytokine production impact disease outcome, and these findings demonstrate that not all DVG generation reduces viral virulence. This study also emphasizes the crucial requirement to examine the quality of virus preparations regarding DVG content to ensure reproducible research.

Comprehensive single cell analysis of pandemic influenza A virus infection in the human airways uncovers cell-type specific host transcriptional signatures relevant for disease progression and pathogenesis

The respiratory epithelium constitutes the first line of defense against invading respiratory pathogens, such as the 2009 pandemic strain of influenza A virus (IAV, H1N1pdm09), and plays a crucial role in the host antiviral response to infection. Despite its importance, however, it remains unknown how individual cell types within the respiratory epithelium respond to IAV infection or how the latter may influence IAV disease progression and pathogenesis. Here, we used single cell RNA sequencing (scRNA-seq) to dissect the host response to IAV infection in its natural target cells. scRNA-seq was performed on human airway epithelial cell (hAEC) cultures infected with either wild-type pandemic IAV (WT) or with a mutant version of IAV (NS1_R38A) that induced a robust innate immune response. We then characterized both the host and viral transcriptomes of more than 19,000 single cells across the 5 major cell types populating the human respiratory epithelium. For all cell types, we observed a wide spectrum of viral burden among single infected cells and a disparate host response between infected and bystander populations. Interestingly, we also identified multiple key differences in the host response to IAV among individual cell types, including high levels of pro-inflammatory cytokines and chemokines in secretory and basal cells and an important role for luminal cells in sensing and restricting incoming virus. Multiple infected cell types were shown to upregulate interferons (IFN), with type III IFNs clearly dominating the antiviral response. Transcriptional changes in genes related to cell differentiation, cell migration, and tissue repair were also identified. Strikingly, we also detected a shift in viral host cell tropism from non-ciliated cells to ciliated cells at later stages of infection and observed major changes in the cellular composition. Microscopic analysis of both WT and NS1_R38A virus-infected hAECs at various stages of IAV infection revealed that the transcriptional changes we observed at 18 hpi were likely driving the downstream histopathological alterations in the airway epithelium. To our knowledge, this is the first study to provide a comprehensive analysis of the cell type-specific host antiviral response to influenza virus infection in its natural target cells – namely, the human respiratory epithelium.

In Vivo Profiling of Individual Multiciliated Cells during Acute Influenza A Virus Infection

Influenza virus infections are thought to be initiated in a small number of cells; however, the heterogeneity across the cellular responses of the epithelial cells during establishment of disease is incompletely understood. Here, we used an H1N1 influenza virus encoding a fluorescent reporter gene, a cell lineage-labeling transgenic mouse line, and single-cell RNA sequencing to explore the range of responses in a susceptible epithelial cell population during an acute influenza A virus (IAV) infection. Focusing on multiciliated cells, we identified a subpopulation that basally expresses interferon-stimulated genes (ISGs), which we hypothesize may be important for the early response to infection. We subsequently found that a population of infected ciliated cells produce most of the ciliated cell-derived inflammatory cytokines, and nearly all bystander ciliated cells induce a broadly antiviral state. From these data together, we propose that variable preexisting gene expression patterns in the initial cells targeted by the virus may ultimately affect the establishment of viral disease. IMPORTANCE Influenza A virus poses a significant threat to public health, and each year, millions of people in the United States alone are exposed to the virus. We do not currently, however, fully understand why some individuals clear the infection asymptomatically and others become severely ill. Understanding how these divergent phenotypes arise could eventually be leveraged to design therapeutics that prevent severe disease. As a first step toward understanding these different infection states, we used a technology that allowed us to determine how thousands of individual murine lung epithelial cells behaved before and during IAV infection. We found that small subsets of epithelial cells exhibited an antiviral state prior to infection, and similarly, some cells made high levels of inflammatory cytokines during infection. We propose that different ratios of these individual cellular responses may contribute to the broader antiviral state of the lung and may ultimately affect disease severity.

A Virus Is a Community: Diversity within Negative-Sense RNA Virus Populations

Negative-sense RNA virus populations are composed of diverse viral components that interact to form a community and shape the outcome of virus infections. At the genomic level, RNA virus populations consist not only of a homogeneous population of standard viral genomes but also of an extremely large number of genome variants, termed viral quasispecies, and nonstandard viral genomes, which include copy-back viral genomes, deletion viral genomes, mini viral RNAs, and hypermutated RNAs. At the particle level, RNA virus populations are composed of pleomorphic particles, particles missing or having additional genomes, and single particles or particle aggregates. As we continue discovering more about the components of negative-sense RNA virus populations and their crucial functions during virus infection, it will become more important to study RNA virus populations as a whole rather than their individual parts. In this review, we will discuss what is known about the components of negative-sense RNA virus communities, speculate how the components of the virus community interact, and summarize what vaccines and antiviral therapies are being currently developed to target or harness these components.

Evidence for Multiple Subpopulations of Herpesvirus-Latently Infected Cells

Kaposi's sarcoma-associated herpesvirus (KSHV)-associated primary effusion lymphomas (PEL) are traditionally viewed as homogenous regarding viral transcription and lineage of origin, but so far this contention has not been explored at the single-cell level. Single-cell RNA sequencing of latently infected PEL supports the existence of multiple subpopulations even within a single cell line. At most 1% of the cells showed evidence of near-complete lytic transcription. The majority of cells only expressed the canonical viral latent transcripts: those originating from the latency locus, the viral interferon regulatory factor locus, and the viral lncRNA nut-1/Pan/T1.1; however, a significant fraction of cells showed various degrees of more permissive transcription, and some showed no evidence of KSHV transcripts whatsoever. Levels of viral interleukin-6 (IL-6)/K2 mRNA emerged as the most distinguishing feature to subset KSHV-infected PEL. One newly uncovered phenotype is the existence of BCBL-1 cells that readily adhered to fibronectin and that displayed mesenchymal lineage-like characteristics. IMPORTANCE Latency is the defining characteristic of the Herpesviridae and central to the tumorigenesis phenotype of Kaposi's sarcoma-associated herpesvirus (KSHV). KSHV-driven primary effusion lymphomas (PEL) rapidly develop resistance to therapy, suggesting tumor instability and plasticity. At any given time, a fraction of PEL cells spontaneously reactivate KSHV, suggesting transcriptional heterogeneity even within a clonal cell line under optimal growth conditions. This study employed single-cell mRNA sequencing to explore the within-population variability of KSHV transcription and how it relates to host cell transcription. Individual clonal PEL cells exhibited differing patterns of viral transcription. Most cells showed the canonical pattern of KSHV latency (LANA, vCyc, vFLIP, Kaposin, and vIRFs), but a significant fraction evidenced extended viral gene transcription, including of the viral IL-6 homolog, open reading frame K2. This study suggests new targets of intervention for PEL. It establishes a conceptual framework to design KSHV cure studies analogous to those for HIV.

Single-Cell and Bulk RNA-Sequencing Reveal Differences in Monocyte Susceptibility to Influenza A Virus Infection Between Africans and Europeans

There is considerable inter-individual and inter-population variability in response to viruses. The potential of monocytes to elicit type-I interferon responses has attracted attention to their role in viral infections. Here, we use single-cell RNA-sequencing to characterize the role of cellular heterogeneity in human variation of monocyte responses to influenza A virus (IAV) exposure. We show widespread inter-individual variability in the percentage of IAV-infected monocytes. Notably, individuals with high cellular susceptibility to IAV are characterized by a lower activation at basal state of an IRF/STAT-induced transcriptional network, which includes antiviral genes such as IFITM3, MX1 and OAS3. Upon IAV challenge, we find that cells escaping viral infection display increased mRNA expression of type-I interferon stimulated genes and decreased expression of ribosomal genes, relative to both infected cells and those never exposed to IAV. We also uncover a stronger resistance of CD16⁺ monocytes to IAV infection, together with CD16⁺ -specific mRNA expression of IL6 and TNF in response to IAV. Finally, using flow cytometry and bulk RNA-sequencing across 200 individuals of African and European ancestry, we observe a higher number of CD16 ⁺ monocytes and lower susceptibility to IAV infection among monocytes from individuals of African-descent. Based on these data, we hypothesize that higher basal monocyte activation, driven by environmental factors and/or weak-effect genetic variants, underlies the lower cellular susceptibility to IAV infection of individuals of African ancestry relative to those of European ancestry. Further studies are now required to investigate how such cellular differences in IAV susceptibility translate into population differences in clinical outcomes and susceptibility to severe influenza.

Transcriptomics of chicken cecal tonsils and intestine after infection with low pathogenic avian influenza virus H9N2

Influenza viruses cause severe respiratory infections in humans and birds, triggering global health concerns and economic burden. Influenza infection is a dynamic process involving complex biological host responses. The objective of this study was to illustrate global biological processes in ileum and cecal tonsils at early time points after chickens were infected with low pathogenic avian influenza virus (LPAIV) H9N2 through transcriptome analysis. Total RNA isolated from ileum and cecal tonsils of non-infected and infected layers at 12-, 24- and 72-h post-infection (hpi) was used for mRNA sequencing analyses to characterize differentially expressed genes and overrepresented pathways. Statistical analysis highlighted transcriptomic signatures significantly occurring 24 and 72 hpi, but not earlier at 12 hpi. Interferon (IFN)-inducible and IFN-stimulated gene (ISG) expression was increased, followed by continued expression of various heat-shock proteins (HSP), including HSP60, HSP70, HSP90 and HSP110. Some upregulated genes involved in innate antiviral responses included DDX60, MX1, RSAD2 and CMPK2. The ISG15 antiviral mechanism pathway was highly enriched in ileum and cecal tonsils at 24 hpi. Overall, most affected pathways were related to interferon production and the heat-shock response. Research on these candidate genes and pathways is warranted to decipher underlying mechanisms of immunity against LPAIV in chickens.

Understanding the pathogenesis of infectious diseases by single-cell RNA sequencing

Infections are highly orchestrated and dynamic processes, which involve both pathogen and host. Transcriptional profiling at the single-cell level enables the analysis of cell diversity, heterogeneity of the immune response, and detailed molecular mechanisms underlying infectious diseases caused by bacteria, viruses, fungi, and parasites. Herein, we highlight recent remarkable advances in single-cell RNA sequencing (scRNA-seq) technologies and their applications in the investigation of host-pathogen interactions, current challenges and potential prospects for disease treatment are discussed as well. We propose that with the aid of scRNA-seq, the mechanism of infectious diseases will be further revealed thus inspiring the development of novel interventions and therapies.

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.

The transcriptional landscape of Venezuelan equine encephalitis virus (TC-83) infection

Venezuelan Equine Encephalitis Virus (VEEV) is a major biothreat agent that naturally causes outbreaks in humans and horses particularly in tropical areas of the western hemisphere, for which no antiviral therapy is currently available. The host response to VEEV and the cellular factors this alphavirus hijacks to support its effective replication or evade cellular immune responses are largely uncharacterized. We have previously demonstrated tremendous cell-to-cell heterogeneity in viral RNA (vRNA) and cellular transcript levels during flaviviral infection using a novel virus-inclusive single-cell RNA-Seq approach. Here, we used this unbiased, genome-wide approach to simultaneously profile the host transcriptome and vRNA in thousands of single cells during infection of human astrocytes with the live-attenuated vaccine strain of VEEV (TC-83). Host transcription was profoundly suppressed, yet “superproducer cells” with extremely high vRNA abundance emerged during the first viral life cycle and demonstrated an altered transcriptome relative to both uninfected cells and cells with high vRNA abundance harvested at later time points. Additionally, cells with increased structural-to-nonstructural transcript ratio exhibited upregulation of intracellular membrane trafficking genes at later time points. Loss- and gain-of-function experiments confirmed pro- and antiviral activities in both vaccine and virulent VEEV infections among the products of transcripts that positively or negatively correlated with vRNA abundance, respectively. Lastly, comparison with single cell transcriptomic data from other viruses highlighted common and unique pathways perturbed by infection across evolutionary scales. This study provides a high-resolution characterization of the VEEV (TC-83)-host interplay, identifies candidate targets for antivirals, and establishes a comparative single-cell approach to study the evolution of virus-host interactions.

Little is known about the host response to Venezuelan Equine Encephalitis Virus (VEEV) and the cellular factors this alphavirus hijacks to support effective replication or evade cellular immune responses. Monitoring dynamics of host and viral RNA (vRNA) during viral infection at a single-cell level can provide insight into the virus-host interplay at a high resolution. Here, a single-cell RNA sequencing technology that detects host and viral RNA was used to investigate the interactions between TC-83, the vaccine strain of VEEV, and the human host during the course of infection of U-87 MG cells (human astrocytoma). Virus abundance and host transcriptome were heterogeneous across cells from the same culture. Subsets of differentially expressed genes, positively or negatively correlating with vRNA abundance, were identified and subsequently in vitro validated as candidate proviral and antiviral factors, respectively, in TC-83 and/or virulent VEEV infections. In the first replication cycle, “superproducer” cells exhibited rapid increase in vRNA abundance and unique gene expression patterns. At later time points, cells with increased structural-to-nonstructural transcript ratio demonstrated upregulation of intracellular membrane trafficking genes. Lastly, comparing the VEEV dataset with published datasets on other RNA viruses revealed unique and overlapping responses across viral clades. Overall, this study improves the understanding of VEEV-host interactions, reveals candidate targets for antiviral approaches, and establishes a comparative single-cell approach to study the evolution of virus-host interactions.

Viral and host heterogeneity and their effects on the viral life cycle

Traditionally, the viral replication cycle is envisioned as a single, well-defined loop with four major steps: attachment and entry into a target cell, replication of the viral genome, maturation of viral proteins and genome packaging into infectious progeny, and egress and dissemination to the next target cell. However, for many viruses, a growing body of evidence points towards extreme heterogeneity in each of these steps. In this Review, we reassess the major steps of the viral replication cycle by highlighting recent advances that show considerable variability during viral infection. First, we discuss heterogeneity in entry receptors, followed by a discussion on error-prone and low-fidelity polymerases and their impact on viral diversity. Next, we cover the implications of heterogeneity in genome packaging and assembly on virion morphology. Last, we explore alternative egress mechanisms, including tunnelling nanotubes and host microvesicles. In summary, we discuss the implications of viral phenotypic, morphological and genetic heterogeneity on pathogenesis and medicine. This Review highlights common themes and unique features that give nuance to the viral replication cycle.

Cell type- and replication stage-specific influenza virus responses in vivo

Influenza A viruses (IAVs) remain a significant global health burden. Activation of the innate immune response is important for controlling early virus replication and spread. It is unclear how early IAV replication events contribute to immune detection. Additionally, while many cell types in the lung can be infected, it is not known if all cell types contribute equally to establish the antiviral state in the host. Here, we use single-cycle influenza A viruses (scIAVs) to characterize the early immune response to IAV in vitro and in vivo. We found that the magnitude of virus replication contributes to antiviral gene expression within infected cells prior to the induction of a global response. We also developed a scIAV that is only capable of undergoing primary transcription, the earliest stage of virus replication. Using this tool, we uncovered replication stage-specific responses in vitro and in vivo. Using several innate immune receptor knockout cell lines, we identify RIG-I as the predominant antiviral detector of primary virus transcription and amplified replication in vitro. Through a Cre-inducible reporter mouse, we used scIAVs expressing Cre-recombinase to characterize cell type-specific responses in vivo. Individual cell types upregulate unique sets of antiviral genes in response to both primary virus transcription and amplified replication. We also identified antiviral genes that are only upregulated in response to direct infection. Altogether, these data offer insight into the early mechanisms of antiviral gene activation during influenza A infection.

Single cell heterogeneity in influenza A virus gene expression shapes the innate antiviral response to infection

Viral infection outcomes are governed by the complex and dynamic interplay between the infecting virus population and the host response. It is increasingly clear that both viral and host cell populations are highly heterogeneous, but little is known about how this heterogeneity influences infection dynamics or viral pathogenicity. To dissect the interactions between influenza A virus (IAV) and host cell heterogeneity, we examined the combined host and viral transcriptomes of thousands of individual cells, each infected with a single IAV virion. We observed complex patterns of viral gene expression and the existence of multiple distinct host transcriptional responses to infection at the single cell level. We show that human H1N1 and H3N2 strains differ significantly in patterns of both viral and host anti-viral gene transcriptional heterogeneity at the single cell level. Our analyses also reveal that semi-infectious particles that fail to express the viral NS can play a dominant role in triggering the innate anti-viral response to infection. Altogether, these data reveal how patterns of viral population heterogeneity can serve as a major determinant of antiviral gene activation.

COVID-19 Is Distinct From SARS-CoV-2-Negative Community-Acquired Pneumonia

Background: Corona virus disease (COVID-19) is an infectious respiratory disease that has spread rapidly across the world. Many studies have already evaluated the clinical features of COVID-19, but how it compares with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-negative community-acquired pneumonia (SN-CAP) is still unclear. Moreover, COVID-19 mortality is correlated with disease severity, but indicators for severity grading have not been specified. We aimed to analyze the clinical characteristics of COVID-19 in comparison with SN-CAP and find indicators for disease severity in COVID-19. Methods: Patients diagnosed with COVID-19 and SN-CAP were enrolled. Clinical, radiological, and laboratory data were analyzed. Results: The numbers of COVID-19 and SN-CAP patients enrolled were 304 and 138, respectively. The age of the patients was not significantly different between the groups. Compared with SN-CAP, COVID-19 patients had more symptoms of fever and dyspnea; and showed significant difference in blood count results. Computed tomography (CT) imaging of COVID-19 patients showed patchy ground-glass opacities that correlated with disease severity, whereas the CT imaging of SN-CAP patients showed patchy high-density shadows. COVID-19 patients were classified into moderate, severe, and critically severe groups. The severe and critically severe groups had elevated levels of white blood cells (WBC), neutrophils, platelets, C-reaction protein (CRP), lymphocyte ratio (NLR), platelet to lymphocyte ratio (PLR), troponin-I, creatinine, and blood urea nitrogen (BUN). However, they had decreased levels of lymphocytes, lymphocyte ratio, and albumin. Compared with the younger patients, the older COVID-19 individuals had more chronic diseases and significantly elevated levels of WBC, neutrophil, and CRP levels. Conclusion: SN-CAP showed more inflammatory reaction than COVID-19. Old people with chronic diseases are more susceptible to COVID-19 and have a high likelihood of developing severe and critically severe infection. Levels of WBC, lymphocytes, neutrophils, CRP, NLR, PLR, troponin-I, creatinine, and BUN are important indicators for severity grading in COVID-19.

Microfluidics in Single-Cell Virology: Technologies and Applications

Microfluidics has proven to be a powerful tool for probing biology at the single-cell level. However, it is only in the past 5 years that single-cell microfluidics has been used in the field of virology. An array of strategies based on microwells, microvalves, and droplets is now available for tracking viral infection dynamics, identifying cell subpopulations with particular phenotypes, as well as high-throughput screening. The insights into the virus-host interactions gained at the single-cell level are unprecedented and usually inaccessible by population-based experiments. Therefore, single-cell microfluidics, which opens new avenues for mechanism elucidation and development of antiviral therapeutics, would be a valuable tool for the study of viral pathogenesis.

Host-Virus Interaction: How Host Cells Defend against Influenza A Virus Infection

Influenza A viruses (IAVs) are highly contagious pathogens infecting human and numerous animals. The viruses cause millions of infection cases and thousands of deaths every year, thus making IAVs a continual threat to global health. Upon IAV infection, host innate immune system is triggered and activated to restrict virus replication and clear pathogens. Subsequently, host adaptive immunity is involved in specific virus clearance. On the other hand, to achieve a successful infection, IAVs also apply multiple strategies to avoid be detected and eliminated by the host immunity. In the current review, we present a general description on recent work regarding different host cells and molecules facilitating antiviral defenses against IAV infection and how IAVs antagonize host immune responses.

Single-Cell Analysis Uncovers a Vast Diversity in Intracellular Viral Defective Interfering RNA Content Affecting the Large Cell-to-Cell Heterogeneity in Influenza A Virus Replication

Virus replication displays a large cell-to-cell heterogeneity; yet, not all sources of this variability are known. Here, we study the effect of defective interfering (DI) particle (DIP) co-infection on cell-to-cell variability in influenza A virus (IAV) replication. DIPs contain a large internal deletion in one of their eight viral RNAs (vRNA) and are, thus, defective in virus replication. Moreover, they interfere with virus replication. Using single-cell isolation and reverse transcription polymerase chain reaction, we uncovered a large between-cell heterogeneity in the DI vRNA content of infected cells, which was confirmed for DI mRNAs by single-cell RNA sequencing. A high load of intracellular DI vRNAs and DI mRNAs was found in low-productive cells, indicating their contribution to the large cell-to-cell variability in virus release. Furthermore, we show that the magnitude of host cell mRNA expression (some factors may inhibit virus replication), but not the ribosome content, may further affect the strength of single-cell virus replication. Finally, we show that the load of viral mRNAs (facilitating viral protein production) and the DI mRNA content are, independently from one another, connected with single-cell virus production. Together, these insights advance single-cell virology research toward the elucidation of the complex multi-parametric origin of the large cell-to-cell heterogeneity in virus infections.