Proofreading-Deficient Coronaviruses Adapt for Increased Fitness over Long-Term Passage without Reversion of Exoribonuclease-Inactivating Mutations

Eight-amino-acid sequence at the N-terminus of SARS-CoV-2 nsp1 is involved in stabilizing viral genome replication

Coronavirus disease 19 is caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) enveloped virus with a single-stranded positive-sense ribonucleic acid (RNA) genome. The CoV non-structural protein (nsp) 1 is a multifunctional protein that undergoes translation shutoff, messenger RNA (mRNA) cleavage, and RNA binding. The C-terminal region is involved in translational shutoff and RNA cleavage. The N-terminal region of SARS-CoV-2 nsp1 is highly conserved among isolated SARS-CoV-2 variants. However, the I-004 variant, isolated during the early SARS-CoV-2 pandemic, lost eight amino acids in the nsp1 region. In this study, we showed that the eight amino acids are important for viral replication in infected interferon-incompetent cells and that the recombinant virus that lost these amino acids had low pathogenicity in the lungs of hamster models. The loss of eight amino acids-induced mutations occurred in the 5' untranslated region (UTR), suggesting that nsp1 contributes to the stability of the viral genome during replication.

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV)

We analyzed the spike protein S1/S2 cleavage of selected strains of a prototype coronavirus, mouse hepatitis virus (MHV) by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of selected MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results-but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

A Proofreading Mutation with an Allosteric Effect Allows a Cluster of SARS-CoV-2 Viruses to Rapidly Evolve

The RNA-dependent RNA polymerase of the severe acute respiratory syndrome coronavirus 2 virus is error prone, with errors being corrected by the exonuclease (NSP14) proofreading mechanism. However, the mutagenesis and subsequent evolutionary trajectory of the virus is mediated by the delicate interplay of replicase fidelity and environmental pressures. Here, we have shown that a single, distal mutation (F60S) in NSP14 can have a profound impact upon proofreading with an increased accumulation of mutations and elevated evolutionary rate being observed. Understanding the implications of these changes is crucial, as these underlying mutational processes may have important implications for understanding the population-wide evolution of the virus. This study underscores the urgent need for continued research into the replicative mechanisms of this virus to combat its continued impact on global health, through the re-emergence of immuno-evasive variants.

SARS-CoV-2 entry into and evolution within a skilled nursing facility

SARS-CoV-2 belongs to the family Coronaviridae which includes multiple human pathogens that have an outsized impact on aging populations. As a novel human pathogen, SARS-CoV-2 is undergoing continuous adaptation to this new host species and there is evidence of this throughout the scientific and public literature. However, most investigations of SARS-CoV-2 evolution have focused on large-scale collections of data across diverse populations and/or living environments. Here we investigate SARS-CoV-2 evolution in epidemiologically linked individuals within a single outbreak at a skilled nursing facility beginning with initial introduction of the pathogen. The data demonstrate that SARS-CoV-2 was introduced to the facility multiple times without establishing an interfacility transmission chain, followed by a single introduction that infected many individuals within a week. This large-scale introduction by a single genotype then persisted in the facility. SARS-CoV-2 sequences were investigated at both the consensus and intra-host variation levels. Understanding the variability in SARS-CoV-2 during transmission chains will assist in understanding the spread of this disease and can ultimately inform best practices for mitigation strategies.

CRISPR techniques and potential for the detection and discrimination of SARS-CoV-2 variants of concern

The continuing evolution of the SARS-CoV-2 virus has led to the emergence of many variants, including variants of concern (VOCs). CRISPR-Cas systems have been used to develop techniques for the detection of variants. These techniques have focused on the detection of variant-specific mutations in the spike protein gene of SARS-CoV-2. These sequences mostly carry single-nucleotide mutations and are difficult to differentiate using a single CRISPR-based assay. Here we discuss the specificity of the Cas9, Cas12, and Cas13 systems, important considerations of mutation sites, design of guide RNA, and recent progress in CRISPR-based assays for SARS-CoV-2 variants. Strategies for discriminating single-nucleotide mutations include optimizing the position of mismatches, modifying nucleotides in the guide RNA, and using two guide RNAs to recognize the specific mutation sequence and a conservative sequence. Further research is needed to confront challenges in the detection and differentiation of variants and sublineages of SARS-CoV-2 in clinical diagnostic and point-of-care applications.

Genomic diversity of SARS-CoV-2 can be accelerated by mutations in the nsp14 gene

Coronaviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), encode a proofreading exonuclease, nonstructural protein 14 (nsp14), that helps ensure replication competence at a low evolutionary rate compared with other RNA viruses. In the current pandemic, SARS-CoV-2 has accumulated diverse genomic mutations including in nsp14. Here, to clarify whether amino acid substitutions in nsp14 affect the genomic diversity and evolution of SARS-CoV-2, we searched for amino acid substitutions in nature that may interfere with nsp14 function. We found that viruses carrying a proline-to-leucine change at position 203 (P203L) have a high evolutionary rate and that a recombinant SARS-CoV-2 virus with the P203L mutation acquired more diverse genomic mutations than wild-type virus during its replication in hamsters. Our findings suggest that substitutions, such as P203L, in nsp14 may accelerate the genomic diversity of SARS-CoV-2, contributing to virus evolution during the pandemic.

COVID-19: A state of art on immunological responses, mutations, and treatment modalities in riposte

Over the last few years, the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) unleashed a global public health catastrophe that had a substantial influence on human physical and mental health, the global economy, and socio-political dynamics. SARS-CoV-2 is a respiratory pathogen and the cause of ongoing COVID-19 pandemic, which testified how unprepared humans are for pandemics. Scientists and policymakers continue to face challenges in developing ideal therapeutic agents and vaccines, while at the same time deciphering the pathology and immunology of SARS-CoV-2. Challenges in the early part of the pandemic included the rapid development of diagnostic assays, vaccines, and therapeutic agents. The ongoing transmission of COVID-19 is coupled with the emergence of viral variants that differ in their transmission efficiency, virulence, and vaccine susceptibility, thus complicating the spread of the pandemic. Our understanding of how the human immune system responds to these viruses as well as the patient groups (such as the elderly and immunocompromised individuals) who are often more susceptible to serious illness have both been aided by this epidemic. COVID-19 causes different symptoms to occur at different stages of infection, making it difficult to determine distinct treatment regimens employed for the various clinical phases of the disease. Unsurprisingly, determining the efficacy of currently available medications and developing novel therapeutic strategies have been a process of trial and error. The global scientific community collaborated to research and develop vaccines at a neck-breaking speed. This review summarises the overall picture of the COVID-19 pandemic, different mutations in SARS-CoV-2, immune response, and the treatment modalities against SARS-CoV-2.

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV)

We have analyzed the spike protein S1/S2 cleavage site of selected strains of MHV by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of the various MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results-but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We thus identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

Cellular APOBEC3A deaminase drives mutations in the SARS-CoV-2 genome

