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

Structural analysis of SARS-CoV-2 genome and predictions of the human interactome

Specific elements of viral genomes regulate interactions within host cells. Here, we calculated the secondary structure content of >2000 coronaviruses and computed >100 000 human protein interactions with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genomic regions display different degrees of conservation. SARS-CoV-2 domain encompassing nucleotides 22 500-23 000 is conserved both at the sequence and structural level. The regions upstream and downstream, however, vary significantly. This part of the viral sequence codes for the Spike S protein that interacts with the human receptor angiotensin-converting enzyme 2 (ACE2). Thus, variability of Spike S is connected to different levels of viral entry in human cells within the population. Our predictions indicate that the 5' end of SARS-CoV-2 is highly structured and interacts with several human proteins. The binding proteins are involved in viral RNA processing, include double-stranded RNA specific editases and ATP-dependent RNA-helicases and have strong propensity to form stress granules and phase-separated assemblies. We propose that these proteins, also implicated in viral infections such as HIV, are selectively recruited by SARS-CoV-2 genome to alter transcriptional and post-transcriptional regulation of host cells and to promote viral replication.

Structural analysis of SARS-CoV-2 genome and predictions of the human interactome

Specific elements of viral genomes regulate interactions within host cells. Here, we calculated the secondary structure content of >2000 coronaviruses and computed >100 000 human protein interactions with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genomic regions display different degrees of conservation. SARS-CoV-2 domain encompassing nucleotides 22 500-23 000 is conserved both at the sequence and structural level. The regions upstream and downstream, however, vary significantly. This part of the viral sequence codes for the Spike S protein that interacts with the human receptor angiotensin-converting enzyme 2 (ACE2). Thus, variability of Spike S is connected to different levels of viral entry in human cells within the population. Our predictions indicate that the 5' end of SARS-CoV-2 is highly structured and interacts with several human proteins. The binding proteins are involved in viral RNA processing, include double-stranded RNA specific editases and ATP-dependent RNA-helicases and have strong propensity to form stress granules and phase-separated assemblies. We propose that these proteins, also implicated in viral infections such as HIV, are selectively recruited by SARS-CoV-2 genome to alter transcriptional and post-transcriptional regulation of host cells and to promote viral replication.

Coronavirus: proteomics analysis of chicken kidney tissue infected with variant 2 (IS-1494)-like avian infectious bronchitis virus

Avian infectious bronchitis virus is one of the most important gammacoronaviruses, which causes a highly contagious disease. In this study, we investigated changes in the proteome of kidney tissue of specific-pathogen-free (SPF) chickens that were infected with an isolate of the nephrotropic variant 2 genotype (IS/1494/06) of avian coronavirus. Twenty 1-day-old SPF White Leghorn chickens were randomly divided into two groups, each comprising 10 chickens, which were kept in separate positive-pressure isolators. Chickens in group A served as a virus-free control group up to the end of the experiment, whereas chickens in group B were inoculated with 0.1 ml of 10^4.5 EID₅₀ of the IBV/chicken/Iran/UTIVO-C/2014 isolate of IBV, and kidney tissue samples were collected at 2 and 7 days post-inoculation (dpi) from both groups. Sequencing of five protein spots at 2 dpi and 22 spots at 7 dpi that showed differential expression by two-dimensional electrophoresis (2DE) along with fold change greater than 2 was done by MS-MALDI/TOF/TOF. Furthermore, the corresponding protein-protein interaction (PPI) networks at 2 and 7 dpi were identified to develop a detailed understanding of the mechanism of molecular pathogenesis. Topological graph analysis of this undirected PPI network revealed the effect of 10 genes in the 2 dpi PPI network and nine genes in the 7 dpi PPI network during virus pathogenesis. Proteins that were found by 2DE analysis and MS/TOF-TOF mass spectrometry to be down- or upregulated were subjected to PPI network analysis to identify interactions with other cellular components. The results show that cellular metabolism was altered due to viral infection. Additionally, multifunctional heat shock proteins with a significant role in host cell survival may be employed circuitously by the virus to reach its target. The data from this study suggest that the process of pathogenesis that occurs during avian coronavirus infection involves the regulation of vital cellular processes and the gradual disruption of critical cellular functions.

Role of SARS-CoV-2 in Altering the RNA-Binding Protein and miRNA-Directed Post-Transcriptional Regulatory Networks in Humans

The outbreak of a novel coronavirus SARS-CoV-2 responsible for the COVID-19 pandemic has caused a worldwide public health emergency. Due to the constantly evolving nature of the coronaviruses, SARS-CoV-2-mediated alterations on post-transcriptional gene regulations across human tissues remain elusive. In this study, we analyzed publicly available genomic datasets to systematically dissect the crosstalk and dysregulation of the human post-transcriptional regulatory networks governed by RNA-binding proteins (RBPs) and micro-RNAs (miRs) due to SARS-CoV-2 infection. We uncovered that 13 out of 29 SARS-CoV-2-encoded proteins directly interacted with 51 human RBPs, of which the majority of them were abundantly expressed in gonadal tissues and immune cells. We further performed a functional analysis of differentially expressed genes in mock-treated versus SARS-CoV-2-infected lung cells that revealed enrichment for the immune response, cytokine-mediated signaling, and metabolism-associated genes. This study also characterized the alternative splicing events in SARS-CoV-2-infected cells compared to the control, demonstrating that skipped exons and mutually exclusive exons were the most abundant events that potentially contributed to differential outcomes in response to the viral infection. A motif enrichment analysis on the RNA genomic sequence of SARS-CoV-2 clearly revealed the enrichment for RBPs such as SRSFs, PCBPs, ELAVs, and HNRNPs, suggesting the sponging of RBPs by the SARS-CoV-2 genome. A similar analysis to study the interactions of miRs with SARS-CoV-2 revealed functionally important miRs that were highly expressed in immune cells, suggesting that these interactions may contribute to the progression of the viral infection and modulate the host immune response across other human tissues. Given the need to understand the interactions of SARS-CoV-2 with key post-transcriptional regulators in the human genome, this study provided a systematic computational analysis to dissect the role of dysregulated post-transcriptional regulatory networks controlled by RBPs and miRs across tissue types during a SARS-CoV-2 infection.

Role of SARS-CoV-2 in altering the RNA binding protein and miRNA directed post-transcriptional regulatory networks in humans

