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

Diverse effects of coronavirus-defective viral genomes on the synthesis of IFNβ and ISG15 mRNAs and coronavirus replication

BACKGROUND

The mechanism by which coronavirus-defective viral genomes (DVGs) affect coronavirus and host cells during infection remains unclear. A variety of DVGs with different RNA structures can be synthesized from coronavirus-infected cells, and these DVGs can also encode proteins. Consequently, in the present study, we first dissected the effects of individual DVGs on the synthesis of IFNβ and ISG15 mRNAs at the RNA, protein and combined levels, and then examined whether different coronavirus-DVGs have different effects on the synthesis of IFNβ and ISG15 mRNAs and coronavirus replication both individually and collectively under different infection conditions.

METHODS

To dissect the effects of individual DVGs on the synthesis of IFNβ and ISG15 mRNAs at the RNA, protein and combined levels, DVG 2.2 and DVG 5.1, which were previously identified in coronavirus-infected cells, and their mutants were constructed followed by transfection. Western blot and RT‒qPCR were used to detect the synthesis of protein and to quantify the synthesis of IFNβ and ISG15 mRNAs, respectively. To examined whether different coronavirus-DVGs have different effects on the synthesis of IFNβ and ISG15 mRNAs and coronavirus replication both individually and collectively under different infection conditions, different naturally occurring DVGs were selected and constructed followed by transfection after or before coronavirus infection and by RT‒qPCR and hemagglutination assay.

RESULTS

These results suggested that (i) coronavirus-DVGs at the RNA, protein and combined levels have different effects on the synthesis of IFNβ and ISG15 mRNAs, (ii) coronavirus-DVGs can inhibit coronavirus replication at least partly through interferon signaling and (iii) different DVGs have different effects on the synthesis of IFNβ and ISG15 mRNAs and coronavirus replication both individually and collectively under different infection conditions.

CONCLUSIONS

Coronavirus replication can be regulated by diverse coronavirus-derived DVGs at least partly through innate immunity. Such regulation may contribute to the pathogenesis of coronavirus. The DVG populations in coronavirus-infected cells with the ability to inhibit coronavirus replication are expected to be potential resources for the identification of antivirals at the level of RNA, protein or in combination, and the methods used in the current study can be used as a platform for this purpose.

Preferential cleavage of the coronavirus defective viral genome by cellular endoribonuclease with characteristics of RNase L

In testing whether coronavirus defective viral genome 12.7 (DVG12.7) with transcription regulating sequence (TRS) can synthesize subgenomic mRNA (sgmRNA) in coronavirus-infected cells, it was unexpectedly found by Northern blot assay that not only sgmRNA (designated sgmDVG 12.7) but also an RNA fragment with a size less than sgmDVG 12.7 was identified. A subsequent study demonstrated that the identified RNA fragment (designated clvDVG) was a cleaved RNA product originating from DVG12.7, and the cleaved sites were located in the loop region of stem‒loop structure and after UU and UA dinucleotides. clvDVG was also identified in mock-infected HRT-18 cells transfected with DVG12.7 transcript, indicating that cellular endoribonuclease is responsible for the cleavage. In addition, the sequence and structure surrounding the cleavage sites can affect the cleavage efficiency of DVG12.7. The cleavage features are therefore consistent with the general criteria for RNA cleavage by cellular RNase L. Furthermore, both the cleavage of rRNA and the synthesis of clvDVG were also identified in A549 cells. Because (i) the cleavage sites occurred predominantly after single-stranded UA and UU dinucleotides, (ii) the sequence and structure surrounding the cleavage sites affected the cleavage efficiency, (iii) the cleavage of rRNA is an index of the activation of RNase L, and (iv) the cleavage of both rRNA and DVG12.7 was identified in A549 cells, the results together indicated that the preferential cleavage of DVG12.7 is correlated with cellular endoribonuclease with the characteristics of RNase L and such cleavage features have not been previously characterized in coronaviruses.

