ASSA: Fast identification of statistically significant interactions between long RNAs

Elucidating the regulatory role of long non-coding RNAs in drought stress response during seed germination in leaf mustard

Leaf mustard (Brassica juncea L. Czern & Coss), an important vegetable crop, experiences pronounced adversity due to seasonal drought stress, particularly at the seed germination stage. Although there is partial comprehension of drought-responsive genes, the role of long non-coding RNAs (lncRNAs) in adjusting mustard's drought stress response is largely unexplored. In this study, we showed that the drought-tolerant cultivar 'Weiliang' manifested a markedly lower base water potential (-1.073 MPa vs -0.437 MPa) and higher germination percentage (41.2% vs 0%) than the drought-susceptible cultivar 'Shuidong' under drought conditions. High throughput RNA sequencing techniques revealed a significant repertoire of lncRNAs from both cultivars during germination under drought stress, resulting in the identification of 2,087 differentially expressed lncRNAs (DELs) and their correspondingly linked 12,433 target genes. It was noted that 84 genes targeted by DEL exhibited enrichment in the photosynthesis pathway. Gene network construction showed that MSTRG.150397, a regulatory lncRNA, was inferred to potentially modulate key photosynthetic genes (Psb27, PetC, PetH, and PsbW), whilst MSTRG.107159 was indicated as an inhibitory regulator of six drought-responsive PIP genes. Further, weighted gene co-expression network analysis (WGCNA) corroborated the involvement of light intensity and stress response genes targeted by the identified DELs. The precision and regulatory impact of lncRNA were verified through qPCR. This study extends our knowledge of the regulatory mechanisms governing drought stress responses in mustard, which will help strategies to augment drought tolerance in this crop.

How to do RNA-RNA Interaction Prediction? A Use-Case Driven Handbook Using IntaRNA

Computational prediction of RNA-RNA interactions (RRI) is a central methodology for the specific investigation of inter-molecular RNA interactions and regulatory effects of non-coding RNAs like eukaryotic microRNAs or prokaryotic small RNAs. Available methods can be classified according to their underlying prediction strategies, each implicating specific capabilities and restrictions often not transparent to the non-expert user. Within this work, we review seven classes of RRI prediction strategies and discuss the advantages and limitations of respective tools, since such knowledge is essential for selecting the right tool in the first place.Among the RRI prediction strategies, accessibility-based approaches have been shown to provide the most reliable predictions. Here, we describe how IntaRNA, as one of the state-of-the-art accessibility-based tools, can be applied in various use cases for the task of computational RRI prediction. Detailed hands-on examples for individual RRI predictions as well as large-scale target prediction scenarios are provided. We illustrate the flexibility and capabilities of IntaRNA through the examples. Each example is designed using real-life data from the literature and is accompanied by instructions on interpreting the respective results from IntaRNA output. Our use-case driven instructions enable non-expert users to comprehensively understand and utilize IntaRNA's features for effective RRI predictions.

LncRTPred: Predicting RNA-RNA mode of interaction mediated by lncRNA

Long non-coding RNAs (lncRNAs) play a significant role in various biological processes. Hence, it is utmost important to elucidate their functions in order to understand the molecular mechanism of a complex biological system. This versatile RNA molecule has diverse modes of interaction, one of which constitutes lncRNA-mRNA interaction. Hence, identifying its target mRNA is essential to understand the function of an lncRNA explicitly. Existing lncRNA target prediction tools mainly adopt thermodynamics approach. Large execution time and inability to perform real-time prediction limit their usage. Further, lack of negative training dataset has been a hindrance in the path of developing machine learning (ML) based lncRNA target prediction tools. In this work, we have developed a ML-based lncRNA-mRNA target prediction model- 'LncRTPred'. Here we have addressed the existing problems by generating reliable negative dataset and creating robust ML models. We have identified the non-interacting lncRNA and mRNAs from the unlabelled dataset using BLAT. It is further filtered to get a reliable set of outliers. LncRTPred provides a cumulative_model_score as the final output against each query. In terms of prediction accuracy, LncRTPred outperforms other popular target prediction protocols like LncTar. Further, we have tested its performance against experimentally validated disease-specific lncRNA-mRNA interactions. Overall, performance of LncRTPred is heavily dependent on the size of the training dataset, which is highly reflected by the difference in its performance for human and mouse species. Its performance for human species shows better as compared to that for mouse when applied on an unknown data due to smaller size of the training dataset in case of mouse compared to that of human. Availability of increased number of lncRNA-mRNA interaction data for mouse will improve the performance of LncRTPred in future. Both webserver and standalone versions of LncRTPred are available. Web server link: http://bicresources.jcbose.ac.in/zhumur/lncrtpred/index.html. Github Link: https://github.com/zglabDIB/LncRTPred.

An in vivo model of glioblastoma radiation resistance identifies long non-coding RNAs and targetable kinases

Key molecular regulators of acquired radiation resistance in recurrent glioblastoma (GBM) are largely unknown, with a dearth of accurate preclinical models. To address this, we generated 8 GBM patient-derived xenograft (PDX) models of acquired radiation therapy–selected (RTS) resistance compared with same-patient, treatment-naive (radiation-sensitive, unselected; RTU) PDXs. These likely unique models mimic the longitudinal evolution of patient recurrent tumors following serial radiation therapy. Indeed, while whole-exome sequencing showed retention of major genomic alterations in the RTS lines, we did detect a chromosome 12q14 amplification that was associated with clinical GBM recurrence in 2 RTS models. A potentially novel bioinformatics pipeline was applied to analyze phenotypic, transcriptomic, and kinomic alterations, which identified long noncoding RNAs (lncRNAs) and targetable, PDX-specific kinases. We observed differential transcriptional enrichment of DNA damage repair pathways in our RTS models, which correlated with several lncRNAs. Global kinomic profiling separated RTU and RTS models, but pairwise analyses indicated that there are multiple molecular routes to acquired radiation resistance. RTS model–specific kinases were identified and targeted with clinically relevant small molecule inhibitors. This cohort of in vivo RTS patient-derived models will enable future preclinical therapeutic testing to help overcome the treatment resistance seen in patients with GBM.