The outbreak of a novel coronavirus SARS-CoV-2 responsible for COVID-19 pandemic has caused worldwide public health emergency. Due to the constantly evolving nature of the coronaviruses, SARS-CoV-2 mediated alteration on post-transcriptional gene regulation across human tissues remains elusive. In this study, we analyze publicly available genomic datasets to systematically dissect the crosstalk and dysregulation of human post-transcriptional regulatory networks governed by RNA binding proteins (RBPs) and micro-RNAs (miRs), due to SARS-CoV-2 infection. We uncovered that 13 out of 29 SARS-CoV-2 encoded proteins directly interact with 51 human RBPs of which majority of them were abundantly expressed in gonadal tissues and immune cells. We further performed a functional analysis of differentially expressed genes in mock-treated versus SARS-CoV-2 infected lung cells that revealed enrichment for immune response, cytokine-mediated signaling, and metabolism associated genes. This study also characterized the alternative splicing events in SARS-CoV-2 infected cells compared to control demonstrating that skipped exons and mutually exclusive exons were the most abundant events that potentially contributed to differential outcomes in response to viral infection. Motif enrichment analysis on the RNA genomic sequence of SARS-CoV-2 clearly revealed the enrichment for RBPs such as SRSFs, PCBPs, ELAVs, and HNRNPs suggesting the sponging of RBPs by SARS-CoV-2 genome. A similar analysis to study the interactions of miRs with SARS-CoV-2 revealed functionally important miRs that were highly expressed in immune cells, suggesting that these interactions may contribute to the progression of the viral infection and modulate host immune response across other human tissues. Given the need to understand the interactions of SARS-CoV-2 with key post-transcriptional regulators in the human genome, this study provides a systematic computational analysis to dissect the role of dysregulated post-transcriptional regulatory networks controlled by RBPs and miRs, across tissues types during SARS-CoV-2 infection.

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.

Coronavirus RNA Proofreading: Molecular Basis and Therapeutic Targeting

The coronavirus disease 2019 (COVID-19) that is wreaking havoc on worldwide public health and economies has heightened awareness about the lack of effective antiviral treatments for human coronaviruses (CoVs). Many current antivirals, notably nucleoside analogs (NAs), exert their effect by incorporation into viral genomes and subsequent disruption of viral replication and fidelity. The development of anti-CoV drugs has long been hindered by the capacity of CoVs to proofread and remove mismatched nucleotides during genome replication and transcription. Here, we review the molecular basis of the CoV proofreading complex and evaluate its potential as a drug target. We also consider existing nucleoside analogs and novel genomic techniques as potential anti-CoV therapeutics that could be used individually or in combination to target the proofreading mechanism.

A comprehensive review about SARS-CoV-2

The coronavirus disease (COVID-19) was first identified in China, December 2019. Since then, it has spread the length and breadth of the world at an unprecedented, alarming rate. Severe acute respiratory syndrome coronavirus (SARS-CoV)-2, which causes COVID-19, has much in common with its closest homologs, SARS-CoV and Middle East respiratory syndrome-CoV. The virus-host interaction of SARS-CoV-2 uses the same receptor, ACE2, which is similar to that of SARS-CoV, which spreads through the respiratory tract. Patients with COVID-19 report symptoms including mild-to-severe fever, cough and fatigue; very few patients report gastrointestinal infections. There are no specific antiviral strategies. A few strong medications are under investigation, so we have to focus on proposals which ought to be taken to forestall this infection in a living host.

The interface between coronaviruses and host cell RNA biology: Novel potential insights for future therapeutic intervention

Coronaviruses, including SARS-Cov-2, are RNA-based pathogens that interface with a large variety of RNA-related cellular processes during infection. These processes include capping, polyadenylation, localization, RNA stability, translation, and regulation by RNA binding proteins or noncoding RNA effectors. The goal of this article is to provide an in-depth perspective on the current state of knowledge of how various coronaviruses interact with, usurp, and/or avoid aspects of these cellular RNA biology machineries. A thorough understanding of how coronaviruses interact with RNA-related posttranscriptional processes in the cell should allow for new insights into aspects of viral pathogenesis as well as identify new potential avenues for the development of anti-coronaviral therapeutics. This article is categorized under: RNA in Disease and Development > RNA in Disease.

RNA-mediated translation regulation in viral genomes: computational advances in the recognition of sequences and structures

RNA structures are widely distributed across all life forms. The global conformation of these structures is defined by a variety of constituent structural units such as helices, hairpin loops, kissing-loop motifs and pseudoknots, which often behave in a modular way. Their ubiquitous distribution is associated with a variety of functions in biological processes. The location of these structures in the genomes of RNA viruses is often coordinated with specific processes in the viral life cycle, where the presence of the structure acts as a checkpoint for deciding the eventual fate of the process. These structures have been found to adopt complex conformations and exert their effects by interacting with ribosomes, multiple host translation factors and small RNA molecules like miRNA. A number of such RNA structures have also been shown to regulate translation in viruses at the level of initiation, elongation or termination. The role of various computational studies in the preliminary identification of such sequences and/or structures and subsequent functional analysis has not been fully appreciated. This review aims to summarize the processes in which viral RNA structures have been found to play an active role in translational regulation, their global conformational features and the bioinformatics/computational tools available for the identification and prediction of these structures.

In silico identification of conserved cis-acting RNA elements in the SARS-CoV-2 genome

Aim:

The aim of this study was to computationally predict conserved RNA sequences and structures known as cis-acting RNA elements (CREs) in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome.

Materials & methods:

Bioinformatics tools were used to analyze and predict CREs by obtaining viral sequences from available databases.

Results:

Computational analysis revealed the presence of RNA stem-loop structures within the 3′ end of the ORF1ab region analogous to previously identified SARS-CoV genomic packaging signals. Alignment-based RNA secondary structure predictions of the 5′ end of the SARS-CoV-2 genome also identified conserved CREs.

Conclusion:

These CREs may be potential vaccine and/or antiviral therapeutic targets; however, further studies are warranted to confirm their roles in the SARS-CoV-2 life cycle.

Establishment of a Virulent Full-Length cDNA Clone for Type I Feline Coronavirus Strain C3663

Feline infectious peritonitis (FIP) is one of the most important infectious diseases in cats and is caused by feline coronavirus (FCoV). Tissue culture-adapted type I FCoV shows reduced FIP induction in experimental infections, which complicates the understanding of FIP pathogenesis caused by type I FCoV. We previously found that the type I FCoV strain C3663 efficiently induces FIP in specific-pathogen-free cats through the naturally infectious route. In this study, we employed a bacterial artificial chromosome-based reverse genetics system to gain more insights into FIP caused by the C3633 strain. We successfully generated recombinant virus (rC3663) from Fcwf-4 cells transfected with infectious cDNA that showed growth kinetics similar to those shown by the parental virus. Next, we constructed a reporter C3663 virus carrying the nanoluciferase (Nluc) gene to measure viral replication with high sensitivity. The inhibitory effects of different compounds against rC3663-Nluc could be measured within 24 h postinfection. Furthermore, we found that A72 cells derived from canine fibroblasts permitted FCoV replication without apparent cytopathic effects. Thus, our reporter virus is useful for uncovering the infectivity of type I FCoV in different cell lines, including canine-derived cells. Surprisingly, we uncovered aberrant viral RNA transcription of rC3663 in A72 cells. Overall, we succeeded in obtaining infectious cDNA clones derived from type I FCoV that retained its virulence. Our recombinant FCoVs are powerful tools for increasing our understanding of the viral life cycle and pathogenesis of FIP-inducing type I FCoV.IMPORTANCE Feline coronavirus (FCoV) is one of the most significant coronaviruses, because this virus induces feline infectious peritonitis (FIP), which is a lethal disease in cats. Tissue culture-adapted type I FCoV often loses pathogenicity, which complicates research on type I FCoV-induced feline infectious peritonitis (FIP). Since we previously found that type I FCoV strain C3663 efficiently induces FIP in specific-pathogen-free cats, we established a reverse genetics system for the C3663 strain to obtain recombinant viruses in the present study. By using a reporter C3663 virus, we were able to examine the inhibitory effect of 68 compounds on C3663 replication in Fcwf-4 cells and infectivity in a canine-derived cell line. Interestingly, one canine cell line, A72, permitted FCoV replication but with low efficiency and aberrant viral gene expression.

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.

