Mutagenic analysis of the coronavirus intergenic consensus sequence

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.

High-Resolution Analysis of Coronavirus Gene Expression by RNA Sequencing and Ribosome Profiling

Members of the family Coronaviridae have the largest genomes of all RNA viruses, typically in the region of 30 kilobases. Several coronaviruses, such as Severe acute respiratory syndrome-related coronavirus (SARS-CoV) and Middle East respiratory syndrome-related coronavirus (MERS-CoV), are of medical importance, with high mortality rates and, in the case of SARS-CoV, significant pandemic potential. Other coronaviruses, such as Porcine epidemic diarrhea virus and Avian coronavirus, are important livestock pathogens. Ribosome profiling is a technique which exploits the capacity of the translating ribosome to protect around 30 nucleotides of mRNA from ribonuclease digestion. Ribosome-protected mRNA fragments are purified, subjected to deep sequencing and mapped back to the transcriptome to give a global "snap-shot" of translation. Parallel RNA sequencing allows normalization by transcript abundance. Here we apply ribosome profiling to cells infected with Murine coronavirus, mouse hepatitis virus, strain A59 (MHV-A59), a model coronavirus in the same genus as SARS-CoV and MERS-CoV. The data obtained allowed us to study the kinetics of virus transcription and translation with exquisite precision. We studied the timecourse of positive and negative-sense genomic and subgenomic viral RNA production and the relative translation efficiencies of the different virus ORFs. Virus mRNAs were not found to be translated more efficiently than host mRNAs; rather, virus translation dominates host translation at later time points due to high levels of virus transcripts. Triplet phasing of the profiling data allowed precise determination of translated reading frames and revealed several translated short open reading frames upstream of, or embedded within, known virus protein-coding regions. Ribosome pause sites were identified in the virus replicase polyprotein pp1a ORF and investigated experimentally. Contrary to expectations, ribosomes were not found to pause at the ribosomal frameshift site. To our knowledge this is the first application of ribosome profiling to an RNA virus.

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.

Coronaviruses as vectors: stability of foreign gene expression

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

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 reverse genetics and development of vectors for gene expression

Knowledge of coronavirus replication, transcription, and virus-host interaction has been recently improved by engineering of coronavirus infectious cDNAs. With the transmissible gastroenteritis virus (TGEV) genome the efficient (>40 microg per 106 cells) and stable (>20 passages) expression of the foreign genes has been shown. Knowledge of the transcription mechanism in coronaviruses has been significantly increased, making possible the fine regulation of foreign gene expression. A new family of vectors based on single coronavirus genomes, in which essential genes have been deleted, has emerged including replication-competent, propagation-deficient vectors. Vector biosafety is being increased by relocating the RNA packaging signal to the position previously occupied by deleted essential genes, to prevent the rescue of fully competent viruses that might arise from recombination events with wild-type field coronaviruses. The large cloning capacity of coronaviruses (>5 kb) and the possibility of engineering the tissue and species tropism to target expression to different organs and animal species, including humans, has increased the potential of coronaviruses as vectors for vaccine development and, possibly, gene therapy.

The severe acute respiratory syndrome coronavirus Nsp15 protein is an endoribonuclease that prefers manganese as a cofactor

Nonstructural protein 15 (Nsp15) of the severe acute respiratory syndrome coronavirus (SARS-CoV) produced in Escherichia coli has endoribonuclease activity that preferentially cleaved 5' of uridylates of RNAs. Blocking either the 5' or 3' terminus did not affect cleavage. Double- and single-stranded RNAs were both substrates for Nsp15 but with different kinetics for cleavage. Mn(2+) at 2 to 10 mM was needed for optimal endoribonuclease activity, but Mg(2+) and several other divalent metals were capable of supporting only a low level of activity. Concentrations of Mn(2+) needed for endoribonuclease activity induced significant conformation change(s) in the protein, as measured by changes in tryptophan fluorescence. A similar endoribonucleolytic activity was detected for the orthologous protein from another coronavirus, demonstrating that the endoribonuclease activity of Nsp15 may be common to coronaviruses. This work presents an initial biochemical characterization of a novel coronavirus endoribonuclease.

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.

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.

Coronaviruses as vectors: position dependence of foreign gene expression

Coronaviruses are the enveloped, positive-stranded RNA viruses with the largest RNA genomes known. Several features make these viruses attractive as vaccine and therapeutic vectors: (i) deletion of their nonessential genes is strongly attenuating; (ii) the genetic space thus created allows insertion of foreign information; and (iii) their tropism can be modified by manipulation of the viral spike. We studied here their ability to serve as expression vectors by inserting two different foreign genes and evaluating systematically the genomic position dependence of their expression, using a murine coronavirus as a model. Renilla and firefly luciferase expression cassettes, each provided with viral transcription regulatory sequences (TRSs), were inserted at several genomic positions, both independently in different viruses and combined within one viral genome. Recombinant viruses were generated by using a convenient method based on targeted recombination and host cell switching. In all cases high expression levels of the foreign genes were observed without severe effects on viral replication in vitro. The expression of the inserted gene appeared to be dependent on its genomic position, as well as on the identity of the gene. Expression levels increased when the luciferase gene was inserted closer to the 3' end of the genome. The foreign gene insertions generally reduced the expression of upstream viral genes. The results are consistent with coronavirus transcription models in which the transcription from upstream TRSs is attenuated by downstream TRSs. Altogether, our observations clearly demonstrate the potential of coronaviruses as (multivalent) expression vectors.

The stability of the duplex between sense and antisense transcription-regulating sequences is a crucial factor in arterivirus subgenomic mRNA synthesis

