Mutational and fitness landscapes of an RNA virus revealed through population sequencing

Difference of Intrahost Dynamics of the Second Human Pegivirus and Hepatitis C Virus in HPgV-2/HCV-Coinfected Patients

Background

The second human pegivirus (HPgV-2) and hepatitis C virus (HCV) belong to the Flaviviridae family and share some common genome features. However, the two viruses exhibit significantly different genetic diversity. The comparison of intrahost dynamics of HPgV-2 and HCV that mainly reflect virus-host interactions is needed to elucidate their intrahost difference of genetic diversity and the possible mechanisms.

Methods

Intrahost single nucleotide variations (iSNVs) were identified by means of next-generation sequencing from both cross-sectional and longitudinal samples from HPgV-2- and HCV-coinfected patients. The levels of human cytokines were quantified in the patient before and after HCV elimination by the treatment of direct-acting antivirals (DAA).

Results

Unlike HCV, the viral sequences of HPgV-2 are highly conserved among HPgV-2-infected patients. However, iSNV analysis confirmed the intrahost variation or quasispecies of HPgV-2. Almost all iSNVs of HPgV-2 did not accumulate or transmit within host over time, which may explain the highly conserved HPgV-2 consensus sequence. Intrahost variation of HPgV-2 mainly causes nucleotide transition in particular at the 3rd codon position and synonymous substitutions, indicating purifying or negative selection posed by host immune system. Cytokine data further indicate that HPgV-2 infection alone may not efficiently stimulate innate immune responses since proinflammatory cytokine expression dramatically decreased with elimination of HCV.

Conclusion

This study provided new insights into the intrahost genomic variations and evolutionary dynamics of HPgV-2 as well as the impact of host immune selection and virus polymerase on virus evolution. The different genetic diversity of HPgV-2 and HCV makes HPgV-2 a potential new model to investigate RNA virus diversity and the mechanism of viral polymerase in modulating virus replication.

High prevalence of tryptophan-truncated S quasispecies in treatment-naïve chronic hepatitis B patients

Hepatitis B virus surface antigen (HBsAg) encoded by the S gene is highly expressed during the replication cycle of hepatitis B virus (HBV). However, the frequent usage of tryptophan in HBsAg, which leads to a high cost of biosynthesis, is inconsistent with the high expression level of this protein. Tryptophan-truncated mutation of HBsAg, that is, a tryptophan to stop codon mutation resulting in truncated HBsAg, might help to maintain its high expression with lower biosynthetic cost. We aimed to investigate the prevalence of tryptophan-truncated S quasispecies in treatment-naïve patients with chronic hepatitis B (CHB) by applying CirSeq as well as a site-by-site algorithm developed by us to identify variants at extremely low frequencies in the carboxyl terminus of HBsAg. A total of 730 mutations were identified in 27 patients with CHB, varying from seven to 56 mutations per sample. The number of synonymous mutations was much higher than that of nonsynonymous mutations in the reverse transcriptase (RT) coding region and vice versa in the S coding region, implying that the evolutionary constraints on the RT and S genes might be different. We showed that 25 (92.6 %) of 27 patients had at least one S-truncated mutation, most of which were derived from tryptophan, indicating a high prevalence of tryptophan-truncated S mutations in treatment-naïve patients with CHB. In terms of the RT gene, 21 (77.8 %) patients had pre-existing drug-resistant mutations, while no truncated mutations were detected. Our findings that tryptophan-truncated S quasispecies and drug-resistant RT mutants were highly prevalent in treatment-naïve patients with CHB provide new insights into the composition of the HBV population, which might help optimize the treatment and management of patients with CHB.

A comparative analysis of parechovirus protein structures with other picornaviruses

Parechoviruses belong to the genus Parechovirus within the family Picornaviridae and are non-enveloped icosahedral viruses with a single-stranded RNA genome. Parechoviruses include human and animal pathogens classified into six species. Those that infect humans belong to the Parechovirus A species and can cause infections ranging from mild gastrointestinal or respiratory illness to severe neonatal sepsis. There are no approved antivirals available to treat parechovirus (nor any other picornavirus) infections. In this parechovirus review, we focus on the cleaved protein products resulting from the polyprotein processing after translation comparing and contrasting their known or predicted structures and functions to those of other picornaviruses. The review also includes our original analysis from sequence and structure prediction. This review highlights significant structural differences between parechoviral and other picornaviral proteins, suggesting that parechovirus drug development should specifically be directed to parechoviral targets.

The Molecular Evolution of Type 2 Vaccine-Derived Polioviruses in Individuals with Primary Immunodeficiency Diseases

The oral poliovirus vaccine (OPV), which prevents person-to-person transmission of poliovirus by inducing robust intestinal immunity, has been a crucial tool for global polio eradication. However, polio outbreaks, mainly caused by type 2 circulating vaccine-derived poliovirus (cVDPV2), are increasing worldwide. Meanwhile, immunodeficiency-associated vaccine-derived poliovirus (iVDPV) is considered another risk factor during the final stage of global polio eradication. Patients with primary immunodeficiency diseases are associated with higher risks for long-term iVDPV infections. Although a limited number of chronic iVDPV excretors were reported, the recent identification of a chronic type 2 iVDPV (iVDPV2) excretor in the Philippines highlights the potential risk of inapparent iVDPV infection for expanding cVDPV outbreaks. Further research on the genetic characterizations and molecular evolution of iVDPV2, based on comprehensive iVDPV surveillance, will be critical for elucidating the remaining risk of iVDPV2 during the post-OPV era.

Fitness Estimation for Viral Variants in the Context of Cellular Coinfection

Animal models are frequently used to characterize the within-host dynamics of emerging zoonotic viruses. More recent studies have also deep-sequenced longitudinal viral samples originating from experimental challenges to gain a better understanding of how these viruses may evolve in vivo and between transmission events. These studies have often identified nucleotide variants that can replicate more efficiently within hosts and also transmit more effectively between hosts. Quantifying the degree to which a mutation impacts viral fitness within a host can improve identification of variants that are of particular epidemiological concern and our ability to anticipate viral adaptation at the population level. While methods have been developed to quantify the fitness effects of mutations using observed changes in allele frequencies over the course of a host’s infection, none of the existing methods account for the possibility of cellular coinfection. Here, we develop mathematical models to project variant allele frequency changes in the context of cellular coinfection and, further, integrate these models with statistical inference approaches to demonstrate how variant fitness can be estimated alongside cellular multiplicity of infection. We apply our approaches to empirical longitudinally sampled H5N1 sequence data from ferrets. Our results indicate that previous studies may have significantly underestimated the within-host fitness advantage of viral variants. These findings underscore the importance of considering the process of cellular coinfection when studying within-host viral evolutionary dynamics.

The Adaptive Potential of the Middle Domain of Yeast Hsp90

