Characterization of an essential RNA secondary structure in the 3' untranslated region of the murine coronavirus genome

Genomic Evolution Strategy in SARS-CoV-2 Lineage B: Coevolution of Cis Elements

In the SARS-CoV-2 lineage, RNA elements essential for its viral life cycle, including genome replication and gene expression, have been identified. Still, the precise structures and functions of these RNA regions in coronaviruses remain poorly understood. This lack of knowledge points out the need for further research to better understand these crucial aspects of viral biology and, in time, prepare for future outbreaks. In this research, the in silico analysis of the cis RNA structures that act in the alpha-, beta-, gamma-, and deltacoronavirus genera has provided a detailed view of the presence and adaptation of the structures of these elements in coronaviruses. The results emphasize the importance of these cis elements in viral biology and their variability between different viral variants. Some coronavirus variants in some groups, depending on the cis element (stem-loop1 and -2; pseudoknot stem-loop1 and -2, and s2m), exhibited functional adaptation. Additionally, the conformation flexibility of the s2m element in the SARS variants was determined, suggesting a coevolution of this element in this viral group. The variability in secondary structures suggests genomic adaptations that may be related to replication processes, genetic regulation, as well as the specific pathogenicity of each variant. The results suggest that RNA structures in coronaviruses can adapt and evolve toward different viral variants, which has important implications for viral adaptation, pathogenicity, and future therapeutic strategies.

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.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

The coronavirus recombination pathway

Recombination is thought to be a mechanism that facilitates cross-species transmission in coronaviruses, thus acting as a driver of coronavirus spillover and emergence. Despite its significance, the mechanism of recombination is poorly understood, limiting our potential to estimate the risk of novel recombinant coronaviruses emerging in the future. As a tool for understanding recombination, here, we outline a framework of the recombination pathway for coronaviruses. We review existing literature on coronavirus recombination, including comparisons of naturally observed recombinant genomes as well as in vitro experiments, and place the findings into the recombination pathway framework. We highlight gaps in our understanding of coronavirus recombination illustrated by the framework and outline how further experimental research is critical for disentangling the molecular mechanism of recombination from external environmental pressures. Finally, we describe how an increased understanding of the mechanism of recombination can inform pandemic predictive intelligence, with a retrospective emphasis on SARS-CoV-2.

Deletion of the s2m RNA Structure in the Avian Coronavirus Infectious Bronchitis Virus and Human Astrovirus Results in Sequence Insertions

Coronaviruses infect a wide variety of host species, resulting in a range of diseases in both humans and animals. The coronavirus genome consists of a large positive-sense single-stranded molecule of RNA containing many RNA structures. One structure, denoted s2m and consisting of 41 nucleotides, is located within the 3' untranslated region (3' UTR) and is shared between some coronavirus species, including infectious bronchitis virus (IBV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2, as well as other pathogens, including human astrovirus. Using a reverse genetic system to generate recombinant viruses, we investigated the requirement of the s2m structure in the replication of IBV, a globally distributed economically important Gammacoronavirus that infects poultry causing respiratory disease. Deletion of three nucleotides predicted to destabilize the canonical structure of the s2m or the deletion of the nucleotides corresponding to s2m impacted viral replication in vitro. In vitro passaging of the recombinant IBV with the s2m sequence deleted resulted in a 36-nucleotide insertion in place of the deletion, which was identified to be composed of a duplication of flanking sequences. A similar result was observed following serial passage of human astrovirus with a deleted s2m sequence. RNA modeling indicated that deletion of the nucleotides corresponding to the s2m impacted other RNA structures present in the IBV 3' UTR. Our results indicated for both IBV and human astrovirus a preference for nucleotide occupation in the genome location corresponding to the s2m, which is independent of the specific s2m sequence. IMPORTANCE Coronaviruses infect many species, including humans and animals, with substantial effects on livestock, particularly with respect to poultry. The coronavirus RNA genome consists of structural elements involved in viral replication whose roles are poorly understood. We investigated the requirement of the RNA structural element s2m in the replication of the Gammacoronavirus infectious bronchitis virus, an economically important viral pathogen of poultry. Using reverse genetics to generate recombinant IBVs with either a disrupted or deleted s2m, we showed that the s2m is not required for viral replication in cell culture; however, replication is decreased in tracheal tissue, suggesting a role for the s2m in the natural host. Passaging of these viruses as well as human astrovirus lacking the s2m sequence demonstrated a preference for nucleotide occupation, independent of the s2m sequence. RNA modeling suggested deletion of the s2m may negatively impact other essential RNA structures.

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.

Computational Pipeline for Reference-Free Comparative Analysis of RNA 3D Structures Applied to SARS-CoV-2 UTR Models

RNA is a unique biomolecule that is involved in a variety of fundamental biological functions, all of which depend solely on its structure and dynamics. Since the experimental determination of crystal RNA structures is laborious, computational 3D structure prediction methods are experiencing an ongoing and thriving development. Such methods can lead to many models; thus, it is necessary to build comparisons and extract common structural motifs for further medical or biological studies. Here, we introduce a computational pipeline dedicated to reference-free high-throughput comparative analysis of 3D RNA structures. We show its application in the RNA-Puzzles challenge, in which five participating groups attempted to predict the three-dimensional structures of 5'- and 3'-untranslated regions (UTRs) of the SARS-CoV-2 genome. We report the results of this puzzle and discuss the structural motifs obtained from the analysis. All simulated models and tools incorporated into the pipeline are open to scientific and academic use.

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.

A Novel Motif in the 3′-UTR of PRRSV-2 Is Critical for Viral Multiplication and Contributes to Enhanced Replication Ability of Highly Pathogenic or L1 PRRSV

Highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) with enhanced replication capability emerged in China and has become dominant epidemic strain since 2006. Up to now, the replication-regulated genes of PRRSV have not been fully clarified. Here, by swapping the genes or elements between HP-PRRSV and classical PRRSV based on infectious clones, NSP1, NSP2, NSP7, NSP9 and 3′-UTR are found to contribute to the high replication efficiency of HP-PRRSV. Further study revealed that mutations at positions 117th or 119th in the 3′-UTR are significantly related to replication efficiency, and the nucleotide at position 120th is critical for viral rescue. The motif composed by 117–120th nucleotides was quite conservative within each lineage of PRRSV; mutations in the motif of HP-PRRSV and currently epidemic lineage 1 (L1) PRRSV showed higher synthesis ability of viral negative genomic RNA, suggesting that those mutations were beneficial for viral replication. RNA structure analysis revealed that this motif maybe involved into a pseudoknot in the 3′-UTR. The results discovered a novel motif, 117–120th nucleotide in the 3′-UTR, that is critical for replication of PRRSV-2, and mutations in the motif contribute to the enhanced replicative ability of HP-PRRSV or L1 PRRSV. Our findings will help to understand the molecular basis of PRRSV replication and find the potential factors resulting in an epidemic strain of PRRSV.

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

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

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.

