Drivers of within-host genetic diversity in acute infections of viruses

Navigating a Fine Balance: Point-Mutant Cheater Viruses Disrupt the Viral Replication Cycle

Cheater viruses cannot replicate on their own yet replicate faster than the wild type (WT) when the 2 viruses coinfect the same cell. Cheaters must possess dual genetic features: a defect, which leads to their inability to infect cells on their own, and a selective advantage over WT during coinfection. Previously, we have discovered 2 point-mutant cheaters of the MS2 bacteriophage. Here, we set out to discover the possible repertoire of cheater MS2 viruses by performing experimental evolution at a very high multiplicity of infection. Our results revealed a third point-mutant cheater that arose in 8 biological replicas. Each of the 3 primary cheaters disrupts the fine balance necessary for phage replication, in different ways that create a defect + advantage. We found that over time, the point-mutant cheaters accumulate additional secondary mutations, which alter other stages of the viral replication cycle, complementing the disruptions created by the original cheater. Intriguingly, cheater and secondary mutations almost always reside in very close proximity on the genome. This region encodes for multiple functions: overlapping reading frames as well as overlapping RNA structures critical for transitioning from one stage to another in the viral replication cycle. This region of overlap explains the dual functions of cheaters, as one mutation can have pleiotropic effects. Overall, these findings underscore how viruses, whose dense genomes often have overlapping functions, can easily evolve point-mutant cheaters, and how cheaters can evolve to alter the intricate balance of the viral replication cycle.

The ageing virus hypothesis: Epigenetic ageing beyond the Tree of Life

A recent thought-provoking theory argues that complex organisms using epigenetic information for their normal development and functioning must irreversibly age as a result of epigenetic signal loss. Importantly, the scope of this theory could be considerably expanded, with scientific benefits, by analyzing epigenetic ageing beyond the borders of the Tree of Life. Viruses that use epigenetic signals for their normal functioning may also age, that is, present an increasing risk of failing to complete their individual life cycle and to disappear with time. As viruses are ancient, abundant, and infect a considerable diversity of hosts, the ageing virus hypothesis, if verified, would have important consequences for many fields of the Life sciences. Uncovering ageing viruses would integrate the most abundant and biologically central entities on Earth into theories of ageing, enhance virology, gerontology, evolutionary biology, molecular ecology, genomics, and possibly medicine through the development of new therapies manipulating viral ageing.

Apollo: A comprehensive GPU-powered within-host simulator for viral evolution and infection dynamics across population, tissue, and cell

Modern sequencing instruments bring unprecedented opportunity to study within-host viral evolution in conjunction with viral transmissions between hosts. However, no computational simulators are available to assist the characterization of within-host dynamics. This limits our ability to interpret epidemiological predictions incorporating within-host evolution and to validate computational inference tools. To fill this need we developed Apollo, a GPU-accelerated, out-of-core tool for within-host simulation of viral evolution and infection dynamics across population, tissue, and cellular levels. Apollo is scalable to hundreds of millions of viral genomes and can handle complex demographic and population genetic models. Apollo can replicate real within-host viral evolution; accurately recapturing observed viral sequences from an HIV cohort derived from initial population-genetic configurations. For practical applications, using Apollo-simulated viral genomes and transmission networks, we validated and uncovered the limitations of a widely used viral transmission inference tool.

Reovirus infection induces transcriptome-wide unique A-to-I editing changes in the murine fibroblasts