Identification of subgenomic mRNAs derived from the coronavirus 1a/1b protein gene: Implications for coronavirus transcription

Synthesis of coronavirus subgenomic mRNA (sgmRNA) is guided by the transcription regulatory sequence (TRS). sgmRNA derived from the body TRS (TRS-B) located at the 1a/1b protein gene is designated 1ab/sgmRNA. In the current study, we comprehensively identified the 1ab/sgmRNAs synthesized from TRS-Bs located at the 1a/1b protein genes of different coronavirus genera both in vitro and in vivo by RT‒PCR and sequencing. The results suggested that the degree of sequence homology between the leader TRS (TRS-L) and TRS-B may not be a decisive factor for 1ab/sgmRNA synthesis. This observation led us to revisit the coronavirus transcription mechanism and to propose that the disassociation of coronavirus polymerase from the viral genome may be a prerequisite for sgmRNA synthesis. Once the polymerase can disassociate at TRS-B, the sequence homology between TRS-L and TRS-B is important for sgmRNA synthesis. The study therefore extends our understanding of transcription mechanisms.

Biological characterization of coronavirus noncanonical transcripts in vitro and in vivo

BACKGROUND

In addition to the well-known coronavirus genomes and subgenomic mRNAs, the existence of other coronavirus RNA species, which are collectively referred to as noncanonical transcripts, has been suggested; however, their biological characteristics have not yet been experimentally validated in vitro and in vivo.

METHODS

To comprehensively determine the amounts, species and structures of noncanonical transcripts for bovine coronavirus in HRT-18 cells and mouse hepatitis virus A59, a mouse coronavirus, in mouse L cells and mice, nanopore direct RNA sequencing was employed. To experimentally validate the synthesis of noncanonical transcripts under regular infection, Northern blotting was performed. Both Northern blotting and nanopore direct RNA sequencing were also applied to examine the reproducibility of noncanonical transcripts. In addition, Northern blotting was also employed to determine the regulatory features of noncanonical transcripts under different infection conditions, including different cells, multiplicities of infection (MOIs) and coronavirus strains.

RESULTS

In the current study, we (i) experimentally determined that coronavirus noncanonical transcripts were abundantly synthesized, (ii) classified the noncanonical transcripts into seven populations based on their structures and potential synthesis mechanisms, (iii) showed that the species and amounts of the noncanonical transcripts were reproducible during regular infection but regulated in altered infection environments, (iv) revealed that coronaviruses may employ various mechanisms to synthesize noncanonical transcripts, and (v) found that the biological characteristics of coronavirus noncanonical transcripts were similar between in vitro and in vivo conditions.

CONCLUSIONS

The biological characteristics of noncanonical coronavirus transcripts were experimentally validated for the first time. The identified features of noncanonical transcripts in terms of abundance, reproducibility and variety extend the current model for coronavirus gene expression. The capability of coronaviruses to regulate the species and amounts of noncanonical transcripts may contribute to the pathogenesis of coronaviruses during infection, posing potential challenges in disease control. Thus, the biology of noncanonical transcripts both in vitro and in vivo revealed here can provide a database for biological research, contributing to the development of antiviral strategies.

Unveiling the biology of defective viral genomes in vitro and in vivo: implications for gene expression and pathogenesis of coronavirus

BACKGROUND

Defective viral genome (DVG) is a truncated version of the full-length virus genome identified in most RNA viruses during infection. The synthesis of DVGs in coronavirus has been suggested; however, the fundamental characteristics of coronavirus DVGs in gene expression and pathogenesis have not been systematically analyzed.

METHODS

Nanopore direct RNA sequencing was used to investigate the characteristics of coronavirus DVGs in gene expression including reproducibility, abundance, species and genome structures for bovine coronavirus in cells, and for mouse hepatitis virus (MHV)-A59 (a mouse coronavirus) in cells and in mice. The MHV-A59 full-length genomic cDNAs (~ 31 kilobases) were in vitro constructed to experimentally validate the origin of coronavirus DVG. The synthesis of DVGs was also experimentally identified by RT-PCR followed by sequencing. In addition, the alterations of DVGs in amounts and species under different infection environments and selection pressures including the treatment of antiviral remdesivir and interferon were evaluated based on the banding patterns by RT-PCR.

