Historical Perspective on the Discovery of the Quasispecies Concept

Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture

BACKGROUND AND PURPOSE

There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis.

EXPERIMENTAL APPROACH

The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured.

KEY RESULTS

Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species.

CONCLUSION AND IMPLICATIONS

SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

Genetic Variability Among U.S.-Sentinel Cotton Plot Cotton Leafroll Dwarf Virus and Globally Available Reference Isolates Based on ORF0 Diversity

The aphid-transmitted polerovirus, cotton leafroll dwarf virus (CLRDV), first characterized from symptomatic cotton plants in South America, has been identified in commercial cotton plantings in the United States. Here, the CLRDV intraspecific diversity was investigated by comparative sequence analysis of the most divergent CLRDV coding region, ORF0/P0. Bayesian analysis of ORF0 sequences for U.S. and reference populations resolved three well-supported sister clades comprising one U.S. and two South American lineages. Principal component analysis (PCA) identified seven statistically supported intraspecific populations. The Bayesian phylogeny and PCA dendrogram-inferred relationships were congruent. Population analysis of ORF0 sequences indicated most lineages have evolved under negative selection, albeit certain sites/isolates evolved under positive selection. Both U.S. and South American isolates exhibited extensive ORF0 diversity. At least two U.S. invasion foci were associated with their founder populations in Alabama-Georgia and eastern Texas. The Alabama-Georgia founder is implicated as the source of recent widespread expansion and establishment of secondary disease foci throughout the southeastern-central United States. Based on the geographically restricted distribution, spread of another extant Texas population appeared impeded by a population bottleneck. Extant CLRDV isolates represent several putative introductions potentially associated with catastrophic weather events dispersing viruliferous cotton aphids of unknown origin(s).

Toward a Categorization of Virus-ncRNA Interactions in the World of RNA to Disentangle the Tiny Secrets of Dengue Virus

In recent years, the function of noncoding RNAs (ncRNAs) as regulatory molecules of cell physiology has begun to be better understood. Advances in viral molecular biology have shown that host ncRNAs, cellular factors, and virus-derived ncRNAs and their interplay are strongly disturbed during viral infections. Nevertheless, the folding of RNA virus genomes has also been identified as a critical factor in regulating canonical and non-canonical functions. Due to the influence of host ncRNAs and the structure of RNA viral genomes, complex molecular and cellular processes in infections are modulated. We propose three main categories to organize the current information about RNA-RNA interactions in some well-known human viruses. The first category shows examples of host ncRNAs associated with the immune response triggered in viral infections. Even though miRNAs introduce a standpoint, they are briefly presented to keep researchers moving forward in uncovering other RNAs. The second category outlines interactions between virus-host ncRNAs, while the third describes how the structure of the RNA viral genome serves as a scaffold for processing virus-derived RNAs. Our grouping may provide a comprehensive framework to classify ncRNA-host-cell interactions for emerging viruses and diseases. In this sense, we introduced them to organize DENV-host-cell interactions.

RNA three-dimensional structure drives the sequence organization of potato spindle tuber viroid quasispecies

RNA viruses and viroids exist and evolve as quasispecies due to error-prone replication. Quasispecies consist of a few dominant master sequences alongside numerous variants that contribute to genetic diversity. Upon environmental changes, certain variants within quasispecies have the potential to become the dominant sequences, leading to the emergence of novel infectious strains. However, the emergence of new infectious variants remains unpredictable. Using mutant pools prepared by saturation mutagenesis of selected stem and loop regions, our study of potato spindle tuber viroid (PSTVd) demonstrates that mutants forming local three-dimensional (3D) structures similar to the wild type (WT) are more likely to accumulate in PSTVd quasispecies. The selection mechanisms underlying this biased accumulation are likely associated with cell-to-cell movement and long-distance trafficking. Moreover, certain trafficking-defective PSTVd mutants can be spread by functional sister genomes in the quasispecies. Our study reveals that the RNA 3D structure of stems and loops constrains the evolution of viroid quasispecies. Mutants with a structure similar to WT have a higher likelihood of being maintained within the quasispecies and can potentially give rise to novel infectious variants. These findings emphasize the potential of targeting RNA 3D structure as a more robust approach to defend against viroid infections.

