The UCUAAAC promoter motif is not required for high-frequency leader recombination in bovine coronavirus defective interfering RNA

Classification, replication, and transcription of Nidovirales

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

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

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

A Conserved Stem-Loop Structure within ORF5 Is a Frequent Recombination Hotspot for Porcine Reproductive and Respiratory Syndrome Virus 1 (PRRSV-1) with a Particular Modified Live Virus (MLV) Strain

The emergence of recombinant PRRSV strains has been observed for more than a decade. These recombinant viruses are characterized by a genome that contains genetic material from at least two different parental strains. Due to the advanced sequencing techniques and a growing number of data bank entries, the role of PRRSV recombinants has become increasingly important since they are sometimes associated with clinical outbreaks. Chimeric viruses observed more recently are products of PRRSV wild-type and vaccine strains. Here, we report on three PRRSV-1 isolates from geographically distant farms with differing clinical manifestations. A sequencing and recombination analysis revealed that these strains are crossovers between different wild-type strains and the same modified live virus vaccine strain. Interestingly, the recombination breakpoint of all analyzed isolates appears at the beginning of open reading frame 5 (ORF5). RNA structure predictions indicate a conserved stem loop in close proximity to the recombination hotspot, which is a plausible cause of a polymerase template switch during RNA replication. Further research into the mechanisms of the stem loop is needed to help understand the PRRSV recombination process and the role of MLVs as parental strains.

Recombination in Coronaviruses, with a Focus on SARS-CoV-2

Recombination is a common evolutionary tool for RNA viruses, and coronaviruses are no exception. We review here the evidence for recombination in SARS-CoV-2 and reconcile nomenclature for recombinants, discuss their origin and fitness, and speculate how recombinants could make a difference in the future of the COVID-19 pandemics.

The strand-biased transcription of SARS-CoV-2 and unbalanced inhibition by remdesivir

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a positive single-stranded RNA virus, causes the coronavirus disease 19 pandemic. During the viral replication and transcription, the RNA-dependent RNA polymerase "jumps" along the genome template, resulting in discontinuous negative-stranded transcripts. Although the sense-mRNA architectures of SARS-CoV-2 were reported, its negative strand was unexplored. Here, we deeply sequenced both strands of RNA and found SARS-CoV-2 transcription is strongly biased to form the sense strand with variable transcription efficiency for different genes. During negative strand synthesis, numerous non-canonical fusion transcripts are also formed, driven by 3-15 nt sequence homology scattered along the genome but more prone to be inhibited by SARS-CoV-2 RNA polymerase inhibitor remdesivir. The drug also represses more of the negative than the positive strand synthesis as supported by a mathematic simulation model and experimental quantifications. Overall, this study opens new sights into SARS-CoV-2 biogenesis and may facilitate the antiviral vaccine development and drug design.

Recent Progress in Torovirus Molecular Biology

Torovirus (ToV) has recently been classified into the new family Tobaniviridae, although it belonged to the Coronavirus (CoV) family historically. ToVs are associated with enteric diseases in animals and humans. In contrast to CoVs, which are recognised as pathogens of veterinary and medical importance, little attention has been paid to ToVs because their infections are usually asymptomatic or not severe; for a long time, only one equine ToV could be propagated in cultured cells. However, bovine ToVs, which predominantly cause diarrhoea in calves, have been detected worldwide, leading to economic losses. Porcine ToVs have also spread globally; although they have not caused serious economic losses, coinfections with other pathogens can exacerbate their symptoms. In addition, frequent inter- or intra-recombination among ToVs can increase pathogenesis or unpredicted host adaptation. These findings have highlighted the importance of ToVs as pathogens and the need for basic ToV research. Here, we review recent progress in the study of ToV molecular biology including reverse genetics, focusing on the similarities and differences between ToVs and CoVs.

Direct RNA nanopore sequencing of full-length coronavirus genomes provides novel insights into structural variants and enables modification analysis

