Stem-loop III in the 5' untranslated region is a cis-acting element in bovine coronavirus defective interfering RNA replication

Alpha- and betacoronavirus cis-acting RNA elements

Coronaviruses have exceptionally large RNA genomes and employ multiprotein replication/transcription complexes to orchestrate specific steps of viral RNA genome replication and expression. Most of these processes involve viral cis-acting RNA elements that are engaged in vital RNA-RNA and/or RNA-protein interactions. Over the past years, a large number of studies provided interesting new insight into the structures and, to a lesser extent, functions of specific RNA elements for representative coronaviruses, and there is evidence to suggest that (a majority of) these RNA elements are conserved across genetically divergent coronavirus genera. It is becoming increasingly clear that at least some of these elements do not function in isolation but operate through complex and highly dynamic RNA-RNA interactions. This article reviews structural and functional aspects of cis-acting RNA elements conserved in alpha- and betacoronavirus 5'- and 3'-terminal genome regions, focusing on their critical roles in viral RNA synthesis and gene expression.

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.

Targeting the conserved coronavirus octamer motif GGAAGAGC is a strategy for the development of coronavirus vaccine

BACKGROUND

Coronaviruses are pathogens of humans and animals that cause widespread and costly diseases. The development of effective strategies to combat the threat of coronaviruses is therefore a top priority. The conserved coronavirus octamer motif 5'GGAAGAGC3' exists in the 3' untranslated region of all identified coronaviruses. In the current study, we aimed to examine whether targeting the coronavirus octamer motif GGAAGAGC is a promising approach to develop coronavirus vaccine.

METHODS

Plaque assays were used to determine the titers of mouse hepatitis virus (MHV)-A59 octamer mutant (MHVoctm) and wild-type (wt) MHV-A59 (MHVwt). Western blotting was used for the determination of translation efficiency of MHVoctm and MHVwt. Plaque assays and RT-qPCR were employed to examine whether MHVoctm was more sensitive to interferon treatment than MHVwt. Weight loss, clinical signs, survival rate, viral RNA detection and histopathological examination were used to evaluate whether MHVoctm was a vaccine candidate against MHVwt infection in BALB/c mice.

RESULTS

In this study, we showed that (i) the MHVoctm with mutation of coronavirus octamer was able to grow to high titers but attenuated in mice, (ii) with the reduced multiplicity of infection (MOI), the difference in gene expression between MHVoctm and MHVwt became more evident in cultured cells, (iii) MHVoctm was more sensitive to interferon treatment than MHVwt and (iv) mice inoculated with MHVoctm were protected from MHVwt infection.

CONCLUSIONS

Based on the results obtained from cultured cells, it was suggested that the synergistic effects of octamer mutation, multiplicity of infection and immune response may be a mechanism explaining the distinct phenotypes of octamer-mutated coronavirus in cell culture and mice. In addition, targeting the conserved coronavirus octamer motif is a strategy for development of coronavirus vaccine. Since the conserved octamer exists in all coronaviruses, this strategy of targeting the conserved octamer motif can also be applied to other human and animal coronaviruses for the development of coronavirus vaccines, especially the emergence of novel coronaviruses such as SARS-CoV-2, saving time and cost for vaccine development and disease control.

High-resolution structure of stem-loop 4 from the 5'-UTR of SARS-CoV-2 solved by solution state NMR

We present the high-resolution structure of stem-loop 4 of the 5'-untranslated region (5_SL4) of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) genome solved by solution state nuclear magnetic resonance spectroscopy. 5_SL4 adopts an extended rod-like structure with a single flexible looped-out nucleotide and two mixed tandem mismatches, each composed of a G•U wobble base pair and a pyrimidine•pyrimidine mismatch, which are incorporated into the stem-loop structure. Both the tandem mismatches and the looped-out residue destabilize the stem-loop structure locally. Their distribution along the 5_SL4 stem-loop suggests a role of these non-canonical elements in retaining functionally important structural plasticity in particular with regard to the accessibility of the start codon of an upstream open reading frame located in the RNA's apical loop. The apical loop-although mostly flexible-harbors residual structural features suggesting an additional role in molecular recognition processes. 5_SL4 is highly conserved among the different variants of SARS-CoV-2 and can be targeted by small molecule ligands, which it binds with intermediate affinity in the vicinity of the non-canonical elements within the stem-loop structure.

Biological characterization of coronavirus noncanonical transcripts in vitro and in vivo

BACKGROUND

In addition to the well-known coronavirus genomes and subgenomic mRNAs, the existence of other coronavirus RNA species, which are collectively referred to as noncanonical transcripts, has been suggested; however, their biological characteristics have not yet been experimentally validated in vitro and in vivo.

METHODS

To comprehensively determine the amounts, species and structures of noncanonical transcripts for bovine coronavirus in HRT-18 cells and mouse hepatitis virus A59, a mouse coronavirus, in mouse L cells and mice, nanopore direct RNA sequencing was employed. To experimentally validate the synthesis of noncanonical transcripts under regular infection, Northern blotting was performed. Both Northern blotting and nanopore direct RNA sequencing were also applied to examine the reproducibility of noncanonical transcripts. In addition, Northern blotting was also employed to determine the regulatory features of noncanonical transcripts under different infection conditions, including different cells, multiplicities of infection (MOIs) and coronavirus strains.

RESULTS

In the current study, we (i) experimentally determined that coronavirus noncanonical transcripts were abundantly synthesized, (ii) classified the noncanonical transcripts into seven populations based on their structures and potential synthesis mechanisms, (iii) showed that the species and amounts of the noncanonical transcripts were reproducible during regular infection but regulated in altered infection environments, (iv) revealed that coronaviruses may employ various mechanisms to synthesize noncanonical transcripts, and (v) found that the biological characteristics of coronavirus noncanonical transcripts were similar between in vitro and in vivo conditions.

CONCLUSIONS

The biological characteristics of noncanonical coronavirus transcripts were experimentally validated for the first time. The identified features of noncanonical transcripts in terms of abundance, reproducibility and variety extend the current model for coronavirus gene expression. The capability of coronaviruses to regulate the species and amounts of noncanonical transcripts may contribute to the pathogenesis of coronaviruses during infection, posing potential challenges in disease control. Thus, the biology of noncanonical transcripts both in vitro and in vivo revealed here can provide a database for biological research, contributing to the development of antiviral strategies.

COVID-19: a wreak havoc across the globe

Coronavirus disease (COVID-19) is an infectious airborne viral pneumonia caused by a novel virus belonging to the family coronaviridae. On February 11, 2019, the Internal Committee on Taxonomy of Virus (ICTV) announced the name of the novel virus as "severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). One of the proteins present on its membrane i.e. the Spike protein is responsible for the attachment of the virus to the host. It spreads through the salivary droplets released when an infected person sneezes or coughs. The best way to slow down the disease is by protecting self by washing hands and using the disinfectant. Most of the infected people experience mild to moderate breathing issues. Serious illness might develop in people with underlying cardiovascular problems, diabetes and other immuno-compromised diseases. To date, there is no effective medicine available in the market which is effective in COVID-19. However, healthcare professionals are using ritonavir, flavipiravir, lopinavir, hydroxychloroquine and remdesivir. Along with the medicines, some countries are using convalescent plasma and mesenchymal stem cells for treatment. Till date, it has claimed millions of death worldwide. In this detailed review, we have discussed the structure of SARS-CoV-2, essential proteins, its lifecycle, transmission, symptoms, pathology, clinical features, diagnosis, prevention, treatment and epidemiology of the disease.

Evolution of naturally arising SARS-CoV-2 defective interfering particles