Cell-cell communication and initial population composition shape the structure of potato spindle tuber viroid quasispecies

RNA viruses and viroids replicate with high mutation rates, forming quasispecies, population of variants centered around dominant sequences. The mechanisms governing quasispecies remain unclear. Plasmodesmata regulate viroid movement and were hypothesized to impact viroid quasispecies. Here, we sequenced the progeny of potato spindle tuber viroid intermediate (PSTVd-I) strain from mature guard cells lacking plasmodesmal connections and from in vitro-cultivated mesophyll cell protoplasts from systemic leaves of early-infected tomato (Solanum lycopersicum) plants. Remarkably, more variants accumulated in guard cells compared to whole leaves. Similarly, after extended cell culture, we observed more variants in cultivated mesophyll protoplasts. Coinfection and single-cell sequencing experiments demonstrated that the same plant cell can be infected multiple times by the same or different PSTVd sequences. To study the impact of initial population composition on PSTVd-I quasispecies, we conducted coinfections with PSTVd-I and variants. Two inoculum ratios (10:1 or 1:10) established quasispecies with or without PSTVd-I as the master sequence. In the absence of the master sequence, the percentage of novel variants initially increased. Moreover, a 1:1 PSTVd-I/variant RNA ratio resulted in PSTVd-I dominating (>50%), while the variants reached 20%. After PSTVd-I-only infection, the variants reached around 10%, while after variant-only infection, the variants were significantly more than 10%. These results emphasize the role of cell-to-cell communication and initial population composition in shaping PSTVd quasispecies.

SARS-CoV-2 mutant spectra as variant of concern nurseries: endless variation?

Introduction

SARS-CoV-2 isolates of a given clade may contain low frequency genomes that encode amino acids or deletions which are typical of a different clade.

Methods

Here we use high resolution ultra-deep sequencing to analyze SARS-CoV-2 mutant spectra.

Results

In 6 out of 11 SARS-CoV-2 isolates from COVID-19 patients, the mutant spectrum of the spike (S)-coding region included two or more amino acids or deletions, that correspond to discordant viral clades. A similar observation is reported for laboratory populations of SARS-CoV-2 USA-WA1/2020, following a cell culture infection in the presence of remdesivir, ribavirin or their combinations. Moreover, some of the clade-discordant genome residues are found in the same haplotype within an amplicon.

Discussion

We evaluate possible interpretations of these findings, and reviewed precedents for rapid selection of genomes with multiple mutations in RNA viruses. These considerations suggest that intra-host evolution may be sufficient to generate minority sequences which are closely related to sequences typical of other clades. The results provide a model for the origin of variants of concern during epidemic spread─in particular Omicron lineages─that does not require prolonged infection, involvement of immunocompromised individuals, or participation of intermediate, non-human hosts.

A Deep Sequencing Strategy for Investigation of Virus Variants within African Swine Fever Virus-Infected Pigs

African swine fever virus (ASFV) is the causative agent of African swine fever, an economically important disease of pigs, often with a high case fatality rate. ASFV has demonstrated low genetic diversity among isolates collected within Eurasia. To explore the influence of viral variants on clinical outcomes and infection dynamics in pigs experimentally infected with ASFV, we have designed a deep sequencing strategy. The variant analysis revealed unique SNPs at <10% frequency in several infected pigs as well as some SNPs that were found in more than one pig. In addition, a deletion of 10,487 bp (resulting in the complete loss of 21 genes) was present at a nearly 100% frequency in the ASFV DNA from one pig at position 6362-16849. This deletion was also found to be present at low levels in the virus inoculum and in two other infected pigs. The current methodology can be used for the currently circulating Eurasian ASFVs and also adapted to other ASFV strains and genotypes. Comprehensive deep sequencing is critical for following ASFV molecular evolution, especially for the identification of modifications that affect virus virulence.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

