The defective component of viral populations

A mathematical model simulating the adaptive immune response in various vaccines and vaccination strategies

Vaccination has been widely recognized as an effective measure for preventing infectious diseases. To facilitate quantitative research into the activation of adaptive immune responses in the human body by vaccines, it is important to develop an appropriate mathematical model, which can provide valuable guidance for vaccine development. In this study, we constructed a novel mathematical model to simulate the dynamics of antibody levels following vaccination, based on principles from immunology. Our model offers a concise and accurate representation of the kinetics of antibody response. We conducted a comparative analysis of antibody dynamics within the body after administering several common vaccines, including traditional inactivated vaccines, mRNA vaccines, and future attenuated vaccines based on defective interfering viral particles (DVG). Our findings suggest that booster shots play a crucial role in enhancing Immunoglobulin G (IgG) antibody levels, and we provide a detailed discussion on the advantages and disadvantages of different vaccine types. From a mathematical standpoint, our model proposes four essential approaches to guide vaccine design: enhancing antigenic T-cell immunogenicity, directing the production of high-affinity antibodies, reducing the rate of IgG decay, and lowering the peak level of vaccine antigen-antibody complexes. Our study contributes to the understanding of vaccine design and its application by explaining various phenomena and providing guidance in comprehending the interactions between antibodies and antigens during the immune process.

Chikungunya Virus RNA Secondary Structures Impact Defective Viral Genome Production

Chikungunya virus (CHIKV) is a mosquito-borne RNA virus that poses an emerging threat to humans. In a manner similar to other RNA viruses, CHIKV encodes an error-prone RNA polymerase which, in addition to producing full-length genomes, gives rise to truncated, non-functional genomes, which have been coined defective viral genomes (DVGs). DVGs have been intensively studied in the context of therapy, as they can inhibit viral replication and dissemination in their hosts. In this work, we interrogate the influence of viral RNA secondary structures on the production of CHIKV DVGs. We experimentally map RNA secondary structures of the CHIKV genome using selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), which couples chemical labelling with next-generation sequencing. We correlate the inferred secondary structure with preferred deletion sites of CHIKV DVGs. We document an increased probability of DVG generation with truncations at unpaired nucleotides within the secondary structure. We then generated a CHIKV mutant bearing synonymous changes at the nucleotide level to disrupt the existing RNA secondary structure (CHIKV-D2S). We show that CHIKV-D2S presents altered DVG generation compared to wild-type virus, correlating with the change in RNA secondary structure obtained by SHAPE-MaP. Our work thus demonstrates that RNA secondary structure impacts CHIKV DVG production during replication.

Expanding Possibilities for Foreign Gene Expression by Cucumber Green Mottle Mosaic Virus Genome-Based Bipartite Vector System

Expanding possibilities for foreign gene expression in cucurbits, we present a novel approach utilising a bipartite vector system based on the cucumber green mottle mosaic virus (CGMMV) genome. Traditional full-length CGMMV vectors face limitations such as a restricted cargo capacity and unstable foreign gene expression. To address these challenges, we developed two 'deconstructed' CGMMV genomes, DG-1 and DG-2. DG-1 features a major internal deletion, resulting in the loss of crucial replicase enzyme domains, rendering it incapable of self-replication. However, a staggered infiltration of DG-1 in CGMMV-infected plants enabled successful replication and movement, facilitating gene-silencing experiments. Conversely, DG-2 was engineered to enhance replication rates and provide multiple cloning sites. Although it exhibited higher replication rates, DG-2 remained localised within infiltrated tissue, displaying trans-replication and restricted movement. Notably, DG-2 demonstrated utility in expressing GFP, with a peak expression observed between 6 and 10 days post-infiltration. Overall, our bipartite system represents a significant advancement in functional genomics, offering a robust tool for foreign gene expression in Nicotiana benthamiana.

Accumulation Dynamics of Defective Genomes during Experimental Evolution of Two Betacoronaviruses

Virus-encoded replicases often generate aberrant RNA genomes, known as defective viral genomes (DVGs). When co-infected with a helper virus providing necessary proteins, DVGs can multiply and spread. While DVGs depend on the helper virus for propagation, they can in some cases disrupt infectious virus replication, impact immune responses, and affect viral persistence or evolution. Understanding the dynamics of DVGs alongside standard viral genomes during infection remains unclear. To address this, we conducted a long-term experimental evolution of two betacoronaviruses, the human coronavirus OC43 (HCoV-OC43) and the murine hepatitis virus (MHV), in cell culture at both high and low multiplicities of infection (MOI). We then performed RNA-seq at regular time intervals, reconstructed DVGs, and analyzed their accumulation dynamics. Our findings indicate that DVGs evolved to exhibit greater diversity and abundance, with deletions and insertions being the most common types. Notably, some high MOI deletions showed very limited temporary existence, while others became prevalent over time. We observed differences in DVG abundance between high and low MOI conditions in HCoV-OC43 samples. The size distribution of HCoV-OC43 genomes with deletions differed between high and low MOI passages. In low MOI lineages, short and long DVGs were the most common, with an additional cluster in high MOI lineages which became more prevalent along evolutionary time. MHV also showed variations in DVG size distribution at different MOI conditions, though they were less pronounced compared to HCoV-OC43, suggesting a more random distribution of DVG sizes. We identified hotspot regions for deletions that evolved at a high MOI, primarily within cistrons encoding structural and accessory proteins. In conclusion, our study illustrates the widespread formation of DVGs during betacoronavirus evolution, influenced by MOI and cell- and virus-specific factors.

Expression dynamics of the aplysia abyssovirus

Two recent studies documented the genome of a novel, extremely large (35.9 kb), nidovirus in RNA sequence databases from the marine neural model Aplysia californica. The goal of the present study was to document the distribution and transcriptional dynamics of this virus, Aplysia abyssovirus 1 (AAbV), in maricultured and wild animals. We confirmed previous findings that AAbV RNA is widespread and reaches extraordinary levels in apparently healthy animals. Transmission electron microscopy identified viral replication factories in ciliated gill epithelial cells but not in neurons where viral RNA is most highly expressed. Viral transcripts do not exhibit evidence of discontinuous RNA synthesis as in coronaviruses but are consistent with production of a single leaderless subgenomic RNA, as in the Gill-associated virus of Penaeus monodon. Splicing patterns in chronically infected adults suggested high levels of defective genomes, possibly explaining the lack of obvious disease signs in high viral load animals.

Shared sequence characteristics identified in non-canonical rearrangements of HSV-1 genomes

IMPORTANCE