The distribution of fitness effects (DFEs) of new mutations across different environments quantifies the potential for adaptation in a given environment and its cost in others. So far, results regarding the cost of adaptation across environments have been mixed, and most studies have sampled random mutations across different genes. Here, we quantify systematically how costs of adaptation vary along a large stretch of protein sequence by studying the distribution of fitness effects of the same ≈2,300 amino-acid changing mutations obtained from deep mutational scanning of 119 amino acids in the middle domain of the heat shock protein Hsp90 in five environments. This region is known to be important for client binding, stabilization of the Hsp90 dimer, stabilization of the N-terminal-Middle and Middle-C-terminal interdomains, and regulation of ATPase–chaperone activity. Interestingly, we find that fitness correlates well across diverse stressful environments, with the exception of one environment, diamide. Consistent with this result, we find little cost of adaptation; on average only one in seven beneficial mutations is deleterious in another environment. We identify a hotspot of beneficial mutations in a region of the protein that is located within an allosteric center. The identified protein regions that are enriched in beneficial, deleterious, and costly mutations coincide with residues that are involved in the stabilization of Hsp90 interdomains and stabilization of client-binding interfaces, or residues that are involved in ATPase–chaperone activity of Hsp90. Thus, our study yields information regarding the role and adaptive potential of a protein sequence that complements and extends known structural information.

From genotypes to organisms: State-of-the-art and perspectives of a cornerstone in evolutionary dynamics

Understanding how genotypes map onto phenotypes, fitness, and eventually organisms is arguably the next major missing piece in a fully predictive theory of evolution. We refer to this generally as the problem of the genotype-phenotype map. Though we are still far from achieving a complete picture of these relationships, our current understanding of simpler questions, such as the structure induced in the space of genotypes by sequences mapped to molecular structures, has revealed important facts that deeply affect the dynamical description of evolutionary processes. Empirical evidence supporting the fundamental relevance of features such as phenotypic bias is mounting as well, while the synthesis of conceptual and experimental progress leads to questioning current assumptions on the nature of evolutionary dynamics-cancer progression models or synthetic biology approaches being notable examples. This work delves with a critical and constructive attitude into our current knowledge of how genotypes map onto molecular phenotypes and organismal functions, and discusses theoretical and empirical avenues to broaden and improve this comprehension. As a final goal, this community should aim at deriving an updated picture of evolutionary processes soundly relying on the structural properties of genotype spaces, as revealed by modern techniques of molecular and functional analysis.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Rescue of Infectious Sindbis Virus by Yeast Spheroplast-Mammalian Cell Fusion

Sindbis virus (SINV), a positive-sense single stranded RNA virus that causes mild symptoms in humans, is transmitted by mosquito bites. SINV reverse genetics have many implications, not only in understanding alphavirus transmission, replication cycle, and virus-host interactions, but also in biotechnology and biomedical applications. The rescue of SINV infectious particles is usually achieved by transfecting susceptible cells (BHK-21) with SINV-infectious mRNA genomes generated from cDNA constructed via in vitro translation (IVT). That procedure is time consuming, costly, and relies heavily on reagent quality. Here, we constructed a novel infectious SINV cDNA construct that expresses its genomic RNA in yeast cells controlled by galactose induction. Using spheroplasts made from this yeast, we established a robust polyethylene glycol-mediated yeast: BHK-21 fusion protocol to rescue infectious SINV particles. Our approach is timesaving and utilizes common lab reagents for SINV rescue. It could be a useful tool for the rescue of large single strand RNA viruses, such as SARS-CoV-2.

Increased RNA virus population diversity improves adaptability

The replication machinery of most RNA viruses lacks proofreading mechanisms. As a result, RNA virus populations harbor a large amount of genetic diversity that confers them the ability to rapidly adapt to changes in their environment. In this work, we investigate whether further increasing the initial population diversity of a model RNA virus can improve adaptation to a single selection pressure, thermal inactivation. For this, we experimentally increased the diversity of coxsackievirus B3 (CVB3) populations across the capsid region. We then compared the ability of these high diversity CVB3 populations to achieve resistance to thermal inactivation relative to standard CVB3 populations in an experimental evolution setting. We find that viral populations with high diversity are better able to achieve resistance to thermal inactivation at both the temperature employed during experimental evolution as well as at a more extreme temperature. Moreover, we identify mutations in the CVB3 capsid that confer resistance to thermal inactivation, finding significant mutational epistasis. Our results indicate that even naturally diverse RNA virus populations can benefit from experimental augmentation of population diversity for optimal adaptation and support the use of such viral populations in directed evolution efforts that aim to select viruses with desired characteristics.

VERSO: A comprehensive framework for the inference of robust phylogenies and the quantification of intra-host genomic diversity of viral samples

We introduce VERSO, a two-step framework for the characterization of viral evolution from sequencing data of viral genomes, which is an improvement on phylogenomic approaches for consensus sequences. VERSO exploits an efficient algorithmic strategy to return robust phylogenies from clonal variant profiles, also in conditions of sampling limitations. It then leverages variant frequency patterns to characterize the intra-host genomic diversity of samples, revealing undetected infection chains and pinpointing variants likely involved in homoplasies. On simulations, VERSO outperforms state-of-the-art tools for phylogenetic inference. Notably, the application to 6,726 amplicon and RNA sequencing samples refines the estimation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution, while co-occurrence patterns of minor variants unveil undetected infection paths, which are validated with contact tracing data. Finally, the analysis of SARS-CoV-2 mutational landscape uncovers a temporal increase of overall genomic diversity and highlights variants transiting from minor to clonal state and homoplastic variants, some of which fall on the spike gene. Available at: https://github.com/BIMIB-DISCo/VERSO.

The evolutionary origin of the universal distribution of mutation fitness effect

An intriguing fact long defying explanation is the observation of a universal exponential distribution of beneficial mutations in fitness effect for different microorganisms. To explain this effect, we use a population model including mutation, directional selection, linkage, and genetic drift. The multiple-mutation regime of adaptation at large population sizes (traveling wave regime) is considered. We demonstrate analytically and by simulation that, regardless of the inherent distribution of mutation fitness effect across genomic sites, an exponential distribution of fitness effects emerges in the long term. This result follows from the exponential statistics of the frequency of the less-fit alleles, f, that we predict to evolve, in the long term, for both polymorphic and monomorphic sites. We map the logarithmic slope of the distribution onto the previously derived fixation probability and demonstrate that it increases linearly in time. Our results demonstrate a striking difference between the distribution of fitness effects observed experimentally for naturally occurring mutations, and the "inherent" distribution obtained in a directed-mutagenesis experiment, which can have any shape depending on the organism. Based on these results, we develop a new method to measure the fitness effect of mutations for each variable residue using DNA sequences sampled from adapting populations. This new method is not sensitive to linkage effects and does not require the one-site model assumptions.

