Role of nucleotides immediately flanking the transcription-regulating sequence core in coronavirus subgenomic mRNA synthesis

Regulation viral RNA transcription and replication by higher-order RNA structures within the nsp1 coding region of MERS coronavirus

Coronavirus (CoV) possesses numerous functional cis-acting elements in its positive-strand genomic RNA. Although most of these RNA structures participate in viral replication, the functions of RNA structures in the genomic RNA of CoV in viral replication remain unclear. In this study, we investigated the functions of the higher-order RNA stem-loop (SL) structures SL5B, SL5C, and SL5D in the ORF1a coding region of Middle East respiratory syndrome coronavirus (MERS-CoV) in viral replication. Our approach, using reverse genetics of a bacterial artificial chromosome system, revealed that SL5B and SL5C play essential roles in the discontinuous transcription of MERS-CoV. In silico analyses predicted that SL5C interacts with a bulged stem-loop (BSL) in the 3' untranslated region, suggesting that the RNA structure of SL5C is important for viral RNA transcription. Conversely, SL5D did not affect transcription, but mediated the synthesis of positive-strand genomic RNA. Additionally, the RNA secondary structure of SL5 in the revertant virus of the SL5D mutant was similar to that of the wild-type, indicating that the RNA structure of SL5D can finely tune RNA replication in MERS-CoV. Our data indicate novel regulatory mechanisms of viral RNA transcription and replication by higher-order RNA structures in the MERS-CoV genomic RNA.

m6A Methylation of Transcription Leader Sequence of SARS-CoV-2 Impacts Discontinuous Transcription of Subgenomic mRNAs

The SARS-CoV-2 genome has been shown to be m6A methylated at several positions in vivo. Strikingly, a DRACH motif, the recognition motif for adenosine methylation, resides in the core of the transcriptional regulatory leader sequence (TRS-L) at position A74, which is highly conserved and essential for viral discontinuous transcription. Methylation at position A74 correlates with viral pathogenicity. Discontinuous transcription produces a set of subgenomic mRNAs that function as templates for translation of all structural and accessory proteins. A74 is base-paired in the short stem-loop structure 5'SL3 that opens during discontinuous transcription to form long-range RNA-RNA interactions with nascent (-)-strand transcripts at complementary TRS-body sequences. A74 can be methylated by the human METTL3/METTL14 complex in vitro. Here, we investigate its impact on the structural stability of 5'SL3 and the long-range TRS-leader:TRS-body duplex formation necessary for synthesis of subgenomic mRNAs of all four viral structural proteins. Methylation uniformly destabilizes 5'SL3 and long-range duplexes and alters their relative equilibrium populations, suggesting that the m6A74 modification acts as a regulator for the abundance of viral structural proteins due to this destabilization.

Classification, replication, and transcription of Nidovirales

Nidovirales is one order of RNA virus, with the largest single-stranded positive sense RNA genome enwrapped with membrane envelope. It comprises four families (Arterividae, Mesoniviridae, Roniviridae, and Coronaviridae) and has been circulating in humans and animals for almost one century, posing great threat to livestock and poultry，as well as to public health. Nidovirales shares similar life cycle: attachment to cell surface, entry, primary translation of replicases, viral RNA replication in cytoplasm, translation of viral proteins, virion assembly, budding, and release. The viral RNA synthesis is the critical step during infection, including genomic RNA (gRNA) replication and subgenomic mRNAs (sg mRNAs) transcription. gRNA replication requires the synthesis of a negative sense full-length RNA intermediate, while the sg mRNAs transcription involves the synthesis of a nested set of negative sense subgenomic intermediates by a discontinuous strategy. This RNA synthesis process is mediated by the viral replication/transcription complex (RTC), which consists of several enzymatic replicases derived from the polyprotein 1a and polyprotein 1ab and several cellular proteins. These replicases and host factors represent the optimal potential therapeutic targets. Hereby, we summarize the Nidovirales classification, associated diseases, "replication organelle," replication and transcription mechanisms, as well as related regulatory factors.

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

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

A recombination-resistant genome for live attenuated and stable PEDV vaccines by engineering the transcriptional regulatory sequences

IMPORTANCE

Coronaviruses are important pathogens of humans and animals, and vaccine developments against them are imperative. Due to the ability to induce broad and prolonged protective immunity and the convenient administration routes, live attenuated vaccines (LAVs) are promising arms for controlling the deadly coronavirus infections. However, potential recombination events between vaccine and field strains raise a safety concern for LAVs. The porcine epidemic diarrhea virus (PEDV) remodeled TRS (RMT) mutant generated in this study replicated efficiently in both cell culture and in pigs and retained protective immunogenicity against PEDV challenge in pigs. Furthermore, the RMT PEDV was resistant to recombination and genetically stable. Therefore, RMT PEDV can be further optimized as a backbone for the development of safe LAVs.

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed "discontinuous transcription" that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

High temporal resolution Nanopore sequencing dataset of SARS-CoV-2 and host cell RNAs

BACKGROUND

Recent studies have disclosed the genome, transcriptome, and epigenetic compositions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the effect of viral infection on gene expression of the host cells. It has been demonstrated that, besides the major canonical transcripts, the viral genome also codes for noncanonical RNA molecules. While the structural characterizations have revealed a detailed transcriptomic architecture of the virus, the kinetic studies provided poor and often misleading results on the dynamics of both the viral and host transcripts due to the low temporal resolution of the infection event and the low virus/cell ratio (multiplicity of infection [MOI] = 0.1) applied for the infection. It has never been tested whether the alteration in the host gene expressions is caused by aging of the cells or by the viral infection.

FINDINGS

In this study, we used Oxford Nanopore's direct cDNA and direct RNA sequencing methods for the generation of a high-coverage, high temporal resolution transcriptomic dataset of SARS-CoV-2 and of the primate host cells, using a high infection titer (MOI = 5). Sixteen sampling time points ranging from 1 to 96 hours with a varying time resolution and 3 biological replicates were used in the experiment. In addition, for each infected sample, corresponding noninfected samples were employed. The raw reads were mapped to the viral and to the host reference genomes, resulting in 49,661,499 mapped reads (54,62 Gbs). The genome of the viral isolate was also sequenced and phylogenetically classified.

CONCLUSIONS