Mutations and genetic rearrangements are the primary driving forces of evolution. Viruses provide valuable model systems for investigating these mechanisms due to their rapid evolutionary rates and vast genetic variability. To investigate genetic rearrangements in the double-stranded DNA genome of herpes simplex virus type 1, the viral population was serially passaged in various cell types. The serial passaging led to formation of defective genomes, resulted from cell-specific non-canonical rearrangements (NCRs). Interestingly, we discovered shared sequence characteristics underlying the formation of these NCRs across all cell types. Moreover, most NCRs identified in clinical samples shared these characteristics. Based on our findings, we propose a model elucidating the formation of NCRs during viral replication within the nucleus of eukaryotic cells.

Metagenome-derived virus-microbe ratios across ecosystems

It is generally assumed that viruses outnumber cells on Earth by at least tenfold. Virus-to-microbe ratios (VMR) are largely based on counts of fluorescently labelled virus-like particles. However, these exclude intracellular viruses and potentially include false positives (DNA-containing vesicles, gene-transfer agents, unspecifically stained inert particles). Here, we develop a metagenome-based VMR estimate (mVRM) that accounts for DNA viruses across all stages of their replication cycles (virion, intracellular lytic and lysogenic) by using normalised RPKM (reads per kilobase of gene sequence per million of mapped metagenome reads) counts of the major capsid protein (MCP) genes and cellular universal single-copy genes (USCGs) as proxies for virus and cell counts, respectively. After benchmarking this strategy using mock metagenomes with increasing VMR, we inferred mVMR across different biomes. To properly estimate mVMR in aquatic ecosystems, we generated metagenomes from co-occurring cellular and viral fractions (>50 kDa-200 µm size-range) in freshwater, seawater and solar saltern ponds (10 metagenomes, 2 control metaviromes). Viruses outnumbered cells in freshwater by ~13 fold and in plankton from marine and saline waters by ~2-4 fold. However, across an additional set of 121 diverse non-aquatic metagenomes including microbial mats, microbialites, soils, freshwater and marine sediments and metazoan-associated microbiomes, viruses, on average, outnumbered cells by barely two-fold. Although viruses likely are the most diverse biological entities on Earth, their global numbers might be closer to those of cells than previously estimated.

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.

Superinfection exclusion creates spatially distinct influenza virus populations

Interactions between viruses during coinfections can influence viral fitness and population diversity, as seen in the generation of reassortant pandemic influenza A virus (IAV) strains. However, opportunities for interactions between closely related viruses are limited by a process known as superinfection exclusion (SIE), which blocks coinfection shortly after primary infection. Using IAVs, we asked whether SIE, an effect which occurs at the level of individual cells, could limit interactions between populations of viruses as they spread across multiple cells within a host. To address this, we first measured the kinetics of SIE in individual cells by infecting them sequentially with 2 isogenic IAVs, each encoding a different fluorophore. By varying the interval between addition of the 2 IAVs, we showed that early in infection SIE does not prevent coinfection, but that after this initial lag phase the potential for coinfection decreases exponentially. We then asked how the kinetics of SIE onset controlled coinfections as IAVs spread asynchronously across monolayers of cells. We observed that viruses at individual coinfected foci continued to coinfect cells as they spread, because all new infections were of cells that had not yet established SIE. In contrast, viruses spreading towards each other from separately infected foci could only establish minimal regions of coinfection before reaching cells where coinfection was blocked. This created a pattern of separate foci of infection, which was recapitulated in the lungs of infected mice, and which is likely to be applicable to many other viruses that induce SIE. We conclude that the kinetics of SIE onset segregate spreading viral infections into discrete regions, within which interactions between virus populations can occur freely, and between which they are blocked.

Genetic complementation fosters evolvability in complex fitness landscapes

The ability of natural selection to optimize traits depends on the topology of the genotype-fitness map (fitness landscape). Epistatic interactions produce rugged fitness landscapes, where adaptation is constrained by the presence of low-fitness intermediates. Here, we used simulations to explore how evolvability in rugged fitness landscapes is influenced by genetic complementation, a process whereby different sequence variants mutually compensate for their deleterious mutations. We designed our model inspired by viral populations, in which genetic variants are known to interact frequently through coinfection. Our simulations indicate that genetic complementation enables a more efficient exploration of rugged fitness landscapes. Although this benefit may be undermined by genetic parasites, its overall effect on evolvability remains positive in populations that exhibit strong relatedness between interacting sequences. Similar processes could operate in contexts other than viral coinfection, such as in the evolution of ploidy.

Arm race between Rift Valley fever virus and host

Rift Valley fever (RVF) is a zoonotic disease caused by Rift Valley fever virus (RVFV), an emerging arbovirus within the Phenuiviridae family of Bunyavirales that has potential to cause severe diseases in both humans and livestock. It increases the incidence of abortion or foetal malformation in ruminants and leads to clinical manifestations like encephalitis or haemorrhagic fever in humans. Upon virus invasion, the innate immune system from the cell or the organism is activated to produce interferon (IFN) and prevent virus proliferation. Meanwhile, RVFV initiates countermeasures to limit antiviral responses at transcriptional and protein levels. RVFV nonstructural proteins (NSs) are the key virulent factors that not only perform immune evasion but also impact the cell replication cycle and has cytopathic effects. In this review, we summarize the innate immunity host cells employ depending on IFN signal transduction pathways, as well as the immune evasion mechanisms developed by RVFV primarily with the inhibitory activity of NSs protein. Clarifying the arms race between host innate immunity and RVFV immune evasion provides new avenues for drug target screening and offers possible solutions to current and future epidemics.

Treatment of Sindbis Virus-Infected Neurons with Antibody to E2 Alters Synthesis of Complete and nsP1-Expressing Defective Viral RNAs

Alphaviruses are positive-sense RNA viruses that are important causes of viral encephalomyelitis. Sindbis virus (SINV), the prototype alphavirus, preferentially infects neurons in mice and is a model system for studying mechanisms of viral clearance from the nervous system. Antibody specific to the SINV E2 glycoprotein plays an important role in SINV clearance, and this effect is reproduced in cultures of infected mature neurons. To determine how anti-E2 antibody affects SINV RNA synthesis, Oxford Nanopore Technologies direct long-read RNA sequencing was used to sequence viral RNAs following antibody treatment of infected neurons. Differentiated AP-7 rat olfactory neuronal cells, an in vitro model for mature neurons, were infected with SINV and treated with anti-E2 antibody. Whole-cell RNA lysates were collected for sequencing of poly(A)-selected RNA 24, 48, and 72 h after infection. Three primary species of viral RNA were produced: genomic, subgenomic, and defective viral genomes (DVGs) encoding the RNA capping protein nsP1. Antibody treatment resulted in overall lower production of SINV RNA, decreased synthesis of subgenomic RNA relative to genomic RNA, and suppressed production of the nsP1 DVG. The nsP1 DVG was packaged into virus particles and could be translated. Because antibody-treated cells released a higher proportion of virions with noncapped genomes and transient transfection to express the nsP1 DVG improved viral RNA capping in antibody-treated cells, we postulate that one mechanism by which antibody inhibits SINV replication in neurons is to suppress DVG synthesis and thus decrease production of infectious virions containing capped genomes. IMPORTANCE Alphaviruses are important causes of viral encephalomyelitis without approved treatments or vaccines. Antibody to the Sindbis virus (SINV) E2 glycoprotein is required for immune-mediated noncytolytic virus clearance from neurons. We used direct RNA nanopore sequencing to evaluate how anti-E2 antibody affects SINV replication at the RNA level. Antibody altered the viral RNAs produced by decreasing the proportion of subgenomic relative to genomic RNA and suppressing production of a previously unrecognized defective viral genome (DVG) encoding nsP1, the viral RNA capping enzyme. Antibody-treated neurons released a lower proportion of SINV particles with capped genomes necessary for translation and infection. Decreased nsP1 DVG production in antibody-treated neurons led to lower expression of nsP1 protein, decreased genome capping efficiency, and release of fewer infectious virus particles. Capping was increased with exogenous expression of the nsP1 DVG. These studies identify a novel alphavirus DVG function and new mechanism for antibody-mediated control of virus replication.