The conversion of Adenosine (A) to Inosine (I), by Adenosine Deaminases Acting on RNA or ADARs, is an essential post-transcriptional modification that contributes to proteome diversity and regulation in metazoans including humans. In addition to its transcriptome-regulating role, ADARs also play a major part in immune response to viral infection, where an interferon response activates interferon-stimulated genes, such as ADARp150, in turn dynamically regulating host-virus interactions. A previous report has shown that infection from reoviruses, despite strong activation of ADARp150, does not influence the editing of some of the major known editing targets, while likely editing others, suggesting a potentially nuanced editing pattern that may depend on different factors. However, the results were based on a handful of selected editing sites and did not cover the entire transcriptome. Thus, to determine whether and how reovirus infection specifically affects host ADAR editing patterns, we analyzed a publicly available deep-sequenced RNA-seq dataset, from murine fibroblasts infected with wild-type and mutant reovirus strains that allowed us to examine changes in editing patterns on a transcriptome-wide scale. To the best of our knowledge, this is the first transcriptome-wide report on host editing changes after reovirus infection. Our results demonstrate that reovirus infection induces unique nuanced editing changes in the host, including introducing sites uniquely edited in infected samples. Genes with edited sites are overrepresented in pathways related to immune regulation, cellular signaling, metabolism, and growth. Moreover, a shift in editing targets has also been observed, where the same genes are edited in infection and control conditions but at different sites, or where the editing rate is increased for some and decreased for other differential targets, supporting the hypothesis of dynamic and condition-specific editing by ADARs.

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.

Open questions in the social lives of viruses

Social interactions among viruses occur whenever multiple viral genomes infect the same cells, hosts, or populations of hosts. Viral social interactions range from cooperation to conflict, occur throughout the viral world, and affect every stage of the viral lifecycle. The ubiquity of these social interactions means that they can determine the population dynamics, evolutionary trajectory, and clinical progression of viral infections. At the same time, social interactions in viruses raise new questions for evolutionary theory, providing opportunities to test and extend existing frameworks within social evolution. Many opportunities exist at this interface: Insights into the evolution of viral social interactions have immediate implications for our understanding of the fundamental biology and clinical manifestation of viral diseases. However, these opportunities are currently limited because evolutionary biologists only rarely study social evolution in viruses. Here, we bridge this gap by (1) summarizing the ways in which viruses can interact socially, including consequences for social evolution and evolvability; (2) outlining some open questions raised by viruses that could challenge concepts within social evolution theory; and (3) providing some illustrative examples, data sources, and conceptual questions, for studying the natural history of social viruses.

Mutation rate, selection, and epistasis inferred from RNA virus haplotypes via neural posterior estimation

RNA viruses are particularly notorious for their high levels of genetic diversity, which is generated through the forces of mutation and natural selection. However, disentangling these two forces is a considerable challenge, and this may lead to widely divergent estimates of viral mutation rates, as well as difficulties in inferring the fitness effects of mutations. Here, we develop, test, and apply an approach aimed at inferring the mutation rate and key parameters that govern natural selection, from haplotype sequences covering full-length genomes of an evolving virus population. Our approach employs neural posterior estimation, a computational technique that applies simulation-based inference with neural networks to jointly infer multiple model parameters. We first tested our approach on synthetic data simulated using different mutation rates and selection parameters while accounting for sequencing errors. Reassuringly, the inferred parameter estimates were accurate and unbiased. We then applied our approach to haplotype sequencing data from a serial passaging experiment with the MS2 bacteriophage, a virus that parasites Escherichia coli. We estimated that the mutation rate of this phage is around 0.2 mutations per genome per replication cycle (95% highest density interval: 0.051-0.56). We validated this finding with two different approaches based on single-locus models that gave similar estimates but with much broader posterior distributions. Furthermore, we found evidence for reciprocal sign epistasis between four strongly beneficial mutations that all reside in an RNA stem loop that controls the expression of the viral lysis protein, responsible for lysing host cells and viral egress. We surmise that there is a fine balance between over- and underexpression of lysis that leads to this pattern of epistasis. To recap, we have developed an approach for joint inference of the mutation rate and selection parameters from full haplotype data with sequencing errors and used it to reveal features governing MS2 evolution.

Host immunological responses facilitate development of SARS-CoV-2 mutations in patients receiving monoclonal antibody treatments

