Generation and selection of coronavirus defective interfering RNA with large open reading frame by RNA recombination and possible editing

Evolution of naturally arising SARS-CoV-2 defective interfering particles

Defective interfering (DI) particles arise during virus propagation, are conditional on parental virus for replication and packaging, and interfere with viral expansion. There is much interest in developing DIs as anti-viral agents. Here we characterize DI particles that arose following serial passaging of SARS-CoV-2 at high multiplicity of infection. The prominent DIs identified have lost ~84% of the SARS-CoV-2 genome and are capable of attenuating parental viral titers. Synthetic variants of the DI genomes also interfere with infection and can be used as conditional, gene delivery vehicles. In addition, the DI genomes encode an Nsp1-10 fusion protein capable of attenuating viral replication. These results identify naturally selected defective viral genomes that emerged and stably propagated in the presence of parental virus.

Genomes from defective interfering (DI) particles following serial passaging of SARS-CoV-2 reveal a fusion protein that attenuates viral replication. Synthetic, recombinant DI genomes are designed to interfere with SARS-CoV-2 replication.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called ‘Tiled-ClickSeq’, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

A synthetic defective interfering SARS-CoV-2

Viruses thrive by exploiting the cells they infect, but in order to replicate and infect other cells they must produce viral proteins. As a result, viruses are also susceptible to exploitation by defective versions of themselves that do not produce such proteins. A defective viral genome with deletions in protein-coding genes could still replicate in cells coinfected with full-length viruses. Such a defective genome could even replicate faster due to its shorter size, interfering with the replication of the virus. We have created a synthetic defective interfering version of SARS-CoV-2, the virus causing the Covid-19 pandemic, assembling parts of the viral genome that do not code for any functional protein but enable the genome to be replicated and packaged. This synthetic defective genome replicates three times faster than SARS-CoV-2 in coinfected cells, and interferes with it, reducing the viral load of infected cells by half in 24 hours. The synthetic genome is transmitted as efficiently as the full-length genome, suggesting the location of the putative packaging signal of SARS-CoV-2. A version of such a synthetic construct could be used as a self-promoting antiviral therapy: by enabling replication of the synthetic genome, the virus would promote its own demise.

Parallel ClickSeq and Nanopore sequencing elucidates the rapid evolution of defective-interfering RNAs in Flock House virus

Defective-Interfering RNAs (DI-RNAs) have long been known to play an important role in virus replication and transmission. DI-RNAs emerge during virus passaging in both cell-culture and their hosts as a result of non-homologous RNA recombination. However, the principles of DI-RNA emergence and their subsequent evolution have remained elusive. Using a combination of long- and short-read Next-Generation Sequencing, we have characterized the formation of DI-RNAs during serial passaging of Flock House virus (FHV) in cell-culture over a period of 30 days in order to elucidate the pathways and potential mechanisms of DI-RNA emergence and evolution. For short-read RNAseq, we employed 'ClickSeq' due to its ability to sensitively and confidently detect RNA recombination events with nucleotide resolution. In parallel, we used the Oxford Nanopore Technologies's (ONT) MinION to resolve full-length defective and wild-type viral genomes. Together, these accurately resolve both rare and common RNA recombination events, determine the correlation between recombination events, and quantifies the relative abundance of different DI-RNAs throughout passaging. We observe the formation of a diverse pool of defective RNAs at each stage of viral passaging. However, many of these 'intermediate' species, while present in early stages of passaging, do not accumulate. After approximately 9 days of passaging we observe the rapid accumulation of DI-RNAs with a correlated reduction in specific infectivity and with the Nanopore data find that DI-RNAs are characterized by multiple RNA recombination events. This suggests that intermediate DI-RNA species are not competitive and that multiple recombination events interact epistatically to confer 'mature' DI-RNAs with their selective advantage allowing for their rapid accumulation. Alternatively, it is possible that mature DI-RNA species are generated in a single event involving multiple RNA rearrangements. These insights have important consequences for our understanding of the mechanisms, determinants and limitations in the emergence and evolution of DI-RNAs.