Sequence analyses of RNA virus genomes remain challenging owing to the exceptional genetic plasticity of these viruses. Because of high mutation and recombination rates, genome replication by viral RNA-dependent RNA polymerases leads to populations of closely related viruses, so-called "quasispecies." Standard (short-read) sequencing technologies are ill-suited to reconstruct large numbers of full-length haplotypes of (1) RNA virus genomes and (2) subgenome-length (sg) RNAs composed of noncontiguous genome regions. Here, we used a full-length, direct RNA sequencing (DRS) approach based on nanopores to characterize viral RNAs produced in cells infected with a human coronavirus. By using DRS, we were able to map the longest (∼26-kb) contiguous read to the viral reference genome. By combining Illumina and Oxford Nanopore sequencing, we reconstructed a highly accurate consensus sequence of the human coronavirus (HCoV)-229E genome (27.3 kb). Furthermore, by using long reads that did not require an assembly step, we were able to identify, in infected cells, diverse and novel HCoV-229E sg RNAs that remain to be characterized. Also, the DRS approach, which circumvents reverse transcription and amplification of RNA, allowed us to detect methylation sites in viral RNAs. Our work paves the way for haplotype-based analyses of viral quasispecies by showing the feasibility of intra-sample haplotype separation. Even though several technical challenges remain to be addressed to exploit the potential of the nanopore technology fully, our work illustrates that DRS may significantly advance genomic studies of complex virus populations, including predictions on long-range interactions in individual full-length viral RNA haplotypes.

Direct RNA nanopore sequencing of full-length coronavirus genomes provides novel insights into structural variants and enables modification analysis

Sequence analyses of RNA virus genomes remain challenging owing to the exceptional genetic plasticity of these viruses. Because of high mutation and recombination rates, genome replication by viral RNA-dependent RNA polymerases leads to populations of closely related viruses, so-called “quasispecies.” Standard (short-read) sequencing technologies are ill-suited to reconstruct large numbers of full-length haplotypes of (1) RNA virus genomes and (2) subgenome-length (sg) RNAs composed of noncontiguous genome regions. Here, we used a full-length, direct RNA sequencing (DRS) approach based on nanopores to characterize viral RNAs produced in cells infected with a human coronavirus. By using DRS, we were able to map the longest (∼26-kb) contiguous read to the viral reference genome. By combining Illumina and Oxford Nanopore sequencing, we reconstructed a highly accurate consensus sequence of the human coronavirus (HCoV)-229E genome (27.3 kb). Furthermore, by using long reads that did not require an assembly step, we were able to identify, in infected cells, diverse and novel HCoV-229E sg RNAs that remain to be characterized. Also, the DRS approach, which circumvents reverse transcription and amplification of RNA, allowed us to detect methylation sites in viral RNAs. Our work paves the way for haplotype-based analyses of viral quasispecies by showing the feasibility of intra-sample haplotype separation. Even though several technical challenges remain to be addressed to exploit the potential of the nanopore technology fully, our work illustrates that DRS may significantly advance genomic studies of complex virus populations, including predictions on long-range interactions in individual full-length viral RNA haplotypes.

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.

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.

Two deletion variants of Middle East respiratory syndrome coronavirus found in a patient with characteristic symptoms

Significant sequence variation of Middle East respiratory syndrome coronavirus (MERS CoV) has never been detected since it was first reported in 2012. A MERS patient came from Korea to China in late May 2015. The patient was 44 years old and had symptoms including high fever, dry cough with a little phlegm, and shortness of breath, which are roughly consistent with those associated with MERS, and had had close contact with individuals with confirmed cases of MERS.After one month of therapy with antiviral, anti-infection, and immune-enhancing agents, the patient recovered in the hospital and was discharged. A nasopharyngeal swab sample was collected for direct sequencing, which revealed two deletion variants of MERS CoV. Deletions of 414 and 419 nt occurred between ORF5 and the E protein, resulting in a partial protein fusion or truncation of ORF5 and the E protein. Functional analysis by bioinformatics and comparison to previous studies implied that the two variants might be defective in their ability to package MERS CoV. However, the mechanism of how these deletions occurred and what effects they have need to be further investigated.

Continuous and Discontinuous RNA Synthesis in Coronaviruses

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

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.

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.

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.

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.

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.

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.

Widespread distribution and a new recombinant species of Brazilian virus associated with cotton blue disease

BACKGROUND

Cotton blue disease (CBD), an important global cotton crop pathology responsible for major economic losses, is prevalent in the major cotton-producing states of Brazil. Typical CBD symptoms include stunting due to internodal shortening, leaf rolling, intense green foliage, and yellowing veins. Atypical CBD symptoms, including reddish and withered leaves, were also observed in Brazilian cotton fields in 2007. Recently, a Polerovirus named Cotton leafroll dwarf virus (CLRDV) was shown to be associated with CBD.

