Gene N proximal and distal RNA motifs regulate coronavirus nucleocapsid mRNA transcription

m6A Methylation of Transcription Leader Sequence of SARS-CoV-2 Impacts Discontinuous Transcription of Subgenomic mRNAs

The SARS-CoV-2 genome has been shown to be m6A methylated at several positions in vivo. Strikingly, a DRACH motif, the recognition motif for adenosine methylation, resides in the core of the transcriptional regulatory leader sequence (TRS-L) at position A74, which is highly conserved and essential for viral discontinuous transcription. Methylation at position A74 correlates with viral pathogenicity. Discontinuous transcription produces a set of subgenomic mRNAs that function as templates for translation of all structural and accessory proteins. A74 is base-paired in the short stem-loop structure 5'SL3 that opens during discontinuous transcription to form long-range RNA-RNA interactions with nascent (-)-strand transcripts at complementary TRS-body sequences. A74 can be methylated by the human METTL3/METTL14 complex in vitro. Here, we investigate its impact on the structural stability of 5'SL3 and the long-range TRS-leader:TRS-body duplex formation necessary for synthesis of subgenomic mRNAs of all four viral structural proteins. Methylation uniformly destabilizes 5'SL3 and long-range duplexes and alters their relative equilibrium populations, suggesting that the m6A74 modification acts as a regulator for the abundance of viral structural proteins due to this destabilization.

Classification, replication, and transcription of Nidovirales

Nidovirales is one order of RNA virus, with the largest single-stranded positive sense RNA genome enwrapped with membrane envelope. It comprises four families (Arterividae, Mesoniviridae, Roniviridae, and Coronaviridae) and has been circulating in humans and animals for almost one century, posing great threat to livestock and poultry，as well as to public health. Nidovirales shares similar life cycle: attachment to cell surface, entry, primary translation of replicases, viral RNA replication in cytoplasm, translation of viral proteins, virion assembly, budding, and release. The viral RNA synthesis is the critical step during infection, including genomic RNA (gRNA) replication and subgenomic mRNAs (sg mRNAs) transcription. gRNA replication requires the synthesis of a negative sense full-length RNA intermediate, while the sg mRNAs transcription involves the synthesis of a nested set of negative sense subgenomic intermediates by a discontinuous strategy. This RNA synthesis process is mediated by the viral replication/transcription complex (RTC), which consists of several enzymatic replicases derived from the polyprotein 1a and polyprotein 1ab and several cellular proteins. These replicases and host factors represent the optimal potential therapeutic targets. Hereby, we summarize the Nidovirales classification, associated diseases, "replication organelle," replication and transcription mechanisms, as well as related regulatory factors.

Molecular mechanisms of human coronavirus NL63 infection and replication

Human coronavirus NL63 (HCoV-NL63) is spread globally, causing upper and lower respiratory tract infections mainly in young children. HCoV-NL63 shares a host receptor (ACE2) with severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 but, unlike them, HCoV-NL63 primarily develops into self-limiting mild to moderate respiratory disease. Although with different efficiency, both HCoV-NL63 and SARS-like CoVs infect ciliated respiratory cells using ACE2 as receptor for binding and cell entry. Working with SARS-like CoVs require access to BSL-3 facilities, while HCoV-NL63 research can be performed at BSL-2 laboratories. Thus, HCoV-NL63 could be used as a safer surrogate for comparative studies on receptor dynamics, infectivity and virus replication, disease mechanism, and potential therapeutic interventions against SARS-like CoVs. This prompted us to review the current knowledge on the infection mechanism and replication of HCoV-NL63. Specifically, after a brief overview on the taxonomy, genomic organization and virus structure, this review compiles the current HCoV-NL63-related research in virus entry and replication mechanism, including virus attachment, endocytosis, genome translation, and replication and transcription. Furthermore, we reviewed cumulative knowledge on the susceptibility of different cells to HCoV-NL63 infection in vitro, which is essential for successful virus isolation and propagation, and contribute to address different scientific questions from basic science to the development and assessment of diagnostic tools, and antiviral therapies. Finally, we discussed different antiviral strategies that have been explored to suppress replication of HCoV-NL63, and other related human coronaviruses, by either targeting the virus or enhancing host antiviral mechanisms.

