Conflicting and ambiguous names of overlapping ORFs in the SARS-CoV-2 genome: A homology-based resolution

Transcriptomic profiles of monocyte-derived macrophages exposed to SARS-CoV-2 VOCs reveal immune-evasion escape driven by delta

BACKGROUND

Since the breakout of COVID-19, the mutated forms of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have shown enhanced rates of transmission and adaptation to humans. The variants of concern (VOC), designated Alpha, Beta, Gamma, Delta, and Omicron emerged independent of one another, and in turn rapidly became dominant. The success of each VOC, as well as the virus fitness, were enabled by altered intrinsic functional properties and, reasonably, to virus antigenicity changes, conferring the ability to evade a primed immune response.

METHODS

We analysed the gene expression profiles of monocyte-derived macrophages (MDM) isolated from whole blood of healthy participants exposed to the 5 different SARS-CoV-2 VOC: D614G, Alpha (B.1.1.7), Gamma (P1), Delta (B.1.617.2), and Omicron BA.1 (B.1.1.529), and to the HCoV-OC43 strain, a coronavirus already present in the population before the SARS-CoV-2 pandemic. Whole transcriptome RNA-Seq, for both coding and non-coding RNAs, was then made.

RESULTS

After exposure to the 5 VOC of MDM, we initially assessed the presence of the viral SARS-CoV-2 transcripts to confirm viral entry. We then analysed the RNA-Seq data and observed a significant deregulation of both coding and non-coding RNAs. In particular, our RNA-Seq analysis showed a significant up-regulation of several genes involved in different immunological processes, such as PARP9/PARP14 axes, in macrophages exposed to D614G, Alpha, and Gamma variants. Surprisingly, our data showed that macrophages exposed to the Delta variant exhibited a transcriptional profile more similar to the naïve control group, while macrophages exposed to the Omicron variant showed intermediate differentially expressed genes (DEGs) between the two groups. By checking the canonical markers for M1/M2 differentiation states, we did not observe any expression in macrophages exposed to the Delta variant, suggesting an M0 status, comparable to the naïve control group. Finally, we observed a significant deregulation of 3 main types of non-coding RNAs (ncRNAs): long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and small nucleolar RNAs (snoRNAs), some of which are common to coronaviruses, and some specific to SARS-CoV-2.

CONCLUSION

The SARS-CoV-2-dependent alteration of the transcriptome of monocyte-derived macrophage (MDM)-infected cells can be linked to the chronological order of the variants' appearance in the human population. Our data suggest an evolution of VOC in modulating the host immune response, with a strong change in pace beginning with the advent of the Delta variant. MDMs exposed to Delta showed a failure in the activation of the adaptive immune response, and this correlates with the more severe symptoms developed by people affected with this SARS-CoV-2 variant.

Examining the feasibility of replacing ORF3a with fluorescent genes to construct SARS-CoV-2 reporter viruses

The SARS-CoV-2 genome encodes at least nine accessory proteins, including innate immune antagonist and putative viroporin ORF3a. ORF3a plays a role in many stages of the viral replication cycle, including immune modulation. We constructed two recombinant (r)SARS-CoV-2 viruses in which the ORF3a gene was replaced with mScarlet (mS) or mNeonGreen (mNG), denoted as rSARS-CoV-2-Δ3a-mS and rSARS-CoV-2-Δ3a-mNG, respectively. rSARS-CoV-2-Δ3a-mNG generated a fluorescent signal after infection in both A549-ACE-2-TMPRSS2 (AAT) and Vero-E6-TMPRSS2 (VTN) cells, unlike rSARS-CoV-2-Δ3a-mS. rSARS-CoV-2-Δ3a-mS mS protein could be detected immunologically in VTN but not AAT cells, indicating the expression of a non-fluorescent mS protein. The analysis of the viral transcriptomes in infected AAT cells by nanopore direct RNA sequencing (dRNAseq) revealed that the level of mS transcript was below the limit of detection in AAT cells. rSARS-CoV-2-Δ3a-mNG virus was found to be genetically stable in AAT and VTN cells, but rSARS-CoV-2-Δ3a-mS acquired partial deletions of the mS gene during sequential passaging in VTN cells, creating the virus rSARS-CoV-2-Δ3a-ΔmS. The mS deletion in VTN cells removes the chromophore coding sequence, and this may explain the presence of a non-fluorescent mS protein detected in VTN cells. The rSARS-CoV-2-Δ3a-mNG, rSARS-CoV-2-Δ3a-mS and rSARS-CoV-2-Δ3a-ΔmS viruses all replicated to a lower titre and produced smaller plaques than the parental rSARS-CoV-2-S-D614G. Interestingly, the rSARS-CoV-2-Δ3a-ΔmS virus produced higher virus titres and larger plaque sizes than rSARS-CoV-2-Δ3a-mS. This suggested that both the insertion of mS coding sequence and the deletion of ORF3a coding sequence contributed to attenuation. In comparison with rSARS-CoV-2, the rSARS-CoV-2-Δ3a-mS and rSARS-CoV-2-Δ3a-mNG viruses showed increased sensitivity to pre-treatment of cells with IFN-α but did not exhibit a dose-dependent increase in replication in the presence of the Janus kinase-signal transducer and activator of transcription signalling pathway inhibitor, ruxolitinib. In conclusion, the replacement of the ORF3a coding sequence with those of fluorescent reporter proteins attenuated the replication of SARS-CoV-2 and its ability to effectively evade the innate immune response in vitro.

