Basic concepts in RNA virus evolution

Binding Evolution of the Dengue Virus Envelope Against DC-SIGN: A Combined Approach of Phylogenetics and Molecular Dynamics Analyses Over 30 Years of Dengue Virus in Brazil

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's "Through the Looking-Glass", postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.

Genetic Diversity of Human Respiratory Syncytial Virus during COVID-19 Pandemic in Yaoundé, Cameroon, 2020-2021

Worldwide, human respiratory syncytial virus (HRSV) is a major cause of severe infections of the lower respiratory system, affecting individuals of all ages. This study investigated the genetic variability of HRSV during the COVID-19 outbreak in Yaoundé; nasopharyngeal samples positive for HRSV were collected from different age groups between July 2020 and October 2021. A semi-nested RT-PCR was performed on the second hypervariable region of the G gene of detected HRSV, followed by sequencing and phylogenetic assessment. Throughout the study, 40 (37.7%) of the 106 HRSV-positive samples successfully underwent G-gene amplification. HRSV A and HRSV B co-circulated at rates of 47.5% and 52.5%, respectively. HRSV A clustered in the GA2.3.5 genetic lineage (ON1) and HRSV B clustered in the GB5.0.5a genetic lineage (BA9). Differences in circulating genotypes were observed between pre- and post-pandemic years for HRSV A. Predictions revealed potential N-glycosylation sites at positions 237-318 of HRSV A and positions 228-232-294 of HRSV B. This study reports the molecular epidemiology of HRSV in Cameroon during the COVID-19 pandemic. It describes the exclusive co-circulation of two genetic lineages. These findings highlight the importance of implementing comprehensive molecular surveillance to prevent the unexpected emergence of other diseases.

Current Trends in RNA Virus Detection via Nucleic Acid Isothermal Amplification-Based Platforms

Ribonucleic acid (RNA) viruses are one of the major classes of pathogens that cause human diseases. The conventional method to detect RNA viruses is real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), but it has some limitations. It is expensive and time-consuming, with infrastructure and trained personnel requirements. Its high throughput requires sophisticated automation and large-scale infrastructure. Isothermal amplification methods have been explored as an alternative to address these challenges. These methods are rapid, user-friendly, low-cost, can be performed in less specialized settings, and are highly accurate for detecting RNA viruses. Microfluidic technology provides an ideal platform for performing virus diagnostic tests, including sample preparation, immunoassays, and nucleic acid-based assays. Among these techniques, nucleic acid isothermal amplification methods have been widely integrated with microfluidic platforms for RNA virus detection owing to their simplicity, sensitivity, selectivity, and short analysis time. This review summarizes some common isothermal amplification methods for RNA viruses. It also describes commercialized devices and kits that use isothermal amplification techniques for SARS-CoV-2 detection. Furthermore, the most recent applications of isothermal amplification-based microfluidic platforms for RNA virus detection are discussed in this article.

Bioinformatic Approaches for Comparative Analysis of Viruses

The field of viral genomic studies has experienced an unprecedented increase in data volume. New strains of known viruses are constantly being added to the GenBank database and so are completely new species with little or no resemblance to our databases of sequences. In addition to this, metagenomic techniques have the potential to further increase the number and rate of sequenced genomes. Besides, it is important to consider that viruses have a set of unique features that often break down molecular biology dogmas, e.g., the flux of information from RNA to DNA in retroviruses and the use of RNA molecules as genomes. As a result, extracting meaningful information from viral genomes remains a challenge and standard methods for comparing the unknown and our databases of characterized sequences may need adaptations. Thus, several bioinformatic approaches and tools have been created to address the challenge of analyzing viral data. This chapter offers descriptions and protocols of some of the most important bioinformatic techniques for comparative analysis of viruses. The authors also provide comments and discussion on how viruses' unique features can affect standard analyses and how to overcome some of the major sources of problems. Protocols and topics emphasize online tools (which are more accessible to users) and give the real experience of what most bioinformaticians do in day-by-day work with command-line pipelines. The topics discussed include (1) clustering related genomes, (2) whole genome multiple sequence alignments for small RNA viruses, (3) protein alignment for marker genes and species affiliation, (4) variant calling and annotation, and (5) virome analyses and pathogen identification.

Further Molecular Diagnosis Determines Lack of Evidence for Real Seed Transmission of Tomato Leaf Curl New Delhi Virus in Cucurbits

Begomoviruses (family Geminiviridae) cause serious diseases in many crop families. Since 2013, the Spanish isolate of tomato leaf curl New Delhi virus (ToLCNDV) has been a limiting factor for cucurbits production in the Mediterranean basin, forcing farmers to adapt new management and control techniques. Although it is well-known that begomoviruses are naturally transmitted by the whitefly Bemisia tabaci, the capacity of these viruses to be vertically transmitted through seeds remains controversial. Clarifying the potential ToLCNDV seed transmission is essential to understand the epidemiology of this threating-for-cucurbits virus and to design appropriate control strategies. We assessed ToLCNDV distribution in the leaves, flowers and seeds of the infected plants of susceptible Cucumis melo accessions and toleration to the infected genotypes of Cucurbita moschata by conventional and quantitative PCR. We analyzed whether the viral particle was transmitted to offspring. We also evaluated ToLCNDV presence in commercial seeds of cucurbits (zucchini (Cucurbita pepo), melon (C. melo), cucumber (Cucumis sativus) and watermelon (Citrullus lanatus)) and in their progenies. As the assayed seedlings remained symptomless, we increased the reliability and accuracy of detection in these samples by searching for replicative forms of ToLCNDV by combining Southern blot hybridization and rolling-circle amplification (RCA). However, integral genomic DNA was not identified in the plants of offspring. Although the seedborne nature of ToLCNDV was confirmed, our results do not support the transmission of this virus from contaminated seeds to progeny.

HMMPolish: a coding region polishing tool for TGS-sequenced RNA viruses

Access to accurate viral genomes is important to downstream data analysis. Third-generation sequencing (TGS) has recently become a popular platform for virus sequencing because of its long read length. However, its per-base error rate, which is higher than next-generation sequencing, can lead to genomes with errors. Polishing tools are thus needed to correct errors either before or after sequence assembly. Despite promising results of available polishing tools, there is still room to improve the error correction performance to perform more accurate genome assembly. The errors, particularly those in coding regions, can hamper analysis such as linage identification and variant monitoring. In this work, we developed a novel pipeline, HMMPolish, for correcting (polishing) errors in protein-coding regions of known RNA viruses. This tool can be applied to either raw TGS reads or the assembled sequences of the target virus. By utilizing profile Hidden Markov Models of protein families/domains in known viruses, HMMPolish can correct errors that are ignored by available polishers. We extensively validated HMMPolish on 34 datasets that covered four clinically important viruses, including HIV-1, influenza-A, norovirus, and severe acute respiratory syndrome coronavirus 2. These datasets contain reads with different properties, such as sequencing depth and platforms (PacBio or Nanopore). The benchmark results against popular/representative polishers show that HMMPolish competes favorably on error correction in coding regions of known RNA viruses.

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation

Intrinsically disordered regions (IDRs) are critical for cellular function, yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. We focussed here on the first two domains of coronavirus nucleocapsid proteins, the disordered N-terminal domain (NTD) followed by the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system to explore the interplay between an IDR adjacent to a folded domain, and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved fuzzy RNA:protein complex is found across nucleocapsid protein orthologs, despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

