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.

The miRNA-target interactions: An underestimated intricacy.

MicroRNAs (miRNAs) play indispensable roles in posttranscriptional gene regulation. Their cellular regulatory impact is determined not solely by their sheer number, which likely amounts to >2000 individual miRNAs in human, than by the regulatory effectiveness of single miRNAs. Although, one begins to develop an understanding of the complex mechanisms underlying miRNA-target interactions (MTIs), the overall knowledge of MTI functionality is still rather patchy. In this critical review, we summarize key features of mammalian MTIs. We especially highlight latest insights on (i) the dynamic make-up of miRNA binding sites including non-canonical binding sites, (ii) the cooperativity between miRNA binding sites, (iii) the adaptivity of MTIs through sequence modifications, (iv) the bearing of intra-cellular miRNA localization changes and (v) the role of cell type and cell status specific miRNA interaction partners. The MTI biology is discussed against the background of state-of-the-art approaches with particular emphasis on experimental strategies for evaluating miRNA functionality.

Emerging concepts of miRNA therapeutics: from cells to clinic.

MicroRNAs (miRNAs) are very powerful genetic regulators, as evidenced by the fact that a single miRNA can direct entire cellular pathways via interacting with a broad spectrum of target genes. This property renders miRNAs as highly interesting therapeutic tools to restore cell functions that are altered as part of a disease phenotype. However, this strength of miRNAs is also a weakness because their cellular effects are so numerous that off-target effects can hardly be avoided. In this review, we point out the main challenges and the strategies to specifically address the problems that need to be surmounted in the push toward a therapeutic application of miRNAs. Particular emphasis is given to approaches that have already found their way into clinical studies.

Extracellular vesicles and PD-L1 suppress macrophages, inducing therapy resistance in TP53-deficient B-cell malignancies.

Genetic alterations in the DNA damage response (DDR) pathway are a frequent mechanism of resistance to chemoimmunotherapy (CIT) in B-cell malignancies. We have previously shown that the synergy of CIT relies on secretory crosstalk elicited by chemotherapy between the tumor cells and macrophages. Here, we show that loss of multiple different members of the DDR pathway inhibits macrophage phagocytic capacity in vitro and in vivo. Particularly, loss of TP53 led to decreased phagocytic capacity ex vivo across multiple B-cell malignancies. We demonstrate via in vivo cyclophosphamide treatment using the Eµ-TCL1 mouse model that loss of macrophage phagocytic capacity in Tp53-deleted leukemia is driven by a significant downregulation of a phagocytic transcriptomic signature using small conditional RNA sequencing. By analyzing the tumor B-cell proteome, we identified a TP53-specific upregulation of proteins associated with extracellular vesicles (EVs). We abrogated EV biogenesis in tumor B-cells via clustered regularly interspaced short palindromic repeats (CRISPR)-knockout (KO) of RAB27A and confirmed that the EVs from TP53-deleted lymphoma cells were responsible for the reduced phagocytic capacity and the in vivo CIT resistance. Furthermore, we observed that TP53 loss led to an upregulation of both PD-L1 cell surface expression and secretion of EVs by lymphoma cells. Disruption of EV bound PD-L1 by anti-PD-L1 antibodies or PD-L1 CRISPR-KO improved macrophage phagocytic capacity and in vivo therapy response. Thus, we demonstrate enhanced EV release and increased PD-L1 expression in TP53-deficient B-cell lymphomas as novel mechanisms of macrophage function alteration in CIT resistance. This study indicates the use of checkpoint inhibition in the combination treatment of B-cell malignancies with TP53 loss.

Skin treatment with non-thermal plasma modulates the immune system through miR-223-3p and its target genes.