Identification and Characterization of a Human Coronavirus 229E Nonstructural Protein 8-Associated RNA 3'-Terminal Adenylyltransferase Activity

Previously, coronavirus nsp8 proteins were suggested to have template-dependent RNA polymerase activities resembling those of RNA primases or even canonical RNA-dependent RNA polymerases, while more recent studies have suggested an essential cofactor function of nsp8 (plus nsp7) for nsp12-mediated RNA-dependent RNA polymerase activity. In an effort to reconcile conflicting data from earlier studies, the study revisits coronavirus nsp8-associated activities using additional controls and proteins. The data obtained for three coronavirus nsp8 proteins provide evidence that the proteins share metal ion-dependent RNA 3′ polyadenylation activities that are greatly stimulated by a short oligo(U) stretch in the template strand. In contrast, nsp8 was found to be unable to select and incorporate appropriate (matching) nucleotides to produce cRNA products from heteropolymeric and other homooligomeric templates. While confirming the critical role of nsp8 in coronavirus replication, the study amends the list of activities mediated by coronavirus nsp8 proteins in the absence of other proteins.

Coronavirus nonstructural protein 8 (nsp8) has been suggested to have diverse activities, including noncanonical template-dependent polymerase activities. Here, we characterized a recombinant form of the human coronavirus 229E (HCoV-229E) nsp8 and found that the protein has metal ion-dependent RNA 3′-terminal adenylyltransferase (TATase) activity, while other nucleotides were not (or very inefficiently) transferred to the 3′ ends of single-stranded and (fully) double-stranded acceptor RNAs. Using partially double-stranded RNAs, very efficient TATase activity was observed if the opposite (template) strand contained a short 5′ oligo(U) sequence, while very little (if any) activity was detected for substrates with other homopolymeric or heteropolymeric sequences in the 5′ overhang. The oligo(U)-assisted/templated TATase activity on partial-duplex RNAs was confirmed for two other coronavirus nsp8 proteins, suggesting that the activity is conserved among coronaviruses. Replacement of a conserved Lys residue with Ala abolished the in vitro RNA-binding and TATase activities of nsp8 and caused a nonviable phenotype when the corresponding mutation was introduced into the HCoV-229E genome, confirming that these activities are mediated by nsp8 and critical for viral replication. In additional experiments, we obtained evidence that nsp8 has a pronounced specificity for adenylate and is unable to incorporate guanylate into RNA products, which strongly argues against the previously proposed template-dependent RNA polymerase activity of this protein. Given the presence of an oligo(U) stretch at the 5′ end of coronavirus minus-strand RNAs, it is tempting to speculate (but remains to be confirmed) that the nsp8-mediated TATase activity is involved in the 3′ polyadenylation of viral plus-strand RNAs.

IMPORTANCE Previously, coronavirus nsp8 proteins were suggested to have template-dependent RNA polymerase activities resembling those of RNA primases or even canonical RNA-dependent RNA polymerases, while more recent studies have suggested an essential cofactor function of nsp8 (plus nsp7) for nsp12-mediated RNA-dependent RNA polymerase activity. In an effort to reconcile conflicting data from earlier studies, the study revisits coronavirus nsp8-associated activities using additional controls and proteins. The data obtained for three coronavirus nsp8 proteins provide evidence that the proteins share metal ion-dependent RNA 3′ polyadenylation activities that are greatly stimulated by a short oligo(U) stretch in the template strand. In contrast, nsp8 was found to be unable to select and incorporate appropriate (matching) nucleotides to produce cRNA products from heteropolymeric and other homooligomeric templates. While confirming the critical role of nsp8 in coronavirus replication, the study amends the list of activities mediated by coronavirus nsp8 proteins in the absence of other proteins.

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.

Prediction of RNA-protein interactions using conjoint triad feature and chaos game representation

RNA-protein interactions (RPIs) play a very important role in a wide range of post-transcriptional regulations, and identifying whether a given RNA-protein pair can form interactions or not is a vital prerequisite for dissecting the regulatory mechanisms of functional RNAs. Currently, expensive and time-consuming biological assays can only determine a very small portion of all RPIs, which calls for computational approaches to help biologists efficiently and correctly find candidate RPIs. Here, we integrated a successful computing algorithm, conjoint triad feature (CTF), and another method, chaos game representation (CGR), for representing RNA-protein pairs and by doing so developed a prediction model based on these representations and random forest (RF) classifiers. When testing two benchmark datasets, RPI369 and RPI2241, the combined method (CTF+CGR) showed some superiority compared with four existing tools. Especially on RPI2241, the CTF+CGR method improved prediction accuracy (ACC) from 0.91 (the best record of all published works) to 0.95. When independently testing a newly constructed dataset, RPI1449, which only contained experimentally validated RPIs released between 2014 and 2016, our method still showed some generalization capability with an ACC of 0.75. Accordingly, we believe that our hybrid CTF+CGR method will be an important tool for predicting RPIs in the future.

Identification and molecular analysis of a lncRNA-HOTAIR transcript from secondary hair follicle of cashmere goat reveal integrated regulatory network with the expression regulated potentially by its promoter methylation

The HOTAIR transcript is transcribed from the antisense strand within the HOXC gene cluster, and it is thought to play a role in regulating the inductive capacity of dermal papilla cells during the reconstruction of hair-follicle. In the current investigation, we firstly isolated and characterized a lncRNA-HOTAIR transcript from the secondary hair follicle of cashmere goat. Also, we analyzed its transcriptional pattern and methylation level of HOTAIR gene promoter in secondary hair follicle of cashmere goat during anagen and telogen stages. Nucleotide composition analysis indicated that the contents of Adenine (A) and Thymine (T) are higher than that of Guanine (G) and Cytosine (C) in lncRNA-HOTAIR transcript of cashmere goat with the highest frequency distribution of AG nucleotide pair (8.06%). The regulatory network analysis showed a directly or indirectly complex regulatory relationships between lncRNA-HOTAIR of cashmere goat and its potential target molecules: miRNAs, mRNAs and proteins. Also, we showed that lncRNA-HOTAIR was properly transcribed at both anagen and telogen stages of secondary hair follicle of cashmere goat with the anagen being significantly higher than telogen in its expression, which suggest that lncRNA-HOTAIR transcript might be involved in the reconstruction of secondary hair follicle with the formation and growth of cashmere fiber. Taken together with methylation analysis of HOTAIR gene promoter, our data suggest that the promoter methylation of HOTAIR gene most likely is involved in its transcriptional suppression in secondary hair follicle of cashmere goat.

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.

Cyclophilins and cyclophilin inhibitors in nidovirus replication