Citrus tristeza virus: Evolution of Complex and Varied Genotypic Groups

Amongst the Closteroviridae, Citrus tristeza virus (CTV) is almost unique in possessing a number of distinct and characterized strains, isolates of which produce a wide range of phenotype combinations among its different hosts. There is little understanding to connect genotypes to phenotypes, and to complicate matters more, these genotypes are found throughout the world as members of mixed populations within a single host plant. There is essentially no understanding of how combinations of genotypes affect symptom expression and disease severity. We know little about the evolution of the genotypes that have been characterized to date, little about the biological role of their diversity and particularly, about the effects of recombination. Additionally, genotype grouping has not been standardized. In this study we utilized an extensive array of CTV genomic information to classify the major genotypes, and to determine the major evolutionary processes that led to their formation and subsequent retention. Our analyses suggest that three major processes act on these genotypes: (1) ancestral diversification of the major CTV lineages, followed by (2) conservation and co-evolution of the major functional domains within, though not between CTV genotypes, and (3) extensive recombination between lineages that have given rise to new genotypes that have subsequently been retained within the global population. The effects of genotype diversity and host-interaction are discussed, as is a proposal for standardizing the classification of existing and novel CTV genotypes.

A point mutation within the replicase gene differentially affects coronavirus genome versus minigenome replication

During the construction of the transmissible gastroenteritis virus (TGEV) full-length cDNA clone, a point mutation at position 637 that was present in the defective minigenome DI-C was maintained as a genetic marker. Sequence analysis of the recovered viruses showed a reversion at this position to the original virus sequence. The effect of point mutations at nucleotide 637 was analyzed by reverse genetics using a TGEV full-length cDNA clone and cDNAs from TGEV-derived minigenomes. The replacement of nucleotide 637 of TGEV genome by a T, as in the DI-C sequence, or an A severely affected virus recovery from the cDNA, yielding mutant viruses with low titers and small plaques compared to those of the wild type. In contrast, T or A at position 637 was required for minigenome rescue in trans by the helper virus. No relationship between these observations and RNA secondary-structure predictions was found, indicating that mutations at nucleotide 637 most likely had an effect at the protein level. Nucleotide 637 occupies the second codon position at amino acid 108 of the pp1a polyprotein. This position is predicted to map in the N-terminal polyprotein papain-like proteinase (PLP-1) cleavage site at the p9/p87 junction. Replacement of G-637 by A, which causes a drastic amino acid change (Gly to Asp) at position 108, affected PLP-1-mediated cleavage in vitro. A correlation was found between predicted cleaving and noncleaving mutations and efficient virus rescue from cDNA and minigenome amplification, respectively.

Coronavirus genome structure and replication

In addition to the SARS coronavirus (treated separately elsewhere in this volume), the complete genome sequences of six species in the coronavirus genus of the coronavirus family [avian infectious bronchitis virus-Beaudette strain (IBV-Beaudette), bovine coronavirus-ENT strain (BCoV-ENT), human coronavirus-229E strain (HCoV-229E), murine hepatitis virus-A59 strain (MHV-A59), porcine transmissible gastroenteritis-Purdue 115 strain (TGEV-Purdue 115), and porcine epidemic diarrhea virus-CV777 strain (PEDV-CV777)] have now been reported. Their lengths range from 27,317 nt for HCoV-229E to 31,357 nt for the murine hepatitis virus-A59, establishing the coronavirus genome as the largest known among RNA viruses. The basic organization of the coronavirus genome is shared with other members of the Nidovirus order (the torovirus genus, also in the family Coronaviridae, and members of the family Arteriviridae) in that the nonstructural proteins involved in proteolytic processing, genome replication, and subgenomic mRNA synthesis (transcription) (an estimated 14–16 end products for coronaviruses) are encoded within the 5′-proximal two-thirds of the genome on gene 1 and the (mostly) structural proteins are encoded within the 3′-proximal one-third of the genome (8–9 genes for coronaviruses). Genes for the major structural proteins in all coronaviruses occur in the 5′ to 3′ order as S, E, M, and N. The precise strategy used by coronaviruses for genome replication is not yet known, but many features have been established. This chapter focuses on some of the known features and presents some current questions regarding genome replication strategy, the cis-acting elements necessary for genome replication [as inferred from defective interfering (DI) RNA molecules], the minimum sequence requirements for autonomous replication of an RNA replicon, and the importance of gene order in genome replication.

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.

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

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