The Past Is Present: Coevolution of Viruses and Host Resistance Within Geographic Centers of Plant Diversity

Understanding the coevolutionary history of plants, pathogens, and disease resistance is vital for plant pathology. Here, I review Francis O. Holmes's work with tobacco mosaic virus (TMV) framed by the foundational work of Nikolai Vavilov on the geographic centers of origin of plants and crop wild relatives (CWRs) and T. Harper Goodspeed's taxonomy of the genus Nicotiana. Holmes developed a hypothesis that the origin of host resistance to viruses was due to coevolution of both at a geographic center. In the 1950s, Holmes proved that genetic resistance to TMV, especially dominant R-genes, was centered in South America for Nicotiana and other solanaceous plants, including Capsicum, potato, and tomato. One seeming exception was eggplant (Solanum melongena). Not until the acceptance of plate tectonics in the 1960s and recent advances in evolutionary taxonomy did it become evident that northeast Africa was the home of eggplant CWRs, far from Holmes's geographic center for TMV-R-gene coevolution. Unbeknownst to most plant pathologists, Holmes's ideas predated those of H.H. Flor, including experimental proof of the gene-for-gene interaction, identification of R-genes, and deployment of dominant host genes to protect crop plants from virus-associated yield losses.

Pangenome graphs in infectious disease: a comprehensive genetic variation analysis of Neisseria meningitidis leveraging Oxford Nanopore long reads

Whole genome sequencing has revolutionized infectious disease surveillance for tracking and monitoring the spread and evolution of pathogens. However, using a linear reference genome for genomic analyses may introduce biases, especially when studies are conducted on highly variable bacterial genomes of the same species. Pangenome graphs provide an efficient model for representing and analyzing multiple genomes and their variants as a graph structure that includes all types of variations. In this study, we present a practical bioinformatics pipeline that employs the PanGenome Graph Builder and the Variation Graph toolkit to build pangenomes from assembled genomes, align whole genome sequencing data and call variants against a graph reference. The pangenome graph enables the identification of structural variants, rearrangements, and small variants (e.g., single nucleotide polymorphisms and insertions/deletions) simultaneously. We demonstrate that using a pangenome graph, instead of a single linear reference genome, improves mapping rates and variant calling for both simulated and real datasets of the pathogen Neisseria meningitidis. Overall, pangenome graphs offer a promising approach for comparative genomics and comprehensive genetic variation analysis in infectious disease. Moreover, this innovative pipeline, leveraging pangenome graphs, can bridge variant analysis, genome assembly, population genetics, and evolutionary biology, expanding the reach of genomic understanding and applications.

Host Cell Targets for Unconventional Antivirals against RNA Viruses

The recent COVID-19 crisis has highlighted the importance of RNA-based viruses. The most prominent members of this group are SARS-CoV-2 (coronavirus), HIV (human immunodeficiency virus), EBOV (Ebola virus), DENV (dengue virus), HCV (hepatitis C virus), ZIKV (Zika virus), CHIKV (chikungunya virus), and influenza A virus. With the exception of retroviruses which produce reverse transcriptase, the majority of RNA viruses encode RNA-dependent RNA polymerases which do not include molecular proofreading tools, underlying the high mutation capacity of these viruses as they multiply in the host cells. Together with their ability to manipulate the immune system of the host in different ways, their high mutation frequency poses a challenge to develop effective and durable vaccination and/or treatments. Consequently, the use of antiviral targeting agents, while an important part of the therapeutic strategy against infection, may lead to the selection of drug-resistant variants. The crucial role of the host cell replicative and processing machinery is essential for the replicative cycle of the viruses and has driven attention to the potential use of drugs directed to the host machinery as therapeutic alternatives to treat viral infections. In this review, we discuss small molecules with antiviral effects that target cellular factors in different steps of the infectious cycle of many RNA viruses. We emphasize the repurposing of FDA-approved drugs with broad-spectrum antiviral activity. Finally, we postulate that the ferruginol analog (18-(phthalimide-2-yl) ferruginol) is a potential host-targeted antiviral.

