Coronavirus mRNA synthesis: identification of novel transcription initiation signals which are differentially regulated by different leader sequences

Comparative analysis of bioinformatics tools to characterize SARS-CoV-2 subgenomic RNAs

During the replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), positive-sense genomic RNA and subgenomic RNAs (sgRNAs) are synthesized by a discontinuous process of transcription characterized by a template switch, regulated by transcription-regulating sequences (TRS). Although poorly known about makeup and dynamics of sgRNAs population and function of its constituents, next-generation sequencing approaches with the help of bioinformatics tools have made a significant contribution to expand the knowledge of sgRNAs in SARS-CoV-2. For this scope to date, Periscope, LeTRS, sgDI-tector, and CORONATATOR have been developed. However, limited number of studies are available to compare the performance of such tools. To this purpose, we compared Periscope, LeTRS, and sgDI-tector in the identification of canonical (c-) and noncanonical (nc-) sgRNA species in the data obtained with the Illumina ARTIC sequencing protocol applied to SARS-CoV-2-infected Caco-2 cells, sampled at different time points. The three software showed a high concordance rate in the identification and in the quantification of c-sgRNA, whereas more differences were observed in nc-sgRNA. Overall, LeTRS and sgDI-tector result to be adequate alternatives to Periscope to analyze Fastq data from sequencing platforms other than Nanopore.

sgRNAs: A SARS-CoV-2 emerging issue

Like for other coronaviruses, SARS-CoV-2 gene expression strategy is based on the synthesis of a nested set of subgenomic mRNA species (sgRNAs). These sgRNA are synthesized using a "discontinuous transcription" mechanism that relies on template switching at Transcription Regulatory Sequences (TRS). Both canonical (c-sgRNA) and non-canonical (nc-sgRNA, less numerous) subgenomic RNA species can be produced. Currently, sgRNAs are investigated on the basis of sequence data obtained through next generation sequencing (NGS), and bioinformatic tools are crucial for their identification, characterization and quantification. To date, few software have been developed to this aim, whose reliability and applicability to all the available NGS platforms need to be established, to build confidence on the information resulting from such tools. In fact, these information may be crucial for the in depth elucidation of viral expression strategy, particularly in respect of the significance of nc-sgRNAs, and for the possible use of sgRNAs as potential markers of virus replicative activity in infected patients.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

New insights about the regulation of Nidovirus subgenomic mRNA synthesis

The members of the Order Nidovirales share a similar genome organization with two overlapping nonstructural polyproteins encoded in the 5' two-thirds and the structural proteins encoded in the 3' third. They also express their 3' region proteins from a nested set of 3' co-terminal subgenomic messenger RNAs (sg mRNAs). Some but not all of the Nidovirus sg mRNAs also have a common 5' leader sequence that is acquired by a discontinuous RNA synthesis mechanism regulated by multiple 3' body transcription regulating sequences (TRSs) and the 5' leader TRS. Initial studies detected a single major body TRS for each 3' sg mRNA with a few alternative functional TRSs reported. The recent application of advanced techniques, such as next generation sequencing and ribosomal profiling, in studies of arteriviruses and coronaviruses has revealed an expanded sg mRNA transcriptome and coding capacity.

Expanded subgenomic mRNA transcriptome and coding capacity of a nidovirus

Members of the order Nidovirales express their structural protein ORFs from a nested set of 3' subgenomic mRNAs (sg mRNAs), and for most of these ORFs, a single genomic transcription regulatory sequence (TRS) was identified. Nine TRSs were previously reported for the arterivirus Simian hemorrhagic fever virus (SHFV). In the present study, which was facilitated by next-generation sequencing, 96 SHFV body TRSs were identified that were functional in both infected MA104 cells and macaque macrophages. The abundance of sg mRNAs produced from individual TRSs was consistent over time in the two different cell types. Most of the TRSs are located in the genomic 3' region, but some are in the 5' ORF1a/1b region and provide alternative sources of nonstructural proteins. Multiple functional TRSs were identified for the majority of the SHFV 3' ORFs, and four previously identified TRSs were found not to be the predominant ones used. A third of the TRSs generated sg mRNAs with variant leader-body junction sequences. Sg mRNAs encoding E', GP2, or ORF5a as their 5' ORF as well as sg mRNAs encoding six previously unreported alternative frame ORFs or 14 previously unreported C-terminal ORFs of known proteins were also identified. Mutation of the start codon of two C-terminal ORFs in an infectious clone reduced virus yield. Mass spectrometry detected one previously unreported protein and suggested translation of some of the C-terminal ORFs. The results reveal the complexity of the transcriptional regulatory mechanism and expanded coding capacity for SHFV, which may also be characteristic of other nidoviruses.

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.

Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus

