RNA recombination in a coronavirus: recombination between viral genomic RNA and transfected RNA fragments

Systematic detection of co-infection and intra-host recombination in more than 2 million global SARS-CoV-2 samples

Systematic monitoring of SARS-CoV-2 co-infections between different lineages and assessing the risk of intra-host recombinant emergence are crucial for forecasting viral evolution. Here we present a comprehensive analysis of more than 2 million SARS-CoV-2 raw read datasets submitted to the European COVID-19 Data Portal to identify co-infections and intra-host recombination. Co-infection was observed in 0.35% of the investigated cases. Two independent procedures were implemented to detect intra-host recombination. We show that sensitivity is predominantly determined by the density of lineage-defining mutations along the genome, thus we used an expanded list of mutually exclusive defining mutations of specific variant combinations to increase statistical power. We call attention to multiple challenges rendering recombinant detection difficult and provide guidelines for the reduction of false positives arising from chimeric sequences produced during PCR amplification. Additionally, we identify three recombination hotspots of Delta - Omicron BA.1 intra-host recombinants.

Self-Replicating RNA Derived from the Genomes of Positive-Strand RNA Viruses

Self-replicating RNA derived from the genomes of positive-strand RNA viruses represents a powerful tool for both molecular studies on virus biology and approaches to novel safe and effective vaccines. The following chapter summarizes the principles how such RNAs can be established and used for design of vaccines. Due to the large variety of strategies needed to circumvent specific pitfalls in the design of such constructs the technical details of the experiments are not described here but can be found in the cited literature.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic analyses suggest recombination was a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing >400,000 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Heterologous RNA Recombination in the Cystoviruses φ6 and φ8: A Mechanism of Viral Variation and Genome Repair

Recombination and mutation of viral genomes represent major mechanisms for viral evolution and, in many cases, moderate pathogenicity. Segmented genome viruses frequently undergo reassortment of the genome via multiple infection of host organisms, with influenza and reoviruses being well-known examples. Specifically, major genomic shifts mediated by reassortment are responsible for radical changes in the influenza antigenic determinants that can result in pandemics requiring rapid preventative responses by vaccine modifications. In contrast, smaller mutational changes brought about by the error-prone viral RNA polymerases that, for the most part, lack a replication base mispairing editing function produce small mutational changes in the RNA genome during replication. Referring again to the influenza example, the accumulated mutations-known as drift-require yearly vaccine updating and rapid worldwide distribution of each new formulation. Coronaviruses with a large positive-sense RNA genome have long been known to undergo intramolecular recombination likely mediated by copy choice of the RNA template by the viral RNA polymerase in addition to the polymerase-based mutations. The current SARS-CoV-2 origin debate underscores the importance of understanding the plasticity of viral genomes, particularly the mechanisms responsible for intramolecular recombination. This review describes the use of the cystovirus bacteriophage as an experimental model for recombination studies in a controlled manner, resulting in the development of a model for intramolecular RNA genome alterations. The review relates the sequence of experimental studies from the laboratory of Leonard Mindich, PhD at the Public Health Research Institute-then in New York City-and covers a period of approximately 12 years. Hence, this is a historical scientific review of research that has the greatest relevance to current studies of emerging RNA virus pathogens.

Detection of SARS-CoV-2 intra-host recombination during superinfection with Alpha and Epsilon variants in New York City

Recombination is an evolutionary process by which many pathogens generate diversity and acquire novel functions. Although a common occurrence during coronavirus replication, detection of recombination is only feasible when genetically distinct viruses contemporaneously infect the same host. Here, we identify an instance of SARS-CoV-2 superinfection, whereby an individual was infected with two distinct viral variants: Alpha (B.1.1.7) and Epsilon (B.1.429). This superinfection was first noted when an Alpha genome sequence failed to exhibit the classic S gene target failure behavior used to track this variant. Full genome sequencing from four independent extracts reveals that Alpha variant alleles comprise around 75% of the genomes, whereas the Epsilon variant alleles comprise around 20% of the sample. Further investigation reveals the presence of numerous recombinant haplotypes spanning the genome, specifically in the spike, nucleocapsid, and ORF 8 coding regions. These findings support the potential for recombination to reshape SARS-CoV-2 genetic diversity.

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.. [PMID: 35043112 PMCID: PMC8764720 DOI: 10.1101/2022.01.04.474799] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

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.

Reduced subgenomic RNA expression is a molecular indicator of asymptomatic SARS-CoV-2 infection

Background

It is estimated that up to 80% of infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are asymptomatic and asymptomatic patients can still effectively transmit the virus and cause disease. While much of the effort has been placed on decoding single nucleotide variation in SARS-CoV-2 genomes, considerably less is known about their transcript variation and any correlation with clinical severity in human hosts, as defined here by the presence or absence of symptoms.

Methods

To assess viral genomic signatures of disease severity, we conducted a systematic characterization of SARS-CoV-2 transcripts and genetic variants in 81 clinical specimens collected from symptomatic and asymptomatic individuals using multi-scale transcriptomic analyses including amplicon-seq, short-read metatranscriptome and long-read Iso-seq.

Results

Here we show a highly coordinated and consistent pattern of sgRNA expression from individuals with robust SARS-CoV-2 symptomatic infection and their expression is significantly repressed in the asymptomatic infections. We also observe widespread inter- and intra-patient variants in viral RNAs, known as quasispecies frequently found in many RNA viruses. We identify unique sets of deletions preferentially found primarily in symptomatic individuals, with many likely to confer changes in SARS-CoV-2 virulence and host responses. Moreover, these frequently occurring structural variants in SARS-CoV-2 genomes serve as a mechanism to further induce SARS-CoV-2 proteome complexity.

Conclusions

Our results indicate that differential sgRNA expression and structural mutational burden are highly correlated with the clinical severity of SARS-CoV-2 infection. Longitudinally monitoring sgRNA expression and structural diversity could further guide treatment responses, testing strategies, and vaccine development.