Norovirus evolution in immunodeficient mice reveals potentiated pathogenicity via a single nucleotide change in the viral capsid

Interferons (IFNs) are key controllers of viral replication, with intact IFN responses suppressing virus growth and spread. Using the murine norovirus (MNoV) system, we show that IFNs exert selective pressure to limit the pathogenic evolutionary potential of this enteric virus. In animals lacking type I IFN signaling, the nonlethal MNoV strain CR6 rapidly acquired enhanced virulence via conversion of a single nucleotide. This nucleotide change resulted in amino acid substitution F514I in the viral capsid, which led to >10,000-fold higher replication in systemic organs including the brain. Pathogenicity was mediated by enhanced recruitment and infection of intestinal myeloid cells and increased extraintestinal dissemination of virus. Interestingly, the trade-off for this mutation was reduced fitness in an IFN-competent host, in which CR6 bearing F514I exhibited decreased intestinal replication and shedding. In an immunodeficient context, a spontaneous amino acid change can thus convert a relatively avirulent viral strain into a lethal pathogen.

Genome Number and Size Polymorphism in Zika Virus Infectious Units

The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero. How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood.

Zika virus (ZIKV; Flaviviridae, Flavivirus) is an arthropod-borne infection that can result in severe outcomes, particularly in fetuses infected in utero. It has been assumed that infection by ZIKV, as well as other viruses, is largely initiated by individual virus particles binding to and entering a cell. However, recent studies have demonstrated that multiple virus particles are frequently delivered to a cell simultaneously and that this collective particle delivery enhances infection. ZIKV is maintained in nature between Aedes aegypti mosquitos and vertebrate hosts, including humans. Human infection is initiated through the injection of a relatively small initial inoculum comprised of a genetically complex virus population. Since most mutations decrease virus fitness, collective particle transmission could benefit ZIKV and other arthropod-borne diseases by facilitating the maintenance of genetic complexity and adaptability during infection or through other mechanisms. Therefore, we utilized a barcoded ZIKV to quantify the number of virus genomes that initiate a plaque. We found that individual plaques contain a mean of 10 infecting viral genomes (range, 1 to 212). Few plaques contained more than two dominant genomes. To determine whether multigenome infectious units consist of collectively transmitting virions, infectious units of ZIKV were then separated mechanically by centrifugation, and heavier fractions were found to contain more genomes per plaque-forming unit, with larger diameters. Finally, larger/heavier infectious units reformed after removal. These data suggest that ZIKV populations consist of a variety of infectious unit sizes, likely mostly made up of aggregates, and only rarely begin with a single virus genome.

IMPORTANCE The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero. How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood. This research demonstrates a novel mechanism of ZIKV transmission through multigenome aggregates, providing insight into ZIKV evolution, immunologic evasion, and better future therapeutic design. This study shows that ZIKV plaques result from collections of genomes rather than individual genomes, increasing the potential for interactions between ZIKV genotypes.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Principles of dengue virus evolvability derived from genotype-fitness maps in human and mosquito cells

Dengue virus (DENV) cycles between mosquito and mammalian hosts. To examine how DENV populations adapt to these different host environments, we used serial passage in human and mosquito cell lines and estimated fitness effects for all single-nucleotide variants in these populations using ultra-deep sequencing. This allowed us to determine the contributions of beneficial and deleterious mutations to the collective fitness of the population. Our analysis revealed that the continuous influx of a large burden of deleterious mutations counterbalances the effect of rare, host-specific beneficial mutations to shape the path of adaptation. Beneficial mutations preferentially map to intrinsically disordered domains in the viral proteome and cluster to defined regions in the genome. These phenotypically redundant adaptive alleles may facilitate host-specific DENV adaptation. Importantly, the evolutionary constraints described in our simple system mirror trends observed across DENV and Zika strains, indicating it recapitulates key biophysical and biological constraints shaping long-term viral evolution.

Addressing Non-linear System Dynamics of Single-Strand RNA Virus-Host Interaction

Positive single-strand ribonucleic acid [(+)ssRNA] viruses can cause multiple outbreaks, for which comprehensive tailored therapeutic strategies are still missing. Virus and host cell dynamics are tightly connected, generating a complex dynamics that conveys in virion assembly to ensure virus spread in the body. Starting from the knowledge of relevant processes in (+ss)RNA virus replication, transcription, translation, virions budding and shedding, and their respective energy costs, we built up a systems thinking (ST)-based diagram of the virus-host interaction, comprehensive of stocks, flows, and processes as well-described in literature. In ST approach, stocks and flows are expressed by a proxy of the energy embedded and transmitted, respectively, whereas processes are referred to the energy required for the system functioning. In this perspective, healthiness is just a particular configuration, in which stocks relevant for the system (equivalent but not limited to proteins, RNA, DNA, and all metabolites required for the survival) are constant, and the system behavior is stationary. At time of infection, the presence of additional stocks (e.g., viral protein and RNA and all metabolites required for virion assembly and spread) confers a complex network of feedbacks leading to new configurations, which can evolve to maximize the virions stock, thus changing the system structure, output, and purpose. The dynamic trajectories will evolve to achieve a new stationary status, a phenomenon described in microbiology as integration and symbiosis when the system is resilient enough to the changes, or the system may stop functioning and die. Application of external driving forces, acting on processes, can affect the dynamic trajectories adding a further degree of complexity, which can be captured by ST approach, used to address these new configurations. Investigation of system configurations in response to external driving forces acting is developed by computational analysis based on ST diagrams, with the aim at designing novel therapeutic approaches.

Globally defining the effects of mutations in a picornavirus capsid

The capsids of non-enveloped viruses are highly multimeric and multifunctional protein assemblies that play key roles in viral biology and pathogenesis. Despite their importance, a comprehensive understanding of how mutations affect viral fitness across different structural and functional attributes of the capsid is lacking. To address this limitation, we globally define the effects of mutations across the capsid of a human picornavirus. Using this resource, we identify structural and sequence determinants that accurately predict mutational fitness effects, refine evolutionary analyses, and define the sequence specificity of key capsid-encoded motifs. Furthermore, capitalizing on the derived sequence requirements for capsid-encoded protease cleavage sites, we implement a bioinformatic approach for identifying novel host proteins targeted by viral proteases. Our findings represent the most comprehensive investigation of mutational fitness effects in a picornavirus capsid to date and illuminate important aspects of viral biology, evolution, and host interactions.

A virus is made up of genetic material that is encased with a protective protein coat called the capsid. The capsid also helps the virus to infect host cells by binding to the host receptor proteins and releasing its genetic material. Inside the cell, the virus hitchhikes the infected cell’s machinery to grow or replicate its own genetic material.

Viral capsids are the main target of the host’s defence system, and therefore, continuously change in an attempt to escape the immune system by introducing alterations (known as mutations) into the genes encoding viral capsid proteins. Mutations occur randomly, and so while some changes to the viral capsid might confer an advantage, others may have no effect at all, or even weaken the virus.