Molecular Population Genetics of Aspen Mosaic-Associated Virus in Finland and Sweden

Aspen mosaic-associated virus (AsMaV) is a newly identified Emaravirus, in the family Fimoviridae, Bunyavirales, associated with mosaic symptoms in aspen trees (Populus tremula). Aspen trees are widely distributed in Europe and understanding the population structure of AsMaV may aid in the development of better management strategies. The virus genome consists of five negative-sense single-stranded RNA (-ssRNA) molecules. To investigate the genetic diversity and population parameters of AsMaV, different regions of the genome were amplified and analyzed and full-length sequence of the divergent isolates were cloned and sequenced. The results show that RNA3 or nucleoprotein is a good representative for studying genetic diversity in AsMaV. Developed RT-PCR-RFLP was able to identify areas with a higher number of haplotypes and could be applied for screening the large number of samples. In general, AsMaV has a conserved genome and based on the phylogenetic studies, geographical structuring was observed in AsMaV isolates from Sweden and Finland, which could be attributed to founder effects. The genome of AsMaV is under purifying selection but not distributed uniformly on genomic RNAs. Distant AsMaV isolates displayed amino acid sequence variations compared to other isolates, and bioinformatic analysis predicted potential post-translational modification sites in some viral proteins.

Pest categorisation of the avocado sunblotch viroid

The EFSA Panel on Plant Health conducted a pest categorisation of the avocado sunblotch viroid (ASBVd) for the EU. The identity of ASBVd, a member of the genus Avsunviroid (family Avsunviroidae) is clearly defined and detection and identification methods are available. The pathogen is not included in the EU Commission Implementing Regulation 2019/2072. ASBVd has been reported in Australia, Ghana, Guatemala, Israel, Mexico, Peru, South Africa, USA (California, Florida) and Venezuela. In the EU, it has been reported in Greece (Crete Island) and Spain. The pathogen could establish in the EU wherever avocado (Persea americana) is grown. The only known natural host of ASBVd is avocado to which it causes the severe 'avocado sunblotch' disease, characterised by white, yellow, red or necrotic depressed areas or scars on the fruit surface, bleached veins and petioles of the leaf, and rectangular cracking patterns in the bark of the old branches. Fruit yield and quality are severely diminished. ASBVd infects under experimental conditions a few more species in the family Lauraceae. The viroid is naturally transmitted at an extremely high rate by seeds (up to 100% in asymptomatically infected trees), but with a low efficiency by pollen (only to the produced seeds), and possibly through root grafts. Plants for planting, including seeds, and fresh avocado fruits were identified as the most relevant pathways for further entry of ASBVd into the EU. Avocado crops are cultivated in southern EU countries. Should the pest further enter and establish in the EU, impact on the production of avocado is expected. Phytosanitary measures are available to prevent entry and spread of the viroid in the EU. ASBVd fulfils the criteria that are within the remit of EFSA to assess for it to be regarded as a potential Union quarantine pest.

Viruses of a key coral symbiont exhibit temperature-driven productivity across a reefscape

Viruses can affect coral health by infecting their symbiotic dinoflagellate partners (Symbiodiniaceae). Yet, viral dynamics in coral colonies exposed to environmental stress have not been studied at the reef scale, particularly within individual viral lineages. We sequenced the viral major capsid protein (mcp) gene of positive-sense single-stranded RNA viruses known to infect symbiotic dinoflagellates ('dinoRNAVs') to analyze their dynamics in the reef-building coral, Porites lobata. We repeatedly sampled 54 colonies harboring Cladocopium C15 dinoflagellates, across three environmentally distinct reef zones (fringing reef, back reef, and forereef) around the island of Moorea, French Polynesia over a 3-year period and spanning a reef-wide thermal stress event. By the end of the sampling period, 28% (5/18) of corals in the fringing reef experienced partial mortality versus 78% (14/18) of corals in the forereef. Over 90% (50/54) of colonies had detectable dinoRNAV infections. Reef zone influenced the composition and richness of viral mcp amino acid types ('aminotypes'), with the fringing reef containing the highest aminotype richness. The reef-wide thermal stress event significantly increased aminotype dispersion, and this pattern was strongest in the colonies that experienced partial mortality. These findings demonstrate that dinoRNAV infections respond to environmental fluctuations experienced in situ on reefs. Further, viral productivity will likely increase as ocean temperatures continue to rise, potentially impacting the foundational symbiosis underpinning coral reef ecosystems.

An alternating transmission model between mice and mosquitoes for genetic study of dengue virus

Rapidly increased incidence and prevalence of dengue virus serotype 2 (DENV-2) in recent decades highlight the need for better understanding of the selective pressures that drive genetic and phenotypic changes in this virus. We simulated the transfer of DENV-2 between human hosts and mosquito vectors by horizontally transmitting the virus between suckling mice and Aedes aegypti (Linnaeus, Diptera: Culicidae). A total of 3 cycles of alternating transmission were performed and 3 passages of virus population were harvested from the infected sucking mice. The viral titer in mice brain and infectivity to mosquitoes of theses viral populations were tested. The genome of the viruses was also sequenced. Results showed that viral titer were similar and infection rate in the mosquitoes were not significantly different among those 3 passages. This in vivo model could be utilized to explore virus evolution and genetic variance in alternating transmission.

Host-Specific Interplay between Foot-and-Mouth Disease Virus 3D Polymerase and the Type-I Interferon Pathway

Foot-and-mouth disease (FMD) is a highly contagious viral disease affecting cloven-hoofed animals. One of the issues related to this disease is the persistence of its causative agent, foot-and-mouth disease virus (FMDV). While the mechanisms of FMDV persistence remain unclear, there are clues that it may be related to protein-protein interactions (PPI) between viral proteins and cellular proteins involved in the interferon (IFN) response. Since FMDV persistence has been described in cattle, sheep and goats but not in swine, we screened PPI involving FMDV proteins and sixteen major type-I IFN pathway proteins from these four species by nanoluciferase-2-hybrid complementation assay, in order to identify new PPI and determine their host specificity. As the results concerning the 3Dpol were the most interesting in view of the limited data concerning its role in immune escape, we decided to focus particularly on this protein. The identified PPI were confirmed by GST pull-down. We identified PPI between 3Dpol and seven IFN pathway proteins, namely, IKKα, IKKε, IRF3, IRF7, NEMO, MDA5 and MAVS. These PPI are conserved among the four studied species, with the exception of the one between 3Dpol and MAVS, which was only found with the swine protein. We also showed, using luciferase reporter assays, that 3Dpol could inhibit the induction phase of the IFN pathway. These results demonstrate, for the first time, a putative role for 3Dpol in FMDV innate immune escape.

Quasispecies Analysis of SARS-CoV-2 of 15 Different Lineages during the First Year of the Pandemic Prompts Scratching under the Surface of Consensus Genome Sequences