RESULTS

The results are as follows: (i) the structures of DVGs are with diversity, (ii) DVGs are overall synthesized with moderate (MHV-A59 in cells) to high (BCoV in cells and MHV-A59 in mice) reproducibility under regular infection with the same virus inoculum, (iii) DVGs can be synthesized from the full-length coronavirus genome, (iv) the sequences flanking the recombination point of DVGs are AU-rich and thus may contribute to the recombination events during gene expression, (v) the species and amounts of DVG are altered under different infection environments, and (vi) the biological nature of DVGs between in vitro and in vivo is similar.

CONCLUSIONS

The identified biological characteristics of coronavirus DVGs in terms of abundance, reproducibility, and variety extend the current model for coronavirus gene expression. In addition, the biological features of alterations in amounts and species of coronavirus DVGs under different infection environments may assist the coronavirus to adapt to the altered environments for virus fitness and may contribute to the coronavirus pathogenesis. Consequently, the unveiled biological features may assist the community to study the gene expression mechanisms of DVGs and their roles in pathogenesis, contributing to the development of antiviral strategy and public health.

Towards Understanding Long COVID: SARS-CoV-2 Strikes the Host Cell Nucleus

Despite what its name suggests, the effects of the COVID-19 pandemic causative agent "Severe Acute Respiratory Syndrome Coronavirus-2" (SARS-CoV-2) were not always confined, neither temporarily (being long-term rather than acute, referred to as Long COVID) nor spatially (affecting several body systems). Moreover, the in-depth study of this ss(+) RNA virus is defying the established scheme according to which it just had a lytic cycle taking place confined to cell membranes and the cytoplasm, leaving the nucleus basically "untouched". Cumulative evidence shows that SARS-CoV-2 components disturb the transport of certain proteins through the nuclear pores. Some SARS-CoV-2 structural proteins such as Spike (S) and Nucleocapsid (N), most non-structural proteins (remarkably, Nsp1 and Nsp3), as well as some accessory proteins (ORF3d, ORF6, ORF9a) can reach the nucleoplasm either due to their nuclear localization signals (NLS) or taking a shuttle with other proteins. A percentage of SARS-CoV-2 RNA can also reach the nucleoplasm. Remarkably, controversy has recently been raised by proving that-at least under certain conditions-, SARS-CoV-2 sequences can be retrotranscribed and inserted as DNA in the host genome, giving rise to chimeric genes. In turn, the expression of viral-host chimeric proteins could potentially create neo-antigens, activate autoimmunity and promote a chronic pro-inflammatory state.

Putative Host-Derived Insertions in the Genomes of Circulating SARS-CoV-2 Variants

Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2-infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template-switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution. IMPORTANCE Throughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show that host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic material, potentially enhancing its adaptability, virulence, and spread.

Putative host-derived insertions in the genomes of circulating SARS-CoV-2 variants

Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2 infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution.

IMPORTANCE

Throughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic materials, potentially enhancing its adaptability, virulence, and spread.

Genetic Engineering Systems to Study Human Viral Pathogens from the Coronaviridae Family

The COVID-19 pandemic caused by the previously unknown SARS-CoV-2 Betacoronavirus made it extremely important to develop simple and safe cellular systems which allow manipulation of the viral genome and high-throughput screening of its potential inhibitors. In this review, we made an attempt at summarizing the currently existing data on genetic engineering systems used to study not only SARS-CoV-2, but also other viruses from the Coronaviridae family. In addition, the review covers the basic knowledge about the structure and the life cycle of coronaviruses.

The SARS-CoV-2 subgenome landscape and its novel regulatory features

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