This dataset can serve as a valuable resource for profiling the SARS-CoV-2 transcriptome dynamics, the virus-host interactions, and the RNA base modifications. Comparison of expression profiles of the host gene in the virally infected and in noninfected cells at different time points allows making a distinction between the effect of the aging of cells in culture and the viral infection. These data can provide useful information for potential novel gene annotations and can also be used for studying the currently available bioinformatics pipelines.

A live-attenuated SARS-CoV-2 vaccine candidate with accessory protein deletions

We report a live-attenuated SARS-CoV-2 vaccine candidate with (i) re-engineered viral transcription regulator sequences and (ii) deleted open-reading-frames (ORF) 3, 6, 7, and 8 (∆3678). The ∆3678 virus replicates about 7,500-fold lower than wild-type SARS-CoV-2 on primary human airway cultures, but restores its replication on interferon-deficient Vero-E6 cells that are approved for vaccine production. The ∆3678 virus is highly attenuated in both hamster and K18-hACE2 mouse models. A single-dose immunization of the ∆3678 virus protects hamsters from wild-type virus challenge and transmission. Among the deleted ORFs in the ∆3678 virus, ORF3a accounts for the most attenuation through antagonizing STAT1 phosphorylation during type-I interferon signaling. We also developed an mNeonGreen reporter ∆3678 virus for high-throughput neutralization and antiviral testing. Altogether, the results suggest that ∆3678 SARS-CoV-2 may serve as a live-attenuated vaccine candidate and a research tool for potential biosafety level-2 use.

Prevention and Control of Porcine Epidemic Diarrhea: The Development of Recombination-Resistant Live Attenuated Vaccines

Porcine epidemic diarrhea (PED), causing up to 100% mortality in neonatal pigs, is a highly contagious enteric disease caused by PED virus (PEDV). The highly virulent genogroup 2 (G2) PEDV emerged in 2010 and has caused huge economic losses to the pork industry globally. It was first reported in the US in 2013, caused country-wide outbreaks, and posed tremendous hardship for many pork producers in 2013–2014. Vaccination of pregnant sows/gilts with live attenuated vaccines (LAVs) is the most effective strategy to induce lactogenic immunity in the sows/gilts and provide a passive protection via the colostrum and milk to suckling piglets against PED. However, there are still no safe and effective vaccines available after about one decade of endeavor. One of the biggest concerns is the potential reversion to virulence of an LAV in the field. In this review, we summarize the status and the major obstacles in PEDV LAV development. We also discuss the function of the transcriptional regulatory sequences in PEDV transcription, contributing to recombination, and possible strategies to prevent the reversion of LAVs. This article provides insights into the rational design of a promising LAV without safety issues.

Accurate Identification of Transcription Regulatory Sequences and Genes in Coronaviruses

Transcription regulatory sequences (TRSs), which occur upstream of structural and accessory genes as well as the 5′ end of a coronavirus genome, play a critical role in discontinuous transcription in coronaviruses. We introduce two problems collectively aimed at identifying these regulatory sequences as well as their associated genes. First, we formulate the TRS Identification problem of identifying TRS sites in a coronavirus genome sequence with prescribed gene locations. We introduce CORSID-A, an algorithm that solves this problem to optimality in polynomial time. We demonstrate that CORSID-A outperforms existing motif-based methods in identifying TRS sites in coronaviruses. Second, we demonstrate for the first time how TRS sites can be leveraged to identify gene locations in the coronavirus genome. To that end, we formulate the TRS and Gene Identification problem of simultaneously identifying TRS sites and gene locations in unannotated coronavirus genomes. We introduce CORSID to solve this problem, which includes a web-based visualization tool to explore the space of near-optimal solutions. We show that CORSID outperforms state-of-the-art gene finding methods in coronavirus genomes. Furthermore, we demonstrate that CORSID enables de novo identification of TRS sites and genes in previously unannotated coronavirus genomes. CORSID is the first method to perform accurate and simultaneous identification of TRS sites and genes in coronavirus genomes without the use of any prior information.

How the Replication and Transcription Complex Functions in Jumping Transcription of SARS-CoV-2

Background: Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although unprecedented efforts are underway to develop therapeutic strategies against this disease, scientists have acquired only a little knowledge regarding the structures and functions of the CoV replication and transcription complex (RTC). Ascertaining all the RTC components and the arrangement of them is an indispensably step for the eventual determination of its global structure, leading to completely understanding all of its functions at the molecular level. Results: The main results include: 1) hairpins containing the canonical and non-canonical NSP15 cleavage motifs are canonical and non-canonical transcription regulatory sequence (TRS) hairpins; 2) TRS hairpins can be used to identify recombination regions in CoV genomes; 3) RNA methylation participates in the determination of the local RNA structures in CoVs by affecting the formation of base pairing; and 4) The eventual determination of the CoV RTC global structure needs to consider METTL3 in the experimental design. Conclusions: In the present study, we proposed the theoretical arrangement of NSP12-15 and METTL3 in the global RTC structure and constructed a model to answer how the RTC functions in the jumping transcription of CoVs. As the most important finding, TRS hairpins were reported for the first time to interpret NSP15 cleavage, RNA methylation of CoVs and their association at the molecular level. Our findings enrich fundamental knowledge in the field of gene expression and its regulation, providing a crucial basis for future studies.

[SARS-CoV-2 with transcription regulatory sequence motif mutation poses a greater threat]

OBJECTIVE

To analyze the mutations in transcription regulatory sequences (TRSs) of coronaviruss (CoV) to provide the basis for exploring the patterns of SARS-CoV-2 transmission and outbreak.

METHODS

A combined evolutionary and molecular functional analysis of all sets of publicly available genomic data of viruses was performed.

RESULTS

A leader transcription regulatory sequence (TRS-L) usually comprises the first 60-70 nts of the 5' UTR in a CoV genome, and the body transcription regulatory sequences (TRS-Bs) are located immediately upstream of the genes other than ORF1a and 1b. In each CoV genome, the TRS-L and TRS-Bs share a specific consensus sequence, namely the TRS motif. Any changes of nucleotide residues in the TRS motifs are defined as TRS motif mutations. Mutations in the TRS-L or multiple TRS-Bs result in superattenuated variants. The spread of super-attenuated variants may cause an increase in asymptomatic or mild infections, prolonged incubation periods and a decreased detection rate of the viruses, thus posing new challenges to SARS-CoV-2 prevention and control. The super-attenuated variants also increase their possibility of long-term coexistence with humans. The Delta variant is significantly different from all the previous variants and may lead to a large-scale transmission. The Delta variant (B.1.617.2) with TRS motif mutation has already appeared and shown signs of spreading in Singapore, which, and even the Southeast Asia, may become the new epicenter of the next wave of SARS-CoV-2 outbreak.