The number of genetic variations in the SARS-CoV-2 genome has been increasing primarily due to continuous viral mutations. Here, we report that the human APOBEC3A (A3A) cytidine deaminase plays a critical role in the induction of C-to-U substitutions in the SARS-CoV-2 genome. Bioinformatic analysis of the chronological genetic changes in a sequence database indicated that the largest UC-to-UU mutation signature, consistent with APOBEC-recognized nucleotide motifs, was predominant in single-stranded RNA regions of the viral genome. In SARS-CoV-2-infected cells, exogenous expression of A3A but not expression of other APOBEC proteins induced UC-to-UU mutations in viral RNA (vRNA). Additionally, the mutated C bases were often located at the tips in bulge or loop regions in the vRNA secondary structure. Interestingly, A3A mRNA expression was drastically increased by interferons (IFNs) and tumour necrosis factor-α (TNF-α) in epithelial cells derived from the respiratory system, a site of efficient SARS-CoV-2 replication. Moreover, the UC-to-UU mutation rate was increased in SARS-CoV-2 produced from lung epithelial cells treated with IFN-ß and TNF-α, but not from CRISPR/Cas9-based A3A knockout cells. Collectively, these findings demonstrate that A3A is a primary host factor that drives mutations in the SARS-CoV-2 RNA genome via RNA editing.

Mechanisms of Coronavirus Genome Stability As Potential Targets for Antiviral Drugs

The COVID-19 pandemic has made it necessary to create antivirals active against the SARS-CoV-2 coronavirus. One of the widely used strategies to fight off viral infections is the use of modified nucleoside analogues that inhibit viral replication by incorporating DNA or RNA into the growing chain, thus stopping its synthesis. The difficulty of using this method of treatment in the case of SARS-CoV-2 is that coronaviruses have an effective mechanism for maintaining genome stability. Its central element is the nsp14 protein, which is characterized by exonuclease activity, due to which incorrectly included and noncanonical nucleotides are removed from the 3' end of the growing RNA chain. Inhibitors of nsp14 exonuclease and nucleoside analogues resistant to its action are viewed as potential targets for anticoronavirus therapy.

Virus Diversity, Abundance, and Evolution in Three Different Bat Colonies in Switzerland

Bats are increasingly recognized as reservoirs for many different viruses that threaten public health, such as Hendravirus, Ebolavirus, Nipahvirus, and SARS- and MERS-coronavirus. To assess spillover risk, viromes of bats from different parts of the world have been investigated in the past. As opposed to most of these prior studies, which determined the bat virome at a single time point, the current work was performed to monitor changes over time. Specifically, fecal samples of three endemic Swiss bat colonies consisting of three different bat species were collected over three years and analyzed using next-generation sequencing. Furthermore, single nucleotide variants of selected DNA and RNA viruses were analyzed to investigate virus genome evolution. In total, sequences of 22 different virus families were found, of which 13 are known to infect vertebrates. Most interestingly, in a Vespertilio murinus colony, sequences from a MERS-related beta-coronavirus were consistently detected over three consecutive years, which allowed us to investigate viral genome evolution in a natural reservoir host.

Proteolytic Processing of the Coronavirus Replicase Nonstructural Protein 14 Exonuclease Is Not Required for Virus Replication but Alters RNA Synthesis and Viral Fitness

Coronaviruses (CoVs) initiate replication by translation of the positive-sense RNA genome into the replicase polyproteins connecting 16 nonstructural protein domains (nsp1-16), which are subsequently processed by viral proteases to yield mature nsp. For the betacoronavirus murine hepatitis virus (MHV), total inhibition of translation or proteolytic processing of replicase polyproteins results in rapid cessation of RNA synthesis. The nsp5-3CLpro (Mpro) processes nsps7-16, which assemble into functional replication-transcription complexes (RTCs), including the enzymatic nsp12-RdRp and nsp14-exoribonuclease (ExoN)/N7-methyltransferase. The nsp14-ExoN activity mediates RNA-dependent RNA proofreading, high-fidelity RNA synthesis, and replication. To date, the solved partial RTC structures, biochemistry, and models use or assume completely processed, mature nsp. Here, we demonstrate that in MHV, engineered deletion of the cleavage sites between nsp13-14 and nsp14-15 allowed recovery of replication-competent virus. Compared to wild-type (WT) MHV, the nsp13-14 and nsp14-15 cleavage deletion mutants demonstrated delayed replication kinetics, impaired genome production, altered abundance and patterns of recombination, and impaired competitive fitness. Further, the nsp13-14 and nsp14-15 mutant viruses demonstrated mutation frequencies that were significantly higher than with the WT. The results demonstrate that cleavage of nsp13-14 or nsp14-15 is not required for MHV viability and that functions of the RTC/nsp14-ExoN are impaired when assembled with noncleaved intermediates. These data will inform future genetic, structural, biochemical, and modeling studies of coronavirus RTCs and nsp 13, 14, and 15 and may reveal new approaches for inhibition or attenuation of CoV infection.

IMPORTANCE Coronavirus replication requires proteolytic maturation of the nonstructural replicase proteins to form the replication-transcription complex. Coronavirus replication-transcription complex models assume mature subunits; however, mechanisms of coronavirus maturation and replicase complex formation have yet to be defined. Here, we show that for the coronavirus murine hepatitis virus, cleavage between the nonstructural replicase proteins nsp13-14 and nsp14-15 is not required for replication but does alter RNA synthesis and recombination. These results shed new light on the requirements for coronavirus maturation and replication-transcription complex assembly, and they may reveal novel therapeutic targets and strategies for attenuation.

Mutations in the SARS-CoV-2 RNA-dependent RNA polymerase confer resistance to remdesivir by distinct mechanisms

The nucleoside analog remdesivir (RDV) is a Food and Drug Administration-approved antiviral for treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections. Thus, it is critical to understand factors that promote or prevent RDV resistance. We passaged SARS-CoV-2 in the presence of increasing concentrations of GS-441524, the parent nucleoside of RDV. After 13 passages, we isolated three viral lineages with phenotypic resistance as defined by increases in half-maximal effective concentration from 2.7- to 10.4-fold. Sequence analysis identified nonsynonymous mutations in nonstructural protein 12 RNA-dependent RNA polymerase (nsp12-RdRp): V166A, N198S, S759A, V792I, and C799F/R. Two lineages encoded the S759A substitution at the RdRp Ser759-Asp-Asp active motif. In one lineage, the V792I substitution emerged first and then combined with S759A. Introduction of S759A and V792I substitutions at homologous nsp12 positions in murine hepatitis virus demonstrated transferability across betacoronaviruses; introduction of these substitutions resulted in up to 38-fold RDV resistance and a replication defect. Biochemical analysis of SARS-CoV-2 RdRp encoding S759A demonstrated a roughly 10-fold decreased preference for RDV-triphosphate (RDV-TP) as a substrate, whereas nsp12-V792I diminished the uridine triphosphate concentration needed to overcome template-dependent inhibition associated with RDV. The in vitro-selected substitutions identified in this study were rare or not detected in the greater than 6 million publicly available nsp12-RdRp consensus sequences in the absence of RDV selection. The results define genetic and biochemical pathways to RDV resistance and emphasize the need for additional studies to define the potential for emergence of these or other RDV resistance mutations in clinical settings.

