Pervasive generation of non-canonical subgenomic RNAs by SARS-CoV-2

Cov-trans: an efficient algorithm for discontinuous transcript assembly in coronaviruses

BACKGROUND

Discontinuous transcription allows coronaviruses to efficiently replicate and transmit within host cells, enhancing their adaptability and survival. Assembling viral transcripts is crucial for virology research and the development of antiviral strategies. However, traditional transcript assembly methods primarily designed for variable alternative splicing events in eukaryotes are not suitable for the viral transcript assembly problem. The current algorithms designed for assembling viral transcripts often struggle with low accuracy in determining the transcript boundaries. There is an urgent need to develop a highly accurate viral transcript assembly algorithm.

RESULTS

In this work, we propose Cov-trans, a reference-based transcript assembler specifically tailored for the discontinuous transcription of coronaviruses. Cov-trans first identifies canonical transcripts based on discontinuous transcription mechanisms, start and stop codons, as well as reads alignment information. Subsequently, it formulates the assembly of non-canonical transcripts as a path extraction problem, and introduces a mixed integer linear programming to recover these non-canonical transcripts.

CONCLUSION

Experimental results show that Cov-trans outperforms other assemblers in both accuracy and recall, with a notable strength in accurately identifying the boundaries of transcripts. Cov-trans is freely available at https://github.com/computer-Bioinfo/Cov-trans.git .

Improved sub-genomic RNA prediction with the ARTIC protocol

Viral subgenomic RNA (sgRNA) plays a major role in SARS-COV2's replication, pathogenicity, and evolution. Recent sequencing protocols, such as the ARTIC protocol, have been established. However, due to the viral-specific biological processes, analyzing sgRNA through viral-specific read sequencing data is a computational challenge. Current methods rely on computational tools designed for eukaryote genomes, resulting in a gap in the tools designed specifically for sgRNA detection. To address this, we make two contributions. Firstly, we present sgENERATE, an evaluation pipeline to study the accuracy and efficacy of sgRNA detection tools using the popular ARTIC sequencing protocol. Using sgENERATE, we evaluate periscope, a recently introduced tool that detects sgRNA from ARTIC sequencing data. We find that periscope has biased predictions and high computational costs. Secondly, using the information produced from sgENERATE, we redesign the algorithm in periscope to use multiple references from canonical sgRNAs to mitigate alignment issues and improve sgRNA and non-canonical sgRNA detection. We evaluate periscope and our algorithm, periscope_multi, on simulated and biological sequencing datasets and demonstrate periscope_multi's enhanced sgRNA detection accuracy. Our contribution advances tools for studying viral sgRNA, paving the way for more accurate and efficient analyses in the context of viral RNA discovery.

Cellular dynamics shape recombination frequency in coronaviruses

Coronavirus genomes have evolutionary histories shaped extensively by recombination. Yet, how often recombination occurs at a cellular level, or the factors that regulate recombination rates, are poorly understood. Utilizing experimental co-infections with pairs of genetically distinct coronaviruses, we found that recombination is both frequent and rare during coinfection. Recombination occurred in every instance of co-infection yet resulted in relatively few recombinant RNAs. By integrating a discrete-time Susceptible-Infected-Removed (SIR) model, we found that rates of recombination are determined primarily by rates of cellular co-infection, rather than other possible barriers such as RNA compartmentalization. By staggering the order and timing of infection with each virus we also found that rates of co-infection are themselves heavily influenced by genetic and ecological mechanisms, including superinfection exclusion and the relative fitness of competing viruses. Our study highlights recombination as a potent yet regulated force: frequent enough to ensure a steady influx of genetic variation but also infrequent enough to maintain genomic integrity. As recombination is thought to be an important driver of host-switching and disease emergence, our study provides new insights into the factors that regulate coronavirus recombination and evolution more broadly.

Adaptive truncation of the S gene in IBV during chicken embryo passaging plays a crucial role in its attenuation

Like all coronaviruses, infectious bronchitis virus, the causative agent of infectious bronchitis in chickens, exhibits a high mutation rate. Adaptive mutations that arise during the production of live attenuated vaccines against IBV often decrease virulence. The specific impact of these mutations on viral pathogenicity, however, has not been fully elucidated. In this study, we identified a mutation at the 3' end of the S gene in an IBV strain that was serially passaged in chicken embryos, and showed that this mutation resulted in a 9-aa truncation of the cytoplasmic tail (CT) of the S protein. This phenomenon of CT truncation has previously been observed in the production of attenuated vaccines against other coronaviruses such as the porcine epidemic diarrhea virus. We next discovered that the 9-aa truncation in the S protein CT resulted in the loss of the endoplasmic-reticulum-retention signal (KKSV). Rescue experiments with recombinant viruses confirmed that the deletion of the KKSV motif impaired the localization of the S protein to the endoplasmic-reticulum-Golgi intermediate compartment (ERGIC) and increased its expression on the cell surface. This significantly reduced the incorporation of the S protein into viral particles, impaired early subgenomic RNA and protein synthesis, and ultimately reduced viral invasion efficiency in CEK cells. In vivo experiments in chickens confirmed the reduced pathogenicity of the mutant IBV strains. Additionally, we showed that the adaptive mutation altered the TRS-B of ORF3 and impacted the transcriptional regulation of this gene. Our findings underscore the significance of this adaptive mutation in the attenuation of IBV infection and provide a novel strategy for the development of live attenuated IBV vaccines.

Interstitial macrophages are a focus of viral takeover and inflammation in COVID-19 initiation in human lung