CONCLUSION

TRS motif mutation will occur in all variants of SARS-CoV-2 and may result in super-attenuated variants. Only super-attenuated variants with TRS motif mutations will eventually lose the abilities of cross-species transmission and causing outbreaks.

Structures and functions of coronavirus replication-transcription complexes and their relevance for SARS-CoV-2 drug design

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed millions of people and continues to cause massive global upheaval. Coronaviruses are positive-strand RNA viruses with an unusually large genome of ~30 kb. They express an RNA-dependent RNA polymerase and a cohort of other replication enzymes and supporting factors to transcribe and replicate their genomes. The proteins performing these essential processes are prime antiviral drug targets, but drug discovery is hindered by our incomplete understanding of coronavirus RNA synthesis and processing. In infected cells, the RNA-dependent RNA polymerase must coordinate with other viral and host factors to produce both viral mRNAs and new genomes. Recent research aiming to decipher and contextualize the structures, functions and interplay of the subunits of the SARS-CoV-2 replication and transcription complex proteins has burgeoned. In this Review, we discuss recent advancements in our understanding of the molecular basis and complexity of the coronavirus RNA-synthesizing machinery. Specifically, we outline the mechanisms and regulation of RNA translation, replication and transcription. We also discuss the composition of the replication and transcription complexes and their suitability as targets for antiviral therapy.

Porcine Respiratory Coronavirus as a Model for Acute Respiratory Coronavirus Disease

In the light of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, we have developed a porcine respiratory coronavirus (PRCV) model for in depth mechanistic evaluation of the pathogenesis, virology and immune responses of this important family of viruses. Pigs are a large animal with similar physiology and immunology to humans and are a natural host for PRCV. Four PRCV strains were investigated and shown to induce different degrees of lung pathology. Importantly, although all four strains replicated equally well in porcine cell lines in vitro and in the upper respiratory tract in vivo, PRCV strains causing more severe lung pathology were also able to replicate in ex vivo tracheal organ cultures as well as in vivo in the trachea and lung. The time course of infection of PRCV 135, which caused the most severe pulmonary pathology, was investigated. Virus was shed from the upper respiratory tract until day 10 post infection, with infection of the respiratory mucosa, as well as olfactory and sustentacular cells, providing an excellent model to study upper respiratory tract disease in addition to the commonly known lower respiratory tract disease from PRCV. Infected animals made antibody and T cell responses that cross reacted with the four PRCV strains and Transmissible Gastroenteritis Virus. The antibody response was reproduced in vitro in organ cultures. Comparison of mechanisms of infection and immune control in pigs infected with PRCVs of differing pathogenicity with human data from SARS-CoV-2 infection and from our in vitro organ cultures, will enable key events in coronavirus infection and disease pathogenesis to be identified.

SARS-CoV-2, Cardiovascular Diseases, and Noncoding RNAs: A Connected Triad

Coronavirus Disease 2019 (COVID-19), caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is characterized by important respiratory impairments frequently associated with severe cardiovascular damages. Moreover, patients with pre-existing comorbidity for cardiovascular diseases (CVD) often present a dramatic increase in inflammatory cytokines release, which increases the severity and adverse outcomes of the infection and, finally, mortality risk. Despite this evident association at the clinical level, the mechanisms linking CVD and COVID-19 are still blurry and unresolved. Noncoding RNAs (ncRNAs) are functional RNA molecules transcribed from DNA but usually not translated into proteins. They play an important role in the regulation of gene expression, either in relatively stable conditions or as a response to different stimuli, including viral infection, and are therefore considered a possible important target in the design of specific drugs. In this review, we introduce known associations and interactions between COVID-19 and CVD, discussing the role of ncRNAs within SARS-CoV-2 infection from the perspective of the development of efficient pharmacological tools to treat COVID-19 patients and taking into account the equally dramatic associated consequences, such as those affecting the cardiovascular system.

Duplex formation between the template and the nascent strand in the transcription-regulating sequences is associated with the site of template switching in SARS - CoV-2

Recently published transcriptomic data of the SARS-CoV-2 coronavirus show that there is a large variation in the frequency and steady state levels of subgenomic mRNA sequences. This variation is derived from discontinuous subgenomic RNA synthesis, where the polymerase switches template from a 3' proximal genome body sequence to a 5' untranslated leader sequence. This leads to a fusion between the common 5' leader sequence and a 3' proximal body sequence in the RNA product. This process revolves around a common core sequence (CS) that is present at both the template sites that make up the fusion junction. Base-pairing between the leader CS and the nascent complementary minus strand body CS, and flanking regions (together called the transcription regulating sequence, TRS) is vital for this template switching event. However, various factors can influence the site of template switching within the same TRS duplex. Here, we model the duplexes formed between the leader and complementary body TRS regions, hypothesizing the role of the stability of the TRS duplex in determining the major sites of template switching for the most abundant mRNAs. We indicate that the stability of secondary structures and the speed of transcription play key roles in determining the probability of template switching in the production of subgenomic RNAs. We speculate on the effect of reported variant nucleotide substitutions on our models.

In vivo structure and dynamics of the SARS-CoV-2 RNA genome

The dynamics of SARS-CoV-2 RNA structure and their functional relevance are largely unknown. Here we develop a simplified SPLASH assay and comprehensively map the in vivo RNA-RNA interactome of SARS-CoV-2 genome across viral life cycle. We report canonical and alternative structures including 5'-UTR and 3'-UTR, frameshifting element (FSE) pseudoknot and genome cyclization in both cells and virions. We provide direct evidence of interactions between Transcription Regulating Sequences, which facilitate discontinuous transcription. In addition, we reveal alternative short and long distance arches around FSE. More importantly, we find that within virions, while SARS-CoV-2 genome RNA undergoes intensive compaction, genome domains remain stable but with strengthened demarcation of local domains and weakened global cyclization. Taken together, our analysis reveals the structural basis for the regulation of replication, discontinuous transcription and translational frameshifting, the alternative conformations and the maintenance of global genome organization during the whole life cycle of SARS-CoV-2, which we anticipate will help develop better antiviral strategies.

Biology and Pathogenesis of SARS-CoV-2: Understandings for Therapeutic Developments against COVID-19