In vitro translational analysis of genomic, defective, and satellite RNAs of Cryphonectria hypovirus 3-GH2

Cryphonectria hypovirus 3-GH2 (CHV3-GH2) is a member of the fungal virus family Hypoviridae that differs from previously characterized members in having a single large open reading frame with the potential to encode a protein of 326 kDa from its 9.8-kb genome. The N-terminal portion of the ORF contains sequence motifs that are somewhat similar to papain-like proteinases identified in other hypoviruses. Translation of the ORF is predicted to release autocatalytically a 32.5-kDa protein. A defective RNA, predicted to encode a 91.6-kDa protein representing most of the N-terminal proteinase fused to the entire putative helicase domain, and two satellite RNAs, predicted to encode very small proteins, also are associated with CHV3-GH2 infected fungal cultures. We performed in vitro translation experiments to examine expression of these RNAs. Translation of three RT-PCR clones representing different lengths of the amino-terminal portion of the ORF of the genomic RNA resulted in autocatalytic release of the predicted 32.5-kDa protein. Site-directed mutagenesis was used to map the processing site between Gly(297) and Thr(298). In vitro translation of multiple independent cDNA clones of CHV3-GH2-defective RNA 2 resulted in protein products of approximately 92 kDa, predicted to be the full-length translation product, 32 kDa, predicted to represent the N-terminal proteinase, and 60 kDa, predicted to represent the C-terminal two-thirds of the full-length product. In vitro translation of cDNA clones representing satellite RNA 4 resulted in products of slightly less than 10 kDa, consistent with the predicted 9.4 kDa product.

Sequence elements involved in the rescue of IBV defective RNA CD-91

Deletion mutagenesis has been used to identify essential regions for rescue of coronavirus defective RNAs (D-RNAs). Using this technique on a cloned IBV D-RNA CD-91, we have identified a region potentially important in its rescue. Comparing the sequence of D-RNAs rescued with those not rescued we have deduced that a 72 base region corresponding to base number 13,824 to 13,896 in the viral genome is required for rescue. This may be an IBV D-RNA packaging signal or a cis-acting element involved in replication. Further experiments and modification of our techniques will be required to differentiate between the two processes.

Deletion analysis of a defective interfering Semliki Forest virus RNA genome defines a region in the nsP2 sequence that is required for efficient packaging of the genome into virus particles