Wong and Ngan et al. characterize the expression and structural variation of SARS-CoV-2 subgenomic RNAs (sgRNAs) in diagnostic specimens from symptomatic and asymptomatic COVID-19 patients. In a series of genomic and transcriptomic analyses, the authors observe reduced sgRNA expression and a distinct set of structural deletions in asymptomatic infections compared with symptomatic infections.

Infection with SARS-CoV-2 can lead to symptoms of different severity, with some individuals remaining entirely asymptomatic but still responsible for much of the viral transmission. Here, we sought to identify markers of the severity of symptoms of SARS-CoV-2 infection, as defined by the presence or absence of symptoms. We used a combination of methods to study SARS-CoV-2 genes and their readouts, known as transcripts, in clinical samples from people with symptomatic and asymptomatic infections. We demonstrate that transcripts responsible for making viral proteins are seen at lower levels in asymptomatic infections. We also identify structural changes in the viral genes and transcripts that potentially influence the host’s response to the virus. Our study defines potential markers of symptom severity that may ultimately guide risk mitigation and testing strategies.

Update on geographical variation and distribution of SARS-nCoV-2: A systematic review

Knowledge of a new mutant strain of SARS-coronavirus (CoV-2) is enormously essential to identify a targeted drug and for the development of the vaccine. In this article, we systematically reviewed the different mutation strains (variant of concern [VOC] and variant of interest [VOI]) which were found in different countries such as the UK, Singapore, China, Germany, Vietnam, Western Africa, Dublin, Ireland, Brazil, Iran, Italy, France, America, and Philippines. We searched four literature databases (PubMed, EMBASE, NATURE, and Willey online library) with suitable keywords and the time filter was November 2019 to June 16, 2021. To understand the worldwide spread of variants of SARS-CoV-2, we included a total of 27 articles of case reports, clinical and observational studies in the systematic review. However, these variants mostly spread because of their ability to increase transmission, virulence, and escape immunity. So, in this paper is we found mutated strains of SARS-CoV-2 like VOCs that are found in different regions across the globe are ALPHA strain in the U.K, BETA strain in South Africa, GAMMA strain in Brazil, Gamma and Beta strains in European Countries, and some VOIs like Theta variant in the Philippines.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium identifies target cells, alterations in gene expression, and cell state changes

There are currently limited Food and Drug Administration (FDA)-approved drugs and vaccines for the treatment or prevention of Coronavirus Disease 2019 (COVID-19). Enhanced understanding of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and pathogenesis is critical for the development of therapeutics. To provide insight into viral replication, cell tropism, and host-viral interactions of SARS-CoV-2, we performed single-cell (sc) RNA sequencing (RNA-seq) of experimentally infected human bronchial epithelial cells (HBECs) in air-liquid interface (ALI) cultures over a time course. This revealed novel polyadenylated viral transcripts and highlighted ciliated cells as a major target at the onset of infection, which we confirmed by electron and immunofluorescence microscopy. Over the course of infection, the cell tropism of SARS-CoV-2 expands to other epithelial cell types including basal and club cells. Infection induces cell-intrinsic expression of type I and type III interferons (IFNs) and interleukin (IL)-6 but not IL-1. This results in expression of interferon-stimulated genes (ISGs) in both infected and bystander cells. This provides a detailed characterization of genes, cell types, and cell state changes associated with SARS-CoV-2 infection in the human airway.

The zoonotic potential of bat-borne coronaviruses

Seven zoonoses - human infections of animal origin - have emerged from the Coronaviridae family in the past century, including three viruses responsible for significant human mortality (SARS-CoV, MERS-CoV, and SARS-CoV-2) in the past twenty years alone. These three viruses, in addition to two older CoV zoonoses (HCoV-229E and HCoV-NL63) are believed to be originally derived from wild bat reservoir species. We review the molecular biology of the bat-derived Alpha- and Betacoronavirus genera, highlighting features that contribute to their potential for cross-species emergence, including the use of well-conserved mammalian host cell machinery for cell entry and a unique capacity for adaptation to novel host environments after host switching. The adaptive capacity of coronaviruses largely results from their large genomes, which reduce the risk of deleterious mutational errors and facilitate range-expanding recombination events by offering heightened redundancy in essential genetic material. Large CoV genomes are made possible by the unique proofreading capacity encoded for their RNA-dependent polymerase. We find that bat-borne SARS-related coronaviruses in the subgenus Sarbecovirus, the source clade for SARS-CoV and SARS-CoV-2, present a particularly poignant pandemic threat, due to the extraordinary viral genetic diversity represented among several sympatric species of their horseshoe bat hosts. To date, Sarbecovirus surveillance has been almost entirely restricted to China. More vigorous field research efforts tracking the circulation of Sarbecoviruses specifically and Betacoronaviruses more generally is needed across a broader global range if we are to avoid future repeats of the COVID-19 pandemic.

Translational control of coronaviruses

Coronaviruses represent a large family of enveloped RNA viruses that infect a large spectrum of animals. In humans, the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic and is genetically related to SARS-CoV and Middle East respiratory syndrome-related coronavirus (MERS-CoV), which caused outbreaks in 2002 and 2012, respectively. All viruses described to date entirely rely on the protein synthesis machinery of the host cells to produce proteins required for their replication and spread. As such, virus often need to control the cellular translational apparatus to avoid the first line of the cellular defense intended to limit the viral propagation. Thus, coronaviruses have developed remarkable strategies to hijack the host translational machinery in order to favor viral protein production. In this review, we will describe some of these strategies and will highlight the role of viral proteins and RNAs in this process.

Coronavirus Detection in the Clinical Microbiology Laboratory: Are We Ready for Identifying and Diagnosing a Novel Virus?