Prevention and Control of Porcine Epidemic Diarrhea: The Development of Recombination-Resistant Live Attenuated Vaccines

Porcine epidemic diarrhea (PED), causing up to 100% mortality in neonatal pigs, is a highly contagious enteric disease caused by PED virus (PEDV). The highly virulent genogroup 2 (G2) PEDV emerged in 2010 and has caused huge economic losses to the pork industry globally. It was first reported in the US in 2013, caused country-wide outbreaks, and posed tremendous hardship for many pork producers in 2013–2014. Vaccination of pregnant sows/gilts with live attenuated vaccines (LAVs) is the most effective strategy to induce lactogenic immunity in the sows/gilts and provide a passive protection via the colostrum and milk to suckling piglets against PED. However, there are still no safe and effective vaccines available after about one decade of endeavor. One of the biggest concerns is the potential reversion to virulence of an LAV in the field. In this review, we summarize the status and the major obstacles in PEDV LAV development. We also discuss the function of the transcriptional regulatory sequences in PEDV transcription, contributing to recombination, and possible strategies to prevent the reversion of LAVs. This article provides insights into the rational design of a promising LAV without safety issues.

SARS-CoV-2 Within-Host and in vitro Genomic Variability and Sub-Genomic RNA Levels Indicate Differences in Viral Expression Between Clinical Cohorts and in vitro Culture

Background

Low frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model.

Methods

Two clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33).

Results

Intrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts.

Discussion and Conclusion

The positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

How hydrolytic exoribonucleases impact human disease: Two sides of the same story

RNAs are extremely important molecules inside the cell which perform many different functions. For example, messenger RNAs, transfer RNAs, and ribosomal RNAs are involved in protein synthesis, whereas non-coding RNAs have numerous regulatory roles. Ribonucleases are the enzymes responsible for the processing and degradation of all types of RNAs, having multiple roles in every aspect of RNA metabolism. However, the involvement of RNases in disease is still not well understood. This review focuses on the involvement of the RNase II/RNB family of 3'-5' exoribonucleases in human disease. This can be attributed to direct effects, whereby mutations in the eukaryotic enzymes of this family (Dis3 (or Rrp44), Dis3L1 (or Dis3L), and Dis3L2) are associated with a disease, or indirect effects, whereby mutations in the prokaryotic counterparts of RNase II/RNB family (RNase II and/or RNase R) affect the physiology and virulence of several human pathogens. In this review, we will compare the structural and biochemical characteristics of the members of the RNase II/RNB family of enzymes. The outcomes of mutations impacting enzymatic function will be revisited, in terms of both the direct and indirect effects on disease. Furthermore, we also describe the SARS-CoV-2 viral exoribonuclease and its importance to combat COVID-19 pandemic. As a result, RNases may be a good therapeutic target to reduce bacterial and viral pathogenicity. These are the two perspectives on RNase II/RNB family enzymes that will be presented in this review.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Origin, phylogeny, variability and epitope conservation of SARS-CoV-2 worldwide

The coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses innumerous challenges, like understanding what triggered the emergence of this new human virus, how this RNA virus is evolving or how the variability of viral genome may impact the primary structure of proteins that are targets for vaccine. We analyzed 19471 SARS-CoV-2 genomes available at the GISAID database from all over the world and 3335 genomes of other Coronoviridae family members available at GenBank, collecting SARS-CoV-2 high-quality genomes and distinct Coronoviridae family genomes. Additionally, we analyzed 199,984 spike glycoprotein sequences. Here, we identify a SARS-CoV-2 emerging cluster containing 13 closely related genomes isolated from bat and pangolin that showed evidence of recombination, which may have contributed to the emergence of SARS-CoV-2. The analyzed SARS-CoV-2 genomes presented 9632 single nucleotide variants (SNVs) corresponding to a variant density of 0.3 over the genome, and a clear geographic distribution. SNVs are unevenly distributed throughout the genome and hotspots for mutations were found for the spike gene and ORF 1ab. We describe a set of predicted spike protein epitopes whose variability is negligible. Additionally, all predicted epitopes for the structural E, M and N proteins are highly conserved. The amino acid changes present in the spike glycoprotein of variables of concern (VOCs) comprise between 3.4% and 20.7% of the predicted epitopes of this protein. These results favors the continuous efficacy of the available vaccines targeting the spike protein, and other structural proteins. Multiple epitopes vaccines should sustain vaccine efficacy since at least some of the epitopes present in variability regions of VOCs are conserved and thus recognizable by antibodies.

Identifying SARS-CoV-2 antiviral compounds by screening for small molecule inhibitors of nsp14/nsp10 exoribonuclease

SARS-CoV-2 is a coronavirus that emerged in 2019 and rapidly spread across the world causing a deadly pandemic with tremendous social and economic costs. Healthcare systems worldwide are under great pressure, and there is an urgent need for effective antiviral treatments. The only currently approved antiviral treatment for COVID-19 is remdesivir, an inhibitor of viral genome replication. SARS-CoV-2 proliferation relies on the enzymatic activities of the non-structural proteins (nsp), which makes them interesting targets for the development of new antiviral treatments. With the aim to identify novel SARS-CoV-2 antivirals, we have purified the exoribonuclease/methyltransferase (nsp14) and its cofactor (nsp10) and developed biochemical assays compatible with high-throughput approaches to screen for exoribonuclease inhibitors. We have screened a library of over 5000 commercial compounds and identified patulin and aurintricarboxylic acid (ATA) as inhibitors of nsp14 exoribonuclease in vitro. We found that patulin and ATA inhibit replication of SARS-CoV-2 in a VERO E6 cell-culture model. These two new antiviral compounds will be valuable tools for further coronavirus research as well as potentially contributing to new therapeutic opportunities for COVID-19.

[Basic information of Coronavirus]

Coronaviruses are pathogens that infect many of animals, resulting in respiratory or enteric diseases. Coronaviruses constitute Nidovirales together with Arteriviridae. Most of human coronaviruses are known to cause mild illness and common cold. However, an epidemic of severe acute respiratory syndrome (SARS) occurred in 2002, ten years after SARS epidemic Middle East respiratory syndrome (MERS) emerged in 2012. Now, we face on a novel coronavirus which emerges in end of 2019. This novel coronavirus is named as SARS-CoV-2. SARS-CoV-2 is spread to worldwide within one to two months and causes coronavirus disease 2019 (COVID-19), respiratory illness. Coronaviruses are enveloped viruses possessing a positive-sense and large single stranded RNA genomes. The 5' two-thirds of the CoV genome consists of two overlapping open reading frames (ORFs 1a and 1b) that encode non-structural proteins (nsps). The other one-third of the genome consists of ORFs encoding structural proteins, including spike (S), membrane (M), envelope (E) and nucleocapsid (N) proteins, and accessory proteins. Upon infection of CoV into host cells, the translation of two precursor polyproteins, pp1a and pp1ab, occurs and these polyproteins are cleaved into 16 nsps by viral proteases. Structural proteins assemble to the vesicles located from ER to Golgi (ER Golgiintermediate compartment) and virions bud into the vesicles. Virions are released from infectedcells via exocytosis.