RESULTS

To understand the distribution and genetic diversity of CLRDV in Brazil, we analyzed 23 CBD-symptomatic plants from susceptible cotton varieties originating from five of the six most important cotton-growing states, from 2004-2007. Here, we report on CLRDV diversity in plants with typical or atypical CBD symptoms by comparing viral coat protein, RNA polymerase (RdRp), and intergenic region genomic sequences.

CONCLUSION

The virus had a widespread distribution with a low genetic diversity; however, three divergent isolates were associated with atypical CBD symptoms. These divergent isolates had a CLRDV-related coat protein but a distinct RdRp sequence, and probably arose from recombination events. Based on the taxonomic rules for the family Luteoviridae, we propose that these three isolates represent isolates of a new species in the genus Polerovirus.

Bovine-like coronaviruses isolated from four species of captive wild ruminants are homologous to bovine coronaviruses, based on complete genomic sequences

We sequenced and analyzed the full-length genomes of four coronaviruses (CoVs), each from a distinct wild-ruminant species in Ohio: sambar deer (Cervus unicolor), a waterbuck (Kobus ellipsiprymnus), a sable antelope (Hippotragus niger), and a white-tailed deer (Odocoileus virginianus). The fecal samples from the sambar deer, the waterbuck, and the white-tailed deer were collected during winter dysentery outbreaks and sporadic diarrhea cases in 1993 and 1994 (H. Tsunemitsu, Z. R. el-Kanawati, D. R. Smith, H. H. Reed, and L. J. Saif, J. Clin. Microbiol. 33:3264-3269, 1995). A fecal sample from a sable antelope was collected in 2003 from an Ohio wild-animal habitat during the same outbreak when a bovine-like CoV from a giraffe (GiCoV) was isolated (M. Hasoksuz, K. Alekseev, A. Vlasova, X. Zhang, D. Spiro, R. Halpin, S. Wang, E. Ghedin, and L. J. Saif, J. Virol. 81:4981-4990, 2007). For two of the CoVs (sambar deer and waterbuck), complete genomes from both the cell culture-adapted and gnotobiotic-calf-passaged strains were also sequenced and analyzed. Phylogenetically, wild-ruminant CoVs belong to group 2a CoVs, with the closest relatedness to recent bovine CoV (BCoV) strains. High nucleotide identities (99.4 to 99.6%) among the wild-ruminant strains and recent BCoV strains (BCoV-LUN and BCoV-ENT, isolated in 1998) further confirm the close relatedness. Comparative genetic analysis of CoVs of captive wild ruminants with BCoV strains suggests that no specific genomic markers are present that allow discrimination between the bovine strains and bovine-like CoVs from captive wild ruminants; furthermore, no specific genetic markers were identified that defined cell cultured or calf-passaged strains or the host origin of strains. The results of this study confirm prior reports of biologic and antigenic similarities between bovine and wild-ruminant CoVs and suggest that cattle may be reservoirs for CoVs that infect captive wild ruminants or vice versa and that these CoVs may represent host range variants of an ancestral CoV.

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.

5'-proximal hot spot for an inducible positive-to-negative-strand template switch by coronavirus RNA-dependent RNA polymerase

Coronaviruses have a positive-strand RNA genome and replicate through the use of a 3' nested set of subgenomic mRNAs each possessing a leader (65 to 90 nucleotides [nt] in length, depending on the viral species) identical to and derived from the genomic leader. One widely supported model for leader acquisition states that a template switch takes place during the generation of negative-strand antileader-containing templates used subsequently for subgenomic mRNA synthesis. In this process, the switch is largely driven by canonical heptameric donor sequences at intergenic sites on the genome that match an acceptor sequence at the 3' end of the genomic leader. With experimentally placed 22-nt-long donor sequences within a bovine coronavirus defective interfering (DI) RNA we have shown that matching sites occurring anywhere within a 65-nt-wide 5'-proximal genomic acceptor hot spot (nt 33 through 97) can be used for production of templates for subgenomic mRNA synthesis from the DI RNA. Here we report that with the same experimental approach, template switches can be induced in trans from an internal site in the DI RNA to the negative-strand antigenome of the helper virus. For these, a 3'-proximal 89-nt acceptor hot spot on the viral antigenome (nt 35 through 123), largely complementary to that described above, was found. Molecules resulting from these switches were not templates for subgenomic mRNA synthesis but, rather, ambisense chimeras potentially exceeding the viral genome in length. The results suggest the existence of a coronavirus 5'-proximal partially double-stranded template switch-facilitating structure of discrete width that contains both the viral genome and antigenome.

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.