The 1,244-nucleotide genome of Semliki Forest virus (SFV) defective interfering (DI) RNA 19 (DI-19) is coterminal with the infectious genome and contains two major deletions. One deletion removes the end of the nsP1 gene and the beginning of the nsP2 gene, and the other removes the end of the nsP2 gene, the nsP3 and nsP4 genes, and all of the structural protein genes (M. Thomson and N. J. Dimmock, Virology 199:354-365, 1994). Like all DI SFV RNAs, DI-19 contains three regions that are conserved. Region a comprises the 5' terminus continuous with part of the nsP1 gene, region b comprises a central part of the nsP2 gene, and region c comprises the 3' terminus and the associated untranslated region. A deletion analysis of the 265-nucleotide b region (nucleotides 679 to 943, inclusive) was undertaken to determine its role in genome replication and packaging into DI virus particles. Deleted plasmids were constructed and transcribed, and the resulting DI RNAs were transfected into SFV-infected BHK cells. Putative progeny DI virus particles that had been released into the tissue culture fluid were then serially passaged in new monolayers together with added high-multiplicity SFV, and cells and tissue culture fluids were tested for the presence of DI RNA by reverse transcription-PCR. DI RNA that had all of the b region deleted was replicated well in BHK-21 cells, as shown by the presence of large amounts of negative-sense DI RNA and an increase in the amount of positive-sense RNA in the cytoplasm, but was packaged very inefficiently, as indicated by very low amounts of DI RNA in the tissue culture fluid. The genome of a deletion mutant that retained the 3' 224 nucleotides of region b was packaged successfully, but one that retained only the 5' 41 nucleotides was not detected in the tissue culture fluid. These and other data suggest that nucleotides 720 to 777 of region b are of particular importance in the packaging process. This finding agrees with data obtained with Ross River virus and contrasts with the well-studied Sindbis alphavirus major packaging signal that is located within the nsP1 gene.

A bulged stem-loop structure in the 3' untranslated region of the genome of the coronavirus mouse hepatitis virus is essential for replication

The 3' untranslated region (UTR) of the positive-sense RNA genome of the coronavirus mouse hepatitis virus (MHV) contains sequences that are necessary for the synthesis of negative-strand viral RNA as well as sequences that may be crucial for both genomic and subgenomic positive-strand RNA synthesis. We have found that the entire 3' UTR of MHV could be replaced by the 3' UTR of bovine coronavirus (BCV), which diverges overall by 31% in nucleotide sequence. This exchange between two viruses that are separated by a species barrier was carried out by targeted RNA recombination. Our results define regions of the two 3' UTRs that are functionally equivalent despite having substantial sequence substitutions, deletions, or insertions with respect to each other. More significantly, our attempts to generate an unallowed substitution of a particular portion of the BCV 3' UTR for the corresponding region of the MHV 3' UTR led to the discovery of a bulged stem-loop RNA secondary structure, adjacent to the stop codon of the nucleocapsid gene, that is essential for MHV viral RNA replication.

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 and packaging of coronavirus infectious bronchitis virus defective RNAs lacking a long open reading frame

The construction of a full-length clone of the avian coronavirus infectious bronchitis virus (IBV) defective RNA (D-RNA), CD-91 (9,080 nucleotides [Z. Penzes et al., Virology 203:286-293]), downstream of the bacteriophage T7 promoter is described. Electroporation of in vitro T7-transcribed CD-91 RNA into IBV helper virus-infected primary chick kidney cells resulted in the production of CD-91 RNA as a replicating D-RNA in subsequent passages. Three CD-91 deletion mutants were constructed--CD-44, CD-58, and CD-61--in which 4,639, 3,236, and 2,953 nucleotides, respectively, were removed from CD-91, resulting in the truncation of the CD-91 long open reading frame (ORF) from 6,465 to 1,311, 1,263, or 2,997 nucleotides in CD-44, CD-58, or CD-61, respectively. Electroporation of in vitro T7-transcribed RNA from the three constructs into IBV helper virus-infected cells resulted in the replication and packaging of CD-58 and CD-61 but not CD-44 RNA. The ORF of CD-61 was further truncated by the insertion of stop codons into the CD-61 sequence by PCR mutagenesis, resulting in constructs CD-61T11 (ORF: nucleotides 996 to 1,058, encoding 20 amino acids), CD-61T22 (ORF: nucleotides 996 to 2,294, encoding 432 amino acids), and CD-61T24 (ORF: nucleotides 996 to 2,450, encoding 484 amino acids), all of which were replicated and packaged to the same levels as observed for either CD-61 or CD-91. Analysis of the D-RNAs showed that the CD-91- or CD-61-specific long ORFs had not been restored. Our data indicate that IBV D-RNAs based on the natural D-RNA, CD-91, do not require a long ORF for efficient replication. In addition, a 1.4-kb sequence, corresponding to IBV sequence at the 5' end of the 1b gene, may be involved in the packaging of IBV D-RNAs or form part of a cis-acting replication element.

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.

