CLIPick: a sensitive peak caller for expression-based deconvolution of HITS-CLIP signals

Widespread 8-oxoguanine modifications of miRNA seeds differentially regulate redox-dependent cancer development

Oxidative stress contributes to tumourigenesis by altering gene expression. One accompanying modification, 8-oxoguanine (o8G) can change RNA-RNA interactions via o8G•A base pairing, but its regulatory roles remain elusive. Here, on the basis of o8G-induced guanine-to-thymine (o8G > T) variations featured in sequencing, we discovered widespread position-specific o8Gs in tumour microRNAs, preferentially oxidized towards 5' end seed regions (positions 2-8) with clustered sequence patterns and clinically associated with patients in lower-grade gliomas and liver hepatocellular carcinoma. We validated that o8G at position 4 of miR-124 (4o8G-miR-124) and 4o8G-let-7 suppress lower-grade gliomas, whereas 3o8G-miR-122 and 4o8G-let-7 promote malignancy of liver hepatocellular carcinoma by redirecting the target transcriptome to oncogenic regulatory pathways. Stepwise oxidation from tumour-promoting 3o8G-miR-122 to tumour-suppressing 2,3o8G-miR-122 occurs and its specific modulation in mouse liver effectively attenuates diethylnitrosamine-induced hepatocarcinogenesis. These findings provide resources and insights into epitranscriptional o8G regulation of microRNA functions, reprogrammed by redox changes, implicating its control for cancer treatment.

Genome-wide cross-cancer analysis illustrates the critical role of bimodal miRNA in patient survival and drug responses to PI3K inhibitors

Heterogeneity of cancer means many tumorigenic genes are only aberrantly expressed in a subset of patients and thus follow a bimodal distribution, having two modes of expression within a single population. Traditional statistical techniques that compare sample means between cancer patients and healthy controls fail to detect bimodally expressed genes. We utilize a mixture modeling approach to identify bimodal microRNA (miRNA) across cancers, find consistent sources of heterogeneity, and identify potential oncogenic miRNA that may be used to guide personalized therapies. Pathway analysis was conducted using target genes of the bimodal miRNA to identify potential functional implications in cancer. In vivo overexpression experiments were conducted to elucidate the clinical importance of bimodal miRNA in chemotherapy treatments. In nine types of cancer, tumors consistently displayed greater bimodality than normal tissue. Specifically, in liver and lung cancers, high expression of miR-105 and miR-767 was indicative of poor prognosis. Functional pathway analysis identified target genes of miR-105 and miR-767 enriched in the phosphoinositide-3-kinase (PI3K) pathway, and analysis of over 200 cancer drugs in vitro showed that drugs targeting the same pathway had greater efficacy in cell lines with high miR-105 and miR-767 levels. Overexpression of the two miRNA facilitated response to PI3K inhibitor treatment. We demonstrate that while cancer is marked by considerable genetic heterogeneity, there is between-cancer concordance regarding the particular miRNA that are more variable. Bimodal miRNA are ideal biomarkers that can be used to stratify patients for prognosis and drug response in certain types of cancer.

Bimodal genes can be defined as those having two modes of expression within the same population. A variety of statistical methodologies have been employed to assess bimodal gene expression, but current methods and their applications have been limited. Given the advances in next-generation sequencing as well as the extensive regulatory role of miRNA, assessing bimodality in miRNA-seq data can greatly broaden our understanding of factors underlying tumor progression. The goal of the current study was to utilize a novel mixture modeling approach to identify bimodal miRNA and then demonstrate their importance in cancer by evaluating their ability to predict overall survival and drug response. Our results showed that high levels of bimodal miRNA expression was characteristic of cancer. Additionally, several bimodal miRNA were common to multiple cancer types, suggesting that certain miRNA consistently account for tumor heterogeneity and may be involved in general oncogenic processes. Our study points to the potential of bimodal miRNA to facilitate precise prognostic evaluation and effective treatment strategies.