Defective interfering (DI) particles arise during virus propagation, are conditional on parental virus for replication and packaging, and interfere with viral expansion. There is much interest in developing DIs as anti-viral agents. Here we characterize DI particles that arose following serial passaging of SARS-CoV-2 at high multiplicity of infection. The prominent DIs identified have lost ~84% of the SARS-CoV-2 genome and are capable of attenuating parental viral titers. Synthetic variants of the DI genomes also interfere with infection and can be used as conditional, gene delivery vehicles. In addition, the DI genomes encode an Nsp1-10 fusion protein capable of attenuating viral replication. These results identify naturally selected defective viral genomes that emerged and stably propagated in the presence of parental virus.

Genomes from defective interfering (DI) particles following serial passaging of SARS-CoV-2 reveal a fusion protein that attenuates viral replication. Synthetic, recombinant DI genomes are designed to interfere with SARS-CoV-2 replication.

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

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

Prevention and Control of Porcine Epidemic Diarrhea: The Development of Recombination-Resistant Live Attenuated Vaccines

Porcine epidemic diarrhea (PED), causing up to 100% mortality in neonatal pigs, is a highly contagious enteric disease caused by PED virus (PEDV). The highly virulent genogroup 2 (G2) PEDV emerged in 2010 and has caused huge economic losses to the pork industry globally. It was first reported in the US in 2013, caused country-wide outbreaks, and posed tremendous hardship for many pork producers in 2013–2014. Vaccination of pregnant sows/gilts with live attenuated vaccines (LAVs) is the most effective strategy to induce lactogenic immunity in the sows/gilts and provide a passive protection via the colostrum and milk to suckling piglets against PED. However, there are still no safe and effective vaccines available after about one decade of endeavor. One of the biggest concerns is the potential reversion to virulence of an LAV in the field. In this review, we summarize the status and the major obstacles in PEDV LAV development. We also discuss the function of the transcriptional regulatory sequences in PEDV transcription, contributing to recombination, and possible strategies to prevent the reversion of LAVs. This article provides insights into the rational design of a promising LAV without safety issues.

A synthetic defective interfering SARS-CoV-2

Viruses thrive by exploiting the cells they infect, but in order to replicate and infect other cells they must produce viral proteins. As a result, viruses are also susceptible to exploitation by defective versions of themselves that do not produce such proteins. A defective viral genome with deletions in protein-coding genes could still replicate in cells coinfected with full-length viruses. Such a defective genome could even replicate faster due to its shorter size, interfering with the replication of the virus. We have created a synthetic defective interfering version of SARS-CoV-2, the virus causing the Covid-19 pandemic, assembling parts of the viral genome that do not code for any functional protein but enable the genome to be replicated and packaged. This synthetic defective genome replicates three times faster than SARS-CoV-2 in coinfected cells, and interferes with it, reducing the viral load of infected cells by half in 24 hours. The synthetic genome is transmitted as efficiently as the full-length genome, suggesting the location of the putative packaging signal of SARS-CoV-2. A version of such a synthetic construct could be used as a self-promoting antiviral therapy: by enabling replication of the synthetic genome, the virus would promote its own demise.

Structural phylogenetic analysis reveals lineage-specific RNA repetitive structural motifs in all coronaviruses and associated variations in SARS-CoV-2

In many single-stranded (ss) RNA viruses, the cis-acting packaging signal that confers selectivity genome packaging usually encompasses short structured RNA repeats. These structural units, termed repetitive structural motifs (RSMs), potentially mediate capsid assembly by specific RNA–protein interactions. However, general knowledge of the conservation and/or the diversity of RSMs in the positive-sense ssRNA coronaviruses (CoVs) is limited. By performing structural phylogenetic analysis, we identified a variety of RSMs in nearly all CoV genomic RNAs, which are exclusively located in the 5′-untranslated regions (UTRs) and/or in the inter-domain regions of poly-protein 1ab coding sequences in a lineage-specific manner. In all alpha- and beta-CoVs, except for Embecovirus spp, two to four copies of 5′-gUUYCGUc-3′ RSMs displaying conserved hexa-loop sequences were generally identified in Stem-loop 5 (SL5) located in the 5′-UTRs of genomic RNAs. In Embecovirus spp., however, two to eight copies of 5′-agc-3′/guAAu RSMs were found in the coding regions of non-structural protein (NSP) 3 and/or NSP15 in open reading frame (ORF) 1ab. In gamma- and delta-CoVs, other types of RSMs were found in several clustered structural elements in 5′-UTRs and/or ORF1ab. The identification of RSM-encompassing structural elements in all CoVs suggests that these RNA elements play fundamental roles in the life cycle of CoVs. In the recently emerged SARS-CoV-2, beta-CoV-specific RSMs are also found in its SL5, displaying two copies of 5′-gUUUCGUc-3′ motifs. However, multiple sequence alignment reveals that the majority of SARS-CoV-2 possesses a variant RSM harboring SL5b C241U, and intriguingly, several variations in the coding sequences of viral proteins, such as Nsp12 P323L, S protein D614G, and N protein R203K-G204R, are concurrently found with such variant RSM. In conclusion, the comprehensive exploration for RSMs reveals phylogenetic insights into the RNA structural elements in CoVs as a whole and provides a new perspective on variations currently found in SARS-CoV-2.

An overview of basic molecular biology of SARS-CoV-2 and current COVID-19 prevention strategies

Coronavirus Disease 2019 (COVID-19) manifests as extreme acute respiratory conditions caused by a novel beta coronavirus named severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) which is reported to be the seventh coronavirus to infect humans. Like other SARS-CoVs it has a large positive-stranded RNA genome. But, specific furin site in the spike protein, mutation prone and phylogenetically mess open reading frame1ab (Orf1ab) separates SARS-CoV-2 from other RNA viruses. Since the outbreak (February-March 2020), researchers, scientists, and medical professionals are inspecting all possible facts and aspects including its replication, detection, and prevention strategies. This led to the prompt identification of its basic biology, genome characterization, structural and expression based functional information of proteins, and utilization of this information in optimizing strategies to prevent its spread. This review summarizes the recent updates on the basic molecular biology of SARS-CoV-2 and prevention strategies undertaken worldwide to tackle COVID-19. This recent information can be implemented for the development and designing of therapeutics against SARS-CoV-2.

Dengue virus-free defective interfering particles have potent and broad anti-dengue virus activity

Dengue virus (DENV) is spread from human to human through the bite of the female Aedes aegypti mosquito and leads to about 100 million clinical infections yearly. Treatment options and vaccine availability for DENV are limited. Defective interfering particles (DIPs) are considered a promising antiviral approach but infectious virus contamination has limited their development. Here, a DENV-derived DIP production cell line was developed that continuously produced DENV-free DIPs. The DIPs contained and could deliver to cells a DENV serotype 2 subgenomic defective-interfering RNA, which was originally discovered in DENV infected patients. The DIPs released into cell culture supernatant were purified and could potently inhibit replication of all DENV serotypes in cells. Antiviral therapeutics are limited for many viral infection. The DIP system described could be re-purposed to make antiviral DIPs for many other RNA viruses such as SARS-CoV-2, yellow fever, West Nile and Zika viruses.

Analysis of Emerging Variants in Structured Regions of the SARS-CoV-2 Genome

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two sets of SARS-CoV-2 genomic sequences to identify emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, 20 variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5' untranslated region (UTR), including a group of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-ss molecular switch in the 3'UTR. Finally, 5 variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M (stem loop 2 m) selfish genetic element, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. Our analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences.

De novo 3D models of SARS-CoV-2 RNA elements from consensus experimental secondary structures