Subgenomic mRNAs of nidoviruses (arteriviruses and coronaviruses) are composed of a common leader sequence and a "body" part of variable size, which are derived from the 5'- and 3'-proximal part of the genome, respectively. Leader-to-body joining has been proposed to occur during minus-strand RNA synthesis and to involve transfer of the nascent RNA strand from one site in the template to another. This discontinuous step in subgenomic RNA synthesis is guided by short transcription-regulating sequences (TRSs) that are present at both these template sites (leader TRS and body TRS). Sense-antisense base pairing between the leader TRS in the plus strand and the body TRS complement in the minus strand is crucial for strand transfer. Here we show that extending the leader TRS-body TRS duplex beyond its wild-type length dramatically enhanced the subgenomic mRNA synthesis of the arterivirus Equine arteritis virus (EAV). Generally, the relative amount of a subgenomic mRNA correlated with the calculated stability of the corresponding leader TRS-body TRS duplex. In addition, various leader TRS mutations induced the generation of minor subgenomic RNA species that were not detected upon infection with wild-type EAV. The synthesis of these RNA species involved leader-body junction events at sites that bear only limited resemblance to the canonical TRS. However, with the mutant leader TRS, but not with the wild-type leader TRS, these sequences could form a duplex that was stable enough to direct subgenomic RNA synthesis, again demonstrating that the stability of the leader TRS-body TRS duplex is a crucial factor in arterivirus subgenomic mRNA synthesis.

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.

Transcription regulatory sequences and mRNA expression levels in the coronavirus transmissible gastroenteritis virus

The transcription regulatory sequences (TRSs) of the coronavirus transmissible gastroenteritis virus (TGEV) have been characterized by using a helper virus-dependent expression system based on coronavirus-derived minigenomes to study the synthesis of subgenomic mRNAs. The TRSs are located at the 5' end of TGEV genes and include a highly conserved core sequence (CS), 5'-CUAAAC-3', that is essential for mediating a 100- to 1,000-fold increase in mRNA synthesis when it is located in the appropriate context. The relevant sequences contributing to TRS activity have been studied by extending the CS 5' upstream and 3' downstream. Sequences from virus genes flanking the CS influenced transcription levels from moderate (10- to 20-fold variation) to complete mRNA synthesis silencing, as shown for a canonical CS at nucleotide (nt) 120 from the initiation codon of the S gene that did not lead to the production of the corresponding mRNA. An optimized TRS has been designed comprising 88 nt from the N gene TRS, the CS, and 3 nt 3' to the M gene CS. Further extension of the 5'-flanking nucleotides (i.e., by 176 nt) decreased subgenomic RNA levels. The expression of a reporter gene (beta-glucuronidase) by using the selected TRS led to the production of 2 to 8 microg of protein per 10(6) cells. The presence of an appropriate Kozak context led to a higher level of protein expression. Virus protein levels were shown to be dependent on transcription and translation regulation.

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.

Coronavirus derived expression systems

Both helper dependent expression systems, based on two components, and single genomes constructed by targeted recombination, or by using infectious cDNA clones, have been developed. The sequences that regulate transcription have been characterized mainly using helper dependent expression systems and it will now be possible to validate them using single genomes. The genome of coronaviruses has been engineered by modification of the infectious cDNA leading to an efficient (>20 microg ml(-1)) and stable (>20 passages) expression of the foreign gene. The possibility of engineering the tissue and species tropism to target expression to different organs and animal species, including humans, increases the potential of coronaviruses as vectors. Thus, coronaviruses are promising virus vectors for vaccine development and, possibly, for gene therapy.

Genetic manipulation of arterivirus alternative mRNA leader-body junction sites reveals tight regulation of structural protein expression

To express its structural proteins, the arterivirus Equine arteritis virus (EAV) produces a nested set of six subgenomic (sg) RNA species. These RNA molecules are generated by a mechanism of discontinuous transcription, during which a common leader sequence, representing the 5' end of the genomic RNA, is attached to the bodies of the sg RNAs. The connection between the leader and body parts of an mRNA is formed by a short, conserved sequence element termed the transcription-regulating sequence (TRS), which is present at the 3' end of the leader as well as upstream of each of the structural protein genes. With the exception of RNA3, only one body TRS was previously assumed to be used to join the leader and body of each EAV sg RNA. Here we show that for the synthesis of two other sg RNAs, RNA4 and RNA5, alternative leader-body junction sites that differ substantially in transcriptional activity are used. By site-directed mutagenesis of an EAV infectious cDNA clone, the alternative TRSs used to generate RNA3, -4, and -5 were inactivated, which strongly influenced the corresponding RNA levels and the production of infectious progeny virus. The relative amounts of RNA produced from alternative TRSs differed significantly and corresponded to the relative infectivities of the virus mutants. This strongly suggested that the structural proteins that are expressed from these RNAs are limiting factors during the viral life cycle and that the discontinuous step in sg RNA synthesis is crucial for the regulation of their expression. On the basis of a theoretical analysis of the predicted RNA structure of the 3' end of the EAV genome, we propose that the local secondary RNA structure of the body TRS regions is an important factor in the regulation of the discontinuous step in EAV sg mRNA synthesis.

Expression of reporter genes from the defective RNA CD-61 of the coronavirus infectious bronchitis virus

The defective RNA (D-RNA) CD-61, derived from the Beaudette strain of the avian coronavirus infectious bronchitis virus (IBV), was used as an RNA vector for the expression of two reporter genes, luciferase and chloramphenicol acetyltransferase (CAT). D-RNAs expressing the CAT gene were demonstrated to be capable of producing CAT protein in a helper-dependent expression system to about 1.6 microgram per 10(6) cells. The reporter genes were expressed from two different sites within the CD-61 sequence and expression was not affected by interruption of the CD-61-specific ORF. Expression of the reporter genes was under the control of a transcription-associated sequence (TAS) derived from the Beaudette gene 5, normally used for the transcription of IBV subgenomic mRNA 5. The Beaudette gene 5 TAS is composed of two tandem repeats of the IBV canonical consensus sequence involved in the acquisition of a leader sequence during the discontinuous transcription of IBV subgenomic mRNAs. It is demonstrated that only one canonical sequence is required for expression of mRNA 5 or for the expression of an mRNA from a D-RNA and that either sequence can function as an acceptor site for acquisition of the leader sequence.