Cyclophilins (Cyps) belong to the family of peptidyl-prolyl isomerases (PPIases). The PPIase activity of most Cyps is inhibited by the immunosuppressive drug cyclosporin A and several of its non-immunosuppressive analogs, which can also block the replication of nidoviruses (arteriviruses and coronaviruses). Cyclophilins have been reported to play an essential role in the replication of several other RNA viruses, including human immunodeficiency virus-1, hepatitis C virus, and influenza A virus. Likewise, the replication of various nidoviruses was reported to depend on Cyps or other PPIases. This review summarizes our current understanding of this class of nidovirus-host interactions, including the potential function of in particular CypA and the inhibitory effect of Cyp inhibitors. Also the involvement of the FK-506-binding proteins and parvulins is discussed. The nidovirus data are placed in a broader perspective by summarizing the most relevant data on Cyp interactions and Cyp inhibitors for other RNA viruses.

•

Nidovirus replication is inhibited by cyclophilin inhibitors.

•

Arterivirus replication depends on cyclophilin A.

•

Cyclosporin A blocks arterivirus RNA synthesis.

•

Using cyclophilin inhibitors against nidoviruses in vivo needs more investigation.

RPiRLS: Quantitative Predictions of RNA Interacting with Any Protein of Known Sequence

RNA-protein interactions (RPIs) have critical roles in numerous fundamental biological processes, such as post-transcriptional gene regulation, viral assembly, cellular defence and protein synthesis. As the number of available RNA-protein binding experimental data has increased rapidly due to high-throughput sequencing methods, it is now possible to measure and understand RNA-protein interactions by computational methods. In this study, we integrate a sequence-based derived kernel with regularized least squares to perform prediction. The derived kernel exploits the contextual information around an amino acid or a nucleic acid as well as the repetitive conserved motif information. We propose a novel machine learning method, called RPiRLS to predict the interaction between any RNA and protein of known sequences. For the RPiRLS classifier, each protein sequence comprises up to 20 diverse amino acids but for the RPiRLS-7G classifier, each protein sequence is represented by using 7-letter reduced alphabets based on their physiochemical properties. We evaluated both methods on a number of benchmark data sets and compared their performances with two newly developed and state-of-the-art methods, RPI-Pred and IPMiner. On the non-redundant benchmark test sets extracted from the PRIDB, the RPiRLS method outperformed RPI-Pred and IPMiner in terms of accuracy, specificity and sensitivity. Further, RPiRLS achieved an accuracy of 92% on the prediction of lncRNA-protein interactions. The proposed method can also be extended to construct RNA-protein interaction networks. The RPiRLS web server is freely available at http://bmc.med.stu.edu.cn/RPiRLS.

Prediction of lncRNA-protein interactions using HeteSim scores based on heterogeneous networks

Massive studies have indicated that long non-coding RNAs (lncRNAs) are critical for the regulation of cellular biological processes by binding with RNA-related proteins. However, only a few experimentally supported lncRNA-protein associations have been reported. Existing network-based methods are typically focused on intrinsic features of lncRNA and protein but ignore the information implicit in the topologies of biological networks associated with lncRNAs. Considering the limitations in previous methods, we propose PLPIHS, an effective computational method for Predicting lncRNA-Protein Interactions using HeteSim Scores. PLPIHS uses the HeteSim measure to calculate the relatedness score for each lncRNA-protein pair in the heterogeneous network, which consists of lncRNA-lncRNA similarity network, lncRNA-protein association network and protein-protein interaction network. An SVM classifier to predict lncRNA-protein interactions is built with the HeteSim scores. The results show that PLPIHS performs significantly better than the existing state-of-the-art approaches and achieves an AUC score of 0.97 in the leave-one-out validation test. We also compare the performances of networks with different connectivity density and find that PLPIHS performs well across all the networks. Furthermore, we use the proposed method to identify the related proteins for lncRNA MALAT1. Highly-ranked proteins are verified by the biological studies and demonstrate the effectiveness of our method.

Nidovirus RNA polymerases: Complex enzymes handling exceptional RNA genomes

Coronaviruses and arteriviruses are distantly related human and animal pathogens that belong to the order Nidovirales. Nidoviruses are characterized by their polycistronic plus-stranded RNA genome, the production of subgenomic mRNAs and the conservation of a specific array of replicase domains, including key RNA-synthesizing enzymes. Coronaviruses (26-34 kilobases) have the largest known RNA genomes and their replication presumably requires a processive RNA-dependent RNA polymerase (RdRp) and enzymatic functions that suppress the consequences of the typically high error rate of viral RdRps. The arteriviruses have significantly smaller genomes and form an intriguing package with the coronaviruses to analyse viral RdRp evolution and function. The RdRp domain of nidoviruses resides in a cleavage product of the replicase polyprotein named non-structural protein (nsp) 12 in coronaviruses and nsp9 in arteriviruses. In all nidoviruses, the C-terminal RdRp domain is linked to a conserved N-terminal domain, which has been coined NiRAN (nidovirus RdRp-associated nucleotidyl transferase). Although no structural information is available, the functional characterization of the nidovirus RdRp and the larger enzyme complex of which it is part, has progressed significantly over the past decade. In coronaviruses several smaller, non-enzymatic nsps were characterized that direct RdRp function, while a 3'-to-5' exoribonuclease activity in nsp14 was implicated in fidelity. In arteriviruses, the nsp1 subunit was found to maintain the balance between genome replication and subgenomic mRNA production. Understanding RdRp behaviour and interactions during RNA synthesis and subsequent processing will be key to rationalising the evolutionary success of nidoviruses and the development of antiviral strategies.

Long-Range Communication between Different Functional Sites in the Picornaviral 3C Protein

The 3C protein is a master regulator of the picornaviral infection cycle, responsible for both cleaving viral and host proteins, and interacting with genomic RNA replication elements. Here we use nuclear magnetic resonance spectroscopy and molecular dynamics simulations to show that 3C is conformationally dynamic across multiple timescales. Binding of peptide and RNA lead to structural dynamics changes at both the protease active site and the RNA-binding site, consistent with these sites being dynamically coupled. Indeed, binding of RNA influences protease activity, and likewise, interactions at the active site affect RNA binding. We propose that RNA and peptide binding re-shapes the conformational energy landscape of 3C to regulate subsequent functions, including formation of complexes with other viral proteins. The observed channeling of the 3C energy landscape may be important for regulation of the viral infection cycle.

The Nonstructural Proteins Directing Coronavirus RNA Synthesis and Processing