Is Autophagy a Friend or Foe in SARS-CoV-2 Infection?

As obligate parasites, viruses need to hijack resources from infected cells to complete their lifecycle. The interaction between the virus and host determines the viral infection process, including viral propagation and the disease's outcome. Understanding the interaction between the virus and host factors is a basis for unraveling the intricate biological processes in the infected cells and thereby developing more efficient and targeted antivirals. Among the various fundamental virus-host interactions, autophagy plays vital and also complicated roles by directly engaging in the viral lifecycle and functioning as an anti- and/or pro-viral factor. Autophagy thus becomes a promising target against virus infection. Since the COVID-19 pandemic, there has been an accumulation of studies aiming to investigate the roles of autophagy in SARS-CoV-2 infection by using different models and from distinct angles, providing valuable information for systematically and comprehensively dissecting the interplay between autophagy and SARS-CoV-2. In this review, we summarize the advancements in the studies of the interaction between SARS-CoV-2 and autophagy, as well as detailed molecular mechanisms. We also update the current knowledge on the pharmacological strategies used to suppress SARS-CoV-2 replication through remodeling autophagy. These extensive studies on SARS-CoV-2 and autophagy can advance our understanding of virus-autophagy interaction and provide insights into developing efficient antiviral therapeutics by regulating autophagy.

Computational investigation of the inhibitory interaction of IRF3 and SARS-CoV-2 accessory protein ORF3b

ORF3b is one of the SARS-CoV-2 accessory proteins. Previous experimental study suggested that ORF3b prevents IRF3 translocating to nucleus. However, the biophysical mechanism of ORF3b-IRF3 interaction is elusive. Here, we explored the conformation ensemble of ORF3b using all-atom replica exchange molecular dynamics simulation. Disordered ORF3b has mixed α-helix, β-turn and loop conformers. The potential ORF3b-IRF3 binding modes were searched by docking representative ORF3b conformers with IRF3, and 50 ORF3b-IRF3 complex poses were screened using molecular dynamics simulations ranging from 500 to 1000 ns. We found that ORF3b binds IRF3 predominantly on its CBP binding and phosphorylated pLxIS motifs, with CBP binding site has the highest binding affinity. The ORF3b-IRF3 binding residues are highly conserved in SARS-CoV-2. Our results provided biophysics insights into ORF3b-IRF3 interaction and explained its interferon antagonism mechanism.

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

SARS-CoV-2 biology and host interactions

The zoonotic emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the ensuing coronavirus disease 2019 (COVID-19) pandemic have profoundly affected our society. The rapid spread and continuous evolution of new SARS-CoV-2 variants continue to threaten global public health. Recent scientific advances have dissected many of the molecular and cellular mechanisms involved in coronavirus infections, and large-scale screens have uncovered novel host-cell factors that are vitally important for the virus life cycle. In this Review, we provide an updated summary of the SARS-CoV-2 life cycle, gene function and virus-host interactions, including recent landmark findings on general aspects of coronavirus biology and newly discovered host factors necessary for virus replication.