Variant-Specific Analysis Reveals a Novel Long-Range RNA-RNA Interaction in SARS-CoV-2 Orf1a

Since the start of the COVID-19 pandemic, understanding the pathology of the SARS-CoV-2 RNA virus and its life cycle has been the priority of many researchers. Currently, new variants of the virus have emerged with various levels of pathogenicity and abundance within the human-host population. Although much of viral pathogenicity is attributed to the viral Spike protein’s binding affinity to human lung cells’ ACE2 receptor, comprehensive knowledge on the distinctive features of viral variants that might affect their life cycle and pathogenicity is yet to be attained. Recent in vivo studies into the RNA structure of the SARS-CoV-2 genome have revealed certain long-range RNA-RNA interactions. Using in silico predictions and a large population of SARS-CoV-2 sequences, we observed variant-specific evolutionary changes for certain long-range RRIs. We also found statistical evidence for the existence of one of the thermodynamic-based RRI predictions, namely Comp1, in the Beta variant sequences. A similar test that disregarded sequence variant information did not, however, lead to significant results. When performing population-based analyses, aggregate tests may fail to identify novel interactions due to variant-specific changes. Variant-specific analyses can result in de novo RRI identification.

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

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

Thermal inactivation of COVID-19 specimens improves RNA quality and quantity

The rapid spread of coronavirus disease 2019 (COVID-19), a disease caused by severe acute respiratory syndrome coronavirus 2, poses a huge demand for immediate diagnosis. Real-time reverse transcriptase-polymerase chain reaction (rRT-PCR) of nasopharyngeal (NP) and oropharyngeal (OP) swabs have been used to confirm the clinical diagnosis. To avoid the risk of viral-exposure of laboratory workers, thermal inactivation is currently recommended but has unknown effects on the accuracy of the rRT-PCR results. Thirty-six NP/OP specimens were collected from COVID-19 patients and subjected to thermal inactivation (60°C for 30 min) or the RNA extraction processes to activate the form. Here, our data showed that the concentration of extracted-RNA increases upon thermal inactivation compared to the active form (p = .028).  Significantly higher levels of RNA copy number were obtained in inactivated compared to the active samples for both N and ORF1ab genes (p = .009, p = .032, respectively). Thermal inactivation elevated concentration and copy number of extracted-RNA, possibly through viral-capsid degradation and/or nucleoprotein denaturation.

Activation of viral transcription by stepwise largescale folding of an RNA virus genome

The genomes of RNA viruses contain regulatory elements of varying complexity. Many plus-strand RNA viruses employ largescale intra-genomic RNA-RNA interactions as a means to control viral processes. Here, we describe an elaborate RNA structure formed by multiple distant regions in a tombusvirus genome that activates transcription of a viral subgenomic mRNA. The initial step in assembly of this intramolecular RNA complex involves the folding of a large viral RNA domain, which generates a discontinuous binding pocket. Next, a distally-located protracted stem-loop RNA structure docks, via base-pairing, into the binding site and acts as a linchpin that stabilizes the RNA complex and activates transcription. A multi-step RNA folding pathway is proposed in which rate-limiting steps contribute to a delay in transcription of the capsid protein-encoding viral subgenomic mRNA. This study provides an exceptional example of the complexity of genome-scale viral regulation and offers new insights into the assembly schemes utilized by large intra-genomic RNA structures.

Evaluation of a recombination-resistant coronavirus as a broadly applicable, rapidly implementable vaccine platform