Depletion of CD4 and CD8 T Cells Reduces Acute Disease and Is Not Associated with Hearing Loss in ML29-Infected STAT1-/- Mice

Lassa virus (LASV) is a zoonotic virus endemic to western Africa that can cause a potentially lethal and hemorrhagic disease, Lassa fever (LF). Survivors suffer a myriad of sequelae, most notably sudden onset sensorineural hearing loss (SNHL), the mechanism of which remains unclear. Unfortunately, studies aiming to identify the mechanism of these sequelae are limited due to the biosafety level 4 (BSL4) requirements of LASV itself. ML29, a reassortant virus proposed as an experimental vaccine candidate against LASV, is potentially an ideal surrogate model of LF in STAT1-/- mice due to similar phenotype in these animals. We intended to better characterize ML29 pathogenesis and potential sequelae in this animal model. Our results indicate that while both CD4 and CD8 T cells are responsible for acute disease in ML29 infection, ML29 induces significant hearing loss in a mechanism independent of either CD4 or CD8 T cells. We believe that this model could provide valuable information for viral-associated hearing loss in general.

The simple emergence of complex molecular function

At odds with a traditional view of molecular evolution that seeks a descent-with-modification relationship between functional sequences, new functions can emerge de novo with relative ease. At early times of molecular evolution, random polymers could have sufficed for the appearance of incipient chemical activity, while the cellular environment harbours a myriad of proto-functional molecules. The emergence of function is facilitated by several mechanisms intrinsic to molecular organization, such as redundant mapping of sequences into structures, phenotypic plasticity, modularity or cooperative associations between genomic sequences. It is the availability of niches in the molecular ecology that filters new potentially functional proposals. New phenotypes and subsequent levels of molecular complexity could be attained through combinatorial explorations of currently available molecular variants. Natural selection does the rest. This article is part of the theme issue 'Emergent phenomena in complex physical and socio-technical systems: from cells to societies'.

Droplet Microfluidics for High-Resolution Virology

Microfluidics has enabled a new era of cellular and molecular assays due to the small length scales, parallelization, and the modularity of various analysis and actuation functions. Droplet microfluidics, in particular, has been instrumental in providing new tools for biology with its ability to quickly and reproducibly generate drops that act as individual reactors. A notable beneficiary of this technology has been single-cell RNA sequencing, which has revealed new heterogeneities and interactions for the fundamental unit of life. However, viruses far surpass the diversity of cellular life, affect the dynamics of all ecosystems, and are a chronic source of global health crises. Despite their impact on the world, high-throughput and high-resolution viral profiling has been difficult, with conventional methods being limited to population-level averaging, large sample volumes, and few cultivable hosts. Consequently, most viruses have not been identified and studied. Droplet microfluidics holds the potential to address many of these limitations and offers new levels of sensitivity and throughput for virology. This Feature highlights recent efforts that have applied droplet microfluidics to the detection and study of viruses, including for diagnostics, virus-host interactions, and cell-independent virus assays. In combination with traditional virology methods, droplet microfluidics should prove a potent tool toward achieving a better understanding of the most abundant biological species on Earth.

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.

Therapeutic interfering particles exploiting viral replication and assembly mechanisms show promising performance: a modelling study

Defective interfering particles arise spontaneously during a viral infection as mutants lacking essential parts of the viral genome. Their ability to replicate in the presence of the wild-type (WT) virus (at the expense of viable viral particles) is mimicked and exploited by therapeutic interfering particles. We propose a strategy for the design of therapeutic interfering RNAs (tiRNAs) against positive-sense single-stranded RNA viruses that assemble via packaging signal-mediated assembly. These tiRNAs contain both an optimised version of the virus assembly manual that is encoded by multiple dispersed RNA packaging signals and a replication signal for viral polymerase, but lack any protein coding information. We use an intracellular model for hepatitis C viral (HCV) infection that captures key aspects of the competition dynamics between tiRNAs and viral genomes for virally produced capsid protein and polymerase. We show that only a small increase in the assembly and replication efficiency of the tiRNAs compared with WT virus is required in order to achieve a treatment efficacy greater than 99%. This demonstrates that the proposed tiRNA design could be a promising treatment option for RNA viral infections.

The evolution of cheating in viruses

The success of many viruses depends upon cooperative interactions between viral genomes. However, whenever cooperation occurs, there is the potential for 'cheats' to exploit that cooperation. We suggest that: (1) the biology of viruses makes viral cooperation particularly susceptible to cheating; (2) cheats are common across a wide range of viruses, including viral entities that are already well studied, such as defective interfering genomes, and satellite viruses. Consequently, the evolutionary theory of cheating could help us understand and manipulate viral dynamics, while viruses also offer new opportunities to study the evolution of cheating.

Organ-specific genome diversity of replication-competent SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is not always confined to the respiratory system, as it impacts people on a broad clinical spectrum from asymptomatic to severe systemic manifestations resulting in death. Further, accumulation of intra-host single nucleotide variants during prolonged SARS-CoV-2 infection may lead to emergence of variants of concern (VOCs). Still, information on virus infectivity and intra-host evolution across organs is sparse. We report a detailed virological analysis of thirteen postmortem coronavirus disease 2019 (COVID-19) cases that provides proof of viremia and presence of replication-competent SARS-CoV-2 in extrapulmonary organs of immunocompromised patients, including heart, kidney, liver, and spleen (NCT04366882). In parallel, we identify organ-specific SARS-CoV-2 genome diversity and mutations of concern N501Y, T1027I, and Y453F, while the patient had died long before reported emergence of VOCs. These mutations appear in multiple organs and replicate in Vero E6 cells, highlighting their infectivity. Finally, we show two stages of fatal disease evolution based on disease duration and viral loads in lungs and plasma. Our results provide insights about the pathogenesis and intra-host evolution of SARS-CoV-2 and show that COVID-19 treatment and hygiene measures need to be tailored to specific needs of immunocompromised patients, even when respiratory symptoms cease.