The tremendous majority of SARS-CoV-2 genomic data so far neglected intra-host genetic diversity. Here, we studied SARS-CoV-2 quasispecies based on data generated by next-generation sequencing (NGS) of complete genomes. SARS-CoV-2 raw NGS data had been generated for nasopharyngeal samples collected between March 2020 and February 2021 by the Illumina technology on a MiSeq instrument, without prior PCR amplification. To analyze viral quasispecies, we designed and implemented an in-house Excel file (“QuasiS”) that can characterize intra-sample nucleotide diversity along the genomes using data of the mapping of NGS reads. We compared intra-sample genetic diversity and global genetic diversity available from Nextstrain. Hierarchical clustering of all samples based on the intra-sample genetic diversity was performed and visualized with the Morpheus web application. NGS mapping data from 110 SARS-CoV-2-positive respiratory samples characterized by a mean depth of 169 NGS reads/nucleotide position and for which consensus genomes that had been obtained were classified into 15 viral lineages were analyzed. Mean intra-sample nucleotide diversity was 0.21 ± 0.65%, and 5357 positions (17.9%) exhibited significant (>4%) diversity, in ≥2 genomes for 1730 (5.8%) of them. ORF10, spike, and N genes had the highest number of positions exhibiting diversity (0.56%, 0.34%, and 0.24%, respectively). Nine hot spots of intra-sample diversity were identified in the SARS-CoV-2 NSP6, NSP12, ORF8, and N genes. Hierarchical clustering delineated a set of six genomes of different lineages characterized by 920 positions exhibiting intra-sample diversity. In addition, 118 nucleotide positions (0.4%) exhibited diversity at both intra- and inter-patient levels. Overall, the present study illustrates that the SARS-CoV-2 consensus genome sequences are only an incomplete and imperfect representation of the entire viral population infecting a patient, and that quasispecies analysis may allow deciphering more accurately the viral evolutionary pathways.

COVID-19 vs. Cancer Immunosurveillance: A Game of Thrones within an Inflamed Microenviroment

The COVID-19 pandemic accounts for more than 500 million confirmed infections and over 6 million deaths worldwide in the last 2 years. SARS-CoV-2 causes a highly complex form of inflammation that affects the human organism both acutely and chronically. In the same line, cancer as an inflammation-induced and immune-editing disease appears to cross-react with immune system at different levels including early interactions during carcinogenesis and later cross-talks within the tumor microenvironment. With all that in mind, a reasonable question one might address is whether the SARS-CoV-2 infection and the derived "long lasting inflammatory status" that is frequently observed in patients, might affect the cancer immunosurveillance mechanisms and consequently their risk of developing cancer, as well as the tumor and immune cell behaviors within the inflamed microenvironment. On this context, this review intends to outline and discuss the existing knowledge on SARS-CoV-2-mediated immunomodulation under the prism of changes that might be able to interfere with cancer cell immunoescape and the overall tumor progression and response to conventional therapeutics. Our goal is to highlight a potential interplay between the COVID-19 immunopathology and cancer immune-microenvironment that may pave the way for thorough investigation in the future.

COVIDOA: a novel evolutionary optimization algorithm based on coronavirus disease replication lifecycle

This paper presents a novel bio-inspired optimization algorithm called Coronavirus Optimization Algorithm (COVIDOA). COVIDOA is an evolutionary search strategy that mimics the mechanism of coronavirus when hijacking human cells. COVIDOA is inspired by the frameshifting technique used by the coronavirus for replication. The proposed algorithm is tested using 20 standard benchmark optimization functions with different parameter values. Besides, we utilized five IEEE Congress of Evolutionary Computation (CEC) benchmark test functions (CECC06, 2019 Competition) and five CEC 2011 real-world problems to prove the proposed algorithm's efficiency. The proposed algorithm is compared to eight of the most popular and recent metaheuristic algorithms from the state-of-the-art in terms of best cost, average cost (AVG), corresponding standard deviation (STD), and convergence speed. The results demonstrate that COVIDOA is superior to most existing metaheuristics.

Ranking the effectiveness of non-pharmaceutical interventions to counter COVID-19 in UK universities with vaccinated population

Several universities around the world have resumed in-person teaching after successful vaccination campaigns have covered 70/80% of the population. In this study, we combine a new compartmental model with an optimal control formulation to discover, among different non-pharmaceutical interventions, the best prevention strategy to maximize on-campus activities while keeping spread under control. Composed of two interconnected Susceptible-Exposed-Infected-Quarantined-Recovered (SEIQR) structures, the model enables staff-to-staff infections, student-to-staff cross infections, student-to-student infections, and environment-to-individual infections. Then, we model input variables representing the implementation of different non-pharmaceutical interventions and formulate and solve optimal control problems for four desired scenarios: minimum number of cases, minimum intervention, minimum non-quarantine intervention, and minimum quarantine intervention. Our results reveal the particular significance of mask wearing and social distancing in universities with vaccinated population (with proportions according to UK data). The study also reveals that quarantining infected students has a higher importance than quarantining staff. In contrast, other measures such as environmental disinfection seems to be less important.

Remission of Liquid Tumors and SARS-CoV-2 Infection: a Literature Review

The coronavirus disease 2019 (COVID-19) pandemic has produced a new global challenge for patients with cancer. The disease and the immunosuppression induced by cancer therapies have generated a perfect storm of conditions to increase the severity of the symptoms and worsen the prognosis. However, a few clinical reports showcased the power of viruses to induce remission in some patients suffering from liquid tumors. Here, we reviewed six cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that resulted in cancer remission, simultaneously highlighting the strengths and the unique challenges of oncolytic virotherapy. Virotherapy has become a special case of cancer immunotherapy. This paradigm-shifting concept suggests that oncolytic viruses are not only promising agents to combat particularly immunologically suppressed, immunotherapy-resistant tumors but also that the trigger of local inflammation, such as SARS-CoV-2 infection of the respiratory pathways, may trigger an abscopal effect sufficient to induce the remission of systemic cancer.

Viruses as Living Systems—A Metacybernetic View

The debate over whether viruses are living organisms tends to be paradigmatically determined. The metabolic paradigm denies that they are, while new research evidences the opposite. The purpose of this paper is to deliver a generic model for viral contexts that explains why viruses are alive. It will take a systems biology approach, with a qualitative part (using metacybernetics) to provide deeper explanations of viral contexts, and a quantitative part (using Fisher Information deriving from the variational principle of Extreme Physical Information) which is in principle able to take measurements and predict outcomes. The modelling process provides an extended view of the epigenetic processes of viruses. The generic systems biology model will depict viruses as autonomous entities with metaphysical processes of autopoietic self-organisation and adaptation, enabling them to maintain their physical viability and hence, within their populations, mutate and evolve. The autopoietic epigenetic processes are shown to describe their capability to change, and these are both qualitatively and quantitatively explored, the latter providing an approach to make measurements of physical phenomena under uncertainty. Viruses maintain their fitness when they are able to maintain their stability, and this is indicated by information flow efficacy. A brief case study is presented on the COVID-19 virus from the perspective that it is a living system, and this includes outcome predictions given Fisher Information conditions for known contexts.

The Intrapopulation Genetic Diversity of RNA Virus May Influence the Sensitivity of Chlorine Disinfection

