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

Assembly of SARS-CoV-2 nucleocapsid protein with nucleic acid

The viral genome of SARS-CoV-2 is packaged by the nucleocapsid (N-)protein into ribonucleoprotein particles (RNPs), 38 ± 10 of which are contained in each virion. Their architecture has remained unclear due to the pleomorphism of RNPs, the high flexibility of N-protein intrinsically disordered regions, and highly multivalent interactions between viral RNA and N-protein binding sites in both N-terminal (NTD) and C-terminal domain (CTD). Here we explore critical interaction motifs of RNPs by applying a combination of biophysical techniques to ancestral and mutant proteins binding different nucleic acids in an in vitro assay for RNP formation, and by examining nucleocapsid protein variants in a viral assembly assay. We find that nucleic acid-bound N-protein dimers oligomerize via a recently described protein-protein interface presented by a transient helix in its long disordered linker region between NTD and CTD. The resulting hexameric complexes are stabilized by multivalent protein-nucleic acid interactions that establish crosslinks between dimeric subunits. Assemblies are stabilized by the dimeric CTD of N-protein offering more than one binding site for stem-loop RNA. Our study suggests a model for RNP assembly where N-protein scaffolding at high density on viral RNA is followed by cooperative multimerization through protein-protein interactions in the disordered linker.

Classification, replication, and transcription of Nidovirales

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

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

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

3'UTR of ALV-J can affect viral replication through promoting transcription and mRNA nuclear export

IMPORTANCE

3'UTRs can affect gene transcription and post-transcriptional regulation in multiple ways, further influencing the function of proteins in a unique manner. Recently, ALV-J has been mutating and evolving rapidly, especially the 3'UTR of viral genome. Meanwhile, clinical symptoms caused by ALV-J have changed significantly. In this study, we found that the ALV-J strains containing △-r-TM-type 3'UTR are the most abundant. By constructing ALV-J infectious clones and subgenomic vectors containing different 3'UTRs, we prove that 3'UTRs directly affect viral tissue preference and can promote virus replication as an enhancer. ALV-J strain containing 3'UTR of △-r-TM proliferated fastest in primary cells. All five forms of 3'UTRs can assist intron-containing viral mRNA nuclear export, with similar efficiency. ALV-J mRNA half-life is not influenced by different 3'UTRs. Our results dissect the roles of 3'UTR on regulating viral replication and pathogenicity, providing novel insights into potential anti-viral strategies.

RNA structure and multiple weak interactions balance the interplay between RNA binding and phase separation of SARS-CoV-2 nucleocapsid

The nucleocapsid (N) protein of SARS-CoV-2 binds viral RNA, condensing it inside the virion, and phase separating with RNA to form liquid-liquid condensates. There is little consensus on what differentiates sequence-independent N-RNA interactions in the virion or in liquid droplets from those with specific genomic RNA (gRNA) motifs necessary for viral function inside infected cells. To identify the RNA structures and the N domains responsible for specific interactions and phase separation, we use the first 1,000 nt of viral RNA and short RNA segments designed as models for single-stranded and paired RNA. Binding affinities estimated from fluorescence anisotropy of these RNAs to the two-folded domains of N (the NTD and CTD) and comparison to full-length N demonstrate that the NTD binds preferentially to single-stranded RNA, and while it is the primary RNA-binding site, it is not essential to phase separation. Nuclear magnetic resonance spectroscopy identifies two RNA-binding sites on the NTD: a previously characterized site and an additional although weaker RNA-binding face that becomes prominent when binding to the primary site is weak, such as with dsRNA or a binding-impaired mutant. Phase separation assays of nucleocapsid domains with double-stranded and single-stranded RNA structures support a model where multiple weak interactions, such as with the CTD or the NTD's secondary face promote phase separation, while strong, specific interactions do not. These studies indicate that both strong and multivalent weak N-RNA interactions underlie the multifunctional abilities of N.

Inhibiting mechanism of Alpiniae oxyphyllae fructus polysaccharide 3 against the replication of porcine epidemic diarrhea virus

Porcine epidemic diarrhea virus (PEDV) causes diarrhea, vomiting, and death in piglets. Our previous study has revealed the anti-PEDV activity of Alpiniae oxyphyllae fructus polysaccharide 3 (AOFP3). However, it is still unknown whether AOFP3 can inhibit the replication of PEDV. Therefore, the effect of AOFP3 on PEDV replication was investigated in the present study, along with analysis of viral RdRp activity and expression of hnRNP A1 by RNA polymerase activity assay in vitro, RIP assay, and Western blotting. The results showed that both the PEDV gene and protein levels in IPEC-J2 cells decreased with AOFP3 treatment. In addition, AOFP3 significantly reduced PEDV's replication by down-regulating the activity of PEDV RdRp and reducing the expression of hnRNP A1, whereas only the bind of RdRp to PEDV 3'UTR was inhibited in AOFP3 treated cells.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Recombination between coronaviruses and synthetic RNAs and biorisk implications motivated by a SARS-CoV-2 FCS origin controversy

