A putative new SARS-CoV protein, 3c, encoded in an ORF overlapping ORF3a

The Main Protease of Middle East Respiratory Syndrome Coronavirus Induces Cleavage of Mitochondrial Antiviral Signaling Protein to Antagonize the Innate Immune Response

Mitochondrial antiviral signaling protein (MAVS) is a crucial signaling adaptor in the sensing of positive-sense RNA viruses and the subsequent induction of the innate immune response. Coronaviruses have evolved multiple mechanisms to evade this response, amongst others, through their main protease (Mpro), which is responsible for the proteolytic cleavage of the largest part of the viral replicase polyproteins pp1a and pp1ab. Additionally, it can cleave cellular substrates, such as innate immune signaling factors, to dampen the immune response. Here, we show that MAVS is cleaved in cells infected with Middle East respiratory syndrome coronavirus (MERS-CoV), but not in cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This cleavage was independent of cellular negative feedback mechanisms that regulate MAVS activation. Furthermore, MERS-CoV Mpro expression induced MAVS cleavage upon overexpression and suppressed the activation of the interferon-β (IFN-β) and nuclear factor-κB (NF-κB) response. We conclude that we have uncovered a novel mechanism by which MERS-CoV downregulates the innate immune response, which is not observed among other highly pathogenic coronaviruses.

The HLA-II immunopeptidome of SARS-CoV-2

Targeted synthetic vaccines have the potential to transform our response to viral outbreaks, yet the design of these vaccines requires a comprehensive knowledge of viral immunogens. Here, we report severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) peptides that are naturally processed and loaded onto human leukocyte antigen-II (HLA-II) complexes in infected cells. We identify over 500 unique viral peptides from canonical proteins as well as from overlapping internal open reading frames. Most HLA-II peptides colocalize with known CD4+ T cell epitopes in coronavirus disease 2019 patients, including 2 reported immunodominant regions in the SARS-CoV-2 membrane protein. Overall, our analyses show that HLA-I and HLA-II pathways target distinct viral proteins, with the structural proteins accounting for most of the HLA-II peptidome and nonstructural and noncanonical proteins accounting for the majority of the HLA-I peptidome. These findings highlight the need for a vaccine design that incorporates multiple viral elements harboring CD4+ and CD8+ T cell epitopes to maximize vaccine effectiveness.

ORF3c is expressed in SARS-CoV-2-infected cells and inhibits innate sensing by targeting MAVS

Most SARS-CoV-2 proteins are translated from subgenomic RNAs (sgRNAs). While the majority of these sgRNAs are monocistronic, some viral mRNAs encode more than one protein. One example is the ORF3a sgRNA that also encodes ORF3c, an enigmatic 41-amino-acid peptide. Here, we show that ORF3c is expressed in SARS-CoV-2-infected cells and suppresses RIG-I- and MDA5-mediated IFN-β induction. ORF3c interacts with the signaling adaptor MAVS, induces its C-terminal cleavage, and inhibits the interaction of RIG-I with MAVS. The immunosuppressive activity of ORF3c is conserved among members of the subgenus sarbecovirus, including SARS-CoV and coronaviruses isolated from bats. Notably, however, the SARS-CoV-2 delta and kappa variants harbor premature stop codons in ORF3c, demonstrating that this reading frame is not essential for efficient viral replication in vivo and is likely compensated by other viral proteins. In agreement with this, disruption of ORF3c does not significantly affect SARS-CoV-2 replication in CaCo-2, CaLu-3, or Rhinolophus alcyone cells. In summary, we here identify ORF3c as an immune evasion factor of SARS-CoV-2 that suppresses innate sensing in infected cells.

The SARS-CoV-2 protein ORF3c is a mitochondrial modulator of innate immunity

The SARS-CoV-2 genome encodes a multitude of accessory proteins. Using comparative genomic approaches, an additional accessory protein, ORF3c, has been predicted to be encoded within the ORF3a sgmRNA. Expression of ORF3c during infection has been confirmed independently by ribosome profiling. Despite ORF3c also being present in the 2002-2003 SARS-CoV, its function has remained unexplored. Here we show that ORF3c localizes to mitochondria, where it inhibits innate immunity by restricting IFN-β production, but not NF-κB activation or JAK-STAT signaling downstream of type I IFN stimulation. We find that ORF3c is inhibitory after stimulation with cytoplasmic RNA helicases RIG-I or MDA5 or adaptor protein MAVS, but not after TRIF, TBK1 or phospho-IRF3 stimulation. ORF3c co-immunoprecipitates with the antiviral proteins MAVS and PGAM5 and induces MAVS cleavage by caspase-3. Together, these data provide insight into an uncharacterized mechanism of innate immune evasion by this important human pathogen.