Coronaviruses are animal and human pathogens that can cause lethal zoonotic infections like SARS and MERS. They have polycistronic plus-stranded RNA genomes and belong to the order Nidovirales, a diverse group of viruses for which common ancestry was inferred from the common principles underlying their genome organization and expression, and from the conservation of an array of core replicase domains, including key RNA-synthesizing enzymes. Coronavirus genomes (~ 26–32 kilobases) are the largest RNA genomes known to date and their expansion was likely enabled by acquiring enzyme functions that counter the commonly high error frequency of viral RNA polymerases. The primary functions that direct coronavirus RNA synthesis and processing reside in nonstructural protein (nsp) 7 to nsp16, which are cleavage products of two large replicase polyproteins translated from the coronavirus genome. Significant progress has now been made regarding their structural and functional characterization, stimulated by technical advances like improved methods for bioinformatics and structural biology, in vitro enzyme characterization, and site-directed mutagenesis of coronavirus genomes. Coronavirus replicase functions include more or less universal activities of plus-stranded RNA viruses, like an RNA polymerase (nsp12) and helicase (nsp13), but also a number of rare or even unique domains involved in mRNA capping (nsp14, nsp16) and fidelity control (nsp14). Several smaller subunits (nsp7–nsp10) act as crucial cofactors of these enzymes and contribute to the emerging “nsp interactome.” Understanding the structure, function, and interactions of the RNA-synthesizing machinery of coronaviruses will be key to rationalizing their evolutionary success and the development of improved control strategies.

Coronavirus cis-Acting RNA Elements

Coronaviruses have exceptionally large RNA genomes of approximately 30 kilobases. Genome replication and transcription is mediated by a multisubunit protein complex comprised of more than a dozen virus-encoded proteins. The protein complex is thought to bind specific cis-acting RNA elements primarily located in the 5'- and 3'-terminal genome regions and upstream of the open reading frames located in the 3'-proximal one-third of the genome. Here, we review our current understanding of coronavirus cis-acting RNA elements, focusing on elements required for genome replication and packaging. Recent bioinformatic, biochemical, and genetic studies suggest a previously unknown level of conservation of cis-acting RNA structures among different coronavirus genera and, in some cases, even beyond genus boundaries. Also, there is increasing evidence to suggest that individual cis-acting elements may be part of higher-order RNA structures involving long-range and dynamic RNA-RNA interactions between RNA structural elements separated by thousands of nucleotides in the viral genome. We discuss the structural and functional features of these cis-acting RNA elements and their specific functions in coronavirus RNA synthesis.

Identification of RNA Binding Proteins Associated with Dengue Virus RNA in Infected Cells Reveals Temporally Distinct Host Factor Requirements

BACKGROUND

There are currently no vaccines or antivirals available for dengue virus infection, which can cause dengue hemorrhagic fever and death. A better understanding of the host pathogen interaction is required to develop effective therapies to treat DENV. In particular, very little is known about how cellular RNA binding proteins interact with viral RNAs. RNAs within cells are not naked; rather they are coated with proteins that affect localization, stability, translation and (for viruses) replication.

METHODOLOGY/PRINCIPAL FINDINGS

Seventy-nine novel RNA binding proteins for dengue virus (DENV) were identified by cross-linking proteins to dengue viral RNA during a live infection in human cells. These cellular proteins were specific and distinct from those previously identified for poliovirus, suggesting a specialized role for these factors in DENV amplification. Knockdown of these proteins demonstrated their function as viral host factors, with evidence for some factors acting early, while others late in infection. Their requirement by DENV for efficient amplification is likely specific, since protein knockdown did not impair the cell fitness for viral amplification of an unrelated virus. The protein abundances of these host factors were not significantly altered during DENV infection, suggesting their interaction with DENV RNA was due to specific recruitment mechanisms. However, at the global proteome level, DENV altered the abundances of proteins in particular classes, including transporter proteins, which were down regulated, and proteins in the ubiquitin proteasome pathway, which were up regulated.

CONCLUSIONS/SIGNIFICANCE

The method for identification of host factors described here is robust and broadly applicable to all RNA viruses, providing an avenue to determine the conserved or distinct mechanisms through which diverse viruses manage the viral RNA within cells. This study significantly increases the number of cellular factors known to interact with DENV and reveals how DENV modulates and usurps cellular proteins for efficient amplification.

Evidence supporting a premature termination mechanism for subgenomic RNA transcription in Pelargonium line pattern virus: identification of a critical long-range RNA-RNA interaction and functional variants through mutagenesis

Pelargonium line pattern virus (PLPV) is a plus-strand RNA virus that has been proposed as type species of a tentative new genus, Pelarspovirus, in the family Tombusviridae. One of the singular traits of members of this prospective genus is the production of a unique subgenomic (sg) mRNA that is structurally and functionally tricistronic. Here, we have aimed to get insights into the mechanism that governs PLPV sg mRNA transcription. A long-range RNA-RNA interaction that is critical for the process has been identified through RNA folding predictions and mutational analysis of the viral genome. Such interaction seems to occur in the plus-strand, likely acts in cis, and specifically mediates the synthesis of sg RNA-sized minus-strand. The accumulation of this RNA species is easily detectable in plants and its generation can be uncoupled from that of the plus-strand sg mRNA. All these data together with the observation that 5' ends of PLPV genomic and sg mRNAs have sequence resemblances (as expected if both act as promoters in the corresponding minus-strand), support that premature termination is the mechanism underlying PLPV sg mRNA formation.

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.

Sequestration of RNA by grass carp Ctenopharyngodon idella TIA1 is associated with its positive role in facilitating grass carp reovirus infection

Previous report demonstrated that grass carp reovirus (GCRV) infection resulted in unlinking cellular stress granule formation from aggregation of grass carp Ctenopharyngodon idella TIA1 (CiTIA1). Here, we provided evidence to show that CiTIA1 bound to synthesized ssRNA and dsRNA in vitro. Both GST-pull down assay and RNA immunoprecipitation analysis confirmed the association between GCRV-specific RNA and GST-tagged CiTIA1 in C. idella kidney (CIK) cells. Furthermore, CiTIA1 was shown to protect dsRNA of virus-origin from degradation in CIK cells through Northern blot analysis. Finally, transient overexpression of CiTIA1 enhanced the replication efficiency of GCRV in CIK cells. Taken together, our results suggested that cellular CiTIA1 might facilitate GCRV replication through sequestrating and protecting viral RNA from degradation.

Design of Potential RNAi (miRNA and siRNA) Molecules for Middle East Respiratory Syndrome Coronavirus (MERS-CoV) Gene Silencing by Computational Method

The Middle East respiratory syndrome coronavirus (MERS-CoV) is a virus that manifests itself in viral infection with fever, cough, shortness of breath, renal failure and severe acute pneumonia, which often result in a fatal outcome. MERS-CoV has been shown to spread between people who are in close contact. Transmission from infected patients to healthcare personnel has also been observed and is irredeemable with present technology. Genetic studies on MERS-CoV have shown that ORF1ab encodes replicase polyproteins and play a foremost role in viral infection. Therefore, ORF1ab replicase polyprotein may be used as a suitable target for disease control. Viral activity can be controlled by RNA interference (RNAi) technology, a leading method for post transcriptional gene silencing in a sequence-specific manner. However, there is a genetic inconsistency in different viral isolates; it is a great challenge to design potential RNAi (miRNA and siRNA) molecules which can silence the respective target genes rather than any other viral gene simultaneously. In the current study, four effective miRNA and five siRNA molecules for silencing of nine different strains of MERS-CoV were rationally designed and corroborated using computational methods, which might lead to knockdown the activity of virus. siRNA and miRNA molecules were predicted against ORF1ab gene of different strains of MERS-CoV as effective candidate using computational methods. Thus, this method may provide an insight for the chemical synthesis of antiviral RNA molecule for the treatment of MERS-CoV, at genomic level.