Characterizing Transcriptional Regulatory Sequences in Coronaviruses and Their Role in Recombination

Novel coronaviruses, including SARS-CoV-2, SARS, and MERS, often originate from recombination events. The mechanism of recombination in RNA viruses is template switching. Coronavirus transcription also involves template switching at specific regions, called transcriptional regulatory sequences (TRS). It is hypothesized but not yet verified that TRS sites are prone to recombination events. Here, we developed a tool called SuPER to systematically identify TRS in coronavirus genomes and then investigated whether recombination is more common at TRS. We ran SuPER on 506 coronavirus genomes and identified 465 TRS-L and 3,509 TRS-B. We found that the TRS-L core sequence (CS) and the secondary structure of the leader sequence are generally conserved within coronavirus genera but different between genera. By examining the location of recombination breakpoints with respect to TRS-B CS, we observed that recombination hotspots are more frequently colocated with TRS-B sites than expected.

Effects of Coronavirus Persistence on the Genome Structure and Subsequent Gene Expression, Pathogenicity and Adaptation Capability

Coronaviruses are able to establish persistence. However, how coronaviruses react to persistence and whether the selected viruses have altered their characteristics remain unclear. In this study, we found that the persistent infection of bovine coronavirus (BCoV), which is in the same genus as SARS-COV-2, led to alterations of genome structure, attenuation of gene expression, and the synthesis of subgenomic mRNA (sgmRNA) with a previously unidentified pattern. Subsequent analyses revealed that the altered genome structures were associated with the attenuation of gene expression. In addition, the genome structure at the 5' terminus and the cellular environment during the persistence were responsible for the sgmRNA synthesis, solving the previously unanswered question regarding the selection of transcription regulatory sequence for synthesis of BCoV sgmRNA 12.7. Although the BCoV variants (BCoV-p95) selected under the persistence replicated efficiently in cells without persistent infection, its pathogenicity was still lower than that of wild-type (wt) BCoV. Furthermore, in comparison with wt BCoV, the variant BCoV-p95 was not able to efficiently adapt to the challenges of alternative environments, suggesting wt BCoV is genetically robust. We anticipate that the findings derived from this fundamental research can contribute to the disease control and treatments against coronavirus infection including SARS-CoV-2.

Characterizing transcriptional regulatory sequences in coronaviruses and their role in recombination

Novel coronaviruses, including SARS-CoV-2, SARS, and MERS, often originate from recombination events. The mechanism of recombination in RNA viruses is template switching. Coronavirus transcription also involves template switching at specific regions, called transcriptional regulatory sequences (TRS). It is hypothesized but not yet verified that TRS sites are prone to recombination events. Here, we developed a tool called SuPER to systematically identify TRS in coronavirus genomes and then investigated whether recombination is more common at TRS. We ran SuPER on 506 coronavirus genomes and identified 465 TRS-L and 3509 TRS-B. We found that the TRS-L core sequence (CS) and the secondary structure of the leader sequence are generally conserved within coronavirus genera but different between genera. By examining the location of recombination breakpoints with respect to TRS-B CS, we observed that recombination hotspots are more frequently co-located with TRS-B sites than expected.

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

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

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

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

Characterization of the Role of Hexamer AGUAAA and Poly(A) Tail in Coronavirus Polyadenylation

Similar to eukaryotic mRNA, the positive-strand coronavirus genome of ~30 kilobases is 5’-capped and 3’-polyadenylated. It has been demonstrated that the length of the coronaviral poly(A) tail is not static but regulated during infection; however, little is known regarding the factors involved in coronaviral polyadenylation and its regulation. Here, we show that during infection, the level of coronavirus poly(A) tail lengthening depends on the initial length upon infection and that the minimum length to initiate lengthening may lie between 5 and 9 nucleotides. By mutagenesis analysis, it was found that (i) the hexamer AGUAAA and poly(A) tail are two important elements responsible for synthesis of the coronavirus poly(A) tail and may function in concert to accomplish polyadenylation and (ii) the function of the hexamer AGUAAA in coronaviral polyadenylation is position dependent. Based on these findings, we propose a process for how the coronaviral poly(A) tail is synthesized and undergoes variation. Our results provide the first genetic evidence to gain insight into coronaviral polyadenylation.