cis Requirement for N-specific protein sequence in bovine coronavirus defective interfering RNA replication

A naturally occurring 2.2-kb defective interfering (DI) RNA of the bovine coronavirus, structurally a simple fusion of the genomic termini, contains a single contiguous open reading frame (ORF) or 1.7 kb composed of the 5'-terminal 288 nucleotides of polymerase gene 1a and all 1,344 nucleotides of the nucleocapsid protein (N) gene. The ORF must remain open throughout most of its sequence for replication to occur. To determine the qualitative importance of the N portion of the chimeric ORF in DI RNA replication, transcripts of mutated reporter-containing constructs were tested for replication in helper virus-infected cells. It was determined that the N ORF could not be replaced by the naturally occurring internal I protein ORF, accomplished by deleting the first base in the N start codon which leads to a +1 frameshift, nor could it be replaced by the chloramphenicol acetyltransferase ORF. Furthermore, 3'-terminal truncations of the N gene leaving less than 85% of its total length were likewise not tolerated. Small in-frame deletions and in-frame foreign sequence insertions of up to 99 nucleotides within certain regions of the N ORF were tolerated, however, but the rate of DI RNA accumulation in these cases was lower. These results indicate that there is a requirement for translation of most if not all of the N protein in cis for optimal replication of the bovine coronavirus DI RNA and suggest that a similar requirement may exist for viral genome replication.

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.

Translation but not the encoded sequence is essential for the efficient propagation of the defective interfering RNAs of the coronavirus mouse hepatitis virus

The defective interfering (DI) RNA MIDI of mouse hepatitis virus strain A59 (MHV-A59) contains a large open reading frame (ORF) spanning almost its entire genome. This ORF consists of sequences derived from ORF1a, ORF1b, and the nucleocapsid gene. We have previously demonstrated that mutations that disrupt the ORF decrease the fitness of MIDI and its derivatives (R. J. de Groot, R. G. van der Most, and W. J. M. Spaan, J. Virol. 66:5898-5905, 1992). To determine whether translation of the ORF per se is required or whether the encoded polypeptide or a specific sequence is involved, we analyzed sets of related DI RNAs containing different ORFs. After partial deletion of ORF1b and nucleocapsid gene sequences, disruption of the remaining ORF is still lethal; translation of the entire ORF is not essential, however. When a large fragment of the MHV-A59 spike gene, which is not present in any of the MHV-A59 DI RNAs identified so far, was inserted in-frame into a MIDI derivative, translation across this sequence was vital to DI RNA survival. Thus, the translated sequence is irrelevant, indicating that translation per se plays a crucial role in DI virus propagation. Next, it was examined during which step of the viral life cycle translation plays its role. Since the requirement for translation also exists in DI RNA-transfected and MHV-infected cells, it follows that either the synthesis or degradation of DI RNAs is affected by translation.

Specific binding of host cellular proteins to multiple sites within the 3' end of mouse hepatitis virus genomic RNA

