Transcription regulatory sequences and mRNA expression levels in the coronavirus transmissible gastroenteritis virus

Strenuous expression of porcine epidemic diarrhea virus ORF3 protein suggests host resistance

Porcine epidemic diarrhea virus is attenuated upon adaptation to cell culture. Exclusively genomic mutations have been traced to the ORF3 gene of the laboratory strains. Previous attempts to express the protein were unsuccessful. We sought to express the ORF3 protein in both mammalian and bacteria cells as a prerequisite for investigation of the protein's role. For prokaryotic expression, two vector systems, pET28-a(+) and pGEX-4T-1 were constructed and expressed in Escherichia coli cells. For eukaryotic analyses, ORF3/pEGFP-C1 vector constructs were expressed in human embryonic, green monkey kidney and mouse fibrous cells. Intriguingly, there was minimal expression of the ORF3 gene. Following a documented hint that truncated ORF3 revealed higher expression, ORF3 gene was truncated. The simple modular architecture research tool analysis predicted two transmembrane domains between amino acid (aa) 41-63 and aa 76-98. Consequently, we generated two fragments; ORF-N (aa 1-98) inclusive of transmembrane domains and ORF3-C (aa 99-224). These truncated sequences were constructed as the whole gene here referred to as ORF3 wild type (wt). Coomassie blue stained gels revealed bands of ORF3-C expressed as a fusion protein of 17.5 and 39 kDa in pET28-a(+) and pGEX-4T-1 vectors, respectively. In contrast, ORF3-N was not. Additionally, ORF3-N induction decreased total cellular proteins suggesting inhibition of protein synthesis or metabolism. Solubility tests carried out at 30 °C, 25 °C and 18 °C showed that ORF3 formed inclusion bodies. Similar findings were observed in mammalian cells. Noteworthy, morphological distortions appeared in mammalian cells expressing ORF3 protein or its truncated mutants suggesting significance in host viability.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Porcine respiratory coronavirus genome sequences; comparisons and relationships to transmissible gastroenteritis viruses

Porcine respiratory coronavirus (PRCV) was initially detected in Europe, and later in the United States of America (US), in the 1980s. In this study we obtained and compared PRCV sequences from Europe and the US, and investigated how these are related to transmissible gastroenteritis virus (TGEV) sequences. The whole genome sequences of Danish (1/90-DK), Italian (PRCV15087/12 III NPTV Parma), and Belgian PRCV (91V44) strains are presented. These sequences were aligned with nine other PRCV sequences from Europe and the US, and 43 TGEV sequences. Following alignment of the PRCV sequences, it was apparent that multiple amino acid variations in the structural proteins were distinct between the European and US strains. The alignments were used to build phylogenetic trees to infer the evolutionary relationships between the strains. In these trees, the European PRCV strains clustered as a separate group, whereas the US strains of PRCV all clustered with TGEVs.

Development of a colloidal gold immunochromatographic strip for the simultaneous detection of porcine epidemic diarrhea virus and transmissible gastroenteritis virus

In recent years, porcine diarrhea-associated viruses have caused significant economic losses globally. These viruses present similar clinical symptoms, such as watery diarrhea, dehydration, and vomiting. Co-infections with porcine epidemic diarrhea virus (PEDV) and transmissible gastroenteritis virus (TGEV) are common. For the rapid and on-site preliminary diagnosis on the pig farms, this study aimed to develop a colloidal gold immunochromatography assay (GICA) strip for the detection of PEDV and TGEV simultaneously. The GICA kit showed that there was no cross-reactivity with the other five common porcine viruses. With visual observation, the lower limits were approximately 104 TCID50/mL and 104 TCID50/mL for PEDV and TGEV, respectively. The GICA strip could be stored at 4°C or 25°C for 12 months without affecting its efficacy. To validate the GICA strip, 121 clinical samples were tested. The positive rates of PEDV and TGEV were 42.9 and 9.9%, respectively, and the co-infection rate of the two viruses was 5.8% based on the duplex GICA strip. Thus, the established GICA strip is a rapid, specific, and stable tool for on-site preliminary diagnosis of PEDV- and TGEV-associated diarrhea.

Establishment and application of a quadruplex real-time RT-qPCR assay for differentiation of TGEV, PEDV, PDCoV, and PoRVA