Dependence of coronavirus RNA replication on an NH2-terminal partial nonstructural protein 1 in cis

Genomes of positive (+)-strand RNA viruses use cis-acting signals to direct both translation and replication. Here we examine two 5'-proximal cis-replication signals of different character in a defective interfering (DI) RNA of the bovine coronavirus (BCoV) that map within a 322-nucleotide (nt) sequence (136 nt from the genomic 5' untranslated region and 186 nt from the nonstructural protein 1 [nsp1]-coding region) not found in the otherwise-identical nonreplicating subgenomic mRNA7 (sgmRNA7). The natural DI RNA is structurally a fusion of the two ends of the BCoV genome that results in a single open reading frame between a partial nsp1-coding region and the entire N gene. (i) In the first examination, mutation analyses of a recently discovered long-range RNA-RNA base-paired structure between the 5' untranslated region and the partial nsp1-coding region showed that it, possibly in concert with adjacent stem-loops, is a cis-acting replication signal in the (+) strand. We postulate that the higher-order structure promotes (+)-strand synthesis. (ii) In the second examination, analyses of multiple frame shifts, truncations, and point mutations within the partial nsp1-coding region showed that synthesis of a PEFP core amino acid sequence within a group A lineage betacoronavirus-conserved NH2-proximal WAPEFPWM domain is required in cis for DI RNA replication. We postulate that the nascent protein, as part of an RNA-associated translating complex, acts to direct the DI RNA to a critical site, enabling RNA replication. We suggest that these results have implications for viral genome replication and explain, in part, why coronavirus sgmRNAs fail to replicate.

IMPORTANCE

cis-Acting RNA and protein structures that regulate (+)-strand RNA virus genome synthesis are potential sites for blocking virus replication. Here we describe two: a previously suspected 5'-proximal long-range higher-order RNA structure and a novel nascent NH2-terminal protein component of nsp1 that are common among betacoronaviruses of group A lineage.

Identification of cis-acting elements on positive-strand subgenomic mRNA required for the synthesis of negative-strand counterpart in bovine coronavirus

It has been demonstrated that, in addition to genomic RNA, sgmRNA is able to serve as a template for the synthesis of the negative-strand [(−)-strand] complement. However, the cis-acting elements on the positive-strand [(+)-strand] sgmRNA required for (−)-strand sgmRNA synthesis have not yet been systematically identified. In this study, we employed real-time quantitative reverse transcription polymerase chain reaction to analyze the cis-acting elements on bovine coronavirus (BCoV) sgmRNA 7 required for the synthesis of its (−)-strand counterpart by deletion mutagenesis. The major findings are as follows. (1) Deletion of the 5'-terminal leader sequence on sgmRNA 7 decreased the synthesis of the (−)-strand sgmRNA complement. (2) Deletions of the 3' untranslated region (UTR) bulged stem-loop showed no effect on (−)-strand sgmRNA synthesis; however, deletion of the 3' UTR pseudoknot decreased the yield of (−)-strand sgmRNA. (3) Nucleotides positioned from −15 to −34 of the sgmRNA 7 3'-terminal region are required for efficient (−)-strand sgmRNA synthesis. (4) Nucleotide species at the 3'-most position (−1) of sgmRNA 7 is correlated to the efficiency of (−)-strand sgmRNA synthesis. These results together suggest, in principle, that the 5'- and 3'-terminal sequences on sgmRNA 7 harbor cis-acting elements are critical for efficient (−)-strand sgmRNA synthesis in BCoV.

Regulation of coronaviral poly(A) tail length during infection