Coronaviruses are positive sense, single-stranded, enveloped, and non-segmented RNA viruses that belong to the Coronaviridae family within the order Nidovirales and suborder Coronavirinae. Two Alphacoronavirus strains: HCoV-229E and HCoV-NL63 and five Betacoronaviruses: HCoV-HKU1, HCoV-OC43, SARS-CoV, MERS-CoV, and SARS-CoV-2 have so far been recognized as Human Coronaviruses (HCoVs). Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 is currently the greatest concern for humanity. Despite the overflow of research on SARS-CoV-2 and other HCoVs published every week, existing knowledge in this area is insufficient for the complete understanding of the viruses and the diseases caused by them. This review is based on the analysis of 210 published works, and it attempts to cover the basic biology of coronaviruses, including the genetic characteristics, life cycle, and host-pathogen interaction, pathogenesis, the antiviral drugs, and vaccines against HCoVs, especially focusing on SARS-CoV-2. Furthermore, we will briefly discuss the potential link between extracellular vesicles (EVs) and SARS-CoV-2/COVID-19 pathophysiology.

Generation of a Novel SARS-CoV-2 Sub-genomic RNA Due to the R203K/G204R Variant in Nucleocapsid: Homologous Recombination has Potential to Change SARS-CoV-2 at Both Protein and RNA Level

BACKGROUND

Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host's anti-viral immune response, in turn affecting the frequency of variants over time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally.

METHODS

Deep-sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. Results: Sequence analysis suggests that the 3 adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep-sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames.

CONCLUSIONS

The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans, resulting in both coding changes and novel sub-genomic RNA transcripts, suggests this as a mechanism for diversification and adaptation within its new host.

Generation of a novel SARS-CoV-2 sub-genomic RNA due to the R203K/G204R variant in nucleocapsid: homologous recombination has potential to change SARS-CoV-2 at both protein and RNA level

BACKGROUND

Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host’s anti-viral immune response, in turn affecting the frequency of variants over-time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally.

METHODS

Deep sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome.

RESULTS

Sequence analysis suggests that the three adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence (CS) of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames.

CONCLUSIONS

The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans resulting in both coding changes and novel sub-genomic RNA transcripts suggests this as a mechanism for diversification and adaptation within its new host.

The strand-biased transcription of SARS-CoV-2 and unbalanced inhibition by remdesivir

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a positive single-stranded RNA virus, causes the coronavirus disease 19 pandemic. During the viral replication and transcription, the RNA-dependent RNA polymerase "jumps" along the genome template, resulting in discontinuous negative-stranded transcripts. Although the sense-mRNA architectures of SARS-CoV-2 were reported, its negative strand was unexplored. Here, we deeply sequenced both strands of RNA and found SARS-CoV-2 transcription is strongly biased to form the sense strand with variable transcription efficiency for different genes. During negative strand synthesis, numerous non-canonical fusion transcripts are also formed, driven by 3-15 nt sequence homology scattered along the genome but more prone to be inhibited by SARS-CoV-2 RNA polymerase inhibitor remdesivir. The drug also represses more of the negative than the positive strand synthesis as supported by a mathematic simulation model and experimental quantifications. Overall, this study opens new sights into SARS-CoV-2 biogenesis and may facilitate the antiviral vaccine development and drug design.

Analysis of Emerging Variants in Structured Regions of the SARS-CoV-2 Genome

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two sets of SARS-CoV-2 genomic sequences to identify emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, 20 variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5' untranslated region (UTR), including a group of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-ss molecular switch in the 3'UTR. Finally, 5 variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M (stem loop 2 m) selfish genetic element, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. Our analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences.

Characterizing Transcriptional Regulatory Sequences in Coronaviruses and Their Role in Recombination

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

A review of novel coronavirus disease (COVID-19): based on genomic structure, phylogeny, current shreds of evidence, candidate vaccines, and drug repurposing

Coronavirus disease (COVID-19) pandemic is instigated by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). As of March 13, 2021, more than 118.9 million cases were infected with COVID-19 worldwide. SARS-CoV-2 is a positive-sense single-stranded RNA beta-CoV. Most COVID-19 infected individuals recover within 1-3 weeks. Nevertheless, approximately 5% of patients develop acute respiratory distress syndrome and other systemic complications, leading to death. Structural genetic analyses of SARS-CoV-2 have shown genomic resemblances but a low evolutionary correlation to SARS-CoV-1 responsible for the 2002-2004 outbreak. The S glycoprotein is critical for cell adhesion and the entrance of the virus into the host. The process of cell entry uses the cellular receptor named angiotensin-converting enzyme 2. Recent evidence proposed that the CD147 as a SARS-CoV-2's potential receptor. The viral genome is mainly held by two non-structural proteins (NSPs), ORF1a and ORF1ab, along with structural proteins. Although NSPs are conserved among the βCoVs, mutations in NSP2 and NSP3 may play critical roles in transmitting the virus and cell tropism. To date, no specific/targeted anti-viral treatments exist. Notably, more than 50 COVID-19 candidate vaccines in clinical trials, and a few being administered. Preventive precautions are the primary strategy to limit the viral load transmission and spread, emphasizing the urgent need for developing significant drug targets and vaccines against COVID-19. This review provides a cumulative overview of the genomic structure, transmission, phylogeny of SARS-CoV-2 from Indian clusters, treatment options, updated discoveries, and future standpoints for COVID-19.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02749-0.

Clinical and molecular aspects of veterinary coronaviruses

Coronaviruses are a large group of RNA viruses that infect a wide range of animal species. The replication strategy of coronaviruses involves recombination and mutation events that lead to the possibility of cross-species transmission. The high plasticity of the viral receptor due to a continuous modification of the host species habitat may be the cause of cross-species transmission that can turn into a threat to other species including the human population. The successive emergence of highly pathogenic coronaviruses such as the Severe Acute Respiratory Syndrome (SARS) in 2003, the Middle East Respiratory Syndrome Coronavirus in 2012, and the recent SARS-CoV-2 has incentivized a number of studies on the molecular basis of the coronavirus and its pathogenesis. The high degree of interrelatedness between humans and wild and domestic animals and the modification of animal habitats by human urbanization, has favored new viral spreads. Hence, knowledge on the main clinical signs of coronavirus infection in the different hosts and the distinctive molecular characteristics of each coronavirus is essential to prevent the emergence of new coronavirus diseases. The coronavirus infections routinely studied in veterinary medicine must be properly recognized and diagnosed not only to prevent animal disease but also to promote public health.