Quasispecies Nature of RNA Viruses: Lessons from the Past

Viral quasispecies are distinct but closely related mutants formed by the disparity in viral genomes due to recombination, mutations, competition, and selection pressure. Theoretical derivation for the origin of a quasispecies is owed to the error-prone replication by polymerase and mutants of RNA replicators. Here, we briefly addressed the theoretical and mathematical origin of quasispecies and their dynamics. The impact of quasispecies for major salient human pathogens is reviewed. In the current global scenario, rapid changes in geographical landscapes favor the origin and selection of mutants. It comes as no surprise that a cauldron of mutants poses a significant risk to public health, capable of causing pandemics. Mutation rates in RNA viruses are magnitudes higher than in DNA organisms, explaining their enhanced virulence and evolvability. RNA viruses cause the most devastating pandemics; for example, members of the Orthomyxoviridae family caused the great influenza pandemic (1918 flu or Spanish flu), the SARS (severe acute respiratory syndrome) and MERS (Middle East respiratory syndrome) outbreak, and the human immunodeficiency viruses (HIV), lentiviruses of the Retroviridae family, caused worldwide devastation. Rapidly evolving RNA virus populations are a daunting challenge for the designing of effective control measures like vaccines. Developing awareness of the evolutionary dispositions of RNA viral mutant spectra and what influences their adaptation and virulence will help curtail outbreaks of past and future pathogens.

The Global Impact of COVID-19: Historical Development, Molecular Characterization, Drug Discovery and Future Directions

In December 2019, an outbreak of a respiratory disease called the coronavirus disease 2019 (COVID-19) caused by a new coronavirus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) began in Wuhan, China. The SARS-CoV-2, an encapsulated positive-stranded RNA virus, spread worldwide with disastrous consequences for people's health, economies, and quality of life. The disease has had far-reaching impacts on society, including economic disruption, school closures, and increased stress and anxiety. It has also highlighted disparities in healthcare access and outcomes, with marginalized communities disproportionately affected by the SARS-CoV-2. The symptoms of COVID-19 range from mild to severe. There is presently no effective cure. Nevertheless, significant progress has been made in developing COVID-19 vaccine for different therapeutic targets. For instance, scientists developed multifold vaccine candidates shortly after the COVID-19 outbreak after Pfizer and AstraZeneca discovered the initial COVID-19 vaccines. These vaccines reduce disease spread, severity, and mortality. The addition of rapid diagnostics to microscopy for COVID-19 diagnosis has proven crucial. Our review provides a thorough overview of the historical development of COVID-19 and molecular and biochemical characterization of the SARS-CoV-2. We highlight the potential contributions from insect and plant sources as anti-SARS-CoV-2 and present directions for future research.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

Influence of viral genome properties on polymerase fidelity

The first step of viral evolution takes place during genome replication via the error-prone viral polymerase. Among the mutants that arise through this process, only a few well-adapted variants will be selected by natural selection, renewing the viral genome population. Viral polymerase-mediated errors are thought to occur stochastically. However, accumulating evidence suggests that viral polymerase-mediated mutations are heterogeneously distributed throughout the viral genome. Here, we review work that supports this concept and provides mechanistic insights into how specific features of the viral genome could modulate viral polymerase-mediated errors. A predisposition to accumulate viral polymerase-mediated errors at specific loci in the viral genome may guide evolution to specific pathways, thus opening new directions of research to better understand viral evolutionary dynamics.

Molecular characterization and evolution dynamics of influenza B viruses circulating in Germany from season 1996/1997 to 2019/2020