The Curious Strategy of Multipartite Viruses

Multipartite virus genomes are composed of several segments, each packaged in a distinct viral particle. Although this puzzling genome architecture is found in ∼17% of known viral species, its distribution among hosts or among distinct types of genome-composing nucleic acid remains poorly understood. No convincing advantage of multipartitism has been identified, yet the maintenance of genomic integrity appears problematic. Here we review recent studies shedding light on these issues. Multipartite viruses rapidly modify the copy number of each segment/gene from one host species to another, a putative benefit if host switches are common. One multipartite virus functions in a multicellular way: The segments do not all need to be present in the same cell and can functionally complement across cells, maintaining genome integrity within hosts. The genomic integrity maintenance during host-to-host transmission needs further elucidation. These features challenge several virology foundations and could apply to other multicomponent viral systems.

The Social Life of Viruses

Despite their simplicity, viruses exhibit certain types of social interactions. Situations in which a given virus achieves higher fitness in combination with other members of the viral population have been described at the level of transmission, replication, suppression of host immune responses, and host killing, enabling the evolution of viral cooperation. Although cellular coinfection with multiple viral particles is the typical playground for these interactions, cooperation between viruses infecting different cells is also established through cellular and viral-encoded communication systems. In general, the stability of cooperation is compromised by cheater genotypes, as best exemplified by defective interfering particles. As predicted by social evolution theory, cheater invasion can be avoided when cooperators interact preferentially with other cooperators, a situation that is promoted in spatially structured populations. Processes such as transmission bottlenecks, organ compartmentalization, localized spread of infection foci, superinfection exclusion, and even discrete intracellular replication centers promote multilevel spatial structuring in viruses. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Reovirus RNA recombination is sequence directed and generates internally deleted defective genome segments during passage

For viruses with segmented genomes, genetic diversity is generated by genetic drift, reassortment, and recombination. Recombination produces RNA populations distinct from full-length gene segments and can influence viral population dynamics, persistence, and host immune responses. Viruses in the Reoviridae family, including rotavirus and mammalian orthoreovirus (reovirus), have been reported to package segments containing rearrangements or internal deletions. Rotaviruses with RNA segments containing rearrangements have been isolated from immunocompromised and immunocompetent children and in vitro following serial passage at relatively high multiplicity. Reoviruses that package small, defective RNA segments have established chronic infections in cells and in mice. However, the mechanism and extent of Reoviridae RNA recombination are undefined. Towards filling this gap in knowledge, we determined the titers and RNA segment profiles for reovirus and rotavirus following serial passage in cultured cells. The viruses exhibited occasional titer reductions characteristic of interference. Reovirus strains frequently accumulated segments that retained 5' and 3' terminal sequences and featured large internal deletions, while similarly fragmented segments were rarely detected in rotavirus populations. Using next-generation RNA-sequencing to analyze RNA molecules packaged in purified reovirus particles, we identified distinct recombination sites within individual viral genome segments. Recombination junctions were frequently but not always characterized by short direct sequence repeats upstream and downstream that spanned junction sites. Taken together, these findings suggest that reovirus accumulates defective gene segments featuring internal deletions during passage and undergoes sequence-directed recombination at distinct sites.IMPORTANCE Viruses in the Reoviridae family include important pathogens of humans and other animals and have segmented RNA genomes. Recombination in RNA virus populations can facilitate novel host exploration and increased disease severity. The extent, patterns, and mechanisms of Reoviridae recombination and the functions and effects of recombined RNA products are poorly understood. Here, we provide evidence that mammalian orthoreovirus regularly synthesizes RNA recombination products that retain terminal sequences but contain internal deletions, while rotavirus rarely synthesizes such products. Recombination occurs more frequently at specific sites in the mammalian orthoreovirus genome, and short regions of identical sequence are often detected at junction sites. These findings suggest that mammalian orthoreovirus recombination events are directed in part by RNA sequences. An improved understanding of recombined viral RNA synthesis may enhance our capacity to engineer improved vaccines and virotherapies in the future.

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.

Genomic diversity in a population of Spodoptera frugiperda nucleopolyhedrovirus

Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV) represents a strong candidate to develop environmental-friendly pesticides against the fall armyworm (Spodoptera frugiperda), a widespread pest that poses a severe threat to different crops around the world. To date, SfMNPV genomic diversity of different isolates has been mainly studied by means of restriction pattern analyses and by sequencing of the egt region. Here, the genomic diversity present inside an isolate of SfMNPV was explored using high-throughput sequencing for the first time. We identified 704 intrahost single nucleotide variants, from which 184 are nonsynonymous mutations distributed among 82 different coding sequences. We detected several structural variants affecting SfMNPV genome, including two previously reported deletions inside the egt region. A comparative analysis between polymorphisms present in different SfMNPV isolates and our intraisolate diversity data suggests that coding regions with higher genetic diversity are associated with oral infectivity or unknown functions. In this context, through molecular evolution studies we provide evidence of diversifying selection acting on sf29, a putative collagenase which could contribute to the oral infectivity of SfMNPV. Overall, our results contribute to deepen our understanding of the coevolution between SfMNPV and the fall armyworm and will be useful to improve the applicability of this virus as a biological control agent.

Defective viral genomes from chikungunya virus are broad-spectrum antivirals and prevent virus dissemination in mosquitoes

Defective viral genomes (DVGs) are truncated and/or rearranged viral genomes produced during virus replication. Described in many RNA virus families, some of them have interfering activity on their parental virus and/or strong immunostimulatory potential, and are being considered in antiviral approaches. Chikungunya virus (CHIKV) is an alphavirus transmitted by Aedes spp. that infected millions of humans in the last 15 years. Here, we describe the DVGs arising during CHIKV infection in vitro in mammalian and mosquito cells, and in vivo in experimentally infected Aedes aegypti mosquitoes. We combined experimental and computational approaches to select DVG candidates most likely to have inhibitory activity and showed that, indeed, they strongly interfere with CHIKV replication both in mammalian and mosquito cells. We further demonstrated that some DVGs present broad-spectrum activity, inhibiting several CHIKV strains and other alphaviruses. Finally, we showed that pre-treating Aedes aegypti with DVGs prevented viral dissemination in vivo.

Lassa Virus Vaccine Candidate ML29 Generates Truncated Viral RNAs Which Contribute to Interfering Activity and Attenuation