Porcine viral diarrhea is a common ailment in clinical settings, causing significant economic losses to the swine industry. Notable culprits behind porcine viral diarrhea encompass transmissible gastroenteritis virus (TGEV), porcine epidemic diarrhea virus (PEDV), porcine deltacoronavirus (PDCoV), and porcine rotavirus-A (PoRVA). Co-infections involving the viruses are a common occurrence in clinical settings, thereby amplifying the complexities associated with differential diagnosis. As a consequence, it is therefore necessary to develop a method that can detect and differentiate all four porcine diarrhea viruses (TGEV, PEDV, PDCoV, and PoRVA) with a high sensitivity and specificity. Presently, polymerase chain reaction (PCR) is the go-to method for pathogen detection. In comparison to conventional PCR, TaqMan real-time PCR offers heightened sensitivity, superior specificity, and enhanced accuracy. This study aimed to develop a quadruplex real-time RT-qPCR assay, utilizing TaqMan probes, for the distinctive detection of TGEV, PEDV, PDCoV, and PoRVA. The quadruplex real-time RT-qPCR assay, as devised in this study, exhibited the capacity to avoid the detection of unrelated pathogens and demonstrated commendable specificity, sensitivity, repeatability, and reproducibility, boasting a limit of detection (LOD) of 27 copies/μL. In a comparative analysis involving 5483 clinical samples, the results from the commercial RT-qPCR kit and the quadruplex RT-qPCR for TGEV, PEDV, PDCoV, and PoRVA detection were entirely consistent. Following sample collection from October to March in Guangxi Zhuang Autonomous Region, we assessed the prevalence of TGEV, PEDV, PDCoV, and PoRVA in piglet diarrhea samples, revealing positive detection rates of 0.2 % (11/5483), 8.82 % (485/5483), 1.22 % (67/5483), and 4.94 % (271/5483), respectively. The co-infection rates of PEDV/PoRVA, PEDV/PDCoV, TGEV/PED/PoRVA, and PDCoV/PoRVA were 0.39 %, 0.11 %, 0.01 %, and 0.03 %, respectively, with no detection of other co-infections, as determined by the quadruplex real-time RT-qPCR. This research not only established a valuable tool for the simultaneous differentiation of TGEV, PEDV, PDCoV, and PoRVA in practical applications but also provided crucial insights into the prevalence of these viral pathogens causing diarrhea in Guangxi.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

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.

Regulatory Non-Coding RNAs during Porcine Viral Infections: Potential Targets for Antiviral Therapy

Pigs play important roles in agriculture and bio-medicine; however, porcine viral infections have caused huge losses to the pig industry and severely affected the animal welfare and social public safety. During viral infections, many non-coding RNAs are induced or repressed by viruses and regulate viral infection. Many viruses have, therefore, developed a number of mechanisms that use ncRNAs to evade the host immune system. Understanding how ncRNAs regulate host immunity during porcine viral infections is critical for the development of antiviral therapies. In this review, we provide a summary of the classification, production and function of ncRNAs involved in regulating porcine viral infections. Additionally, we outline pathways and modes of action by which ncRNAs regulate viral infections and highlight the therapeutic potential of artificial microRNA. Our hope is that this information will aid in the development of antiviral therapies based on ncRNAs for the pig industry.

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.

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

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

Integrated metabolomics and transcriptomics analyses reveal metabolic responses to TGEV infection in porcine intestinal epithelial cells

Transmissible gastroenteritis virus (TGEV) is a coronavirus that infects piglets with severe diarrhoea, vomiting, dehydration, and even death, causing huge economic losses to the pig industry. The underlying pathogenesis of TGEV infection and the effects of TGEV infection on host metabolites remain poorly understood. To investigate the critical metabolites and regulatory factors during TGEV infection in intestinal porcine epithelial cells (IPEC-J2), we performed metabolomic and transcriptomic analyses of TGEV-infected IPEC-J2 cells by LC/MS and RNA-seq techniques. A total of 87 differential metabolites and 489 differentially expressed genes were detected. A series of metabolites and candidate genes from glutathione metabolism and AMPK signalling pathway were examined through combined analysis of metabolome and transcriptome. We found glutathione peroxidase 3 (GPX3) is markedly reduced after TGEV infection, and a significant negative correlation between AMPK signalling pathway and TGEV infection. Exogenous addition of the AMPK activator COH-SR4 significantly downregulates stearoyl coenzyme A (SCD1) mRNA and inhibits TGEV replication; while exogenous GSK-690693 significantly promotes TGEV infection by inhibiting AMPK signalling pathway. In summary, our study provides insights into the key metabolites and regulators for TGEV infection from the metabolome and transcriptome perspective, which will offer promising antiviral metabolic and molecular targets and enrich the understanding of the existence of a similar mechanism in the host.

Development of a Multiplex Quantitative PCR for Detecting Porcine Epidemic Diarrhea Virus, Transmissible Gastroenteritis Virus, and Porcine Deltacoronavirus Simultaneously in China

Porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis virus (TGEV), and porcine deltacoronavirus (PDCoV) belong to the category of swine enteric coronavirus that cause acute diarrhea in piglets, which has resulted in massive losses to the pig husbandry. Therefore, a sensitive and rapid detection method which can differentially detect these viruses that lead to mixed infections in clinical cases, is urgently needed. According to the conserved regions of the PEDV M gene, TGEV S gene, and PDCoV N gene, and the reference gene of porcine (β-Actin), we designed new specific primers and probes for the multiplex qPCR assay capable of simultaneously detecting three RNA viruses. This method, with a great specificity, did not cross-react with the common porcine virus. Moreover, the limit of detection of the method we developed could reach 10 copies/μL ,and the intra- and inter-group coefficients of variation of it below 3%. Applying this assay to detect 462 clinical samples which were collected in 2022-2023, indicated that the discrete positive rates of PEDV, TGEV, and PDCoV were 19.70%, 0.87%, and 10.17%, respectively. The mixed infection rates of PEDV/TGEV, PEDV/PDCoV, TGEV/PDCoV, and PEDV/TGEV/PDCoV were 3.25%, 23.16%, 0.22%, and 11.90%, respectively. All in all, the multiplex qPCR assay we developed as a tool for differential and rapid diagnosing can be put on the active prevention and control of PEDV, TGEV, and PDCoV, , which can create great value in the diagnosis of swine diarrhea diseases.