Alternative Reading Frames are an Underappreciated Source of Protein Sequence Novelty

Protein-coding DNA sequences can be translated into completely different amino acid sequences if the nucleotide triplets used are shifted by a non-triplet amount on the same DNA strand or by translating codons from the opposite strand. Such "alternative reading frames" of protein-coding genes are a major contributor to the evolution of novel protein products. Recent studies demonstrating this include examples across the three domains of cellular life and in viruses. These sequences increase the number of trials potentially available for the evolutionary invention of new genes and also have unusual properties which may facilitate gene origin. There is evidence that the structure of the standard genetic code contributes to the features and gene-likeness of some alternative frame sequences. These findings have important implications across diverse areas of molecular biology, including for genome annotation, structural biology, and evolutionary genomics.

A novel approach to finding conserved features in low-variability gene alignments characterises RNA motifs in SARS-CoV and SARS-CoV-2

Collections of genetic sequences belonging to related organisms contain information on the evolutionary constraints to which the organisms have been subjected. Heavily constrained regions can be investigated to understand their roles in an organism's life cycle, and drugs can be sought to disrupt these roles. In organisms with low genetic diversity, such as newly-emerged pathogens, it is key to obtain this information early to develop new treatments. Here, we present methods that ensure we can leverage all the information available in a low-signal, low-noise set of sequences, to find contiguous regions of relatively conserved nucleic acid. We demonstrate the application of these methods by analysing over 5 million genome sequences of the recently-emerged RNA virus SARS-CoV-2 and correlating these results with an analysis of 119 genome sequences of SARS-CoV. We propose the precise location of a previously described packaging signal, and discuss explanations for other regions of high conservation.

SARS-CoV-2 ORF3c impairs mitochondrial respiratory metabolism, oxidative stress, and autophagic flux

Coronaviruses encode a variable number of accessory proteins that are involved in host-virus interaction, suppression of immune responses, or immune evasion. SARS-CoV-2 encodes at least twelve accessory proteins, whose roles during infection have been studied. Nevertheless, the role of the ORF3c accessory protein, an alternative open reading frame of ORF3a, has remained elusive. Herein, we show that the ORF3c protein has a mitochondrial localization and alters mitochondrial metabolism, inducing a shift from glucose to fatty acids oxidation and enhanced oxidative phosphorylation. These effects result in increased ROS production and block of the autophagic flux. In particular, ORF3c affects lysosomal acidification, blocking the normal autophagic degradation process and leading to autolysosome accumulation. We also observed different effect on autophagy for SARS-CoV-2 and batCoV RaTG13 ORF3c proteins; the 36R and 40K sites are necessary and sufficient to determine these effects.

The HLA-II immunopeptidome of SARS-CoV-2

Targeted synthetic vaccines have the potential to transform our response to viral outbreaks; yet the design of these vaccines requires a comprehensive knowledge of viral immunogens, including T-cell epitopes. Having previously mapped the SARS-CoV-2 HLA-I landscape, here we report viral peptides that are naturally processed and loaded onto HLA-II complexes in infected cells. We identified over 500 unique viral peptides from canonical proteins, as well as from overlapping internal open reading frames (ORFs), revealing, for the first time, the contribution of internal ORFs to the HLA-II peptide repertoire. Most HLA-II peptides co-localized with the known CD4+ T cell epitopes in COVID-19 patients. We also observed that two reported immunodominant regions in the SARS-CoV-2 membrane protein are formed at the level of HLA-II presentation. Overall, our analyses show that HLA-I and HLA-II pathways target distinct viral proteins, with the structural proteins accounting for most of the HLA-II peptidome and non-structural and non-canonical proteins accounting for the majority of the HLA-I peptidome. These findings highlight the need for a vaccine design that incorporates multiple viral elements harboring CD4+ and CD8+ T cell epitopes to maximize the vaccine effectiveness.

Immunosuppression as a Hub for SARS-CoV-2 Mutational Drift