AGO CLIP-based imputation of potent siRNA sequences targeting SARS-CoV-2 with antifibrotic miRNA-like activity

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is associated with fatal pulmonary fibrosis. Small interfering RNAs (siRNAs) can be developed to induce RNA interference against SARS-CoV-2, and their susceptible target sites can be inferred by Argonaute crosslinking immunoprecipitation sequencing (AGO CLIP). Here, by reanalysing AGO CLIP data in RNA viruses, we delineated putative AGO binding in the conserved non-structural protein 12 (nsp12) region encoding RNA-dependent RNA polymerase (RdRP) in SARS-CoV-2. We utilised the inferred AGO binding to optimise the local RNA folding parameter to calculate target accessibility and predict all potent siRNA target sites in the SARS-CoV-2 genome, avoiding sequence variants. siRNAs loaded onto AGO also repressed seed (positions 2–8)-matched transcripts by acting as microRNAs (miRNAs). To utilise this, we further screened 13 potential siRNAs whose seed sequences were matched to known antifibrotic miRNAs and confirmed their miRNA-like activity. A miR-27-mimicking siRNA designed to target the nsp12 region (27/RdRP) was validated to silence a synthesised nsp12 RNA mimic in lung cell lines and function as an antifibrotic miR-27 in regulating target transcriptomes related to TGF-β signalling. siRNA sequences with an antifibrotic miRNA-like activity that could synergistically treat COVID-19 are available online (http://clip.korea.ac.kr/covid19).

AGO-accessible anticancer siRNAs designed with synergistic miRNA-like activity

Small interfering RNAs (siRNAs) therapeutically induce RNA interference (RNAi) of disease-causing genes, but they also silence hundreds of seed-matched off-targets as behaving similar to microRNAs (miRNAs). miRNAs control the pathophysiology of tumors, wherein their accessible binding sites can be sequenced by Argonaute crosslinking immunoprecipitation (AGO CLIP). Herein, based on AGO CLIP, we develop potent anticancer siRNAs utilizing miRNA-like activity (mi/siRNAs). The mi/siRNAs contain seed sequences (positions 2-7) of tumor-suppressive miRNAs while maintaining perfect sequence complementarity to the AGO-accessible tumor target sites. Initially, host miRNA interactions with human papillomavirus 18 (HPV18) were identified in cervical cancer by AGO CLIP, revealing tumor-suppressive activity of miR-1/206 and miR-218. Based on the AGO-miRNA binding sites, mi/siRNAs were designed to target E6 and E7 (E6/E7) transcript with seed sequences of miR-1/206 (206/E7) and miR-218 (218/E7). Synergistic anticancer activity of 206/E7 and 218/E7 was functionally validated and confirmed via RNA sequencing and in vivo xenograft models (206/E7). Other mi/siRNA sequences were additionally designed for cervical, ovarian, and breast cancer, and available as an online tool (http://ago.korea.ac.kr/misiRNA); some of the mi/siRNAs were validated for their augmented anticancer activity (206/EphA2 and 206/Her2). mi/siRNAs could coordinate miRNA-like activity with robust siRNA function, demonstrating the potential of AGO CLIP analysis for RNAi therapeutics.

Galaxy CLIP-Explorer: a web server for CLIP-Seq data analysis

BACKGROUND

Post-transcriptional regulation via RNA-binding proteins plays a fundamental role in every organism, but the regulatory mechanisms lack important understanding. Nevertheless, they can be elucidated by cross-linking immunoprecipitation in combination with high-throughput sequencing (CLIP-Seq). CLIP-Seq answers questions about the functional role of an RNA-binding protein and its targets by determining binding sites on a nucleotide level and associated sequence and structural binding patterns. In recent years the amount of CLIP-Seq data skyrocketed, urging the need for an automatic data analysis that can deal with different experimental set-ups. However, noncanonical data, new protocols, and a huge variety of tools, especially for peak calling, made it difficult to define a standard.

FINDINGS

CLIP-Explorer is a flexible and reproducible data analysis pipeline for iCLIP data that supports for the first time eCLIP, FLASH, and uvCLAP data. Individual steps like peak calling can be changed to adapt to different experimental settings. We validate CLIP-Explorer on eCLIP data, finding similar or nearly identical motifs for various proteins in comparison with other databases. In addition, we detect new sequence motifs for PTBP1 and U2AF2. Finally, we optimize the peak calling with 3 different peak callers on RBFOX2 data, discuss the difficulty of the peak-calling step, and give advice for different experimental set-ups.

CONCLUSION

CLIP-Explorer finally fills the demand for a flexible CLIP-Seq data analysis pipeline that is applicable to the up-to-date CLIP protocols. The article further shows the limitations of current peak-calling algorithms and the importance of a robust peak detection.

Position-specific oxidation of miR-1 encodes cardiac hypertrophy

In pathophysiology, reactive oxygen species oxidize biomolecules that contribute to disease phenotypes¹. One such modification, 8-oxoguanine² (o⁸G), is abundant in RNA³ but its epitranscriptional role has not been investigated for microRNAs (miRNAs). Here we specifically sequence oxidized miRNAs in a rat model of the redox-associated condition cardiac hypertrophy⁴. We find that position-specific o⁸G modifications are generated in seed regions (positions 2-8) of selective miRNAs, and function to regulate other mRNAs through o⁸G•A base pairing. o⁸G is induced predominantly at position 7 of miR-1 (7o⁸G-miR-1) by treatment with an adrenergic agonist. Introducing 7o⁸G-miR-1 or 7U-miR-1 (in which G at position 7 is substituted with U) alone is sufficient to cause cardiac hypertrophy in mice, and the mRNA targets of o⁸G-miR-1 function in affected phenotypes; the specific inhibition of 7o⁸G-miR-1 in mouse cardiomyocytes was found to attenuate cardiac hypertrophy. o⁸G-miR-1 is also implicated in patients with cardiomyopathy. Our findings show that the position-specific oxidation of miRNAs could serve as an epitranscriptional mechanism to coordinate pathophysiological redox-mediated gene expression.

LNRLMI: Linear neighbour representation for predicting lncRNA-miRNA interactions

LncRNA and miRNA are key molecules in mechanism of competing endogenous RNAs(ceRNA), and their interactions have been discovered with important roles in gene regulation. As supplementary to the identification of lncRNA‐miRNA interactions from CLIP‐seq experiments, in silico prediction can select the most potential candidates for experimental validation. Although developing computational tool for predicting lncRNA‐miRNA interaction is of great importance for deciphering the ceRNA mechanism, little effort has been made towards this direction. In this paper, we propose an approach based on linear neighbour representation to predict lncRNA‐miRNA interactions (LNRLMI). Specifically, we first constructed a bipartite network by combining the known interaction network and similarities based on expression profiles of lncRNAs and miRNAs. Based on such a data integration, linear neighbour representation method was introduced to construct a prediction model. To evaluate the prediction performance of the proposed model, k‐fold cross validations were implemented. As a result, LNRLMI yielded the average AUCs of 0.8475 ± 0.0032, 0.8960 ± 0.0015 and 0.9069 ± 0.0014 on 2‐fold, 5‐fold and 10‐fold cross validation, respectively. A series of comparison experiments with other methods were also conducted, and the results showed that our method was feasible and effective to predict lncRNA‐miRNA interactions via a combination of different types of useful side information. It is anticipated that LNRLMI could be a useful tool for predicting non‐coding RNA regulation network that lncRNA and miRNA are involved in.