sgRNAs: A SARS-CoV-2 emerging issue

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

Transmissible Gastroenteritis Virus: An Update Review and Perspective

Transmissible gastroenteritis virus (TGEV) is a member of the alphacoronavirus genus, which has caused huge threats and losses to pig husbandry with a 100% mortality in infected piglets. TGEV is observed to be recombining and evolving unstoppably in recent years, with some of these recombinant strains spreading across species, which makes the detection and prevention of TGEV more complex. This paper reviews and discusses the basic biological properties of TGEV, factors affecting virulence, viral receptors, and the latest research advances in TGEV infection-induced apoptosis and autophagy to improve understanding of the current status of TGEV and related research processes. We also highlight a possible risk of TGEV being zoonotic, which could be evidenced by the detection of CCoV-HuPn-2018 in humans.

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.

Conserved recombination patterns across coronavirus subgenera

Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination break points at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination break points across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination break-point hotspot locations. We find that while the locations of recombination break points are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination break points most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination break-point distributions in coronavirus genomes sampled from nature.

RNA polymerase inaccuracy underlies SARS-CoV-2 variants and vaccine heterogeneity

Both the SARS-CoV-2 virus and its mRNA vaccines depend on RNA polymerases (RNAP)1,2; however, these enzymes are inherently error-prone and can introduce variants into the RNA3. To understand SARS-CoV-2 evolution and vaccine efficacy, it is critical to identify the extent and distribution of errors introduced by the RNAPs involved in each process. Current methods lack the sensitivity and specificity to measure de novo RNA variants in low input samples like viral isolates3. Here, we determine the frequency and nature of RNA errors in both SARS-CoV-2 and its vaccine using a targeted Accurate RNA Consensus sequencing method (tARC-seq). We found that the viral RNA-dependent RNAP (RdRp) makes ~1 error every 10,000 nucleotides - higher than previous estimates4. We also observed that RNA variants are not randomly distributed across the genome but are associated with certain genomic features and genes, such as S (Spike). tARC-seq captured a number of large insertions, deletions and complex mutations that can be modeled through non-programmed RdRp template switching. This template switching feature of RdRp explains many key genetic changes observed during the evolution of different lineages worldwide, including Omicron. Further sequencing of the Pfizer-BioNTech COVID-19 vaccine revealed an RNA variant frequency of ~1 in 5,000, meaning most of the vaccine transcripts produced in vitro by T7 phage RNAP harbor a variant. These results demonstrate the extraordinary genetic diversity of viral populations and the heterogeneous nature of an mRNA vaccine fueled by RNAP inaccuracy. Along with functional studies and pandemic data, tARC-seq variant spectra can inform models to predict how SARS-CoV-2 may evolve. Finally, our results may help improve future vaccine development and study design as mRNA therapies continue to gain traction.

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

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

Molecular characterization of the FCoV-like canine coronavirus HLJ-071 in China

Background

According to the differences of antigen and genetic composition, canine coronavirus (CCoV) consists of two genotypes, CCoV-I and CCoV-II. Since 2004, CCoVs with point mutations or deletions of NSPs are contributing to the changes in tropism and virulence in dogs.

Results

In this study, we isolated a CCoV, designated HLJ-071, from a dead 5-week-old female Welsh Corgi with severe diarrhea and vomit. Sequence analysis suggested that HLJ-071 bearing a complete ORF3abc compared with classic CCoV isolates (1-71, K378 and S378). In addition, a variable region was located between S gene and ORF 3a gene, in which a deletion with 104 nts for HLJ-071 when compared with classic CCoV strains 1-71, S378 and K378. Phylogenetic analysis based on the S gene and complete sequences showed that HLJ-071 was closely related to FCoV II. Recombination analysis suggested that HLJ-071 originated from the recombination of FCoV 79-1683, FCoV DF2 and CCoV A76. Finally, according to cell tropism experiments, it suggested that HLJ-071 could replicate in canine macrophages/monocytes cells.

Conclusion

The present study involved the isolation and genetic characterization of a variant CCoV strain and spike protein and ORF3abc of CCoV might play a key role in viral tropism, which could affect the replication in monocyte/macrophage cells. It will provide essential information for further understanding the evolution in China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-03073-8.

A brief molecular insight of COVID-19: epidemiology, clinical manifestation, molecular mechanism, cellular tropism and immuno-pathogenesis

In December 2019, the emergence and expansion of novel and infectious respiratory virus SARS-CoV-2 originated from Wuhan, China caused an unprecedented threat to the public health and became a global pandemic. SARS-CoV-2 is an enveloped, positive sense and single stranded RNA virus belonging to genera betacoronavirus, of Coronaviridae family. The viral genome sequencing studies revealed 75-80% similarity with SARS-CoV. SARS-CoV-2 mainly affects the lower respiratory system and may progress to pneumonia and Acute Respiratory Distress Syndrome (ARDS). Apart from life-threatening situations and burden on the global healthcare system, the COVID-19 pandemic has imposed several challenges on the worldwide economics and livelihood. The novel pathogen is highly virulent, rapidly mutating and has a tendency to cross the species boundaries such as from bats to humans through the evolution and natural selection from intermediate host. In this review we tried to summarize the overall picture of SARS-CoV-2 including origin/ emergence, epidemiology, pathogenesis, genome organization, comparative analysis with other CoVs, infection and replication mechanism along with cellular tropism and immunopathogenesis which will provide a brief panoramic view about the virus and disease.