Nucleocapsid phosphorylation and RNA helicase DDX1 recruitment enables coronavirus transition from discontinuous to continuous transcription

Coronaviruses contain a positive-sense single-stranded genomic (g) RNA, which encodes nonstructural proteins. Several subgenomic mRNAs (sgmRNAs) encoding structural proteins are generated by template switching from the body transcription regulatory sequence (TRS) to the leader TRS. The process preferentially generates shorter sgmRNA. Appropriate readthrough of body TRSs is required to produce longer sgmRNAs and full-length gRNA. We find that phosphorylation of the viral nucleocapsid (N) by host glycogen synthase kinase-3 (GSK-3) is required for template switching. GSK-3 inhibition selectively reduces the generation of gRNA and longer sgmRNAs, but not shorter sgmRNAs. N phosphorylation allows recruitment of the RNA helicase DDX1 to the phosphorylated-N-containing complex, which facilitates template readthrough and enables longer sgmRNA synthesis. DDX1 knockdown or loss of helicase activity markedly reduces the levels of longer sgmRNAs. Thus, coronaviruses employ a unique strategy for the transition from discontinuous to continuous transcription to ensure balanced sgmRNAs and full-length gRNA synthesis.

Genomic and single nucleotide polymorphism analysis of infectious bronchitis coronavirus

Infectious bronchitis virus (IBV) is a Gammacoronavirus that causes a highly contagious respiratory disease in chickens. A QX-like strain was analysed by high-throughput Illumina sequencing and genetic variation across the entire viral genome was explored at the sub-consensus level by single nucleotide polymorphism (SNP) analysis. Thirteen open reading frames (ORFs) in the order 5'-UTR-1a-1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3'UTR were predicted. The relative frequencies of missense: silent SNPs were calculated to obtain a comparative measure of variability in specific genes. The most variable ORFs in descending order were E, 3b, 5'UTR, N, 1a, S, 1ab, M, 4c, 5a, 6b. The E and 3b protein products play key roles in coronavirus virulence, and RNA folding demonstrated that the mutations in the 5'UTR did not alter the predicted secondary structure. The frequency of SNPs in the Spike (S) protein ORF of 0.67% was below the genomic average of 0.76%. Only three SNPS were identified in the S1 subunit, none of which were located in hypervariable region (HVR) 1 or HVR2. The S2 subunit was considerably more variable containing 87% of the polymorphisms detected across the entire S protein. The S2 subunit also contained a previously unreported multi-A insertion site and a stretch of four consecutive mutated amino acids, which mapped to the stalk region of the spike protein. Template-based protein structure modelling produced the first theoretical model of the IBV spike monomer. Given the lack of diversity observed at the sub-consensus level, the tenet that the HVRs in the S1 subunit are very tolerant of amino acid changes produced by genetic drift is questioned.

Nuclear proteins hijacked by mammalian cytoplasmic plus strand RNA viruses

Plus strand RNA viruses that replicate in the cytoplasm face challenges in supporting the numerous biosynthetic functions required for replication and propagation. Most of these viruses are genetically simple and rely heavily on co-opting cellular proteins, particularly cellular RNA-binding proteins, into new roles for support of virus infection at the level of virus-specific translation, and building RNA replication complexes. In the course of infectious cycles many nuclear-cytoplasmic shuttling proteins of mostly nuclear distribution are detained in the cytoplasm by viruses and re-purposed for their own gain. Many mammalian viruses hijack a common group of the same factors. This review summarizes recent gains in our knowledge of how cytoplasmic RNA viruses use these co-opted host nuclear factors in new functional roles supporting virus translation and virus RNA replication and common themes employed between different virus groups.

•

Nuclear shuttling host proteins are commonly hijacked by RNA viruses to support replication.

•

A limited group of ubiquitous RNA binding proteins are commonly hijacked by a broad range of viruses.

•

Key virus proteins alter roles of RNA binding proteins in different stages of virus replication.

Identification of a gamma interferon-activated inhibitor of translation-like RNA motif at the 3' end of the transmissible gastroenteritis coronavirus genome modulating innate immune response

A 32-nucleotide (nt) RNA motif located at the 3′ end of the transmissible gastroenteritis coronavirus (TGEV) genome was found to specifically interact with the host proteins glutamyl-prolyl-tRNA synthetase (EPRS) and arginyl-tRNA synthetase (RRS). This RNA motif has high homology in sequence and secondary structure with the gamma interferon-activated inhibitor of translation (GAIT) element, which is located at the 3′ end of several mRNAs encoding proinflammatory proteins. The GAIT element is involved in the translation silencing of these mRNAs through its interaction with the GAIT complex (EPRS, heterogeneous nuclear ribonucleoprotein Q, ribosomal protein L13a, and glyceraldehyde 3-phosphate dehydrogenase) to favor the resolution of inflammation. Interestingly, we showed that the viral RNA motif bound the GAIT complex and inhibited the in vitro translation of a chimeric mRNA containing this RNA motif. To our knowledge, this is the first GAIT-like motif described in a positive RNA virus. To test the functional role of the GAIT-like RNA motif during TGEV infection, a recombinant coronavirus harboring mutations in this motif was engineered and characterized. Mutations of the GAIT-like RNA motif did not affect virus growth in cell cultures. However, an exacerbated innate immune response, mediated by the melanoma differentiation-associated gene 5 (MDA5) pathway, was observed in cells infected with the mutant virus compared with the response observed in cells infected with the parental virus. Furthermore, the mutant virus was more sensitive to beta interferon than the parental virus. All together, these data strongly suggested that the viral GAIT-like RNA motif modulates the host innate immune response.

The innate immune response is the first line of antiviral defense that culminates with the synthesis of interferon and proinflammatory cytokines to limit virus replication. Coronaviruses encode several proteins that interfere with the innate immune response at different levels, but to date, no viral RNA counteracting antiviral response has been described. In this work, we have characterized a 32-nt RNA motif located at the 3′ end of the TGEV genome that specifically interacted with EPRS and RRS. This RNA motif presented high homology with the GAIT element, involved in the modulation of the inflammatory response. Moreover, the disruption of the viral GAIT-like RNA motif led to an exacerbated innate immune response triggered by MDA5, indicating that the GAIT-like RNA motif counteracts the host innate immune response. These novel findings may be of relevance for other coronaviruses and could serve as the basis for the development of novel antiviral strategies.

RPI-Pred: predicting ncRNA-protein interaction using sequence and structural information