Influenza B viruses are responsible for significant disease burden caused by viruses of both the Yamagata- and Victoria-lineage. Since the circulating patterns of influenza B viruses in different countries vary we investigated molecular properties and evolution dynamics of influenza B viruses circulating in Germany between 1996 and 2020. A change of the dominant lineage occurred in Germany in seven seasons in over past 25 years. A total of 676 sequences of hemagglutinin coding domain 1 (HA1) and 516 sequences of neuraminidase (NA) genes of Yamagata- and Victoria-lineage viruses were analyzed using time-scaled phylogenetic tree. Phylogenetic analysis demonstrated that Yamagata-lineage viruses are more diverse than the Victoria-lineage viruses and could be divided into nine genetic groups whereas Victoria-lineage viruses presented six genetic groups. Comparative phylogenetic analyses of both the HA and NA segments together revealed a number of inter-lineage as well as inter- and intra-clade reassortants. We identified key amino acid substitutions in major HA epitopes such as in four antigenic sites and receptor-binding sites (RBS) and in the regions close to them, with most substitutions in the 120-loop of both lineage viruses. Altogether, seventeen substitutions were fixed over time within the Yamagata-lineage with twelve of them in the antigenic sites. Thirteen substitutions were identified within the Victoria-lineage, with eleven of them in the antigenic sites. Moreover, all Victoria-lineage viruses of the 2017/2018 season were characterized by a deletion of two amino acids at the position 162-163 in the antigenic site of HA1. The viruses with triple deletion Δ162-164 were found in Germany since season 2018/2019. We highlighted the interplay between substitutions in the glycosylation sites and RBS and antigenic epitope during HA evolution. The results obtained underscore the need for continuous monitoring of circulating influenza B viruses. Early detection of strains with genetic and antigenic variation is essential to predict the circulation patterns for the following season. Such information is important for the development of optimal vaccines and strategies for prevention and control of influenza.

Determinants of Virus Variation, Evolution, and Host Adaptation

Virus evolution is the change in the genetic structure of a viral population over time and results in the emergence of new viral variants, strains, and species with novel biological properties, including adaptation to new hosts. There are host, vector, environmental, and viral factors that contribute to virus evolution. To achieve or fine tune compatibility and successfully establish infection, viruses adapt to a particular host species or to a group of species. However, some viruses are better able to adapt to diverse hosts, vectors, and environments. Viruses generate genetic diversity through mutation, reassortment, and recombination. Plant viruses are exposed to genetic drift and selection pressures by host and vector factors, and random variants or those with a competitive advantage are fixed in the population and mediate the emergence of new viral strains or species with novel biological properties. This process creates a footprint in the virus genome evident as the preferential accumulation of substitutions, insertions, or deletions in areas of the genome that function as determinants of host adaptation. Here, with respect to plant viruses, we review the current understanding of the sources of variation, the effect of selection, and its role in virus evolution and host adaptation.

Functional cure of a chronic virus infection by shifting the virus - host equilibrium state

The clinical handling of chronic virus infections remains a challenge. Here we describe recent progress in the understanding of virus - host interaction dynamics. Based on the systems biology concept of multi-stability and the prediction of multiplicative cooperativity between virus-specific cytotoxic T cells and neutralising antibodies, we argue for the requirements to engage multiple immune system components for functional cure strategies. Our arguments are derived from LCMV model system studies and are translated to HIV-1 infection.

Impact of Intrahost NS5 Nucleotide Variations on Dengue Virus Replication

Due to the nature of RNA viruses, their high mutation rates produce a population of closely related but genetically diverse viruses, termed quasispecies. To determine the role of quasispecies in DENV disease severity, 22 isolates (10 from mild cases, 12 from fatal cases) were obtained, amplified, and sequenced with Next Generation Sequencing using the Illumina MiSeq platform. Using variation calling, unique wildtype nucleotide positions were selected and analyzed for variant nucleotides between mild and fatal cases. The analysis of variant nucleotides between mild and fatal cases showed 6 positions with a significant difference of p < 0.05 with 1 position in the structural region, and 5 positions in the non-structural (NS) regions. All variations were found to have a higher percentage in fatal cases. To further investigate the genetic changes that affect the virus’s properties, reverse genetics (rg) viruses containing substitutions with the variations were generated and viral growth properties were examined. We found that the virus variant rgNS5-T7812G (G81G) had higher replication rates in both Baby hamster kidney cells (BHK-21) and Vero cells while rgNS5-C9420A (A617A) had a higher replication rate only in BHK-21 cells compared to wildtype virus. Both variants were considered temperature sensitive whereby the viral titers of the variants were relatively lower at 39°C, but was higher at 35 and 37°C. Additionally, the variants were thermally stable compared to wildtype at temperatures of 29, 37, and 39°C. In conclusion, viral quasispecies found in isolates from the 2015 DENV epidemic, resulted in variations with significant difference between mild and fatal cases. These variations, NS5-T7812G (G81G) and NS5-C9420A (A617A), affect viral properties which may play a role in the virulence of DENV.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

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

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

