Long-range RNA-RNA Interactions Research Articles

M6A methylation of transcription leader sequence impacts discontinuous transcription of subgenomic mRNAs.

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

Dengue Virus Capsid Protein Facilitates Genome Compaction and Packaging.

Dengue virus (DENV) is a single-stranded (+)-sense RNA virus that infects humans and mosquitoes, posing a significant health risk in tropical and subtropical regions. Mature virions are composed of an icosahedral shell of envelope (E) and membrane (M) proteins circumscribing a lipid bilayer, which in turn contains a complex of the approximately 11 kb genomic RNA with capsid (C) proteins. Whereas the structure of the envelope is clearly defined, the structure of the packaged genome in complex with C proteins remains elusive. Here, we investigated the interactions of C proteins with viral RNA, in solution and inside mature virions, via footprinting and cross-linking experiments. We demonstrated that C protein interaction with DENV genomes saturates at an RNA:C protein ratio below 1:250. Moreover, we also showed that the length of the RNA genome interaction sites varies, in a multimodal distribution, consistent with the C protein binding to each RNA site mostly in singlets or pairs (and, in some instances, higher numbers). We showed that interaction sites are preferentially sites with low base pairing, as previously measured by 2'-acetylation analyzed by primer extension (SHAPE) reactivity indicating structuredness. We found a clear association pattern emerged: RNA-C protein binding sites are strongly associated with long-range RNA-RNA interaction sites, particularly inside virions. This, in turn, explains the need for C protein in viral genome packaging: the protein has a chief role in coordinating these key interactions, promoting proper packaging of viral RNA. Such sites are, thus, highly consequential for viral assembly, and, as such, may be targeted in future drug development strategies against these and related viruses.

A potential long-range RNA-RNA interaction in the HIV-1 RNA

It is well-established that viral and cellular mRNAs alike harbour functional long-range intra-molecular RNA-RNA interactions. Despite the biological importance of such interactions, their identification and characterization remain challenging. Here we present a computational method for the identification of certain kinds of long-range intra-molecular RNA-RNA interactions involving the loop nucleotides of a hairpin loop. Using the computational method, we analysed 4272 HIV-1 genomic mRNAs. A potential long-range intra-molecular RNA-RNA interaction within the HIV-1 genomic RNA was identified. The long-range interaction is mediated by a kissing loop structure between two stem-loops of the previously reported SHAPE-based secondary structure of the entire HIV-1 genome. Structural modelling studies were carried out to show that the kissing loop structure not only is sterically feasible, but also contains a conserved RNA structural motif often found in compact RNA pseudoknots. The computational method should be generally applicable to the identification of potential long-range intra-molecular RNA-RNA interactions in any viral or cellular mRNA sequence. Communicated by Ramaswamy H. Sarma

Variant-Specific Analysis Reveals a Novel Long-Range RNA-RNA Interaction in SARS-CoV-2 Orf1a.

Since the start of the COVID-19 pandemic, understanding the pathology of the SARS-CoV-2 RNA virus and its life cycle has been the priority of many researchers. Currently, new variants of the virus have emerged with various levels of pathogenicity and abundance within the human-host population. Although much of viral pathogenicity is attributed to the viral Spike protein’s binding affinity to human lung cells’ ACE2 receptor, comprehensive knowledge on the distinctive features of viral variants that might affect their life cycle and pathogenicity is yet to be attained. Recent in vivo studies into the RNA structure of the SARS-CoV-2 genome have revealed certain long-range RNA-RNA interactions. Using in silico predictions and a large population of SARS-CoV-2 sequences, we observed variant-specific evolutionary changes for certain long-range RRIs. We also found statistical evidence for the existence of one of the thermodynamic-based RRI predictions, namely Comp1, in the Beta variant sequences. A similar test that disregarded sequence variant information did not, however, lead to significant results. When performing population-based analyses, aggregate tests may fail to identify novel interactions due to variant-specific changes. Variant-specific analyses can result in de novo RRI identification.