Endemic species of coronavirus (HCoV-OC43, HCoV-229E, HCoV-NL63, and HCoV-HKU1) are frequent causes of upper respiratory tract infections. Three highly pathogenic coronaviruses have been associated with outbreaks and epidemics and have challenged clinical microbiology laboratories to quickly develop assays for diagnosis. Their initial characterization was achieved by molecular methods. With the great advance in metagenomic whole-genome sequencing directly from clinical specimens, diagnosis of novel coronaviruses could be quickly implemented into the workflow of managing cases of pneumonia of unknown cause, which will markedly affect the time of the initial characterization and accelerate the initiation of outbreak control measures.

Critical Roles of N6-Methyladenosine (m6A) in Cancer and Virus Infection

Studies have shown that epigenetic abnormalities are involved in various diseases, including cancer. In particular, in order to realize precision medicine, the integrated analysis of genetics and epigenetics is considered to be important; detailed epigenetic analysis in the medical field has been becoming increasingly important. In the epigenetics analysis, DNA methylation and histone modification analyses have been actively studied for a long time, and many important findings were accumulated. On the other hand, recently, attention has also been focused on RNA modification in the field of epigenetics; now it is known that RNA modification is associated with various biological functions, such as regulation of gene expression. Among RNA modifications, functional analysis of N⁶-methyladenosine (m⁶A), the most abundant RNA modification found from humans to plants is actively progressing, and it has also been known that m⁶A abnormality is involved in cancer and other diseases. Importantly, recent studies have shown that m⁶A is related to viral infections. Considering the current world situation under threat of viral infections, it is important to deepen knowledge of RNA modification from the viewpoint of viral diseases. Hence, in this review, we have summarized the recent findings regarding the roles of RNA modifications in biological functions, cancer biology, and virus infection, particularly focusing on m⁶A in mRNA.

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium

SARS-CoV-2, the causative agent of COVID-19, has tragically burdened individuals and institutions around the world. There are currently no approved drugs or vaccines for the treatment or prevention of COVID-19. Enhanced understanding of SARS-CoV-2 infection and pathogenesis is critical for the development of therapeutics. To reveal insight into viral replication, cell tropism, and host-viral interactions of SARS-CoV-2 we performed single-cell RNA sequencing of experimentally infected human bronchial epithelial cells (HBECs) in air-liquid interface cultures over a time-course. This revealed novel polyadenylated viral transcripts and highlighted ciliated cells as a major target of infection, which we confirmed by electron microscopy. Over the course of infection, cell tropism of SARS-CoV-2 expands to other epithelial cell types including basal and club cells. Infection induces cell-intrinsic expression of type I and type III IFNs and IL6 but not IL1. This results in expression of interferon-stimulated genes in both infected and bystander cells. We observe similar gene expression changes from a COVID-19 patient ex vivo. In addition, we developed a new computational method termed CONditional DENSity Embedding (CONDENSE) to characterize and compare temporal gene dynamics in response to infection, which revealed genes relating to endothelin, angio-genesis, interferon, and inflammation-causing signaling pathways. In this study, we conducted an in-depth analysis of SARS-CoV-2 infection in HBECs and a COVID-19 patient and revealed genes, cell types, and cell state changes associated with infection.

Emerging coronavirus diseases and future perspectives

Coronavirus related infectious diseases seems to be biggest challenge of 21 century that have been constantly emerging and threating public health around the globe. Coronavirus disease-19 (COVID-19) that was detected as cause of respiratory tract infection in China by end the December 2019 impelled World Health Organization to declare in January 2020 public health emergency of international concern and consequently pandemic in March 2020. Over a past six months COVID-19 pandemic has wrapped up all continents except Antarctica. Scientists around the globe are finding way to tackle and reduce the ultimate risk and size of pandemic with lower morbidity and mortality rates. In this context, technologies such as sequencing, Crispr and artificial intelligence are playing vital role in diagnosis and management of infectious disease in contrast to conventional methods. Despite of this, there is a need to have rapid and early diagnostic tools and systems that recognize infectious disease in asymptotic condition. Here we provide an overview on the recent CoV outbreak and contribution of technologies with the emphasis on the future management for detection of such infectious diseases.

The Architecture of SARS-CoV-2 Transcriptome

SARS-CoV-2 is a betacoronavirus responsible for the COVID-19 pandemic. Although the SARS-CoV-2 genome was reported recently, its transcriptomic architecture is unknown. Utilizing two complementary sequencing techniques, we present a high-resolution map of the SARS-CoV-2 transcriptome and epitranscriptome. DNA nanoball sequencing shows that the transcriptome is highly complex owing to numerous discontinuous transcription events. In addition to the canonical genomic and 9 subgenomic RNAs, SARS-CoV-2 produces transcripts encoding unknown ORFs with fusion, deletion, and/or frameshift. Using nanopore direct RNA sequencing, we further find at least 41 RNA modification sites on viral transcripts, with the most frequent motif, AAGAA. Modified RNAs have shorter poly(A) tails than unmodified RNAs, suggesting a link between the modification and the 3' tail. Functional investigation of the unknown transcripts and RNA modifications discovered in this study will open new directions to our understanding of the life cycle and pathogenicity of SARS-CoV-2.

Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses

Over the past 20 years, several coronaviruses have crossed the species barrier into humans, causing outbreaks of severe, and often fatal, respiratory illness. Since SARS-CoV was first identified in animal markets, global viromics projects have discovered thousands of coronavirus sequences in diverse animals and geographic regions. Unfortunately, there are few tools available to functionally test these viruses for their ability to infect humans, which has severely hampered efforts to predict the next zoonotic viral outbreak. Here, we developed an approach to rapidly screen lineage B betacoronaviruses, such as SARS-CoV and the recent SARS-CoV-2, for receptor usage and their ability to infect cell types from different species. We show that host protease processing during viral entry is a significant barrier for several lineage B viruses and that bypassing this barrier allows several lineage B viruses to enter human cells through an unknown receptor. We also demonstrate how different lineage B viruses can recombine to gain entry into human cells, and confirm that human ACE2 is the receptor for the recently emerging SARS-CoV-2.