The clinical course of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), is largely determined by host factors, with a wide range of outcomes. Despite an extensive vaccination campaign and high rates of infection worldwide, the pandemic persists, adapting to overcome antiviral immunity acquired through prior exposure. The source of many such major adaptations is variants of concern (VOCs), novel SARS-CoV-2 variants produced by extraordinary evolutionary leaps whose origins remain mostly unknown. In this study, we tested the influence of factors on the evolutionary course of SARS-CoV-2. Electronic health records of individuals infected with SARS-CoV-2 were paired to viral whole-genome sequences to assess the effects of host clinical parameters and immunity on the intra-host evolution of SARS-CoV-2. We found slight, albeit significant, differences in SARS-CoV-2 intra-host diversity, which depended on host parameters such as vaccination status and smoking. Only one viral genome had significant alterations as a result of host parameters; it was found in an immunocompromised, chronically infected woman in her 70s. We highlight the unusual viral genome obtained from this woman, which had an accelerated mutational rate and an excess of rare mutations, including near-complete truncating of the accessory protein ORF3a. Our findings suggest that the evolutionary capacity of SARS-CoV-2 during acute infection is limited and mostly unaffected by host characteristics. Significant viral evolution is seemingly exclusive to a small subset of COVID-19 cases, which typically prolong infections in immunocompromised patients. In these rare cases, SARS-CoV-2 genomes accumulate many impactful and potentially adaptive mutations; however, the transmissibility of such viruses remains unclear.

Computational Analysis Predicts Correlations among Amino Acids in SARS-CoV-2 Proteomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of -0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

Discovery of Novel Chinese Medicine Compounds Targeting 3CL Protease by Virtual Screening and Molecular Dynamics Simulation

The transmission and infectivity of COVID-19 have caused a pandemic that has lasted for several years. This is due to the constantly changing variants and subvariants that have evolved rapidly from SARS-CoV-2. To discover drugs with therapeutic potential for COVID-19, we focused on the 3CL protease (3CL^pro) of SARS-CoV-2, which has been proven to be an important target for COVID-19 infection. Computational prediction techniques are quick and accurate enough to facilitate the discovery of drugs against the 3CL^pro of SARS-CoV-2. In this paper, we used both ligand-based virtual screening and structure-based virtual screening to screen the traditional Chinese medicine small molecules that have the potential to target the 3CL^pro of SARS-CoV-2. MD simulations were used to confirm these results for future in vitro testing. MCCS was then used to calculate the normalized free energy of each ligand and the residue energy contribution. As a result, we found ZINC15676170, ZINC09033700, and ZINC12530139 to be the most promising antiviral therapies against the 3CL^pro of SARS-CoV-2.

Translation landscape of SARS-CoV-2 noncanonical subgenomic RNAs

The ongoing COVID-19 pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with a positive-stranded RNA genome. Current proteomic studies of SARS-CoV-2 mainly focus on the proteins encoded by its genomic RNA (gRNA) or canonical subgenomic RNAs (sgRNAs). Here, we systematically investigated the translation landscape of SARS-CoV-2, especially its noncanonical sgRNAs. We first constructed a strict pipeline, named vipep, for identifying reliable peptides derived from RNA viruses using RNA-seq and mass spectrometry data. We applied vipep to analyze 24 sets of mass spectrometry data related to SARS-CoV-2 infection. In addition to known canonical proteins, we identified many noncanonical sgRNA-derived peptides, which stably increase after viral infection. Furthermore, we explored the potential functions of those proteins encoded by noncanonical sgRNAs and found that they can bind to viral RNAs and may have immunogenic activity. The generalized vipep pipeline is applicable to any RNA viruses and these results have expanded the SARS-CoV-2 translation map, providing new insights for understanding the functions of SARS-CoV-2 sgRNAs.

D155Y substitution of SARS-CoV-2 ORF3a weakens binding with Caveolin-1

The clinical manifestation of the recent pandemic COVID-19, caused by the novel SARS-CoV-2 virus, varies from mild to severe respiratory illness. Although environmental, demographic and co-morbidity factors have an impact on the severity of the disease, contribution of the mutations in each of the viral genes towards the degree of severity needs a deeper understanding for designing a better therapeutic approach against COVID-19. Open Reading Frame-3a (ORF3a) protein has been found to be mutated at several positions. In this work, we have studied the effect of one of the most frequently occurring mutants, D155Y of ORF3a protein, found in Indian COVID-19 patients. Using computational simulations we demonstrated that the substitution at 155th changed the amino acids involved in salt bridge formation, hydrogen-bond occupancy, interactome clusters, and the stability of the protein compared with the other substitutions found in Indian patients. Protein-protein docking using HADDOCK analysis revealed that substitution D155Y weakened the binding affinity of ORF3a with caveolin-1 compared with the other substitutions, suggesting its importance in the overall stability of ORF3a-caveolin-1 complex, which may modulate the virulence property of SARS-CoV-2.