BackgroundThe role of host immunity in emergence of evasive SARS-CoV-2 Spike mutations under therapeutic monoclonal antibody (mAb) pressure remains to be explored.MethodsIn a prospective, observational, monocentric ORCHESTRA cohort study, conducted between March 2021 and November 2022, mild-to-moderately ill COVID-19 patients (n = 204) receiving bamlanivimab, bamlanivimab/etesevimab, casirivimab/imdevimab, or sotrovimab were longitudinally studied over 28 days for viral loads, de novo Spike mutations, mAb kinetics, seroneutralization against infecting variants of concern, and T cell immunity. Additionally, a machine learning-based circulating immune-related biomarker (CIB) profile predictive of evasive Spike mutations was constructed and confirmed in an independent data set (n = 19) that included patients receiving sotrovimab or tixagevimab/cilgavimab.ResultsPatients treated with various mAbs developed evasive Spike mutations with remarkable speed and high specificity to the targeted mAb-binding sites. Immunocompromised patients receiving mAb therapy not only continued to display significantly higher viral loads, but also showed higher likelihood of developing de novo Spike mutations. Development of escape mutants also strongly correlated with neutralizing capacity of the therapeutic mAbs and T cell immunity, suggesting immune pressure as an important driver of escape mutations. Lastly, we showed that an antiinflammatory and healing-promoting host milieu facilitates Spike mutations, where 4 CIBs identified patients at high risk of developing escape mutations against therapeutic mAbs with high accuracy.ConclusionsOur data demonstrate that host-driven immune and nonimmune responses are essential for development of mutant SARS-CoV-2. These data also support point-of-care decision making in reducing the risk of mAb treatment failure and improving mitigation strategies for possible dissemination of escape SARS-CoV-2 mutants.FundingThe ORCHESTRA project/European Union's Horizon 2020 research and innovation program.

Global genetic diversity and evolutionary patterns among Potato leafroll virus populations

Potato leafroll virus (PLRV) is a widespread and one of the most damaging viral pathogens causing significant quantitative and qualitative losses in potato worldwide. The current knowledge of the geographical distribution, standing genetic diversity and the evolutionary patterns existing among global PLRV populations is limited. Here, we employed several bioinformatics tools and comprehensively analyzed the diversity, genomic variability, and the dynamics of key evolutionary factors governing the global spread of this viral pathogen. To date, a total of 84 full-genomic sequences of PLRV isolates have been reported from 22 countries with most genomes documented from Kenya. Among all PLRV-encoded major proteins, RTD and P0 displayed the highest level of nucleotide variability. The highest percentage of mutations were associated with RTD (38.81%) and P1 (31.66%) in the coding sequences. We detected a total of 10 significantly supported recombination events while the most frequently detected ones were associated with PLRV genome sequences reported from Kenya. Notably, the distribution patterns of recombination breakpoints across different genomic regions of PLRV isolates remained variable. Further analysis revealed that with exception of a few positively selected codons, a major part of the PLRV genome is evolving under strong purifying selection. Protein disorder prediction analysis revealed that CP-RTD had the highest percentage (48%) of disordered amino acids and the majority (27%) of disordered residues were positioned at the C-terminus. These findings will extend our current knowledge of the PLRV geographical prevalence, genetic diversity, and evolutionary factors that are presumably shaping the global spread and successful adaptation of PLRV as a destructive potato pathogen to geographically isolated regions of the world.

A general approach to identify low-frequency variants within influenza samples collected during routine surveillance