Rapid Spread of Mutant Alleles in Worldwide SARS-CoV-2 Strains Revealed by Genome-Wide Single Nucleotide Polymorphism and Variation Analysis

The novel coronavirus (SARS-CoV-2) has become a pandemic and is threatening human health globally. Here, we report nine newly evolved SARS-CoV-2 single nucleotide polymorphism (SNP) alleles those underwent a rapid increase (seven cases) or decrease (two cases) in their frequency for 30-80% in the initial four months, which are further confirmed by intrahost single nucleotide variation analysis using raw sequence data including 8,217 samples. The nine SNPs are mostly (8/9) located in the coding region and are mainly (6/9) nonsynonymous substitutions. The nine SNPs show a complete linkage in SNP pairs and belong to three different linkage groups, named LG_1 to LG_3. Analyses in population genetics show signatures of adaptive selection toward the mutants in LG_1, but no signal of selection for LG_2. Population genetic analysis results on LG_3 show geological differentiation. Analyses on geographic COVID-19 cases and published clinical data provide evidence that the mutants in LG_1 and LG_3 benefit virus replication and those in LG_1 have a positive correlation with the disease severity in COVID-19-infected patients. The mutants in LG_2 show a bias toward mildness of the disease based on available public clinical data. Our findings may be instructive for epidemiological surveys and disease control of COVID-19 in the future.

Multiple novel non-canonically transcribed sub-genomic mRNAs produced by avian coronavirus infectious bronchitis virus

Coronavirus sub-genomic mRNA (sgmRNA) synthesis occurs via a process of discontinuous transcription involving complementary transcription regulatory sequences (TRSs), one (TRS-L) encompassing the leader sequence of the 5' untranslated region (UTR), and the other upstream of each structural and accessory gene (TRS-B). Several coronaviruses have an ORF located between the N gene and the 3'-UTR, an area previously thought to be non-coding in the Gammacoronavirus infectious bronchitis virus (IBV) due to a lack of a canonical TRS-B. Here, we identify a non-canonical TRS-B allowing for a novel sgmRNA relating to this ORF to be produced in several strains of IBV: Beaudette, CR88, H120, D1466, Italy-02 and QX. Interestingly, the potential protein produced by this ORF is prematurely truncated in the Beaudette strain. A single nucleotide deletion was made in the Beaudette strain allowing for the generation of a recombinant IBV (rIBV) that had the potential to express a full-length protein. Assessment of this rIBV in vitro demonstrated that restoration of the full-length potential protein had no effect on viral replication. Further assessment of the Beaudette-derived RNA identified a second non-canonically transcribed sgmRNA located within gene 2. Deep sequencing analysis of allantoic fluid from Beaudette-infected embryonated eggs confirmed the presence of both the newly identified non-canonically transcribed sgmRNAs and highlighted the potential for further yet unidentified sgmRNAs. This HiSeq data, alongside the confirmation of non-canonically transcribed sgmRNAs, indicates the potential of the coronavirus genome to encode a larger repertoire of genes than has currently been identified.

Translation-Associated Mutational U-Pressure in the First ORF of SARS-CoV-2 and Other Coronaviruses

Within 4 months of the ongoing COVID-19 pandemic caused by SARS-CoV-2, more than 250 nucleotide mutations have been detected in ORF1ab of the virus isolated from infected persons from different parts of the globe. These observations open up an obvious question about the rate and direction of mutational pressure for further vaccine and therapeutics designing. In this study, we did a comparative analysis of ORF1a and ORF1b by using the first isolate (Wuhan strain) as the parent sequence. We observed that most of the nucleotide mutations are C to U transitions. The rate of synonymous C to U transitions is significantly higher than the rate of non-synonymous ones, indicating negative selection on amino acid substitutions. Further, trends in nucleotide usage bias have been investigated in 49 coronaviruses species. A strong bias in nucleotide usage in fourfold degenerate sites toward uracil residues is seen in ORF1ab of all the studied coronaviruses: both in the ORF1a and in the ORF1b translated thanks to the programmed ribosomal frameshifting that has an efficiency of 14 – 45% in different species. A more substantial mutational U-pressure is observed in ORF1a than in ORF1b perhaps because ORF1a is translated more frequently than ORF1b. Mutational U-pressure is there even in ORFs that are not translated from genomic RNA plus strands, but the bias is weaker than in ORF1ab. Unlike other nucleotide mutations, mutational U-pressure caused by cytosine deamination, mostly occurring during the RNA plus strand replication and also translation, cannot be corrected by the proof-reading machinery of coronaviruses. The knowledge generated on the mutational U-pressure that becomes stronger during translation of viral RNA plus strands has implications for vaccine and nucleoside analog development for treating COVID-19 and other coronavirus infections.

The molecular virology of coronaviruses

Few human pathogens have been the focus of as much concentrated worldwide attention as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of COVID-19. Its emergence into the human population and ensuing pandemic came on the heels of severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV), two other highly pathogenic coronavirus spillovers, which collectively have reshaped our view of a virus family previously associated primarily with the common cold. It has placed intense pressure on the collective scientific community to develop therapeutics and vaccines, whose engineering relies on a detailed understanding of coronavirus biology. Here, we present the molecular virology of coronavirus infection, including its entry into cells, its remarkably sophisticated gene expression and replication mechanisms, its extensive remodeling of the intracellular environment, and its multifaceted immune evasion strategies. We highlight aspects of the viral life cycle that may be amenable to antiviral targeting as well as key features of its biology that await discovery.

Analysis of Rapidly Emerging Variants in Structured Regions of the SARS-CoV-2 Genome

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two cohorts of SARS-CoV-2 genomic sequences to identify rapidly emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, twenty variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5'UTR, including a set of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-stable molecular switch in the 3'UTR. Finally, five variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M stem loop, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. This analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences, as rapidly emerging variations in these regions could lead to drug resistance.

Characterizing transcriptional regulatory sequences in coronaviruses and their role in recombination

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

Comparative analysis of coronavirus genomic RNA structure reveals conservation in SARS-like coronaviruses