RNA virus populations are not clonal; rather, they comprise a mutant swarm in which sequences are slightly different from the master sequence. Genetic diversity within a population (intrapopulation genetic diversity) is critical for RNA viruses to survive under environmental stresses. Disinfection has become an important practice in the control of pathogenic viruses; however, the impact of intrapopulation genetic diversity on the sensitivity of disinfection, defined as -log10 (postdisinfected infectious titer/predisinfected titer), has not been elucidated. In this study, we serially passaged populations of rhesus rotavirus. We demonstrated that populations with reduced chlorine sensitivity emerged at random and independently of chlorine exposure. Sequencing analysis revealed that compared with sensitive populations, less-sensitive ones had higher non-synonymous genetic diversity of the outer capsid protein gene, suggesting that changes in the amino acid sequences of the outer capsid protein were the main factors influencing chlorine sensitivity. No common mutations were found among less-sensitive populations, indicating that rather than specific mutations, the diversity of the outer capsid protein itself was associated with the disinfection sensitivity and that the disinfection sensitivity changed stochastically. Simulation results suggest that the disinfection sensitivity of a genetically diverse population is destabilized if cooperative viral clusters including multiple sequences are formed. These results advocate that any prevention measures leading to low intrapopulation genetic diversity are important to prevent the spread and evolution of pathogenic RNA viruses in society.

Antiretroviral Drug-Resistance Mutations on the Gag Gene: Mutation Dynamics during Analytic Treatment Interruption among Individuals Experiencing Virologic Failure

We describe drug-resistance mutation dynamics of the gag gene among individuals under antiretroviral virologic failure who underwent analytical treatment interruption (ATI). These mutations occur in and around the cleavage sites that form the particles that become the mature HIV-1 virus. The study involved a 12-week interruption in antiretroviral therapy (ART) and sequencing of the gag gene in 38 individuals experiencing virologic failure and harboring triple-class resistant HIV strains. Regions of the gag gene surrounding the NC-p2 and p1-p6 cleavage sites were sequenced at baseline before ATI and after 12 weeks from plasma HIV RNA using population-based Sanger sequencing. Fourteen of the sixteen patients sequenced presented at least one mutation in the gag gene at baseline, with an average of 4.93 mutations per patient. All the mutations had reverted to the wild type by the end of the study. Mutations in the gag gene complement mutations in the pol gene to restore HIV fitness. Those mutations around cleavage sites and within substrates contribute to protease inhibitor resistance and difficulty in re-establishing effective virologic suppression. ART interruption in the presence of antiretroviral resistant HIV strains was used here as a practical measure for more adapted HIV profiles in the absence of ART selective pressure.

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

Strainline: full-length de novo viral haplotype reconstruction from noisy long reads

Haplotype-resolved de novo assembly of highly diverse virus genomes is critical in prevention, control and treatment of viral diseases. Current methods either can handle only relatively accurate short read data, or collapse haplotype-specific variations into consensus sequence. Here, we present Strainline, a novel approach to assemble viral haplotypes from noisy long reads without a reference genome. Strainline is the first approach to provide strain-resolved, full-length de novo assemblies of viral quasispecies from noisy third-generation sequencing data. Benchmarking on simulated and real datasets of varying complexity and diversity confirm this novelty and demonstrate the superiority of Strainline.

A defective viral genome strategy elicits broad protective immunity against respiratory viruses

RNA viruses generate defective viral genomes (DVGs) that can interfere with replication of the parental wild-type virus. To examine their therapeutic potential, we created a DVG by deleting the capsid-coding region of poliovirus. Strikingly, intraperitoneal or intranasal administration of this genome, which we termed eTIP1, elicits an antiviral response, inhibits replication, and protects mice from several RNA viruses, including enteroviruses, influenza, and SARS-CoV-2. While eTIP1 replication following intranasal administration is limited to the nasal cavity, its antiviral action extends non-cell-autonomously to the lungs. eTIP1 broad-spectrum antiviral effects are mediated by both local and distal type I interferon responses. Importantly, while a single eTIP1 dose protects animals from SARS-CoV-2 infection, it also stimulates production of SARS-CoV-2 neutralizing antibodies that afford long-lasting protection from SARS-CoV-2 reinfection. Thus, eTIP1 is a safe and effective broad-spectrum antiviral generating short- and long-term protection against SARS-CoV-2 and other respiratory infections in animal models.

Microorganisms as Shapers of Human Civilization, from Pandemics to Even Our Genomes: Villains or Friends? A Historical Approach

Universal history is characterized by continuous evolution, in which civilizations are born and die. This evolution is associated with multiple factors, among which the role of microorganisms is often overlooked. Viruses and bacteria have written or decisively contributed to terrible episodes of history, such as the Black Death in 14th century Europe, the annihilation of pre-Columbian American civilizations, and pandemics such as the 1918 Spanish flu or the current COVID-19 pandemic caused by the coronavirus SARS-CoV-2. Nevertheless, it is clear that we could not live in a world without these tiny beings. Endogenous retroviruses have been key to our evolution and for the regulation of gene expression, and the gut microbiota helps us digest compounds that we could not otherwise process. In addition, we have used microorganisms to preserve or prepare food for millennia and more recently to obtain drugs such as antibiotics or to develop recombinant DNA technologies. Due to the enormous importance of microorganisms for our survival, they have significantly influenced the population genetics of different human groups. This paper will review the role of microorganisms as "villains" who have been responsible for tremendous mortality throughout history but also as "friends" who help us survive and evolve.

Simulated vector transmission differentially influences dynamics of two viral variants of deformed wing virus in honey bees (Apis mellifera)

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5′ IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) genetic diversity and occurrence of wild type and vaccine-like strains in the United States swine industry

Porcine reproductive and respiratory syndrome virus genotype 2 (PRRSV-2) genetic diversity in the U.S. was assessed using a database comprising 10 years’ worth of sequence data obtained from swine production systems routine monitoring and outbreak investigations. A total of 26,831 ORF5 PRRSV-2 sequences from 34 production systems were included in this analysis. Within group mean genetic distance (i.e. mean proportion of nucleotide differences within ORF5) per year according to herd type was calculated for all PRRSV-2 sequences. The percent nucleotide difference between each sequence and the ORF5 sequences from four commercially available PRRSV-2 vaccines (Ingelvac PRRS MLV, Ingelvac PRRS ATP, Fostera PRRS, and Prevacent PRRS) within the same lineage over time was used to classify sequences in wild-type or vaccine-like. The mean ORF5 genetic distance fluctuated from 0.09 to 0.13, being generally smaller in years in which there was a relative higher frequency of dominant lineage. Vaccine-like sequences comprised about one fourth of sequences obtained through routine monitoring of PRRS. We found that lineage 5 sequences were mostly Ingelvac PRRS MLV-like. Lineage 8 sequences up to 2011 were 62.9% Ingelvac PRRS ATP-like while the remaining were wild-type viruses. From 2012 onwards, 51.9% of lineage 8 sequences were Ingelvac PRRS ATP-like, 45.0% were Fostera PRRS-like, and only 3.2% were wild-type. For lineage 1 sequences, 0.1% and 1.7% of the sequences were Prevacent PRRS-like in 2009–2018 and 2019, respectively. These results suggest that repeated introductions of vaccine-like viruses through use of modified live vaccines might decrease within-lineage viral diversity as vaccine-like strains become more prevalent. Overall, this compilation of private data from routine monitoring provides valuable information on PRRSV viral diversity.

Emergence of a New Lineage 1C Variant of Porcine Reproductive and Respiratory Syndrome Virus 2 in the United States