Coronaviruses are a family of enveloped, single-stranded, positive-strand RNA viruses classified within the Nidovirales order. This coronavirus family consists of pathogens of many animal species and of humans, including the recently isolated severe acute respiratory syndrome coronavirus (SARS-CoV). This review is divided into two main parts; the first concerns the animal coronaviruses and their pathogenesis, with an emphasis on the functions of individual viral genes, and the second discusses the newly described human emerging pathogen, SARS-CoV. The coronavirus part covers (i) a description of a group of coronaviruses and the diseases they cause, including the prototype coronavirus, murine hepatitis virus, which is one of the recognized animal models for multiple sclerosis, as well as viruses of veterinary importance that infect the pig, chicken, and cat and a summary of the human viruses; (ii) a short summary of the replication cycle of coronaviruses in cell culture; (iii) the development and application of reverse genetics systems; and (iv) the roles of individual coronavirus proteins in replication and pathogenesis. The SARS-CoV part covers the pathogenesis of SARS, the developing animal models for infection, and the progress in vaccine development and antiviral therapies. The data gathered on the animal coronaviruses continue to be helpful in understanding SARS-CoV.

Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus

Coronaviruses are a family of enveloped, single-stranded, positive-strand RNA viruses classified within the Nidovirales order. This coronavirus family consists of pathogens of many animal species and of humans, including the recently isolated severe acute respiratory syndrome coronavirus (SARS-CoV). This review is divided into two main parts; the first concerns the animal coronaviruses and their pathogenesis, with an emphasis on the functions of individual viral genes, and the second discusses the newly described human emerging pathogen, SARS-CoV. The coronavirus part covers (i) a description of a group of coronaviruses and the diseases they cause, including the prototype coronavirus, murine hepatitis virus, which is one of the recognized animal models for multiple sclerosis, as well as viruses of veterinary importance that infect the pig, chicken, and cat and a summary of the human viruses; (ii) a short summary of the replication cycle of coronaviruses in cell culture; (iii) the development and application of reverse genetics systems; and (iv) the roles of individual coronavirus proteins in replication and pathogenesis. The SARS-CoV part covers the pathogenesis of SARS, the developing animal models for infection, and the progress in vaccine development and antiviral therapies. The data gathered on the animal coronaviruses continue to be helpful in understanding SARS-CoV.

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.

Coronavirus replication and pathogenesis: Implications for the recent outbreak of severe acute respiratory syndrome (SARS), and the challenge for vaccine development

A novel coronavirus has been recently identified as the causative agent of the severe acute respiratory syndrome (SARS) outbreak that has accounted for more than 8000 infected people worldwide. This review will discuss current knowledge on coronavirus replication, pathogenesis, evolution, and vaccine strategies, as well as the most recent findings on SARS coronavirus.

Down-regulation of transcription of the proapoptotic gene BNip3 in cultured astrocytes by murine coronavirus infection