Structures and functions of coronavirus replication-transcription complexes and their relevance for SARS-CoV-2 drug design

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed millions of people and continues to cause massive global upheaval. Coronaviruses are positive-strand RNA viruses with an unusually large genome of ~30 kb. They express an RNA-dependent RNA polymerase and a cohort of other replication enzymes and supporting factors to transcribe and replicate their genomes. The proteins performing these essential processes are prime antiviral drug targets, but drug discovery is hindered by our incomplete understanding of coronavirus RNA synthesis and processing. In infected cells, the RNA-dependent RNA polymerase must coordinate with other viral and host factors to produce both viral mRNAs and new genomes. Recent research aiming to decipher and contextualize the structures, functions and interplay of the subunits of the SARS-CoV-2 replication and transcription complex proteins has burgeoned. In this Review, we discuss recent advancements in our understanding of the molecular basis and complexity of the coronavirus RNA-synthesizing machinery. Specifically, we outline the mechanisms and regulation of RNA translation, replication and transcription. We also discuss the composition of the replication and transcription complexes and their suitability as targets for antiviral therapy.

Computational methods for inferring location and genealogy of overlapping genes in virus genomes: approaches and applications

Viruses may evolve to increase the amount of encoded genetic information by means of overlapping genes, which utilize several reading frames. Such overlapping genes may be especially impactful for genomes of small size, often serving a source of novel accessory proteins, some of which play a crucial role in viral pathogenicity or in promoting the systemic spread of virus. Diverse genome-based metrics were proposed to facilitate recognition of overlapping genes that otherwise may be overlooked during genome annotation. They can detect the atypical codon bias associated with the overlap (e.g. a statistically significant reduction in variability at synonymous sites) or other sequence-composition features peculiar to overlapping genes. In this review, I compare nine computational methods, discuss their strengths and limitations, and survey how they were applied to detect candidate overlapping genes in the genome of SARS-CoV-2, the etiological agent of COVID-19 pandemic.

Science's Response to CoVID-19

CoVID-19 is a multi-symptomatic disease which has made a global impact due to its ability to spread rapidly, and its relatively high mortality rate. Beyond the heroic efforts to develop vaccines, which we do not discuss herein, the response of scientists and clinicians to this complex problem has reflected the need to detect CoVID-19 rapidly, to diagnose patients likely to show adverse symptoms, and to treat severe and critical CoVID-19. Here we aim to encapsulate these varied and sometimes conflicting approaches and the resulting data in terms of chemistry and biology. In the process we highlight emerging concepts, and potential future applications that may arise out of this immense effort.

Prediction of two novel overlapping ORFs in the genome of SARS-CoV-2

Six candidate overlapping genes have been detected in SARS-CoV-2, yet current methods struggle to detect overlapping genes that recently originated. However, such genes might encode proteins beneficial to the virus, and provide a model system to understand gene birth. To complement existing detection methods, I first demonstrated that selection pressure to avoid stop codons in alternative reading frames is a driving force in the origin and retention of overlapping genes. I then built a detection method, CodScr, based on this selection pressure. Finally, I combined CodScr with methods that detect other properties of overlapping genes, such as a biased nucleotide and amino acid composition. I detected two novel ORFs (ORF-Sh and ORF-Mh), overlapping the spike and membrane genes respectively, which are under selection pressure and may be beneficial to SARS-CoV-2. ORF-Sh and ORF-Mh are present, as ORF uninterrupted by stop codons, in 100% and 95% of the SARS-CoV-2 genomes, respectively.

SARS-CoV-2: from its discovery to genome structure, transcription, and replication