The positive-strand coronavirus genome of ~30 kilobase in length and subgenomic (sg) mRNAs of shorter lengths, are 5’ and 3’-co-terminal by virtue of a common 5’-capped leader and a common 3’-polyadenylated untranslated region. Here, by ligating head-to-tail viral RNAs from bovine coronavirus-infected cells and sequencing across the ligated junctions, it was learned that at the time of peak viral RNA synthesis [6 hours postinfection (hpi)] the 3’ poly(A) tail on genomic and sgmRNAs is ~65 nucleotides (nt) in length. Surprisingly, this length was found to vary throughout infection from ~45 nt immediately after virus entry (at 0 to 4 hpi) to ~65 nt later on (at 6 h to 9 hpi) and from ~65 nt (at 6 h to 9 hpi) to ~30 nt (at 120-144 hpi). With the same method, poly(U) sequences of the same lengths were simultaneously found on the ligated viral negative-strand RNAs. Functional analyses of poly(A) tail length on specific viral RNA species, furthermore, revealed that translation, in vivo, of RNAs with the longer poly(A) tail was enhanced over those with the shorter poly(A). Although the mechanisms by which the tail lengths vary is unknown, experimental results together suggest that the length of the poly(A) and poly(U) tails is regulated. One potential function of regulated poly(A) tail length might be that for the coronavirus genome a longer poly(A) favors translation. The regulation of coronavirus translation by poly(A) tail length resembles that during embryonal development suggesting there may be mechanistic parallels.

An optimal cis-replication stem-loop IV in the 5' untranslated region of the mouse coronavirus genome extends 16 nucleotides into open reading frame 1

The 288-nucleotide (nt) 3' untranslated region (UTR) in the genome of the bovine coronavirus (BCoV) and 339-nt 3' UTR in the severe acute respiratory syndrome (SARS) coronavirus (SCoV) can each replace the 301-nt 3' UTR in the mouse hepatitis coronavirus (MHV) for virus replication, thus demonstrating common 3' cis-replication signals. Here, we show that replacing the 209-nt MHV 5' UTR with the ∼63%-sequence-identical 210-nt BCoV 5' UTR by reverse genetics does not yield viable virus, suggesting 5' end signals are more stringent or possibly are not strictly 5' UTR confined. To identify potential smaller, 5'-common signals, each of three stem-loop (SL) signaling domains and one inter-stem-loop domain from the BCoV 5' UTR was tested by replacing its counterpart in the MHV genome. The SLI/II domain (nucleotides 1 to 84) and SLIII domain (nucleotides 85 to 141) each immediately enabled near-wild-type (wt) MHV-like progeny, thus behaving similarly to comparable 5'-proximal regions of the SCoV 5' UTR as shown by others. The inter-stem-loop domain (nt 142 to 173 between SLs III and IV) enabled small plaques only after genetic adaptation. The SLIV domain (nt 174 to 210) required a 16-nt extension into BCoV open reading frame 1 (ORF1) for apparent stabilization of a longer BCoV SLIV (nt 174 to 226) and optimal virus replication. Surprisingly, pleiomorphic SLIV structures, including a terminal loop deletion, were found among debilitated progeny from intra-SLIV chimeras. The results show the inter-stem-loop domain to be a potential novel species-specific cis-replication element and that cis-acting SLIV in the viral genome extends into ORF1 in a manner that stabilizes its lower stem and is thus not 5' UTR confined.

Subgenomic messenger RNA amplification in coronaviruses

