Search

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.

Background Genetic recombination is a tremendous source of intrahost diversity in viruses and is critical for their ability to rapidly adapt to new environments or fitness challenges. While viruses are routinely characterized using high-throughput sequencing techniques, characterizing the genetic products of recombination in next-generation sequencing data remains a challenge. Viral recombination events can be highly diverse and variable in nature, including simple duplications and deletions, or more complex events such as copy/snap-back recombination, intervirus or intersegment recombination, and insertions of host nucleic acids. Due to the variable mechanisms driving virus recombination and the different selection pressures acting on the progeny, recombination junctions rarely adhere to simple canonical sites or sequences. Furthermore, numerous different events may be present simultaneously in a viral population, yielding a complex mutational landscape. Findings We have previously developed an algorithm called ViReMa (Virus Recombination Mapper) that bootstraps the bowtie short-read aligner to capture and annotate a wide range of recombinant species found within virus populations. Here, we have updated ViReMa to provide an “error density” function designed to accurately detect recombination events in the longer reads now routinely generated by the Illumina platforms and provide output reports for multiple types of recombinant species using standardized formats. We demonstrate the utility and flexibility of ViReMa in different settings to report deletion events in simulated data from Flock House virus, copy-back RNA species in Sendai viruses, short duplication events in HIV, and virus-to-host recombination in an archaeal DNA virus.

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. 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. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (Δ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 glycosylated, and disruption of 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.

SARS-CoV-2 is the etiological agent responsible for the COVID-19 pandemic. It is estimated that only 10 aerosol-borne virus particles are sufficient to establish a secondary infection with SARS-CoV-2. However, the dispersal pattern of SARS-CoV-2 is highly variable and only 10– 20% of cases are responsible for up 80% of secondary infections. The heterogeneous nature of SARS-CoV-2 transmission suggests that super-spreader events play an important role in viral transmission. Super-spreader events occur when a single person is responsible for an unusually high number of secondary infections due to a combination of biological, environmental, and/or behavioral factors. While super-spreader events have been identified as a significant factor driving SARS-CoV-2 transmission, epidemiologic studies have consistently shown that education settings do not play a major role in community transmission. However, an outbreak of SARS-CoV-2 was recently reported among 186 children (aged 10-17) and adults (aged 18 +) after attending an overnight summer camp in Texas in June 2021. To understand the transmission dynamics of the outbreak, RNA was isolated from 36 nasopharyngeal swabs collected from patients that attended the camp and 19 control patients with no known connection to the outbreak. Genome sequencing on the Oxford Nanopore platform was performed using the ARTIC approaches for library preparation and bioinformatic analysis. SARS-CoV-2 amplicons were produced from all RNA samples and >70% of the viral genome was successfully reconstructed with >10X coverage for 46 samples. Phylogenetic methods were used to estimate the transmission history and suggested that the outbreak was the result of a single introduction. We also found evidence for secondary transmission from campers to the community. Together, these findings demonstrate that super-spreader events may occur during large gatherings of children.

Genetic recombination is a tremendous source of intra-host diversity in viruses and is critical for their ability to rapidly adapt to new environments or fitness challenges. While viruses are routinely characterized using high-throughput sequencing techniques, characterizing the genetic products of recombination in next-generation sequencing data remains a challenge. Viral recombination events can be highly diverse and variable in nature, including simple duplications and deletions, or more complex events such as copy/snap-back recombination, inter-virus or inter-segment recombination and insertions of host nucleic acids. Due to the variable mechanisms driving virus recombination and the different selection pressures acting on the progeny, recombination junctions rarely adhere to simple canonical sites or sequences. Furthermore, numerous different events may be present simultaneously in a viral population, yielding a complex mutational landscape. We have previously developed an algorithm called ViReMa (Virus Recombination Mapper) that bootstraps the bowtie short-read aligner to capture and annotate a wide-range of recombinant species found within virus populations. Here, we have updated ViReMa to provide an ‘error-density’ function designed to accurately detect recombination events in the longer reads now routinely generated by the Illumina platforms and provide output reports for multiple types of recombinant species using standardized formats. We demonstrate the utility and flexibility of ViReMa in different settings to report deletion events in simulated data from Flock House virus, copy-back RNA species in Sendai viruses, short duplication events in HIV, and virus to host recombination in an archaeal DNA virus.

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 (Δ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.

Adaptation of viruses to their environments occurs through the acquisition of both novel Single-Nucleotide Variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.

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 5’UTR, 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.