High affinity interaction between nucleocapsid protein and leader/intergenic sequence of mouse hepatitis virus RNA

The nucleocapsid (N) protein of mouse hepatitis virus (MHV) is the major virion structural protein. It associates with both viral genomic RNA and subgenomic mRNAs and has structural and non-structural roles in replication including viral RNA-dependent RNA transcription, genome replication, encapsidation and translation. These processes all involve RNA-protein interactions between the N protein and viral RNAs. To better understand the RNA-binding properties of this multifunctional protein, the N protein was expressed in Escherichia coli as a chimeric protein fused to glutathione-S-transferase (GST). Biochemical analyses of RNA-binding properties were performed on full-length and partial N protein segments to define the RNA-binding domain. The full-length N protein and the GST-N protein fusion product had similar binding activities with a dissociation constant (K(d)) of 14 nM when the MHV 5'-leader sequence was used as ligand. The smallest N protein fragment which retained RNA-binding activity was a 55 aa segment containing residues 177-231 which bound viral RNA with a K(d) of 32 nM. A consensus viral sequence recognized by the N protein was inferred from these studies; AAUCYAAAC was identified to be the potential minimum ligand for the N protein. Although the core UCYAA sequence is often tandemly repeated in viral genomes, ligands containing one or more repeats of UCYAA showed no difference in binding to the N protein. Together these data demonstrate a high-affinity, specific interaction between the N protein and a conserved RNA sequence present at the 5'-ends of MHV mRNA.

Subgenomic mRNA regulation by a distal RNA element in a (+)-strand RNA virus

Subgenomic (sg) mRNAs are synthesized by (+)-strand RNA viruses to allow for efficient translation of products encoded 3' in their genomes. This strategy also provides a means for regulating the expression of such products via modulation of sg mRNA accumulation. We have studied the mechanism by which sg mRNAs levels are controlled in tomato bushy stunt virus, a small (+)-strand RNA virus which synthesizes two sg mRNAs during infections. Neither the viral capsid nor movement proteins were found to play any significant role in modulating the accumulation levels of either sg mRNA. Deletion analysis did, however, identify a 12-nt-long RNA sequence located approximately 1,000 nt upstream from the site of initiation of sg mRNA2 synthesis that was required specifically for accumulation of sg mRNA2. Further analysis revealed a potential base-pairing interaction between this sequence and a sequence located just 5' to the site of initiation for sg mRNA2 synthesis. Mutant genomes in which this interaction was either disrupted or maintained were analyzed and the results indicated a positive correlation between the predicted stability of the base-pairing interaction and the efficiency of sg mRNA2 accumulation. The functional significance of the long-distance interaction was further supported by phylogenetic sequence analysis which revealed conservation of base-pairing interactions of similar stability and relative position in the genomes of different tombusviruses. It is proposed that the upstream sequence represents a cis-acting RNA element which facilitates sg mRNA accumulation by promoting efficient synthesis of sg mRNA2 via a long-distance RNA-RNA interaction.

Polypyrimidine tract-binding protein binds to the leader RNA of mouse hepatitis virus and serves as a regulator of viral transcription

A cellular protein, previously described as p55, binds specifically to the plus strand of the mouse hepatitis virus (MHV) leader RNA. We have purified this protein and determined by partial peptide sequencing that it is polypyrimidine tract-binding protein (PTB) (also known as heterogeneous nuclear ribonucleoprotein [hnRNP] I), a nuclear protein which shuttles between the nucleus and cytoplasm. PTB plays a role in the regulation of alternative splicing of pre-mRNAs in normal cells and translation of several viruses. By UV cross-linking and immunoprecipitation studies using cellular extracts and a recombinant PTB, we have established that PTB binds to the MHV plus-strand leader RNA specifically. Deletion analyses of the leader RNA mapped the PTB-binding site to the UCUAA pentanucleotide repeats. Using a defective-interfering RNA reporter system, we have further shown that the PTB-binding site in the leader RNA is critical for MHV RNA synthesis. This and our previous study (H.-P. Li, X. Zhang, R. Duncan, L. Comai, and M. M. C. Lai, Proc. Natl. Acad. Sci. USA 94:9544-9549, 1997) combined thus show that two cellular hnRNPs, PTB and hnRNP A1, bind to the transcription-regulatory sequences of MHV RNA and may participate in its transcription.

Importance of the positive-strand RNA secondary structure of a murine coronavirus defective interfering RNA internal replication signal in positive-strand RNA synthesis

The RNA elements that are required for replication of defective interfering (DI) RNA of the JHM strain of mouse hepatitis virus (MHV) consist of three discontinuous genomic regions: about 0.46 to 0.47 kb from both terminal sequences and an internal 58-nucleotide (nt)-long sequence (58-nt region) present at about 0.9 kb from the 5' end of the DI genome. The internal region is important for positive-strand DI RNA synthesis (Y. N. Kim and S. Makino, J. Virol. 69:4963-4971, 1995). We further characterized the 58-nt region in the present study and obtained the following results. (i) The positive-strand RNA structure in solution was comparable with that predicted by computer modeling. (ii) Positive-strand RNA secondary structure, but not negative-strand RNA structure, was important for the biological function of the region. (iii) The biological function had a sequence-specific requirement. We discuss possible mechanisms by which the internal cis-acting signal drives MHV positive-strand DI RNA synthesis.

Characterizations of coronavirus cis-acting RNA elements and the transcription step affecting its transcription efficiency