A Putative long-range RNA-RNA interaction between ORF8and Spike of SARS-CoV-2.

SARS-CoV-2 has affected people worldwide as the causative agent of COVID-19. The virus is related to the highly lethal SARS-CoV-1 responsible for the 2002-2003 SARS outbreak in Asia. Research is ongoing to understand why both viruses have different spreading capacities and mortality rates. Like other beta coronaviruses, RNA-RNA interactions occur between different parts of the viral genomic RNA, resulting in discontinuous transcription and production of various sub-genomic RNAs. These sub-genomic RNAs are then translated into other viral proteins. In this work, we performed a comparative analysis for novel long-range RNA-RNA interactions that may involve the Spike region. Comparing in-silico fragment-based predictions between reference sequences of SARS-CoV-1 and SARS-CoV-2 revealed several predictions amongst which a thermodynamically stable long-range RNA-RNA interaction between (23660-23703 Spike) and (28025-28060 ORF8) unique to SARS-CoV-2 was observed. The patterns of sequence variation using data gathered worldwide further supported the predicted stability of the sub-interacting region (23679-23690 Spike) and (28031-28042 ORF8). Such RNA-RNA interactions can potentially impact viral life cycle including sub-genomic RNA production rates.

The Short- and Long-Range RNA-RNA Interactome of SARS-CoV-2

SummaryThe Coronaviridae is a family of positive-strand RNA viruses that includes SARS-CoV-2, the etiologic agent of the COVID-19 pandemic. Bearing the largest single-stranded RNA genomes in nature, coronaviruses are critically dependent on long-distance RNA-RNA interactions to regulate the viral transcription and replication pathways. Here we experimentally mapped the in vivo RNA-RNA interactome of the full-length SARS-CoV-2 genome and subgenomic mRNAs. We uncovered a network of RNA-RNA interactions spanning tens of thousands of nucleotides. These interactions reveal that the viral genome and subgenomes adopt alternative topologies inside cells and engage in different interactions with host RNAs. Notably, we discovered a long-range RNA-RNA interaction, the FSE-arch, that encircles the programmed ribosomal frameshifting element. The FSE-arch is conserved in the related MERS-CoV and is under purifying selection. Our findings illuminate RNA structure-based mechanisms governing replication, discontinuous transcription, and translation of coronaviruses and will aid future efforts to develop antiviral strategies.

The Role of the RNA-RNA Interactome in the Hepatitis C Virus Life Cycle

RNA virus genomes are multifunctional entities endowed with conserved structural elements that control translation, replication and encapsidation, among other processes. The preservation of these structural RNA elements constraints the genomic sequence variability. The hepatitis C virus (HCV) genome is a positive, single-stranded RNA molecule with numerous conserved structural elements that manage different steps during the infection cycle. Their function is ensured by the association of protein factors, but also by the establishment of complex, active, long-range RNA-RNA interaction networks-the so-called HCV RNA interactome. This review describes the RNA genome functions mediated via RNA-RNA contacts, and revisits some canonical ideas regarding the role of functional high-order structures during the HCV infective cycle. By outlining the roles of long-range RNA-RNA interactions from translation to virion budding, and the functional domains involved, this work provides an overview of the HCV genome as a dynamic device that manages the course of viral infection.

The HCV genome domains 5BSL3.1 and 5BSL3.3 act as managers of translation