We report an ongoing regional outbreak of an emerging porcine reproductive and respiratory syndrome virus (PRRSV2) variant within Lineage 1C affecting 154 breeding and grow-finishing sites in the Midwestern U.S. Transmission seemed to have occurred in two waves, with the first peak of weekly cases occurring between October and December 2020 and the second starting in April 2021. Most of cases occurred within a 120 km radius. Both orf5 and whole genome sequencing results suggest that this represents the emergence of a new variant within Lineage 1C distinct from what has been previously circulating. A case-control study was conducted with 50 cases (sites affected with the newly emerged variant) and 58 controls (sites affected with other PRRSV variants) between October and December 2020. Sites that had a market vehicle that was not exclusive to the production system had 0.04 times the odds of being a case than a control. A spatial cluster (81.42 km radius) with 1.68 times higher the number of cases than controls was found. The average finishing mortality within the first 4 weeks after detection was higher amongst cases (4.50%) than controls (0.01%). The transmission of a highly similar virus between different farms carrying on trough spring rises concerns for the next high transmission season of PRRS.

Fitness Cost of Antiretroviral Drug Resistance Mutations on the pol Gene during Analytical Antiretroviral Treatment Interruption among Individuals Experiencing Virological Failure

HIV cure studies require patients to enter an analytical treatment interruption (ATI). Here, we describe previously unanalyzed data that sheds light on ATI dynamics in PLHIV (People Living with HIV). We present drug resistance mutation dynamics on the pol gene among individuals with antiretroviral virological failure who underwent ATI. The study involved a 12-week interruption in antiretroviral therapy (ART), monitoring of viral load, CD4+/CD8+ T cell counts, and sequencing of the pol gene from 38 individuals experiencing virological failure and harboring 3-class resistant HIV strains: nucleoside reverse transcriptase inhibitors (NRTI) non-nucleoside inhibitors (NNRTI), and protease inhibitors (PI). Protease and reverse transcriptase regions of the pol gene were sequenced at baseline before ATI and every four weeks thereafter from PBMCs and at baseline and after 12 weeks from plasma HIV RNA using population-based Sanger sequencing. Average viral load increased 0.559 log₁₀ copies per milliliter. CD4+ T cell count decreased as soon as ART was withdrawn, an average loss of 99.0 cells/mL. Forty-three percent of the mutations associated with antiretroviral resistance in PBMCs disappeared and fifty-seven percent of the mutations in plasma reverted to wild type, which was less than the 100% reversion expected. In PBMC, the PI mutations reverted more slowly than reverse transcriptase mutations. The patients were projected to need an average of 33.7 weeks for PI to revert compared with 20.9 weeks for NRTI and 19.8 weeks for NNRTI. Mutations in the pol gene can cause virological failure and difficulty in re-establishing effective virological suppression.

Reduced subgenomic RNA expression is a molecular indicator of asymptomatic SARS-CoV-2 infection

Background

It is estimated that up to 80% of infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are asymptomatic and asymptomatic patients can still effectively transmit the virus and cause disease. While much of the effort has been placed on decoding single nucleotide variation in SARS-CoV-2 genomes, considerably less is known about their transcript variation and any correlation with clinical severity in human hosts, as defined here by the presence or absence of symptoms.

Methods

To assess viral genomic signatures of disease severity, we conducted a systematic characterization of SARS-CoV-2 transcripts and genetic variants in 81 clinical specimens collected from symptomatic and asymptomatic individuals using multi-scale transcriptomic analyses including amplicon-seq, short-read metatranscriptome and long-read Iso-seq.

Results

Here we show a highly coordinated and consistent pattern of sgRNA expression from individuals with robust SARS-CoV-2 symptomatic infection and their expression is significantly repressed in the asymptomatic infections. We also observe widespread inter- and intra-patient variants in viral RNAs, known as quasispecies frequently found in many RNA viruses. We identify unique sets of deletions preferentially found primarily in symptomatic individuals, with many likely to confer changes in SARS-CoV-2 virulence and host responses. Moreover, these frequently occurring structural variants in SARS-CoV-2 genomes serve as a mechanism to further induce SARS-CoV-2 proteome complexity.

Conclusions

Our results indicate that differential sgRNA expression and structural mutational burden are highly correlated with the clinical severity of SARS-CoV-2 infection. Longitudinally monitoring sgRNA expression and structural diversity could further guide treatment responses, testing strategies, and vaccine development.

Wong and Ngan et al. characterize the expression and structural variation of SARS-CoV-2 subgenomic RNAs (sgRNAs) in diagnostic specimens from symptomatic and asymptomatic COVID-19 patients. In a series of genomic and transcriptomic analyses, the authors observe reduced sgRNA expression and a distinct set of structural deletions in asymptomatic infections compared with symptomatic infections.

Infection with SARS-CoV-2 can lead to symptoms of different severity, with some individuals remaining entirely asymptomatic but still responsible for much of the viral transmission. Here, we sought to identify markers of the severity of symptoms of SARS-CoV-2 infection, as defined by the presence or absence of symptoms. We used a combination of methods to study SARS-CoV-2 genes and their readouts, known as transcripts, in clinical samples from people with symptomatic and asymptomatic infections. We demonstrate that transcripts responsible for making viral proteins are seen at lower levels in asymptomatic infections. We also identify structural changes in the viral genes and transcripts that potentially influence the host’s response to the virus. Our study defines potential markers of symptom severity that may ultimately guide risk mitigation and testing strategies.

Towards an integrative view of virus phenotypes

Understanding how phenotypes emerge from genotypes is a foundational goal in biology. As challenging as this task is when considering cellular life, it is further complicated in the case of viruses. During replication, a virus as a discrete entity (the virion) disappears and manifests itself as a metabolic amalgam between the virus and the host (the virocell). Identifying traits that unambiguously constitute a virus's phenotype is straightforward for the virion, less so for the virocell. Here, we present a framework for categorizing virus phenotypes that encompasses both virion and virocell stages and considers functional and performance traits of viruses in the context of fitness. Such an integrated view of virus phenotype is necessary for comprehensive interpretation of viral genome sequences and will advance our understanding of viral evolution and ecology.

dsmCRISPR: Dual synthetic mismatches CRISPR/Cas12a-based detection of SARS-CoV-2 D614G mutation

Fast evolving of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus has caused the spreading of COVID-19 disease rapidly around the globe. The mutation, especially in the gene encoding spike protein has helped the virus adapt and evade human immune system, as well as affect the efficacy of the immunizations and treatments. SARS-CoV-2 variant carrying D614G amino acid change at the spike protein is the most dominant strain in the pandemic. Therefore, efficient detection of the SARS-CoV-2 variants including D614G mutation is critical to control the COVID-19 pandemic. Herein, we report a dual synthetic mismatches CRISPR/Cas12a (dsmCRISPR) method to detect the SARS-CoV-2 D614G mutation with high sensitivity and specificity. By targeting SARS-CoV-2 D614G mutation, synthetic mismatch crRNAs were designed from -3 to +3 position around the mutation site. To improve the sensitivity and specificity, a synthetic mismatch primer with a 3'-terminal base complementary to the D614G point mutation and a mismatch next to 3'-terminal base was used to specifically amplify the D614G mutation site with higher annealing temperature. Using synthetic mismatch crRNA-(-1), a higher ratio (13.45) of the fluorescence between G614 and D614 was observed. When combined with mismatch primer to amplify D614G mutation, the fluorescence ratio of G614/D164 template detected was increased by 73.53% to 23.12. This method can detect the SARS-CoV-2 D614G mutation nucleic acid with high sensitivity, which was validated with synthetic SARS-CoV-2 D614G RNA. Therefore, the dsmCRIPSR method has significant potential to serve as a sensitive and specific assay for SARS-CoV-2 D614G detection and could be further extended for the detection of other SARS-CoV-2 variants of interest.

