ViReMaShiny: An Interactive Application for Analysis of Viral Recombination Data

Recombination is an essential driver of virus evolution and adaption, giving rise to new chimeric viruses, structural variants, sub-genomic RNAs, and Defective-RNAs. Next-Generation Sequencing of virus samples, either from experimental or clinical settings, has revealed a complex distribution of recombination events that contributes to the intrahost diversity. We and others have previously developed alignment tools to discover and map these diverse recombination events in NGS data. However, there is no standard for data visualization to contextualize events of interest and downstream analysis often requires bespoke coding. To address this, we present ViReMaShiny, a web-based application built using the R Shiny framework to allow interactive exploration and point-and-click visualization of viral recombination data provided in BED format generated by computational pipelines such as ViReMa (Viral-Recombination-Mapper).

The coronavirus proofreading exoribonuclease mediates extensive viral recombination

Coronaviruses (CoVs) emerge as zoonoses and cause severe disease in humans, demonstrated by the SARS-CoV-2 (COVID-19) pandemic. RNA recombination is required during normal CoV replication for subgenomic mRNA (sgmRNA) synthesis and generates defective viral genomes (DVGs) of unknown function. However, the determinants and patterns of CoV recombination are unknown. Here, we show that divergent {beta}-CoVs SARS-CoV-2, MERS-CoV, and murine hepatitis virus (MHV) perform extensive RNA recombination in culture, generating similar patterns of recombination junctions and diverse populations of DVGs and sgmRNAs. We demonstrate that the CoV proofreading nonstructural protein (nsp14) 3-to-5 exoribonuclease (nsp14-ExoN) is required for normal CoV recombination and that its genetic inactivation causes significantly decreased frequency and altered patterns of recombination in both infected cells and released virions. Thus, nsp14-ExoN is a key determinant of both high fidelity CoV replication and recombination, and thereby represents a highly-conserved and vulnerable target for virus inhibition and attenuation.

QTQTN motif upstream of the furin-cleavage site plays key role in SARS-CoV-2 infection and pathogenesis.

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing 1,2. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates 3. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif ({Delta}QTQTN). Here we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS, and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated4, and disruption its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site - the FCS, loop length, and glycosylation - are required for efficient SARS-CoV-2 replication and pathogenesis.

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis.

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203-205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral RG motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2s continued adaptation to human infection. Author SummarySince its emergence, SARS-CoV-2 has continued to adapt for human infection resulting in the emergence of variants with unique genetic profiles. Most studies of genetic variation have focused on spike, the target of currently available vaccines, leaving the importance of variation elsewhere understudied. Here, we characterize a highly variable motif at residues 203-205 in nucleocapsid. Recreating the prominent nucleocapsid R203K+G204R mutation in an early pandemic background, we show that this mutation is alone sufficient to enhance SARS-CoV-2 replication and pathogenesis. We also link augmentation of SARS-CoV-2 infection by the R203K+G204R mutation to its modulation of nucleocapsid phosphorylation. Finally, we characterize an analogous alanine double substitution at positions 203-204. This mutant was found to mimic R203K+G204R, suggesting augmentation of infection occurs by disrupting the ancestral sequence. Together, our findings illustrate that mutations outside of spike are key components of SARS-CoV-2s adaptation to human infection.

SARS-CoV-2 Uses Nonstructural Protein 16 to Evade Restriction by IFIT1 and IFIT3

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2-O methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2-O MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo, using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive to type I interferon (IFN-I) in vitro. Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2-O methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, a methyltransferase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a possible target for future antiviral therapies. ImportanceSimilar to other coronaviruses, disruption of SARS-CoV-2 NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo, our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1, but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2-O methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called Tiled-ClickSeq, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Furin Cleavage Site Is Key to SARS-CoV-2 Pathogenesis

SARS-CoV-2 has resulted in a global pandemic and shutdown economies around the world. Sequence analysis indicates that the novel coronavirus (CoV) has an insertion of a furin cleavage site (PRRAR) in its spike protein. Absent in other group 2B CoVs, the insertion may be a key factor in the replication and virulence of SARS-CoV-2. To explore this question, we generated a SARS-CoV-2 mutant lacking the furin cleavage site ({Delta}PRRA) in the spike protein. This mutant virus replicated with faster kinetics and improved fitness in Vero E6 cells. The mutant virus also had reduced spike protein processing as compared to wild-type SARS-CoV-2. In contrast, the {Delta}PRRA had reduced replication in Calu3 cells, a human respiratory cell line, and had attenuated disease in a hamster pathogenesis model. Despite the reduced disease, the {Delta}PRRA mutant offered robust protection from SARS-CoV-2 rechallenge. Importantly, plaque reduction neutralization tests (PRNT50) with COVID-19 patient sera and monoclonal antibodies against the receptor-binding domain found a shift, with the mutant virus resulting in consistently reduced PRNT50 titers. Together, these results demonstrate a critical role for the furin cleavage site insertion in SARS-CoV-2 replication and pathogenesis. In addition, these findings illustrate the importance of this insertion in evaluating neutralization and other downstream SARS-CoV-2 assays. ImportanceAs COVID-19 has impacted the world, understanding how SARS-CoV-2 replicates and causes virulence offers potential pathways to disrupt its disease. By removing the furin cleavage site, we demonstrate the importance of this insertion to SARS-CoV-2 replication and pathogenesis. In addition, the findings with Vero cells indicate the likelihood of cell culture adaptations in virus stocks that can influence reagent generation and interpretation of a wide range of data including neutralization and drug efficacy. Overall, our work highlights the importance of this key motif in SARS-CoV-2 infection and pathogenesis. Article SummaryA deletion of the furin cleavage site in SARS-CoV-2 amplifies replication in Vero cells, but attenuates replication in respiratory cells and pathogenesis in vivo. Loss of the furin site also reduces susceptibility to neutralization in vitro.