Emerging and re-emerging zoonotic viral diseases are major threats to global health, economic stability, and national security. Vaccines are key for reducing coronaviral disease burden; however, the utility of live-attenuated vaccines is limited by risks of reversion or repair. Because of their history of emergence events due to their prevalence in zoonotic pools, designing live-attenuated coronavirus vaccines that can be rapidly and broadly implemented is essential for outbreak preparedness. Here, we show that coronaviruses with completely rewired transcription regulatory networks (TRNs) are effective vaccines against SARS-CoV. The TRN-rewired viruses are attenuated and protect against lethal SARS-CoV challenge. While a 3-nt rewired TRN reverts via second-site mutation upon serial passage, a 7-nt rewired TRN is more stable, suggesting that a more extensively rewired TRN might be essential for avoiding growth selection. In summary, rewiring the TRN is a feasible strategy for limiting reversion in an effective live-attenuated coronavirus vaccine candidate that is potentially portable across the Nidovirales order.

Continuous and Discontinuous RNA Synthesis in Coronaviruses

Replication of the coronavirus genome requires continuous RNA synthesis, whereas transcription is a discontinuous process unique among RNA viruses. Transcription includes a template switch during the synthesis of subgenomic negative-strand RNAs to add a copy of the leader sequence. Coronavirus transcription is regulated by multiple factors, including the extent of base-pairing between transcription-regulating sequences of positive and negative polarity, viral and cell protein-RNA binding, and high-order RNA-RNA interactions. Coronavirus RNA synthesis is performed by a replication-transcription complex that includes viral and cell proteins that recognize cis-acting RNA elements mainly located in the highly structured 5' and 3' untranslated regions. In addition to many viral nonstructural proteins, the presence of cell nuclear proteins and the viral nucleocapsid protein increases virus amplification efficacy. Coronavirus RNA synthesis is connected with the formation of double-membrane vesicles and convoluted membranes. Coronaviruses encode proofreading machinery, unique in the RNA virus world, to ensure the maintenance of their large genome size.

Nucleocapsid phosphorylation and RNA helicase DDX1 recruitment enables coronavirus transition from discontinuous to continuous transcription

Coronaviruses contain a positive-sense single-stranded genomic (g) RNA, which encodes nonstructural proteins. Several subgenomic mRNAs (sgmRNAs) encoding structural proteins are generated by template switching from the body transcription regulatory sequence (TRS) to the leader TRS. The process preferentially generates shorter sgmRNA. Appropriate readthrough of body TRSs is required to produce longer sgmRNAs and full-length gRNA. We find that phosphorylation of the viral nucleocapsid (N) by host glycogen synthase kinase-3 (GSK-3) is required for template switching. GSK-3 inhibition selectively reduces the generation of gRNA and longer sgmRNAs, but not shorter sgmRNAs. N phosphorylation allows recruitment of the RNA helicase DDX1 to the phosphorylated-N-containing complex, which facilitates template readthrough and enables longer sgmRNA synthesis. DDX1 knockdown or loss of helicase activity markedly reduces the levels of longer sgmRNAs. Thus, coronaviruses employ a unique strategy for the transition from discontinuous to continuous transcription to ensure balanced sgmRNAs and full-length gRNA synthesis.

Transmissible gastroenteritis virus (TGEV) infection alters the expression of cellular microRNA species that affect transcription of TGEV gene 7

Transmissible gastroenteritis virus (TGEV) is a member of Coronaviridae family. TGEV infection has emerged as a major cause of severe gastroenteritis and leads to alterations of many cellular processes. Meanwhile, the pathogenic mechanism of TGEV is still unclear. microRNAs (miRNAs) are a novel class of small non-coding RNAs which are involved in the regulation of numerous biological processes such as viral infection and cell apoptosis. Accumulating data show that miRNAs are involved in the process of coronavirus infection such as replication of severe acute respiratory syndrome coronavirus (SARS-CoV). However, the link between miRNAs and TGEV infection is unknown. In this study, we performed microRNA microarray assay and predicted targets of altered miRNAs. The results showed TGEV infection caused the change of miRNAs profile. Then we selected miR-4331 for further analysis and subsequently identified cell division cycle-associated protein 7 (CDCA7) as the target of miR-4331. Moreover, miR-4331 showed the ability to inhibit transcription of TGEV gene 7 (a non-structure gene) via directly targeting CDCA7. In conclusion, differentially expressed miR-4331 that is caused by TGEV infection can suppress transcription of TGEV gene 7 via targeting cellular CDCA7. Our key finding is that TGEV selectively manipulates the expression of some cellular miRNAs to regulate its subgenomic transcription.

