isomiRTar: a comprehensive portal of pan-cancer 5'-isomiR targeting

Incautious design of shRNAs for stable overexpression of miRNAs could result in generation of undesired isomiRs

shRNA-mediated strategy of miRNA overexpression based on RNA Polymerase (Pol III) expression cassettes is widely used for miRNA functional studies. For some miRNAs, e.g., encoded in the genome as a part of a polycistronic miRNA cluster, it is most likely the only way for their individual stable overexpression. Here we have revealed that expression of miRNAs longer than 19 nt (e.g. 23 nt in length hsa-miR-93-5p) using such approach could be accompanied by undesired predominant generation of 5' end miRNA isoforms (5'-isomiRs). Extra U residues (up to five) added by Pol III at the 3' end of the transcribed shRNA during transcription termination could cause a shift in the Dicer cleavage position of the shRNA. This results in the formation of 5'-isomiRs, which have a significantly altered seed region compared to the initially encoded canonical hsa-miR-93-5p. We demonstrated that the commonly used qPCR method is insensitive to the formation of 5'-isomiRs and cannot be used to confirm miRNA overexpression. However, the predominant expression of 5'-isomiRs without three or four first nucleotides instead of the canonical isoform could be disclosed based on miRNA-Seq analysis. Moreover, mRNA sequencing data showed that the 5'-isomiRs of hsa-miR-93-5p presumably regulate their own mRNA targets. Thus, omitting miRNA-Seq analysis may lead to erroneous conclusions regarding revealed mRNA targets and possible molecular mechanisms in which studied miRNA is involved. Overall, the presented results show that structures of shRNAs for stable overexpression of miRNAs requires careful design to avoid generation of undesired 5'-isomiRs.

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

Integrated microRNA and proteome analysis of cancer datasets with MoPC

MicroRNAs (miRNAs) are small molecules that play an essential role in regulating gene expression by post-transcriptional gene silencing. Their study is crucial in revealing the fundamental processes underlying pathologies and, in particular, cancer. To date, most studies on miRNA regulation consider the effect of specific miRNAs on specific target mRNAs, providing wet-lab validation. However, few tools have been developed to explain the miRNA-mediated regulation at the protein level. In this paper, the MoPC computational tool is presented, that relies on the partial correlation between mRNAs and proteins conditioned on the miRNA expression to predict miRNA-target interactions in multi-omic datasets. MoPC returns the list of significant miRNA-target interactions and plot the significant correlations on the heatmap in which the miRNAs and targets are ordered by the chromosomal location. The software was applied on three TCGA/CPTAC datasets (breast, glioblastoma, and lung cancer), returning enriched results in three independent targets databases.

MicroRNAs Present in Malignant Pleural Fluid Increase the Migration of Normal Mesothelial Cells In Vitro and May Help Discriminate between Benign and Malignant Effusions

The sensitivity of pleural fluid (PF) analyses for the diagnosis of malignant pleural effusions (MPEs) is low to moderate. Knowledge about the pathobiology and molecular characteristics of this condition is limited. In this study, the crosstalk between stromal cells and tumor cells was investigated in vitro in order to reveal factors that are present in PF which can mediate MPE formation and aid in discriminating between benign and malignant etiologies. Eighteen PF samples, in different proportions, were exposed in vitro to mesothelial MeT-5A cells to determine the biological effects on these cells. Treatment of normal mesothelial MeT-5A cells with malignant PF increased cell viability, proliferation, and migration, and activated different survival-related signaling pathways. We identified differentially expressed miRNAs in PF samples that could be responsible for these changes. Consistently, bioinformatics analysis revealed an enrichment of the discovered miRNAs in migration-related processes. Notably, the abundance of three miRNAs (miR-141-3p, miR-203a-3, and miR-200c-3p) correctly classified MPEs with false-negative cytological examination results, indicating the potential of these molecules for improving diagnosis. Malignant PF produces phenotypic and functional changes in normal mesothelial cells. These changes are partly mediated by certain miRNAs, which, in turn, could serve to differentiate malignant from benign effusions.

The signature of SARS-CoV-2 evolution reflects selective pressures within human guts

In somatic cells, microRNAs (miRNAs) bind to the genomes of RNA viruses and influence their translation and replication. In London and Berlin samples represented in GISAID database, we traced severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages and divided these sequenced in two groups, "Ancestral variants" and "Omicrons," and analyzed them through the prism of the tissue-specific binding between host miRNAs and viral messenger RNAs. We demonstrate a significant number of miRNA-binding sites in the NSP4 region of the SARS-CoV-2 genome, with evidence of evolutionary pressure within this region exerted by human intestinal miRNAs. Notably, in infected cells, NSP4 promotes the formation of double-membrane vesicles, which serve as the scaffolds for replication-transcriptional complexes and protect viral RNA from intracellular destruction. In 3 years of selection, the loss of many miRNA-binding sites in general and those within the NSP4 in particular has shaped the SARS-CoV-2 genomes. With that, the descendants of the BA.2 variants were promoted as dominant strains, which define current momentum of the pandemics.

A novel approach for a joint analysis of isomiR and mRNA expression data reveals features of isomiR targeting in breast cancer

A widely used procedure for selecting significant miRNA-mRNA or isomiR-mRNA pairs out of predicted interactions involves calculating the correlation between expression levels of miRNAs/isomiRs and mRNAs in a series of samples. In this manuscript, we aimed to assess the validity of this procedure by comparing isomiR-mRNA correlation profiles in sets of sequence-based predicted target mRNAs and non-target mRNAs (negative controls). Target prediction was carried out using RNA22 and TargetScan algorithms. Spearman's correlation analysis was conducted using miRNA and mRNA sequencing data of The Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA) project. Luminal A, luminal B, basal-like breast cancer subtypes, and adjacent normal tissue samples were analyzed separately. Using the sets of putative targets and non-targets, we introduced adjusted isomiR targeting activity (ITA)-the number of negatively correlated potential isomiR targets adjusted by the background (estimated using non-target mRNAs). We found that for most isomiRs a significant negative correlation between isomiR-mRNA expression levels appeared more often in a set of predicted targets compared to the non-targets. This trend was detected for both classical seed region binding types (8mer, 7mer-m8, 7mer-A1, 6mer) predicted by TargetScan and the non-classical ones (G:U wobbles and up to one mismatch or unpaired nucleotide within seed sequence) predicted by RNA22. Adjusted ITA distributions were similar for target sites located in 3'-UTRs and coding mRNA sequences, while 5'-UTRs had much lower scores. Finally, we observed strong cancer subtype-specific patterns of isomiR activity, highlighting the differences between breast cancer molecular subtypes and normal tissues. Surprisingly, our target prediction- and correlation-based estimates of isomiR activities were practically non-correlated with the average isomiR expression levels neither in cancerous nor in normal samples.