RNA-protein complexes are essential in mediating important fundamental cellular processes, such as transport and localization. In particular, ncRNA-protein interactions play an important role in post-transcriptional gene regulation like mRNA localization, mRNA stabilization, poly-adenylation, splicing and translation. The experimental methods to solve RNA-protein interaction prediction problem remain expensive and time-consuming. Here, we present the RPI-Pred (RNA-protein interaction predictor), a new support-vector machine-based method, to predict protein-RNA interaction pairs, based on both the sequences and structures. The results show that RPI-Pred can correctly predict RNA-protein interaction pairs with ∼94% prediction accuracy when using sequence and experimentally determined protein and RNA structures, and with ∼83% when using sequences and predicted protein and RNA structures. Further, our proposed method RPI-Pred was superior to other existing ones by predicting more experimentally validated ncRNA-protein interaction pairs from different organisms. Motivated by the improved performance of RPI-Pred, we further applied our method for reliable construction of ncRNA-protein interaction networks. The RPI-Pred is publicly available at: http://ctsb.is.wfubmc.edu/projects/rpi-pred.

What we know but do not understand about nidovirus helicases

•

The ubiquitous nidovirus helicase is a multi-functional enzyme of superfamily 1.

•

Its unique N-terminal domain is most similar to the Upf1 multinuclear zinc-binding domain.

•

It has been implicated in replication, transcription, virion biogenesis, translation and post-transcriptional viral RNA processing.

•

Four different classes of antiviral compounds targeting the helicase have been identified.

Helicases are versatile NTP-dependent motor proteins of monophyletic origin that are found in all kingdoms of life. Their functions range from nucleic acid duplex unwinding to protein displacement and double-strand translocation. This explains their participation in virtually every metabolic process that involves nucleic acids, including DNA replication, recombination and repair, transcription, translation, as well as RNA processing. Helicases are encoded by all plant and animal viruses with a positive-sense RNA genome that is larger than 7 kb, indicating a link to genome size evolution in this virus class. Viral helicases belong to three out of the six currently recognized superfamilies, SF1, SF2, and SF3. Despite being omnipresent, highly conserved and essential, only a few viral helicases, mostly from SF2, have been studied extensively. In general, their specific roles in the viral replication cycle remain poorly understood at present. The SF1 helicase protein of viruses classified in the order Nidovirales is encoded in replicase open reading frame 1b (ORF1b), which is translated to give rise to a large polyprotein following a ribosomal frameshift from the upstream ORF1a. Proteolytic processing of the replicase polyprotein yields a dozen or so mature proteins, one of which includes a helicase. Its hallmark is the presence of an N-terminal multi-nuclear zinc-binding domain, the nidoviral genetic marker and one of the most conserved domains across members of the order. This review summarizes biochemical, structural, and genetic data, including drug development studies, obtained using helicases originating from several mammalian nidoviruses, along with the results of the genomics characterization of a much larger number of (putative) helicases of vertebrate and invertebrate nidoviruses. In the context of our knowledge of related helicases of cellular and viral origin, it discusses the implications of these results for the protein's emerging critical function(s) in nidovirus evolution, genome replication and expression, virion biogenesis, and possibly also post-transcriptional processing of viral RNAs. Using our accumulated knowledge and highlighting gaps in our data, concepts and approaches, it concludes with a perspective on future research aimed at elucidating the role of helicases in the nidovirus replication cycle.

Design of potential RNAi (miRNA and siRNA) molecules for Middle East respiratory syndrome coronavirus (MERS-CoV) gene silencing by computational method

The Middle East respiratory syndrome coronavirus (MERS-CoV) is a virus that manifests itself in viral infection with fever, cough, shortness of breath, renal failure and severe acute pneumonia, which often result in a fatal outcome. MERS-CoV has been shown to spread between people who are in close contact. Transmission from infected patients to healthcare personnel has also been observed and is irredeemable with present technology. Genetic studies on MERS-CoV have shown that ORF 1ab encodes replicase polyproteins and play a foremost role in viral infection. Therefore, ORF 1ab replicase polyprotein may be used as suitable target for disease control. Viral activity can be controlled by RNA interference (RNAi) technology, a leading method for post transcriptional gene silencing in a sequence specific manner. However, there is a genetic inconsistency in different viral isolates; it is a great challenge to design potential RNAi (miRNA and siRNA) molecules which can silence the respective target genes rather than any other viral gene simultaneously. In current study four effective miRNA and five siRNA molecules for silencing of nine different strains of MERS-CoV were rationally designed and corroborated using computational methods, which might lead to knockdown the activity of virus. siRNA and miRNA molecules were predicted against ORF1ab gene of different strains of MERS-CoV as effective candidate using computational methods. Thus, this method may provide an insight for the chemical synthesis of antiviral RNA molecule for the treatment of MERS-CoV, at genomic level.

Design of potential RNAi (miRNA and siRNA) molecules for Middle East respiratory syndrome coronavirus (MERS-CoV) gene silencing by computational method

The Middle East respiratory syndrome coronavirus (MERS-CoV) is a virus that manifests itself in viral infection with fever, cough, shortness of breath, renal failure and severe acute pneumonia, which often result in a fatal outcome. MERS-CoV has been shown to spread between people who are in close contact. Transmission from infected patients to healthcare personnel has also been observed and is irredeemable with present technology. Genetic studies on MERS-CoV have shown that ORF 1ab encodes replicase polyproteins and play a foremost role in viral infection. Therefore, ORF 1ab replicase polyprotein may be used as suitable target for disease control. Viral activity can be controlled by RNA interference (RNAi) technology, a leading method for post transcriptional gene silencing in a sequence specific manner. However, there is a genetic inconsistency in different viral isolates; it is a great challenge to design potential RNAi (miRNA and siRNA) molecules which can silence the respective target genes rather than any other viral gene simultaneously. In current study four effective miRNA and five siRNA molecules for silencing of nine different strains of MERS-CoV were rationally designed and corroborated using computational methods, which might lead to knockdown the activity of virus. siRNA and miRNA molecules were predicted against ORF1ab gene of different strains of MERS-CoV as effective candidate using computational methods. Thus, this method may provide an insight for the chemical synthesis of antiviral RNA molecule for the treatment of MERS-CoV, at genomic level.

Engineering the PRRS virus genome: updates and perspectives

•

We review PRRSV infectious clones and their applications.

•

14 infectious clones are available so far for genotypes I and II.

•

Genomic mutations, insertions, deletions, and replacements are successful.

•

We discuss advances and utilization of PRRSV reverse genetics and future potential.

Porcine reproductive and respiratory syndrome virus (PRRSV) is endemic in most pig producing countries worldwide and causes enormous economic losses to the pork industry. Infectious clones for PRRSV have been constructed, and so far at least 14 different infectious clones are available representing both genotypes I and II. Two strategies have been taken for progeny reconstitution: RNA transfection and DNA transfection. Mutations, insertions, deletions, and replacements of the viral genome have been employed to study the structure function relationship, foreign gene expression, functional complementation, and virulence determinants. Essential regions and non-essential regions for viral replication have been identified in both the coding regions and non-encoding regions. Foreign sequences have successfully been inserted into the nsp2 and N regions and in the space between ORF1b and ORF2a. Chimeras between member viruses in the family Arteriviridae have also been constructed and utilized to study cell tropism and functional complementation. This review discusses the advances and utilization of PRRSV reverse genetics and its potential for future research.

RNA structure analysis of alphacoronavirus terminal genome regions

•

Review of current knowledge of cis-acting RNA elements essential to coronavirus replication.