SARS-CoV-2 is an extremely contagious respiratory virus causing adult atypical pneumonia COVID-19 with severe acute respiratory syndrome (SARS). SARS-CoV-2 has a single-stranded, positive-sense RNA (+RNA) genome of ~ 29.9 kb and exhibits significant genetic shift from different isolates. After entering the susceptible cells expressing both ACE2 and TMPRSS2, the SARS-CoV-2 genome directly functions as an mRNA to translate two polyproteins from the ORF1a and ORF1b region, which are cleaved by two viral proteases into sixteen non-structural proteins (nsp1-16) to initiate viral genome replication and transcription. The SARS-CoV-2 genome also encodes four structural (S, E, M and N) and up to six accessory (3a, 6, 7a, 7b, 8, and 9b) proteins, but their translation requires newly synthesized individual subgenomic RNAs (sgRNA) in the infected cells. Synthesis of the full-length viral genomic RNA (gRNA) and sgRNAs are conducted inside double-membrane vesicles (DMVs) by the viral replication and transcription complex (RTC), which comprises nsp7, nsp8, nsp9, nsp12, nsp13 and a short RNA primer. To produce sgRNAs, RTC starts RNA synthesis from the highly structured gRNA 3' end and switches template at various transcription regulatory sequence (TRSB) sites along the gRNA body probably mediated by a long-distance RNA-RNA interaction. The TRS motif in the gRNA 5' leader (TRSL) is responsible for the RNA-RNA interaction with the TRSB upstream of each ORF and skipping of the viral genome in between them to produce individual sgRNAs. Abundance of individual sgRNAs and viral gRNA synthesized in the infected cells depend on the location and read-through efficiency of each TRSB. Although more studies are needed, the unprecedented COVID-19 pandemic has taught the world a painful lesson that is to invest and proactively prepare future emergence of other types of coronaviruses and any other possible biological horrors.

Coronavirus, the King Who Wanted More Than a Crown: From Common to the Highly Pathogenic SARS-CoV-2, Is the Key in the Accessory Genes?

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), that emerged in late 2019, is the etiologic agent of the current "coronavirus disease 2019" (COVID-19) pandemic, which has serious health implications and a significant global economic impact. Of the seven human coronaviruses, all of which have a zoonotic origin, the pandemic SARS-CoV-2, is the third emerging coronavirus, in the 21st century, highly pathogenic to the human population. Previous human coronavirus outbreaks (SARS-CoV-1 and MERS-CoV) have already provided several valuable information on some of the common molecular and cellular mechanisms of coronavirus infections as well as their origin. However, to meet the new challenge caused by the SARS-CoV-2, a detailed understanding of the biological specificities, as well as knowledge of the origin are crucial to provide information on viral pathogenicity, transmission and epidemiology, and to enable strategies for therapeutic interventions and drug discovery. Therefore, in this review, we summarize the current advances in SARS-CoV-2 knowledges, in light of pre-existing information of other recently emerging coronaviruses. We depict the specificity of the immune response of wild bats and discuss current knowledge of the genetic diversity of bat-hosted coronaviruses that promotes viral genome expansion (accessory gene acquisition). In addition, we describe the basic virology of coronaviruses with a special focus SARS-CoV-2. Finally, we highlight, in detail, the current knowledge of genes and accessory proteins which we postulate to be the major keys to promote virus adaptation to specific hosts (bat and human), to contribute to the suppression of immune responses, as well as to pathogenicity.

SARS-CoV-2 Accessory Proteins in Viral Pathogenesis: Knowns and Unknowns

There are still many unanswered questions concerning viral SARS-CoV-2 pathogenesis in COVID-19. Accessory proteins in SARS-CoV-2 consist of eleven viral proteins whose roles during infection are still not completely understood. Here, a review on the current knowledge of SARS-CoV-2 accessory proteins is summarized updating new research that could be critical in understanding SARS-CoV-2 interaction with the host. Some accessory proteins such as ORF3b, ORF6, ORF7a and ORF8 have been shown to be important IFN-I antagonists inducing an impairment in the host immune response. In addition, ORF3a is involved in apoptosis whereas others like ORF9b and ORF9c interact with cellular organelles leading to suppression of the antiviral response in infected cells. However, possible roles of ORF7b and ORF10 are still awaiting to be described. Also, ORF3d has been reassigned. Relevant information on the knowns and the unknowns in these proteins is analyzed, which could be crucial for further understanding of SARS-CoV-2 pathogenesis and to design strategies counteracting their actions evading immune responses in COVID-19.