Structure and regulation of coronavirus genomes: state-of-the-art and novel insights from SARS-CoV-2 studies

Coronaviruses (CoV) are positive-sense single-stranded RNA viruses, harboring the largest viral RNA genomes known to date. Apart from the primary sequence encoding for all the viral proteins needed for the generation of new viral particles, certain regions of CoV genomes are known to fold into stable structures, controlling several aspects of CoV life cycle, from the regulation of the discontinuous transcription of subgenomic mRNAs, to the packaging of the genome into new virions. Here we review the current knowledge on CoV RNA structures, discussing it in light of the most recent discoveries made possible by analyses of the SARS-CoV-2 genome.

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.

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.

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.

Gene Variations in Cis-Acting Elements between the Taiwan and Prototype Strains of Porcine Epidemic Diarrhea Virus Alter Viral Gene Expression

In 2013, the outbreak of porcine epidemic diarrhea (PED) in Taiwan caused serious economic losses. In this study, we examined whether the variations of the cis-acting elements between the porcine epidemic diarrhea virus (PEDV) Taiwan (TW) strain and the prototype strain CV777 alter gene expression. For this aim, we analyzed the variations of the cis-acting elements in the 5' and 3' untranslated regions (UTRs) between the PEDV TW, CV777, and other reference strains. We also determined the previously unidentified transcription regulatory sequence (TRS), a sequence motif required for coronavirus transcription, and found that a nucleotide deletion in the TW strain, in comparison with CV777 strain, immediately downstream of the leader core sequence alters the identity between the leader TRS and the body TRS. Functional analyses using coronavirus defective interfering (DI) RNA revealed that such variations in cis-acting elements for the TW strain compared with the CV777 strain have an influence on the efficiency of gene expression. The current data show for the first time the evolution of PEDV in terms of cis-acting elements and their effects on gene expression, and thus may contribute to our understanding of recent PED outbreaks worldwide.

Genomic characterization and pathogenicity of a porcine hemagglutinating encephalomyelitis virus strain isolated in China

Porcine hemagglutinating encephalomyelitis virus (PHEV) is a member of the genus betacoronavirus within the family coronaviridae, which invades the central nervous system (CNS) via peripheral nervous system and causes encephalomyelitis or vomiting and wasting disease (VWD) in sucking piglets. Up to now, although few complete nucleotide sequences of PHEV have been reported, they are not annotated. This study aimed to illuminate genome characterization, phylogenesis and pathogenicity of the PHEV/2008 strain. The full length of the PHEV/2008 strain genome was 30,684 bp, with a G + C content of 37.27%. The genome included at a minimum of 11 predicted open reading frames (ORFs) flanked by 5' and 3' untranslated regions (UTR) of 211 and 289 nucleotides. The replicase polyproteins pp1a and pp1ab, which had 4382 and 7094 amino acid residues, respectively, were predicted to be cleaved into 16 subunits by two viral proteinases. Phylogenetic analysis based on the complete genome sequence revealed that PHEV/2008 strain was genetically different from other known PHEV types, which represented a novel genotype (GI-1). In addition, we found that PHEV/2008 was neurotropic and highly pathogenic to 4-week-old BALB/c mice. Taken together, this is the first detailed annotated, complete genomic sequence of a new genotype PHEV strain in China.

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.

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).

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.

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.