The most exposed regions of SARS-CoV-2 structural proteins are subject to strong positive selection and gene overlap may locally modify this behavior

The SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) pandemic that emerged in 2019 has been an unprecedented event in international science, as it has been possible to sequence millions of genomes, tracking their evolution very closely. This has enabled various types of secondary analyses of these genomes, including the measurement of their sequence selection pressure. In this work, we have been able to measure the selective pressure of all the described SARS-CoV-2 genes, even analyzed by sequence regions, and we show how this type of analysis allows us to separate the genes between those subject to positive selection (usually those that code for surface proteins or those exposed to the host immune system) and those subject to negative selection because they require greater conservation of their structure and function. We have also seen that when another gene with an overlapping reading frame appears within a gene sequence, the overlapping sequence between the two genes evolves under a stronger purifying selection than the average of the non-overlapping regions of the main gene. We propose this type of analysis as a useful tool for locating and analyzing all the genes of a viral genome when an adequate number of sequences are available.IMPORTANCEWe have analyzed the selection pressure of all severe acute respiratory syndrome coronavirus 2 genes by means of the nonsynonymous (Ka) to synonymous (Ks) substitution rate. We found that protein-coding genes are exposed to strong positive selection, especially in the regions of interaction with other molecules (host receptor and genome of the virus itself). However, overlapping coding regions are more protected and show negative selection. This suggests that this measure could be used to study viral gene function as well as overlapping genes.

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.

Fitness effects of mutations to SARS-CoV-2 proteins