This study describes the development of an approach to rapidly screen lineage B betacoronaviruses, such as SARS-CoV and the recently emerged SARS-CoV-2, for receptor usage and their ability to infect cell types from different species. Using it, they confirm human ACE2 as the receptor for SARs-CoV-2 and show that host protease processing during viral entry is a significant barrier for viral entry.

Functional assessment of cell entry and receptor usage for lineage B β-coronaviruses, including 2019-nCoV

Over the past 20 years, several coronaviruses have crossed the species barrier into humans, causing outbreaks of severe, and often fatal, respiratory illness. Since SARS-CoV was first identified in animal markets, global viromics projects have discovered thousands of coronavirus sequences in diverse animals and geographic regions. Unfortunately, there are few tools available to functionally test these novel viruses for their ability to infect humans, which has severely hampered efforts to predict the next zoonotic viral outbreak. Here we developed an approach to rapidly screen lineage B betacoronaviruses, such as SARS-CoV and the recent 2019-nCoV, for receptor usage and their ability to infect cell types from different species. We show that host protease processing during viral entry is a significant barrier for several lineage B viruses and that bypassing this barrier allows several lineage B viruses to enter human cells through an unknown receptor. We also demonstrate how different lineage B viruses can recombine to gain entry into human cells and confirm that human ACE2 is the receptor for the recently emerging 2019-nCoV.

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

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

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

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

Biological species in the viral world

Due to their dependence on cellular organisms for metabolism and replication, viruses are typically named and assigned to species according to their genome structure and the original host that they infect. But because viruses often infect multiple hosts and the numbers of distinct lineages within a host can be vast, their delineation into species is often dictated by arbitrary sequence thresholds, which are highly inconsistent across lineages. Here we apply an approach to determine the boundaries of viral species based on the detection of gene flow within populations, thereby defining viral species according to the biological species concept (BSC). Despite the potential for gene transfer between highly divergent genomes, viruses, like the cellular organisms they infect, assort into reproductively isolated groups and can be organized into biological species. This approach revealed that BSC-defined viral species are often congruent with the taxonomic partitioning based on shared gene contents and host tropism, and that bacteriophages can similarly be classified in biological species. These results open the possibility to use a single, universal definition of species that is applicable across cellular and acellular lifeforms.

Deciphering the biology of porcine epidemic diarrhea virus in the era of reverse genetics

Emergence of the porcine epidemic diarrhea virus (PEDV) as a global threat to the swine industry underlies the urgent need for deeper understanding of this virus. To date, we have yet to identify functions for all the major gene products, much less grasp their implications for the viral life cycle and pathogenic mechanisms. A major reason is the lack of genetic tools for studying PEDV. In this review, we discuss the reverse genetics approaches that have been successfully used to engineer infectious clones of PEDV as well as other potential and complementary methods that have yet to be applied to PEDV. The importance of proper cell culture for successful PEDV propagation and maintenance of disease phenotype are addressed in our survey of permissive cell lines. We also highlight areas of particular relevance to PEDV pathogenesis and disease that have benefited from reverse genetics studies and pressing questions that await resolution by such studies. In particular, we examine the spike protein as a determinant of viral tropism, entry and virulence, ORF3 and its association with cell culture adaptation, and the nucleocapsid protein and its potential role in modulating PEDV pathogenicity. Finally, we conclude with an exploration of how reverse genetics can help mitigate the global impact of PEDV by addressing the challenges of vaccine development.

Self-Replicating RNA

Self-replicating RNA derived from the genomes of positive strand RNA viruses represents a powerful tool for both molecular studies on virus biology and approaches to novel safe and effective vaccines. The following chapter summarizes the principles how such RNAs can be established and used for design of vaccines. Due to the large variety of strategies needed to circumvent specific pitfalls in the design of such constructs the technical details of the experiments are not described here but can be found in the cited literature.

Identification of the Mechanisms Causing Reversion to Virulence in an Attenuated SARS-CoV for the Design of a Genetically Stable Vaccine

A SARS-CoV lacking the full-length E gene (SARS-CoV-∆E) was attenuated and an effective vaccine. Here, we show that this mutant virus regained fitness after serial passages in cell culture or in vivo, resulting in the partial duplication of the membrane gene or in the insertion of a new sequence in gene 8a, respectively. The chimeric proteins generated in cell culture increased virus fitness in vitro but remained attenuated in mice. In contrast, during SARS-CoV-∆E passage in mice, the virus incorporated a mutated variant of 8a protein, resulting in reversion to a virulent phenotype. When the full-length E protein was deleted or its PDZ-binding motif (PBM) was mutated, the revertant viruses either incorporated a novel chimeric protein with a PBM or restored the sequence of the PBM on the E protein, respectively. Similarly, after passage in mice, SARS-CoV-∆E protein 8a mutated, to now encode a PBM, and also regained virulence. These data indicated that the virus requires a PBM on a transmembrane protein to compensate for removal of this motif from the E protein. To increase the genetic stability of the vaccine candidate, we introduced small attenuating deletions in E gene that did not affect the endogenous PBM, preventing the incorporation of novel chimeric proteins in the virus genome. In addition, to increase vaccine biosafety, we introduced additional attenuating mutations into the nsp1 protein. Deletions in the carboxy-terminal region of nsp1 protein led to higher host interferon responses and virus attenuation. Recombinant viruses including attenuating mutations in E and nsp1 genes maintained their attenuation after passage in vitro and in vivo. Further, these viruses fully protected mice against challenge with the lethal parental virus, and are therefore safe and stable vaccine candidates for protection against SARS-CoV.