The nucleocapsid protein of human coronavirus NL63

Human coronavirus (HCoV) NL63 was first described in 2004 and is associated with respiratory tract disease of varying severity. At the genetic and structural level, HCoV-NL63 is similar to other members of the Coronavirinae subfamily, especially human coronavirus 229E (HCoV-229E). Detailed analysis, however, reveals several unique features of the pathogen. The coronaviral nucleocapsid protein is abundantly present in infected cells. It is a multi-domain, multi-functional protein important for viral replication and a number of cellular processes. The aim of the present study was to characterize the HCoV-NL63 nucleocapsid protein. Biochemical analyses revealed that the protein shares characteristics with homologous proteins encoded in other coronaviral genomes, with the N-terminal domain responsible for nucleic acid binding and the C-terminal domain involved in protein oligomerization. Surprisingly, analysis of the subcellular localization of the N protein of HCoV-NL63 revealed that, differently than homologous proteins from other coronaviral species except for SARS-CoV, it is not present in the nucleus of infected or transfected cells. Furthermore, no significant alteration in cell cycle progression in cells expressing the protein was observed. This is in stark contrast with results obtained for other coronaviruses, except for the SARS-CoV.

Antigenic structures stably expressed by recombinant TGEV-derived vectors

Coronaviruses (CoVs) are positive-stranded RNA viruses with potential as immunization vectors, expressing high levels of heterologous genes and eliciting both secretory and systemic immune responses. Nevertheless, its high recombination rate may result in the loss of the full-length foreign gene, limiting their use as vectors. Transmissible gastroenteritis virus (TGEV) was engineered to express porcine reproductive and respiratory syndrome virus (PRRSV) small protein domains, as a strategy to improve heterologous gene stability. After serial passage in tissue cultures, stable expression of small PRRSV protein antigenic domains was achieved. Therefore, size reduction of the heterologous genes inserted in CoV-derived vectors led to the stable expression of antigenic domains. Immunization of piglets with these TGEV vectors led to partial protection against a challenge with a virulent PRRSV strain, as immunized animals showed reduced clinical signs and lung damage. Further improvement of TGEV-derived vectors will require the engineering of vectors with decreased recombination rate.

Functional long-range RNA-RNA interactions in positive-strand RNA viruses

Long-range RNA–RNA interactions, many of which span several thousands of nucleotides, have been discovered within the genomes of positive-strand RNA viruses. These interactions mediate fundamental viral processes, including translation, replication and transcription.

In certain plant viruses that have uncapped, non-polyadenylated RNA genomes, translation initiation is facilitated by 3′ cap-independent translational enhancers (3′ CITEs) that are located in or near to their 3′ UTRs. These RNA elements function by binding to either the ribosome-recruiting eukaryotic translation initiation factor 4F (eIF4F) complex or ribosomal subunits, and they enhance translation initiation by engaging the 5′ end of the genome via a 5′-to-3′ RNA-based bridge.

The activities of the internal ribosome entry sites (IRESs) in the 5′ UTRs of various viruses are modulated by RNA-based interactions between the IRESs and elements near to the 3′ ends of their genomes.

In several plant viruses, translational recoding events, including ribosomal frameshifting and stop codon readthrough, have been found to rely on long-range RNA–RNA interactions.

Multiple 5′-to-3′ base-pairing interactions facilitate genome circularization in flaviviruses, which has been proposed to reposition the 5′-bound RNA-dependent RNA polymerase (RdRp) to the initiation site of negative-strand synthesis at the 3′ terminus.

The long-distance interaction between two cis-acting replication elements in tombusviruses generates a bipartite RNA platform for the assembly of the replicase complex and repositions the internally bound RdRp to the 3′ terminus.

Tombusviruses also rely on several long-range interactions that mediate the premature termination of the RdRp during negative-strand synthesis that leads to transcription of subgenomic mRNAs (sgmRNAs).