The initial step in mouse hepatitis virus (MHV) RNA replication is the synthesis of negative-strand RNA from a positive-strand genomic RNA template. Our approach to begin studying MHV RNA replication is to identify the cis-acting signals for RNA synthesis and the proteins which recognize these signals at the 3' end of genomic RNA of MHV. To determine whether host cellular and/or viral proteins interact with the 3' end of the coronavirus genome, an RNase T1 protection/gel mobility shift electrophoresis assay was used to examine cytoplasmic extracts from mock- and MHV-JHM-infected 17Cl-1 murine cells for the ability to form complexes with defined regions of the genomic RNA. We demonstrated the specific binding of host cell proteins to multiple sites within the 3' end of MHV-JHM genomic RNA. By using a set of RNA probes with deletions at either the 5' or 3' end or both ends, two distinct binding sites were located. The first protein-binding element was mapped in the 3'-most 42 nucleotides of the genomic RNA [3' (+42) RNA], and the second element was mapped within an 86-nucleotide sequence encompassing nucleotides 171 to 85 from the 3' end of the genome (171-85 RNA). A single potential stem-loop structure is predicted for the 3' (+)42 RNA, and two stem-loop structures are predicted for the 171-85 RNA. Proteins interacting with these two elements were identified by UV-induced covalent cross-linking to labeled RNAs followed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis analysis. The RNA-protein complex formed with the 3'-most 42 nucleotides contains approximately five host polypeptides, a highly labeled protein of 120 kDa and four minor species with sizes of 103, 81, 70, and 55 kDa. The second protein-binding element, contained within a probe representing nucleotides 487 to 85 from the 3' end of the genome, also appears to bind five host polypeptides, 142, 120, 100, 55, and 33 kDa in size, with the 120-kDa protein being the most abundant. The RNA-protein complexes observed with MHV-infected cells in both RNase protection/gel mobility shift and UV cross-linking assays were identical to those observed with uninfected cells. The possible involvement of the interaction of host proteins with the viral genome during MHV replication is discussed.

Evidence for a pseudoknot in the 3' untranslated region of the bovine coronavirus genome

A potential pseudoknot was found in the 3' untranslated region of the bovine coronavirus genome beginning 63 nt downstream from the stop codon of the N gene. Mutation analysis of the pseudoknot in a cloned defective interfering RNA indicated that this structural element is necessary for defective interfering RNA replication.

Structure and encapsidation of transmissible gastroenteritis coronavirus (TGEV) defective interfering genomes

Serial undiluted passages were performed with the PUR46 strain of TGEV in swine testis (ST) cells. Total cellular RNA was analyzed at different passages after orthophosphate metabolic labeling. Three new defective RNA species of 24, 10.5, and 9.5 kb (DI-A, DI-B, and DI-C respectively) were detected at passage 30, which were highly stable and significantly interfered with helper mRNA synthesis in subsequent passages. By Northern hybridization DIs A, B, and C were detected in purified virions at amounts similar to those of helper RNA. Standard and defective TGEV virions could be sorted in sucrose gradients, indicating that defective and full-length genomes are independently packaged. cDNA synthesis of DI-B and DI-C RNAs was performed by the reverse transcription-polymerase chain reaction (RT-PCR) to give four fragments in each case. Cloning and sequencing of the DI-C PCR products showed that the smallest DI particle comprises 9.5 kb and has 4 discontinuous regions of the genome. It contains 2.1 kb from the 5'-end of the genome, about 7 kb from gene 1b, the first 24 nucleotides of the S gene, 12 nucleotides of ORF 7, and the 0.4 kb of the UTR at the 3'-end.

A cis-acting function for the coronavirus leader in defective interfering RNA replication

To test the hypothesis that the 65-nucleotide (nt) leader on subgenomic mRNAs suffices as a 5'-terminal cis-acting signal for RNA replication, a corollary to the notion that coronavirus mRNAs behave as replicons, synthetic RNA transcripts of a cloned, reporter-containing N mRNA (mRNA 7) of the bovine coronavirus with a precise 5' terminus and a 3' poly(A) of 68 nt were tested for replication after being transfected into helper virus-infected cells. No replication was observed, but synthetic transcripts of a cloned reporter-containing defective interfering (DI) RNA differing from the N mRNA construct by 433 nt of continuous 5'-proximal genomic sequence between the leader and the N open reading frame did replicate and become packaged, indicating the insufficiency of the leader alone as a 5' signal for replication of transfected RNA molecules. The leader was shown to be a necessary part of the cis-acting signal for DI RNA replication, however, since removal of terminal bases that destroyed a predicted intraleader stem-loop also destroyed replicating ability. Surprisingly, when the same stem-loop was disrupted by base substitutions, replication appeared only minimally impaired and the leader was found to have rapidly reverted to wild type during DI RNA replication, a phenomenon reminiscent of high-frequency leader switching in the mouse hepatitis coronavirus. These results suggest that once a minimal structural requirement for leader is fulfilled for initiation of DI RNA replication, the wild-type leader is strongly preferred for subsequent replication. They also demonstrate that, in contrast to reported natural mouse hepatitis coronavirus DI RNAs, the DI RNA of the bovine coronavirus does not require sequence elements originating from discontinuous downstream regions within the polymerase gene for replication or for packaging.

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.

Nonhomologous RNA recombination in tombusviruses: generation and evolution of defective interfering RNAs by stepwise deletions

We used a protoplast system to study the mechanisms involved in the generation and evolution of defective interfering (DI) RNAs of tomato bushy stunt tombusvirus (TBSV). Synthetic transcripts corresponding to different naturally occurring TBSV DI RNAs, or to various artificially constructed TBSV defective RNAs, were analyzed. The relative levels of competitiveness of different DI RNAs were determined by coinoculating their corresponding transcripts into protoplasts along with helper genomic RNA transcripts and monitoring the level of DI RNA accumulation. Further studies were performed to assess the contribution of naked DI RNA stability and DI RNA encapsidation efficiency to the observed levels of competitiveness. In addition, the ability of various defective RNAs to evolve to alternative forms was tested by serially passaging protoplast infections initiated with transcripts corresponding to helper genomic RNA and a single type of defective RNA. These studies, and the analysis of the sequences of observed recombinants, indicate that (i) replication competence is a major factor dictating DI RNA competitiveness and is likely a primary determinant in DI RNA evolution, (ii) DI RNAs are capable of evolving to both smaller and larger forms, and the rates at which various transitions occur differ, (iii) DI RNA-DI RNA recombination and/or rearrangement is responsible for the formation of the evolved RNA molecules which were examined, and (iv) sequence complementarities between positive- and negative-sense strands in the regions of the junctions suggest that, in some cases, base pairing between an incomplete replicase-associated nascent strand and acceptor template may mediate selection of recombination sites. On the basis of our data, we propose a stepwise deletion model to describe the temporal order of events leading to the formation of tombusvirus DI RNAs.

Deletion mapping of a mouse hepatitis virus defective interfering RNA reveals the requirement of an internal and discontiguous sequence for replication

All of the defective interfering (DI) RNAs of mouse hepatitis virus (MHV) contain both the 5' and 3' ends of the viral genomic RNA, which presumably include the cis sequences required for RNA replication. To define the replication signal of MHV RNA, we have used a vaccinia virus-T7 polymerase-transcribed MHV DI RNA to study the effects of sequence deletion on DI RNA replication. Following infection of susceptible cells with a recombinant vaccinia virus expressing T7 RNA polymerase, various cDNA clones derived from a DI RNA (DIssF) of the JHM strain of MHV, which is a 3.5-kb naturally occurring DI RNA, behind a T7 promoter were transfected. On superinfection with a helper MHV, the ability of various DI RNAs to replicate was determined. Serial deletions from the middle of the RNA toward both the 5' and 3' ends demonstrated that 859 nucleotides from the 5' end and 436 nucleotides from the 3' end of the MHV RNA genome were necessary for RNA replication. Surprisingly, an additional stretch of 135 nucleotides located at 3.1 to 3.3 kb from the 5' end of the genome was also required. This stretch is discontiguous from the 5'-end cis replication signal and is present in all of the naturally occurring DI RNAs studied so far. The requirement for a long stretch of 5'- and 3'-end sequences predicts that the subgenomic MHV mRNAs cannot replicate. The efficiency of RNA replication varied with different cDNA constructs, suggesting possible interaction between different regions of DI RNA. The identification of MHV RNA replication signals allowed the construction of an MHV DI-based expression vector, which can express foreign genes, such as the chloramphenicol acetyltransferase gene.