The urgent need for improved policy, regulation, and oversight of research with potential pandemic pathogens (PPPs) has been widely acknowledged. A 2022 article in Frontiers in Virology raises questions, reporting on a 100% sequence homology between the SARS-CoV-2 furin cleavage site (FCS) and the negative strand of a 2017 patented sequence. Even though Ambati and collaborators suspect a possible inadvertent or intentional cause leading to the FCS insert, the related underpinnings have not been studied from the perspective of potential biorisk policy gaps. A commentary on their article contests the low coincidence likelihood that was calculated by Ambati et al., arguing that the sequence match could have been a chance occurrence alone. Additionally, it has been suggested that the odds of the recombination event may be low. These considerations seem to have put many speculations related to any implied viral beginnings, notably from a research setting likely outside the Wuhan Institute of Virology, to rest. However, potential implications for future disasters in terms of biosafety and biosecurity have not been addressed. To demonstrate the feasibility of the Ambati et al. postulate, a theoretical framework is developed that substantially extends the research orientations implicated by these authors and the related patent. It is argued that specific experimental conditions, in combination, could significantly increase the implied recombination profile between coronaviruses and synthetic RNAs. Consequently, this article scrutinizes these largely unrecognized vulnerabilities to discuss implications across the spectrum of the biological risk landscape, with special attention to a potential "crime harvest." Focusing on insufficiently understood features of interaction between the natural and man-made world, vulnerabilities related to contaminants, camouflaging, and various misuse potentials fostered by the digitization and computerization of synthetic biology, it highlights novel biorisk gaps not covered by existing PPP policy. Even though this work does not aim to provide proof of the viral origin, it will make the point that, in theory, a convergence of under-appreciated lab experiments and technologies could have led to the SARS-CoV-2 FCS insert, which analogously could be exploited by various threat actors for the clandestine genesis of similar or even worse pathogens.

Mutagenesis and structural studies reveal the basis for the specific binding of SARS-CoV-2 SL3 RNA element with human TIA1 protein

Viral RNA-host protein interactions are indispensable during RNA virus transcription and replication, but their detailed structural and dynamical features remain largely elusive. Here, we characterize the binding interface for the SARS-CoV-2 stem-loop 3 (SL3) cis-acting element to human TIA1 protein with a combined theoretical and experimental approaches. The highly structured SARS-CoV-2 SL3 has a high binding affinity to TIA1 protein, in which the aromatic stacking, hydrogen bonds, and hydrophobic interactions collectively direct this specific binding. Further mutagenesis studies validate our proposed 3D binding model and reveal two SL3 variants have enhanced binding affinities to TIA1. And disruptions of the identified RNA-protein interactions with designed antisense oligonucleotides dramatically reduce SARS-CoV-2 infection in cells. Finally, TIA1 protein could interact with conserved SL3 RNA elements within other betacoronavirus lineages. These findings open an avenue to explore the viral RNA-host protein interactions and provide a pioneering structural basis for RNA-targeting antiviral drug design.

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

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

The preference signature of the SARS-CoV-2 Nucleocapsid NTD for its 5'-genomic RNA elements

The nucleocapsid protein (N) of SARS-CoV-2 plays a pivotal role during the viral life cycle. It is involved in RNA transcription and accounts for packaging of the large genome into virus particles. N manages the enigmatic balance of bulk RNA-coating versus precise RNA-binding to designated cis-regulatory elements. Numerous studies report the involvement of its disordered segments in non-selective RNA-recognition, but how N organizes the inevitable recognition of specific motifs remains unanswered. We here use NMR spectroscopy to systematically analyze the interactions of N's N-terminal RNA-binding domain (NTD) with individual cis RNA elements clustering in the SARS-CoV-2 regulatory 5'-genomic end. Supported by broad solution-based biophysical data, we unravel the NTD RNA-binding preferences in the natural genome context. We show that the domain's flexible regions read the intrinsic signature of preferred RNA elements for selective and stable complex formation within the large pool of available motifs.

Host microRNA interactions with the SARS-CoV-2 viral genome 3'-untranslated region

The 2019 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has marked the spread of a novel human coronavirus. While the viral life cycle is well understood, most of the interactions at the virus-host interface remain elusive. Furthermore, the molecular mechanisms behind disease severity and immune evasion are still largely unknown. Conserved elements of the viral genome such as secondary structures within the 5'- and 3'-untranslated regions (UTRs) serve as attractive targets of interest and could prove crucial in furthering our understanding of virus-host interactions. It has been proposed that microRNA (miR) interactions with viral components could be used by both the virus and host for their own benefit. Analysis of the SARS-CoV-2 viral genome 3'-UTR has revealed the potential for host cellular miR binding sites, providing sites for specific interactions with the virus. In this study, we demonstrate that the SARS-CoV-2 genome 3'-UTR binds the host cellular miRNAs miR-760-3p, miR-34a-5p, and miR-34b-5p, which have been shown to influence translation of interleukin-6 (IL-6), the IL-6 receptor (IL-6R), as well as progranulin (PGRN), respectively, proteins that have roles in the host immune response and inflammatory pathways. Furthermore, recent work suggests the potential of miR-34a-5p and miR-34b-5p to target and inhibit translation of viral proteins. Native gel electrophoresis and steady-state fluorescence spectroscopy were utilized to characterize the binding of these miRs to their predicted sites within the SARS-CoV-2 genome 3'-UTR. Additionally, we investigated 2'-fluoro-D-arabinonucleic acid (FANA) analogs of these miRNAs as competitive binding inhibitors for these miR binding interactions. The mechanisms detailed in this study have the potential to drive the development of antiviral treatments for SARS-CoV-2 infection, and provide a potential molecular basis for cytokine release syndrome and immune evasion which could implicate the host-virus interface.

The 5'UTR of HCoV-OC43 adopts a topologically constrained structure to intrinsically repress translation