The RNA genome of the hepatitis C virus (HCV) encodes a single open reading frame (ORF) containing numerous functional elements. Among these, the cis-acting replication element (CRE) at the 3′ end of the viral ORF, has become of increasing interest given its dual role as a viral translation repressor and replication enhancer. Long-range RNA-RNA contacts mediated by the CRE build the structural scaffold required for its proper functioning. The recruitment of different cellular factors, many related to the functioning of the translation machinery, might aid in the CRE-exerted downregulation of viral translation. The present data show that the CRE promotes a defect in polysome production, and hinders the assembly of the 80S complex, likely through the direct, high affinity recruitment of the 40S ribosomal subunit. This interaction involves the highly conserved 5BSL3.1 and 5BSL3.3 domains of the CRE, and is strictly dependent on RNA-protein contacts, particularly with the ribosomal proteins RPSA and RPS29. These observations support a model in which the CRE-mediated inhibition of viral translation is a multifactorial process defined by the establishment of long-range RNA-RNA interactions between the 5′ and 3′ ends of the viral genome, the sequestration of the 40S subunit by the CRE, and the subsequent stalling of polysome elongation at the 3′ end of the ORF, all governed by the highly stable hairpin domains 5BSL3.1 and 5BSL3.3. The present data thus suggest a new managerial role in HCV translation for these 5BSL3.1 and 5BSL3.3 domains.

Form confers function: Case of the 3'X region of the hepatitis C virus genome.

At the 3’ end of genomic hepatitis C virus (HCV) RNA there is a highly conserved untranslated region, the 3’X-tail, which forms part of the 3’UTR. This region plays key functions in regulation of critical processes of the viral life cycle. The 3’X region is essential for viral replication and infectivity. It is also responsible for regulation of switching between translation and transcription of the viral RNA. There is some evidence indicating the contribution of the 3’X region to the translation efficiency of the viral polyprotein and to the encapsidation process. Several different secondary structure models of the 3’X region, based on computer predictions and experimental structure probing, have been proposed. It is likely that the 3’X region adopts more than one structural form in infected cells and that a specific equilibrium between the various forms regulates several aspects of the viral life cycle. The most intriguing explanations of the structural heterogeneity problem of the 3’X region came with the discovery of its involvement in long-range RNA-RNA interactions and the potential for homodimer formation. This article summarizes current knowledge on the structure and function of the 3’X region of hepatitis C genomic RNA, reviews previous opinions, presents new hypotheses and summarizes the questions that still remain unanswered.

A Stem-Loop Structure in Potato Leafroll Virus Open Reading Frame 5 (ORF5) Is Essential for Readthrough Translation of the Coat Protein ORF Stop Codon 700 Bases Upstream.

Translational readthrough of the stop codon of the capsid protein (CP) open reading frame (ORF) is used by members of the Luteoviridae to produce their minor capsid protein as a readthrough protein (RTP). The elements regulating RTP expression are not well understood, but they involve long-distance interactions between RNA domains. Using high-resolution mass spectrometry, glutamine and tyrosine were identified as the primary amino acids inserted at the stop codon of Potato leafroll virus (PLRV) CP ORF. We characterized the contributions of a cytidine-rich domain immediately downstream and a branched stem-loop structure 600 to 700 nucleotides downstream of the CP stop codon. Mutations predicted to disrupt and restore the base of the distal stem-loop structure prevented and restored stop codon readthrough. Motifs in the downstream readthrough element (DRTE) are predicted to base pair to a site within 27 nucleotides (nt) of the CP ORF stop codon. Consistent with a requirement for this base pairing, the DRTE of Cereal yellow dwarf virus was not compatible with the stop codon-proximal element of PLRV in facilitating readthrough. Moreover, deletion of the complementary tract of bases from the stop codon-proximal region or the DRTE of PLRV prevented readthrough. In contrast, the distance and sequence composition between the two domains was flexible. Mutants deficient in RTP translation moved long distances in plants, but fewer infection foci developed in systemically infected leaves. Selective 2'-hydroxyl acylation and primer extension (SHAPE) probing to determine the secondary structure of the mutant DRTEs revealed that the functional mutants were more likely to have bases accessible for long-distance base pairing than the nonfunctional mutants. This study reveals a heretofore unknown combination of RNA structure and sequence that reduces stop codon efficiency, allowing translation of a key viral protein.IMPORTANCE Programmed stop codon readthrough is used by many animal and plant viruses to produce key viral proteins. Moreover, such "leaky" stop codons are used in host mRNAs or can arise from mutations that cause genetic disease. Thus, it is important to understand the mechanism(s) of stop codon readthrough. Here, we shed light on the mechanism of readthrough of the stop codon of the coat protein ORFs of viruses in the Luteoviridae by identifying the amino acids inserted at the stop codon and RNA structures that facilitate this "leakiness" of the stop codon. Members of the Luteoviridae encode a C-terminal extension to the capsid protein known as the readthrough protein (RTP). We characterized two RNA domains in Potato leafroll virus (PLRV), located 600 to 700 nucleotides apart, that are essential for efficient RTP translation. We further determined that the PLRV readthrough process involves both local structures and long-range RNA-RNA interactions. Genetic manipulation of the RNA structure altered the ability of PLRV to translate RTP and systemically infect the plant. This demonstrates that plant virus RNA contains multiple layers of information beyond the primary sequence and extends our understanding of stop codon readthrough. Strategic targets that can be exploited to disrupt the virus life cycle and reduce its ability to move within and between plant hosts were revealed.

