Modulation of RNA Condensation by the DEAD-Box Protein eIF4A.

Stress granules are condensates of non-translating mRNAs and proteins involved in the stress response and neurodegenerative diseases. Stress granules form in part through intermolecular RNA-RNA interactions, and to better understand how RNA-based condensation occurs, we demonstrate that RNA is effectively recruited to the surfaces of RNA or RNP condensates in vitro. We demonstrate that, through ATP-dependent RNA binding, the DEAD-box protein eIF4A reduces RNA condensation in vitro and limits stress granule formation in cells. This defines a function for eIF4A to limit intermolecular RNA-RNA interactions in cells. These results establish an important role for eIF4A, and potentially other DEAD-box proteins, as ATP-dependent RNA chaperones that limit the condensation of RNA, analogous to the function of proteins like HSP70 in combatting protein aggregates.

SARS-CoV-2 RNA stabilizes host mRNAs to elicit immunopathogenesis.

SARS-CoV-2 RNA interacts with host factors to suppress interferon responses and simultaneously induces cytokine release to drive the development of severe coronavirus disease 2019 (COVID-19). However, how SARS-CoV-2 hijacks host RNAs to elicit such imbalanced immune responses remains elusive. Here, we analyzed SARS-CoV-2 RNA in situ structures and interactions in infected cells and patient lung samples using RIC-seq. We discovered that SARS-CoV-2 RNA forms 2,095 potential duplexes with the 3' UTRs of 205 host mRNAs to increase their stability by recruiting RNA-binding protein YBX3 in A549 cells. Disrupting the SARS-CoV-2-to-host RNA duplex or knocking down YBX3 decreased host mRNA stability and reduced viral replication. Among SARS-CoV-2-stabilized host targets, NFKBIZ was crucial for promoting cytokine production and reducing interferon responses, probably contributing to cytokine storm induction. Our study uncovers the crucial roles of RNA-RNA interactions in the immunopathogenesis of RNA viruses such as SARS-CoV-2 and provides valuable host targets for drug development.

Circular RNA ZNF609/CKAP5 mRNA interaction regulates microtubule dynamics and tumorigenicity.

Circular RNAs (circRNAs) are widely expressed in eukaryotes and are regulated in many biological processes. Although several studies indicate their activity as microRNA (miRNA) and protein sponges, little is known about their ability to directly control mRNA homeostasis. We show that the widely expressed circZNF609 directly interacts with several mRNAs and increases their stability and/or translation by favoring the recruitment of the RNA-binding protein ELAVL1. Particularly, the interaction with CKAP5 mRNA, which interestingly overlaps the back-splicing junction, enhances CKAP5 translation, regulating microtubule function in cancer cells and sustaining cell-cycle progression. Finally, we show that circZNF609 downregulation increases the sensitivity of several cancer cell lines to different microtubule-targeting chemotherapeutic drugs and that locked nucleic acid (LNA) protectors against the pairing region on circZNF609 phenocopy such effects. These data set an example of how the small effects tuned by circZNF609/CKAP5 mRNA interaction might have a potent output in tumor growth and drug response.

Unique role for lncRNA HOTAIR in defining depot-specific gene expression patterns in human adipose-derived stem cells.

Accumulation of fat above the waist is an important risk factor in developing obesity-related comorbidities independently of BMI or total fat mass. Deciphering the gene regulatory programs of the adipose tissue precursor cells within upper body or abdominal (ABD) and lower body or gluteofemoral (GF) depots is important to understand their differential capacity for lipid accumulation, maturation, and disease risk. Previous studies identified the HOX transcript antisense intergenic RNA (HOTAIR) as a GF-specific lncRNA; however, its role in adipose tissue biology is still unclear. Using three different approaches (silencing of HOTAIR in GF human adipose-derived stem cells [GF hASCs], overexpression of HOTAIR in ABD hASCs, and ChIRP-seq) to localize HOTAIR binding in GF hASC chromatin, we found that HOTAIR binds and modulates expression, both positively and negatively, of genes involved in adipose tissue-specific pathways, including adipogenesis. We further demonstrate a direct interaction between HOTAIR and genes with high RNAPII binding in their gene bodies, especially at their 3' ends or transcription end sites. Computational analysis suggests HOTAIR binds preferentially to the 3' ends of genes containing predicted strong RNA-RNA interactions with HOTAIR. Together, these results reveal a unique function for HOTAIR in hASC depot-specific regulation of gene expression.