Early stages of deadly respiratory diseases including COVID-19 are challenging to elucidate in humans. Here, we define cellular tropism and transcriptomic effects of SARS-CoV-2 virus by productively infecting healthy human lung tissue and using scRNA-seq to reconstruct the transcriptional program in "infection pseudotime" for individual lung cell types. SARS-CoV-2 predominantly infected activated interstitial macrophages (IMs), which can accumulate thousands of viral RNA molecules, taking over 60% of the cell transcriptome and forming dense viral RNA bodies while inducing host profibrotic (TGFB1, SPP1) and inflammatory (early interferon response, CCL2/7/8/13, CXCL10, and IL6/10) programs and destroying host cell architecture. Infected alveolar macrophages (AMs) showed none of these extreme responses. Spike-dependent viral entry into AMs used ACE2 and Sialoadhesin/CD169, whereas IM entry used DC-SIGN/CD209. These results identify activated IMs as a prominent site of viral takeover, the focus of inflammation and fibrosis, and suggest targeting CD209 to prevent early pathology in COVID-19 pneumonia. This approach can be generalized to any human lung infection and to evaluate therapeutics.

Identification of Hammerhead-variant ribozyme sequences in SARS-CoV-2

The SARS-CoV-2 RNA virus and variants, responsible for the COVID-19 pandemic has become endemic, raised a need for further understanding of the viral genome and biology. Despite vast research on SARS-CoV-2, no ribozymes have been found in the virus genome. Here we report the identification of 39 Hammerhead-variant ribozyme sequences (CoV-HHRz) in SARS-CoV-2. These sequences are highly conserved within SARS-CoV-2 variants but show large diversity among other coronaviruses. In vitro CoV-HHRz sequences possess the characteristics of typical ribozymes; cleavage is pH and ion dependent, although their activity is relatively low and Mn2+ is required for cleavage. The cleavage sites of four CoV-HHRz coincide with the breakpoint of expressed subgenomic RNA (sgRNAs) in SARS-CoV-2 transcriptome data suggesting in vivo activity. The CoV-HHRz are involved in processing sgRNAs for ORF7b, ORF 10 and ORF1ab nsp13 which are essential for viral packaging and life cycle.

Experimental and analytical pipeline for sub-genomic RNA landscape of coronavirus by Nanopore sequencer

Coronaviruses (CoVs), including severe acute respiratory syndrome coronavirus 2, can infect a variety of mammalian and avian hosts with significant medical and economic consequences. During the life cycle of CoV, a coordinated series of subgenomic RNAs, including canonical subgenomic messenger RNA and non-canonical defective viral genomes (DVGs), are generated with different biological implications. Studies that adopted the Nanopore sequencer (ONT) to investigate the landscape and dynamics of viral RNA subgenomic transcriptomes applied arbitrary bioinformatics parameters without justification or experimental validation. The current study used bovine coronavirus (BCoV), which can be performed under biosafety level 2 for library construction and experimental validation using traditional colony polymerase chain reaction and Sanger sequencing. Four different ONT protocols, including RNA direct and cDNA direct sequencing with or without exonuclease treatment, were used to generate RNA transcriptomic libraries from BCoV-infected cell lysates. Through rigorously examining the k-mer, gap size, segment size, and bin size, the optimal cutoffs for the bioinformatic pipeline were determined to remove the sequence noise while keeping the informative DVG reads. The sensitivity and specificity of identifying DVG reads using the proposed pipeline can reach 82.6% and 99.6% under the k-mer size cutoff of 15. Exonuclease treatment reduced the abundance of RNA transcripts; however, it was not necessary for future library preparation. Additional recovery of clipped BCoV nucleotide sequences with experimental validation expands the landscape of the CoV discontinuous RNA transcriptome, whose biological function requires future investigation. The results of this study provide the benchmarks for library construction and bioinformatic parameters for studying the discontinuous CoV RNA transcriptome.IMPORTANCEFunctional defective viral genomic RNA, containing all the cis-acting elements required for translation or replication, may play different roles in triggering cell innate immune signaling, interfering with the canonical subgenomic messenger RNA transcription/translation or assisting in establishing persistence infection. This study does not only provide benchmarks for library construction and bioinformatic parameters for studying the discontinuous coronavirus RNA transcriptome but also reveals the complexity of the bovine coronavirus transcriptome, whose functional assays will be critical in future studies.

An RNA-Seq analysis of coronavirus in the skin of the Pangolin

Protection of the Critically Endangered East Asian Pangolin species is hampered by the vulnerability of captive individuals to infection. Studies have previously shown the pangolin to have a unique pseudogenisation of many immunity genes (including IFNE, IFIH1, cGAS, STING, TLR5, and TLR11), and we suspected that these losses could account for this vulnerability. Here we used RNA-Seq data to show the effect of these gene losses on the transcriptional response to a viral skin infection in a deceased pangolin. This virus is very closely related to the one causing the current COVID-19 pandemic in the human population (SARS-CoV2), and we found the most upregulated pathway was the same one previously identified in the lungs of SARS-CoV2-infected humans. As predicted, we found that the pathways downstream of the lost genes were not upregulated. For example, the pseudogenised interferon epsilon (IFNE) is known to be particularly important in epithelial immunity, and we show that interferon-related responses were not upregulated in the infected pangolin skin. We suggest that the pangolin's innate gene pseudogenisation is indeed likely to be responsible for the animal's vulnerability to infection.

Expression dynamics of the aplysia abyssovirus

Two recent studies documented the genome of a novel, extremely large (35.9 kb), nidovirus in RNA sequence databases from the marine neural model Aplysia californica. The goal of the present study was to document the distribution and transcriptional dynamics of this virus, Aplysia abyssovirus 1 (AAbV), in maricultured and wild animals. We confirmed previous findings that AAbV RNA is widespread and reaches extraordinary levels in apparently healthy animals. Transmission electron microscopy identified viral replication factories in ciliated gill epithelial cells but not in neurons where viral RNA is most highly expressed. Viral transcripts do not exhibit evidence of discontinuous RNA synthesis as in coronaviruses but are consistent with production of a single leaderless subgenomic RNA, as in the Gill-associated virus of Penaeus monodon. Splicing patterns in chronically infected adults suggested high levels of defective genomes, possibly explaining the lack of obvious disease signs in high viral load animals.

