Splice-switching of the insulin receptor pre-mRNA alleviates tumorigenic hallmarks in rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is an aggressive pediatric tumor with a poor prognosis for metastasis and recurrent disease. Large-scale sequencing endeavors demonstrate that Rhabdomyosarcomas have a dearth of precisely targetable driver mutations. However, IGF-2 signaling is known to be grossly altered in RMS. The insulin receptor (IR) exists in two alternatively spliced isoforms, IR-A and IR-B. The IGF-2 signaling molecule binds both its innate IGF-1 receptor as well as the insulin receptor variant A (IR-A) with high affinity. Mitogenic and proliferative signaling via the canonical IGF-2 pathway is, therefore, augmented by IR-A. This study shows that RMS patients express increased IR-A levels compared to control tissues that predominantly express the IR-B isoform. We also found that Hif-1α is significantly increased in RMS tumors, portraying their hypoxic phenotype. Concordantly, the alternative splicing of IR adapts to produce more IR-A in response to hypoxic stress. Upon examining the pre-mRNA structure of the gene, we identified a potential hypoxia-responsive element, which is also the binding site for the RNA-binding protein CUG-BP1 (CELF1). We designed Splice Switching Oligonucleotides (SSO) against this binding site to decrease IR-A levels in RMS cell lines and, consequently, rescue the IR-B expression levels. SSO treatment resulted in a significant reduction in cell proliferation, migration, and angiogenesis. Our data shows promising insight into how impeding the IGF-2 pathway by reducing IR-A expression mitigates tumor growth. It is evident that Rhabdomyosarcomas use IR alternative splicing as yet another survival strategy that can be exploited as a therapeutic intervention in conjunction with already established anti-IGF-1 receptor therapies.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry

Background Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry. Results Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat. Conclusions We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Modulation of pre-mRNA structure by hnRNP proteins regulates alternative splicing of MALT1

Alternative splicing is controlled by differential binding of trans-acting RNA binding proteins (RBPs) to cis-regulatory elements in intronic and exonic pre-mRNA regions. How secondary structure in the pre-mRNA transcripts affects recognition by RBPs and determines alternative exon usage is poorly understood. The MALT1 paracaspase is a key component of signaling pathways that mediate innate and adaptive immune responses. Alternative splicing of MALT1 exon7 is critical for controlling optimal T cell activation. Here, we demonstrate that processing of the MALT1 pre-mRNA depends on RNA structural elements that shield the 5′ and 3′ splice sites of the alternatively spliced exon7. By combining biochemical analyses with chemical probing and NMR we show that the RBPs hnRNP U and hnRNP L bind competitively and with comparable affinities to identical stem-loop RNA structures flanking the 5′ and 3′ splice sites of MALT1 exon7. While hnRNP U stabilizes RNA stem-loop conformations that maintain exon7 skipping, hnRNP L unwinds these RNA elements to facilitate recruitment of the essential splicing factor U2AF2 to promote exon7 inclusion. Our data represent a paradigm for the control of splice site selection by differential RBP binding and modulation of pre-mRNA structure.

Quantitative prediction of variant effects on alternative splicing using endogenous pre-messenger RNA structure probing

Splicing is a highly regulated process that depends on numerous factors. It is particularly challenging to quantitatively predict how a mutation will affect precursor messenger RNA (mRNA) structure and the subsequent functional consequences. Here we use a novel Mutational Profiling (-MaP) methodology to obtain highly reproducible endogenous precursor and mature mRNA structural probing data in vivo. We use these data to estimate Boltzmann suboptimal ensembles, and predict the structural consequences of mutations on precursor mRNA structure. Together with a structural analysis of recent cryo-EM spliceosome structures at different stages of the splicing cycle, we determined that the footprint of the Bact complex on precursor mRNA is best able to predict splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene. However, structure alone only achieves 74% accuracy. We therefore developed a β-regression weighting framework that incorporates splice site strength, structure and exonic/intronic splicing regulatory elements which together achieves 90% accuracy for 47 known and six newly discovered splice-altering variants. This combined experimental/computational framework represents a path forward for accurate prediction of splicing related disease-causing variants.

Ribosomal profiling as a tool for studying translation in plants: main results, problems and future prospects