Antiviral therapeutics directed against RNA dependent RNA polymerases from positive-sense viruses

Viruses with positive-sense single stranded RNA (+ssRNA) genomes are responsible for different diseases and represent a global health problem. In addition to developing new vaccines that protect against severe illness on infection, it is imperative to identify new antiviral molecules to treat infected patients. The genome of these RNA viruses generally codes for an enzyme with RNA dependent RNA polymerase (RdRP) activity. This molecule is centrally involved in the duplication of the RNA genome. Inhibition of this enzyme by small molecules will prevent duplication of the RNA genome and thus reduce the viral titer. An overview of the different therapeutic strategies used to inhibit RdRPs from +ssRNA viruses is provided, along with an analysis of these enzymes to highlight new binding sites for inhibitors.

tRNA-like structures and their functions

tRNA-like structures (TLSs) were first identified in the RNA genomes of turnip yellow mosaic virus. Since then, TLSs have been found in many other species including mammals, and the RNAs harboring these structures range from viral genomic RNAs to mRNAs and noncoding RNAs. Some progress has also been made on understanding their functions that include regulation of RNA replication, translation enhancement, RNA-protein interaction, and more. In this review, we summarize the current knowledge about the regulations and functions of these TLSs. Possible future directions of the field are also briefly discussed.

Genetic diversity through social heterosis can increase virulence in RNA viral infections and cancer progression

In viral infections and cancer tumours, negative health outcomes often correlate with increasing genetic diversity. Possible evolutionary processes for such relationships include mutant lineages escaping host control or diversity, per se, creating too many immune system targets. Another possibility is social heterosis where mutations and replicative errors create clonal lineages varying in intrinsic capability for successful dispersal; improved environmental buffering; resource extraction or effective defence against immune systems. Rather than these capabilities existing in one genome, social heterosis proposes complementary synergies occur across lineages in close proximity. Diverse groups overcome host defences as interacting 'social genomes' with group genetic tool kits exceeding limited individual plasticity. To assess the possibility of social heterosis in viral infections and cancer progression, we conducted extensive literature searches for examples consistent with general and specific predictions from the social heterosis hypothesis. Numerous studies found supportive patterns in cancers across multiple tissues and in several families of RNA viruses. In viruses, social heterosis mechanisms probably result from long coevolutionary histories of competition between pathogen and host. Conversely, in cancers, social heterosis is a by-product of recent mutations. Investigating how social genomes arise and function in viral quasi-species swarms and cancer tumours may lead to new therapeutic approaches.

Broad-Spectrum Antiviral Strategies and Nucleoside Analogues

The emergence or re-emergence of viruses with epidemic and/or pandemic potential, such as Ebola, Zika, Middle East Respiratory Syndrome (MERS-CoV), Severe Acute Respiratory Syndrome Coronavirus 1 and 2 (SARS and SARS-CoV-2) viruses, or new strains of influenza represents significant human health threats due to the absence of available treatments. Vaccines represent a key answer to control these viruses. However, in the case of a public health emergency, vaccine development, safety, and partial efficacy concerns may hinder their prompt deployment. Thus, developing broad-spectrum antiviral molecules for a fast response is essential to face an outbreak crisis as well as for bioweapon countermeasures. So far, broad-spectrum antivirals include two main categories: the family of drugs targeting the host-cell machinery essential for virus infection and replication, and the family of drugs directly targeting viruses. Among the molecules directly targeting viruses, nucleoside analogues form an essential class of broad-spectrum antiviral drugs. In this review, we will discuss the interest for broad-spectrum antiviral strategies and their limitations, with an emphasis on virus-targeted, broad-spectrum, antiviral nucleoside analogues and their mechanisms of action.

Uniqueness of RNA Coliphage Qβ Display System in Directed Evolutionary Biotechnology

Phage display technology involves the surface genetic engineering of phages to expose desirable proteins or peptides whose gene sequences are packaged within phage genomes, thereby rendering direct linkage between genotype with phenotype feasible. This has resulted in phage display systems becoming invaluable components of directed evolutionary biotechnology. The M13 is a DNA phage display system which dominates this technology and usually involves selected proteins or peptides being displayed through surface engineering of its minor coat proteins. The displayed protein or peptide’s functionality is often highly reduced due to harsh treatment of M13 variants. Recently, we developed a novel phage display system using the coliphage Qβ as a nano-biotechnology platform. The coliphage Qβ is an RNA phage belonging to the family of Leviviridae, a long investigated virus. Qβ phages exist as a quasispecies and possess features making them comparatively more suitable and unique for directed evolutionary biotechnology. As a quasispecies, Qβ benefits from the promiscuity of its RNA dependent RNA polymerase replicase, which lacks proofreading activity, and thereby permits rapid variant generation, mutation, and adaptation. The minor coat protein of Qβ is the readthrough protein, A₁. It shares the same initiation codon with the major coat protein and is produced each time the ribosome translates the UGA stop codon of the major coat protein with the of misincorporation of tryptophan. This misincorporation occurs at a low level (1/15). Per convention and definition, A₁ is the target for display technology, as this minor coat protein does not play a role in initiating the life cycle of Qβ phage like the pIII of M13. The maturation protein A₂ of Qβ initiates the life cycle by binding to the pilus of the F⁺ host bacteria. The extension of the A₁ protein with a foreign peptide probe recognizes and binds to the target freely, while the A₂ initiates the infection. This avoids any disturbance of the complex and the necessity for acidic elution and neutralization prior to infection. The combined use of both the A₁ and A₂ proteins of Qβ in this display system allows for novel bio-panning, in vitro maturation, and evolution. Additionally, methods for large library size construction have been improved with our directed evolutionary phage display system. This novel phage display technology allows 12 copies of a specific desired peptide to be displayed on the exterior surface of Qβ in uniform distribution at the corners of the phage icosahedron. Through the recently optimized subtractive bio-panning strategy, fusion probes containing up to 80 amino acids altogether with linkers, can be displayed for target selection. Thus, combined uniqueness of its genome, structure, and proteins make the Qβ phage a desirable suitable innovation applicable in affinity maturation and directed evolutionary biotechnology. The evolutionary adaptability of the Qβ phage display strategy is still in its infancy. However, it has the potential to evolve functional domains of the desirable proteins, glycoproteins, and lipoproteins, rendering them superior to their natural counterparts.

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.

Inhibition of viral RNA-dependent RNA polymerases with clinically relevant nucleotide analogs

The treatment of viral infections remains challenging, in particular in the face of emerging pathogens. Broad-spectrum antiviral drugs could potentially be used as a first line of defense. The RNA-dependent RNA polymerase (RdRp) of RNA viruses serves as a logical target for drug discovery and development efforts. Herein we discuss compounds that target RdRp of poliovirus, hepatitis C virus, influenza viruses, respiratory syncytial virus, and the growing data on coronaviruses. We focus on nucleotide analogs and mechanisms of action and resistance.

Polymerase Fidelity Contributes to Foot-and-Mouth Disease Virus Pathogenicity and Transmissibility In Vivo