Knowledge of the fitness effects of mutations to SARS-CoV-2 can inform assessment of new variants, design of therapeutics resistant to escape, and understanding of the functions of viral proteins. However, experimentally measuring effects of mutations is challenging: we lack tractable lab assays for many SARS-CoV-2 proteins, and comprehensive deep mutational scanning has been applied to only two SARS-CoV-2 proteins. Here, we develop an approach that leverages millions of publicly available SARS-CoV-2 sequences to estimate effects of mutations. We first calculate how many independent occurrences of each mutation are expected to be observed along the SARS-CoV-2 phylogeny in the absence of selection. We then compare these expected observations to the actual observations to estimate the effect of each mutation. These estimates correlate well with deep mutational scanning measurements. For most genes, synonymous mutations are nearly neutral, stop-codon mutations are deleterious, and amino acid mutations have a range of effects. However, some viral accessory proteins are under little to no selection. We provide interactive visualizations of effects of mutations to all SARS-CoV-2 proteins (https://jbloomlab.github.io/SARS2-mut-fitness/). The framework we describe is applicable to any virus for which the number of available sequences is sufficiently large that many independent occurrences of each neutral mutation are observed.

Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.

SARS-CoV-2 accessory proteins involvement in inflammatory and profibrotic processes through IL11 signaling

SARS-CoV-2, the cause of the COVID-19 pandemic, possesses eleven accessory proteins encoded in its genome. Their roles during infection are still not completely understood. In this study, transcriptomics analysis revealed that both WNT5A and IL11 were significantly up-regulated in A549 cells expressing individual accessory proteins ORF6, ORF8, ORF9b or ORF9c from SARS-CoV-2 (Wuhan-Hu-1 isolate). IL11 is a member of the IL6 family of cytokines. IL11 signaling-related genes were also differentially expressed. Bioinformatics analysis disclosed that both WNT5A and IL11 were involved in pulmonary fibrosis idiopathic disease and functional assays confirmed their association with profibrotic cell responses. Subsequently, data comparison with lung cell lines infected with SARS-CoV-2 or lung biopsies from patients with COVID-19, evidenced altered profibrotic gene expression that matched those obtained in this study. Our results show ORF6, ORF8, ORF9b and ORF9c involvement in inflammatory and profibrotic responses. Thus, these accessory proteins could be targeted by new therapies against COVID-19 disease.

Interferon antagonists encoded by SARS-CoV-2 at a glance

The innate immune system is a powerful barrier against invading pathogens. Interferons (IFNs) are a major part of the cytokine-mediated anti-viral innate immune response. After recognition of a pathogen by immune sensors, signaling cascades are activated that culminate in the release of IFNs. These activate cells in an autocrine or paracrine fashion eventually setting cells in an anti-viral state via upregulation of hundreds of interferon-stimulated genes (ISGs). To evade the anti-viral effect of the IFN system, successful viruses like the pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolved strategies to counteract both IFN induction and signaling. In fact, more than half of the about 30 proteins encoded by SARS-CoV-2 target the IFN system at multiple levels to escape IFN-mediated restriction. Here, we review recent insights into the molecular mechanisms used by SARS-CoV-2 proteins to suppress IFN production and the establishment of an anti-viral state.

GENCODE: reference annotation for the human and mouse genomes in 2023

GENCODE produces high quality gene and transcript annotation for the human and mouse genomes. All GENCODE annotation is supported by experimental data and serves as a reference for genome biology and clinical genomics. The GENCODE consortium generates targeted experimental data, develops bioinformatic tools and carries out analyses that, along with externally produced data and methods, support the identification and annotation of transcript structures and the determination of their function. Here, we present an update on the annotation of human and mouse genes, including developments in the tools, data, analyses and major collaborations which underpin this progress. For example, we report the creation of a set of non-canonical ORFs identified in GENCODE transcripts, the LRGASP collaboration to assess the use of long transcriptomic data to build transcript models, the progress in collaborations with RefSeq and UniProt to increase convergence in the annotation of human and mouse protein-coding genes, the propagation of GENCODE across the human pan-genome and the development of new tools to support annotation of regulatory features by GENCODE. Our annotation is accessible via Ensembl, the UCSC Genome Browser and https://www.gencodegenes.org.

A proteome-scale map of the SARS-CoV-2-human contactome

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

A proof-of-concept study on the genomic evolution of Sars-Cov-2 in molnupiravir-treated, paxlovid-treated and drug-naïve patients

Little is known about SARS-CoV-2 evolution under Molnupiravir and Paxlovid, the only antivirals approved for COVID-19 treatment. By investigating SARS-CoV-2 variability in 8 Molnupiravir-treated, 7 Paxlovid-treated and 5 drug-naïve individuals at 4 time-points (Days 0-2-5-7), a higher genetic distance is found under Molnupiravir pressure compared to Paxlovid and no-drug pressure (nucleotide-substitutions/site mean±Standard error: 18.7 × 10-4 ± 2.1 × 10-4 vs. 3.3 × 10-4 ± 0.8 × 10-4 vs. 3.1 × 10-4 ± 0.8 × 10-4, P = 0.0003), peaking between Day 2 and 5. Molnupiravir drives the emergence of more G-A and C-T transitions than other mutations (P = 0.031). SARS-CoV-2 selective evolution under Molnupiravir pressure does not differ from that under Paxlovid or no-drug pressure, except for orf8 (dN > dS, P = 0.001); few amino acid mutations are enriched at specific sites. No RNA-dependent RNA polymerase (RdRp) or main proteases (Mpro) mutations conferring resistance to Molnupiravir or Paxlovid are found. This proof-of-concept study defines the SARS-CoV-2 within-host evolution during antiviral treatment, confirming higher in vivo variability induced by Molnupiravir compared to Paxlovid and drug-naive, albeit not resulting in apparent mutation selection.

Cell Entry and Unusual Replication of SARS-CoV-2

BACKGROUND

SARS-CoV-2 is the causative virus for the CoVID-19 pandemic that has frequently mutated to continue to infect and resist available vaccines. Emerging new variants of the virus have complicated notions of immunity conferred by vaccines versus immunity that results from infection. While we continue to progress from epidemic to endemic as a result of this collective immunity, the pandemic remains a morbid and mortal problem.

OBJECTIVE

The SARS-CoV-2 virus has a very complex manner of replication. The spike protein, one of the four structural proteins of the encapsulated virus, is central to the ability of the virus to penetrate cells to replicate. The objective of this review is to summarize these complex features of viral replication.

METHODS

A review of the recent literature was performed on the biology of SARS-CoV-2 infection from published work from PubMed and works reported to preprint servers, e.g., bioRxiv and medRxiv.

RESULTS AND CONCLUSION

The complex molecular and cellular biology involved in SARS-CoV-2 replication and the origination of >30 proteins from a single open reading frame (ORF) have been summarized, as well as the structural biology of spike protein, a critical factor in the cellular entry of the virus, which is a necessary feature for it to replicate and cause disease.

Focus on Marine Animal Safety and Marine Bioresources in Response to the SARS-CoV-2 Crisis

SARS-CoV-2 as a zoonotic virus has significantly affected daily life and social behavior since its outbreak in late 2019. The concerns over its transmission through different media directly or indirectly have evoked great attention about the survival of SARS-CoV-2 virions in the environment and its potential infection of other animals. To evaluate the risk of infection by SARS-CoV-2 and to counteract the COVID-19 disease, extensive studies have been performed to understand SARS-CoV-2 biogenesis and its pathogenesis. This review mainly focuses on the molecular architecture of SARS-CoV-2, its potential for infecting marine animals, and the prospect of drug discovery using marine natural products to combat SARS-CoV-2. The main purposes of this review are to piece together progress in SARS-CoV-2 functional genomic studies and antiviral drug development, and to raise our awareness of marine animal safety on exposure to SARS-CoV-2.

Hydrophobic Alpha-Helical Short Peptides in Overlapping Reading Frames of the Coronavirus Genome

In this study, we show that the coronavirus (CoV) genome may encode many functional hydrophobic alpha-helical peptides (HAHPs) in overlapping reading frames of major coronaviral proteins throughout the entire viral genome. These HAHPs can theoretically be expressed from non-canonical sub-genomic (sg)RNAs that are synthesized in substantial amounts in infected cells. We selected and analyzed five and six HAHPs encoded in the S gene regions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and Middle East respiratory syndrome coronavirus (MERS-CoV), respectively. Two and three HAHPs derived from SARS-CoV-2 and MERS-CoV, respectively, specifically interacted with both the SARS-CoV-2 and MERS-CoV S proteins and inhibited their membrane fusion activity. Furthermore, one of the SARS-CoV-2 HAHPs specifically inhibited viral RNA synthesis by accumulating at the site of viral RNA synthesis. Our data show that a group of HAHPs in the coronaviral genome potentially has a regulatory role in viral propagation.

SARS-CoV-2 ORF6 disrupts innate immune signalling by inhibiting cellular mRNA export

SARS-CoV-2 is a betacoronavirus and the etiological agent of COVID-19, a devastating infectious disease. Due to its far-reaching effect on human health, there is an urgent and growing need to understand the viral molecular biology of SARS-CoV-2 and its interaction with the host cell. SARS-CoV-2 encodes 9 predicted accessory proteins, which are presumed to be dispensable for in vitro replication, most likely having a role in modulating the host cell environment to aid viral replication. Here we show that the ORF6 accessory protein interacts with cellular Rae1 to inhibit cellular protein production by blocking mRNA export. We utilised cell fractionation coupled with mRNAseq to explore which cellular mRNA species are affected by ORF6 expression and show that ORF6 can inhibit the export of many mRNA including those encoding antiviral factors such as IRF1 and RIG-I. We also show that export of these mRNA is blocked in the context of SARS-CoV-2 infection. Together, our studies identify a novel mechanism by which SARS-CoV-2 can manipulate the host cell environment to supress antiviral responses, providing further understanding to the replication strategies of a virus that has caused an unprecedented global health crisis.

Convolutional Neural Network Applied to SARS-CoV-2 Sequence Classification

COVID-19, the illness caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus belonging to the Coronaviridade family, a single-strand positive-sense RNA genome, has been spreading around the world and has been declared a pandemic by the World Health Organization. On 17 January 2022, there were more than 329 million cases, with more than 5.5 million deaths. Although COVID-19 has a low mortality rate, its high capacities for contamination, spread, and mutation worry the authorities, especially after the emergence of the Omicron variant, which has a high transmission capacity and can more easily contaminate even vaccinated people. Such outbreaks require elucidation of the taxonomic classification and origin of the virus (SARS-CoV-2) from the genomic sequence for strategic planning, containment, and treatment of the disease. Thus, this work proposes a high-accuracy technique to classify viruses and other organisms from a genome sequence using a deep learning convolutional neural network (CNN). Unlike the other literature, the proposed approach does not limit the length of the genome sequence. The results show that the novel proposal accurately distinguishes SARS-CoV-2 from the sequences of other viruses. The results were obtained from 1557 instances of SARS-CoV-2 from the National Center for Biotechnology Information (NCBI) and 14,684 different viruses from the Virus-Host DB. As a CNN has several changeable parameters, the tests were performed with forty-eight different architectures; the best of these had an accuracy of 91.94 ± 2.62% in classifying viruses into their realms correctly, in addition to 100% accuracy in classifying SARS-CoV-2 into its respective realm, Riboviria. For the subsequent classifications (family, genera, and subgenus), this accuracy increased, which shows that the proposed architecture may be viable in the classification of the virus that causes COVID-19.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

Homology-based classification of accessory proteins in coronavirus genomes uncovers extremely dynamic evolution of gene content

Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super‐groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike‐like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion‐associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross‐communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.

SARS-CoV-2 Nsp14 mediates the effects of viral infection on the host cell transcriptome

Viral infection involves complex set of events orchestrated by multiple viral proteins. To identify functions of SARS-CoV-2 proteins, we performed transcriptomic analyses of cells expressing individual viral proteins. Expression of Nsp14, a protein involved in viral RNA replication, provoked a dramatic remodeling of the transcriptome that strongly resembled that observed following SARS-CoV-2 infection. Moreover, Nsp14 expression altered the splicing of more than 1000 genes and resulted in a dramatic increase in the number of circRNAs, which are linked to innate immunity. These effects were independent of the Nsp14 exonuclease activity and required the N7-guanine-methyltransferase domain of the protein. Activation of the NFkB pathway and increased expression of CXCL8 occurred early upon Nsp14 expression. We identified IMPDH2, which catalyzes the rate-limiting step of guanine nucleotides biosynthesis, as a key mediator of these effects. Nsp14 expression caused an increase in GTP cellular levels, and the effect of Nsp14 was strongly decreased in the presence of IMPDH2 inhibitors. Together, our data demonstrate an unknown role for Nsp14 with implications for therapy.

Viruses are parasites, relying on the cells they infect to make more of themselves. In doing so they change how an infected cell turns its genes on and off, forcing it to build new virus particles and turning off the immune surveillance that would allow the body to intervene. This is how SARS-CoV-2, the virus that causes COVID, survives with a genome that carries instructions to make just 29 proteins.

One of these proteins, known as Nsp14, is involved in both virus reproduction and immune escape. Previous work has shown that it interacts with IMPDH2, the cellular enzyme that controls the production of the building blocks of the genetic code. The impact of this interaction is not clear.

To find out more, Zaffagni et al. introduced 26 of the SARS-CoV-2 proteins into human cells one at a time. Nsp14 had the most dramatic effect, dialing around 4,000 genes up or down and changing how the cell interprets over 1,000 genes. Despite being just one protein, it mimicked the genetic changes seen during real SARS-CoV-2 infection. Blocking IMPDH2 partially reversed the effects, which suggests that the interaction of Nsp14 with the enzyme might be responsible for the effects of SARS-CoV-2 on the genes of the cell.

Understanding how viral proteins affect cells can explain what happens during infection. This could lead to the discovery of new treatments designed to counteract the effects of the virus. Further work could investigate whether interfering with Nsp14 helps cells to overcome infection.

Understanding the Role of SARS-CoV-2 ORF3a in Viral Pathogenesis and COVID-19

The ongoing SARS-CoV-2 pandemic has shocked the world due to its persistence, COVID-19-related morbidity and mortality, and the high mutability of the virus. One of the major concerns is the emergence of new viral variants that may increase viral transmission and disease severity. In addition to mutations of spike protein, mutations of viral proteins that affect virulence, such as ORF3a, also must be considered. The purpose of this article is to review the current literature on ORF3a, to summarize the molecular actions of SARS-CoV-2 ORF3a, and its role in viral pathogenesis and COVID-19. ORF3a is a polymorphic, multifunctional viral protein that is specific to SARS-CoV/SARS-CoV-2. It was acquired from β-CoV lineage and likely originated from bats through viral evolution. SARS-CoV-2 ORF3a is a viroporin that interferes with ion channel activities in host plasma and endomembranes. It is likely a virion-associated protein that exerts its effect on the viral life cycle during viral entry through endocytosis, endomembrane-associated viral transcription and replication, and viral release through exocytosis. ORF3a induces cellular innate and pro-inflammatory immune responses that can trigger a cytokine storm, especially under hypoxic conditions, by activating NLRP3 inflammasomes, HMGB1, and HIF-1α to promote the production of pro-inflammatory cytokines and chemokines. ORF3a induces cell death through apoptosis, necrosis, and pyroptosis, which leads to tissue damage that affects the severity of COVID-19. ORF3a continues to evolve along with spike and other viral proteins to adapt in the human cellular environment. How the emerging ORF3a mutations alter the function of SARS-CoV-2 ORF3a and its role in viral pathogenesis and COVID-19 is largely unknown. This review provides an in-depth analysis of ORF3a protein's structure, origin, evolution, and mutant variants, and how these characteristics affect its functional role in viral pathogenesis and COVID-19.

Antiviral effects of azithromycin: A narrative review

Viral infections have a great impact on human health. The urgent need to find a cure against different viruses led us to investigations in a vast range of drugs. Azithromycin (AZT), classified as a macrolide, showed various effects on different known viruses such as severe acute respiratory syndrome coronavirus (SARS-CoV), Zika, Ebola, Enterovirus (EVs) and Rhinoviruses (RVs), and Influenza A previously; namely, these viruses, which caused global concerns, are considered as targets for AZT different actions. Due to AZT background in the treatment of known viral infections mentioned above (which is described in this study), in the early stages of COVID-19 (a new zoonotic disease caused by a novel coronavirus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)) development, AZT drew attention to itself due to its antiviral and immunomodulatory effects as a valuable candidate for COVID-19 treatment. AZT usage instructions for treating different viral infections have always been under observation, and COVID-19 is no exception. There are still debates about the use of AZT in COVID-19 treatment. However, eventually, novel researches convinced WHO to announce the discontinuation of AZT use (alone or in combination with hydroxychloroquine) in treating SARS-CoV-2 infection. This research aims to study the structure of all of the viruses mentioned above and the molecular and clinical effects of AZT against the virus.