Guide RNA acrobatics: the one-for-two shuffle.

RNA modifications are crucial for the proper function of the RNAs. The sites of pseudouridines are often specified by dual hairpin guide RNAs, with one or both hairpins identifying a target uridine. In this issue of Genes & Development, Jády and colleagues (pp. 70-83) identify a novel mechanism by which a single guide RNA hairpin can specify two uridines adjacent to each other or separated by 1 nt; i.e., one for two or guide RNA acrobatics.

Guide RNA acrobatics: positioning consecutive uridines for pseudouridylation by H/ACA pseudouridylation loops with dual guide capacity.

Site-specific pseudouridylation of human ribosomal and spliceosomal RNAs is directed by H/ACA guide RNAs composed of two hairpins carrying internal pseudouridylation guide loops. The distal "antisense" sequences of the pseudouridylation loop base-pair with the target RNA to position two unpaired target nucleotides 5'-UN-3', including the 5' substrate U, under the base of the distal stem topping the guide loop. Therefore, each pseudouridylation loop is expected to direct synthesis of a single pseudouridine (Ψ) in the target sequence. However, in this study, genetic depletion and restoration and RNA mutational analyses demonstrate that at least four human H/ACA RNAs (SNORA53, SNORA57, SCARNA8, and SCARNA1) carry pseudouridylation loops supporting efficient and specific synthesis of two consecutive pseudouridines (ΨΨ or ΨNΨ) in the 28S (Ψ3747/Ψ3749), 18S (Ψ1045/Ψ1046), and U2 (Ψ43/Ψ44 and Ψ89/Ψ91) RNAs, respectively. In order to position two substrate Us for pseudouridylation, the dual guide loops form alternative base-pairing interactions with their target RNAs. This remarkable structural flexibility of dual pseudouridylation loops provides an unexpected versatility for RNA-directed pseudouridylation without compromising its efficiency and accuracy. Besides supporting synthesis of at least 6% of human ribosomal and spliceosomal Ψs, evidence indicates that dual pseudouridylation loops also participate in pseudouridylation of yeast and archaeal rRNAs.

Stress-Induced Translation Inhibition through Rapid Displacement of Scanning Initiation Factors.

Cellular responses to environmental stress are frequently mediated by RNA-binding proteins (RBPs). Here, we examined global RBP dynamics in Saccharomyces cerevisiae in response to glucose starvation and heat shock. Each stress induced rapid remodeling of the RNA-protein interactome without corresponding changes in RBP abundance. Consistent with general translation shutdown, ribosomal proteins contacting the mRNA showed decreased RNA association. Among translation components, RNA association was most reduced for initiation factors involved in 40S scanning (eukaryotic initiation factor 4A [eIF4A], eIF4B, and Ded1), indicating a common mechanism of translational repression. In unstressed cells, eIF4A, eIF4B, and Ded1 primarily targeted the 5' ends of mRNAs. Following glucose withdrawal, 5' binding was abolished within 30 s, explaining the rapid translation shutdown, but mRNAs remained stable. Heat shock induced progressive loss of 5' RNA binding by initiation factors over ∼16 min and provoked mRNA degradation, particularly for translation-related factors, mediated by Xrn1. Taken together, these results reveal mechanisms underlying translational control of gene expression during stress.

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

The 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.

Loss of Dynamic RNA Interaction and Aberrant Phase Separation Induced by Two Distinct Types of ALS/FTD-Linked FUS Mutations.