CoV-er all the bases: Structural perspectives of SARS-CoV-2 RNA synthesis

The ongoing Covid-19 pandemic has spurred research in the biology of the nidovirus severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Much focus has been on the viral RNA synthesis machinery due to its fundamental role in viral propagation. The central and essential enzyme of the RNA synthesis process, the RNA-dependent RNA polymerase (RdRp), functions in conjunction with a coterie of viral-encoded enzymes that mediate crucial nucleic acid transactions. Some of these enzymes share common features with other RNA viruses, while others play roles unique to nidoviruses or CoVs. The RdRps are proven targets for viral pathogens, and many of the other nucleic acid processing enzymes are promising targets. The purpose of this review is to summarize recent advances in our understanding of the mechanisms of RNA synthesis in CoVs. By reflecting on these studies, we hope to emphasize the remaining gaps in our knowledge. The recent onslaught of structural information related to SARS-CoV-2 RNA synthesis, in combination with previous structural, genetic and biochemical studies, have vastly improved our understanding of how CoVs replicate and process their genomic RNA. Structural biology not only provides a blueprint for understanding the function of the enzymes and cofactors in molecular detail, but also provides a basis for drug design and optimization. The concerted efforts of researchers around the world, in combination with the renewed urgency toward understanding this deadly family of viruses, may eventually yield new and improved antivirals that provide relief to the current global devastation.

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.

An overview of SARS-COV-2 epidemiology, mutant variants, vaccines, and management strategies

BACKGROUND

Over the last two decades, humanity has observed the extraordinary anomaly caused by novel, weird coronavirus strains, such as severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS). As the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) virus has made its entry into the world, it has dramatically affected life in every domain by continuously producing new variants. The vaccine development is an ongoing process, although some vaccines got marketed. The big challenge is now whether the vaccine candidates can provide long-lasting protection or prevention against mutant variants.

METHODS

The information was gathered from various journals, electronic searches via Internet-based information such as PubMed, Google Scholar, Science Direct, online electronic journals, WHO landscape, world meters, WHO website, and News.

RESULTS

This review will present and discuss some coronavirus disease 19 (COVID-19) related aspects including: the pathophysiology, epidemiology, mutant variants vaccine candidates, vaccine efficacy, and management strategies. Due to the high death rate, continuous spread, an inadequate workforce, lack of required therapeutics, and incomplete understanding of the viral strain, it becomes crucial to build the knowledge of its biological characteristics and make available the rapid diagnostic and vital therapeutic machinery for the combat and management of an infection.

CONCLUSION

The data summarizes current research on the COVID 19 infection and therapeutic interventions, which will direct future decision-making on the effort-worthy phases of the COVID 19 and the development of critical therapeutics. The only possible solution is the vaccine development targeting against all variant strains to halt its progress; the identified theoretical and practical knowledge can eliminate the gaps to improve a better understanding of the novel coronavirus structure and its design of a vaccine. In addition, to that the long-lasting protection is another challenging objective that need to be looked into.

The structure of a dimeric form of SARS-CoV-2 polymerase

The coronavirus SARS-CoV-2 uses an RNA-dependent RNA polymerase (RdRp) to replicate and transcribe its genome. Previous structures of the RdRp revealed a monomeric enzyme composed of the catalytic subunit nsp12, two copies of subunit nsp8, and one copy of subunit nsp7. Here we report an alternative, dimeric form of the enzyme and resolve its structure at 5.5 Å resolution. In this structure, the two RdRps contain only one copy of nsp8 each and dimerize via their nsp7 subunits to adopt an antiparallel arrangement. We speculate that the RdRp dimer facilitates template switching during production of sub-genomic RNAs.

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.

Porcine enteric coronaviruses: an updated overview of the pathogenesis, prevalence, and diagnosis

The recent prevalence of coronavirus (CoV) poses a serious threat to animal and human health. Currently, porcine enteric coronaviruses (PECs), including the transmissible gastroenteritis virus (TGEV), the novel emerging swine acute diarrhoea syndrome coronavirus (SADS-CoV), porcine delta coronavirus (PDCoV), and re-emerging porcine epidemic diarrhoea virus (PEDV), which infect pigs of different ages, have caused more frequent occurrences of diarrhoea, vomiting, and dehydration with high morbidity and mortality in piglets. PECs have the potential for cross-species transmission and are causing huge economic losses in the pig industry in China and the world, which therefore needs to be urgently addressed. Accordingly, this article summarises the pathogenicity, prevalence, and diagnostic methods of PECs and provides an important reference for their improved diagnosis, prevention, and control.

The Recombination Potential between SARS-CoV-2 and MERS-CoV from Cross-Species Spill-over Infections