To better understand the effect of capsid mutations on the virus’ ability to infect host cells, Mattenberger et al. studied the Coxsackievirus B3, which is linked to heart problems and acute heart failure in humans. The researchers analysed around 90% of possible amino acid mutations (over 14,800 mutations) and correlated each mutation to how it influenced the virus’ ability to replicate in human cells grown in the laboratory.

Based on these results, Mattenberger et al. developed a computer model to predict how a particular mutation might affect the virus. The analysis also identified specific amino acid sequences of capsid proteins that are essential for certain tasks, such as building the capsid. It also included an analysis of sequences in the capsid that allow it to be recognized by another viral protein, which cuts the capsid proteins into the right size from a larger precursor. By looking for similar sequences in human genes, the researchers identified several ones that the virus may attack and inactivate to support its own replication.

These findings may help identify potential drug targets to develop new antiviral therapies. For example, proteins of the capsid that are less likely to mutate will provide a better target as they lower the possibility of the virus to become resistant to the treatment. They also highlight new proteins in human cells that could potentially block the virus in cells.

Investigating Mitochondrial Transcriptomes and RNA Processing Using Circular RNA Sequencing

Transcriptomic technologies have revolutionized the study of gene expression and RNA biology. Different RNA sequencing methods enable the analyses of diverse species of transcripts, including their abundance, processing, stability, and other specific features. Mitochondrial transcriptomics has benefited from these technologies that have revealed the surprising complexity of its RNAs. Here we describe a method based upon cyclization of mitochondrial RNAs and next generation sequencing to analyze the steady-state levels and sizes of mitochondrial RNAs, their degradation products, as well as their processing intermediates by capturing both 5' and 3' ends of transcripts.

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240–300 passages, corresponding to approximately 400–600 generations or almost 3 years. The sequences show signs of extensive parallel evolution—the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

Analysis of the longest evolutionary experiment with HIV-1 to-date reveals continuous viral adaptation over several years. The authors quantify the environment-specific mutations that arise and determine the fraction of mutations that co-occur with significantly different frequencies than expected by chance.

Biased Pol II fidelity contributes to conservation of functional domains in the Potato spindle tuber viroid genome

Accurate calculation of mutation rates for viruses and viroids is necessary for evolutionary studies and to evaluate adaptation potential. However, estimation of in vivo mutation rates is complicated by selection, which leads to loss or proliferation of certain mutations. To minimize this concern, lethal mutations, including nonsense and non-synonymous mutations, have been used to determine mutation rates for several viruses and viroids, including Potato spindle tuber viroid (PSTVd). However, this approach has limitations, including focus on a relatively small number of genome sites and the possibility that mutations may not actually be lethal or may be maintained by wild type individuals. To avoid selection bias altogether, we sequenced minus-strand PSTVd dimers from concatemeric replication intermediates. The underlying rationale is that mutations found in only one of the monomers were likely generated de novo during RNA polymerase II (Pol II) transcription of the circular plus-strand RNA genome. This approach yielded an apparent Pol II error rate of ~1/1837 nucleotides per transcription cycle, and an estimated mutation rate of ~1/919 nucleotides for a single replication cycle. Remarkably, de novo mutations were nearly absent from the most conserved, replication-critical regions of the PSTVd genome, suggesting that sequence conservation is a consequence of both essential function and template optimization for greater Pol II fidelity. Such biased fidelity may constitute a novel strategy to ensure population success while allowing abundant sampling of sequence space in other genome regions. Comparison with variants in progeny populations derived from a cloned, wild type PSTVd master sequence revealed that most de novo mutations were lost through selection.

Polymerase errors are the major source of variation in virus and viroid genomes, and as a consequence polymerase error rates are major determinants of adaptation potential. Accurate calculation of in vivo mutation rates is complicated by selection. To circumvent this issue, dimeric PSTVd minus-strand replication intermediates generated in vivo by host RNA polymerase II (Pol II) were sequenced to identify de novo mutations. This analysis revealed a very high error rate for Pol II transcribing genomic PSTVd RNA, leading to an extremely high mutation rate. Remarkably, however, de novo mutations were rare in the most highly conserved, replication-critical genome regions, suggesting these sequences are selected for both function and enhanced transcription fidelity. This biased fidelity may reveal a novel strategy to ensure population survival while maximizing adaptation potential. Further, comparison of mutations identified by minus-strand dimer sequencing with mutations observed in progeny variants derived from wild type PSTVd showed that most de novo mutations were lost through selection.

The Early Evolution of Oral Poliovirus Vaccine Is Shaped by Strong Positive Selection and Tight Transmission Bottlenecks

The emergence of circulating vaccine-derived polioviruses through evolution of the oral polio vaccine (OPV) poses a significant obstacle to polio eradication. Understanding the early genetic changes that occur as OPV evolves and transmits is important for preventing future outbreaks. Here, we use deep sequencing to define the evolutionary trajectories of type 2 OPV in a vaccine trial. By sequencing 497 longitudinal stool samples from 271 OPV2 recipients and household contacts, we were able to examine the extent of convergent evolution in vaccinated individuals and the amount of viral diversity that is transmitted. In addition to rapid reversion of key attenuating mutations, we identify strong selection at 19 sites across the genome. We find that a tight transmission bottleneck limits the onward transmission of these early adaptive mutations. Our results highlight the distinct evolutionary dynamics of live attenuated virus vaccines and have important implications for the success of next-generation OPV.

A Bayesian Framework for Inferring the Influence of Sequence Context on Point Mutations

The probability of point mutations is expected to be highly influenced by the flanking nucleotides that surround them, known as the sequence context. This phenomenon may be mainly attributed to the enzyme that modifies or mutates the genetic material, because most enzymes tend to have specific sequence contexts that dictate their activity. Here, we develop a statistical model that allows for the detection and evaluation of the effects of different sequence contexts on mutation rates from deep population sequencing data. This task is computationally challenging, as the complexity of the model increases exponentially as the context size increases. We established our novel Bayesian method based on sparse model selection methods, with the leading assumption that the number of actual sequence contexts that directly influence mutation rates is minuscule compared with the number of possible sequence contexts. We show that our method is highly accurate on simulated data using pentanucleotide contexts, even when accounting for noisy data. We next analyze empirical population sequencing data from polioviruses and HIV-1 and detect a significant enrichment in sequence contexts associated with deamination by the cellular deaminases ADAR 1/2 and APOBEC3G, respectively. In the current era, where next-generation sequencing data are highly abundant, our approach can be used on any population sequencing data to reveal context-dependent base alterations and may assist in the discovery of novel mutable sites or editing sites.

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