Nidovirus transcription: how to make sense...?

Many positive-stranded RNA viruses use subgenomic mRNAs to express part of their genetic information. To produce structural and accessory proteins, members of the order Nidovirales (corona-, toro-, arteri- and roniviruses) generate a 3' co-terminal nested set of at least three and often seven to nine mRNAs. Coronavirus and arterivirus subgenomic transcripts are not only 3' co-terminal but also contain a common 5' leader sequence, which is derived from the genomic 5' end. Their synthesis involves a process of discontinuous RNA synthesis that resembles similarity-assisted RNA recombination. Most models proposed over the past 25 years assume co-transcriptional fusion of subgenomic RNA leader and body sequences, but there has been controversy over the question of whether this occurs during plus- or minus-strand synthesis. In the latter model, which has now gained considerable support, subgenomic mRNA synthesis takes place from a complementary set of subgenome-size minus-strand RNAs, produced by discontinuous minus-strand synthesis. Sense-antisense base-pairing interactions between short conserved sequences play a key regulatory role in this process. In view of the presumed common ancestry of nidoviruses, the recent finding that ronivirus and torovirus mRNAs do not contain a common 5' leader sequence is surprising. Apparently, major mechanistic differences must exist between nidoviruses, which raises questions about the functions of the common leader sequence and nidovirus transcriptase proteins and the evolution of nidovirus transcription. In this review, nidovirus transcription mechanisms are compared, the experimental systems used are critically assessed and, in particular, the impact of recently developed reverse genetic systems is discussed.

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.

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

The 210-nucleotide (nt) 5' untranslated region (UTR) in the positive-strand bovine coronavirus (BCoV) genome is predicted to contain four higher-order structures identified as stem-loops I to IV, which may function as cis-acting elements in genomic RNA replication. Here, we describe evidence that stem-loop IV, a bulged stem-loop mapping at nt 186 through 215, (i) is phylogenetically conserved among group 2 coronaviruses and may have a homolog in groups 1 and 3, (ii) exists as a higher-order structure on the basis of enzyme probing, (iii) is required as a higher-order element for replication of a BCoV defective interfering (DI) RNA in the positive but not the negative strand, and (iv) as a higher-order structure in wild-type (wt) and mutant molecules that replicate, specifically binds six cellular proteins in the molecular mass range of 25 to 58 kDa as determined by electrophoretic mobility shift and UV cross-linking assays; binding to viral proteins was not detected. Interestingly, the predicted stem-loop IV homolog in the severe acute respiratory syndrome (SARS) coronavirus appears to be group 1-like in that it is in part duplicated with a group 1-like conserved loop sequence and is not group 2-like, as would be expected by the SARS coronavirus group 2-like 3' UTR structure. These results together indicate that stem-loop IV in the BCoV 5' UTR is a cis-acting element for DI RNA replication and that it might function through interactions with cellular proteins. It is postulated that stem-loop IV functions similarly in the virus genome.

Discontinuous subgenomic RNA synthesis in arteriviruses is guided by an RNA hairpin structure located in the genomic leader region

Nidoviruses produce an extensive 3'-coterminal nested set of subgenomic (sg) mRNAs, which are used to express structural proteins and sometimes accessory proteins. In arteriviruses and coronaviruses, these mRNAs contain a common 5' leader sequence, derived from the genomic 5' end. The joining of the leader sequence to different segments derived from the 3'-proximal part of the genome (mRNA bodies) presumably involves a unique mechanism of discontinuous minus-strand RNA synthesis in which base pairing between sense and antisense transcription-regulating sequences (TRSs) plays an essential role. The leader TRS is present in the loop of a hairpin structure that functions in sg mRNA synthesis. In this study, the minimal sequences in the 5'-proximal region of the Equine arteritis virus genome that are required for sg RNA synthesis were delimited through mutagenesis. A full-length cDNA clone was engineered in which this domain was duplicated, allowing us to make mutations and monitor their effects on sg RNA synthesis without seriously affecting genome replication and translation. The leader TRS present in the duplicated sequence was used and yielded novel sg mRNAs with significantly extended leaders. Our combined findings suggest that the leader TRS hairpin (LTH) and its immediate flanking sequences are essential for efficient sg RNA synthesis and form an independent functional entity that could be moved 300 nucleotides downstream of its original position in the genome. We hypothesize that a conformational switch in the LTH region regulates the role of the 5'-proximal region of the arterivirus genome in subgenomic RNA synthesis.

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.