Countries in the Middle-East (ME) are tackling two corona virus outbreaks simultaneously, Middle-Eastern Respiratory Syndrome Coronavirus (MERS-CoV) and the current Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Both viruses infect the same host (humans) and the same cell (type-II alveolar cells) causing lower respiratory illnesses such as pneumonia. Molecularly, MERS-CoV and SARS-CoV-2 enter alveolar cells via spike proteins recognizing dipeptidyl peptidase-4 and angiotensin converting enzyme-II, respectively. Intracellularly, both viruses hide in organelles to generate negative RNA strands and initiate replication using very similar mechanisms. At the transcription level, both viruses utilise identical Transcription Regulatory Sequences (TRSs), which are known recombination cross-over points during replication, to transcribe genes. Using whole genome alignments of both viruses, we identify clusters of high sequence homology at ORF1a and ORF1b. Given the high recombination rates detected in SARS-CoV-2, we speculate that in co-infections recombination is feasible via TRS and/or clusters of homologies. Accordingly, here we recommend mitigation measure and testing for both MERS-CoV and SARS-CoV-2 in ME countries.

Transcriptome & viral growth analysis of SARS-CoV-2-infected Vero CCL-81 cells

Background & objectives:

The genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), belonging to the family Coronaviridae, encodes for structural, non-structural, and accessory proteins, which are required for replication of the virus. These proteins are encoded by different genes present on the SARS-CoV-2 genome. The expression pattern of these genes in the host cells needs to be assessed. This study was undertaken to understand the transcription pattern of the SARS-CoV-2 genes in the Vero CCL-81 cells during the course of infection.

Methods:

Vero CCL-81 cells were infected with the SARS-CoV-2 virus inoculum having a 0.1 multiplicity of infection. The supernatants and cell pellets were harvested after centrifugation at different time points, post-infection. The 50% tissue culture infective dose (TCID₅₀) and cycle threshold (C_t) values of the E and the RdRp-2 genes were calculated. Next-generation sequencing of the harvested sample was carried out to observe the expression pattern of the virus by mapping to the SARS-CoV-2 Wuhan HU-1 reference sequence. The expressions were in terms of the reads per kilobase million (RPKM) values.

Results:

In the inital six hours post-infection, the copy numbers of E and RdRp-2 genes were approximately constant, which raised 10 log-fold and continued to increase till the 12 h post-infection (hpi). The TCID₅₀ was observed in the supernatant after 7 hpi, indicating the release of the viral progeny. ORF8 and ORF7a, along with the nucleocapsid transcript, were found to express at higher levels.

Interpretation & conclusions:

This study was a step towards understanding the growth kinetics of the SARS-CoV-2 replication cycle. The findings indicated that ORF8 and ORF7b gene transcripts were expressed in higher amounts indicating their essential role in viral replication. Future studies need to be conducted to explore their role in the SARS-CoV-2 replication.

Structural Basis for Helicase-Polymerase Coupling in the SARS-CoV-2 Replication-Transcription Complex

SARS-CoV-2 is the causative agent of the 2019-2020 pandemic. The SARS-CoV-2 genome is replicated and transcribed by the RNA-dependent RNA polymerase holoenzyme (subunits nsp7/nsp8₂/nsp12) along with a cast of accessory factors. One of these factors is the nsp13 helicase. Both the holo-RdRp and nsp13 are essential for viral replication and are targets for treating the disease COVID-19. Here we present cryoelectron microscopic structures of the SARS-CoV-2 holo-RdRp with an RNA template product in complex with two molecules of the nsp13 helicase. The Nidovirales order-specific N-terminal domains of each nsp13 interact with the N-terminal extension of each copy of nsp8. One nsp13 also contacts the nsp12 thumb. The structure places the nucleic acid-binding ATPase domains of the helicase directly in front of the replicating-transcribing holo-RdRp, constraining models for nsp13 function. We also observe ADP-Mg²⁺ bound in the nsp12 N-terminal nidovirus RdRp-associated nucleotidyltransferase domain, detailing a new pocket for anti-viral therapy development.

Structural basis for helicase-polymerase coupling in the SARS-CoV-2 replication-transcription complex

SARS-CoV-2 is the causative agent of the 2019-2020 pandemic. The SARS-CoV-2 genome is replicated-transcribed by the RNA-dependent RNA polymerase holoenzyme (subunits nsp7/nsp82/nsp12) along with a cast of accessory factors. One of these factors is the nsp13 helicase. Both the holo-RdRp and nsp13 are essential for viral replication and are targets for treating the disease COVID-19. Here we present cryo-electron microscopic structures of the SARS-CoV-2 holo-RdRp with an RNA template-product in complex with two molecules of the nsp13 helicase. The Nidovirus-order-specific N-terminal domains of each nsp13 interact with the N-terminal extension of each copy of nsp8. One nsp13 also contacts the nsp12-thumb. The structure places the nucleic acid-binding ATPase domains of the helicase directly in front of the replicating-transcribing holo-RdRp, constraining models for nsp13 function. We also observe ADP-Mg2+ bound in the nsp12 N-terminal nidovirus RdRp-associated nucleotidyltransferase domain, detailing a new pocket for anti-viral therapeutic development.