Influenza viruses exhibit considerable diversity between hosts. Additionally, different quasispecies can be found within the same host. High-throughput sequencing technologies can be used to sequence a patient-derived virus population at sufficient depths to identify low-frequency variants (LFV) present in a quasispecies, but many challenges remain for reliable LFV detection because of experimental errors introduced during sample preparation and sequencing. High genomic copy numbers and extensive sequencing depths are required to differentiate false positive from real LFV, especially at low allelic frequencies (AFs). This study proposes a general approach for identifying LFV in patient-derived samples obtained during routine surveillance. Firstly, validated thresholds were determined for LFV detection, whilst balancing both the cost and feasibility of reliable LFV detection in clinical samples. Using a genetically well-defined population of influenza A viruses, thresholds of at least 10⁴ genomes per microlitre and AF of ≥5 % were established as detection limits. Secondly, a subset of 59 retained influenza A (H3N2) samples from the 2016-2017 Belgian influenza season was composed. Thirdly, as a proof of concept for the added value of LFV for routine influenza monitoring, potential associations between patient data and whole genome sequencing data were investigated. A significant association was found between a high prevalence of LFV and disease severity. This study provides a general methodology for influenza LFV detection, which can also be adopted by other national influenza reference centres and for other viruses such as SARS-CoV-2. Additionally, this study suggests that the current relevance of LFV for routine influenza surveillance programmes might be undervalued.

Drivers of adaptive evolution during chronic SARS-CoV-2 infections

In some immunocompromised patients with chronic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, considerable adaptive evolution occurs. Some substitutions found in chronic infections are lineage-defining mutations in variants of concern (VOCs), which has led to the hypothesis that VOCs emerged from chronic infections. In this study, we searched for drivers of VOC-like emergence by consolidating sequencing results from a set of 27 chronic infections. Most substitutions in this set reflected lineage-defining VOC mutations; however, a subset of mutations associated with successful global transmission was absent from chronic infections. We further tested the ability to associate antibody evasion mutations with patient-specific and virus-specific features and found that viral rebound is strongly correlated with the emergence of antibody evasion. We found evidence for dynamic polymorphic viral populations in most patients, suggesting that a compromised immune system selects for antibody evasion in particular niches in a patient’s body. We suggest that a tradeoff exists between antibody evasion and transmissibility and that extensive monitoring of chronic infections is necessary to further understanding of VOC emergence.

Analysis of mutations that arise in chronic SARS-CoV-2 infections shows both overlap and differences with mutations present in pandemic viral variants of concern, highlighting their distinct drivers of evolution.

Thermal stress triggers productive viral infection of a key coral reef symbiont

Climate change-driven ocean warming is increasing the frequency and severity of bleaching events, in which corals appear whitened after losing their dinoflagellate endosymbionts (family Symbiodiniaceae). Viral infections of Symbiodiniaceae may contribute to some bleaching signs, but little empirical evidence exists to support this hypothesis. We present the first temporal analysis of a lineage of Symbiodiniaceae-infecting positive-sense single-stranded RNA viruses ("dinoRNAVs") in coral colonies, which were exposed to a 5-day heat treatment (+2.1 °C). A total of 124 dinoRNAV major capsid protein gene "aminotypes" (unique amino acid sequences) were detected from five colonies of two closely related Pocillopora-Cladocopium (coral-symbiont) combinations in the experiment; most dinoRNAV aminotypes were shared between the two coral-symbiont combinations (64%) and among multiple colonies (82%). Throughout the experiment, seventeen dinoRNAV aminotypes were found only in heat-treated fragments, and 22 aminotypes were detected at higher relative abundances in heat-treated fragments. DinoRNAVs in fragments of some colonies exhibited higher alpha diversity and dispersion under heat stress. Together, these findings provide the first empirical evidence that exposure to high temperatures triggers some dinoRNAVs to switch from a persistent to a productive infection mode within heat-stressed corals. Over extended time frames, we hypothesize that cumulative dinoRNAV production in the Pocillopora-Cladocopium system could affect colony symbiotic status, for example, by decreasing Symbiodiniaceae densities within corals. This study sets the stage for reef-scale investigations of dinoRNAV dynamics during bleaching events.

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

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.

Facilitating Antiviral Drug Discovery Using Genetic and Evolutionary Knowledge