The low fidelity of foot-and-mouth disease virus (FMDV) RNA-dependent RNA polymerase allows FMDV to exhibit high genetic diversity. Previously, we showed that the genetic diversity of FMDV plays an important role in virulence in suckling mice. Here, we mutated the amino acid residue Phe257, located in the finger domain of FMDV polymerase and conserved across FMDV serotypes, to a cysteine (F257C) to study the relationship between viral genetic diversity, virulence, and transmissibility in natural hosts. The single amino acid substitution in FMDV polymerase resulted in a high-fidelity virus variant, rF257C, with growth kinetics indistinguishable from those of wild-type (WT) virus in cell culture, but it displayed smaller plaques and impaired fitness in direct competition assays. Furthermore, we found that rF257C was attenuated in vivo in both suckling mice and pigs (one of its natural hosts). Importantly, contact exposure experiments showed that the rF257C virus exhibited reduced transmissibility compared to that of wild-type FMDV in the porcine model. This study provides evidence that FMDV genetic diversity is important for viral virulence and transmissibility in susceptible animals. Given that type O FMDV exhibits the highest genetic diversity among all seven serotypes of FMDV, we propose that the lower polymerase fidelity of the type O FMDV could contribute to its dominance worldwide.IMPORTANCE Among the seven serotypes of FMDV, serotype O FMDV have the broadest distribution worldwide, which could be due to their high virulence and transmissibility induced by high genetic diversity. In this paper, we generated a single amino acid substitution FMDV variant with a high-fidelity polymerase associated with viral fitness, virulence, and transmissibility in a natural host. The results highlight that maintenance of viral population diversity is essential for interhost viral spread. This study provides evidence that higher genetic diversity of type O FMDV could increase both virulence and transmissibility, thus leading to their dominance in the global epidemic.

Naturally occurring SARS-CoV-2 gene deletions close to the spike S1/S2 cleavage site in the viral quasispecies of COVID19 patients

The SARS-CoV-2 spike (S) protein, the viral mediator for binding and entry into the host cell, has sparked great interest as a target for vaccine development and treatments with neutralizing antibodies. Initial data suggest that the virus has low mutation rates, but its large genome could facilitate recombination, insertions, and deletions, as has been described in other coronaviruses. Here, we deep-sequenced the complete SARS-CoV-2 S gene from 18 patients (10 with mild and 8 with severe COVID-19), and found that the virus accumulates deletions upstream and very close to the S1/S2 cleavage site (PRRAR/S), generating a frameshift with appearance of a stop codon. These deletions were found in a small percentage of the viral quasispecies (2.2%) in samples from all the mild and only half the severe COVID-19 patients. Our results suggest that the virus may generate free S1 protein released to the circulation. We suggest that natural selection has favoured a "Don't burn down the house" strategy, in which free S1 protein may compete with viral particles for the ACE2 receptor, thus reducing the severity of the infection and tissue damage without losing transmission capability.

Papain-Like Proteases as Coronaviral Drug Targets: Current Inhibitors, Opportunities, and Limitations

Papain-like proteases (PLpro) of coronaviruses (CoVs) support viral reproduction and suppress the immune response of the host, which makes CoV PLpro perspective pharmaceutical targets. Their inhibition could both prevent viral replication and boost the immune system of the host, leading to the speedy recovery of the patient. Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is the third CoV outbreak in the last 20 years. Frequent mutations of the viral genome likely lead to the emergence of more CoVs. Inhibitors for CoV PLpro can be broad-spectrum and can diminish present and prevent future CoV outbreaks as PLpro from different CoVs have conservative structures. Several inhibitors have been developed to withstand SARS-CoV and Middle East respiratory syndrome CoV (MERS-CoV). This review summarizes the structural features of CoV PLpro, the inhibitors that have been identified over the last 20 years, and the compounds that have the potential to become novel effective therapeutics against CoVs in the near future.

Viruses That Can and Cannot Coexist With Humans and the Future of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a worldwide pandemic. Many projections concerning the outbreak, such as the estimated number of cases and deaths in upcoming months, have been made available. However, what happens to the virus after the pandemic subsides has not been fully explored. In this article, we discuss the ways that past and present human viruses have emerged via zoonotic transmission, the mechanisms that they have acquired the ability for effective transmission among humans, the process to sustain a chain of transmission to coexist with humans, and the factors important for complete containment leading to eradication of viruses. These aspects of viral disease may provide clues for the future path that SARS-CoV-2 might take in relation to human infection.

Bottleneck Size-Dependent Changes in the Genetic Diversity and Specific Growth Rate of a Rotavirus A Strain

RNA viruses form a dynamic distribution of mutant swarms (termed "quasispecies") due to the accumulation of mutations in the viral genome. The genetic diversity of a viral population is affected by several factors, including a bottleneck effect. Human-to-human transmission exemplifies a bottleneck effect, in that only part of a viral population can reach the next susceptible hosts. In the present study, two lineages of the rhesus rotavirus (RRV) strain of rotavirus A were serially passaged five times at a multiplicity of infection (MOI) of 0.1 or 0.001, and three phenotypes (infectious titer, cell binding ability, and specific growth rate) were used to evaluate the impact of a bottleneck effect on the RRV population. The specific growth rate values of lineages passaged under the stronger bottleneck (MOI of 0.001) were higher after five passages. The nucleotide diversity also increased, which indicated that the mutant swarms of the lineages under the stronger bottleneck effect were expanded through the serial passages. The random distribution of synonymous and nonsynonymous substitutions on rotavirus genome segments indicated that almost all mutations were selectively neutral. Simple simulations revealed that the presence of minor mutants could influence the specific growth rate of a population in a mutant frequency-dependent manner. These results indicate a stronger bottleneck effect can create more sequence spaces for minor sequences.IMPORTANCE In this study, we investigated a bottleneck effect on an RRV population that may drastically affect the viral population structure. RRV populations were serially passaged under two levels of a bottleneck effect, which exemplified human-to-human transmission. As a result, the genetic diversity and specific growth rate of RRV populations increased under the stronger bottleneck effect, which implied that a bottleneck created a new space in a population for minor mutants originally existing in a hidden layer, which includes minor mutations that cannot be distinguished from a sequencing error. The results of this study suggest that the genetic drift caused by a bottleneck in human-to-human transmission explains the random appearance of new genetic lineages causing viral outbreaks, which can be expected according to molecular epidemiology using next-generation sequencing in which the viral genetic diversity within a viral population is investigated.

Limited evolution of the yellow fever virus 17d in a mouse infection model

By infecting mice with the yellow fever virus vaccine strain 17D (YFV-17D; Stamaril®), the dose dependence and evolutionary consequences of neurotropic yellow fever infection was assessed. Highly susceptible AG129 mice were used to allow for a maximal/unlimited expansion of the viral populations. Infected mice uniformly developed neurotropic disease; the virus was isolated from their brains, plaque purified and sequenced. Viral RNA populations were overall rather homogenous [Shannon entropies 0−0.15]. The remaining, yet limited intra-host population diversity (0−11 nucleotide exchanges per genome) appeared to be a consequence of pre-existing clonal heterogeneities (quasispecies) of Stamaril®. In parallel, mice were infected with a molecular clone of YFV-17D which was in vivo launched from a plasmid. Such plasmid-launched YFV-17D had a further reduced and almost clonal evolution. The limited intra-host evolution during unrestricted expansion in a highly susceptible host is relevant for vaccine and drug development against flaviviruses in general. Firstly, a propensity for limited evolution even upon infection with a (very) low inoculum suggests that fractional dosing as implemented in current YF-outbreak control may pose only a limited risk of reversion to pathogenic vaccine-derived virus variants. Secondly, it also largely lowers the chance of antigenic drift and development of resistance to antivirals.

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.