Approaches to Identify and Characterise the Post-Transcriptional Roles of lncRNAs in Cancer

It is becoming increasingly evident that the non-coding genome and transcriptome exert great influence over their coding counterparts through complex molecular interactions. Among non-coding RNAs (ncRNA), long non-coding RNAs (lncRNAs) in particular present increased potential to participate in dysregulation of post-transcriptional processes through both RNA and protein interactions. Since such processes can play key roles in contributing to cancer progression, it is desirable to continue expanding the search for lncRNAs impacting cancer through post-transcriptional mechanisms. The sheer diversity of mechanisms requires diverse resources and methods that have been developed and refined over the past decade. We provide an overview of computational resources as well as proven low-to-high throughput techniques to enable identification and characterisation of lncRNAs in their complex interactive contexts. As more cancer research strategies evolve to explore the non-coding genome and transcriptome, we anticipate this will provide a valuable primer and perspective of how these technologies have matured and will continue to evolve to assist researchers in elucidating post-transcriptional roles of lncRNAs in cancer.

Direct Interactions with Nascent Transcripts Is Potentially a Common Targeting Mechanism of Long Non-Coding RNAs

Although thousands of mammalian long non-coding RNAs (lncRNAs) have been reported in the last decade, their functional annotation remains limited. A wet-lab approach to detect functions of a novel lncRNA usually includes its knockdown followed by RNA sequencing and identification of the deferentially expressed genes. However, identification of the molecular mechanism(s) used by the lncRNA to regulate its targets frequently becomes a challenge. Previously, we developed the ASSA algorithm that detects statistically significant inter-molecular RNA-RNA interactions. Here we designed a workflow that uses ASSA predictions to estimate the ability of an lncRNA to function via direct base pairing with the target transcripts (co- or post-transcriptionally). The workflow was applied to 300+ lncRNA knockdown experiments from the FANTOM6 pilot project producing statistically significant predictions for 71 unique lncRNAs (104 knockdowns). Surprisingly, the majority of these lncRNAs were likely to function co-transcriptionally, i.e., hybridize with the nascent transcripts of the target genes. Moreover, a number of the obtained predictions were supported by independent iMARGI experimental data on co-localization of lncRNA and chromatin. We detected an evolutionarily conserved lncRNA CHASERR (AC013394.2 or LINC01578) that could regulate target genes co-transcriptionally via interaction with a nascent transcript by directing CHD2 helicase. The obtained results suggested that this nuclear lncRNA may be able to activate expression of the target genes in trans by base-pairing with the nascent transcripts and directing the CHD2 helicase to the regulated promoters leading to open the chromatin and active transcription. Our study highlights the possible importance of base-pairing between nuclear lncRNAs and nascent transcripts for the regulation of gene expression.

Prediction of lncRNAs and their interactions with nucleic acids: benchmarking bioinformatics tools

The genomes of mammalian species are pervasively transcribed producing as many noncoding as protein-coding RNAs. There is a growing body of evidence supporting their functional role. Long noncoding RNA (lncRNA) can bind both nucleic acids and proteins through several mechanisms. A reliable computational prediction of the most probable mechanism of lncRNA interaction can facilitate experimental validation of its function. In this study, we benchmarked computational tools capable to discriminate lncRNA from mRNA and predict lncRNA interactions with other nucleic acids. We assessed the performance of 9 tools for distinguishing protein-coding from noncoding RNAs, as well as 19 tools for prediction of RNA-RNA and RNA-DNA interactions. Our conclusions about the considered tools were based on their performances on the entire genome/transcriptome level, as it is the most common task nowadays. We found that FEELnc and CPAT distinguish between coding and noncoding mammalian transcripts in the most accurate manner. ASSA, RIBlast and LASTAL, as well as Triplexator, turned out to be the best predictors of RNA-RNA and RNA-DNA interactions, respectively. We showed that the normalization of the predicted interaction strength to the transcript length and GC content may improve the accuracy of inferring RNA interactions. Yet, all the current tools have difficulties to make accurate predictions of short-trans RNA-RNA interactions-stretches of sparse contacts. All over, there is still room for improvement in each category, especially for predictions of RNA interactions.

Gibbs/MCMC Sampling for Multiple RNA Interaction with Sub-optimal Solutions

Multiple RNA interaction can be modeled as a problem in combinatorial optimization, where the "optimal" structure is driven by an energy-minimization-like algorithm. However, the actual structure may not be optimal in this computational sense. Moreover, it is not necessarily unique. Therefore, alternative sub-optimal solutions are needed to cover the biological ground. We present a combinatorial formulation for the Multiple RNA Interaction problem with approximation algorithms to handle various interaction patterns, which when combined with Gibbs sampling and MCMC (Markov Chain Monte Carlo), can efficiently generate a reasonable number of optimal and sub-optimal solutions. When viable structures are far from an optimal solution, exploring dependence among different parts of the interaction can increase their score and boost their candidacy for the sampling algorithm. By clustering the solutions, we identify few representatives that are distinct enough to suggest possible alternative structures.