The rapid spread of COVID-19 is motivating development of antivirals targeting conserved SARS-CoV-2 molecular machinery. The SARS-CoV-2 genome includes conserved RNA elements that offer potential small-molecule drug targets, but most of their 3D structures have not been experimentally characterized. Here, we provide a compilation of chemical mapping data from our and other labs, secondary structure models, and 3D model ensembles based on Rosetta's FARFAR2 algorithm for SARS-CoV-2 RNA regions including the individual stems SL1-8 in the extended 5' UTR; the reverse complement of the 5' UTR SL1-4; the frameshift stimulating element (FSE); and the extended pseudoknot, hypervariable region, and s2m of the 3' UTR. For eleven of these elements (the stems in SL1-8, reverse complement of SL1-4, FSE, s2m and 3' UTR pseudoknot), modeling convergence supports the accuracy of predicted low energy states; subsequent cryo-EM characterization of the FSE confirms modeling accuracy. To aid efforts to discover small molecule RNA binders guided by computational models, we provide a second set of similarly prepared models for RNA riboswitches that bind small molecules. Both datasets ('FARFAR2-SARS-CoV-2', https://github.com/DasLab/FARFAR2-SARS-CoV-2; and 'FARFAR2-Apo-Riboswitch', at https://github.com/DasLab/FARFAR2-Apo-Riboswitch') include up to 400 models for each RNA element, which may facilitate drug discovery approaches targeting dynamic ensembles of RNA molecules.

Secondary structure of the SARS-CoV-2 5'-UTR

The SARS-CoV-2, a positive-sense single-stranded RNA Coronavirus, is a global threat to human health. Thus, understanding its life cycle mechanistically would be important to facilitate the design of antiviral drugs. A key aspect of viral progression is the synthesis of viral proteins by the ribosome of the human host. In Coronaviruses, this process is regulated by the viral 5' and 3' untranslated regions (UTRs), but the precise regulatory mechanism has not yet been well understood. In particular, the 5'-UTR of the viral genome is most likely involved in translation initiation of viral proteins. Here, we performed inline probing and RNase V1 probing to establish a model of the secondary structure of SARS-CoV-2 5'-UTR. We found that the 5'-UTR contains stable structures including a very stable four-way junction close to the AUG start codon. Sequence alignment analysis of SARS-CoV-2 variants 5'-UTRs revealed a highly conserved structure with few co-variations that confirmed our secondary structure model based on probing experiments.

A Transcription Regulatory Sequence in the 5' Untranslated Region of SARS-CoV-2 Is Vital for Virus Replication with an Altered Evolutionary Pattern against Human Inhibitory MicroRNAs

Our knowledge of the evolution and the role of untranslated region (UTR) in SARS-CoV-2 pathogenicity is very limited. Leader sequence, originated from UTR, is found at the 5' ends of all encoded SARS-CoV-2 transcripts, highlighting its importance. Here, evolution of leader sequence was compared between human pathogenic and non-pathogenic coronaviruses. Then, profiling of microRNAs that can inactivate the key UTR regions of coronaviruses was carried out. A distinguished pattern of evolution in leader sequence of SARS-CoV-2 was found. Mining all available microRNA families against leader sequences of coronaviruses resulted in discovery of 39 microRNAs with a stable thermodynamic binding energy. Notably, SARS-CoV-2 had a lower binding stability against microRNAs. hsa-MIR-5004-3p was the only human microRNA able to target the leader sequence of SARS and to a lesser extent, also SARS-CoV-2. However, its binding stability decreased remarkably in SARS-COV-2. We found some plant microRNAs with low and stable binding energy against SARS-COV-2. Meta-analysis documented a significant (p < 0.01) decline in the expression of MIR-5004-3p after SARS-COV-2 infection in trachea, lung biopsy, and bronchial organoids as well as lung-derived Calu-3 and A549 cells. The paucity of the innate human inhibitory microRNAs to bind to leader sequence of SARS-CoV-2 can contribute to its high replication in infected human cells.

SARS-CoV-2: A critical review of its history, pathogenesis, transmission, diagnosis and treatment

The outbreak of the deadly virus (novel coronavirus or Severe Acute Respiratory Syndrome Coronavirus-2) that emerged in December 2019, remained a controversial subject of intense speculations regarding its origin, became a worldwide health problem resulting in serious coronavirus disease 2019 (acronym COVID-19). The concern regarding this new viral strain "Severe Acute Respiratory Syndrome Coronavirus-2" (acronym SARS-CoV-2) and diseases it causes (COVID-19) is well deserved at all levels. The incidence of COVID-19 infection and infectious patients are increasing at a high rate. Coronaviruses (CoVs), enclosed positive-sense RNA viruses, are distinguished by club-like spikes extending from their surface, an exceptionally large genome of RNA, and a special mechanism for replication. Coronaviruses are associated with a broad variety of human and other animal diseases spanning from enteritis in cattle and pigs and upper chicken respiratory disease to extremely lethal human respiratory infections. With World Health Organization (WHO) declaring COVID-19 as pandemic, we deemed it necessary to provide a detailed review of coronaviruses discussing their history, current situation, coronavirus classification, pathogenesis, structure, mode of action, diagnosis and treatment, the effect of environmental factors, risk reduction and guidelines to understand the virus and develop ways to control it.

Coronaviruses: An Updated Overview of Their Replication and Pathogenesis

Coronaviruses (CoVs), enveloped positive-sense RNA viruses, are characterized by club-like spikes that project from their surface, an unusually large RNA genome, and a unique replication strategy. CoVs cause a variety of diseases in mammals and birds ranging from enteritis in cows and pigs, and upper respiratory tract and kidney disease in chickens to lethal human respiratory infections. Most recently, the novel coronavirus, SARS-CoV-2, which was first identified in Wuhan, China in December 2019, is the cause of a catastrophic pandemic, COVID-19, with more than 8 million infections diagnosed worldwide by mid-June 2020. Here we provide a brief introduction to CoVs discussing their replication, pathogenicity, and current prevention and treatment strategies. We will also discuss the outbreaks of the highly pathogenic Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and Middle Eastern Respiratory Syndrome Coronavirus (MERS-CoV), which are relevant for understanding COVID-19.

Comprehensive in-vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms

SARS-CoV-2 is the positive-sense RNA virus that causes COVID-19, a disease that has triggered a major human health and economic crisis. The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form stable RNA structures and yet, as much as 97% of its 30 kilobases have not been structurally explored in the context of a viral infection. Our limited knowledge of SARS-CoV-2 genomic architecture is a fundamental limitation to both our mechanistic understanding of coronavirus life cycle and the development of COVID-19 RNA-based therapeutics. Here, we apply a novel long amplicon strategy to determine for the first time the secondary structure of the SARS-CoV-2 RNA genome probed in infected cells. In addition to the conserved structural motifs at the viral termini, we report new structural features like a conformationally flexible programmed ribosomal frameshifting pseudoknot, and a host of novel RNA structures, each of which highlights the importance of studying viral structures in their native genomic context. Our in-depth structural analysis reveals extensive networks of well-folded RNA structures throughout Orf1ab and reveals new aspects of SARS-CoV-2 genome architecture that distinguish it from other single-stranded, positive-sense RNA viruses. Evolutionary analysis of RNA structures in SARS-CoV-2 shows that several features of its genomic structure are conserved across beta coronaviruses and we pinpoint individual regions of well-folded RNA structure that merit downstream functional analysis. The native, complete secondary structure of SAR-CoV-2 presented here is a roadmap that will facilitate focused studies on mechanisms of replication, translation and packaging, and guide the identification of new RNA drug targets against COVID-19.