•

Identification of RNA structural elements in alphacoronavirus terminal genome regions.

•

Discussion of intra- and intergeneric conservation of genomic cis-acting RNA elements in alpha- and betacoronaviruses.

Coronavirus genome replication is mediated by a multi-subunit protein complex that is comprised of more than a dozen virally encoded and several cellular proteins. Interactions of the viral replicase complex with cis-acting RNA elements located in the 5′ and 3′-terminal genome regions ensure the specific replication of viral RNA. Over the past years, boundaries and structures of cis-acting RNA elements required for coronavirus genome replication have been extensively characterized in betacoronaviruses and, to a lesser extent, other coronavirus genera. Here, we review our current understanding of coronavirus cis-acting elements located in the terminal genome regions and use a combination of bioinformatic and RNA structure probing studies to identify and characterize putative cis-acting RNA elements in alphacoronaviruses. The study suggests significant RNA structure conservation among members of the genus Alphacoronavirus but also across genus boundaries. Overall, the conservation pattern identified for 5′ and 3′-terminal RNA structural elements in the genomes of alpha- and betacoronaviruses is in agreement with the widely used replicase polyprotein-based classification of the Coronavirinae, suggesting co-evolution of the coronavirus replication machinery with cognate cis-acting RNA elements.

Structural basis for the regulatory function of a complex zinc-binding domain in a replicative arterivirus helicase resembling a nonsense-mediated mRNA decay helicase

All positive-stranded RNA viruses with genomes >∼7 kb encode helicases, which generally are poorly characterized. The core of the nidovirus superfamily 1 helicase (HEL1) is associated with a unique N-terminal zinc-binding domain (ZBD) that was previously implicated in helicase regulation, genome replication and subgenomic mRNA synthesis. The high-resolution structure of the arterivirus helicase (nsp10), alone and in complex with a polynucleotide substrate, now provides first insights into the structural basis for nidovirus helicase function. A previously uncharacterized domain 1B connects HEL1 domains 1A and 2A to a long linker of ZBD, which further consists of a novel RING-like module and treble-clef zinc finger, together coordinating three Zn atoms. On substrate binding, major conformational changes were evident outside the HEL1 domains, notably in domain 1B. Structural characterization, mutagenesis and biochemistry revealed that helicase activity depends on the extensive relay of interactions between the ZBD and HEL1 domains. The arterivirus helicase structurally resembles the cellular Upf1 helicase, suggesting that nidoviruses may also use their helicases for post-transcriptional quality control of their large RNA genomes.

Unmasking the information encoded as structural motifs of viral RNA genomes: a potential antiviral target

RNA viruses show enormous capacity to evolve and adapt to new cellular and molecular contexts, a consequence of mutations arising from errors made by viral RNA-dependent RNA polymerase during replication. Sequence variation must occur, however, without compromising functions essential for the completion of the viral cycle. RNA viruses are safeguarded in this respect by their genome carrying conserved information that does not code only for proteins but also for the formation of structurally conserved RNA domains that directly perform these critical functions. Functional RNA domains can interact with other regions of the viral genome and/or proteins to direct viral translation, replication and encapsidation. They are therefore potential targets for novel therapeutic strategies. This review summarises our knowledge of the functional RNA domains of human RNA viruses and examines the achievements made in the design of antiviral compounds that interfere with their folding and therefore their function.

Arterivirus molecular biology and pathogenesis

Arteriviruses are positive-stranded RNA viruses that infect mammals. They can cause persistent or asymptomatic infections, but also acute disease associated with a respiratory syndrome, abortion or lethal haemorrhagic fever. During the past two decades, porcine reproductive and respiratory syndrome virus (PRRSV) and, to a lesser extent, equine arteritis virus (EAV) have attracted attention as veterinary pathogens with significant economic impact. Particularly noteworthy were the 'porcine high fever disease' outbreaks in South-East Asia and the emergence of new virulent PRRSV strains in the USA. Recently, the family was expanded with several previously unknown arteriviruses isolated from different African monkey species. At the molecular level, arteriviruses share an intriguing but distant evolutionary relationship with coronaviruses and other members of the order Nidovirales. Nevertheless, several of their characteristics are unique, including virion composition and structure, and the conservation of only a subset of the replicase domains encountered in nidoviruses with larger genomes. During the past 15 years, the advent of reverse genetics systems for EAV and PRRSV has changed and accelerated the structure-function analysis of arterivirus RNA and protein sequences. These systems now also facilitate studies into host immune responses and arterivirus immune evasion and pathogenesis. In this review, we have summarized recent advances in the areas of arterivirus genome expression, RNA and protein functions, virion architecture, virus-host interactions, immunity, and pathogenesis. We have also briefly reviewed the impact of these advances on disease management, the engineering of novel candidate live vaccines and the diagnosis of arterivirus infection.

The cellular interactome of the coronavirus infectious bronchitis virus nucleocapsid protein and functional implications for virus biology

The coronavirus nucleocapsid (N) protein plays a multifunctional role in the virus life cycle, from regulation of replication and transcription and genome packaging to modulation of host cell processes. These functions are likely to be facilitated by interactions with host cell proteins. The potential interactome of the infectious bronchitis virus (IBV) N protein was mapped using stable isotope labeling with amino acids in cell culture (SILAC) coupled to a green fluorescent protein-nanotrap pulldown methodology and liquid chromatography-tandem mass spectrometry. The addition of the SILAC label allowed discrimination of proteins that were likely to specifically bind to the N protein over background binding. Overall, 142 cellular proteins were selected as potentially binding to the N protein, many as part of larger possible complexes. These included ribosomal proteins, nucleolar proteins, translation initiation factors, helicases, and hnRNPs. The association of selected cellular proteins with IBV N protein was confirmed by immunoblotting, cosedimentation, and confocal microscopy. Further, the localization of selected proteins in IBV-infected cells as well as their activity during virus infection was assessed by small interfering RNA-mediated depletion, demonstrating the functional importance of cellular proteins in the biology of IBV. This interactome not only confirms previous observations made with other coronavirus and IBV N proteins with both overexpressed proteins and infectious virus but also provides novel data that can be exploited to understand the interaction between the virus and the host cell.

Altering SARS coronavirus frameshift efficiency affects genomic and subgenomic RNA production

In previous studies, differences in the amount of genomic and subgenomic RNA produced by coronaviruses with mutations in the programmed ribosomal frameshift signal of ORF1a/b were observed. It was not clear if these differences were due to changes in genomic sequence, the protein sequence or the frequency of frameshifting. Here, viruses with synonymous codon changes are shown to produce different ratios of genomic and subgenomic RNA. These findings demonstrate that the protein sequence is not the primary cause of altered genomic and subgenomic RNA production. The synonymous codon changes affect both the structure of the frameshift signal and frameshifting efficiency. Small differences in frameshifting efficiency result in dramatic differences in genomic RNA production and TCID₅₀ suggesting that the frameshifting frequency must stay above a certain threshold for optimal virus production. The data suggest that either the RNA sequence or the ratio of viral proteins resulting from different levels of frameshifting affects viral replication.

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.