Signals Involved in Regulation of Hepatitis C Virus RNA Genome Translation and Replication.

Hepatitis C virus (HCV) preferentially replicates in the human liver and frequently causes chronic infection, often leading to cirrhosis and liver cancer. HCV is an enveloped virus classified in the genus Hepacivirus in the family Flaviviridae and has a single-stranded RNA genome of positive orientation. The HCV RNA genome is translated and replicated in the cytoplasm. Translation is controlled by the Internal Ribosome Entry Site (IRES) in the 5′ untranslated region (5′ UTR), while also downstream elements like the cis-replication element (CRE) in the coding region and the 3′ UTR are involved in translation regulation. The cis-elements controlling replication of the viral RNA genome are located mainly in the 5′- and 3′-UTRs at the genome ends but also in the protein coding region, and in part these signals overlap with the signals controlling RNA translation. Many long-range RNA–RNA interactions (LRIs) are predicted between different regions of the HCV RNA genome, and several such LRIs are actually involved in HCV translation and replication regulation. A number of RNA cis-elements recruit cellular RNA-binding proteins that are involved in the regulation of HCV translation and replication. In addition, the liver-specific microRNA-122 (miR-122) binds to two target sites at the 5′ end of the viral RNA genome as well as to at least three additional target sites in the coding region and the 3′ UTR. It is involved in the regulation of HCV RNA stability, translation and replication, thereby largely contributing to the hepatotropism of HCV. However, we are still far from completely understanding all interactions that regulate HCV RNA genome translation, stability, replication and encapsidation. In particular, many conclusions on the function of cis-elements in HCV replication have been obtained using full-length HCV genomes or near-full-length replicon systems. These include both genome ends, making it difficult to decide if a cis-element in question acts on HCV replication when physically present in the plus strand genome or in the minus strand antigenome. Therefore, it may be required to use reduced systems that selectively focus on the analysis of HCV minus strand initiation and/or plus strand initiation.

The 5BSL3.2 Functional RNA Domain Connects Distant Regions in the Hepatitis C Virus Genome.

Viral genomes are complexly folded entities that carry all the information required for the infective cycle. The nucleotide sequence of the RNA virus genome encodes proteins and functional information contained in discrete, highly conserved structural units. These so-called functional RNA domains play essential roles in the progression of infection, which requires their preservation from one generation to the next. Numerous functional RNA domains exist in the genome of the hepatitis C virus (HCV). Among them, the 5BSL3.2 domain in the cis-acting replication element (CRE) at the 3′ end of the viral open reading frame has become of particular interest given its role in HCV RNA replication and as a regulator of viral protein synthesis. These functionalities are achieved via the establishment of a complex network of long-distance RNA–RNA contacts involving (at least as known to date) the highly conserved 3′X tail, the apical loop of domain IIId in the internal ribosome entry site, and/or the so-called Alt region upstream of the CRE. Changing contacts promotes the execution of different stages of the viral cycle. The 5BSL3.2 domain thus operates at the core of a system that governs the progression of HCV infection. This review summarizes our knowledge of the long-range RNA–RNA interaction network in the HCV genome, with special attention paid to the structural and functional consequences derived from the establishment of different contacts. The potential implications of such interactions in switching between the different stages of the viral cycle are discussed.