Crucial mutation in the exoribonuclease domain of nsp14 of PEDV leads to high genetic instability during viral replication

Background

Coronavirus (CoV) nonstructural protein 14 (nsp14) has exoribonuclease (ExoN) activity, responsible for proofreading and contributing to replication fidelity. It has been reported that CoVs exhibit variable sensitivity to nsp14-ExoN deficiency. Betacoronavirus murine hepatitis virus (MHV) and severe acute respiratory syndrome (SARS)-CoV were viable upon nsp14-ExoN deficiency. While betacoronavirus Middle East respiratory syndrome (MERS)-CoV and SARS-CoV-2 were non-viable with disabled nsp14-ExoN. In this study, we investigated the nsp14-ExoN deficiency of alphacoronavirus porcine epidemic diarrhea virus (PEDV) in viral pathogenesis using reverse genetics.

Results

Eight nsp14-ExoN deficient mutants, targeting the predicted active sites and the Zinc finger or mental-coordinating sites, of PEDV were designed. Only one mutant E191A with a mutation in the Mg²⁺-binding site was rescued using the infectious clone of PEDV PC22A strain (icPC22A). The passage no.1–3 (P1-3) of E191A grew to very low titers in Vero cells. To evaluate the pathogenesis of the E191A, 4 or 5-day-old gnotobiotic pigs were inoculated orally with 100 TCID₅₀/pig of the E191A-P1, icPC22A, or mock. All mock pigs did not shed virus in feces or show clinical signs. All pigs inoculated with icPC22A shed high viral RNA levels, had severe diarrhea, and died by 6 days post-inoculation (dpi). In contrast, only 3 pigs (3/4, 75%) in the E191A-P1 group shed low levels of viral RNA and 2 pigs had moderate diarrhea at acute infection phase. At 22 dpi, each pig was challenged orally with 10⁶ plaque forming unit of virulent icPC22A. All pigs in the mock group developed severe diarrhea and 2 of the 5 pigs died. Pigs in the E191A-P1 group had less severe diarrhea and no pigs died. Sanger sequencing analysis revealed that the viral genome in the fecal sample of one E191A-P1-inoculated pig and the P4 virus passaged in vitro lost the E191A mutation, suggesting the genetic instability of the E191A mutant.

Conclusion

The recombinant PEDV variants carrying mutations at the essential functional sites within nsp14-ExoN were either lethal or genetically unstable. Our finding further confirmed the critical role of nsp14-ExoN in CoV life cycle, suggesting that it may be a target for the design of universal anti-CoV drugs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00598-1.

New targets for drug design: importance of nsp14/nsp10 complex formation for the 3'-5' exoribonucleolytic activity on SARS-CoV-2

SARS‐CoV‐2 virus has triggered a global pandemic with devastating consequences. The understanding of fundamental aspects of this virus is of extreme importance. In this work, we studied the viral ribonuclease nsp14, one of the most interferon antagonists from SARS‐CoV‐2. Nsp14 is a multifunctional protein with two distinct activities, an N‐terminal 3’‐to‐5’ exoribonuclease (ExoN) and a C‐terminal N7‐methyltransferase (N7‐MTase), both critical for coronaviruses life cycle, indicating nsp14 as a prominent target for the development of antiviral drugs. In coronaviruses, nsp14 ExoN activity is stimulated through the interaction with the nsp10 protein. We have performed a biochemical characterization of nsp14‐nsp10 complex from SARS‐CoV‐2. We confirm the 3’‐5’ exoribonuclease and MTase activities of nsp14 and the critical role of nsp10 in upregulating the nsp14 ExoN activity. Furthermore, we demonstrate that SARS‐CoV‐2 nsp14 N7‐MTase activity is functionally independent of the ExoN activity and nsp10. A model from SARS‐CoV‐2 nsp14‐nsp10 complex allowed mapping key nsp10 residues involved in this interaction. Our results show that a stable interaction between nsp10 and nsp14 is required for the nsp14‐mediated ExoN activity of SARS‐CoV‐2. We studied the role of conserved DEDD catalytic residues of SARS‐CoV‐2 nsp14 ExoN. Our results show that motif I of ExoN domain is essential for the nsp14 function, contrasting to the functionality of these residues in other coronaviruses, which can have important implications regarding the specific pathogenesis of SARS‐CoV‐2. This work unraveled a basis for discovering inhibitors targeting specific amino acids in order to disrupt the assembly of this complex and interfere with coronaviruses replication.

Identification of multipotent drugs for COVID-19 therapeutics with the evaluation of their SARS-CoV2 inhibitory activity

The SARS-CoV2 is a highly contagious pathogen that causes COVID-19 disease. It has affected millions of people globally with an average lethality of ~3%. There is an urgent need of drugs for the treatment of COVID-19. In the current studies, we have used bioinformatics techniques to screen the FDA approved drugs against nine SARS-CoV2 proteins to identify drugs for repurposing. Additionally, we analyzed if the identified molecules can also affect the human proteins whose expression in lung changed during SARS-CoV2 infection. Targeting such genes may also be a beneficial strategy to curb disease manifestation. We have identified 74 molecules that can bind to various SARS-CoV2 and human host proteins. We experimentally validated our in-silico predictions using vero E6 cells infected with SARS-CoV2 virus. Interestingly, many of our predicted molecules viz. capreomycin, celecoxib, mefloquine, montelukast, and nebivolol showed good activity (IC50) against SARS-CoV2. We hope that these studies may help in the development of new therapeutic options for the treatment of COVID-19.

Intra-host variation and evolutionary dynamics of SARS-CoV-2 populations in COVID-19 patients

BACKGROUND

Since early February 2021, the causative agent of COVID-19, SARS-CoV-2, has infected over 104 million people with more than 2 million deaths according to official reports. The key to understanding the biology and virus-host interactions of SARS-CoV-2 requires the knowledge of mutation and evolution of this virus at both inter- and intra-host levels. However, despite quite a few polymorphic sites identified among SARS-CoV-2 populations, intra-host variant spectra and their evolutionary dynamics remain mostly unknown.

METHODS

Using high-throughput sequencing of metatranscriptomic and hybrid captured libraries, we characterized consensus genomes and intra-host single nucleotide variations (iSNVs) of serial samples collected from eight patients with COVID-19. The distribution of iSNVs along the SARS-CoV-2 genome was analyzed and co-occurring iSNVs among COVID-19 patients were identified. We also compared the evolutionary dynamics of SARS-CoV-2 population in the respiratory tract (RT) and gastrointestinal tract (GIT).

RESULTS