Unmasking the information encoded as structural motifs of viral RNA genomes: a potential antiviral target

RNA viruses show enormous capacity to evolve and adapt to new cellular and molecular contexts, a consequence of mutations arising from errors made by viral RNA-dependent RNA polymerase during replication. Sequence variation must occur, however, without compromising functions essential for the completion of the viral cycle. RNA viruses are safeguarded in this respect by their genome carrying conserved information that does not code only for proteins but also for the formation of structurally conserved RNA domains that directly perform these critical functions. Functional RNA domains can interact with other regions of the viral genome and/or proteins to direct viral translation, replication and encapsidation. They are therefore potential targets for novel therapeutic strategies. This review summarises our knowledge of the functional RNA domains of human RNA viruses and examines the achievements made in the design of antiviral compounds that interfere with their folding and therefore their function.

Functional analysis of the stem loop S3 and S4 structures in the coronavirus 3'UTR

We designed a series of mutations to separately destabilize two helical stems (designated S3 and S4) predicted by a covariation-based model of the coronavirus 3′UTR (Zust et al., 2008). Mouse hepatitis virus genomes containing three or four nucleotide mutations that destabilize either S3 or S4 were viable, whereas genomes carrying these mutations in both S3 and S4 were not viable. A genome carrying these mutations in S3 and S4 plus compensatory mutations restoring base-pairing yielded a virus with wild type phenotype. Larger mutations which completely disrupt S3 or S4 generated various phenotypes. Mutations opening up S3 were lethal. Disruptions of S4 generated both viable and lethal mutants. Genomes carrying the original mutations in S3 or S4 plus compensatory mutations restoring base pairing were viable and had robust growth phenotypes. These results support the Zust model for the coronavirus 3′UTR and suggest that the S3 stem is required for virus viability.

•

Genetic studies of MHV 3′UTR RNA secondary structure support the Zust et al. model.

•

The 3′UTR S3 helical stem is required for virus viability.

•

The 3′UTR S4 helical stem is not required for virus viability.

The first complete genome sequences of clinical isolates of human coronavirus 229E

Human coronavirus 229E has been identified in the mid-1960s, yet still only one full-genome sequence is available. This full-length sequence has been determined from the cDNA-clone Inf-1 that is based on the lab-adapted strain VR-740. Lab-adaptation might have resulted in genomic changes, due to insufficient pressure to maintain gene integrity of non-essential genes. We present here the first full-length genome sequence of two clinical isolates. Each encoded gene was compared to Inf-1. In general, little sequence changes were noted, most could be attributed to genetic drift, since the clinical isolates originate from 2009 to 2010 and VR740 from 1962. Hot spots of substitutions were situated in the S1 region of the Spike, the nucleocapsid gene, and the non-structural protein 3 gene, whereas several deletions were detected in the 3′UTR. Most notable was the difference in genome organization: instead of an ORF4A and ORF4B, an intact ORF4 was present in clinical isolates.

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.

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.

Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria

Bats are reservoirs for emerging zoonotic viruses that can have a profound impact on human and animal health, including lyssaviruses, filoviruses, paramyxoviruses, and severe acute respiratory syndrome coronaviruses (SARS-CoVs). In the course of a project focused on pathogen discovery in contexts where human-bat contact might facilitate more efficient interspecies transmission of viruses, we surveyed gastrointestinal tissue obtained from bats collected in caves in Nigeria that are frequented by humans. Coronavirus consensus PCR and unbiased high-throughput pyrosequencing revealed the presence of coronavirus sequences related to those of SARS-CoV in a Commerson's leaf-nosed bat (Hipposideros commersoni). Additional genomic sequencing indicated that this virus, unlike subgroup 2b CoVs, which includes SARS-CoV, is unique, comprising three overlapping open reading frames between the M and N genes and two conserved stem-loop II motifs. Phylogenetic analyses in conjunction with these features suggest that this virus represents a new subgroup within group 2 CoVs.

RNA Higher-Order Structures Within the Coronavirus 5′ and 3′ Untranslated Regions and Their Roles in Viral Replication

The 5′ and 3′ untranslated regions (UTRs) of all coronaviruses contain RNA higher-order structures which play essential roles in viral transcription and replication. In this chapter we present our current knowledge of how those cis-acting elements were defined and their functional roles in viral transcription and replication. Cellular proteins which have been shown binding to those cis-acting elements and potentially support the RNA discontinuous synthesis model are also discussed. A conserved RNA structure model for the 5′ and 3′ UTRs of group 2 coronaviruses is presented with the known cellular protein binding sites.

RNA conformational changes in the life cycles of RNA viruses, viroids, and virus-associated RNAs