Identification of the protein coding capability of coronavirus defective viral genomes by mass spectrometry

During coronavirus infection, in addition to the well-known coronavirus genomes and subgenomic mRNAs, an abundance of defective viral genomes (DVGs) can also be synthesized. In this study, we aimed to examine whether DVGs can encode proteins in infected cells. Nanopore direct RNA sequencing and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis were employed. With the protein databases generated by nanopore direct RNA sequencing and the cell lysates derived from the RNA-protein pull-down assay, six DVG-encoded proteins were identified by LC-MS/MS based on the featured fusion peptides caused by recombination during DVG synthesis. The results suggest that the coronavirus DVGs have the capability to encode proteins. Consequently, future studies determining the biological function of DVG-encoded proteins may contribute to the understanding of their roles in coronavirus pathogenesis and the development of antiviral strategies.

Comparative analysis of bioinformatics tools to characterize SARS-CoV-2 subgenomic RNAs

During the replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), positive-sense genomic RNA and subgenomic RNAs (sgRNAs) are synthesized by a discontinuous process of transcription characterized by a template switch, regulated by transcription-regulating sequences (TRS). Although poorly known about makeup and dynamics of sgRNAs population and function of its constituents, next-generation sequencing approaches with the help of bioinformatics tools have made a significant contribution to expand the knowledge of sgRNAs in SARS-CoV-2. For this scope to date, Periscope, LeTRS, sgDI-tector, and CORONATATOR have been developed. However, limited number of studies are available to compare the performance of such tools. To this purpose, we compared Periscope, LeTRS, and sgDI-tector in the identification of canonical (c-) and noncanonical (nc-) sgRNA species in the data obtained with the Illumina ARTIC sequencing protocol applied to SARS-CoV-2-infected Caco-2 cells, sampled at different time points. The three software showed a high concordance rate in the identification and in the quantification of c-sgRNA, whereas more differences were observed in nc-sgRNA. Overall, LeTRS and sgDI-tector result to be adequate alternatives to Periscope to analyze Fastq data from sequencing platforms other than Nanopore.

Unveiling the biology of defective viral genomes in vitro and in vivo: implications for gene expression and pathogenesis of coronavirus

BACKGROUND

Defective viral genome (DVG) is a truncated version of the full-length virus genome identified in most RNA viruses during infection. The synthesis of DVGs in coronavirus has been suggested; however, the fundamental characteristics of coronavirus DVGs in gene expression and pathogenesis have not been systematically analyzed.

METHODS

Nanopore direct RNA sequencing was used to investigate the characteristics of coronavirus DVGs in gene expression including reproducibility, abundance, species and genome structures for bovine coronavirus in cells, and for mouse hepatitis virus (MHV)-A59 (a mouse coronavirus) in cells and in mice. The MHV-A59 full-length genomic cDNAs (~ 31 kilobases) were in vitro constructed to experimentally validate the origin of coronavirus DVG. The synthesis of DVGs was also experimentally identified by RT-PCR followed by sequencing. In addition, the alterations of DVGs in amounts and species under different infection environments and selection pressures including the treatment of antiviral remdesivir and interferon were evaluated based on the banding patterns by RT-PCR.

RESULTS

The results are as follows: (i) the structures of DVGs are with diversity, (ii) DVGs are overall synthesized with moderate (MHV-A59 in cells) to high (BCoV in cells and MHV-A59 in mice) reproducibility under regular infection with the same virus inoculum, (iii) DVGs can be synthesized from the full-length coronavirus genome, (iv) the sequences flanking the recombination point of DVGs are AU-rich and thus may contribute to the recombination events during gene expression, (v) the species and amounts of DVG are altered under different infection environments, and (vi) the biological nature of DVGs between in vitro and in vivo is similar.

CONCLUSIONS

The identified biological characteristics of coronavirus DVGs in terms of abundance, reproducibility, and variety extend the current model for coronavirus gene expression. In addition, the biological features of alterations in amounts and species of coronavirus DVGs under different infection environments may assist the coronavirus to adapt to the altered environments for virus fitness and may contribute to the coronavirus pathogenesis. Consequently, the unveiled biological features may assist the community to study the gene expression mechanisms of DVGs and their roles in pathogenesis, contributing to the development of antiviral strategy and public health.

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed "discontinuous transcription" that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

Emerging role of microRNAs and long non-coding RNAs in COVID-19 with implications to therapeutics

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection which is commonly known as COVID-19 (COronaVIrus Disease 2019) has creeped into the human population taking tolls of life and causing tremendous economic crisis. It is indeed crucial to gain knowledge about their characteristics and interactions with human host cells. It has been shown that the majority of our genome consists of non-coding RNAs. Non-coding RNAs including micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs) display significant roles in regulating gene expression in almost all cancers and viral diseases. It is intriguing that miRNAs and lncRNAs remarkably regulate the function and expression of major immune components of SARS-CoV-2. MiRNAs act via RNA interference mechanism in which they bind to the complementary sequences of the viral RNA strand, inducing the formation of silencing complex that eventually degrades or inhibits the viral RNA and viral protein expression. LncRNAs have been extensively shown to regulate gene expression in cytokine storm and thus emerges as a critical target for COVID-19 treatment. These lncRNAs also act as competing endogenous RNAs (ceRNAs) by sponging miRNAs and thus affecting the expression of downstream targets during SARS-CoV-2 infection. In this review, we extensively discuss the role of miRNAs and lncRNAs, describe their mechanism of action and their different interacting human targets cells during SARS-CoV-2 infection. Finally, we discuss possible ways how an interference with their molecular function could be exploited for new therapies against SARS-CoV-2.

sgRNAs: A SARS-CoV-2 emerging issue