Mass spectroscopic characterization of the coronavirus infectious bronchitis virus nucleoprotein and elucidation of the role of phosphorylation in RNA binding by using surface plasmon resonance

Phosphorylation of the coronavirus nucleoprotein (N protein) has been predicted to play a role in RNA binding. To investigate this hypothesis, we examined the kinetics of RNA binding between nonphosphorylated and phosphorylated infectious bronchitis virus N protein with nonviral and viral RNA by surface plasmon resonance (Biacore). Mass spectroscopic analysis of N protein identified phosphorylation sites that were proximal to RNA binding domains. Kinetic analysis, by surface plasmon resonance, indicated that nonphosphorylated N protein bound with the same affinity to viral RNA as phosphorylated N protein. However, phosphorylated N protein bound to viral RNA with a higher binding affinity than nonviral RNA, suggesting that phosphorylation of N protein determined the recognition of virus RNA. The data also indicated that a known N protein binding site (involved in transcriptional regulation) consisting of a conserved core sequence present near the 5' end of the genome (in the leader sequence) functioned by promoting high association rates of N protein binding. Further analysis of the leader sequence indicated that the core element was not the only binding site for N protein and that other regions functioned to promote high-affinity binding.

Regulation of relative abundance of arterivirus subgenomic mRNAs

The subgenomic (sg) mRNAs of arteriviruses (order Nidovirales) form a 5'- and 3'-coterminal nested set with the viral genome. Their 5' common leader sequence is derived from the genomic 5'-proximal region. Fusion of sg RNA leader and "body" segments involves a discontinuous transcription step. Presumably during minus-strand synthesis, the nascent RNA strand is transferred from one site in the genomic template to another, a process guided by conserved transcription-regulating sequences (TRSs) at these template sites. Subgenomic RNA species are produced in different but constant molar ratios, with the smallest RNAs usually being most abundant. Factors thought to influence sg RNA synthesis are size differences between sg RNA species, differences in sequence context between body TRSs, and the mutual influence (or competition) between strand transfer reactions occurring at different body TRSs. Using an Equine arteritis virus infectious cDNA clone, we investigated how body TRS activity affected sg RNA synthesis from neighboring body TRSs. Flanking sequences were standardized by head-to-tail insertion of several copies of an RNA7 body TRS cassette. A perfect gradient of sg RNA abundance, progressively favoring smaller RNA species, was observed. Disruption of body TRS function by mutagenesis did not have a significant effect on the activity of other TRSs. However, deletion of body TRS-containing regions enhanced synthesis of sg RNAs from upstream TRSs but not of those produced from downstream TRSs. The results of this study provide considerable support for the proposed discontinuous extension of minus-strand RNA synthesis as a crucial step in sg RNA synthesis.

Comparison of four regions in the replicase gene of heterologous infectious bronchitis virus strains

Infectious bronchitis virus (IBV) produces six subgenomic (sg) mRNAs, each containing a 64 nucleotide (nt) leader sequence, derived from the 5' end of the genome by a discontinuous process. Several putative functional domains such as a papain-like proteinase (PL(pro)), main protease (M(pro)), RNA-dependent RNA polymerase (RdRp), and RNA helicase encoded by the replicase gene are important for virus replication. We have sequenced four regions of the replicase genes corresponding to the 5'-terminal sequence, PL(pro), M(pro), and RdRp domains from 20 heterologous IBV strains, and compared them with previously published coronavirus sequences. All the coronavirus 5'-termini and PL(pro) domains were divergent, unlike the M(pro) and the RdRp domains that were highly conserved with 28% and 48% conserved residues, respectively. Among IBV strains, the 5' untranslated region including the leader sequence was highly conserved (>94% identical); whereas, the N-terminal coding region and the PL(pro) domains were highly variable ranging from 84.6% to 100%, and 77.6% to 100% identity, respectively. The IBV M(pro) and RdRp domains were highly conserved with 82.7% and 92.7% conserved residues, respectively. The BJ strain was the most different from other IBVs in all four regions of the replicase. Phylogeny-based clustering based on replicase genes was identical to the antigen-based classification of coronaviruses into three groups. However, the IBV strain classification based on replicase gene domains did not correlate with that of the type-specific antigenic groups. The replicase gene sequences of many IBVs recovered from infected chickens were identical to those of vaccine viruses irrespective of serotype, suggesting that either there has been an exchange of genetic material among vaccine and field isolates or that there is a convergent evolution to a specific replicase genotype. There was no correlation between the genotype of any region of the replicase gene and pathotype, suggesting that the replicase is not the sole determinant of IBV pathogenicity.