Evidence for Internal Initiation of RNA Synthesis by the Hepatitis C Virus RNA-Dependent RNA Polymerase NS5B In Cellulo

Initiation of RNA synthesis by the hepatitis C virus (HCV) RNA-dependent RNA polymerase (RdRp) NS5B has been extensively studied in vitro and in cellulo Intracellular replication is thought to rely exclusively on terminal de novo initiation, as it conserves all genetic information of the genome. In vitro, however, additional modes of initiation have been observed. In this study, we aimed to clarify whether the intracellular environment allows for internal initiation of RNA replication by the HCV replicase. We used a dual luciferase replicon harboring a terminal and an internal copy of the viral genomic 5' untranslated region, which was anticipated to support noncanonical initiation. Indeed, a shorter RNA species was detected by Northern blotting with low frequency, depending on the length and sequence composition upstream of the internal initiation site. By introducing mutations at either site, we furthermore established that internal and terminal initiation shared identical sequence requirements. Importantly, lethal point mutations at the terminal site resulted exclusively in truncated replicons. In contrast, the same mutations at the internal site abrogated internal initiation, suggesting a competitive selection of initiation sites, rather than recombination or template-switching events. In conclusion, our data indicate that the HCV replicase is capable of internal initiation in its natural environment, although functional replication likely requires only terminal initiation. Since many other positive-strand RNA viruses generate subgenomic messenger RNAs during their replication cycle, we surmise that their capability for internal initiation is a common and conserved feature of viral RdRps.IMPORTANCE Many aspects of viral RNA replication of hepatitis C virus (HCV) are still poorly understood. The process of RNA synthesis is driven by the RNA-dependent RNA polymerase (RdRp) NS5B. Most mechanistic studies on NS5B so far were performed with in vitro systems using isolated recombinant polymerase. In this study, we present a replicon model, which allows the intracellular assessment of noncanonical modes of initiation by the full HCV replicase. Our results add to the understanding of the biochemical processes underlying initiation of RNA synthesis by NS5B by the discovery of internal initiation in cellulo Moreover, they validate observations made in vitro, showing that the viral polymerase acts very similarly in isolation and in complex with other viral and host proteins. Finally, these observations provide clues about the evolution of RdRps of positive-strand RNA viruses, which might contain the intrinsic ability to initiate internally.

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.

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.

Characterization of Chinese Porcine Epidemic Diarrhea Virus with Novel Insertions and Deletions in Genome

Outbreaks of porcine epidemic diarrhoea virus (PEDV) have caused great economic losses to the global pig industry. PEDV strains with variants in the spike (S) gene have been reported in several countries. To better understand the molecular epidemiology and genetic diversity of PEDV field isolates, in this study, we characterised the complete genome sequence of a novel PEDV variant JSCZ1601 from a outbreak in China in 2016. The PEDV isolate was 28,033 nucleotides (nt) in length without the polyadenylated sequences. Phylogenetic analysis based on the full-length genome sequence of JSCZ1601 grouped it with the pandemic variants determined post-2010 into group 2 (G2). However, the S gene of JSCZ1601 formed a new subgroup separated from the subgroups containing the other G2 strains. Comparative analysis of the amino acids encoded by the S genes revealed the N-terminal of the deduced JSCZ1601 S protein had a novel two-amino-acid deletion (N58 and S59) compared with all identified genogroups. Further, compared with the reference strains, a 'G' insertion was detected in the 5' terminal of the 5'UTR of the JSCZ1601. The animal experiment revealed that this strain was high pathogenic to neonatal pigs. Taken together, a PEDV strain with the new molecular characterizations and phylogenies was found in mainland China. It is necessary to strengthen the monitoring of PEDV variations.

Immune evasion of porcine enteric coronaviruses and viral modulation of antiviral innate signaling

•

Enteric coronaviruses have evolved to modulate the host innate immunity.

•

Viral IFN antagonists have been identified and they are mostly redundant.

•

For protection of intestinal epithelia from enteric viruses, type III IFN plays a major role.

Porcine epidemic diarrhea virus (PEDV) and porcine deltacoronavirus (PDCoV) are emerged and reemerging viruses in pigs, and together with transmissible gastroenteritis virus (TGEV), pose significant economic concerns to the swine industry. These viruses infect epithelial cells of the small intestine and cause watery diarrhea, dehydration, and a high mortality in neonatal piglets. Type I interferons (IFN-α/β) are major antiviral cytokines forming host innate immunity, and in turn, these enteric coronaviruses have evolved to modulate the host innate immune signaling during infection. Accumulating evidence however suggests that IFN induction and signaling in the intestinal epithelial cells differ from other epithelial cells, largely due to distinct features of the gut epithelial mucosal surface and commensal microflora, and it appears that type III interferon (IFN-λ) plays a key role to maintain the antiviral state in the gut. This review describes the recent understanding on the immune evasion strategies of porcine enteric coronaviruses and the role of different types of IFNs for intestinal antiviral innate immunity.

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.

Antigenic structures stably expressed by recombinant TGEV-derived vectors