Detection of Plant Viruses and Disease Management: Relevance of Genetic Diversity and Evolution

Plant viruses cause considerable economic losses and are a threat for sustainable agriculture. The frequent emergence of new viral diseases is mainly due to international trade, climate change, and the ability of viruses for rapid evolution. Disease control is based on two strategies: i) immunization (genetic resistance obtained by plant breeding, plant transformation, cross-protection, or others), and ii) prophylaxis to restrain virus dispersion (using quarantine, certification, removal of infected plants, control of natural vectors, or other procedures). Disease management relies strongly on a fast and accurate identification of the causal agent. For known viruses, diagnosis consists in assigning a virus infecting a plant sample to a group of viruses sharing common characteristics, which is usually referred to as species. However, the specificity of diagnosis can also reach higher taxonomic levels, as genus or family, or lower levels, as strain or variant. Diagnostic procedures must be optimized for accuracy by detecting the maximum number of members within the group (sensitivity as the true positive rate) and distinguishing them from outgroup viruses (specificity as the true negative rate). This requires information on the genetic relationships within-group and with members of other groups. The influence of the genetic diversity of virus populations in diagnosis and disease management is well documented, but information on how to integrate the genetic diversity in the detection methods is still scarce. Here we review the techniques used for plant virus diagnosis and disease control, including characteristics such as accuracy, detection level, multiplexing, quantification, portability, and designability. The effect of genetic diversity and evolution of plant viruses in the design and performance of some detection and disease control techniques are also discussed. High-throughput or next-generation sequencing provides broad-spectrum and accurate identification of viruses enabling multiplex detection, quantification, and the discovery of new viruses. Likely, this technique will be the future standard in diagnostics as its cost will be dropping and becoming more affordable.

Temporal Dynamics of Co-circulating Lineages of Porcine Reproductive and Respiratory Syndrome Virus

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) is the most important endemic pathogen in the U.S. swine industry. Despite control efforts involving improved biosecurity and different vaccination protocols, the virus continues to circulate and evolve. One of the foremost challenges in its control is high levels of genetic and antigenic diversity. Here, we quantify the co-circulation, emergence and sequential turnover of multiple PRRSV lineages in a single swine-producing region in the United States over a span of 9 years (2009–2017). By classifying over 4,000 PRRSV sequences (open-reading frame 5) into phylogenetic lineages and sub-lineages, we document the ongoing diversification and temporal dynamics of the PRRSV population, including the rapid emergence of a novel sub-lineage that appeared to be absent globally pre-2008. In addition, lineage 9 was the most prevalent lineage from 2009 to 2010, but its occurrence fell to 0.5% of all sequences identified per year after 2014, coinciding with the emergence or re-emergence of lineage 1 as the dominant lineage. The sequential dominance of different lineages, as well as three different sub-lineages within lineage 1, is consistent with the immune-mediated selection hypothesis for the sequential turnover in the dominant lineage. As host populations build immunity through natural infection or vaccination toward the most common variant, this dominant (sub-) lineage may be replaced by an emerging variant to which the population is more susceptible. An analysis of patterns of non- synonymous and synonymous mutations revealed evidence of positive selection on immunologically important regions of the genome, further supporting the potential that immune-mediated selection shapes the evolutionary and epidemiological dynamics for this virus. This has important implications for patterns of emergence and re-emergence of genetic variants of PRRSV that have negative impacts on the swine industry. Constant surveillance on PRRSV occurrence is crucial to a better understanding of the epidemiological and evolutionary dynamics of co-circulating viral lineages. Further studies utilizing whole genome sequencing and exploring the extent of cross-immunity between heterologous PRRS viruses could shed further light on PRRSV immunological response and aid in developing strategies that might be able to diminish disease impact.

Tyr82 Amino Acid Mutation in PB1 Polymerase Induces an Influenza Virus Mutator Phenotype

In various positive-sense single-stranded RNA viruses, a low-fidelity viral RNA-dependent RNA polymerase (RdRp) confers attenuated phenotypes by increasing the mutation frequency. We report a negative-sense single-stranded RNA virus RdRp mutant strain with a mutator phenotype. Based on structural data of RdRp, rational targeting of key residues, and screening of fidelity variants, we isolated a novel low-fidelity mutator strain of influenza virus that harbors a Tyr82-to-Cys (Y82C) single-amino-acid substitution in the PB1 polymerase subunit. The purified PB1-Y82C polymerase indeed showed an increased frequency of misincorporation compared with the wild-type PB1 in an in vitro biochemical assay. To further investigate the effects of position 82 on PB1 polymerase fidelity, we substituted various amino acids at this position. As a result, we isolated various novel mutators other than PB1-Y82C with higher mutation frequencies. The structural model of influenza virus polymerase complex suggested that the Tyr82 residue, which is located at the nucleoside triphosphate entrance tunnel, may influence a fidelity checkpoint. Interestingly, although the PB1-Y82C variant replicated with wild-type PB1-like kinetics in tissue culture, the 50% lethal dose of the PB1-Y82C mutant was 10 times lower than that of wild-type PB1 in embryonated chicken eggs. In conclusion, our data indicate that the Tyr82 residue of PB1 has a crucial role in regulating polymerase fidelity of influenza virus and is closely related to attenuated pathogenic phenotypes in vivo IMPORTANCE Influenza A virus rapidly acquires antigenic changes and antiviral drug resistance, which limit the effectiveness of vaccines and drug treatments, primarily owing to its high rate of evolution. Virus populations formed by quasispecies can contain resistance mutations even before a selective pressure is applied. To study the effects of the viral mutation spectrum and quasispecies, high- and low-fidelity variants have been isolated for several RNA viruses. Here, we report the discovery of a low-fidelity RdRp variant of influenza A virus that contains a substitution at Tyr82 in PB1. Viruses containing the PB1-Y82C substitution showed growth kinetics and viral RNA synthesis levels similar to those of the wild-type virus in cell culture; however, they had significantly attenuated phenotypes in a chicken egg infection experiment. These data demonstrated that decreased RdRp fidelity attenuates influenza A virus in vivo, which is a desirable feature for the development of safer live attenuated vaccine candidates.

Predictive Water Virology: Hierarchical Bayesian Modeling for Estimating Virus Inactivation Curve

Hazard analysis and critical control point (HACCP) are a series of actions to be taken to ensure product consumption safety. In food poisoning risk management, researchers in the field of predictive microbiology calculate the values that provide minimum stress (e.g., temperature and contact time in heating) for sufficient microbe inactivation based on mathematical models. HACCP has also been employed for health risk management in sanitation safety planning (SSP), but the application of predictive microbiology to water-related pathogens is difficult because the variety of pathogen types and the complex composition of the wastewater matrix does not allow us to make a simple mathematical model to predict inactivation efficiency. In this study, we performed a systematic review and meta-analysis to construct predictive inactivation curves using free chlorine for enteric viruses based on a hierarchical Bayesian model using parameters such as water quality. Our model considered uncertainty among virus disinfection tests and difference in genotype-dependent sensitivity of a virus to disinfectant. The proposed model makes it possible to identify critical disinfection stress capable of reducing virus concentration that is below the tolerable concentration to ensure human health.