Defective interfering particles (DIPs) are naturally occurring products during virus replication in infected cells. DIPs contain defective viral genomes (DVGs) and interfere with replication and propagation of their corresponding standard viral genomes by competing for viral and cellular resources, as well as promoting innate immune antiviral responses. Consequently, for many different viruses, including mammarenaviruses, DIPs play key roles in the outcome of infection. Due to their ability to broadly interfere with viral replication, DIPs are attractive tools for the development of a new generation of biologics to target genetically diverse and rapidly evolving viruses. Here, we provide evidence that in cells infected with the Lassa fever (LF) vaccine candidate ML29, a reassortant that carries the nucleoprotein (NP) and glycoprotein (GP) dominant antigens of the pathogenic Lassa virus (LASV) together with the L polymerase and Z matrix protein of the non-pathogenic genetically related Mopeia virus (MOPV), L-derived truncated RNA species are readily detected following infection at low multiplicity of infection (MOI) or in persistently-infected cells originally infected at high MOI. In the present study, we show that expression of green fluorescent protein (GFP) driven by a tri-segmented form of the mammarenavirus lymphocytic choriomeningitis virus (r3LCMV-GFP/GFP) was strongly inhibited in ML29-persistently infected cells, and that the magnitude of GFP suppression was dependent on the passage history of the ML29-persistently infected cells. In addition, we found that DIP-enriched ML29 was highly attenuated in immunocompetent CBA/J mice and in Hartley guinea pigs. Likewise, STAT-1^-/- mice, a validated small animal model for human LF associated hearing loss sequelae, infected with DIP-enriched ML29 did not exhibit any hearing abnormalities throughout the observation period (62 days).

Experimental Evolution Reveals a Genetic Basis for Membrane-Associated Virus Release

Many animal viruses replicate and are released from cells in close association to membranes. However, whether this is a passive process or is controlled by the virus remains poorly understood. Importantly, the genetic basis and evolvability of membrane-associated viral shedding have not been investigated. To address this, we performed a directed evolution experiment using coxsackievirus B3, a model enterovirus, in which we repeatedly selected the free-virion or the fast-sedimenting membrane-associated viral subpopulations. The virus responded to this selection regime by reproducibly fixing a series of mutations that altered the extent of membrane-associated viral shedding, as revealed by full-genome ultra-deep sequencing. Specifically, using site-directed mutagenesis, we showed that substitution N63H in the viral capsid protein VP3 reduced the ratio of membrane-associated to free viral particles by 2 orders of magnitude. These findings open new avenues for understanding the mechanisms and implications of membrane-associated viral transmission.

RanDeL-Seq: a High-Throughput Method to Map Viral cis- and trans-Acting Elements

Recent studies have renewed interest in developing novel antiviral therapeutics and vaccines based on defective interfering particles (DIPs)—a subset of viral deletion mutants that conditionally replicate. Identifying and engineering DIPs require that viral cis- and trans-acting elements be accurately mapped.

It has long been known that noncoding genomic regions can be obligate cis elements acted upon in trans by gene products. In viruses, cis elements regulate gene expression, encapsidation, and other maturation processes, but mapping these elements relies on targeted iterative deletion or laborious prospecting for rare spontaneously occurring mutants. Here, we introduce a method to comprehensively map viral cis and trans elements at single-nucleotide resolution by high-throughput random deletion. Variable-size deletions are randomly generated by transposon integration, excision, and exonuclease chewback and then barcoded for tracking via sequencing (i.e., random deletion library sequencing [RanDeL-seq]). Using RanDeL-seq, we generated and screened >23,000 HIV-1 variants to generate a single-base resolution map of HIV-1’s cis and trans elements. The resulting landscape recapitulated HIV-1’s known cis-acting elements (i.e., long terminal repeat [LTR], Ψ, and Rev response element [RRE]) and, surprisingly, indicated that HIV-1’s central DNA flap (i.e., central polypurine tract [cPPT] to central termination sequence [CTS]) is as critical as the LTR, Ψ, and RRE for long-term passage. Strikingly, RanDeL-seq identified a previously unreported ∼300-bp region downstream of RRE extending to splice acceptor 7 that is equally critical for sustained viral passage. RanDeL-seq was also used to construct and screen a library of >90,000 variants of Zika virus (ZIKV). Unexpectedly, RanDeL-seq indicated that ZIKV’s cis-acting regions are larger than the untranscribed (UTR) termini, encompassing a large fraction of the nonstructural genes. Collectively, RanDeL-seq provides a versatile framework for generating viral deletion mutants, enabling discovery of replication mechanisms and development of novel antiviral therapeutics, particularly for emerging viral infections.

Mixtures of Insect-Pathogenic Viruses in a Single Virion: towards the Development of Custom-Designed Insecticides

Alphabaculoviruses (Baculoviridae) are pathogenic DNA viruses of Lepidoptera that have applications as the basis for biological insecticides and expression vectors in biotechnological processes. These viruses have a characteristic physical structure that facilitates the transmission of groups of genomes. We demonstrate that coinfection of a susceptible insect by two different alphabaculovirus species results in the production of mixed-virus occlusion bodies containing the parental viruses. This occurred between closely related and phylogenetically more distant alphabaculoviruses. Approximately half the virions present in proteinaceous viral occlusion bodies produced following coinfection of insects with a mixture of two alphabaculoviruses contained both viruses, indicating that the viruses coinfected and replicated in a single cell and were coenveloped within the same virion. This observation was confirmed by endpoint dilution assay. Moreover, both viruses persisted in the mixed-virus population by coinfection of insects during several rounds of insect-to-insect transmission. Coinfection by viruses that differed in genome size had unexpected results on the length of viral nucleocapsids, which differed from those of both parental viruses. These results have unique implications for the development of alphabaculoviruses as biological control agents of insect pests.IMPORTANCE Alphabaculoviruses are used as biological insecticides and expression vectors in biotechnology and medical applications. We demonstrate that in caterpillars infected with particular mixtures of viruses, the genomes of different baculovirus species can be enveloped together within individual virions and occluded within proteinaceous occlusion bodies. This results in the transmission of mixed-virus populations to the caterpillar stages of moth species. Once established, mixed-virus populations persist by coinfection of insect cells during several rounds of insect-to-insect transmission. Mixed-virus production technology opens the way to the development of custom-designed insecticides for control of different combinations of caterpillar pest species.

Quantification of Ebola virus replication kinetics in vitro

Mathematical modelling has successfully been used to provide quantitative descriptions of many viral infections, but for the Ebola virus, which requires biosafety level 4 facilities for experimentation, modelling can play a crucial role. Ebola virus modelling efforts have primarily focused on in vivo virus kinetics, e.g., in animal models, to aid the development of antivirals and vaccines. But, thus far, these studies have not yielded a detailed specification of the infection cycle, which could provide a foundational description of the virus kinetics and thus a deeper understanding of their clinical manifestation. Here, we obtain a diverse experimental data set of the Ebola virus infection in vitro, and then make use of Bayesian inference methods to fully identify parameters in a mathematical model of the infection. Our results provide insights into the distribution of time an infected cell spends in the eclipse phase (the period between infection and the start of virus production), as well as the rate at which infectious virions lose infectivity. We suggest how these results can be used in future models to describe co-infection with defective interfering particles, which are an emerging alternative therapeutic.