Coronaviruses, including SARS-CoV-2 the etiological agent of COVID-19 disease, have caused multiple epidemic and pandemic outbreaks in the past 20 years1-3. With no vaccines, and only recently developed antiviral therapeutics, we are ill equipped to handle coronavirus outbreaks4. A better understanding of the molecular mechanisms that regulate coronavirus replication and pathogenesis is needed to guide the development of new antiviral therapeutics and vaccines. RNA secondary structures play critical roles in multiple aspects of coronavirus replication, but the extent and conservation of RNA secondary structure across coronavirus genomes is unknown5. Here, we define highly structured RNA regions throughout the MERS-CoV, SARS-CoV, and SARS-CoV-2 genomes. We find that highly stable RNA structures are pervasive throughout coronavirus genomes, and are conserved between the SARS-like CoV. Our data suggests that selective pressure helps preserve RNA secondary structure in coronavirus genomes, suggesting that these structures may play important roles in virus replication and pathogenesis. Thus, disruption of conserved RNA secondary structures could be a novel strategy for the generation of attenuated SARS-CoV-2 vaccines for use against the current COVID-19 pandemic.

Comparative Analysis of Gene Expression in Virulent and Attenuated Strains of Infectious Bronchitis Virus at Subcodon Resolution

Like all coronaviruses, avian infectious bronchitis virus (IBV) possesses a long, single-stranded, positive-sense RNA genome (∼27 kb) and has a complex replication strategy that includes the production of a nested set of subgenomic mRNAs (sgmRNAs). Here, we used whole-transcriptome sequencing (RNASeq) and ribosome profiling (RiboSeq) to delineate gene expression in the IBV M41-CK and Beau-R strains at subcodon resolution. RNASeq facilitated a comparative analysis of viral RNA synthesis and revealed two novel transcription junction sites in the attenuated Beau-R strain, one of which would generate a sgmRNA encoding a ribosomally occupied open reading frame (dORF) located downstream of the nucleocapsid coding region. RiboSeq permitted quantification of the translational efficiency of virus gene expression and identified, for the first time, sites of ribosomal pausing on the genome. Quantification of reads flanking the programmed ribosomal frameshifting (PRF) signal at the genomic RNA ORF1a/ORF1b junction revealed that PRF in IBV is highly efficient (33 to 40%). Triplet phasing of RiboSeq data allowed precise determination of reading frames and revealed the translation of two ORFs (ORF4b and ORF4c on sgmRNA IR), which are widely conserved across IBV isolates. Analysis of differential gene expression in infected primary chick kidney cells indicated that the host cell response to IBV occurs primarily at the level of transcription, with global upregulation of immune-related mRNA transcripts following infection and comparatively modest changes in the translation efficiencies of host genes. Cellular genes and gene networks differentially expressed during virus infection were also identified, giving insights into the host cell response to IBV infection.IMPORTANCE IBV is a major avian pathogen and presents a substantial economic burden to the poultry industry. Improved vaccination strategies are urgently needed to curb the global spread of this virus, and the development of suitable vaccine candidates will be aided by an improved understanding of IBV molecular biology. Our high-resolution data have enabled a precise study of transcription and translation in cells infected with both pathogenic and attenuated forms of IBV and expand our understanding of gammacoronaviral gene expression. We demonstrate that gene expression shows considerable intraspecies variation, with single nucleotide polymorphisms being associated with altered production of sgmRNA transcripts, and our RiboSeq data sets enabled us to uncover novel ribosomally occupied ORFs in both strains. The numerous cellular genes and gene networks found to be differentially expressed during virus infection provide insights into the host cell response to IBV infection.

Interaction of coronavirus nucleocapsid protein with the 5'- and 3'-ends of the coronavirus genome is involved in genome circularization and negative-strand RNA synthesis

Synthesis of the negative‐strand ((−)‐strand) counterpart is the first step of coronavirus (CoV) replication; however, the detailed mechanism of the early event and the factors involved remain to be determined. Here, using bovine coronavirus (BCoV)‐defective interfering (DI) RNA, we showed that (a) a poly(A) tail with a length of 15 nucleotides (nt) was sufficient to initiate efficient (−)‐strand RNA synthesis and (b) substitution of the poly(A) tail with poly(U), (C) or (G) only slightly decreased the efficiency of (−)‐strand synthesis. The findings indicate that in addition to the poly(A) tail, other factors acting in trans may also participate in (−)‐strand synthesis. The BCoV nucleocapsid (N) protein, an RNA‐binding protein, was therefore tested as a candidate. Based on dissociation constant (K_d) values, it was found that the binding affinity between N protein, but not poly(A)‐binding protein, and the 3′‐terminal 55 nt plus a poly(A), poly(U), poly(C) or poly(G) tail correlates with the efficiency of (−)‐strand synthesis. Such an association was also evidenced by the binding affinity between the N protein and 5′‐ and 3′‐terminal cis‐acting elements important for (−)‐strand synthesis. Further analysis demonstrated that N protein can act as a bridge to facilitate interaction between the 5′‐ and 3′‐ends of the CoV genome, leading to circularization of the genome. Together, the current study extends our understanding of the mechanism of CoV (−)‐strand RNA synthesis through involvement of N protein and genome circularization and thus may explain why the addition of N protein in trans is required for efficient CoV replication.

Identification and formation mechanism of a novel noncoding RNA produced by avian infectious bronchitis virus

Viral noncoding (nc) RNAs have been shown to play important roles in viral life cycle. Many viruses employ different mechanism to produce ncRNAs. Here, we report that coronavirus infectious bronchitis virus (IBV) produces a novel ncRNA in virus-infected cells. This ncRNA consists of 563 nucleotides excluding a poly(A) tail, is mainly derived from the 3'-untranslated region of IBV genome, and contains a 63-nt-long of terminal leader sequence derived from the 5' end of the viral genome. Using mutagenesis and reverse genetics, we reveal that this ncRNA is a subgenomic RNA generated by discontinuous transcription mechanism.

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

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

Molecular characteristic and pathogenicity analysis of a virulent recombinant avain infectious bronchitis virus isolated in China