Seven to eight species of viral subgenomic mRNAs are produced in coronavirus-infected cells. These mRNAs are produced in different quantities, and their molar ratios remain constant during viral replication. We studied RNA elements that affect coronavirus transcription efficiency by characterizing a series of cloned coronavirus mouse hepatitis virus (MHV) defective interfering (DI) RNAs containing an inserted intergenic sequence, from which subgenomic DI RNA is transcribed in MHV-infected cells. Certain combinations of upstream and downstream flanking sequences of the intergenic sequence suppressed subgenomic DI RNA transcription, yet changing one of the flanking sequences to a different sequence eliminated transcription suppression. The suppressive effect of certain combinations of flanking sequences, but not all combinations, could be counteracted by altering the intergenic sequence. Thus, the combination of intergenic sequence and flanking sequence affected transcription efficiency. We also characterized another set of DI RNAs designed to clarify which transcription step determines the relative molar ratios of coronavirus mRNAs. Our study indicated that if subgenomic mRNAs were exclusively synthesized from negative-strand genomic RNA, then the relative molar ratios of coronavirus mRNAs were most likely determined after synthesis of the genomic-sized template RNA. If negative-strand subgenomic RNAs were templates for subgenomic mRNAs, then the relative molar ratios of coronavirus mRNAs probably were determined after synthesis of the genomic-sized template RNA used for subgenomic-sized RNA transcription but prior to the completion of the synthesis of subgenomic-sized RNAs containing the leader sequence. The relative molar ratios of coronavirus mRNAs, therefore, seem to have been established prior to a putative replicon-type amplification of subgenomic mRNAs.

The molecular biology of coronaviruses

This chapter discusses the manipulation of clones of coronavirus and of complementary DNAs (cDNAs) of defective-interfering (DI) RNAs to study coronavirus RNA replication, transcription, recombination, processing and transport of proteins, virion assembly, identification of cell receptors for coronaviruses, and processing of the polymerase. The nature of the coronavirus genome is nonsegmented, single-stranded, and positive-sense RNA. Its size ranges from 27 to 32 kb, which is significantly larger when compared with other RNA viruses. The gene encoding the large surface glycoprotein is up to 4.4 kb, encoding an imposing trimeric, highly glycosylated protein. This soars some 20 nm above the virion envelope, giving the virus the appearance-with a little imagination-of a crown or coronet. Coronavirus research has contributed to the understanding of many aspects of molecular biology in general, such as the mechanism of RNA synthesis, translational control, and protein transport and processing. It remains a treasure capable of generating unexpected insights.

The molecular biology of coronaviruses

This chapter discusses the manipulation of clones of coronavirus and of complementary DNAs (cDNAs) of defective-interfering (DI) RNAs to study coronavirus RNA replication, transcription, recombination, processing and transport of proteins, virion assembly, identification of cell receptors for coronaviruses, and processing of the polymerase. The nature of the coronavirus genome is nonsegmented, single-stranded, and positive-sense RNA. Its size ranges from 27 to 32 kb, which is significantly larger when compared with other RNA viruses. The gene encoding the large surface glycoprotein is up to 4.4 kb, encoding an imposing trimeric, highly glycosylated protein. This soars some 20 nm above the virion envelope, giving the virus the appearance-with a little imagination-of a crown or coronet. Coronavirus research has contributed to the understanding of many aspects of molecular biology in general, such as the mechanism of RNA synthesis, translational control, and protein transport and processing. It remains a treasure capable of generating unexpected insights.

Replication of murine coronavirus defective interfering RNA from negative-strand transcripts

The positive-strand defective interfering (DI) RNA of the murine coronavirus mouse hepatitis virus (MHV), when introduced into MHV-infected cells, results in DI RNA replication and accumulation. We studied whether the introduction of negative-strand transcripts of MHV DI RNA would also result in replication. At a location downstream of the T7 promoter and upstream of the human hepatitis delta virus ribozyme domain, we inserted a complete cDNA clone of MHV DI RNA in reverse orientation; in vitro-synthesized RNA from this plasmid yielded a negative-strand RNA copy of the MHV DI RNA. When the negative-strand transcripts of the DI RNA were expressed in MHV-infected cells by a vaccinia virus T7 expression system, positive-strand DI RNAs accumulated in the plasmid-transfected cells. DI RNA replication depended on the expression of T7 polymerase and on the presence of the T7 promoter. Transfection of in vitro-synthesized negative-strand transcripts into MHV-infected cells and serial passage of virus samples from RNA-transfected cells also resulted in accumulation of the DI RNA. Positive-strand DI RNA transcripts were undetectable in sample preparations of the in vitro-synthesized negative-strand DI RNA transcripts, and DI RNA did not accumulate after cotransfection of a small amount of positive-strand DI RNA and truncated-replication-disabled negative-strand transcripts; clearly, the DI RNA replicated from the transfected negative-strand transcripts and not from minute amounts of positive-strand DI RNAs that might be envisioned as artifacts of T7 transcription. Sequence analysis of positive-strand DI RNAs in the cells transfected with negative-strand transcripts showed that DI RNAs maintained the DI-specific unique sequences introduced within the leader sequence. These data indicated that positive-strand DI RNA synthesis occurred from introduced negative-strand transcripts in the MHV-infected cells; this demonstration, using MHV, of DI RNA replication from transfected negative-strand DI RNA transcripts is the first such demonstration among all positive-stranded RNA viruses.

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

The 65-nucleotide leader on the cloned bovine coronavirus defective interfering (DI) RNA, when marked by mutations, has been shown to rapidly convert to the wild-type leader of the helper virus following DI RNA transfection into helper virus-infected cells. A model of leader-primed transcription in which free leader supplied in trans by the helper virus interacts by way of its flanking 5'UCUAAAC3' sequence element with the 3'-proximal 3'AGAUUUG5' promoter on the DI RNA minus strand to prime RNA replication has been used to explain this phenomenon. To test this model, the UCUAAAC element which occurs only once in the BCV 5' untranslated region was either deleted or completely substituted in input DI RNA template, and evidence of leader conversion was sought. In both cases, leader conversion occurred rapidly, indicating that this element is not required on input RNA for the conversion event. Substitution mutations mapped the crossover region to a 24-nucleotide segment that begins within the UCUAAAC sequence and extends downstream. Although structure probing of the bovine coronavirus 5' untranslated region indicated that the UCUAAAC element is in the loop of a prominent stem and thus theoretically available for base pair-directed priming, no evidence of an unattached leader early in infection that might have served as a primer for transcription was found by RNase protection studies. These results together suggest that leader conversion on the DI RNA 5' terminus is not guided by the UCUAAAC element and might arise instead from a high-frequency, region-specific, homologous recombination event perhaps during minus-strand synthesis rather than by leader priming during plus-strand synthesis.