The 32 consensus genomes revealed the co-existence of different genotypes within the same patient. We further identified 40 intra-host single nucleotide variants (iSNVs). Most (30/40) iSNVs presented in a single patient, while ten iSNVs were found in at least two patients or identical to consensus variants. Comparing allele frequencies of the iSNVs revealed a clear genetic differentiation between intra-host populations from the respiratory tract (RT) and gastrointestinal tract (GIT), mostly driven by bottleneck events during intra-host migrations. Compared to RT populations, the GIT populations showed a better maintenance and rapid development of viral genetic diversity following the suspected intra-host bottlenecks.

CONCLUSIONS

Our findings here illustrate the intra-host bottlenecks and evolutionary dynamics of SARS-CoV-2 in different anatomic sites and may provide new insights to understand the virus-host interactions of coronaviruses and other RNA viruses.

Pattern of genomic variation in SARS-CoV-2 (COVID-19) suggests restricted nonrandom changes: Analysis using Shewhart control charts

SARS-CoV-2 is a member of the Coronavirus family which recently originated from the Wuhan province of China and spread very rapidly through the world infecting more than 4 million people. In the past, other Coronaviruses have also been found to cause human infection, but not as widespread as COVID-19. Since Coronavirus sequences constantly change due to mutation and recombination, it is important to understand the pattern of changes and likely path the virus can take in the future. In this study, we have used the Shewhart control chart to identify and analyze hypervariable (hotspots) and hypovariable (coldspots) regions of the virus. Our analysis shows that SARS-CoV-2 has changed in a few regions of the genome. Analysis of SARS-CoV-1 and MERS sequences suggests that over time, mutations start accumulating in different regions and most likely SARS-CoV-2 may also follow a similar path. The results suggest a possible emergence of modified viruses over some time.

Analysis of Indian SARS-CoV-2 Genomes Reveals Prevalence of D614G Mutation in Spike Protein Predicting an Increase in Interaction With TMPRSS2 and Virus Infectivity

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, has emerged as a global pandemic worldwide. In this study, we used ARTIC primers-based amplicon sequencing to profile 225 SARS-CoV-2 genomes from India. Phylogenetic analysis of 202 high-quality assemblies identified the presence of all the five reported clades 19A, 19B, 20A, 20B, and 20C in the population. The analyses revealed Europe and Southeast Asia as two major routes for introduction of the disease in India followed by local transmission. Interestingly, the19B clade was found to be more prevalent in our sequenced genomes (17%) compared to other genomes reported so far from India. Haplotype network analysis showed evolution of 19A and 19B clades in parallel from predominantly Gujarat state in India, suggesting it to be one of the major routes of disease transmission in India during the months of March and April, whereas 20B and 20C appeared to evolve from 20A. At the same time, 20A and 20B clades depicted prevalence of four common mutations 241 C > T in 5' UTR, P4715L, F942F along with D614G in the Spike protein. D614G mutation has been reported to increase virus shedding and infectivity. Our molecular modeling and docking analysis identified that D614G mutation resulted in enhanced affinity of Spike S1-S2 hinge region with TMPRSS2 protease, possibly the reason for increased shedding of S1 domain in G614 as compared to D614. Moreover, we also observed an increased concordance of G614 mutation with the viral load, as evident from decreased Ct value of Spike and the ORF1ab gene.

The Enzymatic Activity of the nsp14 Exoribonuclease Is Critical for Replication of MERS-CoV and SARS-CoV-2

Coronaviruses (CoVs) stand out for their large RNA genome and complex RNA-synthesizing machinery comprising 16 nonstructural proteins (nsps). The bifunctional nsp14 contains 3'-to-5' exoribonuclease (ExoN) and guanine-N7-methyltransferase (N7-MTase) domains. While the latter presumably supports mRNA capping, ExoN is thought to mediate proofreading during genome replication. In line with such a role, ExoN knockout mutants of mouse hepatitis virus (MHV) and severe acute respiratory syndrome coronavirus (SARS-CoV) were previously reported to have crippled but viable hypermutation phenotypes. Remarkably, using reverse genetics, a large set of corresponding ExoN knockout mutations has now been found to be lethal for another betacoronavirus, Middle East respiratory syndrome coronavirus (MERS-CoV). For 13 mutants, viral progeny could not be recovered, unless-as happened occasionally-reversion had first occurred. Only a single mutant was viable, likely because its E191D substitution is highly conservative. Remarkably, a SARS-CoV-2 ExoN knockout mutant was found to be unable to replicate, resembling observations previously made for alpha- and gammacoronaviruses, but starkly contrasting with the documented phenotype of ExoN knockout mutants of the closely related SARS-CoV. Subsequently, we established in vitro assays with purified recombinant MERS-CoV nsp14 to monitor its ExoN and N7-MTase activities. All ExoN knockout mutations that proved lethal in reverse genetics were found to severely decrease ExoN activity while not affecting N7-MTase activity. Our study strongly suggests that CoV nsp14 ExoN has an additional function, which apparently is critical for primary viral RNA synthesis and thus differs from the proofreading function that, based on previous MHV and SARS-CoV studies, was proposed to boost longer-term replication fidelity.IMPORTANCE The bifunctional nsp14 subunit of the coronavirus replicase contains 3'-to-5' exoribonuclease (ExoN) and guanine-N7-methyltransferase domains. For the betacoronaviruses MHV and SARS-CoV, ExoN was reported to promote the fidelity of genome replication, presumably by mediating a form of proofreading. For these viruses, ExoN knockout mutants are viable while displaying an increased mutation frequency. Strikingly, we have now established that the equivalent ExoN knockout mutants of two other betacoronaviruses, MERS-CoV and SARS-CoV-2, are nonviable, suggesting an additional and critical ExoN function in their replication. This is remarkable in light of the very limited genetic distance between SARS-CoV and SARS-CoV-2, which is highlighted, for example, by 95% amino acid sequence identity in their nsp14 sequences. For (recombinant) MERS-CoV nsp14, both its enzymatic activities were evaluated using newly developed in vitro assays that can be used to characterize these key replicative enzymes in more detail and explore their potential as target for antiviral drug development.

Imposed mutational meltdown as an antiviral strategy

Following widespread infections of the most recent coronavirus known to infect humans, SARS‐CoV‐2, attention has turned to potential therapeutic options. With no drug or vaccine yet approved, one focal point of research is to evaluate the potential value of repurposing existing antiviral treatments, with the logical strategy being to identify at least a short‐term intervention to prevent within‐patient progression, while long‐term vaccine strategies unfold. Here, we offer an evolutionary/population‐genetic perspective on one approach that may overwhelm the capacity for pathogen defense (i.e., adaptation) – induced mutational meltdown – providing an overview of key concepts, review of previous theoretical and experimental work of relevance, and guidance for future research. Applied with appropriate care, including target specificity, induced mutational meltdown may provide a general, rapidly implemented approach for the within‐patient eradication of a wide range of pathogens or other undesirable microorganisms.

The molecular virology of coronaviruses