The rugged nature of the RNA structural free energy landscape allows cellular RNAs to respond to environmental conditions or fluctuating levels of effector molecules by undergoing dynamic conformational changes that switch on or off activities such as catalysis, transcription or translation. Infectious RNAs must also temporally control incompatible activities and rapidly complete their life cycle before being targeted by cellular defenses. Viral genomic RNAs must switch between translation and replication, and untranslated subviral RNAs must control other activities such as RNA editing or self-cleavage. Unlike well characterized riboswitches in cellular RNAs, the control of infectious RNA activities by altering the configuration of functional RNA domains has only recently been recognized. In this review, we will present some of these molecular rearrangements found in RNA viruses, viroids and virus-associated RNAs, relating how these dynamic regions were discovered, the activities that might be regulated, and what factors or conditions might cause a switch between conformations.

Bovine coronavirus nonstructural protein 1 (p28) is an RNA binding protein that binds terminal genomic cis-replication elements

Nonstructural protein 1 (nsp1), a 28-kDa protein in the bovine coronavirus (BCoV) and closely related mouse hepatitis coronavirus, is the first protein cleaved from the open reading frame 1 (ORF 1) polyprotein product of genome translation. Recently, a 30-nucleotide (nt) cis-replication stem-loop VI (SLVI) has been mapped at nt 101 to 130 within a 288-nt 5'-terminal segment of the 738-nt nsp1 cistron in a BCoV defective interfering (DI) RNA. Since a similar nsp1 coding region appears in all characterized groups 1 and 2 coronavirus DI RNAs and must be translated in cis for BCoV DI RNA replication, we hypothesized that nsp1 might regulate ORF 1 expression by binding this intra-nsp1 cistronic element. Here, we (i) establish by mutation analysis that the 72-nt intracistronic SLV immediately upstream of SLVI is also a DI RNA cis-replication signal, (ii) show by gel shift and UV-cross-linking analyses that cellular proteins of approximately 60 and 100 kDa, but not viral proteins, bind SLV and SLVI, (SLV-VI) and (iii) demonstrate by gel shift analysis that nsp1 purified from Escherichia coli does not bind SLV-VI but does bind three 5' untranslated region (UTR)- and one 3' UTR-located cis-replication SLs. Notably, nsp1 specifically binds SLIII and its flanking sequences in the 5' UTR with approximately 2.5 muM affinity. Additionally, under conditions enabling expression of nsp1 from DI RNA-encoded subgenomic mRNA, DI RNA levels were greatly reduced, but there was only a slight transient reduction in viral RNA levels. These results together indicate that nsp1 is an RNA-binding protein that may function to regulate viral genome translation or replication but not by binding SLV-VI within its own coding region.

Genetic interactions between an essential 3' cis-acting RNA pseudoknot, replicase gene products, and the extreme 3' end of the mouse coronavirus genome

The upstream end of the 3' untranslated region (UTR) of the mouse hepatitis virus genome contains two essential and overlapping RNA secondary structures, a bulged stem-loop and a pseudoknot, which have been proposed to be elements of a molecular switch that is critical for viral RNA synthesis. It has previously been shown that a particular six-base insertion in loop 1 of the pseudoknot is extremely deleterious to the virus. We have now isolated multiple independent second-site revertants of the loop 1 insertion mutant, and we used reverse-genetics methods to confirm the identities of suppressor mutations that could compensate for the original insertion. The suppressors were localized to two separate regions of the genome. Members of one class of suppressor were mapped to the portions of gene 1 that encode nsp8 and nsp9, thereby providing the first evidence for specific interactions between coronavirus replicase gene products and a cis-acting genomic RNA element. The second class of suppressor was mapped to the extreme 3' end of the genome, a result which pointed to the existence of a direct base-pairing interaction between loop 1 of the pseudoknot and the genomic terminus. The latter finding was strongly supported by phylogenetic evidence and by the construction of a deletion mutant that reduced the 3' UTR to its minimal essential elements. Taken together, the interactions revealed by the two classes of suppressors suggest a model for the initiation of coronavirus negative-strand RNA synthesis.

An RNA pseudoknot in the 3' end of the arterivirus genome has a critical role in regulating viral RNA synthesis

In the life cycle of plus-strand RNA viruses, the genome initially serves as the template for both translation of the viral replicase gene and synthesis of minus-strand RNA and is ultimately packaged into progeny virions. These various processes must be properly balanced to ensure efficient viral proliferation. To achieve this, higher-order RNA structures near the termini of a variety of RNA virus genomes are thought to play a key role in regulating the specificity and efficiency of viral RNA synthesis. In this study, we have analyzed the signals for minus-strand RNA synthesis in the prototype of the arterivirus family, equine arteritis virus (EAV). Using site-directed mutagenesis and an EAV reverse genetics system, we have demonstrated that a stem-loop structure near the 3' terminus of the EAV genome is required for RNA synthesis. We have also obtained evidence for an essential pseudoknot interaction between the loop region of this stem-loop structure and an upstream hairpin residing in the gene encoding the nucleocapsid protein. We propose that the formation of this pseudoknot interaction may constitute a molecular switch that could regulate the specificity or timing of viral RNA synthesis. This hypothesis is supported by the fact that phylogenetic analysis predicted the formation of similar pseudoknot interactions near the 3' end of all known arterivirus genomes, suggesting that this interaction has been conserved in evolution.