Murine coronavirus mouse hepatitis virus (MHV) causes encephalitis and demyelination in the central nervous system of susceptible rodents. Astrocytes are the major target for MHV persistence. However, the mechanisms by which astrocytes survive MHV infection and permit viral persistence are not known. Here we performed DNA microarray analysis on differential gene expression in astrocyte DBT cells by MHV infection and found that the mRNA of the proapoptotic gene BNip3 was significantly decreased following MHV infection. This finding was further confirmed by quantitative reverse transcription-polymerase chain reaction, Western blot analysis, and BNip3-promoter-luciferase reporter system. Interestingly, infection with live and ultraviolet light-inactivated viruses equally repressed BNip3 expression, indicating that the down-regulation of BNip3 expression does not require virus replication and is mediated during cell entry. Furthermore, treatment of cells with chloroquine, which blocks the acidification of endosomes, significantly inhibited the repression of the BNip3 promoter activity induced by the acidic pH-dependent MHV mutant OBLV60, which enters cells via endocytosis, indicating that the down-regulation of BNip3 expression is mediated by fusion between viral envelope and cell membranes during entry. Deletion analysis showed that the sequence between nucleotides 262 and 550 of the 588-base-pair BNip3 promoter is necessary and sufficient for driving the BNip3 expression and that it contains signals that are responsible for MHV-induced down-regulation of BNip3 expression in DBT cells. These results may provide insights into the mechanisms by which MHV evades host antiviral defense and promotes cell survival, thereby allowing its persistence in the host astrocytes.

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.

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.

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.

Expression of hemagglutinin/esterase by a mouse hepatitis virus coronavirus defective-interfering RNA alters viral pathogenesis

A defective-interfering (DI) RNA of mouse hepatitis virus (MHV) was developed as a vector for expressing MHV hemagglutinin/esterase (HE) protein. The virus containing an expressed HE protein (A59-DE-HE) was generated by infecting cells with MHV-A59, which does not express HE, and transfecting the in vitro-transcribed DI RNA containing the HE gene. A similar virus (A59-DE-CAT) expressing the chloramphenicol acetyltransferase (CAT) was used as a control. These viruses were inoculated intracerebrally into mice, and the role of the HE protein in viral pathogenesis was evaluated. Results showed that all mice infected with parental A59 or A59-DE-CAT succumbed to infection by 9 days postinfection (p.i.), demonstrating that inclusion of the DI did not by itself alter pathogenesis. In contrast, 60% of mice infected with A59-DE-HE survived infection. HE- or CAT-specific subgenomic mRNAs were detected in the brains at days 1 and 2 p.i. but not later, indicating that the genes in the DI vector were expressed only in the early stage of viral infection. No significant difference in virus titer or viral antigen expression in brains was observed between A59-DE-HE- and A59-DE-CAT-infected mice, suggesting that virus replication in brain was not affected by the expression of HE. However, at day 3 p.i. there was a slight increase in the extent of inflammatory cell infiltration in the brains of the A59-DE-HE-infected mice. Surprisingly, virus titers in the livers of A59-DE-HE-infected mice were 3 log10 lower than that of the A59-DE-CAT-infected mice at day 6 p.i. Also, substantially less necrosis and viral antigen were detected in the livers of the A59-DE-HE-infected mice. This may account for the reduced mortality of these mice. The possible contribution of the host immune system to this difference in pathogenesis was analyzed by comparing the expression of four cytokines. Results showed that both tumor necrosis factor-alpha and interleukin-6 mRNAs increased in the brains of the A59-DE-HE-infected mice at day 2 p.i., whereas interferon-gamma and interleukin-1 alpha mRNAs were similar between A59-DE-HE- and A59-DE-CAT-infected mice. These data suggest that the transient expression of HE protein enhances an early innate immune response, possibly contributing to the eventual clearance of virus from the liver. This study indicates the feasibility of the DI expression system for studying roles of viral proteins during MHV infection.

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.

Analysis of a recombinant mouse hepatitis virus expressing a foreign gene reveals a novel aspect of coronavirus transcription

We have inserted heterologous genetic material into the nonessential gene 4 of the coronavirus mouse hepatitis virus (MHV) in order to test the applicability of targeted RNA recombination for site-directed mutagenesis of the MHV genome upstream of the nucleocapsid (N) gene and to develop further genetic tools for site-directed mutagenesis of structural genes other than N. Initially, a 19-nucleotide tag was inserted into the start of gene 4a of MHV strain A59 with the N gene deletion mutant Alb4 as the recipient virus. In further work, the entire gene for the green fluorescent protein (GFP) was inserted in place of gene 4, creating the currently largest known RNA virus. The expression of GFP was demonstrated by Western blot analysis of infected cell lysates; however, the level of GFP expression was not sufficient to allow detection of fluorescence of viral plaques. Northern blot analysis of transcripts of GFP recombinants showed the expected alteration of the pattern of the nested MHV subgenomic mRNAs. Surprisingly, though, GFP recombinants also produced an RNA species that was the same size as wild-type mRNA4. Analysis of the 5' end of this species revealed that it was actually a collection of mRNAs originating from 10 different genomic fusion sites, none possessing a canonical intergenic sequence. The finding of these aberrant mRNAs suggests that long-range RNA or the ribonucleoprotein structure of the MHV genome can sometimes be the sole determinant of the site of initiation of transcription.

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

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.