The Role of Temperature in Shaping Mosquito-Borne Viruses Transmission

Mosquito-borne diseases having the greatest impact on human health are typically prevalent in the tropical belt of the world. However, these diseases are conquering temperate regions, raising the question of the role of temperature on their dynamics and expansion. Temperature is one of the most significant abiotic factors affecting, in many ways, insect vectors and the pathogens they transmit. Here, we debate the veracity of this claim by synthesizing current knowledge on the effects of temperature on arboviruses and their vectors, as well as the outcome of their interactions.

Low pathogenic avian influenza virus isolates with different levels of defective genome segments vary in pathogenicity and transmission efficiency

Defective interfering particles (DIPs) of influenza virus are generated through incorporation of highly truncated forms of genome segments, mostly those coding polymerase complex proteins (PB2, PB1, PA). Such particles are able to replicate only in the presence of a virus with the complete genome, thus DIPs may alter the infection outcome by suppressing production of standard virus particles, but also by stimulating the immune response. In the present study we compared the clinical outcome, mortality and transmission in chickens and turkeys infected with the same infectious doses of H7N7 low pathogenic avian influenza virus containing different levels of defective gene segments (95/95(DVG-high) and 95/95(DVG-low)). No clinical signs, mortality or transmission were noted in SPF chickens inoculated with neither virus stock. Turkeys infected with 95/95(DVG-high) showed only slight clinical signs with no mortality, and the virus was transmitted only to birds in direct contact. In contrast, more severe disease, mortality and transmission to direct and indirect contact birds was observed in turkeys infected with 95/95(DVG-low). Apathy, lower water and food intake, respiratory system disorders and a total mortality of 60% were noted. Shedding patterns in contact turkeys indicated more efficient within- and between-host spread of the virus than in 95/95(DVG-high) group. Sequencing of virus genomes showed no mutations that could account for the observed differences in pathogenicity. The results suggest that the abundance of DIPs in the inoculum was the factor responsible for the mild course of infection and disrupted virus transmission.

Competition between social cheater viruses is driven by mechanistically different cheating strategies

Cheater viruses, also known as defective interfering viruses, cannot replicate on their own yet replicate faster than the wild type upon coinfection. While there is growing interest in using cheaters as antiviral therapeutics, the mechanisms underlying cheating have been rarely explored. During experimental evolution of MS2 phage, we observed the parallel emergence of two independent cheater mutants. The first, a point deletion mutant, lacked polymerase activity but was advantageous in viral packaging. The second synonymous mutant cheater displayed a completely different cheating mechanism, involving an altered RNA structure. Continued evolution revealed the demise of the deletion cheater and rise of the synonymous cheater. A mathematical model inferred that while a single cheater is expected to reach an equilibrium with the wild type, cheater demise arises from antagonistic interactions between coinfecting cheaters. These findings highlight layers of parasitism: viruses parasitizing cells, cheaters parasitizing intact viruses, and cheaters may parasitize other cheaters.

Immunopathology and immune homeostasis during viral infection in insects

Organisms clear infections by mounting an immune response that is normally turned off once the pathogens have been cleared. However, sometimes this immune response is not properly or timely arrested, resulting in the host damaging itself. This immune dysregulation may be referred to as immunopathology. While our knowledge of immune and metabolic pathways in insects, particularly in response to viral infections, is growing, little is known about the mechanisms that regulate this immune response and hence little is known about immunopathology in this important and diverse group of organisms. In this chapter we focus both on documenting the molecular mechanisms described involved in restoring immune homeostasis in insects after viral infections and on identifying potential mechanisms for future investigation. We argue that learning about the immunopathological consequences of an improperly regulated immune response in insects will benefit both insect and human health.

Defective Interfering Particles of Negative-Strand RNA Viruses

Viral defective interfering particles (DIPs) were intensely studied several decades ago but research waned leaving open many critical questions. New technologies and other advances led to a resurgence in DIP studies for negative-strand RNA viruses. While DIPs have long been recognized, their exact contribution to the outcome of acute or persistent viral infections has remained elusive. Recent studies have identified defective viral genomes (DVGs) in human infections, including respiratory syncytial virus and influenza, and growing evidence indicates that DVGs influence disease severity and may contribute to viral persistence. Further, several studies have advanced our understanding of key viral and host factors that regulate DIP formation and activity. Here we review these discoveries and highlight key questions moving forward.

Interplay between SARS-CoV-2 and the type I interferon response

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic. An unbalanced immune response, characterized by a weak production of type I interferons (IFN-Is) and an exacerbated release of proinflammatory cytokines, contributes to the severe forms of the disease. SARS-CoV-2 is genetically related to SARS-CoV and Middle East respiratory syndrome-related coronavirus (MERS-CoV), which caused outbreaks in 2003 and 2013, respectively. Although IFN treatment gave some encouraging results against SARS-CoV and MERS-CoV in animal models, its potential as a therapeutic against COVID-19 awaits validation. Here, we describe our current knowledge of the complex interplay between SARS-CoV-2 infection and the IFN system, highlighting some of the gaps that need to be filled for a better understanding of the underlying molecular mechanisms. In addition to the conserved IFN evasion strategies that are likely shared with SARS-CoV and MERS-CoV, novel counteraction mechanisms are being discovered in SARS-CoV-2-infected cells. Since the last coronavirus epidemic, we have made considerable progress in understanding the IFN-I response, including its spatiotemporal regulation and the prominent role of plasmacytoid dendritic cells (pDCs), which are the main IFN-I-producing cells. While awaiting the results of the many clinical trials that are evaluating the efficacy of IFN-I alone or in combination with antiviral molecules, we discuss the potential benefits of a well-timed IFN-I treatment and propose strategies to boost pDC-mediated IFN responses during the early stages of viral infection.

Social Bacteriophages

Despite their simplicity, viruses can display social-like interactions such as cooperation, communication, and cheating. Focusing on bacteriophages, here we review features including viral product sharing, cooperative evasion of antiviral defenses, prudent host exploitation, superinfection exclusion, and inter-phage peptide-mediated signaling. We argue that, in order to achieve a better understanding of these processes, their mechanisms of action need to be considered in the context of social evolution theory, paying special attention to key population-level factors such as genetic relatedness and spatial structure.

A new implication of quasispecies dynamics: Broad virus diversification in absence of external perturbations

RNA genetic elements include many important animal and plant pathogens. They share high mutability, a trait that has multiple implications for the interactions with their host organisms. Here we review evidence of a new adaptive feature of RNA viruses that we term "broadly diversifying selection". It constitutes a new type of positive selection without participation of any external selective agent, and which is built upon a progressive increase of the number of different genomes that dominate the population. The evidence was provided by analyses of mutant spectrum composition of two important viral pathogens, foot-and-mouth disease virus (FMDV) and hepatitis C virus (HCV) after prolonged replication in their respective cell culture environment. Despite being fueled by mutations that arise randomly and in absence of an external guiding selective force, this type of selection prepares the viral population for a response to selective forces still to occur. Since current evidence suggests that broadly diversifying selection is favored by elevated mutation rates and population sizes, it may constitute a more general behavior, relevant also to the adaptive dynamics of microbial populations and cancer cells.