Few human pathogens have been the focus of as much concentrated worldwide attention as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of COVID-19. Its emergence into the human population and ensuing pandemic came on the heels of severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV), two other highly pathogenic coronavirus spillovers, which collectively have reshaped our view of a virus family previously associated primarily with the common cold. It has placed intense pressure on the collective scientific community to develop therapeutics and vaccines, whose engineering relies on a detailed understanding of coronavirus biology. Here, we present the molecular virology of coronavirus infection, including its entry into cells, its remarkably sophisticated gene expression and replication mechanisms, its extensive remodeling of the intracellular environment, and its multifaceted immune evasion strategies. We highlight aspects of the viral life cycle that may be amenable to antiviral targeting as well as key features of its biology that await discovery.

Molecular Diagnosis of COVID-19: Challenges and Research Needs

Molecular diagnosis of COVID-19 primarily relies on the detection of RNA of the SARS-CoV-2 virus, the causative infectious agent of the pandemic. Reverse transcription polymerase chain reaction (RT-PCR) enables sensitive detection of specific sequences of genes that encode the RNA dependent RNA polymerase (RdRP), nucleocapsid (N), envelope (E), and spike (S) proteins of the virus. Although RT-PCR tests have been widely used and many alternative assays have been developed, the current testing capacity and availability cannot meet the unprecedented global demands for rapid, reliable, and widely accessible molecular diagnosis. Challenges remain throughout the entire analytical process, from the collection and treatment of specimens to the amplification and detection of viral RNA and the validation of clinical sensitivity and specificity. We highlight the main issues surrounding molecular diagnosis of COVID-19, including false negatives from the detection of viral RNA, temporal variations of viral loads, selection and treatment of specimens, and limiting factors in detecting viral proteins. We discuss critical research needs, such as improvements in RT-PCR, development of alternative nucleic acid amplification techniques, incorporating CRISPR technology for point-of-care (POC) applications, validation of POC tests, and sequencing of viral RNA and its mutations. Improved assays are also needed for environmental surveillance or wastewater-based epidemiology, which gauges infection on the community level through analyses of viral components in the community's wastewater. Public health surveillance benefits from large-scale analyses of antibodies in serum, although the current serological tests do not quantify neutralizing antibodies. Further advances in analytical technology and research through multidisciplinary collaboration will contribute to the development of mitigation strategies, therapeutics, and vaccines. Lessons learned from molecular diagnosis of COVID-19 are valuable for better preparedness in response to other infectious diseases.

The interface between coronaviruses and host cell RNA biology: Novel potential insights for future therapeutic intervention

Coronaviruses, including SARS-Cov-2, are RNA-based pathogens that interface with a large variety of RNA-related cellular processes during infection. These processes include capping, polyadenylation, localization, RNA stability, translation, and regulation by RNA binding proteins or noncoding RNA effectors. The goal of this article is to provide an in-depth perspective on the current state of knowledge of how various coronaviruses interact with, usurp, and/or avoid aspects of these cellular RNA biology machineries. A thorough understanding of how coronaviruses interact with RNA-related posttranscriptional processes in the cell should allow for new insights into aspects of viral pathogenesis as well as identify new potential avenues for the development of anti-coronaviral therapeutics. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Reverse genetic systems: Rational design of coronavirus live attenuated vaccines with immune sequelae

Since the end of 2019, the global COVID-19 outbreak has once again made coronaviruses a hot topic. Vaccines are hoped to be an effective way to stop the spread of the virus. However, there are no clinically approved vaccines available for coronavirus infections. Reverse genetics technology can realize the operation of RNA virus genomes at the DNA level and provide new ideas and strategies for the development of new vaccines. In this review, we systematically describe the role of reverse genetics technology in studying the effects of coronavirus proteins on viral virulence and innate immunity, cell and tissue tropism and antiviral drug screening. An efficient reverse genetics platform is useful for obtaining the ideal attenuated strain to prepare an attenuated live vaccine.

Selection of viral variants during persistent infection of insectivorous bat cells with Middle East respiratory syndrome coronavirus

Coronaviruses that cause severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) are speculated to have originated in bats. The mechanisms by which these viruses are maintained in individuals or populations of reservoir bats remain an enigma. Mathematical models have predicted long-term persistent infection with low levels of periodic shedding as a likely route for virus maintenance and spillover from bats. In this study, we tested the hypothesis that bat cells and MERS coronavirus (CoV) can co-exist in vitro. To test our hypothesis, we established a long-term coronavirus infection model of bat cells that are persistently infected with MERS-CoV. We infected cells from Eptesicus fuscus with MERS-CoV and maintained them in culture for at least 126 days. We characterized the persistently infected cells by detecting virus particles, protein and transcripts. Basal levels of type I interferon in the long-term infected bat cells were higher, relative to uninfected cells, and disrupting the interferon response in persistently infected bat cells increased virus replication. By sequencing the whole genome of MERS-CoV from persistently infected bat cells, we identified that bat cells repeatedly selected for viral variants that contained mutations in the viral open reading frame 5 (ORF5) protein. Furthermore, bat cells that were persistently infected with ΔORF5 MERS-CoV were resistant to superinfection by wildtype virus, likely due to reduced levels of the virus receptor, dipeptidyl peptidase 4 (DPP4) and higher basal levels of interferon in these cells. In summary, our study provides evidence for a model of coronavirus persistence in bats, along with the establishment of a unique persistently infected cell culture model to study MERS-CoV-bat interactions.

Genome-Wide Mutagenesis of Hepatitis C Virus Reveals Ability of Genome To Overcome Detrimental Mutations

To gain insight into the impact of mutations on the viability of the hepatitis C virus (HCV) genome, we created a set of full-genome mutant libraries, differing from the parent sequence as well as each other, by using a random mutagenesis approach; the proportion of mutations increased across these libraries with declining template amount or dATP concentration. The replication efficiencies of full-genome mutant libraries ranged between 71 and 329 focus-forming units (FFU) per 10⁵ Huh7.5 cells. Mutant libraries with low proportions of mutations demonstrated low replication capabilities, whereas those with high proportions of mutations had their replication capabilities restored. Hepatoma cells transfected with selected mutant libraries, with low (4 mutations per 10,000 bp copied), moderate (33 mutations), and high (66 mutations) proportions of mutations, and their progeny were subjected to serial passage. Predominant virus variants (mutants) from these mutant libraries (Mutant^l, Mutant^m, and Mutant^h, respectively) were evaluated for changes in growth kinetics and particle-to-FFU unit ratio, virus protein expression, and modulation of host cell protein synthesis. Mutant^m and Mutant^l variants produced >3.0-log-higher extracellular progeny per ml than the parent, and Mutant^h produced progeny at a rate 1.0-log lower. More than 80% of the mutations were in a nonstructural part of the mutant genomes, the majority were nonsynonymous, and a moderate to large proportion were in the conserved regions. Our results suggest that the HCV genome has the ability to overcome lethal/deleterious mutations because of the high reproduction rate but highly selects for random, beneficial mutations.IMPORTANCE Hepatitis C virus (HCV) in vivo displays high genetic heterogeneity, which is partly due to the high reproduction and random substitutions during error-prone genome replication. It is difficult to introduce random substitutions in vitro because of limitations in inducing mutagenesis from the 5' end to the 3' end of the genome. Our study has overcome this limitation. We synthesized full-length genomes with few to several random mutations in the background of an HCV clone that can recapitulate all steps of the life cycle. Our study provides evidence of the capability of the HCV genome to overcome deleterious mutations and remain viable. Mutants that emerged from the libraries had diverse phenotype profiles compared to the parent, and putative adaptive mutations mapped to segments of the conserved nonstructural genome. We demonstrate the potential utility of our system for the study of sequence variation that ensures the survival and adaptation of HCV.