Like for other coronaviruses, SARS-CoV-2 gene expression strategy is based on the synthesis of a nested set of subgenomic mRNA species (sgRNAs). These sgRNA are synthesized using a "discontinuous transcription" mechanism that relies on template switching at Transcription Regulatory Sequences (TRS). Both canonical (c-sgRNA) and non-canonical (nc-sgRNA, less numerous) subgenomic RNA species can be produced. Currently, sgRNAs are investigated on the basis of sequence data obtained through next generation sequencing (NGS), and bioinformatic tools are crucial for their identification, characterization and quantification. To date, few software have been developed to this aim, whose reliability and applicability to all the available NGS platforms need to be established, to build confidence on the information resulting from such tools. In fact, these information may be crucial for the in depth elucidation of viral expression strategy, particularly in respect of the significance of nc-sgRNAs, and for the possible use of sgRNAs as potential markers of virus replicative activity in infected patients.

Intragenomic rearrangements involving 5'-untranslated region segments in SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses

BACKGROUND

Variation of the betacoronavirus SARS-CoV-2 has been the bane of COVID-19 control. Documented variation includes point mutations, deletions, insertions, and recombination among closely or distantly related coronaviruses. Here, we describe yet another aspect of genome variation by beta- and alphacoronaviruses that was first documented in an infectious isolate of the betacoronavirus SARS-CoV-2, obtained from 3 patients in Hong Kong that had a 5'-untranslated region segment at the end of the ORF6 gene that in its new location translated into an ORF6 protein with a predicted modified carboxyl terminus. While comparing the amino acid sequences of translated ORF8 genes in the GenBank database, we found a subsegment of the same 5'-UTR-derived amino acid sequence modifying the distal end of ORF8 of an isolate from the United States and decided to carry out a systematic search.

METHODS

Using the nucleotide and in the case of SARS-CoV-2 also the translated amino acid sequence in three reading frames of the genomic termini of coronaviruses as query sequences, we searched for 5'-UTR sequences in regions other than the 5'-UTR in SARS-CoV-2 and reference strains of alpha-, beta-, gamma-, and delta-coronaviruses.

RESULTS

We here report numerous genomic insertions of 5'-untranslated region sequences into coding regions of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses, but not delta- or gammacoronaviruses. To our knowledge this is the first systematic description of such insertions. In many cases, these insertions would change viral protein sequences and further foster genomic flexibility and viral adaptability through insertion of transcription regulatory sequences in novel positions within the genome. Among human Embecorivus betacoronaviruses, for instance, from 65% to all of the surveyed sequences in publicly available databases contain inserted 5'-UTR sequences.

CONCLUSION

The intragenomic rearrangements involving 5'-untranslated region sequences described here, which in several cases affect highly conserved genes with a low propensity for recombination, may underlie the generation of variants homotypic with those of concern or interest and with potentially differing pathogenic profiles. Intragenomic rearrangements thus add to our appreciation of how variants of SARS-CoV-2 and other beta- and alphacoronaviruses may arise.

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.

Subgenomic RNAs and Their Encoded Proteins Contribute to the Rapid Duplication of SARS-CoV-2 and COVID-19 Progression

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently widespread throughout the world, accompanied by a rising number of people infected and breakthrough infection of variants, which make the virus highly transmissible and replicable. A comprehensive understanding of the molecular virological events and induced immunological features during SARS-CoV-2 replication can provide reliable targets for vaccine and drug development. Among the potential targets, subgenomic RNAs and their encoded proteins involved in the life cycle of SARS-CoV-2 are extremely important in viral duplication and pathogenesis. Subgenomic RNAs employ a range of coping strategies to evade immune surveillance from replication to translation, which allows RNAs to synthesize quickly, encode structural proteins efficiently and complete the entire process of virus replication and assembly successfully. This review focuses on the characteristics and functions of SARS-CoV-2 subgenomic RNAs and their encoded proteins and explores in depth the role of subgenomic RNAs in the replication and infection of host cells to provide important clues to the mechanism of COVID-19 pathogenesis.

The roles of APOBEC-mediated RNA editing in SARS-CoV-2 mutations, replication and fitness

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.

On the Origins of Omicron's Unique Spike Gene Insertion

The emergence of a heavily mutated SARS-CoV-2 variant (Omicron; Pango lineage B.1.1.529 and BA sublineages) and its rapid spread to over 75 countries raised a global public health alarm. Characterizing the mutational profile of Omicron is necessary to interpret its clinical phenotypes which are shared with or distinctive from those of other SARS-CoV-2 variants. We compared the mutations of the initially circulating Omicron variant (now known as BA.1) with prior variants of concern (Alpha, Beta, Gamma, and Delta), variants of interest (Lambda, Mu, Eta, Iota, and Kappa), and ~1500 SARS-CoV-2 lineages constituting ~5.8 million SARS-CoV-2 genomes. Omicron's Spike protein harbors 26 amino acid mutations (23 substitutions, 2 deletions, and 1 insertion) that are distinct compared to other variants of concern. While the substitution and deletion mutations appeared in previous SARS-CoV-2 lineages, the insertion mutation (ins214EPE) was not previously observed in any other SARS-CoV-2 lineage. Here, we consider and discuss various mechanisms through which the nucleotide sequence encoding for ins214EPE could have been acquired, including local duplication, polymerase slippage, and template switching. Although we are not able to definitively determine the mechanism, we highlight the plausibility of template switching. Analysis of the homology of the inserted nucleotide sequence and flanking regions suggests that this template-switching event could have involved the genomes of SARS-CoV-2 variants (e.g., the B.1.1 strain), other human coronaviruses that infect the same host cells as SARS-CoV-2 (e.g., HCoV-OC43 or HCoV-229E), or a human transcript expressed in a host cell that was infected by the Omicron precursor.

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.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

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

Background

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

Methods

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

Results

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

Discussion and Conclusion

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