Broad and Dynamic Diversification of Infectious Hepatitis C Virus in a Cell Culture Environment

Previous studies documented that long-term hepatitis C virus (HCV) replication in human hepatoma Huh-7.5 cells resulted in viral fitness gain, expansion of the mutant spectrum, and several phenotypic alterations. In the present work, we show that mutational waves (changes in frequency of individual mutations) occurred continuously and became more prominent as the virus gained fitness. They were accompanied by an increasing proportion of heterogeneous genomic sites that affected 1 position in the initial HCV population and 19 and 69 positions at passages 100 and 200, respectively. Analysis of biological clones of HCV showed that these dynamic events affected infectious genomes, since part of the fluctuating mutations became incorporated into viable genomes. While 17 mutations were scored in 3 biological clones isolated from the initial population, the number reached 72 in 3 biological clones from the population at passage 200. Biological clones differed in their responses to antiviral inhibitors, indicating a phenotypic impact of viral dynamics. Thus, HCV adaptation to a specific constant environment (cell culture without external influences) broadens the mutant repertoire and does not focus the population toward a limited number of dominant genomes. A retrospective examination of mutant spectra of foot-and-mouth disease virus passaged in cell cultures suggests a parallel behavior here described for HCV. We propose that virus diversification in a constant environment has its basis in the availability of multiple alternative mutational pathways for fitness gain. This mechanism of broad diversification should also apply to other replicative systems characterized by high mutation rates and large population sizes.IMPORTANCE The study shows that extensive replication of an RNA virus in a constant biological environment does not limit exploration of sequence space and adaptive options. There was no convergence toward a restricted set of adapted genomes. Mutational waves and mutant spectrum broadening affected infectious genomes. Therefore, profound modifications of mutant spectrum composition and consensus sequence diversification are not exclusively dependent on environmental alterations or the intervention of population bottlenecks.

Viral fitness: history and relevance for viral pathogenesis and antiviral interventions

The quasispecies dynamics of viral populations (continuous generation of variant genomes and competition among them) has as one of its frequent consequences variations in overall multiplication capacity, a major component of viral fitness. This parameter has multiple implications for viral pathogenesis and viral disease control, some of them unveiled thanks to deep sequencing of viral populations. Darwinian fitness is an old concept whose quantification dates back to the early developments of population genetics. It was later applied to viruses (mainly to RNA viruses) to quantify relative multiplication capacities of individual mutant clones or complex populations. The present article reviews the fitness concept and its relevance for the understanding of the adaptive dynamics of viruses in constant and changing environments. Many studies have addressed the fitness cost of escape mutations (to antibodies, cytotoxic T cells or inhibitors) as an influence on the efficacy of antiviral interventions. Here, we summarize the evidence that the basal fitness level can be a determinant of inhibitor resistance.

Population- and Variant-Based Genome Analyses of Viruses from Vaccine-Derived Rabies Cases Demonstrate Product Specific Clusters and Unique Patterns

Rabies in wildlife has been successfully controlled in parts of Europe and North America using oral rabies vaccination, i.e., the distribution of baits containing live-attenuated virus strains. Occasionally, these vaccines caused vaccine virus-induced rabies cases. To elucidate the mechanisms of genetic selection and the effect of viral populations on these rabies cases, a next generation sequencing approach as well as comprehensive data analyses of the genetic diversity of Street Alabama Dufferin (SAD) and ERA vaccine virus strains and vaccine-induced rabies cases from Canada and several European countries were conducted. As a result, twelve newly generated sets of sequencing data from Canada and Poland were added to a pool of previously investigated samples. While the population-based analysis showed a segregation of viruses of ERA vaccine-induced rabies cases from those of SAD Bern original (SAD Bern_orig)-derived rabies cases, the in-depth variant analysis revealed three distinct combinations of selected variants for the ERA vaccine-induced cases, suggesting the presence of multiple replication-competent haplotypes in the investigated ERA-BHK21 vaccine. Our findings demonstrate the potential of a deep sequencing approach in combination with comprehensive analyses on the consensus, population, and variant level.

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

Defective influenza virus particles generated during viral replication carry incomplete viral genomes and can interfere with the replication of competent viruses. These defective genomes are thought to modulate the disease severity and pathogenicity of an influenza virus infection. Different defective viral genomes also introduce another source of variation across a heterogeneous cell population. Evaluating the impact of defective virus genomes on host cell responses cannot be fully resolved at the population level, requiring single-cell transcriptional profiling. Here, we characterized virus and host transcriptomes in individual influenza virus-infected cells, including those of defective viruses that arise during influenza A virus infection. We established an association between defective virus transcription and host responses and validated interfering and immunostimulatory functions of identified dominant defective viral genome species in vitro. This study demonstrates the intricate effects of defective viral genomes on host transcriptional responses and highlights the importance of capturing host-virus interactions at the single-cell level.

Virus and host factors contribute to cell-to-cell variation in viral infections and determine the outcome of the overall infection. However, the extent of the variability at the single-cell level and how it impacts virus-host interactions at a system level are not well understood. To characterize the dynamics of viral transcription and host responses, we used single-cell RNA sequencing to quantify at multiple time points the host and viral transcriptomes of human A549 cells and primary bronchial epithelial cells infected with influenza A virus. We observed substantial variability in viral transcription between cells, including the accumulation of defective viral genomes (DVGs) that impact viral replication. We show (i) a correlation between DVGs and virus-induced variation of the host transcriptional program and (ii) an association between differential inductions of innate immune response genes and attenuated viral transcription in subpopulations of cells. These observations at the single-cell level improve our understanding of the complex virus-host interplay during influenza virus infection.

Collective Viral Spread Mediated by Virion Aggregates Promotes the Evolution of Defective Interfering Particles