The International Virus Bioinformatics Meeting 2022

The International Virus Bioinformatics Meeting 2022 took place online, on 23-25 March 2022, and has attracted about 380 participants from all over the world. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The participants created a highly interactive scientific environment even without physical face-to-face interactions. This meeting is a focal point to gain an insight into the state-of-the-art of the virus bioinformatics research landscape and to interact with researchers in the forefront as well as aspiring young scientists. The meeting featured eight invited and 18 contributed talks in eight sessions on three days, as well as 52 posters, which were presented during three virtual poster sessions. The main topics were: SARS-CoV-2, viral emergence and surveillance, virus-host interactions, viral sequence analysis, virus identification and annotation, phages, and viral diversity. This report summarizes the main research findings and highlights presented at the meeting.

Plasticity in structure and assembly of SARS-CoV-2 nucleocapsid protein

Worldwide SARS-CoV-2 sequencing efforts track emerging mutations in its spike protein, as well as characteristic mutations in other viral proteins. Besides their epidemiological importance, the observed SARS-CoV-2 sequences present an ensemble of viable protein variants, and thereby a source of information on viral protein structure and function. Charting the mutational landscape of the nucleocapsid (N) protein that facilitates viral assembly, we observe variability exceeding that of the spike protein, with more than 86% of residues that can be substituted, on average by three to four different amino acids. However, mutations exhibit an uneven distribution that tracks known structural features but also reveals highly protected stretches of unknown function. One of these conserved regions is in the central disordered linker proximal to the N-G215C mutation that has become dominant in the Delta variant, outcompeting G215 variants without further spike or N-protein substitutions. Structural models suggest that the G215C mutation stabilizes conserved transient helices in the disordered linker serving as protein-protein interaction interfaces. Comparing Delta variant N-protein to its ancestral version in biophysical experiments, we find a significantly more compact and less disordered structure. N-G215C exhibits substantially stronger self-association, shifting the unliganded protein from a dimeric to a tetrameric oligomeric state, which leads to enhanced coassembly with nucleic acids. This suggests that the sequence variability of N-protein is mirrored by high plasticity of N-protein biophysical properties, which we hypothesize can be exploited by SARS-CoV-2 to achieve greater efficiency of viral assembly, and thereby enhanced infectivity.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Lethal Mutagenesis of RNA Viruses and Approved Drugs with Antiviral Mutagenic Activity

In RNA viruses, a small increase in their mutation rates can be sufficient to exceed their threshold of viability. Lethal mutagenesis is a therapeutic strategy based on the use of mutagens, driving viral populations to extinction. Extinction catastrophe can be experimentally induced by promutagenic nucleosides in cell culture models. The loss of HIV infectivity has been observed after passage in 5-hydroxydeoxycytidine or 5,6-dihydro-5-aza-2′-deoxycytidine while producing a two-fold increase in the viral mutation frequency. Among approved nucleoside analogs, experiments with polioviruses and other RNA viruses suggested that ribavirin can be mutagenic, although its mechanism of action is not clear. Favipiravir and molnupiravir exert an antiviral effect through lethal mutagenesis. Both drugs are broad-spectrum antiviral agents active against RNA viruses. Favipiravir incorporates into viral RNA, affecting the G→A and C→U transition rates. Molnupiravir (a prodrug of β-d-N⁴-hydroxycytidine) has been recently approved for the treatment of SARS-CoV-2 infection. Its triphosphate derivative can be incorporated into viral RNA and extended by the coronavirus RNA polymerase. Incorrect base pairing and inefficient extension by the polymerase promote mutagenesis by increasing the G→A and C→U transition frequencies. Despite having remarkable antiviral action and resilience to drug resistance, carcinogenic risks and genotoxicity are important concerns limiting their extended use in antiviral therapy.