Genomic characterization of equine coronavirus

The complete genome sequence of the first equine coronavirus (ECoV) isolate, NC99 strain was accomplished by directly sequencing 11 overlapping fragments which were RT–PCR amplified from viral RNA. The ECoV genome is 30,992 nucleotides in length, excluding the polyA tail. Analysis of the sequence identified 11 open reading frames which encode two replicase polyproteins, five structural proteins (hemagglutinin esterase, spike, envelope, membrane, and nucleocapsid) and four accessory proteins (NS2, p4.7, p12.7, and I). The two replicase polyproteins are predicted to be proteolytically processed by three virus-encoded proteases into 16 non-structural proteins (nsp1–16). The ECoV nsp3 protein had considerable amino acid deletions and insertions compared to the nsp3 proteins of bovine coronavirus, human coronavirus OC43, and porcine hemagglutinating encephalomyelitis virus, three group 2 coronaviruses phylogenetically most closely related to ECoV. The structure of subgenomic mRNAs was analyzed by Northern blot analysis and sequencing of the leader–body junction in each sg mRNA.

An RNA stem-loop within the bovine coronavirus nsp1 coding region is a cis-acting element in defective interfering RNA replication

Higher-order cis-acting RNA replication structures have been identified in the 3'- and 5'-terminal untranslated regions (UTRs) of a bovine coronavirus (BCoV) defective interfering (DI) RNA. The UTRs are identical to those in the viral genome, since the 2.2-kb DI RNA is composed of only the two ends of the genome fused between an internal site within the 738-nucleotide (nt) 5'-most coding region (the nsp1, or p28, coding region) and a site just 4 nt upstream of the 3'-most open reading frame (ORF) (the N gene). The joined ends of the viral genome in the DI RNA create a single continuous 1,635-nt ORF, 288 nt of which come from the 738-nt nsp1 coding region. Here, we have analyzed features of the 5'-terminal 288-nt portion of the nsp1 coding region within the continuous ORF that are required for DI RNA replication. We observed that (i) the 5'-terminal 186 nt of the nsp1 coding region are necessary and sufficient for DI RNA replication, (ii) two Mfold-predicted stem-loops within the 186-nt sequence, named SLV (nt 239 to 310) and SLVI (nt 311 to 340), are supported by RNase structure probing and by nucleotide covariation among closely related group 2 coronaviruses, and (iii) SLVI is a required higher-order structure for DI RNA replication based on mutation analyses. The function of SLV has not been evaluated. We conclude that SLVI within the BCoV nsp1 coding region is a higher-order cis-replication element for DI RNA and postulate that it functions similarly in the viral genome.

A U-turn motif-containing stem-loop in the coronavirus 5' untranslated region plays a functional role in replication

The 5' untranslated region (UTR) of the mouse hepatitis virus (MHV) genome contains cis-acting sequences necessary for transcription and replication. A consensus secondary structural model of the 5' 140 nucleotides of the 5' UTRs of nine coronaviruses (CoVs) derived from all three major CoV groups is presented and characterized by three major stem-loops, SL1, SL2, and SL4. NMR spectroscopy provides structural support for SL1 and SL2 in three group 2 CoVs, including MHV, BCoV, and HCoV-OC43. SL2 is conserved in all CoVs, typically containing a pentaloop (C47-U48-U49-G50-U51 in MHV) stacked on a 5 base-pair stem, with some sequences containing an additional U 3' to U51; SL2 therefore possesses sequence features consistent with a U-turn-like conformation. The imino protons of U48 in the wild-type RNA, and G48 in the U48G SL2 mutant RNA, are significantly protected from exchange with solvent, consistent with a hydrogen bonding interaction critical to the hairpin loop architecture. SL2 is required for MHV replication; MHV genomes containing point substitutions predicted to perturb the SL2 structure (U48C, U48A) were not viable, while those that maintain the structure (U48G and U49A) were viable. The U48C MHV mutant supports both positive- and negative-sense genome-sized RNA synthesis, but fails to direct the synthesis of positive- or negative-sense subgenomic RNAs. These data support the existence of the SL2 in our models, and further suggest a critical role in coronavirus replication.

A hypervariable region within the 3' cis-acting element of the murine coronavirus genome is nonessential for RNA synthesis but affects pathogenesis