Replication Kinetics of B.1.351 and B.1.1.7 SARS-CoV-2 Variants of Concern Including Assessment of a B.1.1.7 Mutant Carrying a Defective ORF7a Gene

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of COVID-19, is a readily transmissible and potentially deadly pathogen which is currently re-defining human susceptibility to pandemic viruses in the modern world. The recent emergence of several genetically distinct descendants known as variants of concern (VOCs) is further challenging public health disease management, due to increased rates of virus transmission and potential constraints on vaccine effectiveness. We report the isolation of SARS-CoV-2 VOCs imported into Australia belonging to the B.1.351 lineage, first described in the Republic of South Africa (RSA), and the B.1.1.7 lineage originally reported in the United Kingdom, and directly compare the replication kinetics of these two VOCs in Vero E6 cells. In this analysis, we also investigated a B.1.1.7 VOC (QLD1516/2021) carrying a 7-nucleotide deletion in the open reading frame 7a (ORF7a) gene, likely truncating and rendering the ORF7a protein of this virus defective. We demonstrate that the replication of the B.1.351 VOC (QLD1520/2020) in Vero E6 cells can be detected earlier than the B.1.1.7 VOCs (QLD1516/2021 and QLD1517/2021), before peaking at 48 h post infection (p.i.), with significantly higher levels of virus progeny. Whilst replication of the ORF7a defective isolate QLD1516/2021 was delayed longer than the other viruses, slightly more viral progeny was produced by the mutant compared to the unmutated isolate QLD1517/2021 at 72 h p.i. Collectively, these findings contribute to our understanding of SARS-CoV-2 replication and evolutionary dynamics, which have important implications in the development of future vaccination, antiviral therapies, and epidemiological control strategies for COVID-19.

SARS-CoV-2 gene content and COVID-19 mutation impact by comparing 44 Sarbecovirus genomes

Despite its clinical importance, the SARS-CoV-2 gene set remains unresolved, hindering dissection of COVID-19 biology. We use comparative genomics to provide a high-confidence protein-coding gene set, characterize evolutionary constraint, and prioritize functional mutations. We select 44 Sarbecovirus genomes at ideally-suited evolutionary distances, and quantify protein-coding evolutionary signatures and overlapping constraint. We find strong protein-coding signatures for ORFs 3a, 6, 7a, 7b, 8, 9b, and a novel alternate-frame gene, ORF3c, whereas ORFs 2b, 3d/3d-2, 3b, 9c, and 10 lack protein-coding signatures or convincing experimental evidence of protein-coding function. Furthermore, we show no other conserved protein-coding genes remain to be discovered. Mutation analysis suggests ORF8 contributes to within-individual fitness but not person-to-person transmission. Cross-strain and within-strain evolutionary pressures agree, except for fewer-than-expected within-strain mutations in nsp3 and S1, and more-than-expected in nucleocapsid, which shows a cluster of mutations in a predicted B-cell epitope, suggesting immune-avoidance selection. Evolutionary histories of residues disrupted by spike-protein substitutions D614G, N501Y, E484K, and K417N/T provide clues about their biology, and we catalog likely-functional co-inherited mutations. Previously reported RNA-modification sites show no enrichment for conservation. Here we report a high-confidence gene set and evolutionary-history annotations providing valuable resources and insights on SARS-CoV-2 biology, mutations, and evolution.

SARS-CoV-2 variants reveal features critical for replication in primary human cells

Since entering the human population, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2; the causative agent of Coronavirus Disease 2019 [COVID-19]) has spread worldwide, causing >100 million infections and >2 million deaths. While large-scale sequencing efforts have identified numerous genetic variants in SARS-CoV-2 during its circulation, it remains largely unclear whether many of these changes impact adaptation, replication, or transmission of the virus. Here, we characterized 14 different low-passage replication-competent human SARS-CoV-2 isolates representing all major European clades observed during the first pandemic wave in early 2020. By integrating viral sequencing data from patient material, virus stocks, and passaging experiments, together with kinetic virus replication data from nonhuman Vero-CCL81 cells and primary differentiated human bronchial epithelial cells (BEpCs), we observed several SARS-CoV-2 features that associate with distinct phenotypes. Notably, naturally occurring variants in Orf3a (Q57H) and nsp2 (T85I) were associated with poor replication in Vero-CCL81 cells but not in BEpCs, while SARS-CoV-2 isolates expressing the Spike D614G variant generally exhibited enhanced replication abilities in BEpCs. Strikingly, low-passage Vero-derived stock preparation of 3 SARS-CoV-2 isolates selected for substitutions at positions 5/6 of E and were highly attenuated in BEpCs, revealing a key cell-specific function to this region. Rare isolate-specific deletions were also observed in the Spike furin cleavage site during Vero-CCL81 passage, but these were rapidly selected against in BEpCs, underscoring the importance of this site for SARS-CoV-2 replication in primary human cells. Overall, our study uncovers sequence features in SARS-CoV-2 variants that determine cell-specific replication and highlights the need to monitor SARS-CoV-2 stocks carefully when phenotyping newly emerging variants or potential variants of concern.

Ipomoeassin-F inhibits the in vitro biogenesis of the SARS-CoV-2 spike protein and its host cell membrane receptor