Genetics of mouse hepatitis virus transcription: evidence that subgenomic negative strands are functional templates

Mouse hepatitis virus (MHV) A59 temperature-sensitive (ts) mutants belonging to complementation group C were characterized and mapped by standard genetic recombination techniques. Temperature shift experiments early in infection suggested that the group C allele can be divided into two phenotypically distinct subgroups, designated C1 and C2. Since previous data indicated that the group C1 mutants probably contained an early defect which affects negative-strand synthesis, RNA synthesis was further examined by analyzing replicative-form (RF) RNA. Full-length as well as subgenomic-length RF RNAs were radiolabeled from 3 to 12 h postinfection (p.i.) and labeled late in infection after shift to the nonpermissive temperature (39.5 degrees C). The relative percent molar ratios of each mRNA and corresponding RF RNA were roughly equivalent throughout infection. Temperature shift experiments at 5.5 or 6.0 h p.i. resulted in an 83 to 92% reduction in the amount of total RF RNA at 39.5 degrees C. Radiolabeling time course experiments after temperature shift to 39.5 degrees C also demonstrated incorporation (6 to 9 h p.i.) into both subgenomic-length and full-length RF RNAs, suggesting that previously transcribed negative strands were functional templates throughout infection. To determine if the reduction in RF RNA was due to a decrease in positive- or negative-strand RNA synthesis, rates of mRNA synthesis were calculated from both full-length and subgenomic-length templates. The rate of mRNA synthesis after the shift was increased at 39.5 degrees C compared with that at 32 degrees C regardless of the template used; however, transcription rates calculated from subgenomic-length templates were similar to those of other viral and eukaryotic polymerases. These findings support the notion that the group C1 allele regulates negative-strand RNA synthesis and strongly suggest that the subgenomic negative-strand RNAs are probably the predominant functional templates for the synthesis of positive-strand RNAs late in infection.

Requirement of the 5'-end genomic sequence as an upstream cis-acting element for coronavirus subgenomic mRNA transcription

We have developed a defective interfering (DI) RNA containing a chloramphenicol acetyltransferase reporter gene, placed behind an intergenic sequence, for studying subgenomic mRNA transcription of mouse hepatitis virus (MHV), a prototype coronavirus. Using this system, we have identified the sequence requirement for MHV subgenomic mRNA transcription. We show that this sequence requirement differs from that for RNA replication. In addition to the previously identified requirement for an intergenic (promoter) sequence, additional sequences from the 5' end of genomic RNA are required for subgenomic mRNA transcription. These upstream sequences include the leader RNA and a spacer sequence between the leader and intergenic sequence, which is derived from the 5' untranslated region and part of gene 1. The spacer sequence requirement is specific, since only the sequence derived from the 5' end of RNA genome, but not from other MHV genomic regions or heterologous sequences, could initiate subgenomic transcription from the intergenic sequence. These results strongly suggest that the wild-type viral subgenomic mRNAs (mRNA2 to mRNA7) and probably their counterpart subgenomic negative-sense RNAs cannot be utilized for mRNA amplification. Furthermore, we have demonstrated that a partial leader sequence present at the 5' end of genome, which lacks the leader-mRNA fusion sequence, could still support subgenomic mRNA transcription. In this case, the leader sequences of the subgenomic transcripts were derived exclusively from the wild-type helper virus, indicating that the MHV leader RNA initiates in trans subgenomic mRNA transcription. Thus, the leader sequence can enhance subgenomic transcription even when it cannot serve as a primer for mRNA synthesis. These results taken together suggest that the 5'-end leader sequence of MHV not only provides a trans-acting primer for mRNA initiation but also serves as a cis-acting element required for the transcription of subgenomic mRNAs. The identification of an upstream cis-acting element for MHV subgenomic mRNA synthesis defines a novel sequence requirement for regulating mRNA synthesis in RNA viruses.

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.