The emergence of SARS-CoV-2, which is responsible for the COVID-19 pandemic, has highlighted the need for rapid characterization of viral mechanisms associated with cellular pathogenesis. Viral UTRs represent conserved genomic elements that contribute to such mechanisms. Structural details of most CoV UTRs are not available, however. Experimental approaches are needed to allow for the facile generation of high-quality viral RNA tertiary structural models, which can facilitate comparative mechanistic efforts. By integrating experimental and computational techniques, we herein report the efficient characterization of conserved RNA structures within the 5'UTR of the HCoV-OC43 genome, a lab-tractable model coronavirus. We provide evidence that the 5'UTR folds into a structure with well-defined stem-loops (SLs) as determined by chemical probing and direct detection of hydrogen bonds by NMR. We combine experimental base-pair restraints with global structural information from SAXS to generate a 3D model that reveals that SL1-4 adopts a topologically constrained structure wherein SLs 3 and 4 coaxially stack. Coaxial stacking is mediated by short linker nucleotides and allows SLs 1 to 2 to sample different cojoint orientations by pivoting about the SL3,4 helical axis. To evaluate the functional relevance of the SL3,4 coaxial helix, we engineered luciferase reporter constructs harboring the HCoV-OC43 5'UTR with mutations designed to abrogate coaxial stacking. Our results reveal that the SL3,4 helix intrinsically represses translation efficiency since the destabilizing mutations correlate with increased luciferase expression relative to wildtype without affecting reporter mRNA levels, thus highlighting how the 5'UTR structure contributes to the viral mechanism.

Deletion of the s2m RNA Structure in the Avian Coronavirus Infectious Bronchitis Virus and Human Astrovirus Results in Sequence Insertions