Genetic diversity is the fuel of evolution and facilitates adaptation to novel environments. However, our understanding of what drives differences in the genetic diversity during the early stages of viral infection is somewhat limited. Here, we use ultra-deep sequencing to interrogate 43 clinical samples taken from early infections of the human-infecting viruses HIV, RSV and CMV. Hundreds to thousands of virus templates were sequenced per sample, allowing us to reveal dramatic differences in within-host genetic diversity among virus populations. We found that increased diversity was mostly driven by presence of multiple divergent genotypes in HIV and CMV samples, which we suggest reflect multiple transmitted/founder viruses. Conversely, we detected an abundance of low frequency hyper-edited genomes in RSV samples, presumably reflecting defective virus genomes (DVGs). We suggest that RSV is characterized by higher levels of cellular co-infection, which allow for complementation and hence elevated levels of DVGs.

The few days or weeks following infection with a virus, termed acute infection, are critical for virus establishment. Here we sought to characterize what leads to differences in the genetic diversity of different viruses sampled during acute infection. We performed ultra-deep sequencing of hundreds to thousands viral genomes from forty-three samples spanning three pathogenic human viruses: HIV, RSV and CMV. We found major differences in the genetic diversity of these different viruses, and in different patients infected with the same virus. We investigated the factors responsible for these differences. We found that the DNA virus CMV was less diverse, most likely since it has a lower mutation rate than the RNA viruses HIV and RSV. We also found that the samples with the highest genetic diversity, which included one CMV sample and two HIV samples, bore evidence for multiple genotype infection. In other words, patients from whom these samples were taken were infected with two different “strains” of the virus. Finally, we also found evidence that viral genomes of HIV, and in particular RSV, are edited by the innate immune system of the host, leading to the presence of defective virus genomes.

An Evolutionary Model of Progression to AIDS

The time to the onset of AIDS symptoms in an HIV infected individual is known to correlate inversely with viremia and the level of immune activation. The correlation exists against the background of strong individual fluctuations demonstrating the existence of hidden variables depending on patient and virus parameters. At the moment, prognosis of the time to AIDS based on patient parameters is not possible. In addition, it is of paramount importance to understand the reason of progression to AIDS in untreated patients to be able to learn to control it by means other than anti-retroviral therapy. Here we develop a mechanistic mathematical model to predict the speed of progression to AIDS in individual untreated patients and patients treated with suboptimal therapy, based on a single-time measurement of several virological and immunological parameters. We show that the gradual increase in virus fitness during a chronic infection causes slow gradual depletion of CD4 T cells. Using the existing evolution models of HIV, we obtain general expressions predicting the time to the onset of AIDS symptoms in terms of the patient parameters, for low-viremia and high-viremia patients separately. We show that the evolution model of AIDS fits the existing data on virus-time correlations better than the alternative model of the deregulation of homeostatic response.

Extended interaction networks with HCV protease NS3-4A substrates explain the lack of adaptive capability against protease inhibitors

Inhibitors against the NS3-4A protease of hepatitis C virus (HCV) have proven to be useful drugs in the treatment of HCV infection. Although variants have been identified with mutations that confer resistance to these inhibitors, the mutations do not restore replicative fitness and no secondary mutations that rescue fitness have been found. To gain insight into the molecular mechanisms underlying the lack of fitness compensation, we screened known resistance mutations in infectious HCV cell culture with different genomic backgrounds. We observed that the Q41R mutation of NS3-4A efficiently rescues the replicative fitness in cell culture for virus variants containing mutations at NS3-Asp168 To understand how the Q41R mutation rescues activity, we performed protease activity assays complemented by molecular dynamics simulations, which showed that protease-peptide interactions far outside the targeted peptide cleavage sites mediate substrate recognition by NS3-4A and support protease cleavage kinetics. These interactions shed new light on the mechanisms by which NS3-4A cleaves its substrates, viral polyproteins and a prime cellular antiviral adaptor protein, the mitochondrial antiviral signaling protein MAVS. Peptide binding is mediated by an extended hydrogen-bond network in NS3-4A that was effectively optimized for protease-MAVS binding in Asp168 variants with rescued replicative fitness from NS3-Q41R. In the protease harboring NS3-Q41R, the N-terminal cleavage products of MAVS retained high affinity to the active site, rendering the protease susceptible for potential product inhibition. Our findings reveal delicately balanced protease-peptide interactions in viral replication and immune escape that likely restrict the protease adaptive capability and narrow the virus evolutionary space.

The mechanism and detection of alternative splicing events in circular RNAs

Circular RNAs (circRNAs) are considered as functional biomolecules with tissue/development-specific expression patterns. Generally, a single gene may generate multiple circRNA variants by alternative splicing, which contain different combinations of exons and/or introns. Due to the low abundance of circRNAs as well as overlapped with their linear counterparts, circRNA enrichment protocol is needed prior to sequencing. Compared with numerous algorithms, which use back-splicing reads for detection and functional characterization of circRNAs, original bioinformatic analyzing tools have been developed to large-scale determination of full-length circRNAs and accurate quantification. This review provides insights into the complexity of circRNA biogenesis and surveys the recent progresses in the experimental and bioinformatic methodologies that focus on accurately full-length circRNAs identification.

Species-Specific Evolution of Ebola Virus during Replication in Human and Bat Cells

Ebola virus (EBOV) causes a severe, often fatal disease in humans and nonhuman primates. Within the past decade, EBOV has caused two large and difficult-to-control outbreaks, one of which recently ended in the Democratic Republic of the Congo. Bats are the likely reservoir of EBOV, but little is known of their relationship with the virus. We perform serial passages of EBOV in human and bat cells and use circular sequencing to compare the short-term evolution of the virus. Virus populations passaged in bat cells have sequence markers indicative of host RNA editing enzyme activity, including evidence for ADAR editing of the EBOV glycoprotein. Multiple regions in the EBOV genome appear to have undergone adaptive evolution when passaged in bat and human cells. Individual mutated viruses are rescued and characterized. Our results provide insight into the host species-specific evolution of EBOV and highlight the adaptive flexibility of the virus.

Applying next-generation sequencing to unravel the mutational landscape in viral quasispecies

Next-generation sequencing (NGS) has revolutionized the scale and depth of biomedical sciences. Because of its unique ability for the detection of sub-clonal variants within genetically diverse populations, NGS has been successfully applied to analyze and quantify the exceptionally-high diversity within viral quasispecies, and many low-frequency drug- or vaccine-resistant mutations of therapeutic importance have been discovered. Although many works have intensively discussed the latest NGS approaches and applications in general, none of them has focused on applying NGS in viral quasispecies studies, mostly due to the limited ability of current NGS technologies to accurately detect and quantify rare viral variants. Here, we summarize several error-correction strategies that have been developed to enhance the detection accuracy of minority variants. We also discuss critical considerations for preparing a sequencing library from viral RNAs and for analyzing NGS data to unravel the mutational landscape.

Identifying Transcription Error-Enriched Genomic Loci Using Nuclear Run-on Circular-Sequencing Coupled with Background Error Modeling