Updated HIV-1 Consensus Sequences Change but Stay Within Similar Distance From Worldwide Samples

HIV consensus sequences are used in various bioinformatic, evolutionary, and vaccine related research. Since the previous HIV-1 subtype and CRF consensus sequences were constructed in 2002, the number of publicly available HIV-1 sequences have grown exponentially, especially from non-EU and US countries. Here, we reconstruct 90 new HIV-1 subtype and CRF consensus sequences from 3,470 high-quality, representative, full genome sequences in the LANL HIV database. While subtypes and CRFs are unevenly spread across the world, in total 89 countries were represented. For consensus sequences that were based on at least 20 genomes, we found that on average 2.3% (range 0.8–10%) of the consensus genome site states changed from 2002 to 2021, of which about half were nucleotide state differences and the rest insertions and deletions. Interestingly, the 2021 consensus sequences were shorter than in 2002, and compared to 4,674 HIV-1 worldwide genome sequences, the 2021 consensuses were somewhat closer to the worldwide genome sequences, i.e., showing on average fewer nucleotide state differences. Some subtypes/CRFs have had limited geographical spread, and thus sampling of subtypes/CRFs is uneven, at least in part, due to the epidemiological dynamics. Thus, taken as a whole, the 2021 consensus sequences likely are good representations of the typical subtype/CRF genome nucleotide states. The new consensus sequences are available at the LANL HIV database.

Plasticity in structure and assembly of SARS-CoV-2 nucleocapsid protein

Worldwide SARS-CoV-2 sequencing efforts track emerging mutations in its spike protein, as well as characteristic mutations in other viral proteins. Besides their epidemiological importance, the observed SARS-CoV-2 sequences present an ensemble of viable protein variants, and thereby a source of information on viral protein structure and function. Charting the mutational landscape of the nucleocapsid (N) protein that facilitates viral assembly, we observe variability exceeding that of the spike protein, with more than 86% of residues that can be substituted, on average by 3-4 different amino acids. However, mutations exhibit an uneven distribution that tracks known structural features but also reveals highly protected stretches of unknown function. One of these conserved regions is in the central disordered linker proximal to the N-G215C mutation that has become dominant in the Delta variant, outcompeting G215 variants without further spike or N-protein substitutions. Structural models suggest that the G215C mutation stabilizes conserved transient helices in the disordered linker serving as protein-protein interaction interfaces. Comparing Delta variant N-protein to its ancestral version in biophysical experiments, we find a significantly more compact and less disordered structure. N-G215C exhibits substantially stronger self-association, shifting the unliganded protein from a dimeric to a tetrameric oligomeric state, which leads to enhanced co-assembly with nucleic acids. This suggests that the sequence variability of N-protein is mirrored by high plasticity of N-protein biophysical properties, which we hypothesize can be exploited by SARS-CoV-2 to achieve greater efficiency of viral assembly, and thereby enhanced infectivity.

RNA-Binding Proteins as Regulators of Internal Initiation of Viral mRNA Translation

Viruses are obligate intracellular parasites that depend on the host’s protein synthesis machinery for translating their mRNAs. The viral mRNA (vRNA) competes with the host mRNA to recruit the translational machinery, including ribosomes, tRNAs, and the limited eukaryotic translation initiation factor (eIFs) pool. Many viruses utilize non-canonical strategies such as targeting host eIFs and RNA elements known as internal ribosome entry sites (IRESs) to reprogram cellular gene expression, ensuring preferential translation of vRNAs. In this review, we discuss vRNA IRES-mediated translation initiation, highlighting the role of RNA-binding proteins (RBPs), other than the canonical translation initiation factors, in regulating their activity.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.