Coronaviruses possess the largest known RNA genome, a 27- to 32-kb (+)-strand molecule that replicates in the cytoplasm. During virus replication, a 3' coterminal nested set of five to eight subgenomic (sg) mRNAs are made that are also 5' coterminal with the genome, because they carry the genomic leader as the result of discontinuous transcription at intergenic donor signals during (-)-strand synthesis when templates for sgmRNA synthesis are made. An unanswered question is whether the sgmRNAs, which appear rapidly and abundantly, undergo posttranscriptional amplification. Here, using RT-PCR and sequence analyses of head-to-tail-ligated (-) strands, we show that after transfection of an in vitro-generated marked sgmRNA into virus-infected cells, the sgmRNA, like the genome, can function as a template for (-)-strand synthesis. Furthermore, when the transfected sgmRNA contains an internally placed RNA-dependent RNA polymerase template-switching donor signal, discontinuous transcription occurs at this site, and a shorter, 3' terminally nested leader-containing sgmRNA is made, as evidenced by its leader-body junction and by the expression of a GFP gene. Thus, in principle, the longer-nested sgmRNAs in a natural infection, all of which contain potential internal template-switching donor signals, can function to increase the number of the shorter 3'-nested sgmRNAs. One predicted advantage of this behavior for coronavirus survivability is an increased chance of maintaining genome fitness in the 3' one-third of the genome via a homologous recombination between the (now independently abundant) WT sgmRNA and a defective genome.

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

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

Coronavirus Genome Replication

Viruses belonging to the family Coronaviridae are unique among RNA viruses because of the unusually large size of their genome, which is of messenger- or positive- or plus-sense. It is ∼30,000 bases or 2–3 times larger than the genomes of most other RNA viruses. Coronaviruses belong to the order Nidovirales, the other three families being the Arteriviridae, Toroviridae and Roniviridae. (For a review of classification and evolutionary relatedness of Nidovirales see Gorbalenya et al. 2006.) This grouping is based on the arrangement and relatedness of open reading frames within their genomes and on the presence in infected cells of multiple subgenomic mRNAs that form a 3'-co-terminal, nested set with the genome. Among the Nidovirales, coronaviruses (and toroviruses) are unique in their possession of a helical nucleocapsid, which is unusual for plus-stranded but not minus-stranded RNA viruses; plus-stranded RNA-containing plant viruses in the Closteroviridae and in the Tobamovirus genus also possess helical capsids. Coronaviruses are very successful and have infected many species of animals, including bats, birds (poultry) and mammals, such as humans and livestock. Coronavirus species are classified into three groups, which were based originally on cross-reacting antibodies and more recently on nucleotide sequence relatedness (Gonzalez et al. 2003). There have been several reviews of coronaviruses published recently and the reader is referred to them for more extensive references (Enjuanes et al. 2006; Masters 2006; Pasternak et al. 2006; Sawicki and Sawicki 2005; Sawicki et al. 2007; Ziebuhr 2005).

Proteomics analysis unravels the functional repertoire of coronavirus nonstructural protein 3

Severe acute respiratory syndrome (SARS) coronavirus infection and growth are dependent on initiating signaling and enzyme actions upon viral entry into the host cell. Proteins packaged during virus assembly may subsequently form the first line of attack and host manipulation upon infection. A complete characterization of virion components is therefore important to understanding the dynamics of early stages of infection. Mass spectrometry and kinase profiling techniques identified nearly 200 incorporated host and viral proteins. We used published interaction data to identify hubs of connectivity with potential significance for virion formation. Surprisingly, the hub with the most potential connections was not the viral M protein but the nonstructural protein 3 (nsp3), which is one of the novel virion components identified by mass spectrometry. Based on new experimental data and a bioinformatics analysis across the Coronaviridae, we propose a higher-resolution functional domain architecture for nsp3 that determines the interaction capacity of this protein. Using recombinant protein domains expressed in Escherichia coli, we identified two additional RNA-binding domains of nsp3. One of these domains is located within the previously described SARS-unique domain, and there is a nucleic acid chaperone-like domain located immediately downstream of the papain-like proteinase domain. We also identified a novel cysteine-coordinated metal ion-binding domain. Analyses of interdomain interactions and provisional functional annotation of the remaining, so-far-uncharacterized domains are presented. Overall, the ensemble of data surveyed here paint a more complete picture of nsp3 as a conserved component of the viral protein processing machinery, which is intimately associated with viral RNA in its role as a virion component.