RNA polymerase transcribes certain genomic loci with higher errors rates. These transcription error-enriched genomic loci (TEELs) have implications in disease. Current deep-sequencing methods cannot distinguish TEELs from post-transcriptional modifications, stochastic transcription errors, and technical noise, impeding efforts to elucidate the mechanisms linking TEELs to disease. Here, we describe background error model-coupled precision nuclear run-on circular-sequencing (EmPC-seq) to discern genomic regions enriched for transcription misincorporations. EmPC-seq innovatively combines a nuclear run-on assay for capturing nascent RNA before post-transcriptional modifications, a circular-sequencing step that sequences the same nascent RNA molecules multiple times to improve accuracy, and a statistical model for distinguishing error-enriched regions among stochastic polymerase errors. Applying EmPC-seq to the ribosomal RNA transcriptome, we show that TEELs of RNA polymerase I are not randomly distributed but clustered together, with higher error frequencies at nascent transcript 3' ends. Our study establishes a reliable method of identifying TEELs with nucleotide precision, which can help elucidate their molecular origins.

Long-term context-dependent genetic adaptation of the viral genetic cloud

MOTIVATION

RNA viruses generate a cloud of genetic variants within each host. This cloud contains high-frequency genotypes, and many rare variants. The dynamics of these variants is crucial to understand viral evolution and their effect on their host.

RESULTS

We use an experimental evolution system to show that the genetic cloud surrounding the Coxsackie virus master sequence slowly, but steadily, evolves over hundreds of generations. This movement is determined by strong context-dependent mutations, where the frequency and type of mutations are affected by neighboring positions, even in silent mutations. This context-dependent mutation pattern serves as a spearhead for the viral population's movement within the adaptive landscape and affects which new dominant variants will emerge. The non-local mutation patterns affect the mutated dinucleotide distribution, and eventually lead to a non-uniform dinucleotide distribution in the main viral sequence. We tested these results on other RNA viruses with similar conclusions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Universally high transcript error rates in bacteria

Errors can occur at any level during the replication and transcription of genetic information. Genetic mutations derived mainly from replication errors have been extensively studied. However, fundamental details of transcript errors, such as their rate, molecular spectrum, and functional effects, remain largely unknown. To globally identify transcript errors, we applied an adapted rolling-circle sequencing approach to Escherichia coli, Bacillus subtilis, Agrobacterium tumefaciens, and Mesoplasma florum, revealing transcript-error rates 3 to 4 orders of magnitude higher than the corresponding genetic mutation rates. The majority of detected errors would result in amino-acid changes, if translated. With errors identified from 9929 loci, the molecular spectrum and distribution of errors were uncovered in great detail. A G→A substitution bias was observed in M. florum, which apparently has an error-prone RNA polymerase. Surprisingly, an increased frequency of nonsense errors towards the 3' end of mRNAs was observed, suggesting a Nonsense-Mediated Decay-like quality-control mechanism in prokaryotes.

Engineering the Live-Attenuated Polio Vaccine to Prevent Reversion to Virulence

The live-attenuated oral poliovirus vaccine (OPV or Sabin vaccine) replicates in gut-associated tissues, eliciting mucosa and systemic immunity. OPV protects from disease and limits poliovirus spread. Accordingly, vaccination with OPV is the primary strategy used to end the circulation of all polioviruses. However, the ability of OPV to regain replication fitness and establish new epidemics represents a significant risk of polio re-emergence should immunization cease. Here, we report the development of a poliovirus type 2 vaccine strain (nOPV2) that is genetically more stable and less likely to regain virulence than the original Sabin2 strain. We introduced modifications within at the 5' untranslated region of the Sabin2 genome to stabilize attenuation determinants, 2C coding region to prevent recombination, and 3D polymerase to limit viral adaptability. Prior work established that nOPV2 is immunogenic in preclinical and clinical studies, and thus may enable complete poliovirus eradication.

The honeybee (Apis mellifera) developmental state shapes the genetic composition of the deformed wing virus-A quasispecies during serial transmission

The main biological threat to the western honeybee (Apis mellifera) is the parasitic mite Varroa destructor, largely because it vectors lethal epidemics of honeybee viruses that, in the absence of this mite, are relatively innocuous. The severe pathology is a direct consequence of excessive virus titres caused by this novel transmission route. However, little is known about how the virus adapts genetically during transmission and whether this influences the pathology. Here, we show that upon injection into honeybee pupae, the deformed wing virus type-A (DWV-A) quasispecies undergoes a rapid, extensive expansion of its sequence space, followed by strong negative selection towards a uniform, common shape by the time the pupae have completed their development, with no difference between symptomatic and asymptomatic adults in either DWV titre or genetic composition. This suggests that the physiological and molecular environment during pupal development has a strong, conservative influence on shaping the DWV-A quasispecies in emerging adults. There was furthermore no evidence of any progressive adaptation of the DWV-A quasispecies to serial intra-abdominal injection, simulating mite transmission, despite the generation of ample variation immediately following each transmission, suggesting that the virus either had already adapted to transmission by injection, or was unaffected by it.

Comparative analysis of CpG islands in equine infectious anemia virus strains

Increasing evidence suggests that DNA methylation has key roles in the replication of retroviruses, including lentiviruses, and pathogenesis of diseases. However, the precise characteristics of CpG islands are not known for many retroviruses. In this study, we compared the distribution of CpG islands among strains of equine infectious anemia virus (EIAV), a lentivirus in the family Retroviridae and a model for HIV research. We identified CpG islands in 32 full-length EIAV genomic sequences obtained from the GenBank database using MethPrimer. Only one CpG island, from 100 to 120 bp, was identified in the genomes of EIAV strains DV10, DLV3-A, and DLV5-10 from China, V26 and V70 from Japan, and IRE H3, IRE F2, IRE F3, and IRE F4 from Ireland. Importantly, the CpG island was located within the Rev gene, which is required for the expression of viral cis-acting elements and the production of new virions. These results suggest that the distribution, length, and genetic properties of CpG islands differ among EIAV strains. Future research should focus on the biological significance of this CpG island within rev to improve our understanding of the precise roles of CpG islands in epigenetic regulation in the species.

Into the Deep (Sequence) of the Foot-and-Mouth Disease Virus Gene Pool: Bottlenecks and Adaptation during Infection in Naïve and Vaccinated Cattle