Analysis of Rapidly Emerging Variants in Structured Regions of the SARS-CoV-2 Genome

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two cohorts of SARS-CoV-2 genomic sequences to identify rapidly emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, twenty variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5'UTR, including a set of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-stable molecular switch in the 3'UTR. Finally, five variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M stem loop, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. This analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences, as rapidly emerging variations in these regions could lead to drug resistance.

Coronavirus Pandemic-Therapy and Vaccines

The current coronavirus COVID-19 pandemic, which originated in Wuhan, China, has raised significant social, psychological and economic concerns in addition to direct medical issues. The rapid spread of severe acute respiratory syndrome-coronavirus (SARS-CoV)-2 to almost every country on the globe and the failure to contain the infections have contributed to fear and panic worldwide. The lack of available and efficient antiviral drugs or vaccines has further worsened the situation. For these reasons, it cannot be overstated that an accelerated effort for the development of novel drugs and vaccines is needed. In this context, novel approaches in both gene therapy and vaccine development are essential. Previous experience from SARS- and MERS-coronavirus vaccine and drug development projects have targeted glycoprotein epitopes, monoclonal antibodies, angiotensin receptor blockers and gene silencing technologies, which may be useful for COVID-19 too. Moreover, existing antivirals used for other types of viral infections have been considered as urgent action is necessary. This review aims at providing a background of coronavirus genetics and biology, examples of therapeutic and vaccine strategies taken and potential innovative novel approaches in progress.

Interaction of coronavirus nucleocapsid protein with the 5'- and 3'-ends of the coronavirus genome is involved in genome circularization and negative-strand RNA synthesis

Synthesis of the negative‐strand ((−)‐strand) counterpart is the first step of coronavirus (CoV) replication; however, the detailed mechanism of the early event and the factors involved remain to be determined. Here, using bovine coronavirus (BCoV)‐defective interfering (DI) RNA, we showed that (a) a poly(A) tail with a length of 15 nucleotides (nt) was sufficient to initiate efficient (−)‐strand RNA synthesis and (b) substitution of the poly(A) tail with poly(U), (C) or (G) only slightly decreased the efficiency of (−)‐strand synthesis. The findings indicate that in addition to the poly(A) tail, other factors acting in trans may also participate in (−)‐strand synthesis. The BCoV nucleocapsid (N) protein, an RNA‐binding protein, was therefore tested as a candidate. Based on dissociation constant (K_d) values, it was found that the binding affinity between N protein, but not poly(A)‐binding protein, and the 3′‐terminal 55 nt plus a poly(A), poly(U), poly(C) or poly(G) tail correlates with the efficiency of (−)‐strand synthesis. Such an association was also evidenced by the binding affinity between the N protein and 5′‐ and 3′‐terminal cis‐acting elements important for (−)‐strand synthesis. Further analysis demonstrated that N protein can act as a bridge to facilitate interaction between the 5′‐ and 3′‐ends of the CoV genome, leading to circularization of the genome. Together, the current study extends our understanding of the mechanism of CoV (−)‐strand RNA synthesis through involvement of N protein and genome circularization and thus may explain why the addition of N protein in trans is required for efficient CoV replication.

Structural and functional conservation of cis-acting RNA elements in coronavirus 5'-terminal genome regions

Structure predictions suggest a partial conservation of RNA structure elements in coronavirus terminal genome regions. Here, we determined the structures of stem-loops (SL) 1 and 2 of two alphacoronaviruses, human coronavirus (HCoV) 229E and NL63, by RNA structure probing and studied the functional relevance of these putative cis-acting elements. HCoV-229E SL1 and SL2 mutants generated by reverse genetics were used to study the effects on viral replication of single-nucleotide substitutions predicted to destabilize the SL1 and SL2 structures. The data provide conclusive evidence for the critical role of SL1 and SL2 in HCoV-229E replication and, in some cases, revealed parallels with previously characterized betacoronavirus SL1 and SL2 elements. Also, we were able to rescue viable HCoV-229E mutants carrying replacements of SL2 with equivalent betacoronavirus structural elements. The data obtained in this study reveal a remarkable degree of structural and functional conservation of 5′-terminal RNA structural elements across coronavirus genus boundaries.

•

Conservation of 5'-terminal SL1 and SL2 elements in alpha- and betacoronaviruses.

•

HCoV-229E reverse genetics data suggest critical role for SL1/2 in viral replication.

•

Replacement of SL2 in HCoV-229E with betacoronavirus equivalents is tolerated.

MERS coronavirus nsp1 participates in an efficient propagation through a specific interaction with viral RNA

MERS-CoV is the only lethal human CoV still endemic in the Arabian Peninsula and neither vaccine nor therapeutics against MERS-CoV infection is available. The nsp1 of CoV is thought to be a major virulence factor because it suppresses protein synthesis through the degradation of host mRNA. In contrast, viral RNA circumvents the nsp1-mediated translational shutoff for an efficient propagation. In this study, we identified amino acid residue in MERS-CoV nsp1 that differ from those of SARS-CoV nsp1, and that appear to be crucial for circumventing the translational shutoff. In addition, reverse genetics analysis suggested the presence of a cis-acting element at the 5'-terminus of the nsp1-coding region, which contributes to the specific recognition of viral RNA that is required for an efficient viral replication. Our results suggest the CoVs share a common mechanism for circumventing the nsp1-mediated translational shutoff.

Characterization of Chinese Porcine Epidemic Diarrhea Virus with Novel Insertions and Deletions in Genome

Outbreaks of porcine epidemic diarrhoea virus (PEDV) have caused great economic losses to the global pig industry. PEDV strains with variants in the spike (S) gene have been reported in several countries. To better understand the molecular epidemiology and genetic diversity of PEDV field isolates, in this study, we characterised the complete genome sequence of a novel PEDV variant JSCZ1601 from a outbreak in China in 2016. The PEDV isolate was 28,033 nucleotides (nt) in length without the polyadenylated sequences. Phylogenetic analysis based on the full-length genome sequence of JSCZ1601 grouped it with the pandemic variants determined post-2010 into group 2 (G2). However, the S gene of JSCZ1601 formed a new subgroup separated from the subgroups containing the other G2 strains. Comparative analysis of the amino acids encoded by the S genes revealed the N-terminal of the deduced JSCZ1601 S protein had a novel two-amino-acid deletion (N58 and S59) compared with all identified genogroups. Further, compared with the reference strains, a 'G' insertion was detected in the 5' terminal of the 5'UTR of the JSCZ1601. The animal experiment revealed that this strain was high pathogenic to neonatal pigs. Taken together, a PEDV strain with the new molecular characterizations and phylogenies was found in mainland China. It is necessary to strengthen the monitoring of PEDV variations.

Coronavirus cis-Acting RNA Elements

Coronaviruses have exceptionally large RNA genomes of approximately 30 kilobases. Genome replication and transcription is mediated by a multisubunit protein complex comprised of more than a dozen virus-encoded proteins. The protein complex is thought to bind specific cis-acting RNA elements primarily located in the 5'- and 3'-terminal genome regions and upstream of the open reading frames located in the 3'-proximal one-third of the genome. Here, we review our current understanding of coronavirus cis-acting RNA elements, focusing on elements required for genome replication and packaging. Recent bioinformatic, biochemical, and genetic studies suggest a previously unknown level of conservation of cis-acting RNA structures among different coronavirus genera and, in some cases, even beyond genus boundaries. Also, there is increasing evidence to suggest that individual cis-acting elements may be part of higher-order RNA structures involving long-range and dynamic RNA-RNA interactions between RNA structural elements separated by thousands of nucleotides in the viral genome. We discuss the structural and functional features of these cis-acting RNA elements and their specific functions in coronavirus RNA synthesis.

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.

Coronaviruses: an overview of their replication and pathogenesis