Non-thermal plasma, a partially ionized gas, holds significant potential for clinical applications, including wound-healing support, oral therapies, and anti-tumour treatments. While its applications showed promising outcomes, the underlying molecular mechanisms remain incompletely understood. We thus apply non-thermal plasma to mouse auricular skin and conducted non-coding RNA sequencing, as well as single-cell blood sequencing. In a time-series analysis (five timepoints spanning 2 hours), we compare the expression of microRNAs in the plasma-treated left ears to the unexposed right ears of the same mice as well as to the ears of unexposed control mice. Our findings indicate specific effects in the treated ears for a set of five miRNAs: mmu-miR-144-5p, mmu-miR-144-3p, mmu-miR-142a-5p, mmu-miR-223-3p, and mmu-miR-451a. Interestingly, mmu-miR-223-3p also exhibits an increase over time in the right non-treated ear of the exposed mice, suggesting systemic effects. Notably, this miRNA, along with mmu-miR-142a-5p and mmu-miR-144-3p, regulates genes and pathways associated with wound healing and tissue regeneration (namely ErbB, FoxO, Hippo, and PI3K-Akt signalling). This co-regulation is particularly remarkable considering the significant seed dissimilarities among the miRNAs. Finally, single-cell sequencing of PBMCs reveals the downregulation of 12 from 15 target genes in B-cells, Cd4+ and Cd8+ T-cells. Collectively, our data provide evidence for a systemic effect of non-thermal plasma.

miRNATissueAtlas2: an update to the human miRNA tissue atlas.