Foot-and-mouth disease virus (FMDV) infects hosts as a population of closely related viruses referred to as a quasispecies. The behavior of this quasispecies has not been described in detail in natural host species. In this study, virus samples collected from vaccinated and non-vaccinated cattle up to 35 days post-experimental infection with FMDV A24-Cruzeiro were analyzed by deep-sequencing. Vaccination induced significant differences compared to viruses from non-vaccinated cattle in substitution rates, entropy, and evidence for adaptation. Genomic variation detected during early infection reflected the diversity inherited from the source virus (inoculum), whereas by 12 days post infection, dominant viruses were defined by newly acquired mutations. Mutations conferring recognized fitness gain occurred and were associated with selective sweeps. Persistent infections always included multiple FMDV subpopulations, suggesting distinct foci of infection within the nasopharyngeal mucosa. Subclinical infection in vaccinated cattle included very early bottlenecks associated with reduced diversity within virus populations. Viruses from both animal cohorts contained putative antigenic escape mutations. However, these mutations occurred during later stages of infection, at which time transmission is less likely to occur. This study improves upon previously published work by analyzing deep sequences of samples, allowing for detailed characterization of FMDV populations over time within multiple hosts.

Sequencing of serially passaged measles virus affirms its genomic stability and reveals a nonrandom distribution of consensus mutations

Oncolytic virotherapy is an emerging treatment option for numerous cancers, with several virus families currently being evaluated in clinical trials. More specifically, vaccine-strain measles virus has arisen as a promising candidate for the treatment of different tumour types in several early clinical trials. Replicating viruses, and especially RNA viruses without proofreading polymerases, can rapidly adapt to varying environments by selecting quasispecies with advantageous genetic mutations. Subsequently, these genetic alterations could potentially weaken the safety profile of virotherapy. In this study, we demonstrate that, following an extended period of virus replication in producer or cancer cell lines, the quasispecies consensus sequence of vaccine strain-derived measles virus accrues a remarkably small number of mutations throughout the nonsegmented negative-stranded RNA genome. Interestingly, we detected a nonrandom distribution of genetic alterations within the genome, with an overall decreasing frequency of mutations from the 3' genome start to its 5' end. Comparing the serially passaged viruses to the parental virus on producer cells, we found that the acquired consensus mutations did not drastically change viral replication kinetics or cytolytic potency. Collectively, our data corroborate the genomic stability and excellent safety profile of oncolytic measles virus, thus supporting its continued development and clinical translation as a promising viro-immunotherapeutic.

Wide spectrum and high frequency of genomic structural variation, including transposable elements, in large double-stranded DNA viruses

Our knowledge of the diversity and frequency of genomic structural variation segregating in populations of large double-stranded (ds) DNA viruses is limited. Here, we sequenced the genome of a baculovirus (Autographa californica multiple nucleopolyhedrovirus [AcMNPV]) purified from beet armyworm (Spodoptera exigua) larvae at depths >195,000× using both short- (Illumina) and long-read (PacBio) technologies. Using a pipeline relying on hierarchical clustering of structural variants (SVs) detected in individual short- and long-reads by six variant callers, we identified a total of 1,141 SVs in AcMNPV, including 464 deletions, 443 inversions, 160 duplications, and 74 insertions. These variants are considered robust and unlikely to result from technical artifacts because they were independently detected in at least three long reads as well as at least three short reads. SVs are distributed along the entire AcMNPV genome and may involve large genomic regions (30,496 bp on average). We show that no less than 39.9 per cent of genomes carry at least one SV in AcMNPV populations, that the vast majority of SVs (75%) segregate at very low frequency (<0.01%) and that very few SVs persist after ten replication cycles, consistent with a negative impact of most SVs on AcMNPV fitness. Using short-read sequencing datasets, we then show that populations of two iridoviruses and one herpesvirus are also full of SVs, as they contain between 426 and 1,102 SVs carried by 52.4–80.1 per cent of genomes. Finally, AcMNPV long reads allowed us to identify 1,757 transposable elements (TEs) insertions, 895 of which are truncated and occur at one extremity of the reads. This further supports the role of baculoviruses as possible vectors of horizontal transfer of TEs. Altogether, we found that SVs, which evolve mostly under rapid dynamics of gain and loss in viral populations, represent an important feature in the biology of large dsDNA viruses.

Virus population dynamics examined with experimental model systems

Experimental evolution permits exploring the effect of controlled environmental variables in virus evolution. Several designs in cell culture and in vivo have established basic concepts that can assist in the interpretation of evolutionary events in the field. Important information has come from cytolytic and persistent infections in cell culture that have unveiled the power of virus-cell coevolution in virus and cell diversification. Equally informative are comparisons of the response of viral populations when subjected to different passage régimens. In particular, plaque-to-plaque transfers in cell culture have revealed unusual genotypes and phenotypes that populate minority layers of viral quasispecies. Some of these viruses display properties that contradict features established in virology textbooks. Several hypotheses and principles of population genetics have found experimental confirmation in experimental designs with viruses. The possibilities of using experimental evolution to understand virus behavior are still largely unexploited.

In vitro evolution of herpes simplex virus 1 (HSV-1) reveals selection for syncytia and other minor variants in cell culture

The large dsDNA virus herpes simplex virus 1 (HSV-1) is considered to be genetically stable, yet it can rapidly evolve in response to strong selective pressures such as antiviral treatment. Deep sequencing has revealed that clinical and laboratory isolates of this virus exist as populations that contain a mixture of minor alleles or variants, similar to many RNA viruses. The classic virology approach of plaque purifying virus creates a genetically homogenous population, but it is not clear how closely this represents the mixed virus populations found in nature. We sought to study the evolution of mixed versus highly purified HSV-1 populations in controlled cell culture conditions, to examine the impact of this genetic diversity on evolution. We found that a mixed population of HSV-1 acquired more genetic diversity and underwent a more dramatic phenotypic shift than a plaque-purified population, producing a viral population that was almost entirely syncytial after just ten passages. At the genomic level, adaptation and genetic diversification occurred at the level of minor alleles or variants in the viral population. Certain genetic variants in the mixed viral population appeared to be positively selected in cell culture, and this shift was also observed in clinical samples during their first passages in vitro. In contrast, the plaque-purified viral population did not appear to change substantially in phenotype or overall quantity of minor allele diversity. These data indicate that HSV-1 is capable of evolving rapidly in a given environment, and that this evolution is facilitated by diversity in the viral population.

Increase in the genetic polymorphism of varicella-zoster virus after passaging in in vitro cell culture

Primary infections with the varicella-zoster virus (VZV) result in varicella, while latent reactivation leads to herpes zoster. Both varicella and zoster can be prevented by live attenuated vaccines. There have been reports suggesting that both clinical VZV strains and those in vaccine preparations are genetically polymorphic, containing mixtures of both wild-type and vaccine-type sequences at certain vaccine-specific sites. In this study, the genetic polymorphism of the VZV genome was examined by analyzing the frequencies of minor alleles at each nucleotide position. Next-generation sequencing of the clinical VZV strain YC02 passaged in an in vitro cell culture was used to identify genetically polymorphic sites (GPS), where the minor allele frequency (MAF) exceeded 5%. The number of GPS increased by 7.3-fold at high passages (p100) when compared to low passages (p17), although the average MAF remained similar. GPS were found in 6 open reading frames (ORFs) in p17, 35, and 54 ORFs in p60 and p100, respectively. GPS were found more frequently in the dispensable gene group than the essential gene group, but the average MAF was greater in the essential gene group. The most common two major/minor base pairs were A/g and T/c. GPS were found in all three passages at 16 positions, all located in the reiterated (R) region. The population diversity as measured by Shannon entropy increased in p60 and p100. However, the entropy remained unchanged in the R regions.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Identification of Loci Associated with Enhanced Virulence in Spodoptera litura Nucleopolyhedrovirus Isolates Using Deep Sequencing