Prediction of conserved long-range RNA-RNA interactions in full viral genomes.

Motivation: Long-range RNA-RNA interactions (LRIs) play an important role in viral replication, however, only a few of these interactions are known and only for a small number of viral species. Up to now, it has been impossible to screen a full viral genome for LRIs experimentally or in silico. Most known LRIs are cross-reacting structures (pseudoknots) undetectable by most bioinformatical tools. Results: We present LRIscan, a tool for the LRI prediction in full viral genomes based on a multiple genome alignment. We confirmed 14 out of 16 experimentally known and evolutionary conserved LRIs in genome alignments of HCV, Tombusviruses, Flaviviruses and HIV-1. We provide several promising new interactions, which include compensatory mutations and are highly conserved in all considered viral sequences. Furthermore, we provide reactivity plots highlighting the hot spots of predicted LRIs. Availability and Implementation: Source code and binaries of LRIscan freely available for download at http://www.rna.uni-jena.de/en/supplements/lriscan/, implemented in Ruby/C ++ and supported on Linux and Windows. Contact: manja@uni-jena.de Supplementary information: Supplementary data are available at Bioinformatics online.

RNA Duplex Map in Living Cells Reveals Higher-Order Transcriptome Structure

RNA has the intrinsic property to base pair, forming complex structures fundamental to its diverse functions. Here, we develop PARIS, a method based on reversible psoralen crosslinking for global mapping of RNA duplexes with near base-pair resolution in living cells. PARIS analysis in three human and mouse cell types reveals frequent long-range structures, higher-order architectures, and RNA-RNA interactions in trans across the transcriptome. PARIS determines base-pairing interactions on an individual-molecule level, revealing pervasive alternative conformations. We used PARIS-determined helices to guide phylogenetic analysis of RNA structures and discovered conserved long-range and alternative structures. XIST, a long noncoding RNA (lncRNA) essential for X chromosome inactivation, folds into evolutionarily conserved RNA structural domains that span many kilobases. XIST A-repeat forms complex inter-repeat duplexes that nucleate higher-order assembly of the key epigenetic silencing protein SPEN. PARIS is a generally applicable and versatile method that provides novel insights into the RNA structurome and interactome. VIDEO ABSTRACT.

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.

Overlapping local and long-range RNA-RNA interactions modulate dengue virus genome cyclization and replication.

The dengue virus genome is a dynamic molecule that adopts different conformations in the infected cell. Here, using RNA folding predictions, chemical probing analysis, RNA binding assays, and functional studies, we identified new cis-acting elements present in the capsid coding sequence that facilitate cyclization of the viral RNA by hybridization with a sequence involved in a local dumbbell structure at the viral 3' untranslated region (UTR). The identified interaction differentially enhances viral replication in mosquito and mammalian cells.

Functional long-range RNA-RNA interactions in positive-strand RNA viruses.

Positive-strand RNA viruses are important human, animal and plant pathogens that are defined by their single-stranded positive-sense RNA genomes. In recent years, it has become increasingly evident that interactions that occur between distantly positioned RNA sequences within these genomes can mediate important viral activities. These long-range intragenomic RNA-RNA interactions involve direct nucleotide base pairing and can span distances of thousands of nucleotides. In this Review, we discuss recent insights into the structure and function of these intriguing genomic features and highlight their diverse roles in the gene expression and genome replication of positive-strand RNA viruses.

Translational readthrough in Tobacco necrosis virus-D