Coronaviruses (CoVs), enveloped positive-sense RNA viruses, are characterized by club-like spikes that project from their surface, an unusually large RNA genome, and a unique replication strategy. Coronaviruses cause a variety of diseases in mammals and birds ranging from enteritis in cows and pigs and upper respiratory disease in chickens to potentially lethal human respiratory infections. Here we provide a brief introduction to coronaviruses discussing their replication and pathogenicity, and current prevention and treatment strategies. We also discuss the outbreaks of the highly pathogenic Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and the recently identified Middle Eastern Respiratory Syndrome Coronavirus (MERS-CoV).

The structure and functions of coronavirus genomic 3' and 5' ends

Coronaviruses (CoVs) are an important cause of illness in humans and animals. Most human coronaviruses commonly cause relatively mild respiratory illnesses; however two zoonotic coronaviruses, SARS-CoV and MERS-CoV, can cause severe illness and death. Investigations over the past 35 years have illuminated many aspects of coronavirus replication. The focus of this review is the functional analysis of conserved RNA secondary structures in the 5' and 3' of the betacoronavirus genomes. The 5' 350 nucleotides folds into a set of RNA secondary structures which are well conserved, and reverse genetic studies indicate that these structures play an important role in the discontinuous synthesis of subgenomic RNAs in the betacoronaviruses. These cis-acting elements extend 3' of the 5'UTR into ORF1a. The 3'UTR is similarly conserved and contains all of the cis-acting sequences necessary for viral replication. Two competing conformations near the 5' end of the 3'UTR have been shown to make up a potential molecular switch. There is some evidence that an association between the 3' and 5'UTRs is necessary for subgenomic RNA synthesis, but the basis for this association is not yet clear. A number of host RNA proteins have been shown to bind to the 5' and 3' cis-acting regions, but the significance of these in viral replication is not clear. Two viral proteins have been identified as binding to the 5' cis-acting region, nsp1 and N protein. A genetic interaction between nsp8 and nsp9 and the region of the 3'UTR that contains the putative molecular switch suggests that these two proteins bind to this region.

Coronaviruses: an overview of their replication and pathogenesis

Coronaviruses (CoVs), enveloped positive-sense RNA viruses, are characterized by club-like spikes that project from their surface, an unusually large RNA genome, and a unique replication strategy. Coronaviruses cause a variety of diseases in mammals and birds ranging from enteritis in cows and pigs and upper respiratory disease in chickens to potentially lethal human respiratory infections. Here we provide a brief introduction to coronaviruses discussing their replication and pathogenicity, and current prevention and treatment strategies. We also discuss the outbreaks of the highly pathogenic Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and the recently identified Middle Eastern Respiratory Syndrome Coronavirus (MERS-CoV).

SHAPE analysis of the RNA secondary structure of the Mouse Hepatitis Virus 5' untranslated region and N-terminal nsp1 coding sequences

SHAPE technology was used to analyze RNA secondary structure of the 5' most 474 nts of the MHV-A59 genome encompassing the minimal 5' cis-acting region required for defective interfering RNA replication. The structures generated were in agreement with previous characterizations of SL1 through SL4 and two recently predicted secondary structure elements, S5 and SL5A. SHAPE provided biochemical support for four additional stem-loops not previously functionally investigated in MHV. Secondary structure predictions for 5' regions of MHV-A59, BCoV and SARS-CoV were similar despite high sequence divergence. The pattern of SHAPE reactivity of in virio genomic RNA, ex virio genomic RNA, and in vitro synthesized RNA was similar, suggesting that binding of N protein or other proteins to virion RNA fails to protect the RNA from reaction with lipid permeable SHAPE reagent. Reverse genetic experiments suggested that SL5C and SL6 within the nsp1 coding sequence are not required for viral replication.

RNA structure analysis of alphacoronavirus terminal genome regions

•

Review of current knowledge of cis-acting RNA elements essential to coronavirus replication.

•

Identification of RNA structural elements in alphacoronavirus terminal genome regions.

•

Discussion of intra- and intergeneric conservation of genomic cis-acting RNA elements in alpha- and betacoronaviruses.

Coronavirus genome replication is mediated by a multi-subunit protein complex that is comprised of more than a dozen virally encoded and several cellular proteins. Interactions of the viral replicase complex with cis-acting RNA elements located in the 5′ and 3′-terminal genome regions ensure the specific replication of viral RNA. Over the past years, boundaries and structures of cis-acting RNA elements required for coronavirus genome replication have been extensively characterized in betacoronaviruses and, to a lesser extent, other coronavirus genera. Here, we review our current understanding of coronavirus cis-acting elements located in the terminal genome regions and use a combination of bioinformatic and RNA structure probing studies to identify and characterize putative cis-acting RNA elements in alphacoronaviruses. The study suggests significant RNA structure conservation among members of the genus Alphacoronavirus but also across genus boundaries. Overall, the conservation pattern identified for 5′ and 3′-terminal RNA structural elements in the genomes of alpha- and betacoronaviruses is in agreement with the widely used replicase polyprotein-based classification of the Coronavirinae, suggesting co-evolution of the coronavirus replication machinery with cognate cis-acting RNA elements.

Genomic characteristics and changes of avian infectious bronchitis virus strain CK/CH/LDL/97I after serial passages in chicken embryos

Background

We previously attenuated the infectious bronchitis virus (IBV) strain CK/CH/LDL/97I and found that it can convey protection against the homologous pathogenic virus.

Objective

To compare the full-length genome sequences of the Chinese IBV strain CK/CH/LDL/97I and its embryo-passaged, attenuated level to identify sequence substitutions responsible for the attenuation and define markers of attenuation.

Methods

The full-length genomes of CK/CH/LDL/97I P5 and P115 were amplified and sequenced. The sequences were assembled and compared using the MEGALIGN program (DNAStar) and a phylogenetic tree was constructed using MEGA4 software.

Results

The CK/CH/LDL/97I virus population contained subpopulations with a mixture of genetic mutants. Changes were observed in nsp4, nsp9, nsp11/12, nsp14, nsp15, nsp16, and ORF3a, but these did not result in amino acid substitutions or did not show functional variations. Amino acid substitutions occurred in the remaining genes between P5 and P115; most were found in the S region, and some of the nucleotide mutations resulted in amino acid substitutions. Among the 9 nsps in the ORF1 region, nsp3 contained the most nucleotide substitutions.

Conclusions

Sequence variations in different genes, especially the S gene and nsp3, in the genomes of CK/CH/LDL/97I viruses might contribute to differences in viral replication, pathogenicity, antigenicity, immunogenicity, and tissue tropism.

Dependence of coronavirus RNA replication on an NH2-terminal partial nonstructural protein 1 in cis

Genomes of positive (+)-strand RNA viruses use cis-acting signals to direct both translation and replication. Here we examine two 5'-proximal cis-replication signals of different character in a defective interfering (DI) RNA of the bovine coronavirus (BCoV) that map within a 322-nucleotide (nt) sequence (136 nt from the genomic 5' untranslated region and 186 nt from the nonstructural protein 1 [nsp1]-coding region) not found in the otherwise-identical nonreplicating subgenomic mRNA7 (sgmRNA7). The natural DI RNA is structurally a fusion of the two ends of the BCoV genome that results in a single open reading frame between a partial nsp1-coding region and the entire N gene. (i) In the first examination, mutation analyses of a recently discovered long-range RNA-RNA base-paired structure between the 5' untranslated region and the partial nsp1-coding region showed that it, possibly in concert with adjacent stem-loops, is a cis-acting replication signal in the (+) strand. We postulate that the higher-order structure promotes (+)-strand synthesis. (ii) In the second examination, analyses of multiple frame shifts, truncations, and point mutations within the partial nsp1-coding region showed that synthesis of a PEFP core amino acid sequence within a group A lineage betacoronavirus-conserved NH2-proximal WAPEFPWM domain is required in cis for DI RNA replication. We postulate that the nascent protein, as part of an RNA-associated translating complex, acts to direct the DI RNA to a critical site, enabling RNA replication. We suggest that these results have implications for viral genome replication and explain, in part, why coronavirus sgmRNAs fail to replicate.