Small-Molecule Antiviral β-d-N4-Hydroxycytidine Inhibits a Proofreading-Intact Coronavirus with a High Genetic Barrier to Resistance

Coronaviruses (CoVs) have emerged from animal reservoirs to cause severe and lethal disease in humans, but there are currently no FDA-approved antivirals to treat the infections. One class of antiviral compounds, nucleoside analogues, mimics naturally occurring nucleosides to inhibit viral replication. While these compounds have been successful therapeutics for several viral infections, mutagenic nucleoside analogues, such as ribavirin and 5-fluorouracil, have been ineffective at inhibiting CoVs. This has been attributed to the proofreading activity of the viral 3'-5' exoribonuclease (ExoN). β-d-N4-Hydroxycytidine (NHC) (EIDD-1931; Emory Institute for Drug Development) has recently been reported to inhibit multiple viruses. Here, we demonstrate that NHC inhibits both murine hepatitis virus (MHV) (50% effective concentration [EC50] = 0.17 μM) and Middle East respiratory syndrome CoV (MERS-CoV) (EC50 = 0.56 μM) with minimal cytotoxicity. NHC inhibited MHV lacking ExoN proofreading activity similarly to wild-type (WT) MHV, suggesting an ability to evade or overcome ExoN activity. NHC inhibited MHV only when added early during infection, decreased viral specific infectivity, and increased the number and proportion of G:A and C:U transition mutations present after a single infection. Low-level NHC resistance was difficult to achieve and was associated with multiple transition mutations across the genome in both MHV and MERS-CoV. These results point to a virus-mutagenic mechanism of NHC inhibition in CoVs and indicate a high genetic barrier to NHC resistance. Together, the data support further development of NHC for treatment of CoVs and suggest a novel mechanism of NHC interaction with the CoV replication complex that may shed light on critical aspects of replication.IMPORTANCE The emergence of coronaviruses (CoVs) into human populations from animal reservoirs has demonstrated their epidemic capability, pandemic potential, and ability to cause severe disease. However, no antivirals have been approved to treat these infections. Here, we demonstrate the potent antiviral activity of a broad-spectrum ribonucleoside analogue, β-d-N4-hydroxycytidine (NHC), against two divergent CoVs. Viral proofreading activity does not markedly impact sensitivity to NHC inhibition, suggesting a novel interaction between a nucleoside analogue inhibitor and the CoV replicase. Further, passage in the presence of NHC generates only low-level resistance, likely due to the accumulation of multiple potentially deleterious transition mutations. Together, these data support a mutagenic mechanism of inhibition by NHC and further support the development of NHC for treatment of CoV infections.

Fitness Barriers Limit Reversion of a Proofreading-Deficient Coronavirus

Coronaviruses encode an exoribonuclease (ExoN) that is important for viral replication, fitness, and virulence, yet coronaviruses with a defective ExoN (ExoN-AA) have not reverted under diverse experimental conditions. In this study, we identify multiple impediments to MHV-ExoN-AA reversion. We show that ExoN-AA reversion is possible but evolutionarily unfavorable. Instead, compensatory mutations outside ExoN-AA motif I are more accessible and beneficial than partial reversion. We also show that coevolution between replicase proteins over long-term passage partially compensates for ExoN-AA motif I but renders the virus inhospitable to a reverted ExoN. Our results reveal the evolutionary basis for the genetic stability of ExoN-inactivating mutations, illuminate complex functional and evolutionary relationships between coronavirus replicase proteins, and identify potential mechanisms for stabilization of ExoN-AA coronavirus mutants.

The 3′-to-5′ exoribonuclease in coronavirus (CoV) nonstructural protein 14 (nsp14-ExoN) mediates RNA proofreading during genome replication. ExoN catalytic residues are arranged in three motifs: I (DE), II (E), and III (D). Alanine replacement of the motif I residues (AA-E-D; four nucleotide substitutions) in murine hepatitis virus (MHV) and severe acute respiratory syndrome (SARS)-CoV yields viable mutants with impaired replication and fitness, increased mutation rates, and attenuated virulence in vivo. Despite these impairments, MHV- and SARS-CoV ExoN motif I AA mutants (ExoN-AA) have not reverted at motif I in diverse in vitro and in vivo environments, suggesting that profound fitness barriers prevent motif I reversion. To test this hypothesis, we engineered MHV-ExoN-AA with 1, 2, or 3 nucleotide mutations along genetic pathways to AA-to-DE reversion. We show that engineered intermediate revertants were viable but had no increased replication or competitive fitness compared to that of MHV-ExoN-AA. In contrast, a low-passage-number (passage 10 [P10]) MHV-ExoN-AA showed increased replication and competitive fitness without reversion of ExoN-AA. Finally, engineered reversion of ExoN-AA to ExoN-DE in the presence of ExoN-AA passage-adaptive mutations resulted in significant fitness loss. These results demonstrate that while reversion is possible, at least one alternative adaptive pathway is more rapidly advantageous than intermediate revertants and may alter the genetic background to render reversion detrimental to fitness. Our results provide an evolutionary rationale for lack of ExoN-AA reversion, illuminate potential multiprotein replicase interactions and coevolution, and support future studies aimed at stabilizing attenuated CoV ExoN-AA mutants.

IMPORTANCE Coronaviruses encode an exoribonuclease (ExoN) that is important for viral replication, fitness, and virulence, yet coronaviruses with a defective ExoN (ExoN-AA) have not reverted under diverse experimental conditions. In this study, we identify multiple impediments to MHV-ExoN-AA reversion. We show that ExoN-AA reversion is possible but evolutionarily unfavorable. Instead, compensatory mutations outside ExoN-AA motif I are more accessible and beneficial than partial reversion. We also show that coevolution between replicase proteins over long-term passage partially compensates for ExoN-AA motif I but renders the virus inhospitable to a reverted ExoN. Our results reveal the evolutionary basis for the genetic stability of ExoN-inactivating mutations, illuminate complex functional and evolutionary relationships between coronavirus replicase proteins, and identify potential mechanisms for stabilization of ExoN-AA coronavirus mutants.

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

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

Chikungunya Virus Fidelity Variants Exhibit Differential Attenuation and Population Diversity in Cell Culture and Adult Mice