The plus-strand RNA genome of Tobacco necrosis virus-D (TNV-D) expresses its polymerase via translational readthrough. The RNA signals involved in this readthrough process were characterized in vitro using a wheat germ extract translation system and in vivo via protoplast infections. The results indicate that (i) TNV-D requires a long-range RNA-RNA interaction between an extended stem-loop (SL) structure proximal to the readthrough site and a sequence in the 3'-untranslated region of its genome; (ii) stability of the extended SL structure is important for its function; (iii) TNV-D readthrough elements are compatible with UAG and UGA, but not UAA; (iv) a readthrough defect can be rescued by a heterologous readthrough element in vitro, but not in vivo; and (v) readthrough elements can also mediate translational frameshifting. These results provide new information on determinants of readthrough in TNV-D and further support the concept of a common general mechanism for readthrough in Tombusviridae.

Norovirus Genome Circularization and Efficient Replication Are Facilitated by Binding of PCBP2 and hnRNP A1

Sequences and structures within the terminal genomic regions of plus-strand RNA viruses are targets for the binding of host proteins that modulate functions such as translation, RNA replication, and encapsidation. Using murine norovirus 1 (MNV-1), we describe the presence of long-range RNA-RNA interactions that were stabilized by cellular proteins. The proteins potentially responsible for the stabilization were selected based on their ability to bind the MNV-1 genome and/or having been reported to be involved in the stabilization of RNA-RNA interactions. Cell extracts were preincubated with antibodies against the selected proteins and used for coprecipitation reactions. Extracts treated with antibodies to poly(C) binding protein 2 (PCBP2) and heterogeneous nuclear ribonucleoprotein (hnRNP) A1 significantly reduced the 5'-3' interaction. Both PCBP2 and hnRNP A1 recombinant proteins stabilized the 5'-3' interactions and formed ribonucleoprotein complexes with the 5' and 3' ends of the MNV-1 genomic RNA. Mutations within the 3' complementary sequences (CS) that disrupt the 5'-3'-end interactions resulted in a significant reduction of the viral titer, suggesting that the integrity of the 3'-end sequence and/or the lack of complementarity with the 5' end is important for efficient virus replication. Small interfering RNA-mediated knockdown of PCBP2 or hnRNP A1 resulted in a reduction in virus yield, confirming a role for the observed interactions in efficient viral replication. PCBP2 and hnRNP A1 induced the circularization of MNV-1 RNA, as revealed by electron microscopy. This study provides evidence that PCBP2 and hnRNP A1 bind to the 5' and 3' ends of the MNV-1 viral RNA and contribute to RNA circularization, playing a role in the virus life cycle.

Computational identification of biologically functional non-hairpin GC-helices in human Argonaute mRNA

BackgroundPerfectly formed duplex elements in RNA occur within folding units, often as a part of hairpin motifs which can be reliably predicted by various RNA folding algorithms. Double helices with consecutive Watson-Crick base-pairing may also be formed between distant RNA segments thereby facilitating long-range interactions of long-chain RNA that may be biologically functional. Here we addressed the potential formation of RNA duplex motifs by long-range RNA-RNA interactions of distantly located matching sequence elements of a single long-chain RNA.ResultsWe generated a Python-based software tool that identifies consecutive RNA duplex elements at any given length and nucleotide content formed by distant sequences. The software tool, dubbed RNAslider, is built on the theoretical RNA structure prediction algorithm Mfold. Source code and sample data sets are available on demand. We found that a small ratio of human genes including the Argonaute (Ago)-like gene family encode mRNAs containing highly GC-rich non-hairpin duplex elements (GC-helix) of equal to or more than 8 base pairs in length and we provide experimental evidence for their biological significance.ConclusionGC-helices are observed preferentially within the 5′-region of mRNAs in an evolutionarily conserved fashion indicating their potential biological role. This view is supported experimentally by post-transcriptional regulation of gene expression of a fusion transcript containing 5′-sequences of human mRNAAgo2 harbouring GC-helices and down-stream coding sequences of Renilla luciferase.