Zoonotic coronaviruses, including SARS-CoV, Middle East respiratory syndrome (MERS-CoV), porcine epidemic diarrhea virus (PEDV) and swine delta coronavirus (SDCoV) have recently emerged causing high morbidity and mortality in human or piglets. No fully protective therapy is still available for these CoVs. Therefore, the development of efficient vaccines is a high priority. Live attenuated vaccines are considered most effective compared to other types of vaccines, as they induce a long-lived, balanced immune response. However, safety is the main concern of this type of vaccines because attenuated viruses can eventually revert to a virulent phenotype. Therefore, an essential feature of any live attenuated vaccine candidate is its stability. In addition, introduction of several safety guards is advisable to increase vaccine safety. In this manuscript, we analyzed the mechanisms by which an attenuated SARS-CoV reverted to a virulent phenotype and describe the introduction of attenuating deletions that maintained virus stability. The virus, engineered with two safety guards, provided full protection against challenge with a lethal SARS-CoV. Understanding the molecular mechanisms leading to pathogenicity and the in vivo evaluation of vaccine genetic stability contributed to a rational design of a promising SARS-CoV vaccine.

PRRSV structure, replication and recombination: Origin of phenotype and genotype diversity

Porcine reproductive and respiratory disease virus (PRRSV) has the intrinsic ability to adapt and evolve. After 25 years of study, this persistent pathogen has continued to frustrate efforts to eliminate infection of herds through vaccination or other elimination strategies. The purpose of this review is to summarize the research on the virion structure, replication and recombination properties of PRRSV that have led to the extraordinary phenotype and genotype diversity that exists worldwide.

•

Review of structure, replication and recombination of porcine reproductive and respiratory syndrome virus.

•

Homologous recombination to produce conventional subgenomic messenger RNA as well as heteroclite RNA.

•

Discussion of structure, replication and recombination mechanisms that have yielded genotypic and phenotypic diversity.

Comparative analysis of four Massachusetts type infectious bronchitis coronavirus genomes reveals a novel Massachusetts type strain and evidence of natural recombination in the genome

Four Massachusetts-type (Mass-type) strains of infectious bronchitis coronavirus (IBV) were compared genetically with the pathogenic M41 and H120 vaccine strains using the complete genomic sequences. The results revealed that strains ck/CH/LNM/091017 and ck/CH/LDL/101212 were closely related to the H120 vaccine, which suggests that they might represent re-isolations of vaccine strains or variants of vaccine strains that have resulted from the accumulated point mutations after several passages in chickens. In contrast, strains ck/CH/LHLJ/07VII and ck/CH/LHLJ/100902 had a close genetic relationship with the pathogenic M41 strain. In addition, molecular markers have been identified that distinguish between field and vaccine (or vaccine-like) Mass-type viruses, which may be able to differentiate between field and vaccine strains for diagnostic purposes. Phylogenetic analysis, and pairwise comparison of full-length genomes and the nine genes, identified the occurrence of recombination events in the genome of strain CK/VH/LHLJ/07VII, which suggests that this virus originated from recombination events between M41- and H120-like strains at the switch site located at the 3' end of the nucleocapsid (N) genes. To our knowledge, this is the first time that evidence for the evolution and natural recombination under field conditions between Mass-type pathogenic and vaccinal IBV strains has been documented. These findings provide insights into the emergence and evolution of the Mass-type IB coronaviruses and may help to explain the emergence of Mass-type IBV in chicken flocks all over the world.

Evolution of infectious bronchitis virus in Taiwan: characterisation of RNA recombination in the nucleocapsid gene

Avian infectious bronchitis virus (IBV) belongs to the Coronaviridae family and causes significant economic loss in Taiwan (TW), even in flocks that have been extensively immunised with Massachusetts (Mass)-serotype vaccines. Phylogenetic analysis of all non-structural and most structural genes shows that TW IBV is genetically distinct from the US strain and more similar to Chinese (CH) IBV. In contrast, the nucleocapsid (N) gene of TW IBV presents phylogenetic incongruence. RNA recombination at the 5′ end of the N gene between TW and US IBV is shown to be responsible for this discordance. Surprisingly, the recombinant N gene is found in all of tested TW IBV isolates, suggesting that a recombination event gave origin to a founder lineage. Our data indicate that RNA recombination in the recombinant 5′ end of the N gene may have caused the emergence of the current IBV population in Taiwan.

A trans-complementing recombination trap demonstrates a low propensity of flaviviruses for intermolecular recombination

Intermolecular recombination between the genomes of closely related RNA viruses can result in the emergence of novel strains with altered pathogenic potential and antigenicity. Although recombination between flavivirus genomes has never been demonstrated experimentally, the potential risk of generating undesirable recombinants has nevertheless been a matter of concern and controversy with respect to the development of live flavivirus vaccines. As an experimental system for investigating the ability of flavivirus genomes to recombine, we developed a "recombination trap," which was designed to allow the products of rare recombination events to be selected and amplified. To do this, we established reciprocal packaging systems consisting of pairs of self-replicating subgenomic RNAs (replicons) derived from tick-borne encephalitis virus (TBEV), West Nile virus (WNV), and Japanese encephalitis virus (JEV) that could complement each other in trans and thus be propagated together in cell culture over multiple passages. Any infectious viruses with intact, full-length genomes that were generated by recombination of the two replicons would be selected and enriched by end point dilution passage, as was demonstrated in a spiking experiment in which a small amount of wild-type virus was mixed with the packaged replicons. Using the recombination trap and the JEV system, we detected two aberrant recombination events, both of which yielded unnatural genomes containing duplications. Infectious clones of both of these genomes yielded viruses with impaired growth properties. Despite the fact that the replicon pairs shared approximately 600 nucleotides of identical sequence where a precise homologous crossover event would have yielded a wild-type genome, this was not observed in any of these systems, and the TBEV and WNV systems did not yield any viable recombinant genomes at all. Our results show that intergenomic recombination can occur in the structural region of flaviviruses but that its frequency appears to be very low and that therefore it probably does not represent a major risk in the use of live, attenuated flavivirus vaccines.