Secondary structure and function of the 5'-proximal region of the equine arteritis virus RNA genome

Nidoviruses produce an extensive 3'-coterminal nested set of subgenomic mRNAs, which are used to express their structural proteins. In addition, arterivirus and coronavirus mRNAs contain a common 5' leader sequence, derived from the genomic 5' end. The joining of this leader sequence to different segments (mRNA bodies) from the genomic 3'-proximal region presumably involves a unique mechanism of discontinuous minus-strand RNA synthesis. Key elements in this process are the so-called transcription-regulating sequences (TRSs), which determine a base-pairing interaction between sense and antisense viral RNA that is essential for leader-to-body joining. To identify RNA structures in the 5'-proximal region of the equine arteritis virus genome that may be involved in subgenomic mRNA synthesis, a detailed secondary RNA structure model was established using bioinformatics, phylogenetic analysis, and RNA structure probing. According to this structure model, the leader TRS is located in the loop of a prominent hairpin (leader TRS hairpin; LTH). The importance of the LTH was supported by the results of a mutagenesis study using an EAV molecular clone. Besides evidence for a direct role of the LTH in subgenomic RNA synthesis, indications for a role of the LTH region in genome replication and/or translation were obtained. Similar LTH structures could be predicted for the 5'-proximal region of all arterivirus genomes and, interestingly, also for most coronaviruses. Thus, we postulate that the LTH is a key structural element in the discontinuous subgenomic RNA synthesis and is likely critical for leader TRS function.

Common RNA replication signals exist among group 2 coronaviruses: evidence for in vivo recombination between animal and human coronavirus molecules

5′ and 3′ UTR sequences on the coronavirus genome are known to carry cis-acting elements for DI RNA replication and presumably also virus genome replication. 5′ UTR-adjacent coding sequences are also thought to harbor cis-acting elements. Here we have determined the 5′ UTR and adjacent 289-nt sequences, and 3′ UTR sequences, for six group 2 coronaviruses and have compared them to each other and to three previously reported group 2 members. Extensive regions of highly similar UTR sequences were found but small regions of divergence were also found indicating group 2 coronaviruses could be subdivided into those that are bovine coronavirus (BCoV)-like (BCoV, human respiratory coronavirus-OC43, human enteric coronavirus, porcine hemagglutinating encephalomyelitis virus, and equine coronavirus) and those that are murine hepatitis virus (MHV)-like (A59, 2, and JHM strains of MHV, puffinosis virus, and rat sialodacryoadenitis virus). The 3′ UTRs of BCoV and MHV have been previously shown to be interchangeable. Here, a reporter-containing BCoV DI RNA was shown to be replicated by all five BCoV-like helper viruses and by MHV-H2 (a human cell-adapted MHV strain), a representative of the MHV-like subgroup, demonstrating group 2 common 5′ and 3′ replication signaling elements. BCoV DI RNA, furthermore, acquired the leader of HCoV-OC43 by leader switching, demonstrating for the first time in vivo recombination between animal and human coronavirus molecules. These results indicate that common replication signaling elements exist among group 2 coronaviruses despite a two-cluster pattern within the group and imply there could exist a high potential for recombination among group members.

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