A 5'-proximal RNA sequence of murine coronavirus as a potential initiation site for genomic-length mRNA transcription

Coronavirus transcription is a discontinuous process, involving interactions between a trans-acting leader and the intergenic transcription initiation sequences. A 9-nucleotide (nt) sequence (UUUAUAAAC), which is located immediately downstream of the leader at the 5' terminus of the mouse hepatitis virus (MHV) genomic RNA, contains a sequence resembling the consensus intergenic sequence (UCUAAAC). It has been shown previously that the presence of the 9-nt sequence facilitates leader RNA switching and may enhance subgenomic mRNA transcription. It is unclear how the 9-nt sequence exerts these functions. In this study, we inserted the 9-nt sequence into a defective interfering (DI) RNA reporter system and demonstrated that mRNA transcription could be initiated from the 9-nt sequence almost as efficiently as from the intergenic sequence between genes 6 and 7. Sequence analysis of the mRNAs showed that the 9-nt sequence served as a site of fusion between the leaders and mRNA. The transcription initiation function of the 9-nt sequence could not be substituted by other 5'-terminal sequences. When the entire 5'-terminal sequence, including four copies of the UCUAA sequence plus the 9-nt sequence, was present, transcription could be initiated from any of the UCUAA copies or the 9-nt sequence, resulting in different copy numbers of the UCUAA sequence and the deletion of the 9-nt sequence in some mRNAs. All of these heterogeneous RNA species were also detected from the 5'-terminal region of the viral genomic-length RNA in MHV-infected cells. These results thus suggest tha the heterogeneity of the copy number of UCUAA sequences at the 5' end, the deletion of the 9-nt sequence in viral and DI RNAs, and the leader RNA switching are the results of transcriptional initiation from the 9-nt site. They also show that an mRNA species (mRNA 1) that lacks the 9-nt sequence can be synthesized during MHV infection. Therefore, MHV genomic RNA replication and mRNA 1 transcription may be distinguishable.

Investigation of the control of coronavirus subgenomic mRNA transcription by using T7-generated negative-sense RNA transcripts

The subgenomic mRNAs of the coronavirus transmissible gastroenteritis virus (TGEV) are not produced in equimolar amounts. We have developed a reporter gene system to investigate the control of this differential subgenomic mRNA synthesis. Transcription of mRNAs by the TGEV polymerase was obtained from negative-sense RNA templates generated in situ from DNA containing a T7 promoter. A series of gene cassettes was produced; these cassettes comprised the reporter chloramphenicol acetyltransferase (CAT) gene downstream of transcription-associated sequences (TASs) (also referred to as intergenic sequences and promoters) believed to be involved in the synthesis of TGEV subgenomic mRNAs 6 and 7. The gene cassettes were designed so that negative-sense RNA copies of the CAT gene with sequences complementary to the TGEV TASs, or modified versions, at the 3' end would be synthesized in situ by T7 RNA polymerase. Using this system, we have demonstrated that CAT was expressed from mRNAs derived from the T7-generated negative-sense RNA transcripts only in TGEV-infected cells and only from transcripts possessing a TGEV negative-sense TAS. Analysis of the CAT mRNAs showed the presence of the TGEV leader RNA sequence at the 5' end, in keeping with observations that all coronavirus mRNAs have a 5' leader sequence corresponding to the 5' end of the genomic RNA. Our results indicated that the CAT mRNAs were transcribed from the in situ-synthesized negative-sense RNA templates without the requirement of TGEV genomic 5' or 3' sequences on the T7-generated negative-sense transcripts (3'-TAS-CAT-5'). Modification of the TGEV TASs indicated (i) that the degree of potential base pairing between the 3' end of the leader RNA and the TGEV negative-sense TAS was not the sole determinant of the amount of subgenomic mRNA transcribed and (ii) that other factors, including nucleotides flanking the TAS, are involved in the regulation of transcription of TGEV subgenomic mRNAs.

Characterization of a murine coronavirus defective interfering RNA internal cis-acting replication signal

The mouse hepatitis virus (MHV) sequences required for replication of the JHM strain of MHV defective interfering (DI) RNA consist of three discontinuous genomic regions: about 0.47 kb from both terminal sequences and a 0.13-kb internal region present at about 0.9 kb from the 5' end of the DI genome. In this study, we investigated the role of the internal 0.13-kb region in MHV RNA replication. Overall sequences of the 0.13-kb regions from various MHV strains were similar to each other, with nucleotide substitutions in some strains; MHV-A59 was exceptional, with three nucleotide deletions. Computer-based secondary-structure analysis of the 0.13-kb region in the positive strand revealed that most of the MHV strains formed the same or a similar main stem-loop structure, whereas only MHV-A59 formed a smaller main stem-loop structure. The RNA secondary structures in the negative strands were much less uniform among the MHV strains. A series of DI RNAs that contained MHV-JHM-derived 5'- and 3'-terminal sequences plus internal 0.13-kb regions derived from various MHV strains were constructed. Most of these DI RNAs replicated in MHV-infected cells, except that MRP-A59, with a 0.13-kb region derived from MHV-A59, failed to replicate. Interestingly, replication of MRP-A59 was temperature dependent; it occurred at 39.5 degrees C but not at 37 or 35 degrees C, whereas a DI RNA with an MHV-JHM-derived 0.13-kb region replicated at all three temperatures. At 37 degrees C, synthesis of MRP-A59 negative-strand RNA was detected in MHV-infected and MRP-A59 RNA-transfected cells. Another DI RNA with the internal 0.13-kb region deleted also synthesized negative-strand RNA in MHV-infected cells. MRP-A59-transfected cells were shifted from 39.5 to 37 degrees C at 5.5 h postinfection, a time when most MHV negative-strand RNAs have already accumulated; after the shift, MRP-A59 positive-strand RNA synthesis ceased. The minimum sequence required for maintenance of the positive-strand major stem-loop structure and biological function of the MHV-JHM 0.13-kb region was about 57 nucleotides. Function was lost in the 50-nucleotide sequence that formed a positive-strand stem-loop structure identical to that of MHV-A59. These studies suggested that the RNA structure made by the internal sequence was important for positive-strand MHV RNA synthesis.