In a coronavirus, an exceptionally long-range interaction, which spans ∼26,000 nucleotides, promotes polymerase repriming during the discontinuous template synthesis step of sgmRNA-N transcription.

A challenge for the future will be to determine how these long-range interactions are integrated and regulated in the complex context of viral RNA genomes.

Long-range intragenomic RNA–RNA interactions in the genomes of positive-strand RNA viruses involve direct nucleotide base pairing and can span distances of thousands of nucleotides. In this Review, Nicholson and White discuss recent insights into the structure and function of these genomic features and highlight their diverse roles in the gene expression and genome replication of positive-strand RNA viruses.

Positive-strand RNA viruses are important human, animal and plant pathogens that are defined by their single-stranded positive-sense RNA genomes. In recent years, it has become increasingly evident that interactions that occur between distantly positioned RNA sequences within these genomes can mediate important viral activities. These long-range intragenomic RNA–RNA interactions involve direct nucleotide base pairing and can span distances of thousands of nucleotides. In this Review, we discuss recent insights into the structure and function of these intriguing genomic features and highlight their diverse roles in the gene expression and genome replication of positive-strand RNA viruses.

Transmissible gastroenteritis coronavirus genome packaging signal is located at the 5' end of the genome and promotes viral RNA incorporation into virions in a replication-independent process

Preferential RNA packaging in coronaviruses involves the recognition of viral genomic RNA, a crucial process for viral particle morphogenesis mediated by RNA-specific sequences, known as packaging signals. An essential packaging signal component of transmissible gastroenteritis coronavirus (TGEV) has been further delimited to the first 598 nucleotides (nt) from the 5' end of its RNA genome, by using recombinant viruses transcribing subgenomic mRNA that included potential packaging signals. The integrity of the entire sequence domain was necessary because deletion of any of the five structural motifs defined within this region abrogated specific packaging of this viral RNA. One of these RNA motifs was the stem-loop SL5, a highly conserved motif in coronaviruses located at nucleotide positions 106 to 136. Partial deletion or point mutations within this motif also abrogated packaging. Using TGEV-derived defective minigenomes replicated in trans by a helper virus, we have shown that TGEV RNA packaging is a replication-independent process. Furthermore, the last 494 nt of the genomic 3' end were not essential for packaging, although this region increased packaging efficiency. TGEV RNA sequences identified as necessary for viral genome packaging were not sufficient to direct packaging of a heterologous sequence derived from the green fluorescent protein gene. These results indicated that TGEV genome packaging is a complex process involving many factors in addition to the identified RNA packaging signal. The identification of well-defined RNA motifs within the TGEV RNA genome that are essential for packaging will be useful for designing packaging-deficient biosafe coronavirus-derived vectors and providing new targets for antiviral therapies.

Global organization of a positive-strand RNA virus genome

The genomes of plus-strand RNA viruses contain many regulatory sequences and structures that direct different viral processes. The traditional view of these RNA elements are as local structures present in non-coding regions. However, this view is changing due to the discovery of regulatory elements in coding regions and functional long-range intra-genomic base pairing interactions. The ∼4.8 kb long RNA genome of the tombusvirus tomato bushy stunt virus (TBSV) contains these types of structural features, including six different functional long-distance interactions. We hypothesized that to achieve these multiple interactions this viral genome must utilize a large-scale organizational strategy and, accordingly, we sought to assess the global conformation of the entire TBSV genome. Atomic force micrographs of the genome indicated a mostly condensed structure composed of interconnected protrusions extending from a central hub. This configuration was consistent with the genomic secondary structure model generated using high-throughput selective 2′-hydroxyl acylation analysed by primer extension (i.e. SHAPE), which predicted different sized RNA domains originating from a central region. Known RNA elements were identified in both domain and inter-domain regions, and novel structural features were predicted and functionally confirmed. Interestingly, only two of the six long-range interactions known to form were present in the structural model. However, for those interactions that did not form, complementary partner sequences were positioned relatively close to each other in the structure, suggesting that the secondary structure level of viral genome structure could provide a basic scaffold for the formation of different long-range interactions. The higher-order structural model for the TBSV RNA genome provides a snapshot of the complex framework that allows multiple functional components to operate in concert within a confined context.

