Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline (http://sourceforge.net/projects/quasr) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro.     

SARS-CoV-2 biology and variants: anticipation of viral evolution and what needs to be done

The global propagation of SARS-CoV-2 and the detection of a large number of variants, some of which have replaced the original clade to become dominant, underscores the fact that the virus is actively exploring its evolutionary space. The longer high levels of viral multiplication occur – permitted by high levels of transmission –, the more the virus can adapt to the human host and find ways to success. The third wave of the COVID-19 pandemic is starting in different parts of the world, emphasizing that transmission containment measures that are being imposed are not adequate. Part of the consideration in determining containment measures is the rationale that vaccination will soon stop transmission and allow a return to normality. However, vaccines themselves represent a selection pressure for evolution of vaccine-resistant variants, so the coupling of a policy of permitting high levels of transmission/virus multiplication during vaccine roll-out with the expectation that vaccines will deal with the pandemic, is unrealistic. In the absence of effective antivirals, it is not improbable that SARS-CoV-2 infection prophylaxis will involve an annual vaccination campaign against ‘dominant’ viral variants, similar to influenza prophylaxis. Living with COVID-19 will be an issue of SARS-CoV-2 variants and evolution. It is therefore crucial to understand how SARS-CoV-2 evolves and what constrains its evolution, in order to anticipate the variants that will emerge. Thus far, the focus has been on the receptor-binding spike protein, but the virus is complex, encoding 26 proteins which interact with a large number of host factors, so the possibilities for evolution are manifold and not predictable a priori. However, if we are to mount the best defence against COVID-19, we must mount it against the variants, and to do this, we must have knowledge about the evolutionary possibilities of the virus. In addition to the generic cellular interactions of the virus, there are extensive polymorphisms in humans (e.g. Lewis, HLA, etc.), some distributed within most or all populations, some restricted to specific ethnic populations and these variations pose additional opportunities for/constraints on viral evolution. We now have the wherewithal – viral genome sequencing, protein structure determination/modelling, protein interaction analysis – to functionally characterize viral variants, but access to comprehensive genome data is extremely uneven. Yet, to develop an understanding of the impacts of such evolution on transmission and disease, we must link it to transmission (viral epidemiology) and disease data (patient clinical data), and the population granularities of these. In this editorial, we explore key facets of viral biology and the influence of relevant aspects of human polymorphisms, human behaviour, geography and climate and, based on this, derive a series of recommendations to monitor viral evolution and predict the types of variants that are likely to arise.     

Phylogeny of <Emphasis Type="Italic">Banana Streak Virus</Emphasis> Reveals Recent and Repetitive Endogenization in the Genome of Its Banana Host (<Emphasis Type="Italic">Musa</Emphasis> sp.)

Banana streak virus (BSV) is a plant dsDNA pararetrovirus (family Caulimoviridae, genus badnavirus). Although integration is not an essential step in the BSV replication cycle, the nuclear genome of banana (Musa sp.) contains BSV endogenous pararetrovirus sequences (BSV EPRVs). Some BSV EPRVs are infectious by reconstituting a functional viral genome. Recent studies revealed a large molecular diversity of episomal BSV viruses (i.e., nonintegrated) while others focused on BSV EPRV sequences only. In this study, the evolutionary history of badnavirus integration in banana was inferred from phylogenetic relationships between BSV and BSV EPRVs. The relative evolution rates and selective pressures (d_N/d_S ratio) were also compared between endogenous and episomal viral sequences. At least 27 recent independent integration events occurred after the divergence of three banana species, indicating that viral integration is a recent and frequent phenomenon. Relaxation of selective pressure on badnaviral sequences that experienced neutral evolution after integration in the plant genome was recorded. Additionally, a significant decrease (35%) in the EPRV evolution rate was observed compared to BSV, reflecting the difference in the evolution rate between episomal dsDNA viruses and plant genome. The comparison of our results with the evolution rate of the Musa genome and other reverse-transcribing viruses suggests that EPRVs play an active role in episomal BSV diversity and evolution.     

CoVrimer: A tool for aligning SARS-CoV-2 primer sequences and selection of conserved/degenerate primers

As mutations in SARS-CoV-2 virus accumulate rapidly, novel primers that amplify this virus sensitively and specifically are in demand. We have developed a webserver named CoVrimer by which users can search for and align existing or newly designed conserved/degenerate primer pair sequences against the viral genome and assess the mutation load of both primers and amplicons. CoVrimer uses mutation data obtained from an online platform established by NGDC-CNCB (12 May 2021) to identify genomic regions, either conserved or with low levels of mutations, from which potential primer pairs are designed and provided to the user for filtering based on generalized and SARS-CoV-2 specific parameters. Alignments of primers and probes can be visualized with respect to the reference genome, indicating variant details and the level of conservation. Consequently, CoVrimer is likely to help researchers with the challenges posed by viral evolution and is freely available at http://konulabapps.bilkent.edu.tr:3838/CoVrimer/.     

Role of spatial patterning of N-protein interactions in SARS-CoV-2 genome packaging

Viruses must efficiently and specifically package their genomes while excluding cellular nucleic acids and viral subgenomic fragments. Some viruses use specific packaging signals, which are conserved sequence or structure motifs present only in the full-length genome. Recent work has shown that viral proteins important for packaging can undergo liquid-liquid phase separation (LLPS), in which one or two viral nucleic acid binding proteins condense with the genome. The compositional simplicity of viral components lends itself well to theoretical modeling compared with more complex cellular organelles. Viral LLPS can be limited to one or two viral proteins and a single genome that is enriched in LLPS-promoting features. In our previous study, we observed that LLPS-promoting sequences of severe acute respiratory syndrome coronavirus 2 are located at the 5′ and 3′ ends of the genome, whereas the middle of the genome is predicted to consist mostly of solubilizing elements. Is this arrangement sufficient to drive single genome packaging, genome compaction, and genome cyclization? We addressed these questions using a coarse-grained polymer model, LASSI, to study the LLPS of nucleocapsid protein with RNA sequences that either promote LLPS or solubilization. With respect to genome cyclization, we find the most optimal arrangement restricts LLPS-promoting elements to the 5′ and 3′ ends of the genome, consistent with the native spatial patterning. Genome compaction is enhanced by clustered LLPS-promoting binding sites, whereas single genome packaging is most efficient when binding sites are distributed throughout the genome. These results suggest that many and variably positioned LLPS-promoting signals can support packaging in the absence of a singular packaging signal which argues against necessity of such a feature. We hypothesize that this model should be generalizable to multiple viruses as well as cellular organelles such as paraspeckles, which enrich specific long RNA sequences in a defined arrangement.     

The Genomic Signature of Human Rhinoviruses A,B and C

Human rhinoviruses are single stranded positive sense RNA viruses that are presented in more than 50% of acute upper respiratory tract infections. Despite extensive studies on the genetic diversity of the virus, little is known about the forces driving it. In order to explain this diversity, many research groups have focused on protein sequence requirements for viable, functional and transmissible virus but have missed out an important aspect of viral evolution such as the genomic ontology of the virus. This study presents for the first time the genomic signature of 111 fully sequenced HRV strains from all three groups HRV-A, HRV-B and HRV-C. We observed an HRV genome tendency to eliminate CpG and UpA dinucleotides, coupling with over-representation of UpG and CpA. We propose a specific mechanism which describes how rapid changes in the HRV genomic sequence can take place under the strict control of conservation of the polypeptide backbone. Moreover, the distribution of the observed under- and over-represented dinucleotides along the HRV genome is presented. Distance matrice tables based on CpG and UpA odds ratios were constructed and viewed as heatmaps and distance trees. None of the suppressions can be attributed to codon usage or in RNA secondary structure requirements. Since viral recognition is dependent on RNA motifs rich in CpG and UpA, it is possible that the overall described genome evolution mechanism acts in order to protect the virus from host recognition.     

Chromosome-protein interactions in polyomavirus virions

In this work, we sought to determine whether the components of the murine polyomavirus capsid establish specific interactions with the minichromosome encapsidated into the mature viral particles by using the cis-diamminedichloroplatinum(II) cross-linking reagent. Our data indicated that VP1, but not minor capsid proteins, interacts with the viral genome in vivo. In addition, semiquantitative PCR assays performed on cross-linked DNA complexes revealed that VP1 binds to all regions of the viral genome but significantly more to the regulatory region. The implications of such an interaction for viral infectivity are discussed.     

Unique pattern of convergent envelope evolution in simian immunodeficiency virus-infected rapid progressor macaques: association with CD4-independent usage of CCR5

The rate of disease development in simian immunodeficiency virus (SIV) infection of macaques varies considerably among individual macaques. While the majority of macaques inoculated with pathogenic SIV develop AIDS within a period of 1 to 2 years, a minority exhibit a rapid disease course characterized by absence or transience of humoral and cellular immune responses and high levels of virus replication with widespread dissemination of SIV in macrophages and multinucleated giant cells. The goal of this study was to examine viral evolution in three SIVsmE543-3-inoculated rapid progressors to determine the contribution of viral evolution to the development of rapid disease and the effect of the absence of immune pressure upon viral evolution. PCR was used to amplify and clone the entire SIV genome from tissues collected at necropsy, and the course of viral evolution was assessed by env sequences cloned from sequential plasma samples of one rapid progressor (RP) macaque. The majority of sequence changes in RP macaques occurred in the envelope gene. Substitutions were observed in all three animals at specific conserved residues in envelope, including loss of a glycosylation site in V1/V2, a D-to-N/V substitution in a highly conserved GDPE motif, and a P-to-V/H/T substitution in the V3 loop analog. A cell-cell fusion assay revealed that representative env clones utilized CCR5 as a coreceptor, independent of CD4. The selection of specific substitutions in envelope in RP macaques suggests novel selection pressures on virus in such animals and suggests that viral variants that evolve in these animals may play a role in disease progression.     

Error thresholds and the constraints to RNA virus evolution

RNA viruses are often thought of as possessing almost limitless adaptability as a result of their extreme mutation rates. However, high mutation rates also put a cap on the size of the viral genome by establishing an error threshold, beyond which lethal numbers of deleterious mutations accumulate. Herein, I argue that a lack of genomic space means that RNA viruses will be subject to important evolutionary constraints because specific sequences are required to encode multiple and often conflicting functions. Empirical evidence for these constraints, and how they limit viral adaptability, is now beginning to accumulate. Documenting the constraints to RNA virus evolution has important implications for predicting the emergence of new viruses and for improving therapeutic procedures.     

Pacing a small cage: mutation and RNA viruses

RNA viruses have an extremely high mutation rate, and we argue that the most plausible explanation for this is a trade-off with replication speed. We suggest that research into further increasing this mutation rate artificially as an antiviral treatment requires a theoretical reevaluation, especially relating to the so-called error threshold. The main evolutionary consequence of a high mutation rate appears to have been to restrict RNA viruses to a small genome; they thus rapidly exploit a limited array of possibilities. Investigating this constraint to their evolution, and how it is occasionally overcome, promises to be fruitful. We explain the many terms used in investigating RNA viral evolution and highlight the specific experimental and comparative work that needs to be done.     

HIV-1 evolution: frustrating therapies,but disclosing molecular mechanisms

Replication of HIV-1 under selective pressure frequently results in the evolution of virus variants that replicate more efficiently under the applied conditions. For example, in patients on antiretroviral therapy, such evolution can result in variants that are resistant to the HIV-1 inhibitors, thus frustrating the therapy. On the other hand, virus evolution can help us to understand the molecular mechanisms that underlie HIV-1 replication. For example, evolution of a defective virus mutant can result in variants that overcome the introduced defect by restoration of the original sequence or by the introduction of additional mutations in the viral genome. Analysis of the evolution pathway can reveal the requirements of the element under study and help to understand its function. Analysis of the escape routes may generate new insight in the viral life cycle and result in the identification of unexpected biological mechanisms. We have developed in vitro HIV-1 evolution into a systematic research tool that allows the study of different aspects of the viral replication cycle. We will briefly review this method of forced virus evolution and provide several examples that illustrate the power of this approach.     

Synethesis and integration of viral DNA in chicken cells at different time after infection with various multiplicities of avian oncornavirus.

To see if integration of the provirus resulting from RNA tumor virus infection is limited to specific sites in the cell DNA, the variation in the number of copies of virus-specific DNA produced and integrated in chicken embryo fibroblasts after RAV-2 infection with different multiplicities has been determined at short times, long times, and several transfers after infection. The number of copies of viral DNA in cells was determined by initial hybridization kinetics of single-stranded viral complementary DNA with a moderate excess of cell DNA. The approach took into account the different sizes of cell DNA and complementary DNA in the hybridization mixture. It was found that uninfected chicken embryo fibroblasts have approximately seven copies, part haploid genome of DNA sequences homologous to part of the Rous-association virus 2 (RAV-2) genome. Infection with RAV-2 adds additional copies, and different sequences, of RAV -2- specific DNA. By 13 h postinfection, there are 3 to 10 additional copies per haploid genome. This number can not be increased by increasing the multiplicity of infection, and stays relatively constant up to 20 h postinfection, when some of the additional viral DNA is integrated. Between 20 and 40 h postinfection, the cells accumulated up to 100 copies per haploid genome of viral DNA. Most of these are unintegrated. This number decreases with cell transfer, until cells are left with one to three copies of additional viral DNA sequences per haploid genome, of which most are integrated. The finding that viral infection causes the permanent addition of one to three copies of integrated viral DNA, despite the cells being confronted with up to 100 copies per haploid genome after infection, is consistent with a hypothesis that chicken cells contain a limited number of specific integration sites for the oncornavirus genome.     

The genomic rate of molecular adaptation of the human influenza A virus

Quantifying adaptive evolution at the genomic scale is an essential yet challenging aspect of evolutionary biology. Here, we develop a method that extends and generalizes previous approaches to estimate the rate of genomic adaptation in rapidly evolving populations and apply it to a large data set of complete human influenza A virus genome sequences. In accord with previous studies, we observe particularly high rates of adaptive evolution in domain 1 of the viral hemagglutinin (HA1). However, our novel approach also reveals previously unseen adaptation in other viral genes. Notably, we find that the rate of adaptation (per codon per year) is higher in surface residues of the viral neuraminidase than in HA1, indicating strong antibody-mediated selection on the former. We also observed high rates of adaptive evolution in several nonstructural proteins, which may relate to viral evasion of T-cell and innate immune responses. Furthermore, our analysis provides strong quantitative support for the hypothesis that human H1N1 influenza experiences weaker antigenic selection than H3N2. As well as shedding new light on the dynamics and determinants of positive Darwinian selection in influenza viruses, the approach introduced here is applicable to other pathogens for which densely sampled genome sequences are available, and hence is ideally suited to the interpretation of next-generation genome sequencing data.