Evolution of enhanced innate immune evasion by SARS-CoV-2

The emergence of SARS-CoV-2 variants of concern suggests viral adaptation to enhance human-to-human transmission1,2. Although much effort has focused on the characterization of changes in the spike protein in variants of concern, mutations outside of spike are likely to contribute to adaptation. Here, using unbiased abundance proteomics, phosphoproteomics, RNA sequencing and viral replication assays, we show that isolates of the Alpha (B.1.1.7) variant3 suppress innate immune responses in airway epithelial cells more effectively than first-wave isolates. We found that the Alpha variant has markedly increased subgenomic RNA and protein levels of the nucleocapsid protein (N), Orf9b and Orf6-all known innate immune antagonists. Expression of Orf9b alone suppressed the innate immune response through interaction with TOM70, a mitochondrial protein that is required for activation of the RNA-sensing adaptor MAVS. Moreover, the activity of Orf9b and its association with TOM70 was regulated by phosphorylation. We propose that more effective innate immune suppression, through enhanced expression of specific viral antagonist proteins, increases the likelihood of successful transmission of the Alpha variant, and may increase in vivo replication and duration of infection4. The importance of mutations outside the spike coding region in the adaptation of SARS-CoV-2 to humans is underscored by the observation that similar mutations exist in the N and Orf9b regulatory regions of the Delta and Omicron variants.

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.