A virulent infectious bronchitis virus (IBV), designated as CK/CH/GD/QY16 (referred as QY16), was isolated from a diseased chicken farm in Guangdong province, China, in 2016. The complete genome of the strain was sequenced and analyzed. The results show that the genome of QY16 consists of 27,670 nucleotides, excluding poly (A) tail, and that its genome organization is 5’ UTR-1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3’ UTR-poly (A) tail. Sequence comparison among QY16 and other IBV strains was conducted and its results demonstrate that the S1 gene of QY16 has the highest nucleotide sequence identity with that of 4/91, and the other part of its genome is highly similar to that of YX10. The results of the phylogenic analysis show that the entire genome of QY16 and most of the QY16 genes are located in the same cluster as those of YX10, except for the S1 gene which is located in the same cluster with that of 4/91. It has been further confirmed by the RDP and SimPlot analysis that QY16 is a recombinant strain deriving from YX10 (as the major parental sequence) and 4/91 (as the minor parental sequence), and that the recombination occurs in a region which includes the 3’-terminal 1b sequence (85 nt) and the 5’-terminal S1 protein gene sequence (1,466 nt). The results of the vaccination-challenge test suggest that QY16 is a nephropathogenic strain of IBV and that the vaccine strains–H120 and 4/91—cannot provide effective protection against it. These results indicate that the continuing evolution of IBV strains by genetic drift and genetic recombination may lead to IBV outbreaks even among the vaccinated chickens in China.

Evaluation of a recombination-resistant coronavirus as a broadly applicable, rapidly implementable vaccine platform

Emerging and re-emerging zoonotic viral diseases are major threats to global health, economic stability, and national security. Vaccines are key for reducing coronaviral disease burden; however, the utility of live-attenuated vaccines is limited by risks of reversion or repair. Because of their history of emergence events due to their prevalence in zoonotic pools, designing live-attenuated coronavirus vaccines that can be rapidly and broadly implemented is essential for outbreak preparedness. Here, we show that coronaviruses with completely rewired transcription regulatory networks (TRNs) are effective vaccines against SARS-CoV. The TRN-rewired viruses are attenuated and protect against lethal SARS-CoV challenge. While a 3-nt rewired TRN reverts via second-site mutation upon serial passage, a 7-nt rewired TRN is more stable, suggesting that a more extensively rewired TRN might be essential for avoiding growth selection. In summary, rewiring the TRN is a feasible strategy for limiting reversion in an effective live-attenuated coronavirus vaccine candidate that is potentially portable across the Nidovirales order.

Attenuation of replication by a 29 nucleotide deletion in SARS-coronavirus acquired during the early stages of human-to-human transmission

A 29 nucleotide deletion in open reading frame 8 (ORF8) is the most obvious genetic change in severe acute respiratory syndrome coronavirus (SARS-CoV) during its emergence in humans. In spite of intense study, it remains unclear whether the deletion actually reflects adaptation to humans. Here we engineered full, partially deleted (-29 nt), and fully deleted ORF8 into a SARS-CoV infectious cDNA clone, strain Frankfurt-1. Replication of the resulting viruses was compared in primate cell cultures as well as Rhinolophus bat cells made permissive for SARS-CoV replication by lentiviral transduction of the human angiotensin-converting enzyme 2 receptor. Cells from cotton rat, goat, and sheep provided control scenarios that represent host systems in which SARS-CoV is neither endemic nor epidemic. Independent of the cell system, the truncation of ORF8 (29 nt deletion) decreased replication up to 23-fold. The effect was independent of the type I interferon response. The 29 nt deletion in SARS-CoV is a deleterious mutation acquired along the initial human-to-human transmission chain. The resulting loss of fitness may be due to a founder effect, which has rarely been documented in processes of viral emergence. These results have important implications for the retrospective assessment of the threat posed by SARS.

Insertion position as well as the inserted TRS and gene sequences differentially affect the retention of foreign gene expression by simian hemorrhagic fever virus (SHFV)

Recombinant SHFV infectious cDNA clones expressing a foreign gene from an additional sg mRNA were constructed. Two 3' genomic region sites, between ORF4' and ORF2b and between ORF4 and ORF5, were utilized for insertion of the myxoma M013 gene with a C-terminal V5 tag followed by one of the three inserted transcription regulatory sequences (TRS), TRS2', TRS4' or TRS7. M013 insertion at the ORF4'/ORF2b site but not the ORF4/ORF5 site generated progeny virus but only the recombinant virus with an inserted TRS2' retained the entire M013 gene through passage four. Insertion of an auto-fluorescent protein gene, iLOV, with an inserted TRS2' at the ORF4'/ORF2b site, generated viable progeny virus. iLOV expression was maintained through passage eight. Although regulation of SHFV subgenomic RNA synthesis is complex, the ORF4'/ORF2b site, which is located between the two sets of minor structural proteins, is able to tolerate foreign gene insertion.

New insights about the regulation of Nidovirus subgenomic mRNA synthesis

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

Expanded subgenomic mRNA transcriptome and coding capacity of a nidovirus

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

Transgene expression in the genome of Middle East respiratory syndrome coronavirus based on a novel reverse genetics system utilizing Red-mediated recombination cloning

Middle East respiratory syndrome coronavirus (MERS-CoV) is a high-priority pathogen in pandemic preparedness research. Reverse genetics systems are a valuable tool to study viral replication and pathogenesis, design attenuated vaccines and create defined viral assay systems for applications such as antiviral screening. Here we present a novel reverse genetics system for MERS-CoV that involves maintenance of the full-length viral genome as a cDNA copy inserted in a bacterial artificial chromosome amenable to manipulation by homologue recombination, based on the bacteriophage λ Red recombination system. Based on a full-length infectious MERS-CoV cDNA clone, optimal genomic insertion sites and expression strategies for GFP were identified and used to generate a reporter MERS-CoV expressing GFP in addition to the complete set of viral proteins. GFP was genetically fused to the N-terminal part of protein 4a, from which it is released during translation via porcine teschovirus 2A peptide activity. The resulting reporter virus achieved titres nearly identical to the wild-type virus 48 h after infection of Vero cells at m.o.i. 0.001 (1×10⁵ p.f.u. ml^-1 and 3×10⁵ p.f.u. ml^-1, respectively), and allowed determination of the 50 % inhibitory concentration for the known MERS-CoV inhibitor cyclosporine A based on fluorescence readout. The resulting value was 2.41 µM, which corresponds to values based on wild-type virus. The reverse genetics system described herein can be efficiently mutated by Red-mediated recombination. The GFP-expressing reporter virus contains the full set of MERS-CoV proteins and achieves wild-type titres in cell culture.

Detection and phylogenetic analysis of porcine epidemic diarrhea virus in central China based on the ORF3 gene and the S1 gene

BACKGROUND

