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

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.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

Generation and functional analysis of defective viral genomes during SARS-CoV-2 infection.. [PMID: 36172120 PMCID: PMC9516852 DOI: 10.1101/2022.09.22.509123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from RNA-seq datasets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hotspots were identified for DVG recombination and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single cell RNA-seq analysis indicated the IFN stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the NGS dataset from a published cohort study and observed significantly higher DVG amount and frequency in symptomatic patients than that in asymptomatic patients. Finally, we observed unusually high DVG frequency in one immunosuppressive patient up to 140 days after admitted to hospital due to COVID-19, first-time suggesting an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection.

Double-stranded RNA drives SARS-CoV-2 nucleocapsid protein to undergo phase separation at specific temperatures

Nucleocapsid protein (N-protein) is required for multiple steps in betacoronaviruses replication. SARS-CoV-2-N-protein condenses with specific viral RNAs at particular temperatures making it a powerful model for deciphering RNA sequence specificity in condensates. We identify two separate and distinct double-stranded, RNA motifs (dsRNA stickers) that promote N-protein condensation. These dsRNA stickers are separately recognized by N-protein's two RNA binding domains (RBDs). RBD1 prefers structured RNA with sequences like the transcription-regulatory sequence (TRS). RBD2 prefers long stretches of dsRNA, independent of sequence. Thus, the two N-protein RBDs interact with distinct dsRNA stickers, and these interactions impart specific droplet physical properties that could support varied viral functions. Specifically, we find that addition of dsRNA lowers the condensation temperature dependent on RBD2 interactions and tunes translational repression. In contrast RBD1 sites are sequences critical for sub-genomic (sg) RNA generation and promote gRNA compression. The density of RBD1 binding motifs in proximity to TRS-L/B sequences is associated with levels of sub-genomic RNA generation. The switch to packaging is likely mediated by RBD1 interactions which generate particles that recapitulate the packaging unit of the virion. Thus, SARS-CoV-2 can achieve biochemical complexity, performing multiple functions in the same cytoplasm, with minimal protein components based on utilizing multiple distinct RNA motifs that control N-protein interactions.

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.

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.

Double-stranded RNA drives SARS-CoV-2 nucleocapsid protein to undergo phase separation at specific temperatures

Betacoronavirus SARS-CoV-2 infections caused the global Covid-19 pandemic. The nucleocapsid protein (N-protein) is required for multiple steps in the betacoronavirus replication cycle. SARS-CoV-2-N-protein is known to undergo liquid-liquid phase separation (LLPS) with specific RNAs at particular temperatures to form condensates. We show that N-protein recognizes at least two separate and distinct RNA motifs, both of which require double-stranded RNA (dsRNA) for LLPS. These motifs are separately recognized by N-protein's two RNA binding domains (RBDs). Addition of dsRNA accelerates and modifies N-protein LLPS in vitro and in cells and controls the temperature condensates form. The abundance of dsRNA tunes N-protein-mediated translational repression and may confer a switch from translation to genome packaging. Thus, N-protein's two RBDs interact with separate dsRNA motifs, and these interactions impart distinct droplet properties that can support multiple viral functions. These experiments demonstrate a paradigm of how RNA structure can control the properties of biomolecular condensates.

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.

In silico identification of conserved cis-acting RNA elements in the SARS-CoV-2 genome

Aim:

The aim of this study was to computationally predict conserved RNA sequences and structures known as cis-acting RNA elements (CREs) in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome.

Materials & methods:

Bioinformatics tools were used to analyze and predict CREs by obtaining viral sequences from available databases.

Results:

Computational analysis revealed the presence of RNA stem-loop structures within the 3′ end of the ORF1ab region analogous to previously identified SARS-CoV genomic packaging signals. Alignment-based RNA secondary structure predictions of the 5′ end of the SARS-CoV-2 genome also identified conserved CREs.

Conclusion:

These CREs may be potential vaccine and/or antiviral therapeutic targets; however, further studies are warranted to confirm their roles in the SARS-CoV-2 life cycle.

The Antiviral and Antitumor Effects of Defective Interfering Particles/Genomes and Their Mechanisms

Defective interfering particles (DIPs), derived naturally from viral particles, are not able to replicate on their own. Several studies indicate that DIPs exert antiviral effects via multiple mechanisms. DIPs are able to activate immune responses and suppress virus replication cycles, such as competing for viral replication products, impeding the packaging, release and invasion of viruses. Other studies show that DIPs can be used as a vaccine against viral infection. Moreover, DIPs/DI genomes display antitumor effects by inducing tumor cell apoptosis and promoting dendritic cell maturation. With genetic modified techniques, it is possible to improve its safety against both viruses and tumors. In this review, a comprehensive discussion on the effects exerted by DIPs is provided. We further highlight the clinical significance of DIPs and propose that DIPs can open up a new platform for antiviral and antitumor therapies.

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.

Interplay between the Poly(A) Tail, Poly(A)-Binding Protein, and Coronavirus Nucleocapsid Protein Regulates Gene Expression of Coronavirus and the Host Cell

In the present study, we investigated the roles of interactions among the poly(A) tail, coronavirus nucleocapsid (N) protein, and poly(A)-binding protein (PABP) in the regulation of coronavirus gene expression. Through dissociation constant (K_d ) comparison, we found that the coronavirus N protein can bind to the poly(A) tail with high affinity, establishing N protein as a PABP. A subsequent analysis with UV cross-linking and immunoprecipitation revealed that the N protein is able to bind to the poly(A) tail in infected cells. Further examination demonstrated that poly(A) tail binding by the N protein negatively regulates translation of coronaviral RNA and host mRNA both in vitro and in cells. Although the N protein can interact with PABP and eukaryotic initiation factor 4G (eIF4G), the poor interaction efficiency between the poly(A)-bound N protein and eIF4E may explain the observed decreased translation efficiency. In addition to interaction with translation factor eIF4G, the N protein is able to interact with coronavirus nonstructural protein 9 (nsp9), a replicase protein required for replication. The study demonstrates interactions among the poly(A) tail, N protein, and PABP both in vitro and in infected cells. Of the interactions, binding of the poly(A) tail to N protein decreases the interaction efficiency between the poly(A) tail and eIF4E, leading to translation inhibition. The poly(A)-dependent translation inhibition by N protein has not been previously demonstrated and thus extends our understanding of coronavirus gene expression.IMPORTANCE Gene expression in coronavirus is a complicated and dynamic process. In this study, we demonstrated that coronavirus N protein is able to bind to the poly(A) tail with high affinity, establishing N protein as a PABP. We also show how the interplay between coronavirus 3' poly(A) tail, PABP, and N protein regulates gene expression of the coronavirus and host cell. Of the interactions, poly(A) tail binding by the N protein negatively regulates translation, and to our knowledge, this inhibition of translation by binding of the N protein to poly(A) tail has not been previously studied. Accordingly, the study provides fundamental molecular details regarding coronavirus infection and expands our knowledge of coronavirus gene expression.

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.

Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea viruses associated with outbreaks of severe diarrhea in piglets in Jiangxi, China 2013

Porcine epidemic diarrhea (PED), caused by porcine epidemic diarrhea virus (PEDV), is a highly contagious, acute enteric viral disease of swine characterized by vomiting, watery diarrhea, dehydration and death. To identify and characterize the field PEDVs associated with the outbreaks of severe diarrhea in piglets in Jiangxi, 2013, the complete genome sequences of two representative strains of PEDV, designated CH/JX-1/2013 and CH/JX-2/2013, were determined and analyzed. The genome sequences of both emergent Jiangxi PEDV strains, CH/JX-1/2013 and CH/JX-2/2013, were 28,038 nucleotides in length excluding 3’ poly (A) tail. Compared to the PEDV CV777 strain, CH/JX-1/2013 and CH/JX-2/2013 had some unique genetic characteristics in the proximal region of the 5´-UTRs. Phylogenetic analysis of the complete genomes and the structural proteins revealed that CH/JX-1/2013 and CH/JX-2/2013 had a close relationship with post-2010 Chinese PEDV strains and US strains identified in 2013. The nucleotide identity between the two Jiangxi strains (CH/JX-1/2013 and CH/JX-2/2013) and 30 strains of PEDV identified ante-2010 and post-2010 ranged from 96.3–97.0% and 97.3–99.7%, respectively. Multiple nucleotide and deduced amino acid mutations were observed in the ORF1a/b, S, ORF3, E, M and N genes among the current field PEDV strains when compared to the CV777 strain. Some of the mutations altered the amino acid charge and hydrophilicity, and notably, there was an amino acid substitution in the middle of one neutralizing epitope (L1371I) of the S gene of both CH/JX-1/2013 and CH/JX-2/2013. Taken together, the accumulated genetic variations of the current field PEDV strains might have led to antigenic changes of the viruses, which might confer the less effectiveness or failure of the CV777-based vaccines currently being widely used in Jiangxi, China.

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.

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.

Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States

Coronaviruses are known to infect humans and other animals and cause respiratory and gastrointestinal diseases. Here we report the emergence of porcine epidemic diarrhea virus (PEDV) in the United States and determination of its origin, evolution, and genotypes based on temporal and geographical evidence. Histological lesions in small intestine sections of affected pigs and the complete genomic sequences of three emergent strains of PEDV isolated from outbreaks in Minnesota and Iowa were characterized. Genetic and phylogenetic analyses of the three U.S. strains revealed a close relationship with Chinese PEDV strains and their likely Chinese origin. The U.S. PEDV strains underwent evolutionary divergence, which can be classified into two sublineages. The three emergent U.S. strains are most closely related to a strain isolated in 2012 from Anhui Province in China, which might be the result of multiple recombination events between different genetic lineages or sublineages of PEDV. Molecular clock analysis of the divergent time based on the complete genomic sequences is consistent with the actual time difference, approximately 2 to 3 years, of the PED outbreaks between China (December 2010) and the United States (May 2013). The finding that the emergent U.S. PEDV strains share unique genetic features at the 5'-untranslated region with a bat coronavirus provided further support of the evolutionary origin of PEDV from bats and potential cross-species transmission. The data from this study have important implications for understanding the ongoing PEDV outbreaks in the United States and will guide future efforts to develop effective preventive and control measures against PEDV.

IMPORTANCE

The sudden emergence of porcine epidemic diarrhea virus (PEDV), a coronavirus, for the first time in the United States causes significant economic and public health concerns. Since its recognition in May 2013, PEDV has rapidly spread across the United States, resulting in high mortality in piglets in more than 17 States now. The ongoing outbreaks of Middle East respiratory syndrome coronavirus in humans from countries in or near the Arabian Peninsula and the historical deadly nature of the 2002 outbreaks of severe acute respiratory syndrome coronavirus create further anxiety over the emergence of PEDV in the United States due to the lack of scientific information about the origin and evolution of this emerging coronavirus. Here we report the detailed genetic characterization, origin, and evolution of emergent PEDV strains in the United States. The results provide much needed information to devise effective preventive and control strategies against PEDV in the United States.

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.

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.

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.

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.

Cis-acting RNA elements in human and animal plus-strand RNA viruses

The RNA genomes of plus-strand RNA viruses have the ability to form secondary and higher-order structures that contribute to their stability and to their participation in inter- and intramolecular interactions. Those structures that are functionally important are called cis-acting RNA elements because their functions cannot be complemented in trans. They can be involved not only in RNA/RNA interactions but also in binding of viral and cellular proteins during the complex processes of translation, RNA replication and encapsidation. Most viral cis-acting RNA elements are located in the highly structured 5'- and 3'-nontranslated regions of the genomes but sometimes they also extend into the adjacent coding sequences. In addition, some cis-acting RNA elements are embedded within the coding sequences far away from the genomic ends. Although the functional importance of many of these structures has been confirmed by genetic and biochemical analyses, their precise roles are not yet fully understood. In this review we have summarized what is known about cis-acting RNA elements in nine families of human and animal plus-strand RNA viruses with an emphasis on the most thoroughly characterized virus families, the Picornaviridae and Flaviviridae.

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.

Structural lability in stem-loop 1 drives a 5' UTR-3' UTR interaction in coronavirus replication

The leader RNA of the 5′ untranslated region (UTR) of coronaviral genomes contains two stem–loop structures denoted SL1 and SL2. Herein, we show that SL1 is functionally and structurally bipartite. While the upper region of SL1 is required to be paired, we observe strong genetic selection against viruses that contain a deletion of A35, an extrahelical nucleotide that destabilizes SL1, in favor of genomes that contain a diverse panel of destabilizing second-site mutations, due to introduction of a noncanonical base pair near A35. Viruses containing destabilizing SL1-ΔA35 mutations also contain one of two specific mutations in the 3′ UTR. Thermal denaturation and imino proton solvent exchange experiments reveal that the lower half of SL1 is unstable and that second-site SL1-ΔA35 substitutions are characterized by one or more features of the wild-type SL1. We propose a “dynamic SL1” model, in which the base of SL1 has an optimized lability required to mediate a physical interaction between the 5′ UTR and the 3′ UTR that stimulates subgenomic RNA synthesis. Although not conserved at the nucleotide sequence level, these general structural characteristics of SL1 appear to be conserved in other coronaviral genomes.

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.

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.

Bovine coronavirus 5'-proximal genomic acceptor hotspot for discontinuous transcription is 65 nucleotides wide

Coronaviruses are positive-strand, RNA-dependent RNA polymerase-utilizing viruses that require a polymerase template switch, characterized as discontinuous transcription, to place a 5'-terminal genomic leader onto subgenomic mRNAs (sgmRNAs). The usually precise switch is thought to occur during the synthesis of negative-strand templates for sgmRNA production and to be directed by heptameric core donor sequences within the genome that match an acceptor core (UCUAAAC in the case of bovine coronavirus) near the 3' end of the 5'-terminal genomic leader. Here it is shown that a 22-nucleotide (nt) donor sequence engineered into a packageable bovine coronavirus defective interfering (DI) RNA and made to match a sequence within the 65-nt virus genomic leader caused a template switch yielding an sgmRNA with only a 33-nt minileader. By changing the donor sequence, acceptor sites between genomic nt 33 and 97 (identical between the DI RNA and the viral genome) could be used to generate sgmRNAs detectable by Northern analysis (approximately 2 to 32 molecules per cell) by 24 h postinfection. Whether the switch was intramolecular only was not determined since a potentially distinguishing acceptor region in the DI RNA rapidly conformed to that in the helper virus genome through a previously described template switch known as leader switching. These results show that crossover acceptor sites for discontinuous transcription (i) need not include the UCUAAAC core and (ii) rest within a surprisingly wide 5'-proximal "hotspot." Overlap of this hotspot with that for leader switching and with elements required for RNA replication suggests that it is part of a larger 5'-proximal multifunctional structure.