Long-read sequencing reveals complex patterns of wraparound transcription in polyomaviruses

Polyomaviruses (PyV) are ubiquitous pathogens that can cause devastating human diseases. Due to the small size of their genomes, PyV utilize complex patterns of RNA splicing to maximize their coding capacity. Despite the importance of PyV to human disease, their transcriptome architecture is poorly characterized. Here, we compare short- and long-read RNA sequencing data from eight human and non-human PyV. We provide a detailed transcriptome atlas for BK polyomavirus (BKPyV), an important human pathogen, and the prototype PyV, simian virus 40 (SV40). We identify pervasive wraparound transcription in PyV, wherein transcription runs through the polyA site and circles the genome multiple times. Comparative analyses identify novel, conserved transcripts that increase PyV coding capacity. One of these conserved transcripts encodes superT, a T antigen containing two RB-binding LxCxE motifs. We find that superT-encoding transcripts are abundant in PyV-associated human cancers. Together, we show that comparative transcriptomic approaches can greatly expand known transcript and coding capacity in one of the simplest and most well-studied viral families.

Virus-like Particles: Measures and Biological Functions

Virus-like particles resemble infectious virus particles in size, shape, and molecular composition; however, they fail to productively infect host cells. Historically, the presence of virus-like particles has been inferred from total particle counts by microscopy, and infectious particle counts or plaque-forming-units (PFUs) by plaque assay; the resulting ratio of particles-to-PFUs is often greater than one, easily 10 or 100, indicating that most particles are non-infectious. Despite their inability to hijack cells for their reproduction, virus-like particles and the defective genomes they carry can exhibit a broad range of behaviors: interference with normal virus growth during co-infections, cell killing, and activation or inhibition of innate immune signaling. In addition, some virus-like particles become productive as their multiplicities of infection increase, a sign of cooperation between particles. Here, we review established and emerging methods to count virus-like particles and characterize their biological functions. We take a critical look at evidence for defective interfering virus genomes in natural and clinical isolates, and we review their potential as antiviral therapeutics. In short, we highlight an urgent need to better understand how virus-like genomes and particles interact with intact functional viruses during co-infection of their hosts, and their impacts on the transmission, severity, and persistence of virus-associated diseases.

Secondary Structure of Subgenomic RNA M of SARS-CoV-2

SARS-CoV-2 belongs to the Coronavirinae family. Like other coronaviruses, SARS-CoV-2 is enveloped and possesses a positive-sense, single-stranded RNA genome of ~30 kb. Genomic RNA is used as the template for replication and transcription. During these processes, positive-sense genomic RNA (gRNA) and subgenomic RNAs (sgRNAs) are created. Several studies presented the importance of the genomic RNA secondary structure in SARS-CoV-2 replication. However, the structure of sgRNAs has remained largely unsolved so far. In this study, we probed the sgRNA M model of SARS-CoV-2 in vitro. The presented model molecule includes 5'UTR and a coding sequence of gene M. This is the first experimentally informed secondary structure model of sgRNA M, which presents features likely to be important in sgRNA M function. The knowledge of sgRNA M structure provides insights to better understand virus biology and could be used for designing new therapeutics.

Biology of the SARS-CoV-2 Coronavirus

New coronavirus infection causing COVID-19, which was first reported in late 2019 in China, initiated severe social and economic crisis that affected the whole world. High frequency of the errors in replication of RNA viruses, zoonotic nature of transmission, and high transmissibility allowed betacoronaviruses to cause the third pandemic in the world since the beginning of 2003: SARS-CoV in 2003, MERS-CoV in 2012, and SARS-CoV-2 in 2019. The latest pandemic united scientific community and served as a powerful impetus in the study of biology of coronaviruses: new routes of virus penetration into the human cells were identified, features of the replication cycle were studied, and new functions of coronavirus proteins were elucidated. It should be recognized that the pandemic was accompanied by the need to obtain and publish results within a short time, which led to the emergence of an array of conflicting data and low reproducibility of research results. We systematized and analyzed scientific literature, filtered the results according to reliability of the methods of analysis used, and prepared a review describing molecular mechanisms of functioning of the SARS-CoV-2 coronavirus. This review considers organization of the genome of the SARS-CoV-2 virus, mechanisms of its gene expression and entry of the virus into the cell, provides information on key mutations that characterize different variants of the virus, and their contribution to pathogenesis of the disease.

Characterization of two SARS-CoV-2 subgenomic RNA dynamics in severe COVID-19 patients

Little is known about Subgenomic RNA (sgRNA) dynamics in patients with Coronavirus diseases 2019 (COVID-19). We collected 147 throat swabs, 74 gut swabs and 46 plasma samples from 117 COVID-19 patients recruited in the LOTUS China trial (ChiCTR2000029308) and compared E and orf7a sgRNA load in patients with different illness duration, outcome, and comorbidities. Both sgRNAs were detected in all the three types of samples, with longest duration of 25, 13, and 17 days for E sgRNA, and 32, 28, and 17 days for orf7a sgRNA in throat, gut, and plasma, respectively. A total of 95% (57/60) of patients had no E sgRNA detected after 10 days post treatment, though 86% of them were still E RNA positive. High correlation on titer was observed between sgRNA encoding E and orf7a gene. sgRNA showed similar variation in the standard care and Lopinavir-Ritonavir group. Patients with diabetes and heart diseases showed higher pharyngeal E sgRNA at the first day (P ​= ​0.016 and 0.013, respectively) but no difference at five days after treatment, compared with patients without such commodities. Patients with hypertension and cerebrovascular diseases showed no difference in the pharyngeal sgRNA levels at both one and five days after treatment, compared with patients without these two commodities. E sgRNA levels in the initial infection showed no correlation with the serum antibody against spike, nucleoprotein, and receptor binding domains at ten days later. sgRNA lasted a long period in COVID-19 patients and might have little effect on humoral response.