IMPORTANCE

cis-Acting RNA and protein structures that regulate (+)-strand RNA virus genome synthesis are potential sites for blocking virus replication. Here we describe two: a previously suspected 5'-proximal long-range higher-order RNA structure and a novel nascent NH2-terminal protein component of nsp1 that are common among betacoronaviruses of group A lineage.

Identification of cis-acting elements on positive-strand subgenomic mRNA required for the synthesis of negative-strand counterpart in bovine coronavirus

It has been demonstrated that, in addition to genomic RNA, sgmRNA is able to serve as a template for the synthesis of the negative-strand [(−)-strand] complement. However, the cis-acting elements on the positive-strand [(+)-strand] sgmRNA required for (−)-strand sgmRNA synthesis have not yet been systematically identified. In this study, we employed real-time quantitative reverse transcription polymerase chain reaction to analyze the cis-acting elements on bovine coronavirus (BCoV) sgmRNA 7 required for the synthesis of its (−)-strand counterpart by deletion mutagenesis. The major findings are as follows. (1) Deletion of the 5'-terminal leader sequence on sgmRNA 7 decreased the synthesis of the (−)-strand sgmRNA complement. (2) Deletions of the 3' untranslated region (UTR) bulged stem-loop showed no effect on (−)-strand sgmRNA synthesis; however, deletion of the 3' UTR pseudoknot decreased the yield of (−)-strand sgmRNA. (3) Nucleotides positioned from −15 to −34 of the sgmRNA 7 3'-terminal region are required for efficient (−)-strand sgmRNA synthesis. (4) Nucleotide species at the 3'-most position (−1) of sgmRNA 7 is correlated to the efficiency of (−)-strand sgmRNA synthesis. These results together suggest, in principle, that the 5'- and 3'-terminal sequences on sgmRNA 7 harbor cis-acting elements are critical for efficient (−)-strand sgmRNA synthesis in BCoV.

The 3'-terminal 55 nucleotides of bovine coronavirus defective interfering RNA harbor cis-acting elements required for both negative- and positive-strand RNA synthesis

The synthesis of the negative-strand [(−)-strand] complement of the ∼30 kilobase, positive-strand [(+)-strand] coronaviral genome is a necessary early step for genome replication. The identification of cis-acting elements required for (−)-strand RNA synthesis in coronaviruses, however, has been hampered due to insufficiencies in the techniques used to detect the (−)-strand RNA species. Here, we employed a method of head-to-tail ligation and real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) to detect and quantitate the synthesis of bovine coronavirus (BCoV) defective interfering (DI) RNA (−) strands. Furthermore, using the aforementioned techniques along with Northern blot assay, we specifically defined the cis-acting RNA elements within the 3′-terminal 55 nucleotides (nts) which function in the synthesis of (−)- or (+)-strand BCoV DI RNA. The major findings are as follows: (i) nts from -5 to -39 within the 3′-terminal 55 nts are the cis-acting elements responsible for (−)-strand BCoV DI RNA synthesis, (ii) nts from −3 to −34 within the 3′-terminal 55 nts are cis-acting elements required for (+)-strand BCoV DI RNA synthesis, and (iii) the nucleotide species at the 3′-most position (−1) is important, but not critical, for both (−)- and (+)-strand BCoV DI RNA synthesis. These results demonstrate that the 3′-terminal 55 nts in BCoV DI RNA harbor cis-acting RNA elements required for both (−)- and (+)-strand DI RNA synthesis and extend our knowledge on the mechanisms of coronavirus replication. The method of head-to-tail ligation and qRT-PCR employed in the study may also be applied to identify other cis-acting elements required for (−)-strand RNA synthesis in coronaviruses.

Origin and characteristics of the recombinant novel avian infectious bronchitis coronavirus isolate ck/CH/LJL/111054

Recombination among infectious bronchitis viruses (IBVs), coupled with point mutations, insertions, and deletions that occur in the genome, is thought to contribute to the emergence of new IBV variants. In this study an IBV, ck/CH/LJL/111054, was isolated from a H120-vaccinated chicken, which presented with a suspected IBV infection. Phylogenetic analysis of the S1 subunit sequence confirmed that strain ck/CH/LJL/111054 is of the Connecticut-type; however, further extensive full-length genomic analysis identified the occurrence of recombination events. Therefore, strain ck/CH/LJL/111054 may have originated from recombination events between Conn- and Mass-like strains at three recombination breakpoints: two located within the nsp3 gene sequence and one in the nsp12 gene sequence. Further, the uptake of the 5' untranslated regions, nsp2, parts of nsp3, nsp4-11, and parts of nsp 12 from Mass-like virus by ck/CH/LJL/111054 might have resulted in changes in viral replication efficiency rather than antigenic changes, via cross-neutralization analysis with the H120 strain. Recombination events coupled with the accumulation of mutations in the ck/CH/LJL/111054 genome may account for its increased virulence in specific-pathogen free chickens.

A leaderless genome identified during persistent bovine coronavirus infection is associated with attenuation of gene expression

The establishment of persistent viral infection is often associated with the selection of one or more mutant viruses. For example, it has been found that an intraleader open reading frame (ORF) in genomic and subgenomic mRNA (sgmRNA) molecules is selected during bovine coronavirus (BCoV) persistence which leads to translation attenuation of the downstream ORF. Here, we report the unexpected identification of leaderless genomes, in addition to leader-containing genomes, in a cell culture persistently infected with BCoV. The discovery was made by using a head-to-tail ligation method that examines genomic 5′-terminal sequences at different times postinfection. Functional analyses of the leaderless genomic RNA in a BCoV defective interfering (DI) RNA revealed that (1) the leaderless genome was able to serve as a template for the synthesis of negative-strand genome, although it cannot perform replicative positive-strand genomic RNA synthesis, and (2) the leaderless genome retained its function in translation and transcription, although the efficiency of these processes was impaired. Therefore, this previously unidentified leaderless genome is associated with the attenuation of genome expression. Whether the leaderless genome contributes to the establishment of persistent infection remains to be determined.

Evidence of pervasive biologically functional secondary structures within the genomes of eukaryotic single-stranded DNA viruses

Single-stranded DNA (ssDNA) viruses have genomes that are potentially capable of forming complex secondary structures through Watson-Crick base pairing between their constituent nucleotides. A few of the structural elements formed by such base pairings are, in fact, known to have important functions during the replication of many ssDNA viruses. Unknown, however, are (i) whether numerous additional ssDNA virus genomic structural elements predicted to exist by computational DNA folding methods actually exist and (ii) whether those structures that do exist have any biological relevance. We therefore computationally inferred lists of the most evolutionarily conserved structures within a diverse selection of animal- and plant-infecting ssDNA viruses drawn from the families Circoviridae, Anelloviridae, Parvoviridae, Nanoviridae, and Geminiviridae and analyzed these for evidence of natural selection favoring the maintenance of these structures. While we find evidence that is consistent with purifying selection being stronger at nucleotide sites that are predicted to be base paired than at sites predicted to be unpaired, we also find strong associations between sites that are predicted to pair with one another and site pairs that are apparently coevolving in a complementary fashion. Collectively, these results indicate that natural selection actively preserves much of the pervasive secondary structure that is evident within eukaryote-infecting ssDNA virus genomes and, therefore, that much of this structure is biologically functional. Lastly, we provide examples of various highly conserved but completely uncharacterized structural elements that likely have important functions within some of the ssDNA virus genomes analyzed here.