The molecular biology of coronaviruses

Coronaviruses are large, enveloped RNA viruses of both medical and veterinary importance. Interest in this viral family has intensified in the past few years as a result of the identification of a newly emerged coronavirus as the causative agent of severe acute respiratory syndrome (SARS). At the molecular level, coronaviruses employ a variety of unusual strategies to accomplish a complex program of gene expression. Coronavirus replication entails ribosome frameshifting during genome translation, the synthesis of both genomic and multiple subgenomic RNA species, and the assembly of progeny virions by a pathway that is unique among enveloped RNA viruses. Progress in the investigation of these processes has been enhanced by the development of reverse genetic systems, an advance that was heretofore obstructed by the enormous size of the coronavirus genome. This review summarizes both classical and contemporary discoveries in the study of the molecular biology of these infectious agents, with particular emphasis on the nature and recognition of viral receptors, viral RNA synthesis, and the molecular interactions governing virion assembly.

Nidovirus transcription: how to make sense...?

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

Bovine coronavirus 5'-proximal genomic acceptor hotspot for discontinuous transcription is 65 nucleotides wide

Coronaviruses are positive-strand, RNA-dependent RNA polymerase-utilizing viruses that require a polymerase template switch, characterized as discontinuous transcription, to place a 5'-terminal genomic leader onto subgenomic mRNAs (sgmRNAs). The usually precise switch is thought to occur during the synthesis of negative-strand templates for sgmRNA production and to be directed by heptameric core donor sequences within the genome that match an acceptor core (UCUAAAC in the case of bovine coronavirus) near the 3' end of the 5'-terminal genomic leader. Here it is shown that a 22-nucleotide (nt) donor sequence engineered into a packageable bovine coronavirus defective interfering (DI) RNA and made to match a sequence within the 65-nt virus genomic leader caused a template switch yielding an sgmRNA with only a 33-nt minileader. By changing the donor sequence, acceptor sites between genomic nt 33 and 97 (identical between the DI RNA and the viral genome) could be used to generate sgmRNAs detectable by Northern analysis (approximately 2 to 32 molecules per cell) by 24 h postinfection. Whether the switch was intramolecular only was not determined since a potentially distinguishing acceptor region in the DI RNA rapidly conformed to that in the helper virus genome through a previously described template switch known as leader switching. These results show that crossover acceptor sites for discontinuous transcription (i) need not include the UCUAAAC core and (ii) rest within a surprisingly wide 5'-proximal "hotspot." Overlap of this hotspot with that for leader switching and with elements required for RNA replication suggests that it is part of a larger 5'-proximal multifunctional structure.

Reconstitution in cultured cells of replicating HDV RNA from pairs of less than full-length RNAs

The genome of hepatitis delta virus (HDV) is a small single-stranded circular RNA that is replicated via RNA-directed RNA synthesis. This makes use of a host RNA polymerase, probably pol II, that normally transcribes DNA templates. In vivo, the host polymerase can initiate replication from transfected linear RNAs using intramolecular template-switching. The present studies report that the polymerase could also achieve intermolecular switching leading to "reconstitution" of full-length HDV RNAs following transfection with two linear RNAs that were less than full length and yet lacking different regions of the genome. These two RNAs were synthesized in vitro, gel purified, pre-annealed, and then transfected into delta293, a cell line conditionally expressing the small delta antigen that is essential for HDV replication. Northern analyses of total RNA harvested from transfected cells detected the accumulation of full-length HDV genomic and antigenomic RNAs. Such reconstitution of full-length replicating HDV RNA was also achieved using nine other pairs of antigenomic RNAs and three pairs of genomic RNAs. Annealing of the RNAs prior to transfection was required for detectable HDV reconstitution. A second cell line, Huh7, also supported reconstitution when a pair of RNAs was cotransfected together with mRNA for the small delta protein. Taken together, these results support a model that observed genome reconstitution is a special form of recombination involving intermolecular template switches and they provide insights into the mechanism of RNA-directed RNA transcription catalyzed by a host RNA polymerase.

Mechanism of RNA recombination in carmo- and tombusviruses: evidence for template switching by the RNA-dependent RNA polymerase in vitro

RNA recombination occurs frequently during replication of tombusviruses and carmoviruses, which are related small plus-sense RNA viruses of plants. The most common recombinants generated by these viruses are either defective interfering (DI) RNAs or chimeric satellite RNAs, which are thought to be generated by template switching of the viral RNA-dependent RNA polymerase (RdRp) during the viral replication process. To test if RNA recombination is mediated by the viral RdRp, we used either a purified recombinant RdRp of Turnip crinkle carmovirus or a partially purified RdRp preparation of Cucumber necrosis tombusvirus. We demonstrated that these RdRp preparations generated RNA recombinants in vitro. The RdRp-driven template switching events occurred between either identical templates or two different RNA templates. The template containing a replication enhancer recombined more efficiently than templates containing artificial sequences. We also observed that AU-rich sequences promote recombination more efficiently than GC-rich sequences. Cloning and sequencing of the generated recombinants revealed that the junction sites were located frequently at the ends of the templates (end-to-end template switching). We also found several recombinants that were generated by template switching involving internal positions in the RNA templates. In contrast, RNA ligation-based RNA recombination was not detected in vitro. Demonstration of the ability of carmo- and tombusvirus RdRps to switch RNA templates in vitro supports the copy-choice models of RNA recombination and DI RNA formation for these viruses.

Phylogenetic analysis reveals a low rate of homologous recombination in negative-sense RNA viruses