Small non-coding RNAs (sncRNAs) are pervasive regulators of physiological and pathological processes. We previously developed the human miRNA Tissue Atlas, detailing the expression of miRNAs across organs in the human body. Here, we present an updated resource containing sequencing data of 188 tissue samples comprising 21 organ types retrieved from six humans. Sampling the organs from the same bodies minimizes intra-individual variability and facilitates the making of a precise high-resolution body map of the non-coding transcriptome. The data allow shedding light on the organ- and organ system-specificity of piwi-interacting RNAs (piRNAs), transfer RNAs (tRNAs), microRNAs (miRNAs) and other non-coding RNAs. As use case of our resource, we describe the identification of highly specific ncRNAs in different organs. The update also contains 58 samples from six tissues of the Tabula Muris collection, allowing to check if the tissue specificity is evolutionary conserved between Homo sapiens and Mus musculus. The updated resource of 87 252 non-coding RNAs from nine non-coding RNA classes for all organs and organ systems is available online without any restrictions (https://www.ccb.uni-saarland.de/tissueatlas2).

miR-34a-5p as molecular hub of pathomechanisms in Huntington's disease.

BACKGROUND: Although a pivotal role of microRNA (miRNA, miR) in the pathogenesis of Huntington's disease (HD) is increasingly recognized, the molecular functions of miRNAs in the pathomechanisms of HD await further elucidation. One of the miRNAs that have been associated with HD is miR-34a-5p, which was deregulated in the mouse R6/2 model and in human HD brain tissues. METHODS: The aim of our study was to demonstrate interactions between miR-34a-5p and HD associated genes. By computational means we predicted 12 801 potential target genes of miR-34a-5p. An in-silico pathway analysis revealed 22 potential miR-34a-5p target genes in the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway "Huntington's disease". RESULTS: Using our high-throughput miRNA interaction reporter assay (HiTmIR) we identified NDUFA9, TAF4B, NRF1, POLR2J2, DNALI1, HIP1, TGM2 and POLR2G as direct miR-34a-5p target genes. Direct binding of miR-34a-5p to target sites in the 3'UTRs of TAF4B, NDUFA9, HIP1 and NRF1 was verified by a mutagenesis HiTmIR assay and by determining endogenous protein levels for HIP1 and NDUFA9. STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) analysis identified protein-protein interaction networks associated with HD like "Glutamine Receptor Signaling Pathway" and "Calcium Ion Transmembrane Import Into Cytosol". CONCLUSION: Our study demonstrates multiple interactions between miR-34a-5p and HD associated target genes and thereby lays the ground for future therapeutic interventions using this miRNA.

Genome-wide DNA methylome and transcriptome changes induced by inorganic nanoparticles in human kidney cells after chronic exposure.

The unique physicochemical properties make inorganic nanoparticles (INPs) an exciting tool in diagnosis and disease management. However, as INPs are relatively difficult to fully degrade and excrete, their unintended accumulation in the tissue might result in adverse health effects. Herein, we provide a methylome-transcriptome framework for chronic effects of INPs, commonly used in biomedical applications, in human kidney TH-1 cells. Renal clearance is one of the most important routes of nanoparticle excretion; therefore, a detailed evaluation of nanoparticle-mediated nephrotoxicity is an important task. Integrated analysis of methylome and transcriptome changes induced by INPs (PEG-AuNPs, Fe3O4NPs, SiO2NPs, and TiO2NPs) revealed significantly deregulated genes with functional classification in immune response, DNA damage, and cancer-related pathways. Although most deregulated genes were unique to individual INPs, a relatively high proportion of them encoded the transcription factors. Interestingly, FOS hypermethylation inversely correlating with gene expression was associated with all INPs exposures. Our study emphasizes the need for a more comprehensive investigation of INPs' biological safety, especially after chronic exposure.

Dynamic and static circulating cancer microRNA biomarkers - a validation study.

For cancers and other pathologies, early diagnosis remains the most promising path to survival. Profiling of longitudinal cohorts facilitates insights into trajectories of biomarkers. We measured microRNA expression in 240 serum samples from patients with colon, lung, and breast cancer and from cancer-free controls. Each patient provided at least two serum samples, one prior to diagnosis and one following diagnosis. The median time interval between the samples was 11.6 years. Using computational models, we evaluated the circulating profiles of 21 microRNAs. The analysis yielded two sets of biomarkers, static ones that show an absolute difference between certain cancer types and controls and dynamic ones where the level over time provided higher diagnostic information content. In the first group, miR-99a-5p stands out for all three cancer types. In the second group, miR-155-5p allows to predict lung cancers and colon cancers. Classification in samples from cancer and non-cancer patients using gradient boosted trees reached an average accuracy of 79.9%. The results suggest that individual change over time or an absolute value at one time point may predict a disease with high specificity and sensitivity.

Ageing-associated small RNA cargo of extracellular vesicles.

Previous work on murine models and humans demonstrated global as well as tissue-specific molecular ageing trajectories of RNAs. Extracellular vesicles (EVs) are membrane vesicles mediating the horizontal transfer of genetic information between different tissues. We sequenced small regulatory RNAs (sncRNAs) in two mouse plasma fractions at five time points across the lifespan from 2-18 months: (1) sncRNAs that are free-circulating (fc-RNA) and (2) sncRNAs bound outside or inside EVs (EV-RNA). Different sncRNA classes exhibit unique ageing patterns that vary between the fcRNA and EV-RNA fractions. While tRNAs showed the highest correlation with ageing in both fractions, rRNAs exhibited inverse correlation trajectories between the EV- and fc-fractions. For miRNAs, the EV-RNA fraction was exceptionally strongly associated with ageing, especially the miR-29 family in adipose tissues. Sequencing of sncRNAs and coding genes in fat tissue of an independent cohort of aged mice up to 27 months highlighted the pivotal role of miR-29a-3p and miR-29b-3p in ageing-related gene regulation that we validated in a third cohort by RT-qPCR.

Outside the limit: questioning the distance restrictions for cooperative miRNA binding sites.

Among the concepts in biology that are widely taken granted is a potentiated cooperative effect of multiple miRNAs on the same target. This strong hypothesis contrasts insufficient experimental evidence. The quantity as well as the quality of required side constraints of cooperative binding remain largely hidden. For miR-21-5p and miR-155-5p, two commonly investigated regulators across diseases, we selected 15 joint target genes. These were chosen to represent various neighboring 3'UTR binding site constellations, partially exceeding the distance rules that have been established for over a decade. We identified different cooperative scenarios with the binding of one miRNA enhancing the binding effects of the other miRNA and vice versa. Using both, reporter assays and whole proteome analyses, we observed these cooperative miRNA effects for genes that bear 3'UTR binding sites at distances greater than the previously defined limits. Astonishingly, the experiments provide even stronger evidence for cooperative miRNA effects than originally postulated. In the light of these findings the definition of targetomes specified for single miRNAs need to be refined by a concept that acknowledges the cooperative effects of miRNAs.

miRTargetLink 2.0-interactive miRNA target gene and target pathway networks.

Which genes, gene sets or pathways are regulated by certain miRNAs? Which miRNAs regulate a particular target gene or target pathway in a certain physiological context? Answering such common research questions can be time consuming and labor intensive. Especially for researchers without computational experience, the integration of different data sources, selection of the right parameters and concise visualization can be demanding. A comprehensive analysis should be central to present adequate answers to complex biological questions. With miRTargetLink 2.0, we develop an all-in-one solution for human, mouse and rat miRNA networks. Users input in the unidirectional search mode either a single gene, gene set or gene pathway, alternatively a single miRNA, a set of miRNAs or an miRNA pathway. Moreover, genes and miRNAs can jointly be provided to the tool in the bidirectional search mode. For the selected entities, interaction graphs are generated from different data sources and dynamically presented. Connected application programming interfaces (APIs) to the tailored enrichment tools miEAA and GeneTrail facilitate downstream analysis of pathways and context-annotated categories of network nodes. MiRTargetLink 2.0 is freely accessible at https://www.ccb.uni-saarland.de/mirtargetlink2.

CoolMPS: evaluation of antibody labeling based massively parallel non-coding RNA sequencing.

Results of massive parallel sequencing-by-synthesis vary depending on the sequencing approach. CoolMPS™ is a new sequencing chemistry that incorporates bases by labeled antibodies. To evaluate the performance, we sequenced 240 human non-coding RNA samples (dementia patients and controls) with and without CoolMPS. The Q30 value as indicator of the per base sequencing quality increased from 91.8 to 94%. The higher quality was reached across the whole read length. Likewise, the percentage of reads mapping to the human genome increased from 84.9 to 86.2%. For both technologies, we computed similar distributions between different RNA classes (miRNA, piRNA, tRNA, snoRNA and yRNA) and within the classes. While standard sequencing-by-synthesis allowed to recover more annotated miRNAs, CoolMPS yielded more novel miRNAs. The correlation between the two methods was 0.97. Evaluating the diagnostic performance, we observed lower minimal P-values for CoolMPS (adjusted P-value of 0.0006 versus 0.0004) and larger effect sizes (Cohen's d of 0.878 versus 0.9). Validating 19 miRNAs resulted in a correlation of 0.852 between CoolMPS and reverse transcriptase-quantitative polymerase chain reaction. Comparison to data generated with Illumina technology confirmed a known shift in the overall RNA composition. With CoolMPS we evaluated a novel sequencing-by-synthesis technology showing high performance for the analysis of non-coding RNAs.

Validation of human microRNA target pathways enables evaluation of target prediction tools.

MicroRNAs are regulators of gene expression. A wide-spread, yet not validated, assumption is that the targetome of miRNAs is non-randomly distributed across the transcriptome and that targets share functional pathways. We developed a computational and experimental strategy termed high-throughput miRNA interaction reporter assay (HiTmIR) to facilitate the validation of target pathways. First, targets and target pathways are predicted and prioritized by computational means to increase the specificity and positive predictive value. Second, the novel webtool miRTaH facilitates guided designs of reporter assay constructs at scale. Third, automated and standardized reporter assays are performed. We evaluated HiTmIR using miR-34a-5p, for which TNF- and TGFB-signaling, and Parkinson's Disease (PD)-related categories were identified and repeated the pipeline for miR-7-5p. HiTmIR validated 58.9% of the target genes for miR-34a-5p and 46.7% for miR-7-5p. We confirmed the targeting by measuring the endogenous protein levels of targets in a neuronal cell model. The standardized positive and negative targets are collected in the new miRATBase database, representing a resource for training, or benchmarking new target predictors. Applied to 88 target predictors with different confidence scores, TargetScan 7.2 and miRanda outperformed other tools. Our experiments demonstrate the efficiency of HiTmIR and provide evidence for an orchestrated miRNA-gene targeting.

The DGCR8 E518K mutation found in Wilms tumors leads to a partial miRNA processing defect that alters gene expression patterns and biological processes.

Wilms tumor (WT) is the most common renal tumor in childhood. We and others have previously identified oncogenic driver mutations affecting the microprocessor genes DROSHA and DGCR8 that lead to altered miRNA expression patterns. In the case of DGCR8, a single recurrent hotspot mutation (E518K) was found in the RNA binding domain. To functionally assess this mutation in vitro, we generated mouse Dgcr8-KO embryonic stem cell (mESC) lines with an inducible expression of wild-type or mutant DGCR8, mirroring the hemizygous mutant expression seen in WT. RNA-seq analysis revealed significant differences of miRNA expression profiles in DGCR8-E518K compared with DGCR8-wild-type mESCs. The E518K mutation only led to a partial rescue of the reported miRNA processing defect in Dgcr8-KO, with selectively reduced expression of numerous canonical miRNAs. Nevertheless, DGCR8-E518K retained significant activity given its ability to still process many miRNAs. Subsequent to altered miRNA levels, the expression of mRNA targets was likewise changed. Functional assays showed that DGCR8-E518K cells still have a partial proliferation and differentiation defect but were able to rescue critical biological processes in embryoid body development. The stem cell program could be shut down and all three germ layers were formed. These findings suggest that the E518K mutation leads to a partial reduction of microprocessor activity and altered specificity with selective impairment only in certain developmental contexts, apparently including nephrogenesis.

What's the target: understanding two decades of in silico microRNA-target prediction.

MOTIVATION: Since the initial discovery of microRNAs as post-transcriptional, regulatory key players in the 1990s, a total number of $2656$ mature microRNAs have been publicly described for Homo sapiens. As discovery of new miRNAs is still on-going, target identification remains to be an essential and challenging step preceding functional annotation analysis. One key challenge for researchers seems to be the selection of the most appropriate tool out of the larger multiverse of published solutions for a given research study set-up. RESULTS: In this review we collectively describe the field of in silico target prediction in the course of time and point out long withstanding principles as well as recent developments. By compiling a catalog of characteristics about the 98 prediction methods and identifying common and exclusive traits, we signpost a simplified mechanism to address the problem of application selection. Going further we devised interpretation strategies for common types of output as generated by frequently used computational methods. To this end, our work specifically aims to make prospective users aware of common mistakes and practical questions that arise during the application of target prediction tools. AVAILABILITY: An interactive implementation of our recommendations including materials shown in the manuscript is freely available at https://www.ccb.uni-saarland.de/mtguide.

miRSwitch: detecting microRNA arm shift and switch events.

Arm selection, the preferential expression of a 3' or 5' mature microRNA (miRNA), is a highly dynamic and tissue-specific process. Time-dependent expression shifts or switches between the arms are also relevant for human diseases. We present miRSwitch, a web server to facilitate the analysis and interpretation of arm selection events. Our species-independent tool evaluates pre-processed small non-coding RNA sequencing (sncRNA-seq) data, i.e. expression matrices or output files from miRNA quantification tools (miRDeep2, miRMaster, sRNAbench). miRSwitch highlights potential changes in the distribution of mature miRNAs from the same precursor. Group comparisons from one or several user-provided annotations (e.g. disease states) are possible. Results can be dynamically adjusted by choosing from a continuous range of highly specific to very sensitive parameters. Users can compare potential arm shifts in the provided data to a human reference map of pre-computed arm shift frequencies. We created this map from 46 tissues and 30 521 samples. As case studies we present novel arm shift information in a Alzheimer's disease biomarker data set and from a comparison of tissues in Homo sapiens and Mus musculus. In summary, miRSwitch offers a broad range of customized arm switch analyses along with comprehensive visualizations, and is freely available at: https://www.ccb.uni-saarland.de/mirswitch/.

miRPathDB 2.0: a novel release of the miRNA Pathway Dictionary Database.

Since the initial release of miRPathDB, tremendous progress has been made in the field of microRNA (miRNA) research. New miRNA reference databases have emerged, a vast amount of new miRNA candidates has been discovered and the number of experimentally validated target genes has increased considerably. Hence, the demand for a major upgrade of miRPathDB, including extended analysis functionality and intuitive visualizations of query results has emerged. Here, we present the novel release 2.0 of the miRNA Pathway Dictionary Database (miRPathDB) that is freely accessible at https://mpd.bioinf.uni-sb.de/. miRPathDB 2.0 comes with a ten-fold increase of pre-processed data. In total, the updated database provides putative associations between 27 452 (candidate) miRNAs, 28 352 targets and 16 833 pathways for Homo sapiens, as well as interactions of 1978 miRNAs, 24 898 targets and 6511 functional categories for Mus musculus. Additionally, we analyzed publications citing miRPathDB to identify common use-cases and further extensions. Based on this evaluation, we added new functionality for interactive visualizations and down-stream analyses of bulk queries. In summary, the updated version of miRPathDB, with its new custom-tailored features, is one of the most comprehensive and advanced resources for miRNAs and their target pathways.

miEAA 2.0: integrating multi-species microRNA enrichment analysis and workflow management systems.

Gene set enrichment analysis has become one of the most frequently used applications in molecular biology research. Originally developed for gene sets, the same statistical principles are now available for all omics types. In 2016, we published the miRNA enrichment analysis and annotation tool (miEAA) for human precursor and mature miRNAs. Here, we present miEAA 2.0, supporting miRNA input from ten frequently investigated organisms. To facilitate inclusion of miEAA in workflow systems, we implemented an Application Programming Interface (API). Users can perform miRNA set enrichment analysis using either the web-interface, a dedicated Python package, or custom remote clients. Moreover, the number of category sets was raised by an order of magnitude. We implemented novel categories like annotation confidence level or localisation in biological compartments. In combination with the miRBase miRNA-version and miRNA-to-precursor converters, miEAA supports research settings where older releases of miRBase are in use. The web server also offers novel comprehensive visualizations such as heatmaps and running sum curves with background distributions. We demonstrate the new features with case studies for human kidney cancer, a biomarker study on Parkinson's disease from the PPMI cohort, and a mouse model for breast cancer. The tool is freely accessible at: https://www.ccb.uni-saarland.de/mieaa2.

Quantitative and time-resolved miRNA pattern of early human T cell activation.

T cells are central to the immune response against various pathogens and cancer cells. Complex networks of transcriptional and post-transcriptional regulators, including microRNAs (miRNAs), coordinate the T cell activation process. Available miRNA datasets, however, do not sufficiently dissolve the dynamic changes of miRNA controlled networks upon T cell activation. Here, we established a quantitative and time-resolved expression pattern for the entire miRNome over a period of 24 h upon human T-cell activation. Based on our time-resolved datasets, we identified central miRNAs and specified common miRNA expression profiles. We found the most prominent quantitative expression changes for miR-155-5p with a range from initially 40 molecules/cell to 1600 molecules/cell upon T-cell activation. We established a comprehensive dynamic regulatory network of both the up- and downstream regulation of miR-155. Upstream, we highlight IRF4 and its complexes with SPI1 and BATF as central for the transcriptional regulation of miR-155. Downstream of miR-155-5p, we verified 17 of its target genes by the time-resolved data recorded after T cell activation. Our data provide comprehensive insights into the range of stimulus induced miRNA abundance changes and lay the ground to identify efficient points of intervention for modifying the T cell response.