Studies into the mechanism for MHV transcription

Previous studies have demonstrated that the MHV genome is divided into seven transcriptional units which are transcribed from highly conserved intergenic start sites (UCU/CAAAC) into mRNA containing a common leader RNA at the 5' end and a coterminal 3' end. In this manuscript, we provide evidence that an additional transcriptional unit is encoded at the 3' end of the MHV genome and is transcribed from a perfect intergenic region into a leader-containing approximately 800 nt mRNA. This mRNA could potentially encode a small 17-18 kDa protein which is identical to the C-terminal third of the nucleocapsid gene.

Mutagenic analysis of the coronavirus intergenic consensus sequence

Previously, a system in which an intergenic region from mouse hepatitis virus (MHV) inserted into an MHV defective interfering (DI) RNA led to transcription of a subgenomic DI RNA in helper virus-infected cells was established. In the present study, a DI cDNA containing one UCUAAAC consensus sequence in the middle of the 0.3-kb-long intergenic region located between genes 6 and 7 was constructed. From this DI cDNA clone, 21 mutant DI RNAs were constructed so that each of the seven consensus sequence nucleotides was changed individually to the three alternative bases. These mutants were used to define how changes in the integrity of MHV transcription consensus sequence UCUAAAC affected mRNA transcription. Except for two mutants with the sequences UGUAAAC and UCGAAAC, all of the mutants supported efficient subgenomic DI RNA transcription. This indicated that MHV transcription regulation was sufficiently flexible to recognize altered consensus sequences. Next, these and other mutants were used to examine the leader-body fusion site on the subgenomic DI RNAs. Sequence analysis demonstrated that all subgenomic DI RNAs analyzed contained two pentanucleotide sequences; the first sequence seemed to be contributed by the leader, and the leader-body fusion most likely took place at either the first or the second nucleotide of the second sequence. This observation was not consistent with the proposed coronavirus transcription model (S. C. Baker and M. M. C. Lai, EMBO J. 9:4173-4179, 1990) which states that nucleotide mismatch can be corrected by RNA polymerase proofreading activity.

Coronavirus mRNA transcription: UV light transcriptional mapping studies suggest an early requirement for a genomic-length template

Mouse hepatitis virus (MHV) synthesizes seven to eight mRNAs, each of which contains a leader RNA derived from the 5' end of the genome. To understand the mechanism of synthesis of these mRNAs, we studied how the synthesis of each mRNA was affected by UV irradiation at different time points after infection. When MHV-infected cells were UV irradiated at a late time in infection (5 h postinfection), the syntheses of the various mRNAs were inhibited to different extents in proportion to the sizes of the mRNAs. Analysis of the UV inactivation kinetics revealed that the UV target size of each mRNA was equivalent to its own physical size. In contrast, when cells were irradiated at 2.5 or 3 h postinfection, there appeared to be two different kinetics of inhibition of mRNA synthesis: the synthesis of every mRNA was inhibited to the same extent by a small UV dose, but the remaining mRNA synthesis was inhibited by additional UV doses at different rates for different mRNAs in proportion to RNA size. The analysis of the UV inactivation kinetics indicated that the UV target sizes for the majority of mRNAs were equivalent to that of the genomic-size RNA early in the infection. These results suggest that MHV mRNA synthesis requires the presence of a genomic-length RNA template at least early in the infection. In contrast, later in the infection, the sizes of the templates used for mRNA synthesis were equivalent to the physical sizes of each mRNA. The possibility that the genomic-length RNA required early in the infection was used only for the synthesis of a polymerase rather than as a template for mRNA synthesis was ruled out by examining the UV sensitivity of a defective interfering (DI) RNA. We found that the UV target size for the DI RNA early in infection was much smaller than that for mRNAs 6 and 7, which are approximately equal to or smaller in size than the DI RNA. This result indicates that even though DI RNA and viral mRNAs are synthesized by the same polymerase, mRNAs are synthesized from a larger (genomic-length) template. We conclude that a genomic-length RNA template is required for MHV subgenomic mRNA synthesis at least early in infection. Several transcription models are proposed.