Reselection of a genomic upstream open reading frame in mouse hepatitis coronavirus 5'-untranslated-region mutants

An AUG-initiated upstream open reading frame (uORF) encoding a potential polypeptide of 3 to 13 amino acids (aa) is found within the 5' untranslated region (UTR) of >75% of coronavirus genomes based on 38 reference strains. Potential CUG-initiated uORFs are also found in many strains. The AUG-initiated uORF is presumably translated following genomic 5'-end cap-dependent ribosomal scanning, but its function is unknown. Here, in a reverse-genetics study with mouse hepatitis coronavirus, the following were observed. (i) When the uORF AUG-initiating codon was replaced with a UAG stop codon along with a U112A mutation to maintain a uORF-harboring stem-loop 4 structure, an unimpaired virus with wild-type (WT) growth kinetics was recovered. However, reversion was found at all mutated sites within five virus passages. (ii) When the uORF was fused with genomic (main) ORF1 by converting three in-frame stop codons to nonstop codons, a uORF-ORF1 fusion protein was made, and virus replicated at WT levels. However, a frameshifting G insertion at virus passage 7 established a slightly 5'-extended original uORF. (iii) When uAUG-eliminating deletions of 20, 30, or 51 nucleotides (nt) were made within stem-loop 4, viable but debilitated virus was recovered. However, a C80U mutation in the first mutant and an A77G mutation in the second appeared by passage 10, which generated alternate uORFs that correlated with restored WT growth kinetics. In vitro, the uORF-disrupting nondeletion mutants showed enhanced translation of the downstream ORF1 compared with the WT. These results together suggest that the uORF represses ORF1 translation yet plays a beneficial but nonessential role in coronavirus replication in cell culture.

Comparative analysis of four Massachusetts type infectious bronchitis coronavirus genomes reveals a novel Massachusetts type strain and evidence of natural recombination in the genome

Four Massachusetts-type (Mass-type) strains of infectious bronchitis coronavirus (IBV) were compared genetically with the pathogenic M41 and H120 vaccine strains using the complete genomic sequences. The results revealed that strains ck/CH/LNM/091017 and ck/CH/LDL/101212 were closely related to the H120 vaccine, which suggests that they might represent re-isolations of vaccine strains or variants of vaccine strains that have resulted from the accumulated point mutations after several passages in chickens. In contrast, strains ck/CH/LHLJ/07VII and ck/CH/LHLJ/100902 had a close genetic relationship with the pathogenic M41 strain. In addition, molecular markers have been identified that distinguish between field and vaccine (or vaccine-like) Mass-type viruses, which may be able to differentiate between field and vaccine strains for diagnostic purposes. Phylogenetic analysis, and pairwise comparison of full-length genomes and the nine genes, identified the occurrence of recombination events in the genome of strain CK/VH/LHLJ/07VII, which suggests that this virus originated from recombination events between M41- and H120-like strains at the switch site located at the 3' end of the nucleocapsid (N) genes. To our knowledge, this is the first time that evidence for the evolution and natural recombination under field conditions between Mass-type pathogenic and vaccinal IBV strains has been documented. These findings provide insights into the emergence and evolution of the Mass-type IB coronaviruses and may help to explain the emergence of Mass-type IBV in chicken flocks all over the world.

Severe acute respiratory syndrome coronavirus nsp1 facilitates efficient propagation in cells through a specific translational shutoff of host mRNA

Severe acute respiratory syndrome (SARS) coronavirus (SCoV) is an enveloped virus containing a single-stranded, positive-sense RNA genome. Nine mRNAs carrying a set of common 5' and 3' untranslated regions (UTR) are synthesized from the incoming viral genomic RNA in cells infected with SCoV. A nonstructural SCoV nsp1 protein causes a severe translational shutoff by binding to the 40S ribosomal subunits. The nsp1-40S ribosome complex further induces an endonucleolytic cleavage near the 5'UTR of host mRNA. However, the mechanism by which SCoV viral proteins are efficiently produced in infected cells in which host protein synthesis is impaired by nsp1 is unknown. In this study, we investigated the role of the viral UTRs in evasion of the nsp1-mediated shutoff. Luciferase activities were significantly suppressed in cells expressing nsp1 together with the mRNA carrying a luciferase gene, while nsp1 failed to suppress luciferase activities of the mRNA flanked by the 5'UTR of SCoV. An RNA-protein binding assay and RNA decay assay revealed that nsp1 bound to stem-loop 1 (SL1) in the 5'UTR of SCoV RNA and that the specific interaction with nsp1 stabilized the mRNA carrying SL1. Furthermore, experiments using an SCoV replicon system showed that the specific interaction enhanced the SCoV replication. The specific interaction of nsp1 with SL1 is an important strategy to facilitate efficient viral gene expression in infected cells, in which nsp1 suppresses host gene expression. Our data indicate a novel mechanism of viral gene expression control by nsp1 and give new insight into understanding the pathogenesis of SARS.

Genetic evidence of a long-range RNA-RNA interaction between the genomic 5' untranslated region and the nonstructural protein 1 coding region in murine and bovine coronaviruses

Higher-order RNA structures in the 5' untranslated regions (UTRs) of the mouse hepatitis coronavirus (MHV) and bovine coronavirus (BCoV), separate species in the betacoronavirus genus, appear to be largely conserved despite an ∼36% nucleotide sequence divergence. In a previous study, each of three 5'-end-proximal cis-acting stem-loop domains in the BCoV genome, I/II, III, and IV, yielded near-wild-type (wt) MHV phenotypes when used by reverse genetics to replace its counterpart in the MHV genome. Replacement with the BCoV 32-nucleotide (nt) inter-stem-loop fourth domain between stem-loops III and IV, however, required blind cell passaging for virus recovery. Here, we describe suppressor mutations within the transplanted BCoV 32-nt domain that along with appearance of potential base pairings identify an RNA-RNA interaction between this domain and a 32-nt region ∼200 nt downstream within the nonstructural protein 1 (Nsp1)-coding region. Mfold and phylogenetic covariation patterns among similarly grouped betacoronaviruses support this interaction, as does cotransplantation of the BCoV 5' UTR and its downstream base-pairing domain. Interestingly, cotransplantation of the BCoV 5' UTR and BCoV Nsp1 coding region directly yielded an MHV wt-like phenotype, which demonstrates a cognate interaction between these two BCoV regions, which in the MHV genome act in a fully interspecies-compliant manner. Surprisingly, the 30-nt inter-stem-loop domain in the MHV genome can be deleted and viral progeny, although debilitated, are still produced. These results together identify a previously undescribed long-range RNA-RNA interaction between the 5' UTR and Nsp1 coding region in MHV-like and BCoV-like betacoronaviruses that is cis acting for viral fitness but is not absolutely required for viral replication in cell culture.

Mouse hepatitis virus stem-loop 4 functions as a spacer element required to drive subgenomic RNA synthesis