Over the course of human history, billions of people worldwide have been infected by various viruses. Despite rapid progress in the development of biomedical techniques, it is still a significant challenge to find promising new antiviral targets and drugs. In the past, antiviral drugs mainly targeted viral proteins when they were used as part of treatment strategies. Since the virus mutation rate is much faster than that of the host, such drugs feature drug resistance and narrow-spectrum antiviral problems. Therefore, the targeting of host molecules has gradually become an important area of research for the development of antiviral drugs. In recent years, rapid advances in high-throughput sequencing techniques have enabled numerous genetic studies (such as genome-wide association studies (GWAS), clustered regularly interspersed short palindromic repeats (CRISPR) screening, etc.) for human diseases, providing valuable genetic and evolutionary resources. Furthermore, it has been revealed that successful drug targets exhibit similar genetic and evolutionary features, which are of great value in identifying promising drug targets and discovering new drugs. Considering these developments, in this article the authors propose a host-targeted antiviral drug discovery strategy based on knowledge of genetics and evolution. We first comprehensively summarized the genetic, subcellular location, and evolutionary features of the human genes that have been successfully used as antiviral targets. Next, the summarized features were used to screen novel druggable antiviral targets and to find potential antiviral drugs, in an attempt to promote the discovery of new antiviral drugs.

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.

Evidence for increased breakthrough rates of SARS-CoV-2 variants of concern in BNT162b2-mRNA-vaccinated individuals

The BNT162b2 mRNA vaccine is highly effective against SARS-CoV-2. However, apprehension exists that variants of concern (VOCs) may evade vaccine protection, due to evidence of reduced neutralization of the VOCs B.1.1.7 and B.1.351 by vaccine sera in laboratory assays. We performed a matched cohort study to examine the distribution of VOCs in infections of BNT162b2 mRNA vaccinees from Clalit Health Services (Israel) using viral genomic sequencing, and hypothesized that if vaccine effectiveness against a VOC is reduced, its proportion among breakthrough cases would be higher than in unvaccinated controls. Analyzing 813 viral genome sequences from nasopharyngeal swabs, we showed that vaccinees who tested positive at least 7 days after the second dose were disproportionally infected with B.1.351, compared with controls. Those who tested positive between 2 weeks after the first dose and 6 days after the second dose were disproportionally infected by B.1.1.7. These findings suggest reduced vaccine effectiveness against both VOCs within particular time windows. Our results emphasize the importance of rigorously tracking viral variants, and of increasing vaccination to prevent the spread of VOCs.

SARS-CoV-2 within-host diversity and transmission

Extensive global sampling and sequencing of the pandemic virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have enabled researchers to monitor its spread and to identify concerning new variants. Two important determinants of variant spread are how frequently they arise within individuals and how likely they are to be transmitted. To characterize within-host diversity and transmission, we deep-sequenced 1313 clinical samples from the United Kingdom. SARS-CoV-2 infections are characterized by low levels of within-host diversity when viral loads are high and by a narrow bottleneck at transmission. Most variants are either lost or occasionally fixed at the point of transmission, with minimal persistence of shared diversity, patterns that are readily observable on the phylogenetic tree. Our results suggest that transmission-enhancing and/or immune-escape SARS-CoV-2 variants are likely to arise infrequently but could spread rapidly if successfully transmitted.

Full genome viral sequences inform patterns of SARS-CoV-2 spread into and within Israel

Full genome sequences are increasingly used to track the geographic spread and transmission dynamics of viral pathogens. Here, with a focus on Israel, we sequence 212 SARS-CoV-2 sequences and use them to perform a comprehensive analysis to trace the origins and spread of the virus. We find that travelers returning from the United States of America significantly contributed to viral spread in Israel, more than their proportion in incoming infected travelers. Using phylodynamic analysis, we estimate that the basic reproduction number of the virus was initially around 2.5, dropping by more than two-thirds following the implementation of social distancing measures. We further report high levels of transmission heterogeneity in SARS-CoV-2 spread, with between 2-10% of infected individuals resulting in 80% of secondary infections. Overall, our findings demonstrate the effectiveness of social distancing measures for reducing viral spread.