Coronaviruses (CoVs) are positive-stranded RNA viruses with potential as immunization vectors, expressing high levels of heterologous genes and eliciting both secretory and systemic immune responses. Nevertheless, its high recombination rate may result in the loss of the full-length foreign gene, limiting their use as vectors. Transmissible gastroenteritis virus (TGEV) was engineered to express porcine reproductive and respiratory syndrome virus (PRRSV) small protein domains, as a strategy to improve heterologous gene stability. After serial passage in tissue cultures, stable expression of small PRRSV protein antigenic domains was achieved. Therefore, size reduction of the heterologous genes inserted in CoV-derived vectors led to the stable expression of antigenic domains. Immunization of piglets with these TGEV vectors led to partial protection against a challenge with a virulent PRRSV strain, as immunized animals showed reduced clinical signs and lung damage. Further improvement of TGEV-derived vectors will require the engineering of vectors with decreased recombination rate.

Infectious Bronchitis Virus as a Vector for the Expression of Heterologous Genes

The avian coronavirus infectious bronchitis virus (IBV) is the causative agent of the respiratory disease infectious bronchitis of domestic fowl, and is controlled by routine vaccination. To explore the potential use of IBV as a vaccine vector a reverse genetics system was utilised to generate infectious recombinant IBVs (rIBVs) expressing the reporter genes enhanced green fluorescent protein (eGFP) or humanised Renilla luciferase (hRluc). Infectious rIBVs were obtained following the replacement of Gene 5 or the intergenic region (IR) with eGFP or hRluc, or the replacement of ORFs 3a and 3b with hRluc. The replacement of Gene 5 with an IBV codon-optimised version of the hRluc gene also resulted in successful rescue of infectious rIBV. Reporter gene expression was confirmed by fluorescence microscopy, or luciferase activity assays, for all successfully rescued rIBVs following infection of primary chick kidney (CK) cells. The genetic stability of rIBVs was analysed by serial passage on CK cells. Recombinant IBV stability varied depending on the genome region being replaced, with the reporter genes maintained up to at least passage 8 (P8) following replacement of Gene 5, P7 for replacement of the IR and P5 for replacement of ORFs 3a and 3b. Codon-optimisation of the hRluc gene, when replacing Gene 5, resulted in an increase in genome stability, with hRluc expression stable up to P10 compared to P8 for standard hRluc. Repeated passaging of rIBVs expressing hRluc at an MOI of 0.01 demonstrated an increase in stability, with hRluc expression stable up to at least P12 following the replacement of Gene 5. This study has demonstrated that heterologous genes can be incorporated into, and expressed from a range of IBV genome locations and that replacement of accessory Gene 5 offers a promising target for realising the potential of IBV as a vaccine vector for other avian pathogens.

Progress in understanding pathogenic mechanisms of porcine transmissible gastroenteritis virus

Porcine transmissible gastroenteritis virus (TG-EV) is an animal coronavirus that causes severe gastroenteritis in young TGEV-seronegative pigs. This review will focus on recent advances in research of the genomic structure, major structural proteins and their function, virus propagation and replication, virus receptors, genetics and pathogenic mechanisms of TGEV. These data will be helpful in understanding the molecular biological characteristics and genetic variation of TGEV and have important theoretical significance for the development of new vaccines and antiviral drugs.

Complete genome of transmissible gastroenteritis virus AYU strain isolated in Shanghai, China

Transmissible gastroenteritis virus strain AYU was isolated in Shanghai. The complete genome has a length of 28,582 bp and contains seven open reading frames. Sequence analysis suggested that Shanghai strain AYU and U.S. strain Purdue P115 are derived from a common ancestor, as they have 99.6% similarity at the nucleotide level.

Long-distance RNA-RNA interactions in the coronavirus genome form high-order structures promoting discontinuous RNA synthesis during transcription

Coronavirus (CoV) transcription requires a high-frequency recombination process that links newly synthesized minus-strand subgenomic RNA copies to the leader region, which is present only once, at the 5' end of the genome. This discontinuous RNA synthesis step is based on the complementarity between the transcription-regulating sequences (TRSs) at the leader region and those preceding each gene in the nascent minus-strand RNA. Furthermore, the template switch requires the physical proximity of RNA genome domains located between 20,000 and 30,000 nucleotides apart. In this report, it is shown that the efficacy of this recombination step is promoted by novel additional long-distance RNA-RNA interactions between RNA motifs located close to the TRSs controlling the expression of each gene and their complementary sequences mapping close to the 5' end of the genome. These interactions would bring together the motifs involved in the recombination process. This finding indicates that the formation of high-order RNA structures in the CoV genome is necessary to control the expression of at least the viral N gene. The requirement of these long-distance interactions for transcription was shown by the engineering of CoV replicons in which the complementarity between the newly identified sequences was disrupted. Furthermore, disruption of complementarity in mutant viruses led to mutations that restored complementarity, wild-type transcription levels, and viral titers by passage in cell cultures. The relevance of these high-order structures for virus transcription is reinforced by the phylogenetic conservation of the involved RNA motifs in CoVs.

In vitro inhibition of transmissible gastroenteritis coronavirus replication in swine testicular cells by short hairpin RNAs targeting the ORF 7 gene

Background