Two murine coronavirus genes suffice for viral RNA synthesis

We identified two mouse hepatitis virus (MHV) genes that suffice for MHV RNA synthesis by using an MHV-JHM-derived defective interfering (DI) RNA, DIssA. DIssA is a naturally occurring self-replicating DI RNA with nearly intact genes 1 and 7. DIssA interferes with most MHV-JHM-specific RNA synthesis, except for synthesis of mRNA 7, which encodes N protein; mRNA 7 synthesis is not inhibited by DIssA. Coinfection of MHV-JHM containing DIssA DI particles and an MHV-A59 RNA- temperature-sensitive mutant followed by subsequent passage of virus at the permissive temperature resulted in elimination of most of the MHV-JHM helper virus. Analysis of intracellular RNAs at the nonpermissive temperature demonstrated efficient synthesis of DIssA and mRNA 7 but not of the helper virus mRNAs. Oligonucleotide fingerprinting analysis demonstrated that the structure of mRNA 7 was MHV-JHM specific and therefore must have been synthesized from the DIssA template RNA. Sequence analysis revealed that DIssA lacks a slightly heterogeneous sequence, which is found in wild-type MHV from the 3' one-third of gene 2-1 to the 3' end of gene 6. Northern (RNA) blot analysis of intracellular RNA species and virus-specific protein analysis confirmed the sequence data. Replication and transcription of another MHV DI RNA were supported in DIssA-replicating cells. Because the products of genes 2 and 2-1 are not essential for MHV replication, we concluded that expression of gene 1 proteins and N protein was sufficient for MHV RNA replication and transcription.

Interactions between the cytoplasmic proteins and the intergenic (promoter) sequence of mouse hepatitis virus RNA: correlation with the amounts of subgenomic mRNA transcribed

Previous studies suggested that coronavirus RNA transcription involves interaction between leader RNA and the intergenic (IG) sequences, probably via protein-RNA interactions (X. M. Zhang, C.-L. Liao, and M. M. C. Lai, J. Virol., 68:4738-4746, 1994; X. M. Zhang and M. M. C. Lai, J. Virol., 68:6626-6633, 1994). To determine whether cellular proteins are involved in this process, we performed UV cross-linking experiments using cytoplasmic extracts of uninfected cells and the IG (promoter) sequence between genes 6 and 7 (IG7) and the 5' untranslational region of mouse hepatitis virus genomic RNA. We demonstrated that three different cellular proteins (p70, p48, and p35/38) bound to the promoter sequence of the template RNA. Deletion analyses of the template RNA mapped the binding site of p35/38 at the consensus transcription initiation signal. In contrast, the binding of p70 and p48 was less specific. p35/38 is the same protein as the one previously identified to bind to the complementary strand of the leader RNA; its binding affinity to the leader was approximately 15 times stronger than that to IG7. Site-directed mutagenesis of the IG sequence revealed that mutations in the consensus sequence of IG7 (UCUAAUCUAAAC to UCGAAAC and GCUAAAG), which resulted in reduced subgenomic mRNA transcription, also caused correspondingly reduced levels of p35/38 binding. These results demonstrated that the extent of protein binding to the IG sequences correlated with the amounts of subgenomic mRNAs transcribed from the IG site. These studies suggest that these RNA-binding proteins are involved in coronavirus RNA transcription and may represent transcription factors.

The effect of two closely inserted transcription consensus sequences on coronavirus transcription

Insertion of an intergenic region from the murine coronavirus mouse hepatitis virus into a mouse hepatitis virus defective interfering (DI) RNA led to transcription of subgenomic DI RNA in helper virus-infected cells. Using this system, we studied how two intergenic regions in close proximity affected subgenomic RNA synthesis. When two intergenic regions were separated by more than 100 nucleotides, slightly less of the larger subgenomic DI RNA (synthesized from the upstream intergenic region) was made; this difference was significant when the intergenic region separation was less than about 35 nucleotides. Deletion of sequences flanking the two intergenic regions inserted in close proximity did not affect transcription. No significant change in the ratio of the two subgenomic DI RNAs was observed when the sequence between the two intergenic regions was altered. Removal of the downstream intergenic region restored transcription of the larger subgenomic DI RNA. The UCUAAAC consensus sequence was needed for efficient suppression of the larger subgenomic DI RNA synthesis. These results demonstrated that the downstream intergenic sequence was suppressing subgenomic DI RNA synthesis from the upstream intergenic region. We discuss possible mechanisms to account for the regulation of this suppression of subgenomic DI RNA synthesis and the ways in which they relate to the general regulation of coronavirus transcription.

Regulation of transcription of coronaviruses

To study factors involved in regulation of transcription of coronaviruses, we constructed defective interfering (DI) RNAs containing sg RNA promoters at multiple positions. Analysis of the amounts of sg DI RNA produced by these DIs resulted in the following observations: (i) a downstream promoter downregulates an upstream promoter; (ii) an upstream promoter has little or no effect on the activity of a downstream promoter. Our data suggest that attenuation of upstream promoter activities by downstream promoter sequences plays an important role in regulating the amounts of sg RNAs produced by coronaviruses. Our observations are in accordance with the models proposed by Konings et al. and Sawicki and Sawicki.

Regulation of coronavirus RNA transcription is likely mediated by protein-RNA interactions