•

147 throat swabs, 74 gut swabs and 46 plasma samples were used for subgenomic RNA (sgRNA) E and orf7a detection and quantification.

•

sgRNA of E and orf7a gene were detected in throat swabs, gut swabs and plasmas and lasted long period.

•

E sgRNA levels in the initial infection showed no correlation with the serum antibody against SARS-CoV-2 ten days later.

The Roles of APOBEC-mediated RNA Editing in SARS-CoV-2 Mutations, Replication and Fitness.. [PMID: 34981048 PMCID: PMC8722585 DOI: 10.1101/2021.12.18.473309]

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.

Efficient Editing of SARS-CoV-2 genomic RNA by Host APOBEC deaminases and Its Potential Impacts on the Viral Replication and Emergence of New Strains in COVID-19 Pandemic

Template switching and duplications in SARS-CoV-2 genomes give rise to insertion variants that merit monitoring

The appearance of multiple new SARS-CoV-2 variants during the COVID-19 pandemic is a matter of grave concern. Some of these variants, such as B.1.617.2, B.1.1.7, and B.1.351, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the Spike (S) protein is thought to be a determinant of SARS-CoV-2 virulence. Here, we identify 346 unique inserts of different lengths in SARS-CoV-2 genomes and present evidence that these inserts reflect actual virus variance rather than sequencing artifacts. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. At least three inserts in the N-terminal domain of the S protein are predicted to lead to escape from neutralizing antibodies, whereas other inserts might result in escape from T-cell immunity. Thus, inserts in the S protein can affect its antigenic properties and merit monitoring.

CoV-er all the bases: Structural perspectives of SARS-CoV-2 RNA synthesis

The ongoing Covid-19 pandemic has spurred research in the biology of the nidovirus severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Much focus has been on the viral RNA synthesis machinery due to its fundamental role in viral propagation. The central and essential enzyme of the RNA synthesis process, the RNA-dependent RNA polymerase (RdRp), functions in conjunction with a coterie of viral-encoded enzymes that mediate crucial nucleic acid transactions. Some of these enzymes share common features with other RNA viruses, while others play roles unique to nidoviruses or CoVs. The RdRps are proven targets for viral pathogens, and many of the other nucleic acid processing enzymes are promising targets. The purpose of this review is to summarize recent advances in our understanding of the mechanisms of RNA synthesis in CoVs. By reflecting on these studies, we hope to emphasize the remaining gaps in our knowledge. The recent onslaught of structural information related to SARS-CoV-2 RNA synthesis, in combination with previous structural, genetic and biochemical studies, have vastly improved our understanding of how CoVs replicate and process their genomic RNA. Structural biology not only provides a blueprint for understanding the function of the enzymes and cofactors in molecular detail, but also provides a basis for drug design and optimization. The concerted efforts of researchers around the world, in combination with the renewed urgency toward understanding this deadly family of viruses, may eventually yield new and improved antivirals that provide relief to the current global devastation.

In vivo structure and dynamics of the SARS-CoV-2 RNA genome

The dynamics of SARS-CoV-2 RNA structure and their functional relevance are largely unknown. Here we develop a simplified SPLASH assay and comprehensively map the in vivo RNA-RNA interactome of SARS-CoV-2 genome across viral life cycle. We report canonical and alternative structures including 5'-UTR and 3'-UTR, frameshifting element (FSE) pseudoknot and genome cyclization in both cells and virions. We provide direct evidence of interactions between Transcription Regulating Sequences, which facilitate discontinuous transcription. In addition, we reveal alternative short and long distance arches around FSE. More importantly, we find that within virions, while SARS-CoV-2 genome RNA undergoes intensive compaction, genome domains remain stable but with strengthened demarcation of local domains and weakened global cyclization. Taken together, our analysis reveals the structural basis for the regulation of replication, discontinuous transcription and translational frameshifting, the alternative conformations and the maintenance of global genome organization during the whole life cycle of SARS-CoV-2, which we anticipate will help develop better antiviral strategies.

SARS-CoV-2 Subgenomic RNAs: Characterization, Utility, and Perspectives

SARS-CoV-2, the etiologic agent at the root of the ongoing COVID-19 pandemic, harbors a large RNA genome from which a tiered ensemble of subgenomic RNAs (sgRNAs) is generated. Comprehensive definition and investigation of these RNA products are important for understanding SARS-CoV-2 pathogenesis. This review summarizes the recent progress on SARS-CoV-2 sgRNA identification, characterization, and application as a viral replication marker. The significance of these findings and potential future research areas of interest are discussed.

Reduced subgenomic RNA expression is a molecular indicator of asymptomatic SARS-CoV-2 infection

Background

It is estimated that up to 80% of infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are asymptomatic and asymptomatic patients can still effectively transmit the virus and cause disease. While much of the effort has been placed on decoding single nucleotide variation in SARS-CoV-2 genomes, considerably less is known about their transcript variation and any correlation with clinical severity in human hosts, as defined here by the presence or absence of symptoms.

Methods

To assess viral genomic signatures of disease severity, we conducted a systematic characterization of SARS-CoV-2 transcripts and genetic variants in 81 clinical specimens collected from symptomatic and asymptomatic individuals using multi-scale transcriptomic analyses including amplicon-seq, short-read metatranscriptome and long-read Iso-seq.

Results

Here we show a highly coordinated and consistent pattern of sgRNA expression from individuals with robust SARS-CoV-2 symptomatic infection and their expression is significantly repressed in the asymptomatic infections. We also observe widespread inter- and intra-patient variants in viral RNAs, known as quasispecies frequently found in many RNA viruses. We identify unique sets of deletions preferentially found primarily in symptomatic individuals, with many likely to confer changes in SARS-CoV-2 virulence and host responses. Moreover, these frequently occurring structural variants in SARS-CoV-2 genomes serve as a mechanism to further induce SARS-CoV-2 proteome complexity.