Spodoptera litura is an emerging pest insect in cotton and arable crops in Central Asia. To explore the possibility of using baculoviruses as biological control agents instead of chemical pesticides, in a previous study we characterized a number of S. litura nucleopolyhedrovirus (SpltNPV) isolates from Pakistan. We found significant differences in speed of kill, an important property of a biological control agent. Here we set out to understand the genetic basis of these differences in speed of kill, by comparing the genome of the fast-killing SpltNPV-Pak-TAX1 isolate with that of the slow-killing SpltNPV-Pak-BNG isolate. These two isolates and the SpltNPV-G2 reference strain from China were deep sequenced with Illumina. As expected, the two Pakistani isolates were closely related with >99% sequence identity, whereas the Chinese isolate was more distantly related. We identified two loci that may be associated with the fast action of the SpltNPV-Pak-TAX1 isolate. First, an analysis of rates of synonymous and non-synonymous mutations identified neutral to positive selection on open reading frame (ORF) 122, encoding a viral fibroblast growth factor (vFGF) that is known to affect virulence in other baculoviruses. Second, the homologous repeat region hr17, a putative enhancer of transcription and origin of replication, is absent in SpltNPV-Pak-TAX1 suggesting it may also affect virulence. Additionally, we found there is little genetic variation within both Pakistani isolates, and we identified four genes under positive selection in both isolates that may have played a role in adaptation of SpltNPV to conditions in Central Asia. Our results contribute to the understanding of the enhanced activity of SpltNPV-Pak-TAX1, and may help to select better SpltNPV isolates for the control of S. litura in Pakistan and elsewhere.

Panicum Mosaic Virus and Its Satellites Acquire RNA Modifications Associated with Host-Mediated Antiviral Degradation

Positive-sense RNA viruses in the Tombusviridae family have genomes lacking a 5' cap structure and prototypical 3' polyadenylation sequence. Instead, these viruses utilize an extensive network of intramolecular RNA-RNA interactions to direct viral replication and gene expression. Here we demonstrate that the genomic RNAs of Panicum mosaic virus (PMV) and its satellites undergo sequence modifications at their 3' ends upon infection of host cells. Changes to the viral and subviral genomes arise de novo within Brachypodium distachyon (herein called Brachypodium) and proso millet, two alternative hosts of PMV, and exist in the infections of a native host, St. Augustinegrass. These modifications are defined by polyadenylation [poly(A)] events and significant truncations of the helper virus 3' untranslated region-a region containing satellite RNA recombination motifs and conserved viral translational enhancer elements. The genomes of PMV and its satellite virus (SPMV) were reconstructed from multiple poly(A)-selected Brachypodium transcriptome data sets. Moreover, the polyadenylated forms of PMV and SPMV RNAs copurify with their respective mature icosahedral virions. The changes to viral and subviral genomes upon infection are discussed in the context of a previously understudied poly(A)-mediated antiviral RNA degradation pathway and the potential impact on virus evolution.IMPORTANCE The genomes of positive-sense RNA viruses have an intrinsic capacity to serve directly as mRNAs upon viral entry into a host cell. These RNAs often lack a 5' cap structure and 3' polyadenylation sequence, requiring unconventional strategies for cap-independent translation and subversion of the cellular RNA degradation machinery. For tombusviruses, critical translational regulatory elements are encoded within the 3' untranslated region of the viral genomes. Here we describe RNA modifications occurring within the genomes of Panicum mosaic virus (PMV), a prototypical tombusvirus, and its satellite agents (i.e., satellite virus and noncoding satellite RNAs), all of which depend on the PMV-encoded RNA polymerase for replication. The atypical RNAs are defined by terminal polyadenylation and truncation events within the 3' untranslated region of the PMV genome. These modifications are reminiscent of host-mediated RNA degradation strategies and likely represent a previously underappreciated defense mechanism against invasive nucleic acids.

A new full-length circular DNA sequencing method for viral-sized genomes reveals that RNAi transgenic plants provoke a shift in geminivirus populations in the field

We present a new method, CIDER-Seq (Circular DNA Enrichment sequencing) for the unbiased enrichment and long-read sequencing of viral-sized circular DNA molecules. We used CIDER-Seq to produce single-read full-length virus genomes for the first time. CIDER-Seq combines PCR-free virus enrichment with Single Molecule Real Time sequencing and a new sequence de-concatenation algorithm. We apply our technique to produce >1200 full-length, highly accurate geminivirus genomes from RNAi-transgenic and control plants in a field trial in Kenya. Using CIDER-Seq we can demonstrate for the first time that the expression of antiviral double-stranded RNA (dsRNA) in transgenic plants causes a consistent shift in virus populations towards species sharing low homology to the transgene derived dsRNA. Our method and its application in an economically important crop plant opens new possibilities in periodic virus sequence surveillance and accurate profiling of diverse circular DNA elements.

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

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

SiNPle: Fast and Sensitive Variant Calling for Deep Sequencing Data

Current high-throughput sequencing technologies can generate sequence data and provide information on the genetic composition of samples at very high coverage. Deep sequencing approaches enable the detection of rare variants in heterogeneous samples, such as viral quasi-species, but also have the undesired effect of amplifying sequencing errors and artefacts. Distinguishing real variants from such noise is not straightforward. Variant callers that can handle pooled samples can be in trouble at extremely high read depths, while at lower depths sensitivity is often sacrificed to specificity. In this paper, we propose SiNPle (Simplified Inference of Novel Polymorphisms from Large coveragE), a fast and effective software for variant calling. SiNPle is based on a simplified Bayesian approach to compute the posterior probability that a variant is not generated by sequencing errors or PCR artefacts. The Bayesian model takes into consideration individual base qualities as well as their distribution, the baseline error rates during both the sequencing and the PCR stage, the prior distribution of variant frequencies and their strandedness. Our approach leads to an approximate but extremely fast computation of posterior probabilities even for very high coverage data, since the expression for the posterior distribution is a simple analytical formula in terms of summary statistics for the variants appearing at each site in the genome. These statistics can be used to filter out putative SNPs and indels according to the required level of sensitivity. We tested SiNPle on several simulated and real-life viral datasets to show that it is faster and more sensitive than existing methods. The source code for SiNPle is freely available to download and compile, or as a Conda/Bioconda package.