Porcine epidemic diarrhea (PED) has increased in severity in China since 2010. To investigate further the infectivity, genetic diversity and molecular epidemiology of its causative agent, the porcine epidemic diarrhea virus (PEDV), we assessed 129 clinical samples, which were the intestinal tissue of piglets with severe diarrhea, from 17 cities in central China. Both the spike (S) glycoprotein (S1, 1-789 amino acids (aa)) and the full-length ORF3 gene of 21 representative field strains from 21 farms in 11 cities were sequenced and analysed.

METHODS

PEDV was detected by reverse transcription-polymerase chain reaction (RT-PCR), and S1 and ORF3 sequences were processed by the Clustal W method via DNAMAN 8 software, and phylogenetic trees were constructed by the neighbor-joining method using MEGA 6 software.

RESULTS

The prevalence of PEDV was 92.25% and was detected in 119 of 129 samples, with 94.03% (63 of 67) of pig farms harbouring the disease. According to the phylogenetic analysis of the S1 genes, our isolates all fell into group G2 (variants) and showed a close relationship to isolates from Chinese (HN1303, CH/ZMDZY/11 and AJ1102), Korean (AD01), American (MN, IA1, IA2 and 13-019349) sources, and these isolates differed genetically from other Chinese (LZC, CH/HNZZ/2011 and SD-M) and Korean (SM98) strains as well Japanese (83-P5 and MK) strains. In addition, our isolates differed from attenuated vaccine strains, CV777 (used in China) and DR13 (used in Korea). According to our derived amino acid sequence analysis, we detected one novel variant PEDV, viz: CH/HNLY, with 4-aa insertion/deletion (RSSS/T) at position 375 and 1-aa (D) deletion at position 430 compared to the CV777 attenuated strain. These mutations were located on the receptor binding domain. Our ORF3 gene analyses showed that the prevalent PEDV isolates were variants, and the isolated strains differed genetically from the vaccine strains.

CONCLUSIONS

These findings illustrated the existence of genetic diversity among geographically distinct PEDV strains, and our study has provided an impetus to conduct further research on the PEDV receptor binding protein and on the new and efficacious vaccines design.

Continuous and Discontinuous RNA Synthesis in Coronaviruses

Replication of the coronavirus genome requires continuous RNA synthesis, whereas transcription is a discontinuous process unique among RNA viruses. Transcription includes a template switch during the synthesis of subgenomic negative-strand RNAs to add a copy of the leader sequence. Coronavirus transcription is regulated by multiple factors, including the extent of base-pairing between transcription-regulating sequences of positive and negative polarity, viral and cell protein-RNA binding, and high-order RNA-RNA interactions. Coronavirus RNA synthesis is performed by a replication-transcription complex that includes viral and cell proteins that recognize cis-acting RNA elements mainly located in the highly structured 5' and 3' untranslated regions. In addition to many viral nonstructural proteins, the presence of cell nuclear proteins and the viral nucleocapsid protein increases virus amplification efficacy. Coronavirus RNA synthesis is connected with the formation of double-membrane vesicles and convoluted membranes. Coronaviruses encode proofreading machinery, unique in the RNA virus world, to ensure the maintenance of their large genome size.

Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea viruses associated with outbreaks of severe diarrhea in piglets in Jiangxi, China 2013

Porcine epidemic diarrhea (PED), caused by porcine epidemic diarrhea virus (PEDV), is a highly contagious, acute enteric viral disease of swine characterized by vomiting, watery diarrhea, dehydration and death. To identify and characterize the field PEDVs associated with the outbreaks of severe diarrhea in piglets in Jiangxi, 2013, the complete genome sequences of two representative strains of PEDV, designated CH/JX-1/2013 and CH/JX-2/2013, were determined and analyzed. The genome sequences of both emergent Jiangxi PEDV strains, CH/JX-1/2013 and CH/JX-2/2013, were 28,038 nucleotides in length excluding 3’ poly (A) tail. Compared to the PEDV CV777 strain, CH/JX-1/2013 and CH/JX-2/2013 had some unique genetic characteristics in the proximal region of the 5´-UTRs. Phylogenetic analysis of the complete genomes and the structural proteins revealed that CH/JX-1/2013 and CH/JX-2/2013 had a close relationship with post-2010 Chinese PEDV strains and US strains identified in 2013. The nucleotide identity between the two Jiangxi strains (CH/JX-1/2013 and CH/JX-2/2013) and 30 strains of PEDV identified ante-2010 and post-2010 ranged from 96.3–97.0% and 97.3–99.7%, respectively. Multiple nucleotide and deduced amino acid mutations were observed in the ORF1a/b, S, ORF3, E, M and N genes among the current field PEDV strains when compared to the CV777 strain. Some of the mutations altered the amino acid charge and hydrophilicity, and notably, there was an amino acid substitution in the middle of one neutralizing epitope (L1371I) of the S gene of both CH/JX-1/2013 and CH/JX-2/2013. Taken together, the accumulated genetic variations of the current field PEDV strains might have led to antigenic changes of the viruses, which might confer the less effectiveness or failure of the CV777-based vaccines currently being widely used in Jiangxi, China.

The structure and functions of coronavirus genomic 3' and 5' ends

Coronaviruses (CoVs) are an important cause of illness in humans and animals. Most human coronaviruses commonly cause relatively mild respiratory illnesses; however two zoonotic coronaviruses, SARS-CoV and MERS-CoV, can cause severe illness and death. Investigations over the past 35 years have illuminated many aspects of coronavirus replication. The focus of this review is the functional analysis of conserved RNA secondary structures in the 5' and 3' of the betacoronavirus genomes. The 5' 350 nucleotides folds into a set of RNA secondary structures which are well conserved, and reverse genetic studies indicate that these structures play an important role in the discontinuous synthesis of subgenomic RNAs in the betacoronaviruses. These cis-acting elements extend 3' of the 5'UTR into ORF1a. The 3'UTR is similarly conserved and contains all of the cis-acting sequences necessary for viral replication. Two competing conformations near the 5' end of the 3'UTR have been shown to make up a potential molecular switch. There is some evidence that an association between the 3' and 5'UTRs is necessary for subgenomic RNA synthesis, but the basis for this association is not yet clear. A number of host RNA proteins have been shown to bind to the 5' and 3' cis-acting regions, but the significance of these in viral replication is not clear. Two viral proteins have been identified as binding to the 5' cis-acting region, nsp1 and N protein. A genetic interaction between nsp8 and nsp9 and the region of the 3'UTR that contains the putative molecular switch suggests that these two proteins bind to this region.

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

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