Coronaviruses infect a wide variety of host species, resulting in a range of diseases in both humans and animals. The coronavirus genome consists of a large positive-sense single-stranded molecule of RNA containing many RNA structures. One structure, denoted s2m and consisting of 41 nucleotides, is located within the 3' untranslated region (3' UTR) and is shared between some coronavirus species, including infectious bronchitis virus (IBV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2, as well as other pathogens, including human astrovirus. Using a reverse genetic system to generate recombinant viruses, we investigated the requirement of the s2m structure in the replication of IBV, a globally distributed economically important Gammacoronavirus that infects poultry causing respiratory disease. Deletion of three nucleotides predicted to destabilize the canonical structure of the s2m or the deletion of the nucleotides corresponding to s2m impacted viral replication in vitro. In vitro passaging of the recombinant IBV with the s2m sequence deleted resulted in a 36-nucleotide insertion in place of the deletion, which was identified to be composed of a duplication of flanking sequences. A similar result was observed following serial passage of human astrovirus with a deleted s2m sequence. RNA modeling indicated that deletion of the nucleotides corresponding to the s2m impacted other RNA structures present in the IBV 3' UTR. Our results indicated for both IBV and human astrovirus a preference for nucleotide occupation in the genome location corresponding to the s2m, which is independent of the specific s2m sequence. IMPORTANCE Coronaviruses infect many species, including humans and animals, with substantial effects on livestock, particularly with respect to poultry. The coronavirus RNA genome consists of structural elements involved in viral replication whose roles are poorly understood. We investigated the requirement of the RNA structural element s2m in the replication of the Gammacoronavirus infectious bronchitis virus, an economically important viral pathogen of poultry. Using reverse genetics to generate recombinant IBVs with either a disrupted or deleted s2m, we showed that the s2m is not required for viral replication in cell culture; however, replication is decreased in tracheal tissue, suggesting a role for the s2m in the natural host. Passaging of these viruses as well as human astrovirus lacking the s2m sequence demonstrated a preference for nucleotide occupation, independent of the s2m sequence. RNA modeling suggested deletion of the s2m may negatively impact other essential RNA structures.

Genome Structure, Life Cycle, and Taxonomy of Coronaviruses and the Evolution of SARS-CoV-2

Coronaviruses have a broad host range and exhibit high zoonotic potential. In this chapter, we describe their genomic organization in terms of encoded proteins and provide an introduction to the peculiar discontinuous transcription mechanism. Further, we present evolutionary conserved genomic RNA secondary structure features, which are involved in the complex replication mechanism. With a focus on computational methods, we review the emergence of SARS-CoV-2 starting with the 2019 strains. In that context, we also discuss the debated hypothesis of whether SARS-CoV-2 was created in a laboratory. We focus on the molecular evolution and the epidemiological dynamics of this recently emerged pathogen and we explain how variants of concern are detected and characterised. COVID-19, the disease caused by SARS-CoV-2, can spread through different transmission routes and also depends on a number of risk factors. We describe how current computational models of viral epidemiology, or more specifically, phylodynamics, have facilitated and will continue to enable a better understanding of the epidemic dynamics of SARS-CoV-2.

Antiviral Activity of Oligonucleotides Targeting the SARS-CoV-2 Genomic RNA Stem-Loop Sequences within the 3'-End of the ORF1b

Increased evidence shows vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exhibited no long-term efficacy and limited worldwide availability, while existing antivirals and treatment options have only limited efficacy. In this study, the main objective was the development of antiviral strategies using nucleic acid-based molecules. To this purpose, partially overlapped 6-19-mer phosphorothioate deoxyoligonucleotides (S-ONs) designed on the SARS-CoV-2 genomic RNA stem-loop packaging sequences within the 3' end of the ORF1b were synthetized using the direct and complementary sequence. Among the S-ONs tested, several oligonucleotides exhibited a fifty percent inhibitory concentration antiviral activity ranging from 0.27 to 34 μM, in the absence of cytotoxicity. The S-ON with a scrambled sequence used in the same conditions was not active. Moreover, selected 10-mer S-ONs were tested using different infectious doses and against different SARS-CoV-2 variants, showing comparable antiviral activity that was abrogated when the central sequence was mutated. Experiments to evaluate the intracellular functional target localization of the S-ON inhibitory activity were also performed. Collectively the data indicate that the SARS-CoV-2 packaging region in the 3' end of the ORF1b may be a promising target candidate for further investigation to develop innovative nucleic-acid-based antiviral therapy.

CovInter: interaction data between coronavirus RNAs and host proteins

Coronavirus has brought about three massive outbreaks in the past two decades. Each step of its life cycle invariably depends on the interactions among virus and host molecules. The interaction between virus RNA and host protein (IVRHP) is unique compared to other virus-host molecular interactions and represents not only an attempt by viruses to promote their translation/replication, but also the host's endeavor to combat viral pathogenicity. In other words, there is an urgent need to develop a database for providing such IVRHP data. In this study, a new database was therefore constructed to describe the interactions between coronavirus RNAs and host proteins (CovInter). This database is unique in (a) unambiguously characterizing the interactions between virus RNA and host protein, (b) comprehensively providing experimentally validated biological function for hundreds of host proteins key in viral infection and (c) systematically quantifying the differential expression patterns (before and after infection) of these key proteins. Given the devastating and persistent threat of coronaviruses, CovInter is highly expected to fill the gap in the whole process of the 'molecular arms race' between viruses and their hosts, which will then aid in the discovery of new antiviral therapies. It's now free and publicly accessible at: https://idrblab.org/covinter/.

Nanopore ReCappable sequencing maps SARS-CoV-2 5' capping sites and provides new insights into the structure of sgRNAs

The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested subgenomic RNAsused to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.

The Neighborhood of the Spike Gene Is a Hotspot for Modular Intertypic Homologous and Nonhomologous Recombination in Coronavirus Genomes

Coronaviruses (CoVs) have very large RNA viral genomes with a distinct genomic architecture of core and accessory open reading frames (ORFs). It is of utmost importance to understand their patterns and limits of homologous and nonhomologous recombination, because such events may affect the emergence of novel CoV strains, alter their host range, infection rate, tissue tropism pathogenicity, and their ability to escape vaccination programs. Intratypic recombination among closely related CoVs of the same subgenus has often been reported; however, the patterns and limits of genomic exchange between more distantly related CoV lineages (intertypic recombination) need further investigation. Here, we report computational/evolutionary analyses that clearly demonstrate a substantial ability for CoVs of different subgenera to recombine. Furthermore, we show that CoVs can obtain-through nonhomologous recombination-accessory ORFs from core ORFs, exchange accessory ORFs with different CoV genera, with other viruses (i.e., toroviruses, influenza C/D, reoviruses, rotaviruses, astroviruses) and even with hosts. Intriguingly, most of these radical events result from double crossovers surrounding the Spike ORF, thus highlighting both the instability and mobile nature of this genomic region. Although many such events have often occurred during the evolution of various CoVs, the genomic architecture of the relatively young SARS-CoV/SARS-CoV-2 lineage so far appears to be stable.

The Remarkable Evolutionary Plasticity of Coronaviruses by Mutation and Recombination: Insights for the COVID-19 Pandemic and the Future Evolutionary Paths of SARS-CoV-2

Coronaviruses (CoVs) constitute a large and diverse subfamily of positive-sense single-stranded RNA viruses. They are found in many mammals and birds and have great importance for the health of humans and farm animals. The current SARS-CoV-2 pandemic, as well as many previous epidemics in humans that were of zoonotic origin, highlights the importance of studying the evolution of the entire CoV subfamily in order to understand how novel strains emerge and which molecular processes affect their adaptation, transmissibility, host/tissue tropism, and patho non-homologous genicity. In this review, we focus on studies over the last two years that reveal the impact of point mutations, insertions/deletions, and intratypic/intertypic homologous and non-homologous recombination events on the evolution of CoVs. We discuss whether the next generations of CoV vaccines should be directed against other CoV proteins in addition to or instead of spike. Based on the observed patterns of molecular evolution for the entire subfamily, we discuss five scenarios for the future evolutionary path of SARS-CoV-2 and the COVID-19 pandemic. Finally, within this evolutionary context, we discuss the recently emerged Omicron (B.1.1.529) VoC.

Amilorides inhibit SARS-CoV-2 replication in vitro by targeting RNA structures

The SARS-CoV-2 pandemic, and the likelihood of future coronavirus pandemics, emphasized the urgent need for development of novel antivirals. Small-molecule chemical probes offer both to reveal aspects of virus replication and to serve as leads for antiviral therapeutic development. Here, we report on the identification of amiloride-based small molecules that potently inhibit OC43 and SARS-CoV-2 replication through targeting of conserved structured elements within the viral 5′-end. Nuclear magnetic resonance–based structural studies revealed specific amiloride interactions with stem loops containing bulge like structures and were predicted to be strongly bound by the lead amilorides in retrospective docking studies. Amilorides represent the first antiviral small molecules that target RNA structures within the 5′ untranslated regions and proximal region of the CoV genomes. These molecules will serve as chemical probes to further understand CoV RNA biology and can pave the way for the development of specific CoV RNA–targeted antivirals.

Metabolomics Signatures of SARS-CoV-2 Infection

For a very long time, viral infections have been considered as one of the most important causes of death and disability around the world. Through the viral infection, viruses as small pathogens enter the host cells and use hosts' biosynthesis machinery to replicate and collect infectious lineages. Moreover, they can modify hosts' metabolic pathways in order to their own purposes. Nowadays (in 2019-2020), the most famous type of viral infection which was caused by a novel type of coronavirus is called COVID-19 disease. It has claimed the lives of many people around the world and is a very serious threat to health. Since investigations of the effects of viruses on host metabolism using metabolomics tools may have given focuses on novel appropriate treatments, in the current review the authors highlighted the virus-host metabolic interactions and metabolomics perspective in COVID-19.

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.

In-Silico analysis reveals lower transcription efficiency of C241T variant of SARS-CoV-2 with host replication factors MADP1 and hnRNP-1

Novel severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has claimed more than 3.3 million lives worldwide and still counting. As per the GISAID database, the genomics of SARS-CoV-2 has been extensively studied, with more than 500 genome submissions per day. Out of several hotspot mutations within the SARS-CoV-2 genome, recent research has focused mainly on the missense variants. Moreover, significantly less attention has been accorded to delineate the role of the untranslated regions (UTRs) of the SARS-CoV-2 genome in the disease progression and etiology. One of the most frequent 5' UTR variants in the SARS-CoV-2 genome is the C241T, with a global frequency of more than 95 %. In the present study, the effect of the C241T mutation has been studied with respect to the changes in RNA structure and its interaction with the host replication factors MADP1 Zinc finger CCHC-type and RNA-binding motif 1 (hnRNP1). The results obtained from molecular docking and molecular dynamics simulation indicated weaker interaction of C241T mutant stem-loops with the host transcription factor MADP1, indicating a reduced replication efficiency. The results are also correlated with increased recovery rates and decreased death rates of global SARS-CoV-2 cases.

RNA-Protein Interaction Analysis of SARS-CoV-2 5' and 3' Untranslated Regions Reveals a Role of Lysosome-Associated Membrane Protein-2a during Viral Infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-strand RNA virus. The viral genome is capped at the 5′ end, followed by an untranslated region (UTR). There is a poly(A) tail at the 3′ end, preceded by a UTR. The self-interaction between the RNA regulatory elements present within the 5′ and 3′ UTRs and their interaction with host/virus-encoded proteins mediate the function of the 5′ and 3′ UTRs. Using an RNA-protein interaction detection (RaPID) assay coupled to liquid chromatography with tandem mass spectrometry, we identified host interaction partners of SARS-CoV-2 5′ and 3′ UTRs and generated an RNA-protein interaction network. By combining these data with the previously known protein-protein interaction data proposed to be involved in virus replication, we generated the RNA-protein-protein interaction (RPPI) network, likely to be essential for controlling SARS-CoV-2 replication. Notably, bioinformatics analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism. Lysosome-associated membrane protein-2a (Lamp2a), the receptor for chaperone-mediated autophagy, is one of the host proteins that interact with the 5′ UTR. Further studies showed that the Lamp2 level is upregulated in SARS-CoV-2-infected cells and that the absence of the Lamp2a isoform enhanced the viral RNA level whereas its overexpression significantly reduced the viral RNA level. Lamp2a and viral RNA colocalize in the infected cells, and there is an increased autophagic flux in infected cells, although there is no change in the formation of autophagolysosomes. In summary, our study provides a useful resource of SARS-CoV-2 5′ and 3′ UTR binding proteins and reveals the role of Lamp2a protein during SARS-CoV-2 infection.

IMPORTANCE Replication of a positive-strand RNA virus involves an RNA-protein complex consisting of viral genomic RNA, host RNA(s), virus-encoded proteins, and host proteins. Dissecting out individual components of the replication complex will help decode the mechanism of viral replication. 5′ and 3′ UTRs in positive-strand RNA viruses play essential regulatory roles in virus replication. Here, we identified the host proteins that associate with the UTRs of SARS-CoV-2, combined those data with the previously known protein-protein interaction data (expected to be involved in virus replication), and generated the RNA-protein-protein interaction (RPPI) network. Analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism, which are important for virus replication. Analysis of one of the interaction partners of the 5′-UTR (Lamp2a) demonstrated its role in reducing the viral RNA level in SARS-CoV-2-infected cells. Collectively, our study provides a resource of SARS-CoV-2 UTR-binding proteins and identifies an important role for host Lamp2a protein during viral infection.

Lantern: an integrative repository of functional annotations for lncRNAs in the human genome

BACKGROUND

With advancements in omics technologies, the range of biological processes where long non-coding RNAs (lncRNAs) are involved, is expanding extensively, thereby generating the need to develop lncRNA annotation resources. Although, there are a plethora of resources for annotating genes, despite the extensive corpus of lncRNA literature, the available resources with lncRNA ontology annotations are rare.

RESULTS

We present a lncRNA annotation extractor and repository (Lantern), developed using PubMed's abstract retrieval engine and NCBO's recommender annotation system. Lantern's annotations were benchmarked against lncRNAdb's manually curated free text. Benchmarking analysis suggested that Lantern has a recall of 0.62 against lncRNAdb for 182 lncRNAs and precision of 0.8. Additionally, we also annotated lncRNAs with multiple omics annotations, including predicted cis-regulatory TFs, interactions with RBPs, tissue-specific expression profiles, protein co-expression networks, coding potential, sub-cellular localization, and SNPs for ~ 11,000 lncRNAs in the human genome, providing a one-stop dynamic visualization platform.

CONCLUSIONS

Lantern integrates a novel, accurate semi-automatic ontology annotation engine derived annotations combined with a variety of multi-omics annotations for lncRNAs, to provide a central web resource for dissecting the functional dynamics of long non-coding RNAs and to facilitate future hypothesis-driven experiments. The annotation pipeline and a web resource with current annotations for human lncRNAs are freely available on sysbio.lab.iupui.edu/lantern.

The SARS-CoV-2 subgenome landscape and its novel regulatory features

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

SARS-CoV-2 main protease suppresses type I interferon production by preventing nuclear translocation of phosphorylated IRF3

Suppression of type I interferon (IFN) response is one pathological outcome of the infection of highly pathogenic human coronaviruses. To effect this, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 encode multiple IFN antagonists. In this study, we reported on the IFN antagonism of SARS-CoV-2 main protease NSP5. NSP5 proteins of both SARS-CoV and SARS-CoV-2 counteracted Sendai virus-induced IFN production. NSP5 variants G15S and K90R commonly seen in circulating strains of SARS-CoV-2 retained the IFN-antagonizing property. The suppressive effect of NSP5 on IFN-β gene transcription induced by RIG-I, MAVS, TBK1 and IKKϵ suggested that NSP5 likely acts at a step downstream of IRF3 phosphorylation in the cytoplasm. NSP5 did not influence steady-state expression or phosphorylation of IRF3, suggesting that IRF3, regardless of its phosphorylation state, might not be the substrate of NSP5 protease. However, nuclear translocation of phosphorylated IRF3 was severely compromised in NSP5-expressing cells. Taken together, our work revealed a new mechanism by which NSP5 proteins encoded by SARS-CoV and SARS-CoV-2 antagonize IFN production by retaining phosphorylated IRF3 in the cytoplasm. Our findings have implications in rational design and development of antiviral agents against SARS-CoV-2.

A novel end-to-end method to predict RNA secondary structure profile based on bidirectional LSTM and residual neural network

Background

Studies have shown that RNA secondary structure, a planar structure formed by paired bases, plays diverse vital roles in fundamental life activities and complex diseases. RNA secondary structure profile can record whether each base is paired with others. Hence, accurate prediction of secondary structure profile can help to deduce the secondary structure and binding site of RNA. RNA secondary structure profile can be obtained through biological experiment and calculation methods. Of them, the biological experiment method involves two ways: chemical reagent and biological crystallization. The chemical reagent method can obtain a large number of prediction data, but its cost is high and always associated with high noise, making it difficult to get results of all bases on RNA due to the limited of sequencing coverage. By contrast, the biological crystallization method can lead to accurate results, yet heavy experimental work and high costs are required. On the other hand, the calculation method is CROSS, which comprises a three-layer fully connected neural network. However, CROSS can not completely learn the features of RNA secondary structure profile since its poor network structure, leading to its low performance.

Results

In this paper, a novel end-to-end method, named as “RPRes, was proposed to predict RNA secondary structure profile based on Bidirectional LSTM and Residual Neural Network.

Conclusions

RPRes utilizes data sets generated by multiple biological experiment methods as the training, validation, and test sets to predict profile, which can compatible with numerous prediction requirements. Compared with the biological experiment method, RPRes has reduced the costs and improved the prediction efficiency. Compared with the state-of-the-art calculation method CROSS, RPRes has significantly improved performance.

Computational strategies to combat COVID-19: useful tools to accelerate SARS-CoV-2 and coronavirus research

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel virus of the family Coronaviridae. The virus causes the infectious disease COVID-19. The biology of coronaviruses has been studied for many years. However, bioinformatics tools designed explicitly for SARS-CoV-2 have only recently been developed as a rapid reaction to the need for fast detection, understanding and treatment of COVID-19. To control the ongoing COVID-19 pandemic, it is of utmost importance to get insight into the evolution and pathogenesis of the virus. In this review, we cover bioinformatics workflows and tools for the routine detection of SARS-CoV-2 infection, the reliable analysis of sequencing data, the tracking of the COVID-19 pandemic and evaluation of containment measures, the study of coronavirus evolution, the discovery of potential drug targets and development of therapeutic strategies. For each tool, we briefly describe its use case and how it advances research specifically for SARS-CoV-2. All tools are free to use and available online, either through web applications or public code repositories. Contact:evbc@unj-jena.de.

Novel perspectives for SARS-CoV-2 genome browsing

SARS-CoV-2 has spread worldwide and caused social, economic, and health turmoil. The first genome assembly of SARS-CoV-2 was produced in Wuhan, and it is widely used as a reference. Subsequently, more than a hundred additional SARS-CoV-2 genomes have been sequenced. While the genomes appear to be mostly identical, there are variations. Therefore, an alignment of all available genomes and the derived consensus sequence could be used as a reference, better serving the science community. Variations are significant, but representing them in a genome browser can become, especially if their sequences are largely identical. Here we summarize the variation in one track. Other information not currently found in genome browsers for SARS-CoV-2, such as predicted miRNAs and predicted TRS as well as secondary structure information, were also added as tracks to the consensus genome. We believe that a genome browser based on the consensus sequence is better suited when considering worldwide effects and can become a valuable resource in the combating of COVID-19. The genome browser is available at http://cov.iaba.online.

The global and local distribution of RNA structure throughout the SARS-CoV-2 genome

SARS-CoV-2 is the causative viral agent of COVID-19, the disease at the center of the current global pandemic. While knowledge of highly structured regions is integral for mechanistic insights into the viral infection cycle, very little is known about the location and folding stability of functional elements within the massive, ∼30kb SARS-CoV-2 RNA genome. In this study, we analyze the folding stability of this RNA genome relative to the structural landscape of other well-known viral RNAs. We present an in-silico pipeline to predict regions of high base pair content across long genomes and to pinpoint hotspots of well-defined RNA structures, a method that allows for direct comparisons of RNA structural complexity within the several domains in SARS-CoV-2 genome. We report that the SARS-CoV-2 genomic propensity for stable RNA folding is exceptional among RNA viruses, superseding even that of HCV, one of the most structured viral RNAs in nature. Furthermore, our analysis suggests varying levels of RNA structure across genomic functional regions, with accessory and structural ORFs containing the highest structural density in the viral genome. Finally, we take a step further to examine how individual RNA structures formed by these ORFs are affected by the differences in genomic and subgenomic contexts, which given the technical difficulty of experimentally separating cellular mixtures of sgRNA from gRNA, is a unique advantage of our in-silico pipeline. The resulting findings provide a useful roadmap for planning focused empirical studies of SARS-CoV-2 RNA biology, and a preliminary guide for exploring potential SARS-CoV-2 RNA drug targets.Importance The RNA genome of SARS-CoV-2 is among the largest and most complex viral genomes, and yet its RNA structural features remain relatively unexplored. Since RNA elements guide function in most RNA viruses, and they represent potential drug targets, it is essential to chart the architectural features of SARS-CoV-2 and pinpoint regions that merit focused study. Here we show that RNA folding stability of SARS-CoV-2 genome is exceptional among viral genomes and we develop a method to directly compare levels of predicted secondary structure across SARS-CoV-2 domains. Remarkably, we find that coding regions display the highest structural propensity in the genome, forming motifs that differ between the genomic and subgenomic contexts. Our approach provides an attractive strategy to rapidly screen for candidate structured regions based on base pairing potential and provides a readily interpretable roadmap to guide functional studies of RNA viruses and other pharmacologically relevant RNA transcripts.

Inhibition of anti-viral stress granule formation by coronavirus endoribonuclease nsp15 ensures efficient virus replication

Cytoplasmic stress granules (SGs) are generally triggered by stress-induced translation arrest for storing mRNAs. Recently, it has been shown that SGs exert anti-viral functions due to their involvement in protein synthesis shut off and recruitment of innate immune signaling intermediates. The largest RNA viruses, coronaviruses, impose great threat to public safety and animal health; however, the significance of SGs in coronavirus infection is largely unknown. Infectious Bronchitis Virus (IBV) is the first identified coronavirus in 1930s and has been prevalent in poultry farm for many years. In this study, we provided evidence that IBV overcomes the host antiviral response by inhibiting SGs formation via the virus-encoded endoribonuclease nsp15. By immunofluorescence analysis, we observed that IBV infection not only did not trigger SGs formation in approximately 80% of the infected cells, but also impaired the formation of SGs triggered by heat shock, sodium arsenite, or NaCl stimuli. We further demonstrated that the intrinsic endoribonuclease activity of nsp15 was responsible for the interference of SGs formation. In fact, nsp15-defective recombinant IBV (rIBV-nsp15-H238A) greatly induced the formation of SGs, along with accumulation of dsRNA and activation of PKR, whereas wild type IBV failed to do so. Consequently, infection with rIBV-nsp15-H238A strongly triggered transcription of IFN-β which in turn greatly affected rIBV-nsp15-H238A replication. Further analysis showed that SGs function as an antiviral hub, as demonstrated by the attenuated IRF3-IFN response and increased production of IBV in SG-defective cells. Additional evidence includes the aggregation of pattern recognition receptors (PRRs) and signaling intermediates to the IBV-induced SGs. Collectively, our data demonstrate that the endoribonuclease nsp15 of IBV interferes with the formation of antiviral hub SGs by regulating the accumulation of viral dsRNA and by antagonizing the activation of PKR, eventually ensuring productive virus replication. We further demonstrated that nsp15s from PEDV, TGEV, SARS-CoV, and SARS-CoV-2 harbor the conserved function to interfere with the formation of chemically-induced SGs. Thus, we speculate that coronaviruses employ similar nsp15-mediated mechanisms to antagonize the host anti-viral SGs formation to ensure efficient virus replication.

Coronavirus encodes the conserved endoribonuclease nsp15, which has been reported to antagonize IFN responses by mediating evasion of recognition by dsRNA sensors. SGs are part of the host cell anti-viral response; not surprisingly, viruses in turn produce an array of antagonists to counteract such host response. Here, we show that IBV prevents the formation of SGs via nsp15, by reducing the accumulation of viral dsRNA, thereby evading the activation of PKR, phosphorylation of eIF2α, and formation of SGs. Depletion of SG scaffold proteins G3BP1/2 decreases IRF3-IFN response and increases the production of IBV. When overexpressed alone, nsp15s from different coronaviruses (IBV, PEDV, TGEV, SARS-CoV, and SARS-CoV-2) interferes with chemically- and physically-induced SGs, probably by targeting essential SGs assembly factors. In this way, coronaviruses antagonize the formation of SGs by nsp15, via reducing the viral dsRNA accumulation and sequestering/depleting critical component of SGs. To our knowledge, this is the first report describing the role of coronavirus nsp15 in the suppression of integral stress response, in crosstalk with anti-innate immune response.

Rfam 14: expanded coverage of metagenomic, viral and microRNA families

Rfam is a database of RNA families where each of the 3444 families is represented by a multiple sequence alignment of known RNA sequences and a covariance model that can be used to search for additional members of the family. Recent developments have involved expert collaborations to improve the quality and coverage of Rfam data, focusing on microRNAs, viral and bacterial RNAs. We have completed the first phase of synchronising microRNA families in Rfam and miRBase, creating 356 new Rfam families and updating 40. We established a procedure for comprehensive annotation of viral RNA families starting with Flavivirus and Coronaviridae RNAs. We have also increased the coverage of bacterial and metagenome-based RNA families from the ZWD database. These developments have enabled a significant growth of the database, with the addition of 759 new families in Rfam 14. To facilitate further community contribution to Rfam, expert users are now able to build and submit new families using the newly developed Rfam Cloud family curation system. New Rfam website features include a new sequence similarity search powered by RNAcentral, as well as search and visualisation of families with pseudoknots. Rfam is freely available at https://rfam.org.

Structural Insight Into the SARS-CoV-2 Nucleocapsid Protein C-Terminal Domain Reveals a Novel Recognition Mechanism for Viral Transcriptional Regulatory Sequences

Coronavirus disease 2019 (COVID-19) has caused massive disruptions to society and the economy, and the transcriptional regulatory mechanisms behind the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are poorly understood. Herein, we determined the crystal structure of the SARS-CoV-2 nucleocapsid protein C-terminal domain (CTD) at a resolution of 2.0 Å, and demonstrated that the CTD has a comparable distinct electrostatic potential surface to equivalent domains of other reported CoVs, suggesting that the CTD has novel roles in viral RNA binding and transcriptional regulation. Further in vitro biochemical assays demonstrated that the viral genomic intergenic transcriptional regulatory sequences (TRSs) interact with the SARS-CoV-2 nucleocapsid protein CTD with a flanking region. The unpaired adeno dinucleotide in the TRS stem-loop structure is a major determining factor for their interactions. Taken together, these results suggested that the nucleocapsid protein CTD is responsible for the discontinuous viral transcription mechanism by recognizing the different patterns of viral TRS during transcription.

Translational control of coronaviruses

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

Amilorides inhibit SARS-CoV-2 replication in vitro by targeting RNA structures

The SARS-CoV-2 pandemic, and the likelihood of future coronavirus pandemics, has rendered our understanding of coronavirus biology more essential than ever. Small molecule chemical probes offer to both reveal novel aspects of virus replication and to serve as leads for antiviral therapeutic development. The RNA-biased amiloride scaffold was recently tuned to target a viral RNA structure critical for translation in enterovirus 71, ultimately uncovering a novel mechanism to modulate positive-sense RNA viral translation and replication. Analysis of CoV RNA genomes reveal many conserved RNA structures in the 5'-UTR and proximal region critical for viral translation and replication, including several containing bulge-like secondary structures suitable for small molecule targeting. Following phylogenetic conservation analysis of this region, we screened an amiloride-based small molecule library against a less virulent human coronavirus, OC43, to identify lead ligands. Amilorides inhibited OC43 replication as seen in viral plaque assays. Select amilorides also potently inhibited replication competent SARS-CoV-2 as evident in the decreased levels of cell free virions in cell culture supernatants of treated cells. Reporter screens confirmed the importance of RNA structures in the 5'-end of the viral genome for small molecule activity. Finally, NMR chemical shift perturbation studies of the first six stem loops of the 5'-end revealed specific amiloride interactions with stem loops 4, 5a, and 6, all of which contain bulge like structures and were predicted to be strongly bound by the lead amilorides in retrospective docking studies. Taken together, the use of multiple orthogonal approaches allowed us to identify the first small molecules aimed at targeting RNA structures within the 5'-UTR and proximal region of the CoV genome. These molecules will serve as chemical probes to further understand CoV RNA biology and can pave the way for the development of specific CoV RNA-targeted antivirals.

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.

SARS-CoV-2 Disrupts Splicing, Translation, and Protein Trafficking to Suppress Host Defenses

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a recently identified coronavirus that causes the respiratory disease known as coronavirus disease 2019 (COVID-19). Despite the urgent need, we still do not fully understand the molecular basis of SARS-CoV-2 pathogenesis. Here, we comprehensively define the interactions between SARS-CoV-2 proteins and human RNAs. NSP16 binds to the mRNA recognition domains of the U1 and U2 splicing RNAs and acts to suppress global mRNA splicing upon SARS-CoV-2 infection. NSP1 binds to 18S ribosomal RNA in the mRNA entry channel of the ribosome and leads to global inhibition of mRNA translation upon infection. Finally, NSP8 and NSP9 bind to the 7SL RNA in the signal recognition particle and interfere with protein trafficking to the cell membrane upon infection. Disruption of each of these essential cellular functions acts to suppress the interferon response to viral infection. Our results uncover a multipronged strategy utilized by SARS-CoV-2 to antagonize essential cellular processes to suppress host defenses.

•

NSP16 binds mRNA recognition domains of U1/U2 snRNAs and disrupts mRNA splicing

•

NSP1 binds in the mRNA entry channel of the ribosome to disrupt protein translation

•

NSP8 and NSP9 bind the signal recognition particle and disrupt protein trafficking

•

These disruptions of protein production suppress the interferon response to 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 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.

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.

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.

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.