FUS is a nuclear RNA-binding protein, and its cytoplasmic aggregation is a pathogenic signature of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). It remains unknown how the FUS-RNA interactions contribute to phase separation and whether its phase behavior is affected by ALS-linked mutations. Here we demonstrate that wild-type FUS binds single-stranded RNA stoichiometrically in a length-dependent manner and that multimers induce highly dynamic interactions with RNA, giving rise to small and fluid condensates. In contrast, mutations in arginine display a severely altered conformation, static binding to RNA, and formation of large condensates, signifying the role of arginine in driving proper RNA interaction. Glycine mutations undergo rapid loss of fluidity, emphasizing the role of glycine in promoting fluidity. Strikingly, the nuclear import receptor Karyopherin-ß2 reverses the mutant defects and recovers the wild-type FUS behavior. We reveal two distinct mechanisms underpinning potentially disparate pathogenic pathways of ALS-linked FUS mutants.

System-wide Profiling of RNA-Binding Proteins Uncovers Key Regulators of Virus Infection.

The compendium of RNA-binding proteins (RBPs) has been greatly expanded by the development of RNA-interactome capture (RIC). However, it remained unknown if the complement of RBPs changes in response to environmental perturbations and whether these rearrangements are important. To answer these questions, we developed "comparative RIC" and applied it to cells challenged with an RNA virus called sindbis (SINV). Over 200 RBPs display differential interaction with RNA upon SINV infection. These alterations are mainly driven by the loss of cellular mRNAs and the emergence of viral RNA. RBPs stimulated by the infection redistribute to viral replication factories and regulate the capacity of the virus to infect. For example, ablation of XRN1 causes cells to be refractory to SINV, while GEMIN5 moonlights as a regulator of SINV gene expression. In summary, RNA availability controls RBP localization and function in SINV-infected cells.

3D feasibility of 2D RNA-RNA interaction paths by stepwise folding simulations.

The structure of an RNA, and even more so its interactions with other RNAs, provide valuable information about its function. Secondary structure-based tools for RNA-RNA interaction predictions provide a quick way to identify possible interaction targets and structures. However, these tools ignore the effect of steric hindrance on the tertiary (3D) structure level, and do not consider whether a suitable folding pathway exists to form the interaction. As a consequence, these tools often predict interactions that are unrealistically long and could be formed (in three dimensions) only by going through highly entangled intermediates. Here, we present a computational pipeline to assess whether a proposed secondary (2D) structure interaction is sterically feasible and reachable along a plausible folding pathway. To this end, we simulate the folding of a series of 3D structures along a given 2D folding path. To avoid the complexity of large-scale atomic resolution simulations, our pipeline uses coarse-grained 3D modeling and breaks up the folding path into small steps, each corresponding to the extension of the interaction by 1 or 2 bp. We apply our pipeline to analyze RNA-RNA interaction formation for three selected RNA-RNA complexes. We find that kissing hairpins, in contrast to interactions in the exterior loop, are difficult to extend and tend to get stuck at an interaction length of 6 bp. Our tool, including source code, documentation, and sample data, is available at www.github.com/irenekb/RRI-3D.

Distinct Gag interaction properties of HIV-1 RNA 5' leader conformers reveal a mechanism for dimeric genome selection.

During HIV-1 assembly, two copies of viral genomic RNAs (gRNAs) are selectively packaged into new viral particles. This process is mediated by specific interactions between HIV-1 Gag and the packaging signals at the 5' leader (5'L) of viral gRNA. 5'L is able to adopt different conformations, which promotes either gRNA dimerization and packaging or Gag translation. Dimerization and packaging are coupled. Whether the selective packaging of the gRNA dimer is due to favorable interactions between Gag and 5'L in the packaging conformation is not known. Here, using RNAs mimicking the two 5'L conformers, we show that the 5'L conformation dramatically affects Gag-RNA interactions. Compared to the RNA in the translation conformation (5'LT), the RNA in the packaging conformation (5'LP) can bind more Gag molecules. Gag associates with 5'LP faster than it binds to 5'LT, whereas Gag dissociates from 5'LP more slowly. The Gag-5'LP complex is more stable at high salt concentrations. The NC-SP2-p6 region of Gag likely accounts for the faster association and slower dissociation kinetics for the Gag-5'LP interaction and for the higher stability. In summary, our data suggest that conformational changes play an important role in the selection of dimeric genomes, probably by affecting the binding kinetics and stability of the Gag-5'L complex.