The 3' cis-acting element for mouse hepatitis virus (MHV) RNA synthesis resides entirely within the 301-nucleotide 3' untranslated region (3' UTR) of the viral genome and consists of three regions. Encompassing the upstream end of the 3' UTR are a bulged stem-loop and an overlapping RNA pseudoknot, both of which are essential to MHV and common to all group 2 coronaviruses. At the downstream end of the genome is the minimal signal for initiation of negative-strand RNA synthesis. Between these two ends is a hypervariable region (HVR) that is only poorly conserved between MHV and other group 2 coronaviruses. Paradoxically, buried within the HVR is an octanucleotide motif (oct), 5'-GGAAGAGC-3', which is almost universally conserved in coronaviruses and is therefore assumed to have a critical biological function. We conducted an extensive mutational analysis of the HVR. Surprisingly, this region tolerated numerous deletions, rearrangements, and point mutations. Most striking, a mutant deleted of the entire HVR was only minimally impaired in tissue culture relative to the wild type. By contrast, the HVR deletion mutant was highly attenuated in mice, causing no signs of clinical disease and minimal weight loss compared to wild-type virus. Correspondingly, replication of the HVR deletion mutant in the brains of mice was greatly reduced compared to that of the wild type. Our results show that neither the HVR nor oct is essential for the basic mechanism of MHV RNA synthesis in tissue culture. However, the HVR appears to play a significant role in viral pathogenesis.

The molecular biology of coronaviruses

Coronaviruses are large, enveloped RNA viruses of both medical and veterinary importance. Interest in this viral family has intensified in the past few years as a result of the identification of a newly emerged coronavirus as the causative agent of severe acute respiratory syndrome (SARS). At the molecular level, coronaviruses employ a variety of unusual strategies to accomplish a complex program of gene expression. Coronavirus replication entails ribosome frameshifting during genome translation, the synthesis of both genomic and multiple subgenomic RNA species, and the assembly of progeny virions by a pathway that is unique among enveloped RNA viruses. Progress in the investigation of these processes has been enhanced by the development of reverse genetic systems, an advance that was heretofore obstructed by the enormous size of the coronavirus genome. This review summarizes both classical and contemporary discoveries in the study of the molecular biology of these infectious agents, with particular emphasis on the nature and recognition of viral receptors, viral RNA synthesis, and the molecular interactions governing virion assembly.

Putative cis-acting stem-loops in the 5' untranslated region of the severe acute respiratory syndrome coronavirus can substitute for their mouse hepatitis virus counterparts

Consensus covariation-based secondary structural models for the 5' 140 nucleotides of the 5' untranslated regions (5'UTRs) from mouse hepatitis virus (MHV) and severe acute respiratory syndrome coronavirus (SCoV) were developed and predicted three major helical stem-loop structures, designated stem-loop 1 (SL1), SL2, and SL4. The SCoV 5'UTR was predicted to contain a fourth stem-loop, named SL3, in which the leader transcriptional regulatory sequence (TRS) is folded into a hairpin loop. cDNAs corresponding to MHV/SCoV chimeric genomes were constructed by replacing the complete MHV 5'UTR with the corresponding SCoV sequence and by separately replacing MHV 5'UTR putative SL1, putative SL2, TRS, and putative SL4 with the corresponding SCoV sequences. Chimeric genomes were transcribed in vitro, and viruses were recovered after electroporation into permissive cells. Genomes in which the MHV 5'UTR SL1, SL2, and SL4 were individually replaced by their SCoV counterparts were viable. Chimeras containing the complete SCoV 5'UTR or the predicted SCoV SL3 were not viable. A chimera containing the SCoV 5'UTR in which the SCoV TRS was replaced with the MHV TRS was also not viable. The chimera containing the entire SCoV 5'UTR failed to direct the synthesis of any virus-specific RNA. Replacing the SCoV TRS with the MHV TRS in the MHV/5'UTR SCoV chimera permitted the synthesis of minus-sense genome-sized RNA but did not support the production of positive- or minus-sense subgenomic RNA7. A similar phenotype was obtained with the MHV/SCoV SL3 chimera. These results suggest a role for the TRS in the replication of minus-sense genomic RNA in addition to its known function in subgenomic RNA synthesis.

RNA signals in the 3' terminus of the genome of Equine arteritis virus are required for viral RNA synthesis

RNA virus genomes contain cis-acting sequences and structural elements involved in virus replication. Both full-length and subgenomic negative-strand RNA synthesis are initiated at the 3' terminus of the positive-strand genomic RNA of Equine arteritis virus (EAV). To investigate the molecular mechanism of EAV RNA synthesis, the RNA secondary structure of the 3'-proximal region of the genome was analysed by chemical and enzymic probing. Based on the RNA secondary structure model derived from this analysis, several deletions were engineered in a full-length cDNA copy of the viral genome. Two RNA domains were identified that are essential for virus replication and most likely play a key role in viral RNA synthesis. The first domain, located directly upstream of the 3' untranslated region (UTR) (nt 12610-12654 of the genome), is mainly single-stranded but contains one small stem-loop structure. The second domain is located within the 3' UTR (nt 12661-12690) and folds into a prominent stem-loop structure with a large loop region. The location of this stem-loop structure near the 3' terminus of the genome suggests that it may act as a recognition signal during the initiation of minus-strand RNA synthesis.

Luxury at a cost? Recombinant mouse hepatitis viruses expressing the accessory hemagglutinin esterase protein display reduced fitness in vitro

Group 2 coronaviruses encode an accessory envelope glycoprotein species, the hemagglutinin esterase (HE), which possesses sialate-O-acetylesterase activity and which, presumably, promotes virus spread and entry in vivo by facilitating reversible virion attachment to O-acetylated sialic acids. While HE may provide a strong selective advantage during natural infection, many laboratory strains of mouse hepatitis virus (MHV) fail to produce the protein. Apparently, their HE genes were inactivated during cell culture adaptation. For this report, we have studied the molecular basis of this phenomenon. By using targeted RNA recombination, we generated isogenic recombinant MHVs which differ exclusively in their expression of HE and produce either the wild-type protein (HE+), an enzymatically inactive HE protein (HE0), or no HE at all. HE expression or the lack thereof did not lead to gross differences in in vitro growth properties. Yet the expression of HE was rapidly lost during serial cell culture passaging. Competition experiments with mixed infections revealed that this was not due to the enzymatic activity: MHVs expressing HE+ or HE0 propagated with equal efficiencies. During the propagation of recombinant MHV-HE+, two types of spontaneous mutants accumulated. One produced an anchorless HE, while the other had a Gly-to-Trp substitution at the predicted C-terminal residue of the HE signal peptide. Neither mutant incorporated HE into virion particles, suggesting that wild-type HE reduces the in vitro propagation efficiency, either at the assembly stage or at a postassembly level. Our findings demonstrate that the expression of "luxury" proteins may come at a fitness penalty. Apparently, under natural conditions the costs of maintaining HE are outweighed by the benefits.

Effect of mutations in the mouse hepatitis virus 3'(+)42 protein binding element on RNA replication

The mouse hepatitis virus (MHV) genome's 3' untranslated region contains cis-acting sequences necessary for replication. Studies of MHV and other coronaviruses have indicated a role for RNA secondary and tertiary elements in replication. Previous work in our laboratory has identified four proteins which form ribonucleoprotein complexes with the 3'-terminal 42 nucleotides [3'(+)42] of the MHV genome. Defective interfering (DI) RNA replication assays have demonstrated a role for the 3'(+)42 host protein binding element in the MHV life cycle. Using gel mobility shift RNase T1 protection assays and secondary structure modeling, we have characterized a possible role for RNA secondary structure in host protein binding to the 3'-terminal 42-nucleotide element. Additionally we have identified a role for the 3'-terminal 42-nucleotide host protein binding element in RNA replication and transcription using DI RNA replication assays and targeted recombination and by directly constructing mutants in this protein binding element using a recently described MHV reverse genetic system. DI RNA replication assays demonstrated that mutations in the 3'(+)42 host protein binding element had a deleterious effect on the accumulation of DI RNA. When the identical mutations were directly inserted into the MHV genome, most mutant genomes were viable but formed smaller plaques than the wild-type parent virus. One mutant was not viable. This mutant directed the synthesis of genome-sized negative-sense RNA approximately as efficiently as the wild type did but had a defect in subgenomic mRNA synthesis. These results point to a potential role for sequences at the extreme 3' end of the MHV genome in subgenomic RNA synthesis.

Coronaviruses as vectors: stability of foreign gene expression

Coronaviruses are enveloped, positive-stranded RNA viruses considered to be promising vectors for vaccine development, as (i) genes can be deleted, resulting in attenuated viruses; (ii) their tropism can be modified by manipulation of their spike protein; and (iii) heterologous genes can be expressed by simply inserting them with appropriate coronaviral transcription signals into the genome. For any live vector, genetic stability is an essential requirement. However, little is known about the genetic stability of recombinant coronaviruses expressing foreign genes. In this study, the Renilla and the firefly luciferase genes were systematically analyzed for their stability after insertion at various genomic positions in the group 1 coronavirus feline infectious peritonitis virus and in the group 2 coronavirus mouse hepatitis virus. It appeared that the two genes exhibit intrinsic differences, the Renilla gene consistently being maintained more stably than the firefly gene. This difference was not caused by genome size restrictions, by different effects of the encoded proteins, or by different consequences of the synthesis of the additional subgenomic mRNAs. The loss of expression of the firefly luciferase was found to result from various, often large deletions of the gene, probably due to RNA recombination. The extent of this process appeared to depend strongly on the coronaviral genomic background, the luciferase gene being much more stable in the feline than in the mouse coronavirus genome. It also depended significantly on the particular genomic location at which the gene was inserted. The data indicate that foreign sequences are more stably maintained when replacing nonessential coronaviral genes.

Inhibition, escape, and attenuated growth of severe acute respiratory syndrome coronavirus treated with antisense morpholino oligomers

The recently emerged severe acute respiratory syndrome coronavirus (SARS-CoV) is a potent pathogen of humans and is capable of rapid global spread. Peptide-conjugated antisense morpholino oligomers (P-PMO) were designed to bind by base pairing to specific sequences in the SARS-CoV (Tor2 strain) genome. The P-PMO were tested for their capacity to inhibit production of infectious virus as well as to probe the function of conserved viral RNA motifs and secondary structures. Several virus-targeted P-PMO and a random-sequence control P-PMO showed low inhibitory activity against SARS coronavirus. Certain other virus-targeted P-PMO reduced virus-induced cytopathology and cell-to-cell spread as a consequence of decreasing viral amplification. Active P-PMO were effective when administered at any time prior to peak viral synthesis and exerted sustained antiviral effects while present in culture medium. P-PMO showed low nonspecific inhibitory activity against translation of nontargeted RNA or growth of the arenavirus lymphocytic choriomeningitis virus. Two P-PMO targeting the viral transcription-regulatory sequence (TRS) region in the 5' untranslated region were the most effective inhibitors tested. After several viral passages in the presence of a TRS-targeted P-PMO, partially drug-resistant SARS-CoV mutants arose which contained three contiguous base point mutations at the binding site of a TRS-targeted P-PMO. Those partially resistant viruses grew more slowly and formed smaller plaques than wild-type SARS-CoV. These results suggest PMO compounds have powerful therapeutic and investigative potential toward coronavirus infection.

Spike protein assembly into the coronavirion: exploring the limits of its sequence requirements

The coronavirus spike (S) protein, required for receptor binding and membrane fusion, is incorporated into the assembling virion by interactions with the viral membrane (M) protein. Earlier we showed that the ectodomain of the S protein is not involved in this process. Here we further defined the requirements of the S protein for virion incorporation. We show that the cytoplasmic domain, not the transmembrane domain, determines the association with the M protein and suffices to effect the incorporation into viral particles of chimeric spikes as well as of foreign viral glycoproteins. The essential sequence was mapped to the membrane-proximal region of the cytoplasmic domain, which is also known to be of critical importance for the fusion function of the S protein. Consistently, only short C-terminal truncations of the S protein were tolerated when introduced into the virus by targeted recombination. The important role of the about 38-residues cytoplasmic domain in the assembly of and membrane fusion by this approximately 1300 amino acids long protein is discussed.

Coronavirus genome structure and replication

In addition to the SARS coronavirus (treated separately elsewhere in this volume), the complete genome sequences of six species in the coronavirus genus of the coronavirus family [avian infectious bronchitis virus-Beaudette strain (IBV-Beaudette), bovine coronavirus-ENT strain (BCoV-ENT), human coronavirus-229E strain (HCoV-229E), murine hepatitis virus-A59 strain (MHV-A59), porcine transmissible gastroenteritis-Purdue 115 strain (TGEV-Purdue 115), and porcine epidemic diarrhea virus-CV777 strain (PEDV-CV777)] have now been reported. Their lengths range from 27,317 nt for HCoV-229E to 31,357 nt for the murine hepatitis virus-A59, establishing the coronavirus genome as the largest known among RNA viruses. The basic organization of the coronavirus genome is shared with other members of the Nidovirus order (the torovirus genus, also in the family Coronaviridae, and members of the family Arteriviridae) in that the nonstructural proteins involved in proteolytic processing, genome replication, and subgenomic mRNA synthesis (transcription) (an estimated 14–16 end products for coronaviruses) are encoded within the 5′-proximal two-thirds of the genome on gene 1 and the (mostly) structural proteins are encoded within the 3′-proximal one-third of the genome (8–9 genes for coronaviruses). Genes for the major structural proteins in all coronaviruses occur in the 5′ to 3′ order as S, E, M, and N. The precise strategy used by coronaviruses for genome replication is not yet known, but many features have been established. This chapter focuses on some of the known features and presents some current questions regarding genome replication strategy, the cis-acting elements necessary for genome replication [as inferred from defective interfering (DI) RNA molecules], the minimum sequence requirements for autonomous replication of an RNA replicon, and the importance of gene order in genome replication.

Viral and cellular proteins involved in coronavirus replication

As the largest RNA virus, coronavirus replication employs complex mechanisms and involves various viral and cellular proteins. The first open reading frame of the coronavirus genome encodes a large polyprotein, which is processed into a number of viral proteins required for viral replication directly or indirectly. These proteins include the RNA-dependent RNA polymerase (RdRp), RNA helicase, proteases, metal-binding proteins, and a number of other proteins of unknown function. Genetic studies suggest that most of these proteins are involved in viral RNA replication. In addition to viral proteins, several cellular proteins, such as heterogeneous nuclear ribonucleoprotein (hnRNP) A1, polypyrimidine-tract-binding (PTB) protein, poly(A)-binding protein (PABP), and mitochondrial aconitase (m-aconitase), have been identified to interact with the critical cis-acting elements of coronavirus replication. Like many other RNA viruses, coronavirus may subvert these cellular proteins from cellular RNA processing or translation machineries to play a role in viral replication.