Coronavirus mRNA transcription was thought to be regulated by the interaction between the leader RNA and the intergenic (IG) sequence, probably involving direct RNA-RNA interactions between complementary sequences. In this study, we found that a 9-nucleotide sequence immediately downstream of the leader RNA up-regulated mRNA transcription and that a particular strain of mouse hepatitis virus (MHV) lacking this 9-nucleotide transcribed subgenomic mRNA species containing unusually heterogeneous leader-fusion sites. These results suggest that the sequence complementarity between the leader and IG is not necessarily required for mRNA transcription. UV cross-linking experiments using cytoplasmic extracts of uninfected cells and the IG sequence showed that three different cellular proteins bound to IG of the template RNA. Deletion analyses and site-directed mutagenesis of IG further demonstrated a correlation between protein-binding and transcription efficiency, suggesting that these RNA-binding proteins are involved in the regulation of coronavirus mRNA transcription. We propose that coronavirus transcription is regulated by RNA-protein and protein-protein interactions.

Detection of negative-stranded subgenomic RNAs but not of free leader in LDV-infected macrophages

The mechanism of synthesis of the seven subgenomic mRNAs of lactate dehydrogenase-elevating virus (LDV) was explored. One proposed mechanism, leader-primed transcription, predicts the formation of free 5'-leader in infected cells which then primes reinitiation of transcription at specific complementary sites on the antigenomic template. No free LDV 5'-leader of 156 nucleotides was detected in LDV-infected macrophages. Another mechanism, independent replication of the subgenomic mRNAs, predicts the presence of negative complements to all subgenomic mRNAs in infected cells which might be generated from subgenomic mRNAs in virions. Full-length antigenomic RNA was detected in LDV-infected macrophages by Northern hybridization at a level of < 1% of that of genomic RNA, but no negative polarity subgenomic RNAs. Negative complements to all subgenomic mRNAs, however, were detected by reverse transcription of total RNA from infected macrophages using as primer an oligonucleotide complementary to the antileader followed by polymerase chain reaction amplification using this sense primer in combination with various oligonucleotide primers complementary to a segment downstream of the junction between the 5' leader and the body of each subgenomic RNA. It is unclear whether these minute amounts of negative subgenomic RNAs function in the replication of the subgenomic mRNAs. They could also be by-products of the RNA replication process. Finally, no subgenomic mRNAs were detected in LDV virions.

Coronavirus leader RNA regulates and initiates subgenomic mRNA transcription both in trans and in cis

Mouse hepatitis virus (MHV), a coronavirus, utilizes a discontinuous transcription mechanism for subgenomic mRNA synthesis. Previous studies (C.-L. Liao and M. C. C. Lai, J. Virol. 68:4727-4737, 1994) have demonstrated that an upstream cis-acting leader sequence serves as a transcriptional enhancer, but the mechanism of transcriptional regulation is not clear. In this study, we constructed a series of defective interfering (DI) RNAs containing the chloramphenicol acetyltransferase (CAT) gene behind a differentially expressed transcription initiation (intergenic) sequence (for mRNA2-1). These DI RNAs had different copy numbers of the UCUAA pentanucleotide sequence at the 3' end of the leader. Transfection of these DI RNA constructs into cells infected with a helper MHV, which contains either two or three UCUAA copies at the 3' end of the leader, resulted in differential expression of CAT activities. We demonstrated that the copy number of UCUAA repeats in the leaders of both helper viral and DI RNAs affected the level of CAT activity, suggesting that MHV leader RNA could regulate both in trans and in cis the transcription of subgenomic mRNAs. The leader RNA of subgenomic mRNAs was derived from either the trans- or the cis-acting leader. Furthermore, insertion of a UA-rich sequence (UUUAUAAAC) immediately downstream of the leader in DI RNA, to match the sequence of helper viral RNA, enhanced the CAT activity by threefold, suggesting that this nine-nucleotide sequence is a cis-acting element. Interestingly, when the nine-nucleotide sequence was absent in DI RNA, the leaders of subgenomic mRNAs were exclusively derived from the helper virus. In contrast, when the nine-nucleotide sequence was present in DI RNA, the leaders were derived from both helper viral and DI RNAs. These results suggest that the nine-nucleotide sequence either is required for the leader RNA to initiate mRNA synthesis or, alternatively, serves as a transcription terminator for the leader RNA synthesis. However, when a constitutively expressed intergenic sequence (for mRNA7) was used, no difference in transcription efficiency was noted, regardless of the copy number of UCUAA in the DI RNA and helper virus. This study thus indicates that MHV subgenomic RNA transcription requires the interaction among the intergenic (promoter) sequence, a trans-acting leader, and a cis-acting leader sequence. A novel model of transcriptional regulation of coronavirus subgenomic mRNAs is presented.

Subgenomic RNA synthesis directed by a synthetic defective interfering RNA of mouse hepatitis virus: a study of coronavirus transcription initiation

We have used a full-length cDNA clone of a mouse hepatitis virus strain A59 defective interfering (DI) RNA, pMIDI-C, and cassette mutagenesis to study the mechanism of coronavirus subgenomic mRNA synthesis. Promoter sequences closely resembling those of subgenomic mRNAs 3 and 7 were inserted into MIDI-C. Both subgenomic RNA promoters gave rise to the synthesis of a subgenomic DI RNA in virus-infected and DI RNA-transfected cells. From a mutagenic analysis of the promoters we concluded the following. (i) The extent of base pairing between the leader RNA and the intergenic promoter sequence does not control subgenomic RNA abundance. (ii) Promoter recognition does not rely on base pairing only. Presumably, transcription initiation requires recognition of the promoter sequence by the transcriptase. (iii) Fusion of leader and body sequences takes place at multiple--possibly random--sites within the intergenic promoter sequence. A model is presented in which, prior to elongation, the leader RNA is trimmed by a processive 3'-->5' nuclease.

Evidence for coronavirus discontinuous transcription