Double-stemmed and split structural variants of fluorescent RNA Mango aptamers.

Aptamers with fluorogenic ligands are emerging as useful tools to quantify and track RNA molecules. The RNA Mango family of aptamers have a useful combination of tight ligand binding, bright fluorescence, and small size. However, the simple structure of these aptamers, with a single base-paired stem capped by a G-quadruplex, can limit the sequence and structural modifications needed for many use-inspired designs. Here we report new structural variants of RNA Mango that have two base-paired stems attached to the quadruplex. Fluorescence saturation analysis of one of the double-stemmed constructs showed a maximum fluorescence that is ∼75% brighter than the original single-stemmed Mango I. A small number of mutations to nucleotides in the tetraloop-like linker of the second stem were subsequently analyzed. The effect of these mutations on the affinity and fluorescence suggested that the nucleobases of the second linker do not directly interact with the fluorogenic ligand (TO1-biotin), but may instead induce higher fluorescence by indirectly altering the ligand properties in the bound state. The effects of the mutations in this second tetraloop-like linker indicate the potential of this second stem for rational design and reselection experiments. Additionally, we demonstrated that a bimolecular mango designed by splitting the double-stemmed Mango can function when two RNA molecules are cotranscribed from different DNA templates in a single in vitro transcription. This bimolecular Mango has potential application in detecting RNA-RNA interactions. Together, these constructs expand the designability of the Mango aptamers to facilitate future applications of RNA imaging.

RNA in situ conformation sequencing reveals novel long-range RNA structures with impact on splicing.

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability.

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools.

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

Mechanisms and Regulation of RNA Condensation in RNP Granule Formation.

Ribonucleoprotein (RNP) granules are RNA-protein assemblies that are involved in multiple aspects of RNA metabolism and are linked to memory, development, and disease. Some RNP granules form, in part, through the formation of intermolecular RNA-RNA interactions. In vitro, such trans RNA condensation occurs readily, suggesting that cells require mechanisms to modulate RNA-based condensation. We assess the mechanisms of RNA condensation and how cells modulate this phenomenon. We propose that cells control RNA condensation through ATP-dependent processes, static RNA buffering, and dynamic post-translational mechanisms. Moreover, perturbations in these mechanisms can be involved in disease. This reveals multiple cellular mechanisms of kinetic and thermodynamic control that maintain the proper distribution of RNA molecules between dispersed and condensed forms.

Multitier regulation of the E. coli extreme acid stress response by CsrA.

CsrA is an RNA-binding protein that regulates processes critical for growth and survival, including central carbon metabolism, motility, biofilm formation, stress responses, and expression of virulence factors in pathogens. Transcriptomics studies in Escherichia coli suggested that CsrA repressed genes involved in surviving extremely acidic conditions. Here, we examine the effects of disrupting CsrA-dependent regulation on the expression of genes and circuitry for acid stress survival and demonstrate CsrA-mediated repression at multiple levels. We show that this repression is critical for managing the trade-off between growth and survival; overexpression of acid stress genes caused by csrA disruption enhances survival under extreme acidity but is detrimental for growth under mildly acidic conditions. In vitro studies confirmed that CsrA binds specifically to mRNAs of structural and regulatory genes for acid stress survival, causing translational repression. We also found that translation of the top-tier acid stress regulator, evgA, is coupled to that of a small leader peptide, evgL, which is repressed by CsrA. Unlike dedicated acid stress response genes, csrA and its sRNA antagonists, csrB and csrC, did not exhibit a substantial response to acid shock. Furthermore, disruption of CsrA regulation of acid stress genes impacted host-microbe interactions in Caenorhabditis elegans, alleviating GABA deficiencies. This study expands the known regulon of CsrA to genes of the extreme acid stress response of E. coli and highlights a new facet of the global role played by CsrA in balancing the opposing physiological demands of stress resistance with the capacity for growth and modulating host interactions.IMPORTANCETo colonize/infect the mammalian intestinal tract, bacteria must survive exposure to the extreme acidity of the stomach. E. coli does this by expressing proteins that neutralize cytoplasmic acidity and cope with molecular damage caused by low pH. Because of the metabolic cost of these processes, genes for surviving acid stress are tightly regulated. Here, we show that CsrA negatively regulates the cascade of expression responsible for the acid stress response. Increased expression of acid response genes due to csrA disruption improved survival at extremely low pH but inhibited growth under mildly acidic conditions. Our findings define a new layer of regulation in the acid stress response of E. coli and a novel physiological function for CsrA.