Chikungunya virus (CHIKV) is a reemerging global health threat that produces debilitating arthritis in people. Like other RNA viruses with high mutation rates, CHIKV produces populations of genetically diverse genomes within a host. While several known CHIKV mutations influence disease severity in vertebrates and transmission by mosquitoes, the role of intrahost diversity in chikungunya arthritic disease has not been studied. In this study, high- and low-fidelity CHIKV variants, previously characterized by altered in vitro population mutation frequencies, were used to evaluate how intrahost diversity influences clinical disease, CHIKV replication, and antibody neutralization in immunocompetent adult mice inoculated in the rear footpads. Both high- and low-fidelity mutations were hypothesized to attenuate CHIKV arthritic disease, replication, and neutralizing antibody levels compared to wild-type (WT) CHIKV. Unexpectedly, high-fidelity mutants elicited more severe arthritic disease than the WT despite comparable CHIKV replication, whereas a low-fidelity mutant produced attenuated disease and replication. Serum antibody developed against both high- and low-fidelity CHIKV exhibited reduced neutralization of WT CHIKV. Using next-generation sequencing (NGS), the high-fidelity mutations were demonstrated to be genetically stable but produced more genetically diverse populations than WT CHIKV in mice. This enhanced diversification was subsequently reproduced after serial in vitro passage. The NGS results contrast with previously reported population diversities for fidelity variants, which focused mainly on part of the E1 gene, and highlight the need for direct measurements of mutation rates to clarify CHIKV fidelity phenotypes.IMPORTANCE CHIKV is a reemerging global health threat that elicits debilitating arthritis in humans. There are currently no commercially available CHIKV vaccines. Like other RNA viruses, CHIKV has a high mutation rate and is capable of rapid intrahost diversification during an infection. In other RNA viruses, virus population diversity associates with disease progression; however, potential impacts of intrahost viral diversity on CHIKV arthritic disease have not been studied. Using previously characterized CHIKV fidelity variants, we addressed whether CHIKV population diversity influences the severity of arthritis and host antibody response in an arthritic mouse model. Our findings show that CHIKV populations with greater genetic diversity can cause more severe disease and stimulate antibody responses with reduced neutralization of low-diversity virus populations in vitro The discordant high-fidelity phenotypes in this study highlight the complexity of inferring replication fidelity indirectly from population diversity.

RNA Virus Fidelity Mutants: A Useful Tool for Evolutionary Biology or a Complex Challenge?

RNA viruses replicate with low fidelity due to the error-prone nature of the RNA-dependent RNA polymerase, which generates approximately one mutation per round of genome replication. Due to the large population sizes produced by RNA viruses during replication, this results in a cloud of closely related virus variants during host infection, of which small increases or decreases in replication fidelity have been shown to result in virus attenuation in vivo, but not typically in vitro. Since the discovery of the first RNA virus fidelity mutants during the mid-aughts, the field has exploded with the identification of over 50 virus fidelity mutants distributed amongst 7 RNA virus families. This review summarizes the current RNA virus fidelity mutant literature, with a focus upon the definition of a fidelity mutant as well as methods to confirm any mutational changes associated with the fidelity mutant. Due to the complexity of such a definition, in addition to reports of unstable virus fidelity phenotypes, the future translational utility of these mutants and applications for basic science are examined.

Low-fidelity Venezuelan equine encephalitis virus polymerase mutants to improve live-attenuated vaccine safety and efficacy

During RNA virus replication, there is the potential to incorporate mutations that affect virulence or pathogenesis. For live-attenuated vaccines, this has implications for stability, as replication may result in mutations that either restore the wild-type phenotype via reversion or compensate for the attenuating mutations by increasing virulence (pseudoreversion). Recent studies have demonstrated that altering the mutation rate of an RNA virus is an effective attenuation tool. To validate the safety of low-fidelity mutations to increase vaccine attenuation, several mutations in the RNA-dependent RNA-polymerase (RdRp) were tested in the live-attenuated Venezuelan equine encephalitis virus vaccine strain, TC-83. Next generation sequencing after passage in the presence of mutagens revealed a mutant containing three mutations in the RdRp, TC-83 3x, to have decreased replication fidelity, while a second mutant, TC-83 4x displayed no change in fidelity, but shared many phenotypic characteristics with TC-83 3x. Both mutants exhibited increased, albeit inconsistent attenuation in an infant mouse model, as well as increased immunogenicity and complete protection against lethal challenge of an adult murine model compared with the parent TC-83. During serial passaging in a highly permissive model, the mutants increased in virulence but remained less virulent than the parent TC-83. These results suggest that the incorporation of low-fidelity mutations into the RdRp of live-attenuated vaccines for RNA viruses can confer increased immunogenicity whilst showing some evidence of increased attenuation. However, while in theory such constructs may result in more effective vaccines, the instability of the vaccine phenotype decreases the likelihood of this being an effective vaccine strategy.

Murine Hepatitis Virus nsp14 Exoribonuclease Activity Is Required for Resistance to Innate Immunity

Coronaviruses (CoVs) are positive-sense RNA viruses that infect numerous mammalian and avian species and are capable of causing severe and lethal disease in humans. CoVs encode several innate immune antagonists that counteract the host innate immune response to facilitate efficient viral replication. CoV nonstructural protein 14 (nsp14) encodes 3'-to-5' exoribonuclease activity (ExoN), which performs a proofreading function and is required for high-fidelity replication. Outside of the order Nidovirales, arenaviruses are the only RNA viruses that encode an ExoN, which functions to degrade double-stranded RNA (dsRNA) replication intermediates. In this study, we tested the hypothesis that CoV ExoN also functions to antagonize the innate immune response. We demonstrate that viruses lacking ExoN activity [ExoN(-)] are sensitive to cellular pretreatment with interferon beta (IFN-β) in a dose-dependent manner. In addition, ExoN(-) virus replication was attenuated in wild-type bone marrow-derived macrophages (BMMs) and partially restored in interferon alpha/beta receptor-deficient (IFNAR^-/-) BMMs. ExoN(-) virus replication did not result in IFN-β gene expression, and in the presence of an IFN-β-mediated antiviral state, ExoN(-) viral RNA levels were not substantially reduced relative to those of untreated samples. However, ExoN(-) virus generated from IFN-β-pretreated cells had reduced specific infectivity and decreased relative fitness, suggesting that ExoN(-) virus generated during an antiviral state is less viable to establish a subsequent infection. Overall, our data suggest murine hepatitis virus (MHV) ExoN activity is required for resistance to the innate immune response, and antiviral mechanisms affecting the viral RNA sequence and/or an RNA modification act on viruses lacking ExoN activity.IMPORTANCE CoVs encode multiple antagonists that prevent or disrupt an efficient innate immune response. Additionally, no specific antiviral therapies or vaccines currently exist for human CoV infections. Therefore, the study of CoV innate immune antagonists is essential for understanding how CoVs overcome host defenses and to maximize potential therapeutic interventions. Here, we sought to determine the contributions of nsp14 ExoN activity in the induction of and resistance to the innate immune response. We show that viruses lacking nsp14 ExoN activity are more sensitive than wild-type MHV to restriction by exogenous IFN-β and that viruses produced in the presence of an antiviral state are less capable of establishing a subsequent viral infection. Our results support the hypothesis that murine hepatitis virus ExoN activity is required for resistance to the innate immune response.