Coronavirus subgenomic mRNA possesses a 5'-end leader sequence which is derived from the 5' end of genomic RNA and is linked to the mRNA body sequence. This study examined whether coronavirus transcription involves a discontinuous transcription step; the possibility that a leader sequence from mouse hepatitis virus (MHV) genomic RNA could be used for MHV subgenomic defective interfering (DI) RNA transcription was examined. This was tested by using helper viruses and DI RNAs that were easily distinguishable. MHV JHM variant JHM(2), which synthesizes a subgenomic mRNA encoding the HE gene, and variant JHM(3-9), which does not synthesize this mRNA, were used. An MHV DI RNA, DI(J3-9), was constructed to contain a JHM(3-9)-derived leader sequence and an inserted intergenic region derived from the region preceding the MHV JHM HE gene. DI(J3-9) replicated efficiently in JHM(2)- or JHM(3-9)-infected cells, whereas synthesis of subgenomic DI RNAs was observed only in JHM(2)-infected cells. Sequence analyses demonstrated that the 5' regions of both helper virus genomic RNAs and genomic DI RNAs maintained their original sequences in DI RNA-replicating cells, indicating that the genomic leader sequences derived from JHM(2) functioned for subgenomic DI RNA transcription. Replication and transcription of DI(J3-9) were observed in cells infected with an MHV A59 strain whose leader sequence was similar to that of JHM(2), except for one nucleotide substitution within the leader sequence. The 5' region of the helper virus genomic RNA and that of the DI RNA were the same as their original structures in virus-infected cells, and the leader sequence of DI(J3-9) subgenomic DI RNA contained the MHV A59-derived leader sequence. The leader sequence of subgenomic DI RNA was derived from that of helper virus; therefore, the genomic leader sequence had a trans-acting property indicative of a discontinuous step in coronavirus transcription.

Genome organization of porcine epidemic diarrhoea virus

In order to study the organization of the genome of porcine epidemic diarrhoea virus (PEDV), we constructed a cDNA library in a phage expression vector by using poly(A) RNA from PEDV-infected Vero cells. An anti-PEDV hyperimmune serum was used to probe the library. The first isolated clone mapped within the N gene and was subsequently used for rescreening the library. The selected clones allowed us to establish the sequence of the 3'-most 7.4 kb of the PEDV genome. Analysis of the cDNA sequences revealed a 3'-coterminal nested structure, which is typical of Coronaviridae and the presence of a hexameric sequence XUA(A/G)AC upstream of each coding region. The amino acid sequences deduced from four of the five ORFs identified showed the characteristic features of the structural proteins S, M, sM and N. Only one ORF located between the S and M genes was found to potentially encode a non-structural polypeptide. Our data lead us to conclude that PEDV is a member of Coronaviridae and belongs to the same genetic subset as TGEV, FIPV and HCV 229E.

Three different cellular proteins bind to complementary sites on the 5'-end-positive and 3'-end-negative strands of mouse hepatitis virus RNA

The termini of viral genomic RNA and its complementary strand are important in the initiation of viral RNA replication, which probably involves both viral and cellular proteins. To detect the possible cellular proteins involved in the replication of mouse hepatitis virus RNA, we performed RNA-protein binding studies with RNAs representing both the 5' and 3' ends of the viral genomic RNA and the 3' end of the negative-strand complementary RNA. Gel-retardation assays showed that both the 5'-end-positive- and 3'-end-negative-strand RNA formed an RNA-protein complex with cellular proteins from the uninfected cells. UV cross-linking experiments further identified a 55-kDa protein bound to the 5' end of the positive-strand viral genomic RNA and two proteins 35 and 38 kDa in size bound to the 3' end of the negative-strand cRNA. The results of the competition assay confirmed the specificity of this RNA-protein binding. No proteins were found to bind to the 3' end of the viral genomic RNA under the same conditions. The binding site of the 55-kDa protein was mapped within the 56-nucleotide region from nucleotides 56 to 112 from the 5' end of the positive-strand RNA, and the 35- and 38-kDa proteins bound to the complementary region on the negative-strand RNA. The 38-kDa protein was detected only in DBT cells but was not detected in HeLa or COS cells, while the 35-kDa protein was found in all three cell types. The juxtaposition of the different cellular proteins on the complementary sites near the ends of the positive- and negative-strand RNAs suggests that these proteins may interact with each other and play a role in mouse hepatitis virus RNA replication.

Effect of intergenic consensus sequence flanking sequences on coronavirus transcription

Insertion of a region, including the 18-nucleotide-long intergenic sequence between genes 6 and 7 of mouse hepatitis virus (MHV) genomic RNA, into an MHV defective interfering (DI) RNA leads to transcription of subgenomic DI RNA in helper virus-infected cells (S. Makino, M. Joo, and J. K. Makino, J. Virol. 66:6031-6041, 1991). In this study, the subgenomic DI RNA system was used to determine how sequences flanking the intergenic region affect MHV RNA transcription and to identify the minimum intergenic sequence required for MHV transcription. DI cDNAs containing the intergenic region between genes 6 and 7, but with different lengths of upstream or downstream flanking sequences, were constructed. All DI cDNAs had an 18-nucleotide-long intergenic region that was identical to the 3' region of the genomic leader sequence, which contains two UCUAA repeat sequences. These constructs included 0 to 1,440 nucleotides of upstream flanking sequence and 0 to 1,671 nucleotides of downstream flanking sequence. An analysis of intracellular genomic DI RNA and subgenomic DI RNA species revealed that there were no significant differences in the ratios of subgenomic to genomic DI RNA for any of the DI RNA constructs. DI cDNAs which lacked the intergenic region flanking sequences and contained a series of deletions within the 18-nucleotide-long intergenic sequence were constructed to determine the minimum sequence necessary for subgenomic DI RNA transcription. Small amounts of subgenomic DI RNA were synthesized from genomic DI RNAs with the intergenic consensus sequences UCUAAAC and GCUAAAC, whereas no subgenomic DI RNA transcription was observed from DI RNAs containing UCUAAAG and GCTAAAG sequences. These analyses demonstrated that the sequences flanking the intergenic sequence between genes 6 and 7 did not play a role in subgenomic DI RNA transcription regulation and that the UCUAAAC consensus sequence was sufficient for subgenomic DI RNA transcription.