A growing number of studies report that viruses can spread in groups in so-called collective infectious units. By increasing the cellular multiplicity of infection, collective dispersal may allow for social-like interactions, such as cooperation or cheating. Yet, little is known about how such interactions evolve. In previous work with vesicular stomatitis virus, we showed that virion aggregation accelerates early infection stages in most cell types, providing a short-term fitness benefit to the virus. Here, we examine the effects of virion aggregation over several infection cycles. Flow cytometry, deep sequencing, infectivity assays, reverse transcription-quantitative PCR, and electron microscopy revealed that virion aggregation rapidly promotes the emergence of defective interfering particles. Therefore, virion aggregation provides immediate fitness benefits to the virus but incurs fitness costs after a few viral generations. This suggests that an optimal strategy for the virus is to undergo virion aggregation only episodically, for instance, during interhost transmission.IMPORTANCE Recent insights have revealed that viruses use a highly diverse set of strategies to release multiple viral genomes into the same target cells, allowing the emergence of beneficial, but also detrimental, interactions among viruses inside infected cells. This has prompted interest among microbial ecologists and evolutionary biologists in studying how collective dispersal impacts the outcome of viral infections. Here, we have used vesicular stomatitis virus as a model system to study the evolutionary implications of collective dissemination mediated by viral aggregates, since this virus can spontaneously aggregate in the presence of saliva. We find that saliva-driven aggregation has a dual effect on viral fitness; whereas aggregation tends to increase infectivity in the very short term, virion aggregates are highly susceptible to invasion by noncooperative defective variants after a few viral generations.

Analysis and Sorting of Individual HSV-1 Particles by Flow Virometry

Flow cytometry has been instrumental in characterizing normal and infected cells. However, until recently, it was not possible to use such an approach to analyze small entities such as bacteria, let alone viruses, owing to the 0.5 μm resolution of most instruments. To circumvent this limitation, some laboratories decorate pathogens with antibodies or nanoparticles. Our laboratory instead exploits an alternative approach that relies on the staining of internal viral constituents with permeable SYTO dyes or the fluorescent tagging of individual viral proteinaceous components, whether capsid, tegument or glycoproteins. This opens up a range of new research avenues and, for example, enabled us to characterize individual herpes simplex virus type 1 particles, discern their different subpopulations, measure the heterogeneity of mature virions in terms of protein content, sort these viral particles with >90% purity and, for the first time, directly address the impact of this heterogeneity on viral fitness. This approach, coined flow virometry or nanoscale flow cytometry, allows for the study of a wide variety of pathogens with high statistical significance and the potential discovery of novel virulence factors.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Direct sequencing of RNA with MinION Nanopore: detecting mutations based on associations

One of the key challenges in the field of genetics is the inference of haplotypes from next generation sequencing data. The MinION Oxford Nanopore sequencer allows sequencing long reads, with the potential of sequencing complete genes, and even complete genomes of viruses, in individual reads. However, MinION suffers from high error rates, rendering the detection of true variants difficult. Here, we propose a new statistical approach named AssociVar, which differentiates between true mutations and sequencing errors from direct RNA/DNA sequencing using MinION. Our strategy relies on the assumption that sequencing errors will be dispersed randomly along sequencing reads, and hence will not be associated with each other, whereas real mutations will display a non-random pattern of association with other mutations. We demonstrate our approach using direct RNA sequencing data from evolved populations of the MS2 bacteriophage, whose small genome makes it ideal for MinION sequencing. AssociVar inferred several mutations in the phage genome, which were corroborated using parallel Illumina sequencing. This allowed us to reconstruct full genome viral haplotypes constituting different strains that were present in the sample. Our approach is applicable to long read sequencing data from any organism for accurate detection of bona fide mutations and inter-strain polymorphisms.

How pirate phage interferes with helper phage: Comparison of the two distinct strategies

Pirate phages use the structural proteins encoded by unrelated helper phages to propagate. The best-studied example is the pirate P4 and helper P2 of coliphages, and it has been known that the Staphylococcus aureus pathogenicity islands (SaPIs) that can encode virulence factors act as pirate phages, too. When alone in the host, the pirate phages act as a prophage, but when the helper phage gene is also in the same host cell, the pirate phage has ability to exploit the helper phages structural proteins to produce pirate phage particles and spread, interfering with the helper phage production. The known helper phages in these systems are temperate phages. Interestingly, the interference of the pirate phage to the helper phage occurs in a different manner between the SaPI-helper system and the P4-P2 system. SaPIs cannot lyse a helper lysogen upon infection, while when a helper phage lyse a SaPI lysogen, most of the phage particles produced are the SaPI particles. On the contrary, in the P4-P2 system, a pirate phage P4 can lyse a helper P2 lysogen to produce mostly the P4 particles, while when P2 phage lyses a P4 lysogen, most of the produced phages are the P2 particles. Here, the consequences of these different strategies in the pirate and helper phage spreading among uninfected host is analyzed by using mathematical models. It is found that SaPI's strategy interferes with the helper phage spreading significantly more than the P4's strategy, because SaPI interferes with the helper phage's main reproduction step, while P4 interferes only by forcing the helper lysogens to lyse. However, the interference is found to be weaker in the spatially structured environment than in the well-mixed environment. This is because, in the spatial setting, the system tends to self-organize so that the helper phages take over the front of propagation due to the need of helper phage for the pirate phage spreading.

A population of bang-bang switches of defective interfering particles makes within-host dynamics of dengue virus controllable

The titre of virus in a dengue patient and the duration of this viraemia has a profound effect on whether or not a mosquito will become infected when it feeds on the patient and this, in turn, is a key driver of the magnitude of a dengue outbreak. The assessment of the heterogeneity of viral dynamics in dengue-infected patients and its precise treatment are still uncertain. Infection onset, patient physiology and immune response are thought to play major roles in the development of the viral load. Research has explored the interference and spontaneous generation of defective virus particles, but have not examined both the antibody and defective particles during natural infection. We explore the intrinsic variability in the within-host dynamics of viraemias for a population of patients using the method of population of models (POMs). A dataset from 208 patients is used to initially calibrate 20,000 models for the infection kinetics for each of the four dengue virus serotypes. The calibrated POMs suggests that naturally generated defective particles may interfere with the viraemia, but the generated defective virus particles are not adequate to reduce high fever and viraemia duration. The effect of adding excess defective dengue virus interfering particles to patients as a therapeutic is evaluated using the calibrated POMs in a bang-bang (on-off or two-step) optimal control setting. Bang-bang control is a class of binary feedback control that turns either 'ON' or 'OFF' at different time points, determined by the system feedback. Here, the bang-bang control estimates the mathematically optimal dose and duration of the intervention for each model in the POM set.

Membrane-Associated Enteroviruses Undergo Intercellular Transmission as Pools of Sibling Viral Genomes

Some viruses are released from cells as pools of membrane-associated virions. By increasing the multiplicity of infection (MOI), this type of collective dispersal could favor viral cooperation, but also the emergence of cheater-like viruses such as defective interfering particles. To better understand this process, we examined the genetic diversity of membrane-associated coxsackievirus infectious units. We find that infected cells release membranous structures (including vesicles) that contain 8-21 infectious particles on average. However, in most cases (62%-93%), these structures do not promote the co-transmission of different viral genetic variants present in a cell. Furthermore, collective dispersal has no effect on viral population sequence diversity. Our results indicate that membrane-associated collective infectious units typically contain viral particles derived from the same parental genome. Hence, if cooperation occurs, it should probably involve sibling viral particles rather than different variants. As shown by social evolution theory, cooperation among siblings should be robust against cheater invasion.