In order to produce proteins essential for their propagation, many pathogenic human viruses, including SARS-CoV-2, the causative agent of COVID-19 respiratory disease, commandeer host biosynthetic machineries and mechanisms. Three major structural proteins, the spike, envelope and membrane proteins, are amongst several SARS-CoV-2 components synthesised at the endoplasmic reticulum (ER) of infected human cells prior to the assembly of new viral particles. Hence, the inhibition of membrane protein synthesis at the ER is an attractive strategy for reducing the pathogenicity of SARS-CoV-2 and other obligate viral pathogens. Using an in vitro system, we demonstrate that the small molecule inhibitor ipomoeassin F (Ipom-F) potently blocks the Sec61-mediated ER membrane translocation and/or insertion of three therapeutic protein targets for SARS-CoV-2 infection; the viral spike and ORF8 proteins together with angiotensin-converting enzyme 2, the host cell plasma membrane receptor. Our findings highlight the potential for using ER protein translocation inhibitors such as Ipom-F as host-targeting, broad-spectrum antiviral agents.This article has an associated First Person interview with the first author of the paper.

Ipomoeassin-F inhibits the in vitro biogenesis of the SARS-CoV-2 spike protein and its host cell membrane receptor

In order to produce proteins essential for their propagation, many pathogenic human viruses, including SARS-CoV-2 the causative agent of COVID-19 respiratory disease, commandeer host biosynthetic machineries and mechanisms. Three major structural proteins, the spike, envelope and membrane proteins, are amongst several SARS-CoV-2 components synthesised at the endoplasmic reticulum (ER) of infected human cells prior to the assembly of new viral particles. Hence, the inhibition of membrane protein synthesis at the ER is an attractive strategy for reducing the pathogenicity of SARS-CoV-2 and other obligate viral pathogens. Using an in vitro system, we demonstrate that the small molecule inhibitor ipomoeassin F (Ipom-F) potently blocks the Sec61-mediated ER membrane translocation/insertion of three therapeutic protein targets for SARS-CoV-2 infection; the viral spike and ORF8 proteins together with angiotensin-converting enzyme 2, the host cell plasma membrane receptor. Our findings highlight the potential for using ER protein translocation inhibitors such as Ipom-F as host-targeting, broad-spectrum, antiviral agents.

The Short- and Long-Range RNA-RNA Interactome of SARS-CoV-2

The Coronaviridae is a family of positive-strand RNA viruses that includes SARS-CoV-2, the etiologic agent of the COVID-19 pandemic. Bearing the largest single-stranded RNA genomes in nature, coronaviruses are critically dependent on long-distance RNA-RNA interactions to regulate the viral transcription and replication pathways. Here we experimentally mapped the in vivo RNA-RNA interactome of the full-length SARS-CoV-2 genome and subgenomic mRNAs. We uncovered a network of RNA-RNA interactions spanning tens of thousands of nucleotides. These interactions reveal that the viral genome and subgenomes adopt alternative topologies inside cells and engage in different interactions with host RNAs. Notably, we discovered a long-range RNA-RNA interaction, the FSE-arch, that encircles the programmed ribosomal frameshifting element. The FSE-arch is conserved in the related MERS-CoV and is under purifying selection. Our findings illuminate RNA structure-based mechanisms governing replication, discontinuous transcription, and translation of coronaviruses and will aid future efforts to develop antiviral strategies.

SARS-CoV-2 gene content and COVID-19 mutation impact by comparing 44 Sarbecovirus genomes

Despite its overwhelming clinical importance, the SARS-CoV-2 gene set remains unresolved, hindering dissection of COVID-19 biology. Here, we use comparative genomics to provide a high-confidence protein-coding gene set, characterize protein-level and nucleotide-level evolutionary constraint, and prioritize functional mutations from the ongoing COVID-19 pandemic. We select 44 complete Sarbecovirus genomes at evolutionary distances ideally-suited for protein-coding and non-coding element identification, create whole-genome alignments, and quantify protein-coding evolutionary signatures and overlapping constraint. We find strong protein-coding signatures for all named genes and for 3a, 6, 7a, 7b, 8, 9b, and also ORF3c, a novel alternate-frame gene. By contrast, ORF10, and overlapping-ORFs 9c, 3b, and 3d lack protein-coding signatures or convincing experimental evidence and are not protein-coding. Furthermore, we show no other protein-coding genes remain to be discovered. Cross-strain and within-strain evolutionary pressures largely agree at the gene, amino-acid, and nucleotide levels, with some notable exceptions, including fewer-than-expected mutations in nsp3 and Spike subunit S1, and more-than-expected mutations in Nucleocapsid. The latter also shows a cluster of amino-acid-changing variants in otherwise-conserved residues in a predicted B-cell epitope, which may indicate positive selection for immune avoidance. Several Spike-protein mutations, including D614G, which has been associated with increased transmission, disrupt otherwise-perfectly-conserved amino acids, and could be novel adaptations to human hosts. The resulting high-confidence gene set and evolutionary-history annotations provide valuable resources and insights on COVID-19 biology, mutations, and evolution.