The expression of eukaryotic genes can be regulated at several stages, including the translation of mRNA. It is known that the structure of mRNA can affect both the efficiency of interaction with the translation apparatus in general and the choice of translation initiation sites. To study the translated fraction of the transcriptome, experimental methods of analysis were developed, the most informative of which is ribosomal profiling (RP, Ribo-seq). Originally developed for use in yeast systems, this method has been adapted for research in translation mechanisms in many plant species. This technology includes the isolation of the polysomal fraction and high-performance sequencing of a pool of mRNA fragments associated with ribosomes. Comparing the results of transcript coverage with reads obtained using the ribosome profiling with the transcriptional efficiency of genes allows the translation efficiency to be evaluated for each transcript. The exact positions of ribosomes determined on mRNA sequences allow determining the translation of open reading frames and switching between the translation of several reading frames – a phenomenon in which two or more overlapping frames are read from one mRNA and different proteins are synthesized. The advantage of this method is that it provides quantitative estimates of ribosome coverage of mRNA and can detect relatively rare translation events. Using this technology, it was possible to identify and classify plant genes by the type of regulation of their expression at the transcription, translation, or both levels. Features of the mRNA structure that affect translation levels have been revealed: the formation of G2 quadruplexes and the presence of specific motifs in the 5’-UTR region, GC content, the presence of alternative translation starts, and the influence of uORFs on the translation of downstream mORFs. In this review, we briefly reviewed the RP methodology and the prospects for its application to study the structural and functional organization and regulation of plant gene expression.

From A to m6A: The Emerging Viral Epitranscriptome

There are over 100 different chemical RNA modifications, collectively known as the epitranscriptome. N6-methyladenosine (m6A) is the most commonly found internal RNA modification in cellular mRNAs where it plays important roles in the regulation of the mRNA structure, stability, translation and nuclear export. This modification is also found in viral RNA genomes and in viral mRNAs derived from both RNA and DNA viruses. A growing body of evidence indicates that m6A modifications play important roles in regulating viral replication by interacting with the cellular m6A machinery. In this review, we will exhaustively detail the current knowledge on m6A modification, with an emphasis on its function in virus biology.

The RNA chaperone protein CspA stimulates translation during cold acclimation by promoting the progression of the ribosomes

CspA is an RNA binding protein expressed during cold-shock in Escherichia coli, capable of stimulating translation of several mRNAs - including its own - at low temperature. We used reconstituted translation systems to monitor the effects of CspA on the different steps of the translation process and probing experiments to analyze the interactions with its target mRNAs. We specifically focused on cspA mRNA which adopts a cold-induced secondary structure at temperatures below 20°C and a more closed conformation at 37°C. We show that at low temperature CspA specifically promotes the translation of the mRNA folded in the conformation less accessible to the ribosome (37°C form). CspA interacts with its mRNA without inducing large structural rearrangement, does not bind the ribosomal subunits and is not able to stimulate the formation of the translation initiation complexes. On the other hand, CspA promotes the progression of the ribosomes during translation of its mRNA at low temperature and this stimulation is mRNA structure-dependent. A similar structure-dependent mechanism may be responsible for the CspA-dependent translation stimulation observed with other probed mRNAs, for which the transition to the elongation phase is progressively facilitated during cold acclimation with the accumulation of CspA.

A translational riboswitch coordinates nascent transcription–translation coupling

Bacterial messenger RNA (mRNA) synthesis by RNA polymerase (RNAP) and first-round translation by the ribosome are often coupled to regulate gene expression, yet how coupling is established and maintained is ill understood. Here, we develop biochemical and single-molecule fluorescence approaches to probe the dynamics of RNAP–ribosome interactions on an mRNA with a translational preQ1-sensing riboswitch in its 5′ untranslated region. Binding of preQ1 leads to the occlusion of the ribosome binding site (RBS), inhibiting translation initiation. We demonstrate that RNAP poised within the mRNA leader region promotes ribosomal 30S subunit binding, antagonizing preQ1-induced RBS occlusion, and that the RNAP–30S bridging transcription factors NusG and RfaH distinctly enhance 30S recruitment and retention, respectively. We further find that, while 30S–mRNA interaction significantly impedes RNAP in the absence of translation, an actively translating ribosome promotes productive transcription. A model emerges wherein mRNA structure and transcription factors coordinate to dynamically modulate the efficiency of transcription–translation coupling.

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population

Background The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of “silent” genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking. Results We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population. Conclusions These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a “Structural Predictivity Index” (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

Combinatorial optimization of mRNA structure, stability, and translation for RNA-based therapeutics

SUMMARYTherapeutic mRNAs and vaccines are being developed for a broad range of human diseases, including COVID-19. However, their optimization is hindered by mRNA instability and inefficient protein expression. Here, we describe design principles that overcome these barriers. We develop a new RNA sequencing-based platform called PERSIST-seq to systematically delineate in-cell mRNA stability, ribosome load, as well as in-solution stability of a library of diverse mRNAs. We find that, surprisingly, in-cell stability is a greater driver of protein output than high ribosome load. We further introduce a method called In-line-seq, applied to thousands of diverse RNAs, that reveals sequence and structure-based rules for mitigating hydrolytic degradation. Our findings show that “superfolder” mRNAs can be designed to improve both stability and expression that are further enhanced through pseudouridine nucleoside modification. Together, our study demonstrates simultaneous improvement of mRNA stability and protein expression and provides a computational-experimental platform for the enhancement of mRNA medicines.