Recombination is increasingly seen as an important means of shaping genetic diversity in RNA viruses. However, observed recombination frequencies vary widely among those viruses studied to date, with only sporadic occurrences reported in RNA viruses with negative-sense genomes. To determine the extent of homologous recombination in negative-sense RNA viruses, phylogenetic analyses of 79 gene sequence alignments from 35 negative-sense RNA viruses (a total of 2154 sequences) were carried out. Powerful evidence was found for recombination, in the form of incongruent phylogenetic trees between different gene regions, in only five sequences from Hantaan virus, Mumps virus and Newcastle disease virus. This is the first report of recombination in these viruses. More tentative evidence for recombination, where conflicting phylogenetic trees were observed (but were without strong bootstrap support) and/or where putative recombinant regions were very short, was found in three alignments from La Crosse virus and Puumala virus. Finally, patterns of sequence variation compatible with the action of recombination, but not definitive evidence for this process, were observed in a further ten viruses: Canine distemper virus, Crimean-Congo haemorrhagic fever virus, Influenza A virus, Influenza B virus, Influenza C virus, Lassa virus, Pirital virus, Rabies virus, Rift Valley Fever virus and Vesicular stomatitis virus. The possibility of recombination in these viruses should be investigated further. Overall, this study reveals that rates of homologous recombination in negative-sense RNA viruses are very much lower than those of mutation, with many viruses seemingly clonal on current data. Consequently, recombination rate is unlikely to be a trait that is set by natural selection to create advantageous or purge deleterious mutations.

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.

Utilizing fowlpox virus recombinants to generate defective RNAs of the coronavirus infectious bronchitis virus

Coronavirus defective RNAs (D-RNAs) have been used as RNA vectors for the expression of heterologous genes and as vehicles for reverse genetics by modifying coronavirus genomes by targetted recombination. D-RNAs based on the avian coronavirus infectious bronchitis virus (IBV) D-RNA CD-61 have been rescued (replicated and packaged into virions) in a helper virus-dependent manner following electroporation of in vitro-generated T7 transcripts into IBV-infected cells. In order to increase the efficiency of rescue of IBV D-RNAs, cDNAs based on CD-61, under the control of a T7 promoter, were integrated into the fowlpox virus (FPV) genome. The 3'-UTR of the D-RNAs was flanked by a hepatitis delta antigenomic ribozyme and T7 terminator sequence to generate suitable 3' ends for rescue by helper IBV. Cells were co-infected simultaneously with IBV, the recombinant FPV (rFPV) containing the D-RNA sequence and a second rFPV expressing T7 RNA polymerase for the initial expression of the D-RNA transcript, subsequently rescued by helper IBV. Rescue of rFPV-derived CD-61 occurred earlier and with higher efficiency than demonstrated previously for electroporation of in vitro T7-generated RNA transcripts in avian cells. Rescue of CD-61 was also demonstrated for the first time in mammalian cells. The rescue of rFPV-derived CD-61 by M41 helper IBV resulted in leader switching, in which the Beaudette-type leader sequence on CD-61 was replaced with the M41 leader sequence, confirming that helper IBV virus replicated the rFPV-derived D-RNA. An rFPV-derived D-RNA containing the luciferase gene under the control of an IBV transcription-associated sequence was also rescued and expressed luciferase on serial passage.

Feline coronavirus type II strains 79-1683 and 79-1146 originate from a double recombination between feline coronavirus type I and canine coronavirus

Recent evidence suggests that the type II feline coronavirus (FCoV) strains 79-1146 and 79-1683 have arisen from a homologous RNA recombination event between FCoV type I and canine coronavirus (CCV). In both cases, the template switch apparently took place between the S and M genes, giving rise to recombinant viruses which encode a CCV-like S protein and the M, N, 7a, and 7b proteins of FCoV type I (K. Motowaka, T. Hohdatsu, H. Hashimoto, and H. Koyama, Microbiol. Immunol. 40:425-433, 1996; H. Vennema, A. Poland, K. Floyd Hawkins, and N. C. Pedersen, Feline Pract. 23:40-44, 1995). In the present study, we have looked for additional FCoV-CCV recombination sites. Four regions in the pol gene were selected for comparative sequence analysis of the type II FCoV strains 79-1683 and 79-1146, the type I FCoV strains TN406 and UCD1, the CCV strain K378, and the TGEV strain Purdue. Our data show that the type II FCoVs have arisen from double recombination events: additional crossover sites were mapped in the ORF1ab frameshifting region of strain 79-1683 and in the 5' half of ORF1b of strain 79-1146.

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.

Generation of coronavirus spike deletion variants by high-frequency recombination at regions of predicted RNA secondary structure

Coronavirus RNA evolves in the central nervous systems (CNS) of mice during persistent infection. This evolution can be monitored by detection of a viral quasispecies of spike deletion variants (SDVs) (C. L. Rowe, S. C. Baker, M. J. Nathan, and J. O. Fleming, J. Virol. 71:2959-2969, 1997). We and others have found that the deletions cluster in the region from 1,200 to 1,800 nucleotides from the 5' end of the spike gene sequence, termed the "hypervariable" region. To address how SDVs might arise, we generated the predicted folding structures of the positive- and negative-strand senses of the entire 4,139-nt spike RNA sequence. We found that a prominent, isolated stem-loop structure is coincident with the hypervariable region in each structure. To determine if this predicted stem-loop is a "hot spot" for RNA recombination, we assessed whether this region of the spike is more frequently deleted than three other selected regions of the spike sequence in a population of viral sequences isolated from the CNS of acutely and persistently infected mice. Using differential colony hybridization of cloned spike reverse transcription-PCR products, we detected SDVs in which the hot spot was deleted but did not detect SDVs in which other regions of the spike sequence were exclusively deleted. Furthermore, sequence analysis and mapping of the crossover sites of 25 distinct patterns of SDVs showed that the majority of crossover sites clustered to two regions at the base of the isolated stem-loop, which we designated as high-frequency recombination sites 1 and 2. Interestingly, the majority of the left and right crossover sites of the SDVs were directly across from or proximal to one another, suggesting that these SDVs are likely generated by intramolecular recombination. Overall, our results are consistent with there being an important role for the spike RNA secondary structure as a contributing factor in the generation of SDVs during persistent infection.