The genomes of many important pathogenic viruses are made of RNA. These genomes encode viral proteins and contain regulatory sequences and structures. In some viruses, distant regions of the RNA genome can interact with each other via base pairing, which suggests that certain genomes may take on well-defined conformations. This concept was investigated using a tombusvirus RNA genome that contains several long-range RNA interactions. The results of microscopic and biochemical analyses indicated a compact genome conformation with structured regions radiating from a central core. The structural model was compatible with some, but not all, long-range interactions, suggesting that the genome is a dynamic molecule that assumes different conformations. The analysis also revealed new structural features of the genome, some of which were shown to be functionally relevant. This study advances our understanding of the role played by global structure in virus genome function and provides a model to further investigate its in role virus reproduction. We anticipate that organizational principles revealed by this investigation will be applicable to other viruses.

Long-distance RNA-RNA interactions in the coronavirus genome form high-order structures promoting discontinuous RNA synthesis during transcription

Coronavirus (CoV) transcription requires a high-frequency recombination process that links newly synthesized minus-strand subgenomic RNA copies to the leader region, which is present only once, at the 5' end of the genome. This discontinuous RNA synthesis step is based on the complementarity between the transcription-regulating sequences (TRSs) at the leader region and those preceding each gene in the nascent minus-strand RNA. Furthermore, the template switch requires the physical proximity of RNA genome domains located between 20,000 and 30,000 nucleotides apart. In this report, it is shown that the efficacy of this recombination step is promoted by novel additional long-distance RNA-RNA interactions between RNA motifs located close to the TRSs controlling the expression of each gene and their complementary sequences mapping close to the 5' end of the genome. These interactions would bring together the motifs involved in the recombination process. This finding indicates that the formation of high-order RNA structures in the CoV genome is necessary to control the expression of at least the viral N gene. The requirement of these long-distance interactions for transcription was shown by the engineering of CoV replicons in which the complementarity between the newly identified sequences was disrupted. Furthermore, disruption of complementarity in mutant viruses led to mutations that restored complementarity, wild-type transcription levels, and viral titers by passage in cell cultures. The relevance of these high-order structures for virus transcription is reinforced by the phylogenetic conservation of the involved RNA motifs in CoVs.

Identification of a noncanonically transcribed subgenomic mRNA of infectious bronchitis virus and other gammacoronaviruses

Coronavirus subgenomic mRNA (sgmRNA) synthesis occurs via a process of discontinuous transcription involving transcription regulatory sequences (TRSs) located in the 5' leader sequence (TRS-L) and upstream of each structural and group-specific gene (TRS-B). Several gammacoronaviruses including infectious bronchitis virus (IBV) contain a putative open reading frame (ORF), localized between the M gene and gene 5, which is controversial due to the perceived absence of a TRS. We have studied the transcription of a novel sgmRNA associated with this potential ORF and found it to be transcribed via a previously unidentified noncanonical TRS-B. Using an IBV reverse genetics system, we demonstrated that the template-switching event during intergenic region (IR) sgmRNA synthesis occurs at the 5' end of the noncanonical TRS-B and recombines between nucleotides 5 and 6 of the 8-nucleotide consensus TRS-L. Introduction of a complete TRS-B showed that higher transcription levels are achieved by increasing the number of nucleotide matches between TRS-L and TRS-B. Translation of a protein from the sgmRNA was demonstrated using enhanced green fluorescent protein, suggesting the translation of a fifth, novel, group-specific protein for IBV. This study has resolved an issue concerning the number of ORFs expressed by members of the Gammacoronavirus genus and proposes the existence of a fifth IBV accessory protein. We confirmed previous reports that coronaviruses can produce sgmRNAs from noncanonical TRS-Bs, which may expand their repertoire of proteins. We also demonstrated that noncanonical TRS-Bs may provide a mechanism by which coronaviruses can control protein expression levels by reducing sgmRNA synthesis.