The 5' 140 nucleotides of the mouse hepatitis virus (MHV) 5' untranslated region (5'UTR) are predicted to contain three secondary structures, stem-loop 1 (SL1), SL2, and SL4. SL1 and SL2 are required for subgenomic RNA synthesis. The current study focuses on SL4, which contains two base-paired regions, SL4a and SL4b. A series of reverse genetic experiments show that SL4a is not required to be base paired. Neither the structure, the sequence, nor the putative 8-amino-acid open reading frame (ORF) in SL4b is required for viral replication. Viruses containing separate deletions of SL4a and SL4b are viable. However, deletion of SL4 is lethal, and genomes carrying this deletion are defective in directing subgenomic RNA synthesis. Deletion of (131)ACA(133) just 3' to SL4 has a profound impact on viral replication. Viruses carrying the (131)ACA(133) deletion were heterogeneous in plaque size. We isolated three viruses with second-site mutations in the 5'UTR which compensated for decreased plaque sizes, delayed growth kinetics, and lower titers associated with the (131)ACA(133) deletion. The second-site mutations are predicted to change either the spacing between SL1 and SL2 or that between SL2 and SL4 or to destabilize the proximal portion of SL4a in our model. A mutant constructed by replacing SL4 with a shorter sequence-unrelated stem-loop was viable. These results suggest that the proposed SL4 in the MHV 5'UTR functions in part as a spacer element that orients SL1, SL2, and the transcriptional regulatory sequence (TRS), and this spacer function may play an important role in directing subgenomic RNA synthesis.

A 5'-proximal stem-loop structure of 5' untranslated region of porcine reproductive and respiratory syndrome virus genome is key for virus replication

Background

It has been well documented that the 5' untranslated region (5' UTR) of many positive-stranded RNA viruses contain key cis-acting regulatory sequences, as well as high-order structural elements. Little is known for such regulatory elements controlling porcine arterivirus replication. We investigated the roles of a conserved stem-loop 2 (SL2) that resides in the 5'UTR of the genome of a type II porcine reproductive and respiratory syndrome virus (PRRSV).

Results

We provided genetic evidences demonstrating that 1) the SL2 in type II PRRSV 5' UTR, N-SL2, could be structurally and functionally substituted by its counterpart in type I PRRSV, E-SL2; 2) the functionality of N-SL2 was dependent upon the G-C rich stem structure, while the ternary-loop size was irrelevant to RNA synthesis; 3) serial deletions showed that the stem integrity of N-SL2 was crucial for subgenomic mRNA synthesis; and 4) when extensive base-pairs in the stem region was deleted, an alternative N-SL2-like structure with different sequence was utilized for virus replication.

Conclusion

Taken together, we concluded that the phylogenetically conserved SL2 in the 5' UTR was crucial for PRRSV virus replication, subgenomic mRNA synthesis in particular.

An optimal cis-replication stem-loop IV in the 5' untranslated region of the mouse coronavirus genome extends 16 nucleotides into open reading frame 1

The 288-nucleotide (nt) 3' untranslated region (UTR) in the genome of the bovine coronavirus (BCoV) and 339-nt 3' UTR in the severe acute respiratory syndrome (SARS) coronavirus (SCoV) can each replace the 301-nt 3' UTR in the mouse hepatitis coronavirus (MHV) for virus replication, thus demonstrating common 3' cis-replication signals. Here, we show that replacing the 209-nt MHV 5' UTR with the ∼63%-sequence-identical 210-nt BCoV 5' UTR by reverse genetics does not yield viable virus, suggesting 5' end signals are more stringent or possibly are not strictly 5' UTR confined. To identify potential smaller, 5'-common signals, each of three stem-loop (SL) signaling domains and one inter-stem-loop domain from the BCoV 5' UTR was tested by replacing its counterpart in the MHV genome. The SLI/II domain (nucleotides 1 to 84) and SLIII domain (nucleotides 85 to 141) each immediately enabled near-wild-type (wt) MHV-like progeny, thus behaving similarly to comparable 5'-proximal regions of the SCoV 5' UTR as shown by others. The inter-stem-loop domain (nt 142 to 173 between SLs III and IV) enabled small plaques only after genetic adaptation. The SLIV domain (nt 174 to 210) required a 16-nt extension into BCoV open reading frame 1 (ORF1) for apparent stabilization of a longer BCoV SLIV (nt 174 to 226) and optimal virus replication. Surprisingly, pleiomorphic SLIV structures, including a terminal loop deletion, were found among debilitated progeny from intra-SLIV chimeras. The results show the inter-stem-loop domain to be a potential novel species-specific cis-replication element and that cis-acting SLIV in the viral genome extends into ORF1 in a manner that stabilizes its lower stem and is thus not 5' UTR confined.

The solution structure of coronaviral stem-loop 2 (SL2) reveals a canonical CUYG tetraloop fold

The transcription and replication of the severe acute respiratory syndrome (SARS) coronavirus (SARS‐CoV) is regulated by specific viral genome sequences within 5′‐ and 3′‐untranslated regions (5′‐UTR and 3′‐UTR). Here we report the solution structure of 5′‐UTR derived stem‐loop 2 (SL2) of SARS‐CoV determined by NMR spectroscopy. The highly conserved pentaloop of SL2 is stacked on 5‐bp stem and adopts a canonical CUYG tetraloop fold with the 3′ nucleotide (U51) flipped out of the stack. The significance of this structure in the context of a previous mutagenesis analysis of SL2 function in replication of the related group 2 coronavirus, mouse hepatitis virus, is discussed.

RNA-RNA and RNA-protein interactions in coronavirus replication and transcription

Coronavirus (CoV) RNA synthesis includes the replication of the viral genome, and the transcription of sgRNAs by a discontinuous mechanism. Both processes are regulated by RNA sequences such as the 5' and 3' untranslated regions (UTRs), and the transcription regulating sequences (TRSs) of the leader (TRS-L) and those preceding each gene (TRS-Bs). These distant RNA regulatory sequences interact with each other directly and probably through protein-RNA and protein-protein interactions involving viral and cellular proteins. By analogy to other plus-stranded RNA viruses, such as polioviruses, in which translation and replication switch involves a cellular factor (PCBP) and a viral protein (3CD) it is conceivable that in CoVs the switch between replication and transcription is also associated with the binding of proteins that are specifically recruited by the replication or transcription complexes. Complexes between RNA motifs such as TRS-L and the TRS-Bs located along the CoV genome are probably formed previously to the transcription start, and most likely promote template-switch of the nascent minus RNA to the TRS-L region. Many cellular proteins interacting with regulatory CoV RNA sequences are members of the heterogeneous nuclear ribonucleoprotein (hnRNP) family of RNA-binding proteins, involved in mRNA processing and transport, which shuttle between the nucleus and the cytoplasm. In the context of CoV RNA synthesis, these cellular ribonucleoproteins might also participate in RNA-protein complexes to bring into physical proximity TRS-L and distant TRS-B, as proposed for CoV discontinuous transcription. In this review, we summarize RNA-RNA and RNA-protein interactions that represent modest examples of complex quaternary RNA-protein structures required for the fine-tuning of virus replication. Design of chemically defined replication and transcription systems will help to clarify the nature and activity of these structures.

Group-specific structural features of the 5'-proximal sequences of coronavirus genomic RNAs

Global predictions of the secondary structure of coronavirus (CoV) 5′ untranslated regions and adjacent coding sequences revealed the presence of conserved structural elements. Stem loops (SL) 1, 2, 4, and 5 were predicted in all CoVs, while the core leader transcription-regulating sequence (L-TRS) forms SL3 in only some CoVs. SL5 in group I and II CoVs, with the exception of group IIa CoVs, is characterized by the presence of a large sequence insertion capable of forming hairpins with the conserved 5′-UUYCGU-3′ loop sequence. Structure probing confirmed the existence of these hairpins in the group I Human coronavirus-229E and the group II Severe acute respiratory syndrome coronavirus (SARS-CoV). In general, the pattern of the 5′ cis-acting elements is highly related to the lineage of CoVs, including features of the conserved hairpins in SL5. The function of these conserved hairpins as a putative packaging signal is discussed.