Structural topological analysis of spike proteins of SARS-CoV-2 variants of concern highlight distinctive amino acid substitution patterns

Since the onset of pandemic in 2019, SARS-CoV-2 has diverged into numerous variants driven by antigenic and infectivity-oriented selection. Some variants have accumulated fitness-enhancing mutations, evaded immunity and spread despite global vaccination campaigns. The spike (S) glycoprotein of SARS-CoV-2 demonstrated the greatest immunogenicity and amino acid substitution diversity owing to its importance in the interaction with human angiotensin receptor 2 (hACE2). The S protein consistently emerges as an amino acid substitution (AAS) hotspot in all six lineages, however, in Omicron this enrichment is significantly higher. This study attempts to design and validate a method of mapping S-protein substitution profile across variants to identify the conserved and AAS regions. A substitution matrix was created based on publicly available databases, and the substitution localization was illustrated on a cryo-electron microscopy generated S-protein model. Our analyses indicated that the diversity of N-terminal (NTD) and receptor-binding (RBD) domains exceeded that of any other regions but still contained extended low substitution density regions particularly considering significantly broader substitution profiles of Omicron BA.2 and BA.4/5. Finally, the substitution matrix was compared to a random sample alignment of variant sequences, revealing discrepancies. Therefore, it was suggested to improve matrix accuracy by processing a large number of S-protein sequences using an automated algorithm. Several critical immunogenic and receptor-interacting residues were identified in the conserved regions within NTD and RBD. In conclusion, the structural and topological analysis of S proteins of SARS-CoV-2 variants highlight distinctive amino acid substitution patterns which may be foundational in predicting future variants.