Dynamically evolving novel overlapping gene as a factor in the SARS-CoV-2 pandemic

Understanding the emergence of novel viruses requires an accurate and comprehensive annotation of their genomes. Overlapping genes (OLGs) are common in viruses and have been associated with pandemics but are still widely overlooked. We identify and characterize ORF3d, a novel OLG in SARS-CoV-2 that is also present in Guangxi pangolin-CoVs but not other closely related pangolin-CoVs or bat-CoVs. We then document evidence of ORF3d translation, characterize its protein sequence, and conduct an evolutionary analysis at three levels: between taxa (21 members of Severe acute respiratory syndrome-related coronavirus), between human hosts (3978 SARS-CoV-2 consensus sequences), and within human hosts (401 deeply sequenced SARS-CoV-2 samples). ORF3d has been independently identified and shown to elicit a strong antibody response in COVID-19 patients. However, it has been misclassified as the unrelated gene ORF3b, leading to confusion. Our results liken ORF3d to other accessory genes in emerging viruses and highlight the importance of OLGs.

The coding capacity of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic¹. To understand the pathogenicity and antigenic potential of SARS-CoV-2 and to develop therapeutic tools, it is essential to profile the full repertoire of its expressed proteins. The current map of SARS-CoV-2 coding capacity is based on computational predictions and relies on homology with other coronaviruses. As the protein complement varies among coronaviruses, especially in regard to the variety of accessory proteins, it is crucial to characterize the specific range of SARS-CoV-2 proteins in an unbiased and open-ended manner. Here, using a suite of ribosome-profiling techniques^2-4, we present a high-resolution map of coding regions in the SARS-CoV-2 genome, which enables us to accurately quantify the expression of canonical viral open reading frames (ORFs) and to identify 23 unannotated viral ORFs. These ORFs include upstream ORFs that are likely to have a regulatory role, several in-frame internal ORFs within existing ORFs, resulting in N-terminally truncated products, as well as internal out-of-frame ORFs, which generate novel polypeptides. We further show that viral mRNAs are not translated more efficiently than host mRNAs; instead, virus translation dominates host translation because of the high levels of viral transcripts. Our work provides a resource that will form the basis of future functional studies.

SARS-CoV-2 gene content and COVID-19 mutation impact by comparing 44 Sarbecovirus genomes

Despite its overwhelming clinical importance, the SARS-CoV-2 gene set remains unresolved, hindering dissection of COVID-19 biology. Here, we use comparative genomics to provide a high-confidence protein-coding gene set, characterize protein-level and nucleotide-level evolutionary constraint, and prioritize functional mutations from the ongoing COVID-19 pandemic. We select 44 complete Sarbecovirus genomes at evolutionary distances ideally-suited for protein-coding and non-coding element identification, create whole-genome alignments, and quantify protein-coding evolutionary signatures and overlapping constraint. We find strong protein-coding signatures for all named genes and for 3a, 6, 7a, 7b, 8, 9b, and also ORF3c, a novel alternate-frame gene. By contrast, ORF10, and overlapping-ORFs 9c, 3b, and 3d lack protein-coding signatures or convincing experimental evidence and are not protein-coding. Furthermore, we show no other protein-coding genes remain to be discovered. Cross-strain and within-strain evolutionary pressures largely agree at the gene, amino-acid, and nucleotide levels, with some notable exceptions, including fewer-than-expected mutations in nsp3 and Spike subunit S1, and more-than-expected mutations in Nucleocapsid. The latter also shows a cluster of amino-acid-changing variants in otherwise-conserved residues in a predicted B-cell epitope, which may indicate positive selection for immune avoidance. Several Spike-protein mutations, including D614G, which has been associated with increased transmission, disrupt otherwise-perfectly-conserved amino acids, and could be novel adaptations to human hosts. The resulting high-confidence gene set and evolutionary-history annotations provide valuable resources and insights on COVID-19 biology, mutations, and evolution.