Transmissible gastroenteritis (TGE) is a highly contagious viral disease of swine, characterized by severe vomiting, diarrhea, and high mortality. Currently, the vaccines for it are only partially effective and no specific drug is available for treatment of TGE virus (TGEV) infection. RNA interference has been confirmed as a new approach for controlling viral infections. In this study, the inhibitory effect of short hairpin RNAs (shRNAs) targeting the ORF 7 gene of TGEV on virus replication was examined.

Results

Four theoretically effective sequences of TGEV ORF 7 gene were designed and selected for construction of shRNA expression plasmids. In the reporter assays, three of four shRNA expression plasmids were able to inhibit significantly the expression of ORF 7 gene and replication of TGEV, as shown by real-time quantitative RT-PCR analysis of viral ORF 7 and N genes and detection of virus titers (TCID₅₀/ml). Stable swine testicular (ST) cells expressing the shRNAs were established. Observation of the cytopathic effect and apoptosis, as well as a cell proliferation assay demonstrated that the three shRNAs were capable of protecting ST cells against TGEV destruction, with high specificity and efficiency.

Conclusions

Our results indicated that plasmid-transcribed shRNAs targeting the ORF 7 gene in the TGEV genome effectively inhibited expression of the viral target gene and viral replication in vitro. These findings provide evidence that the shRNAs have potential therapeutic application for treatment of TGE.

Functional transcriptional regulatory sequence (TRS) RNA binding and helix destabilizing determinants of murine hepatitis virus (MHV) nucleocapsid (N) protein

Coronavirus (CoV) nucleocapsid (N) protein contains two structurally independent RNA binding domains. These are denoted N-terminal domain (NTD) and C-terminal domain and are joined by a charged linker region rich in serine and arginine residues (SR linker). In mouse hepatitis virus (MHV), the NTD binds the transcriptional regulatory sequence (TRS) RNA, a conserved hexanucleotide sequence required for subgenomic RNA synthesis. The NTD is also capable of disrupting a short RNA duplex. We show here that three residues on the β3 (Arg-125 and Tyr-127) and β5 (Tyr-190) strands play key roles in TRS RNA binding and helix destabilization with Ala substitutions of these residues lethal to the virus. NMR studies of the MHV NTD·TRS complex revealed that this region defines a major RNA binding interface in MHV with site-directed spin labeling studies consistent with a model in which the adenosine-rich 3'-region of TRS is anchored by Arg-125, Tyr-127, and Tyr-190 in a way that is critical for efficient subgenomic RNA synthesis in MHV. Characterization of CoV N NTDs from infectious bronchitis virus and from severe acute respiratory syndrome CoV revealed that, although detailed NTD-TRS determinants are distinct from those of MHV NTD, rapid helix destabilization activity of CoV N NTDs is most strongly correlated with CoV function and virus viability.

Subgenomic messenger RNAs: mastering regulation of (+)-strand RNA virus life cycle

Many (+)-strand RNA viruses use subgenomic (SG) RNAs as messengers for protein expression, or to regulate their viral life cycle. Three different mechanisms have been described for the synthesis of SG RNAs. The first mechanism involves internal initiation on a (−)-strand RNA template and requires an internal SGP promoter. The second mechanism makes a prematurely terminated (−)-strand RNA which is used as template to make the SG RNA. The third mechanism uses discontinuous RNA synthesis while making the (−)-strand RNA templates. Most SG RNAs are translated into structural proteins or proteins related to pathogenesis: however other SG RNAs regulate the transition between translation and replication, function as riboregulators of replication or translation, or support RNA–RNA recombination. In this review we discuss these functions of SG RNAs and how they influence viral replication, translation and recombination.

RNA-RNA and RNA-protein interactions in coronavirus replication and transcription

Coronavirus (CoV) RNA synthesis includes the replication of the viral genome, and the transcription of sgRNAs by a discontinuous mechanism. Both processes are regulated by RNA sequences such as the 5' and 3' untranslated regions (UTRs), and the transcription regulating sequences (TRSs) of the leader (TRS-L) and those preceding each gene (TRS-Bs). These distant RNA regulatory sequences interact with each other directly and probably through protein-RNA and protein-protein interactions involving viral and cellular proteins. By analogy to other plus-stranded RNA viruses, such as polioviruses, in which translation and replication switch involves a cellular factor (PCBP) and a viral protein (3CD) it is conceivable that in CoVs the switch between replication and transcription is also associated with the binding of proteins that are specifically recruited by the replication or transcription complexes. Complexes between RNA motifs such as TRS-L and the TRS-Bs located along the CoV genome are probably formed previously to the transcription start, and most likely promote template-switch of the nascent minus RNA to the TRS-L region. Many cellular proteins interacting with regulatory CoV RNA sequences are members of the heterogeneous nuclear ribonucleoprotein (hnRNP) family of RNA-binding proteins, involved in mRNA processing and transport, which shuttle between the nucleus and the cytoplasm. In the context of CoV RNA synthesis, these cellular ribonucleoproteins might also participate in RNA-protein complexes to bring into physical proximity TRS-L and distant TRS-B, as proposed for CoV discontinuous transcription. In this review, we summarize RNA-RNA and RNA-protein interactions that represent modest examples of complex quaternary RNA-protein structures required for the fine-tuning of virus replication. Design of chemically defined replication and transcription systems will help to clarify the nature and activity of these structures.