Experimental confirmation of recombination upstream of the S1 hypervariable region of infectious bronchitis virus

Chimeric infectious bronchitis virus (IBV) genomes with cross-over sites in the S1 gene were generated by co-infection with two distinct IBV strains. Recombinant viruses were collected from chicken embryos, embryonic cultured cells and chickens co-infected with Ark99 and Mass41 strains and purified by differential centrifugation. The recombinant S1 genes were identified by reverse transcription polymerase chain reaction (RTPCR) using heterologous primers and confirmed by nucleotide sequencing. The recombinants with Ark99 5' and Mass41 3' sequences were identified following the in vitro, in ovo and in vivo co-infections. Mixed RNA extracted from Ark99 and Mass41 did not produce RTPCR products with these primers at the PCR conditions used. Cross-over sites within the amplified 580 (Mass41) or 604 (Ark99) bases of the 5' S1 gene could only be detected between nucleotides 50 and 155. While this region, lying upstream of the S1 hypervariable region, corresponded with sites commonly identified in naturally occurring isolates, recombination sites identified in these studies could not be detected within the HVR of S1 of the genomes of chimeric viruses.

A specific host cellular protein binding element near the 3' end of mouse hepatitis virus genomic RNA

A distinct host cellular protein binding element was mapped within a 38-nucleotide (nt) sequence 166-129 nucleotides upstream of the 3' end of the MHV-JHM genome using a RNase T1 protection/gel mobility shift electrophoresis assay. The resultant RNA-protein complex contains six host cellular proteins, one protein of 120-kDa molecular mass, two poorly resolved species approximately 55 kDa in size, a second pair of poorly resolved 40-kDa proteins, and a minor component of 25 kDa. A series of RNA probes containing deletions or clustered transversion mutations were tested for their ability to form complexes with mock- and MHV-JHM-infected cytoplasmic extracts. Three mutant RNA probes (mA, mB, and mC) with deletions at 154-140, 139-129, and 128-118, respectively, expressed 4, 37, and 94% of the host protein binding activity exhibited by the wild-type RNA. Defective interfering (DI) RNAs (DImA, DImB, and DImC) containing corresponding deletions at 154-140, 139-129, 128-118, and another DI RNA (DImD) with a deletion at nucleotides (nts) 112-102, a region which did not affect RNA-protein interactions, were transfected into MHV-JHM-infected 17CL-1 cells to assay the effects of these mutations on DI RNA replication. All of these mutations had an adverse effect on DI RNA replication. However, analysis of negative strand mutant DI RNAs revealed that two mutants (DImC and DImD) carrying deletions having little or no effect on RNA-protein interaction in our RNA-protein binding assays maintained their mutant sequences. In contrast, the other two mutants (DImA and DImB) containing deletions that dramatically decreased RNA-protein binding activity did not maintain their mutations; wild-type sequences were restored in the majority of the progeny negative strand molecules. These data indicate that the 26-nucleotide sequence at positions 154-129 from the 3' end of viral genome is important to both RNA-protein binding and viral replication. This protein binding element contains an 11-nt sequence (UGAGAGAAGUU, positions 139-129) very similar to a more 3' sequence (UGAAUGAAGUU) previously implicated in host protein binding and viral RNA replication (Yu and Leibowitz, 1995a and 1995b).

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.

Genetic engineering of animal RNA viruses

The ability to genetically manipulate viruses has led to extraordinary advances in understanding virus biology and to the establishment of useful vector systems. Initially confined to DNA viruses and retroviruses, RNA viruses have more recently become attractive candidates for expression of heterologous genes and offer promising perspectives for biomedical applications.

Nonhomologous RNA-RNA recombination events at the 3' nontranslated region of the Sindbis virus genome: hot spots and utilization of nonviral sequences

The mechanism of RNA-RNA recombination at the 3' nontranslated region (3'NTR) of the Sindbis virus (SIN) genome was studied by using nonreplicative RNA precursors. The 11.7-kb SIN genome was transcribed in vitro as two nonoverlapping RNA fragments. RNA-1 contained the entire 11.4-kb protein coding sequence of SIN and also carried an additional 1.8-kb nonviral sequence at its 3' end. RNA-2 carried the remaining 0.26 or 0.3 kb of the SIN genome containing the 3'NTR. Transfection of these two fragments into BHK cells resulted in vivo RNA-RNA recombination and release of infectious SIN recombinants. Eighteen plaque-purified recombinant viruses were sequenced to precisely map the RNA-RNA crossover sites at the 3'NTR. Sixteen of the 18 recombinants were found to be genetically heterogeneous at the 3'NTR. Two major clustered sites within the 3'NTR of RNA-2 were found to be fused to multiple locations on the nonviral sequence of RNA-1, resulting in insertions of 10 to 1,085 nucleotides at the 3'NTR. Sequence analysis of crossover sites suggested only limited homology and heteroduplex-forming capability between substrate RNAs. Analysis of additional 23 recombinant viruses generated by mutagenized donor and acceptor templates supports the occurrence of recombination hot spots on donor templates. Introduction of a 17-nucleotide rudimentary replicase recognition signal in the acceptor template alone did not induce the polymerase to reinitiate at the 17-nucleotide signal. Interestingly, deletion of a 24-nucleotide hot spot locus on the donor template abolished crossover events at one of the two sites and allowed the polymerase to reinitiate at the 17-nucleotide replicase recognition signal inserted at the acceptor template. The possible roles of RNA-protein and RNA-RNA interactions in the differential regulation of apparent pausing, template selection, and reinitiation are discussed.

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.