Higher-order structures in the 5' untranslated region (UTR) of plus-strand RNA viruses are known in many cases to function as cis-acting elements in RNA translation, replication, or transcription. Here we describe evidence supporting the structure and a cis-acting function in defective interfering (DI) RNA replication of stem-loop III, the third of four predicted higher-order structures mapping within the 210-nucleotide (nt) 5' UTR of the 32-kb bovine coronavirus (BCoV) genome. Stem-loop III maps at nt 97 through 116, has a calculated free energy of -9.1 kcal/mol in the positive strand and -3.0 kcal/mol in the negative strand, and has associated with it beginning at nt 100 an open reading frame (ORF) potentially encoding an 8-amino-acid peptide. Stem-loop III is presumed to function in the positive strand, but its strand of action has not been established. Stem-loop III (i) shows phylogenetic conservation among group 2 coronaviruses and appears to have a homolog in coronavirus groups 1 and 3, (ii) has in all coronaviruses for which sequence is known a closely associated short, AUG-initiated intra-5' UTR ORF, (iii) is supported by enzyme structure-probing evidence in BCoV RNA, (iv) must maintain stem integrity for DI RNA replication in BCoV DI RNA, and (v) shows a positive correlation between maintenance of the short ORF and maximal DI RNA accumulation in BCoV DI RNA. These results indicate that stem-loop III in the BCoV 5' UTR is a cis-acting element for DI RNA replication and that its associated intra-5' UTR ORF may function to enhance replication. It is postulated that these two elements function similarly in the virus genome.

Discontinuous and non-discontinuous subgenomic RNA transcription in a nidovirus

Arteri-, corona-, toro- and roniviruses are evolutionarily related positive-strand RNA viruses, united in the order Nidovirales. The best studied nidoviruses, the corona- and arteriviruses, employ a unique transcription mechanism, which involves discontinuous RNA synthesis, a process resembling similarity-assisted copy-choice RNA recombination. During infection, multiple subgenomic (sg) mRNAs are transcribed from a mirror set of sg negative-strand RNA templates. The sg mRNAs all possess a short 5' common leader sequence, derived from the 5' end of the genomic RNA. The joining of the non-contiguous 'leader' and 'body' sequences presumably occurs during minus-strand synthesis. To study whether toroviruses use a similar transcription mechanism, we characterized the 5' termini of the genome and the four sg mRNAs of Berne virus (BEV). We show that BEV mRNAs 3-5 lack a leader sequence. Surprisingly, however, RNA 2 does contain a leader, identical to the 5'-terminal 18 residues of the genome. Apparently, BEV combines discontinuous and non-discontinuous RNA synthesis to produce its sg mRNAs. Our findings have important implications for the understanding of the mechanism and evolution of nidovirus transcription.

Recombination and Coronavirus Defective Interfering RNAs

Naturally occurring defective interfering RNAs have been found in 4 of 14 coronavirus species. They range in size from 2.2 kb to approximately 25 kb, or 80% of the 30-kb parent virus genome. The large DI RNAs do not in all cases appear to require helper virus for intracellular replication and it has been postulated that they may on their own function as agents of disease. Coronavirus DI RNAs appear to arise by internal deletions (through nonhomologous recombination events) on the virus genome or on DI RNAs of larger size by a polymerase strand-switching (copy-choice) mechanism. In addition to their use in the study of virus RNA replication and virus assembly, coronavirus DI RNAs are being used in a major way to study the mechanism of a high-frequency, site-specific RNA recombination event that leads to leader acquisition during virus replication (i.e., the leader fusion event that occurs during synthesis of subgenomic mRNAs, and the leader-switching event that can occur during DI RNA replication), a distinguishing feature of coronaviruses (and arteriviruses). Coronavirus DI RNAs are also being engineered as vehicles for the generation of targeted recombinants of the parent virus genome.

Sequence requirements for RNA strand transfer during nidovirus discontinuous subgenomic RNA synthesis

Nidovirus subgenomic mRNAs contain a leader sequence derived from the 5' end of the genome fused to different sequences ('bodies') derived from the 3' end. Their generation involves a unique mechanism of discontinuous subgenomic RNA synthesis that resembles copy-choice RNA recombination. During this process, the nascent RNA strand is transferred from one site in the template to another, during either plus or minus strand synthesis, to yield subgenomic RNA molecules. Central to this process are transcription-regulating sequences (TRSs), which are present at both template sites and ensure the fidelity of strand transfer. Here we present results of a comprehensive co-variation mutagenesis study of equine arteritis virus TRSs, demonstrating that discontinuous RNA synthesis depends not only on base pairing between sense leader TRS and antisense body TRS, but also on the primary sequence of the body TRS. While the leader TRS merely plays a targeting role for strand transfer, the body TRS fulfils multiple functions. The sequences of mRNA leader-body junctions of TRS mutants strongly suggested that the discontinuous step occurs during minus strand synthesis.