Conclusions

Our results indicate that differential sgRNA expression and structural mutational burden are highly correlated with the clinical severity of SARS-CoV-2 infection. Longitudinally monitoring sgRNA expression and structural diversity could further guide treatment responses, testing strategies, and vaccine development.

Age and Sex Modulate SARS-CoV-2 Viral Load Kinetics: A Longitudinal Analysis of 1735 Subjects

The COVID-19 pandemic caused by SARS-CoV-2 represents a public health emergency, which became even more challenging since the detection of highly transmissible variants and strategies against COVID-19 were indistinctly established. We characterized the temporal viral load kinetics in individuals infected by original and variant strains. Naso-oropharyngeal swabs from 33,000 individuals (admitted to the IRCCS Santa Lucia Foundation Drive-in, healthcare professionals and hospitalized patients who underwent routinary screening) from November 2020 to June 2021 were analyzed. Of them, 1735 subjects were selected and grouped according to the viral strain. Diagnostic analyses were performed by CE-IVD RT-PCR-based kits. The subgenomic-RNA component was assessed in 36 subjects using digital PCR. Infection duration, viral load decay speed, effects of age and sex were assessed and compared by extensive statistical analyses. Overall, infection duration and viral load differed between the groups (p < 0.05). Male sex was more present among both original and variant carriers affected with high viral load and showing fast decay speed, whereas original strain carriers with slow decay speed resulted in older (p < 0.05). Subgenomic-RNA was detected in the positive samples, including those with low viral load. This study provides a picture of the viral load kinetics, identifying individuals with similar patterns and showing differential effects of age and sex, thus providing potentially useful information for personalized management of infected subjects.

Insertions in SARS-CoV-2 genome caused by template switch and duplications give rise to new variants that merit monitoring

The appearance of multiple new SARS-CoV-2 variants during the winter of 2020-2021 is a matter of grave concern. Some of these new variants, such as B.1.617.2, B.1.1.7, and B.1.351, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the COVID-19 pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on point nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the spike protein is thought to be a determinant of SARS-CoV-2 virulence and other inserts might have contributed to coronavirus pathogenicity as well. Here, we investigate insertions in SARS-CoV-2 genomes and identify 347 unique inserts of different lengths. We present evidence that these inserts reflect actual virus variance rather than sequencing errors. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. We show that inserts in the Spike glycoprotein can affect its antigenic properties and thus merit monitoring. At least, three inserts in the N-terminal domain of the Spike (ins245IME, ins246DSWG, and ins248SSLT) that were first detected in 2021 are predicted to lead to escape from neutralizing antibodies, whereas other inserts might result in escape from T-cell immunity.

Impact of Full Vaccination with mRNA BNT162b2 on SARS-CoV-2 Infection: Genomic and Subgenomic Viral RNAs Detection in Nasopharyngeal Swab and Saliva of Health Care Workers

SARS-CoV-2 infection was monitored in 1898 health care workers (HCWs) after receiving full vaccination with BNT162b2. Untill 30 June 2021, 10 HCWs tested positive for SARS-CoV-2 using real time RT-PCR, resulting in a 4-month cumulative incidence of 0.005%. The infection was mildly symptomatic in six (60%) and asymptomatic in four (40%) individuals. Among the infected HCWs, eight consenting individuals provided paired NPS and saliva during the course of infection, for the purpose of the analysis performed in the present study. Genomic and subgenomic viral RNAs were investigated using real-time RT-PCR in both biological specimens. The temporal profile of viral load was measured using ddPCR. Viral mutations were also analysed. Subgenomic viral RNA was detected in 8/8 (100%) NPS and in 6/8 (75%) saliva specimens at the baseline. The expression of subgenomic RNA was observed for up to 7 days in 3/8 (38%) symptomatic cases. Moreover, concordance was observed between NPS and saliva in the detection of viral mutations, and both N501Y and 69/70del (associated with the B.1.1.7 variant) were detected in the majority 6/8 (75%) of subjects, while the K417T mutation (associated with the P.1-type variants) was detected in 2/8 (25%) individuals. Overall, our findings report a low frequency of infected HCWs after full vaccination. It is, therefore, important to monitor the vaccinees in order to identify asymptomatic infected individuals. Saliva can be a surrogate for SARS-CoV-2 surveillance, particularly in social settings such as hospitals.

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.

Identification of presented SARS-CoV-2 HLA class I and HLA class II peptides using HLA peptidomics

The human leukocyte antigen (HLA)-bound viral antigens serve as an immunological signature that can be selectively recognized by T cells. As viruses evolve by acquiring mutations, it is essential to identify a range of presented viral antigens. Using HLA peptidomics, we are able to identify severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-derived peptides presented by highly prevalent HLA class I (HLA-I) molecules by using infected cells as well as overexpression of SARS-CoV-2 genes. We find 26 HLA-I peptides and 36 HLA class II (HLA-II) peptides. Among the identified peptides, some are shared between different cells and some are derived from out-of-frame open reading frames (ORFs). Seven of these peptides were previously shown to be immunogenic, and we identify two additional immunoreactive peptides by using HLA multimer staining. These results may aid the development of the next generation of SARS-CoV-2 vaccines based on presented viral-specific antigens that span several of the viral genes.

A colloidal gold-based immunochromatographic strip for rapid detection of SARS-CoV-2 antibodies after vaccination

Vaccination interventions is consideredan important preventive measure to block the transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and protect the organism from pathogen infection effectively. However, a quick and accurate technique to evaluate the immune efficacy of the SARS-CoV-2 inactivated vaccine remains scarce. In this paper, an IgM-IgG antibody combined detection colloidal gold immunochromatography assay kit was optimized and developed, which can assess the efficacy of the inactivated SARS-CoV-2 vaccine. We collected fingertip blood samples from 3 vaccinees and 1 unvaccinated sample. The results showed that the proportion of antibody was high after the second shots immunization. The colloidal gold-based immunochromatographic strip is rapid, convenient and easy to operate. It can be used as an auxiliary method for preliminary evaluation of the antibody effect of vaccine recipients, and provide a reference index for the potential clinical application value of the vaccine.