Interaction between Sars-CoV-2 structural proteins and host cellular receptors: From basic mechanisms to clinical perspectives

Severe acute respiratory syndrome coronavirus 2 (Sars-CoV-2) has caused a global pandemic that has affected the lives of billions of individuals. Sars-CoV-2 primarily infects human cells by binding of the viral spike protein to angiotensin-converting enzyme 2 (ACE2). In addition, novel means of viral entry are currently being investigated, including Neuropillin 1, toll-like receptors (TLRs), cluster of differentiation 147 (CD147), and integrin α5β1. Enriched expression of these proteins across metabolic regulatory organs/tissues, including the circulatory system, liver, pancreas, and intestine contributes to major clinical complications among COVID-19 patients, particularly the development of hypertension, myocardial injury, arrhythmia, acute coronary syndrome and increased coagulation in the circulatory system during and post-infection. Pre-existing metabolic disease, such as cardiovascular disease, obesity, diabetes, and non-alcoholic fatty liver disease, is associated with increased risk of hospitalization, persistent post-infection complications and worse outcomes in patients with COVID-19. This review overviews the biological features of Sars-CoV-2, highlights recent findings that delineate the pathological mechanisms of COVID-19 and the consequent clinical diseases.

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.

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.

Profiling SARS-CoV-2 HLA-I peptidome reveals T cell epitopes from out-of-frame ORFs