Enhanced accumulation of coronavirus defective interfering RNA from expressed negative-strand transcripts by coexpressed positive-strand RNA transcripts

Expression of negative-strand murine coronavirus mouse hepatitis virus (MHV) defective interfering (DI) RNA transcripts in MHV-infected cells results in the accumulation of positive-strand DI RNAs (M. Joo et al., 1996, J. Virol. 70, 5769-5776). However, the expressed negative-strand DI RNA transcripts are poor templates for positive-strand DI RNA synthesis. The present study demonstrated that DI RNA accumulation from the expressed negative-strand DI RNA transcripts in MHV-infected cells was enhanced by the coexpression of complementary RNA transcripts that correspond to the 5' region of positive-strand DI RNA. The positive-strand RNA transcripts corresponding to the 5' end-most 0.7-2.0 kb DI RNA had a similar enhancement effect. The coexpressed positive-strand RNA transcripts lacking the leader sequence or those containing only the leader sequence failed to demonstrate this enhancement effect, demonstrating that the presence of the leader sequence in the coexpressed positive-strand RNA transcripts was necessary, but not sufficient, for the enhancement of DI RNA accumulation from the coexpressed negative-strand DI RNA transcripts. Negative-strand DI RNA transcripts that were coexpressed with the partial-length positive-strand RNA transcripts were no more stable than those expressed alone, suggesting that a higher stability of the expressed negative-strand RNA transcripts was an unlikely reason for the higher DI RNA accumulation in cells coexpressing two complementary DI RNA transcripts. Sequence analyses unexpectedly demonstrated that the leader sequence of the majority of accumulated DI RNAs switched to helper virus derived leader sequence, suggesting that enhancement of DI RNA accumulation was mediated by the efficient utilization of helper virus derived leader sequence for DI RNA synthesis. Furthermore, our data suggested that this leader switching, a type of homologous RNA-RNA recombination, occurred during positive-strand DI RNA synthesis and that MHV positive-strand RNA synthesis mechanism may have a preference toward recognizing double-stranded RNA structures over single-stranded negative-strand RNA to produce positive-strand DI RNAs.

Downstream sequences influence the choice between a naturally occurring noncanonical and closely positioned upstream canonical heptameric fusion motif during bovine coronavirus subgenomic mRNA synthesis

Mechanisms leading to subgenomic mRNA (sgmRNA) synthesis in coronaviruses are poorly understood but are known to involve a heptameric signaling motif, originally called the intergenic sequence. The intergenic sequence is the presumed crossover region (fusion site) for RNA-dependent RNA polymerase (RdRp) during discontinuous transcription, a process leading to sgmRNAs that are both 5' and 3' coterminal. In the bovine coronavirus, the major fusion site for synthesis of mRNA 5 (GGUAGAC) does not conform to the canonical motif (UC[U,C]AAAC) at three positions (underlined), yet it lies just 14 nucleotides downstream from such a sequence (UCCAAAC). The infrequently used canonical sequence, by computer prediction, is buried within the stem of a stable hairpin (-17.2 kcal/mol). Here we document the existence of this stem by enzyme probing and examine its influence and that of neighboring sequences on the unusual choice of fusion sites by analyzing transcripts made in vivo from mutated defective interfering RNA constructs. We learned that (i) mutations that were predicted to unfold the stem-loop in various ways did not switch RdRp crossover to the upstream canonical site, (ii) a totally nonconforming downstream motif resulted in no measurable transcription from either site, (iii) the canonical upstream site does not function ectopically to lend competence to the downstream noncanonical site, and (iv) altering flanking sequences downstream of the downstream noncanonical motif in ways that diminish sequence similarity with the virus genome 5' end caused a dramatic switch to the upstream canonical site. These results show that sequence elements downstream of the noncanonical site can dramatically influence the choice of fusion sites for synthesis of mRNA 5 and are interpreted as being most consistent with a mechanism of similarity-assisted RdRp strand switching during minus-strand synthesis.