A Novel SARS-CoV-2 Viral Sequence Bioinformatic Pipeline Has Found Genetic Evidence That the Viral 3' Untranslated Region (UTR) Is Evolving and Generating Increased Viral Diversity

An unprecedented amount of SARS-CoV-2 sequencing has been performed, however, novel bioinformatic tools to cope with and process these large datasets is needed. Here, we have devised a bioinformatic pipeline that inputs SARS-CoV-2 genome sequencing in FASTA/FASTQ format and outputs a single Variant Calling Format file that can be processed to obtain variant annotations and perform downstream population genetic testing. As proof of concept, we have analyzed over 229,000 SARS-CoV-2 viral sequences up until November 30, 2020. We have identified over 39,000 variants worldwide with increased polymorphisms, spanning the ORF3a gene as well as the 3' untranslated (UTR) regions, specifically in the conserved stem loop region of SARS-CoV-2 which is accumulating greater observed viral diversity relative to chance variation. Our analysis pipeline has also discovered the existence of SARS-CoV-2 hypermutation with low frequency (less than in 2% of genomes) likely arising through host immune responses and not due to sequencing errors. Among annotated non-sense variants with a population frequency over 1%, recurrent inactivation of the ORF8 gene was found. This was found to be present in the newly identified B.1.1.7 SARS-CoV-2 lineage that originated in the United Kingdom. Almost all VOC-containing genomes possess one stop codon in ORF8 gene (Q27∗), however, 13% of these genomes also contains another stop codon (K68∗), suggesting that ORF8 loss does not interfere with SARS-CoV-2 spread and may play a role in its increased virulence. We have developed this computational pipeline to assist researchers in the rapid analysis and characterization of SARS-CoV-2 variation.

Proteomics-Based Insights Into the SARS-CoV-2-Mediated COVID-19 Pandemic: A Review of the First Year of Research

In late 2019, a virus subsequently named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in China and led to a worldwide pandemic of the disease termed coronavirus disease 2019. The global health threat posed by this pandemic led to an extremely rapid and robust mobilization of the scientific and medical communities as evidenced by the publication of more than 10,000 peer-reviewed articles and thousands of preprints in the first year of the pandemic alone. With the publication of the initial genome sequence of SARS-CoV-2, the proteomics community immediately joined this effort publishing, to date, more than 100 peer-reviewed proteomics studies and submitting many more preprints to preprint servers. In this review, we focus on peer-reviewed articles published on the proteome, glycoproteome, and glycome of SARS-CoV-2. At a basic level, proteomic studies provide valuable information on quantitative aspects of viral infection course; information on the identities, sites, and microheterogeneity of post-translational modifications; and, information on protein-protein interactions. At a biological systems level, these studies elucidate host cell and tissue responses, characterize antibodies and other immune system factors in infection, suggest biomarkers that may be useful for diagnosis and disease-course monitoring, and help in the development or repurposing of potential therapeutics. Here, we summarize results from selected early studies to provide a perspective on the current rapidly evolving literature.

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.

Transcriptional and epi-transcriptional dynamics of SARS-CoV-2 during cellular infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses subgenomic RNA (sgRNA) to produce viral proteins for replication and immune evasion. We apply long-read RNA and cDNA sequencing to in vitro human and primate infection models to study transcriptional dynamics. Transcription-regulating sequence (TRS)-dependent sgRNA upregulates earlier in infection than TRS-independent sgRNA. An abundant class of TRS-independent sgRNA consisting of a portion of open reading frame 1ab (ORF1ab) containing nsp1 joins to ORF10, and the 3' untranslated region (UTR) upregulates at 48 h post-infection in human cell lines. We identify double-junction sgRNA containing both TRS-dependent and -independent junctions. We find multiple sites at which the SARS-CoV-2 genome is consistently more modified than sgRNA and that sgRNA modifications are stable across transcript clusters, host cells, and time since infection. Our work highlights the dynamic nature of the SARS-CoV-2 transcriptome during its replication cycle.

Subgenomic RNA identification in SARS-CoV-2 genomic sequencing data

We have developed periscope, a tool for the detection and quantification of subgenomic RNA (sgRNA) in SARS-CoV-2 genomic sequence data. The translation of the SARS-CoV-2 RNA genome for most open reading frames (ORFs) occurs via RNA intermediates termed "subgenomic RNAs." sgRNAs are produced through discontinuous transcription, which relies on homology between transcription regulatory sequences (TRS-B) upstream of the ORF start codons and that of the TRS-L, which is located in the 5' UTR. TRS-L is immediately preceded by a leader sequence. This leader sequence is therefore found at the 5' end of all sgRNA. We applied periscope to 1155 SARS-CoV-2 genomes from Sheffield, United Kingdom, and validated our findings using orthogonal data sets and in vitro cell systems. By using a simple local alignment to detect reads that contain the leader sequence, we were able to identify and quantify reads arising from canonical and noncanonical sgRNA. We were able to detect all canonical sgRNAs at the expected abundances, with the exception of ORF10. A number of recurrent noncanonical sgRNAs are detected. We show that the results are reproducible using technical replicates and determine the optimum number of reads for sgRNA analysis. In VeroE6 ACE2+/- cell lines, periscope can detect the changes in the kinetics of sgRNA in orthogonal sequencing data sets. Finally, variants found in genomic RNA are transmitted to sgRNAs with high fidelity in most cases. This tool can be applied to all sequenced COVID-19 samples worldwide to provide comprehensive analysis of SARS-CoV-2 sgRNA.