Understudied proteins and understudied functions in the model bacterium Bacillus subtilis-A major challenge in current research.

Model organisms such as the Gram-positive bacterium Bacillus subtilis have been studied intensively for decades. However, even for model organisms, no function has been identified for about one fourth of all proteins. It has recently been realized that such understudied proteins as well as poorly studied functions set a limitation to our understanding of the requirements for cellular life, and the Understudied Proteins Initiative has been launched. Of poorly studied proteins, those that are strongly expressed are likely to be important to the cell and should therefore be considered high priority in further studies. Since the functional analysis of unknown proteins can be extremely laborious, a minimal knowledge is required prior to targeted functional studies. In this review, we discuss strategies to obtain such a minimal annotation, for example, from global interaction, expression, or localization studies. We present a set of 41 highly expressed and poorly studied proteins of B. subtilis. Several of these proteins are thought or known to bind RNA and/or the ribosome, some may control the metabolism of B. subtilis, and another subset of particularly small proteins may act as regulatory elements to control the expression of downstream genes. Moreover, we discuss the challenges of poorly studied functions with a focus on RNA-binding proteins, amino acid transport, and the control of metabolic homeostasis. The identification of the functions of the selected proteins not only will strongly advance our knowledge on B. subtilis, but also on other organisms since many of the proteins are conserved in many groups of bacteria.

Determinants of the interaction between the 5'-leader of HIV-1 genome and human lysyl-tRNA synthetase in reverse transcription primer release process.

Reverse transcription of human immunodeficiency virus type 1 (HIV-1) initiates from the 3' end of human tRNALys3. The primer tRNALys3 is selectively packaged into the virus in the form of a complex with human lysyl-tRNA synthetase (LysRS). To facilitate reverse transcription initiation, part of the 5' leader (5'L) of HIV-1 genomic RNA (gRNA) evolves a tRNA anticodon-like element (TLE), which binds LysRS and releases tRNALys3 for primer annealing and reverse transcription initiation. Although TLE has been identified as a key element in 5'L responsible for LysRS binding, how the conformations and various hairpin structures of 5'L regulate 5'L-LysRS interaction is not fully understood. Here, these factors have been individually investigated using direct and competitive fluorescence anisotropy binding experiments. Our data showed that the conformation of 5'L significantly influences its binding affinity with LysRS. The 5'L conformation favoring gRNA dimerization and packaging exhibits much weaker binding affinity with LysRS compared to the alternative 5'L conformation that is not selected for packaging. Additionally, dimerization of 5'L impairs LysRS-5'L interaction. Furthermore, among various regions of 5'L, both the primer binding site/TLE domain and the stem-loop 3 are important for LysRS interaction, whereas the dimerization initiation site and the splicing donor plays a minor role. In contrast, the presence of the transacting responsive and the polyadenylation signal hairpins slightly inhibit LysRS binding. These findings reveal that the conformation and various regions of the 5'L of HIV-1 genome regulate its interaction with human LysRS and the reverse transcription primer release process.