T cell-mediated immunity plays an important role in controlling SARS-CoV-2 infection, but the repertoire of naturally processed and presented viral epitopes on class I human leukocyte antigen (HLA-I) remains uncharacterized. Here, we report the first HLA-I immunopeptidome of SARS-CoV-2 in two cell lines at different times post infection using mass spectrometry. We found HLA-I peptides derived not only from canonical open reading frames (ORFs) but also from internal out-of-frame ORFs in spike and nucleocapsid not captured by current vaccines. Some peptides from out-of-frame ORFs elicited T cell responses in a humanized mouse model and individuals with COVID-19 that exceeded responses to canonical peptides, including some of the strongest epitopes reported to date. Whole-proteome analysis of infected cells revealed that early expressed viral proteins contribute more to HLA-I presentation and immunogenicity. These biological insights, as well as the discovery of out-of-frame ORF epitopes, will facilitate selection of peptides for immune monitoring and vaccine development.

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.

SARS-CoV-2 uses a multipronged strategy to impede host protein synthesis

The coronavirus SARS-CoV-2 is the cause of the ongoing pandemic of COVID-19¹. Coronaviruses have developed a variety of mechanisms to repress host mRNA translation to allow the translation of viral mRNA, and concomitantly block the cellular innate immune response^2,3. Although several different proteins of SARS-CoV-2 have previously been implicated in shutting off host expression^4-7, a comprehensive picture of the effects of SARS-CoV-2 infection on cellular gene expression is lacking. Here we combine RNA sequencing, ribosome profiling and metabolic labelling of newly synthesized RNA to comprehensively define the mechanisms that are used by SARS-CoV-2 to shut off cellular protein synthesis. We show that infection leads to a global reduction in translation, but that viral transcripts are not preferentially translated. Instead, we find that infection leads to the accelerated degradation of cytosolic cellular mRNAs, which facilitates viral takeover of the mRNA pool in infected cells. We reveal that the translation of transcripts that are induced in response to infection (including innate immune genes) is impaired. We demonstrate this impairment is probably mediated by inhibition of nuclear mRNA export, which prevents newly transcribed cellular mRNA from accessing ribosomes. Overall, our results uncover a multipronged strategy that is used by SARS-CoV-2 to take over the translation machinery and to suppress host defences.

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.