The MicroRNA Family Gets Wider: The IsomiRs Classification and Role

miRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different microRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression methods developed for mRNA-seq data. We establish miRglmm, a differential expression method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current differential expression methods in estimating differential expression for miRNA, whether or not there is differential isomiR usage, and simultaneously provides estimates of isomiR-level differential expression.

Product-induced catalytic amplification strategy based on DNA tetrahedron for detection of miRNA-21 in colorectal cancer

A product-induced catalytic amplification (PICA) strategy had been developed for miRNA-21 detection based on DNA tetrahedron module (DTM). The produced DNA fragment could open hairpin structure and increase the concentration of catalyst, accelerating the circular cleavage reaction on DTM by DNAzyme cleavage. The continuously cleavage of DNAzyme on DTM resulted the greatly enhancement of signal. A favorable linear range was achieved from 20 pM to 5 nM with a limit of detection of 7 pM. Furthermore, through the implementation of the PICA strategy, the overall reaction time experienced a noticeable decrease to 30 min. The assessments of the amplification rate and kinetic constant of the PICA strategy were also conducted. These results highlighted the promising potential of the PICA strategy for practical utilization in serum samples.

Comprehensive sequencing profile and functional analysis of IsomiRs in human pancreatic islets and beta cells

AIMS/HYPOTHESIS

MiRNAs regulate gene expression, influencing beta cell function and pathways. Isoforms of miRNA (isomiRs), sequence variants of miRNAs with post-transcriptional modifications, exhibit cell-type-specific expression and functions. Despite their biological significance, a comprehensive isomiR profile in human pancreatic islets and beta cells remains unexplored. This study aims to profile isomiR expression in four beta cell sources: (1) laser capture microdissected human islets (LCM-HI); (2) collagenase-isolated human islets (CI-HI); (3) sorted beta cells; and (4) the EndoC-βH1 beta cell line, and to investigate their potential role in beta cell function.

METHODS

Small RNA-seq and/or small RNA dataset analysis was conducted on human pancreatic islets and beta cells. Data were processed using the sRNAbench bioinformatics pipeline to classify isomiRs based on sequence variations. A beta cell-specific isomiR signature was identified via cross-validation across datasets. Correlations between LCM-HI isomiR expression and in vivo clinical parameters were analysed using regression models. Functional validation of isomiR-411-5p-Ext5p(+1) was performed via overexpression in EndoC-βH1 cells and CI-HI, followed by glucose-stimulated insulin secretion (GSIS) assays and/or transcriptomic analysis.

RESULTS

IsomiRs constituted 59.2 ± 1.9% (LCM-HI), 59.6 ± 2.4% (CI-HI), 42.3 ± 7.2% (sorted beta cells) and 43.8 ± 1.2% (EndoC-βH1) of total miRNA reads (data represented as mean ± SD), with 3' end trimming (Trim3p) being the predominant modification. A beta cell-specific isomiR signature of 30 sequences was identified, with isomiR-411-5p-Ext5p(+1) showing a significant inverse correlation with basal insulin secretion (p=0.0009, partial R2=0.68) and total insulin secretion (p=0.005, partial R2=0.54). Overexpression of isomiR-411-5p-Ext5p(+1), but not of its canonical counterpart, importantly reduced GSIS by 51% ( ± 15.2%; mean ± SD) (p=0.01) in EndoC-βH1 cells. Transcriptomic analysis performed in EndoC-βH1 cells and CI-HI identified 47 genes significantly downregulated by isomiR-411-5p-Ext5p(+1) (false discovery rate [FDR]<0.05) but not by the canonical miRNA, with enriched pathways related to Golgi vesicle biogenesis (FDR=0.017) and trans-Golgi vesicle budding (FDR=0.018). TargetScan analysis confirmed seed sequence-dependent target specificity for 81 genes uniquely regulated by the isomiR (p=1.1 × 10⁻⁹).

CONCLUSIONS/INTERPRETATION

This study provides the first comprehensive isomiR profiling in human islets and beta cells, revealing their substantial contribution to miRNA regulation. IsomiR-411-5p-Ext5p(+1) emerges as a distinct key modulator of insulin secretion and granule dynamics in beta cells. These findings highlight isomiRs as potential biomarkers and therapeutic targets for diabetes, warranting further exploration of their roles in beta cell biology.

Combinatorial Analysis of miRNAs and tRNA Fragments as Potential Biomarkers for Cancer Patients in Liquid Biopsies

Background: Liquid biopsy has gained significant attention as a non-invasive method for cancer detection and monitoring. IsomiRs and tRNA-derived fragments (tRFs) are small non-coding RNAs that arise from non-canonical microRNA (miRNAs) processing and the cleavage of tRNAs, respectively. These small non-coding RNAs have emerged as pro-mising cancer biomarkers, and their distinct expression patterns highlight the need for further exploration of their roles in cancer research. Methods: In this study, we investigated the differential expression profiles of miRNAs, isomiRs, and tRFs in plasma extracellular vesicles (EVs) from colorectal and prostate cancer patients compared to healthy controls. Subsequently, a combinatorial analysis using the CombiROC package was performed to identify a panel of biomarkers with optimal diagnostic accuracy. Results: Our results demonstrate that a combination of miRNAs, isomiRs, and tRFs can effectively di- stinguish cancer patients from healthy controls, achieving accuracy and an area under the curve (AUC) of approximately 80%. Conclusions: These findings highlight the potential of a combinatorial approach to small RNA analysis in liquid biopsies for improved cancer diagnosis and management.

Revealing microRNA regulation in single cells

MicroRNAs (miRNAs) are key regulators of gene expression and control cellular functions in physiological and pathophysiological states. miRNAs play important roles in disease, stress, and development, and are now being investigated for therapeutic approaches. Alternative processing of miRNAs during biogenesis results in the generation of miRNA isoforms (isomiRs) which further diversify miRNA gene regulation. Single-cell RNA-sequencing (scsRNA-seq) technologies, together with computational strategies, enable exploration of miRNAs, isomiRs, and interacting RNAs at the cellular level. By integration with other miRNA-associated single-cell modalities, miRNA roles can be resolved at different stages of processing and regulation. In this review we discuss (i) single-cell experimental assays that measure miRNA and isomiR abundances, and (ii) computational methods for their analysis to investigate the mechanisms of miRNA biogenesis and post-transcriptional regulation.

ncRNA Editing: Functional Characterization and Computational Resources

Non-coding RNAs (ncRNAs) play crucial roles in gene expression regulation, translation, and disease development, including cancer. They are classified by size in short and long non-coding RNAs. This chapter focuses on the functional implications of adenosine-to-inosine (A-to-I) RNA editing in both short (e.g., miRNAs) and long ncRNAs. RNA editing dynamically alters the sequence and structure of primary transcripts, impacting ncRNA biogenesis and function. Notable findings include the role of miRNA editing in promoting glioblastoma invasiveness, characterizing RNA editing hotspots across cancers, and its implications in thyroid cancer and ischemia. This chapter also highlights bioinformatics resources and next-generation sequencing (NGS) technologies that enable comprehensive ncRNAome studies and genome-wide RNA editing detection. Dysregulation of RNA editing machinery has been linked to various human diseases, emphasizing the potential of RNA editing as a biomarker and therapeutic target. This overview integrates current knowledge and computational tools for studying ncRNA editing, providing insights into its biological significance and clinical applications.

Comprehensive profiling of small RNAs and their changes and linkages to mRNAs in schizophrenia and bipolar disorder

We investigated small non-coding RNAs (sncRNAs) from the prefrontal cortex of 93 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) and 77 controls. We uncovered recurring complex sncRNA profiles, with 98% of all sncRNAs being accounted for by miRNA isoforms (60.6%), tRNA-derived fragments (17.8%), rRNA-derived fragments (11.4%), and Y RNA-derived fragments (8.3%). In SCZ, 15% of all sncRNAs exhibit statistically significant changes in their abundance. In BD, the fold changes (FCs) are highly correlated with those in SCZ but less acute. Non-templated nucleotide additions to the 3´-ends of many miRNA isoforms determine their FC independently of miRNA identity or genomic locus of origin. In both SCZ and BD, disease- and age-associated sncRNAs and mRNAs reveal accelerated aging. Co-expression modules between sncRNAs and mRNAs align with the polarities of SCZ changes and implicate sncRNAs in critical processes, including synaptic signaling, neurogenesis, memory, behavior, and cognition.

AmiCa: Atlas of miRNA-gene correlations in cancer

The increasing availability of RNA sequencing data has opened up numerous opportunities to analyze various RNA interactions, including microRNA-target interactions (MTIs). In response to the necessity for a specialized tool to study MTIs in cancer and normal tissues, we developed AmiCa (https://amica.omics.si/), a web server designed for comprehensive analysis of mature microRNA (miRNA) and gene expression in 32 cancer types. Data from 9498 tumor samples and 626 normal samples from The Cancer Genome Atlas were obtained through the Genomic Data Commons and used to calculate differential expression and miRNA-target gene (MTI) correlations. AmiCa provides data on differential expression of miRNAs/genes for cancers for which normal tissue samples were available. In addition, the server calculates and presents correlations separately for tumor and normal samples for cancers for which normal samples are available. Furthermore, it enables the exploration of miRNA/gene expression in all cancer types with different miRNA/gene expression. In addition, AmiCa includes a ranking system for genes and miRNAs that can be used to identify those that are particularly highly expressed in certain cancers compared to other cancers, facilitating targeted and cancer-specific research. Finally, the functionality of AmiCa is illustrated by two case studies.

Distinct brain and lung endothelial miRNA/mRNA profiles after exposure to Plasmodium falciparum-infected red blood cells

MicroRNAs (miRNAs) control 60% of genes expressed in the human body, but their role in malaria pathogenesis is incompletely understood. Here, we demonstrate cell type-specific alterations to the miRNA profiles during the early response to malaria infection in brain and lung endothelial cells (ECs). In brain ECs, incubation with Plasmodium falciparum-infected red blood cells in the ring stage (iRBCs) most significantly affected endocytosis-related miRNAs and mRNAs. Contrastingly, in lung ECs, iRBCs altered electron transport chain-related miRNAs and mRNAs. We present a dataset of inherent differences between microRNA profiles in brain and lung ECs and their extracellular vesicles (EVs). We demonstrated that shear stress affected multiple pathways in brain ECs, which were controlled by numerous human miRNAs. Together, these findings indicate that host miRNAs respond to parasite exposure, accompanied by stimulation of downstream signaling pathways within the ECs. Therefore, we consider miRNAs the initial spark for early host-parasite interaction events.

A three-layer perspective on miRNA regulation in β cell inflammation

MicroRNAs (miRNAs) are noncoding RNA molecules that regulate gene expression post-transcriptionally and influence numerous biological processes. Aberrant miRNA expression is linked to diseases such as diabetes mellitus; indeed, miRNAs regulate pancreatic islet inflammation in both type 1 (T1D) and type 2 diabetes (T2D). Traditionally, miRNA research has focused on canonical sequences and offers a two-layer view - from expression to function. However, advances in RNA sequencing have revealed miRNA variants, called isomiRs, that arise from alternative processing or modifications of canonical sequences. This introduces a three-layer view - from expression, through sequence modifications, to function. We discuss the potential link between cellular stresses and isomiR biogenesis, and how this association could improve our knowledge of islet inflammation and dysfunction.

Quantitative design of cell type-specific mRNA stability from microRNA expression data

Limiting expression to target cell types is a longstanding goal in gene therapy, which could be met by sensing endogenous microRNA. However, an unclear association between microRNA expression and activity currently hampers such an approach. Here, we probe this relationship by measuring the stability of synthetic microRNA-responsive 3'UTRs across 10 cell lines in a library format. By systematically addressing biases in microRNA expression data and confounding factors such as microRNA crosstalk, we demonstrate that a straightforward model can quantitatively predict reporter stability purely from expression data. We use this model to design constructs with previously unattainable response patterns across our cell lines. The rules we derive for microRNA expression data selection and processing should apply to microRNA- responsive devices for any environment with available expression data.

An atlas of small non-coding RNAs in human preimplantation development

Understanding the molecular circuitries that govern early embryogenesis is important, yet our knowledge of these in human preimplantation development remains limited. Small non-coding RNAs (sncRNAs) can regulate gene expression and thus impact blastocyst formation, however, the expression of specific biotypes and their dynamics during preimplantation development remains unknown. Here we identify the abundance of and kinetics of piRNA, rRNA, snoRNA, tRNA, and miRNA from embryonic day (E)3-7 and isolate specific miRNAs and snoRNAs of particular importance in blastocyst formation and pluripotency. These sncRNAs correspond to specific genomic hotspots: an enrichment of the chromosome 19 miRNA cluster (C19MC) in the trophectoderm (TE), and the chromosome 14 miRNA cluster (C14MC) and MEG8-related snoRNAs in the inner cell mass (ICM), which may serve as 'master regulators' of potency and lineage. Additionally, we observe a developmental transition with 21 isomiRs and in tRNA fragment (tRF) codon usage and identify two novel miRNAs. Our analysis provides a comprehensive measure of sncRNA biotypes and their corresponding dynamics throughout human preimplantation development, providing an extensive resource. Better understanding the sncRNA regulatory programmes in human embryogenesis will inform strategies to improve embryo development and outcomes of assisted reproductive technologies. We anticipate broad usage of our data as a resource for studies aimed at understanding embryogenesis, optimising stem cell-based models, assisted reproductive technology, and stem cell biology.

Post-Transcriptional Modifications to miRNAs Undergo Widespread Alterations, Creating a Unique Lung Adenocarcinoma IsomiRome

BACKGROUND

MicroRNAs (miRNAs) modulate the expression of oncogenes and tumor suppressor genes, functioning as significant epigenetic regulators in cancer. IsomiRs are miRNA molecules that have undergone small modifications during miRNA processing. These modifications can alter an isomiR's binding stability with mRNA targets, and certain isomiRs have been implicated in the development of specific cancers. Still, the isomiRomes of many tissues, including the lung, have not been characterized; Methods: In this study, we analyzed small RNA sequencing data for three cohorts of lung adenocarcinoma (LUAD) and adult non-malignant lung (ANL) samples.

RESULTS

We quantified isomiR expression and found 16 A-to-I edited isomiRs expressed in multiple cohorts, as well as 213 5' isomiRs, 128 3' adenylated isomiRs, and 100 3' uridylated isomiRs. Rates of A-to-I editing at editing hotspots correlated with mRNA expression of the editing enzymes ADAR and ADARB1, which were both observed to be deregulated in LUAD. LUAD samples displayed lower overall rates of A-to-I editing and 3' adenylation than ANL samples. Support vector machines and random forest models were trained on one cohort to distinguish ANL and stage I/II LUAD samples using reads per million (RPM) and frequency data for different types of isomiRs. Models trained on A-to-I editing rates at editing hotspots displayed high accuracy when tested on the other two cohorts and compared favorably to classifiers trained on miRNA expression alone; Conclusions: We have identified isomiRs in the human lung and found that their expression differs between non-malignant and tumor tissues, suggesting they hold potential as cancer biomarkers.

Longitudinal associations between microRNAs and weight in the diabetes prevention program

Objective

Circulating microRNAs show cross-sectional associations with overweight and obesity. Few studies provided data to differentiate between a snapshot perspective on these associations versus how microRNAs characterize prodromal risk from disease pathology and complications. This study assessed longitudinal relationships between circulating microRNAs and weight at multiple time-points in the Diabetes Prevention Program trial.

Research design and methods

A subset of participants (n=150) from the Diabetes Prevention Program were included. MicroRNAs were measured from banked plasma using a Fireplex Assay. We used generalized linear mixed models to evaluate relationships between microRNAs and changes in weight at baseline, year-1, and year-2. Logistic regression was used to evaluate whether microRNAs at baseline were associated with weight change after 2 years.

Results

In fully adjusted models that included relevant covariates, seven miRs (i.e., miR-126, miR-15a, miR-192, miR-23a, and miR-27a) were statistically associated with weight over 2 years. MiR-197 and miR-320a remained significant after adjustment for multiple comparisons. Baseline levels of let-7f, miR-17, and miR-320c were significantly associated with 3% weight loss after 2 years in fully adjusted models.

Discussion

This study provided evidence for longitudinal relationships between circulating microRNAs and weight. Because microRNAs characterize the combined effects of genetic determinants and responses to behavioral determinants, they may provide insights about the etiology of overweight and obesity in the context or risk for common, complex diseases. Additional studies are needed to validate the potential genes and biological pathways that might be targeted by these microRNA biomarkers and have mechanistic implications for weight loss and disease prevention.

The roles and mechanisms of coding and noncoding RNA variations in cancer

Functional variations in coding and noncoding RNAs are crucial in tumorigenesis, with cancer-specific alterations often resulting from chemical modifications and posttranscriptional processes mediated by enzymes. These RNA variations have been linked to tumor cell proliferation, growth, metastasis, and drug resistance and are valuable for identifying diagnostic or prognostic cancer biomarkers. The diversity of posttranscriptional RNA modifications, such as splicing, polyadenylation, methylation, and editing, is particularly significant due to their prevalence and impact on cancer progression. Additionally, other modifications, including RNA acetylation, circularization, miRNA isomerization, and pseudouridination, are recognized as key contributors to cancer development. Understanding the mechanisms underlying these RNA modifications in cancer can enhance our knowledge of cancer biology and facilitate the development of innovative therapeutic strategies. Targeting these RNA modifications and their regulatory enzymes may pave the way for novel RNA-based therapies, enabling tailored interventions for specific cancer subtypes. This review provides a comprehensive overview of the roles and mechanisms of various coding and noncoding RNA modifications in cancer progression and highlights recent advancements in RNA-based therapeutic applications.

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.

A systematic review of non-coding RNA therapeutics in early clinical trials: a new perspective against cancer

Targeting non-coding RNAs (ncRNAs), including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), has recently emerged as a promising strategy for treating malignancies and other diseases. In recent years, the development of ncRNA-based therapeutics for targeting protein-coding and non-coding genes has also gained momentum. This review systematically examines ongoing and completed clinical trials to provide a comprehensive overview of the emerging landscape of ncRNA-based therapeutics. Significant efforts have been made to advance ncRNA therapeutics to early clinical studies. The most advanced trials have been conducted with small interfering RNAs (siRNAs), miRNA replacement using nanovector-entrapped miRNA mimics, or miRNA silencing by antisense oligonucleotides. While siRNA-based therapeutics have already received FDA approval, miRNA mimics, inhibitors, and lncRNA-based therapeutics are still under evaluation in preclinical and early clinical studies. We critically discuss the rationale and methodologies of ncRNA targeting strategies to illustrate this rapidly evolving field.

miRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different miRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression (DE) methods developed for mRNA-seq data. We establish miRglmm, a DE method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current DE methods in estimating DE for miRNA, whether or not there is significant isomiR variability, and simultaneously provides estimates of isomiR-level DE.

Short 2'-O-methyl/LNA oligomers as highly-selective inhibitors of miRNA production in vitro and in vivo

MicroRNAs (miRNAs) that share identical or near-identical sequences constitute miRNA families and are predicted to act redundantly. Yet recent evidence suggests that members of the same miRNA family with high sequence similarity might have different roles and that this functional divergence might be rooted in their precursors' sequence. Current knock-down strategies such as antisense oligonucleotides (ASOs) or miRNA sponges cannot distinguish between identical or near identical miRNAs originating from different precursors to allow exploring unique functions of these miRNAs. We here develop a novel strategy based on short 2'-OMe/LNA-modified oligonucleotides to selectively target specific precursor molecules and ablate the production of individual members of miRNA families in vitro and in vivo. Leveraging the highly conserved Xenopus miR-181a family as proof-of-concept, we demonstrate that 2'-OMe/LNA-ASOs targeting the apical region of pre-miRNAs achieve precursor-selective inhibition of mature miRNA-5p production. Furthermore, we extend the applicability of our approach to the human miR-16 family, illustrating its universality in targeting precursors generating identical miRNAs. Overall, our strategy enables efficient manipulation of miRNA expression, offering a powerful tool to dissect the functions of identical or highly similar miRNAs derived from different precursors within miRNA families.

Longitudinal Associations Between MicroRNAs and Weight in the Diabetes Prevention Program

OBJECTIVE

Circulating microRNAs show cross-sectional associations with overweight and obesity. Few studies provided data to differentiate between a snapshot perspective on these associations versus how microRNAs characterize prodromal risk from disease pathology and complications. This study assessed longitudinal relationships between circulating microRNAs and weight at multiple time-points in the Diabetes Prevention Program trial.

RESEARCH DESIGN AND METHODS

A subset of participants (n=150) from the Diabetes Prevention Program were included. MicroRNAs were measured from banked plasma using a Fireplex Assay. We used generalized linear mixed models to evaluate relationships between microRNAs and changes in weight at baseline, year-1, and year-2. Logistic regression was used to evaluate whether microRNAs at baseline were associated with weight change after 2 years.

RESULTS

In fully adjusted models that included relevant covariates, seven miRs (i.e., miR-126, miR-15a, miR-192, miR-23a, and miR-27a) were statistically associated with weight over 2 years. MiR-197 and miR-320a remained significant after adjustment for multiple comparisons. Baseline levels of let-7f, miR-17, and miR-320c were significantly associated with 3% weight loss after 2 years in fully adjusted models.

DISCUSSION

This study provided evidence for longitudinal relationships between circulating microRNAs and weight. Because microRNAs characterize the combined effects of genetic determinants and responses to behavioral determinants, they may provide insights about the etiology of overweight and obesity in the context or risk for common, complex diseases. Additional studies are needed to validate the potential genes and biological pathways that might be targeted by these microRNA biomarkers and have mechanistic implications for weight loss and disease prevention.

MicroRNA-mediated regulation of nonsense-mediated mRNA decay factors: Insights into microRNA prediction tools and profiling techniques

Nonsense-mediated mRNA decay (NMD) stands out as a prominent RNA surveillance mechanism within eukaryotes, meticulously overseeing both RNA abundance and integrity by eliminating aberrant transcripts. These defective transcripts are discerned through the concerted efforts of translating ribosomes, eukaryotic release factors (eRFs), and trans-acting NMD factors, with Up-Frameshift 3 (UPF3) serving as a noteworthy component. Remarkably, in humans, UPF3 exists in two paralogous forms, UPF3A (UPF3) and UPF3B (UPF3X). Beyond its role in quality control, UPF3 wields significant influence over critical cellular processes, including neural development, synaptic plasticity, and axon guidance. However, the precise regulatory mechanisms governing UPF3 remain elusive. MicroRNAs (miRNAs) emerge as pivotal post-transcriptional gene regulators, exerting substantial impact on diverse pathological and physiological pathways. This comprehensive review encapsulates our current understanding of the intricate regulatory nexus between NMD and miRNAs, with particular emphasis on the essential role played by UPF3B in neurodevelopment. Additionally, we bring out the significance of the 3'-untranslated region (3'-UTR) as the molecular bridge connecting NMD and miRNA-mediated gene regulation. Furthermore, we provide an in-depth exploration of diverse computational tools tailored for the prediction of potential miRNA targets. To complement these computational approaches, we delineate experimental techniques designed to validate predicted miRNA-mRNA interactions, empowering readers with the knowledge necessary to select the most appropriate methodology for their specific research objectives.

RNA-RNA interactions between respiratory syncytial virus and miR-26 and miR-27 are associated with regulation of cell cycle and antiviral immunity

microRNAs (miRNAs) regulate nearly all physiological processes but our understanding of exactly how they function remains incomplete, particularly in the context of viral infections. Here, we adapt a biochemical method (CLEAR-CLIP) and analysis pipeline to identify targets of miRNAs in lung cells infected with Respiratory syncytial virus (RSV). We show that RSV binds directly to miR-26 and miR-27 through seed pairing and demonstrate that these miRNAs target distinct gene networks associated with cell cycle and metabolism (miR-27) and antiviral immunity (miR-26). Many of the targets are de-repressed upon infection and we show that the miR-27 targets most sensitive to miRNA inhibition are those associated with cell cycle. Finally, we demonstrate that high confidence chimeras map to long noncoding RNAs (lncRNAs) and pseudogenes in transcriptional regulatory regions. We validate that a proportion of miR-27 and Argonaute 2 (AGO2) is nuclear and identify a long non-coding RNA (lncRNA) as a miR-27 target that is linked to transcriptional regulation of nearby genes. This work expands the target networks of miR-26 and miR-27 to include direct interactions with RSV and lncRNAs and implicate these miRNAs in regulation of key genes that impact the viral life cycle associated with cell cycle, metabolism, and antiviral immunity.

microRNAs and thrombo-inflammation: relationship in sight

PURPOSE OF REVIEW

Thrombo-inflammation is a multifaceted pathologic process involving various cells such as platelets, neutrophils, and monocytes. In recent years, microRNAs have been consistently implicated as regulators of these cells.

RECENT FINDINGS

MicroRNAs play a regulatory role in several platelet receptors that have recently been identified as contributing to thrombo-inflammation and neutrophil extracellular trap (NET) formation. In addition, a growing body of evidence has shown that several intracellular and extracellular microRNAs directly promote NET formation.

SUMMARY

Targeting microRNAs is a promising therapeutic approach to control thrombosis in patients with both infectious and noninfectious inflammatory diseases. Future research efforts should focus on elucidating the specific roles of microRNAs in thrombo-inflammation and translating these findings into tangible benefits for patients.

ExplORRNet: An interactive web tool to explore stage-wise miRNA expression profiles and their interactions with mRNA and lncRNA in human breast and gynecological cancers

Background

MicroRNAs (miRNAs) are key regulators of gene expression that have been implicated in gynecological and breast cancers. Understanding the cancer stage-wise expression patterns of miRNAs and their interactions with other RNA molecules in cancer is crucial to improve cancer diagnosis and treatment planning. Comprehensive web tools that integrate data on the transcriptome, circulating miRNAs, and their validated targets to derive beneficial conclusions in cancer research are lacking.

Methods

Using the Shiny R package, we developed a web tool called ExplORRNet that integrates transcriptomic profiles from The Cancer Genome Atlas and miRNA expression data derived from various sources, including tissues, cell lines, exosomes, serum, and plasma, available in the Gene Expression Omnibus database. Differential expression analyses between normal and tumor tissue samples as well as different stages of cancer, accompanied by gene enrichment and survival analyses, can be performed using specialized R packages. Additionally, a miRNA-messenger RNA (mRNA)-long non-coding RNA (lncRNA) networks are constructed to identify regulatory modules.

Results

Our tool identifies cancer stage-wise differentially regulated miRNAs, mRNAs, and lncRNAs in gynecological and breast cancers. Survival analysis identifies miRNAs associated with patient survival, and functional enrichment analysis provides insights into dysregulated miRNA-related biological processes and pathways. The miRNA-mRNA-lncRNA networks highlight interconnected regulatory molecular modules driving cancer progression. Case studies demonstrate the utility of the ExplORRNet for studying gynecological and breast cancers.

Conclusion

ExplORRNet is an intuitive and user-friendly web tool that provides a deeper understanding of dysregulated miRNAs and their functional implications in gynecological and breast cancers. We hope our ExplORRNet tool has potential utility among the clinical and basic researchers and will be beneficial to the entire cancer genomics community to encourage and facilitate mining the rapidly growing public databases to progress the field of precision oncology. The ExplORRNet is available at https://mirna.cs.ut.ee.

Upregulation of developmentally-downregulated miR-1247-5p promotes neuroprotection and axon regeneration in vivo

Numerous micro-RNAs (miRNAs) affect neurodevelopment and neuroprotection, but potential roles of many miRNAs in regulating these processes are still unknown. Here, we used the retinal ganglion cell (RGC) central nervous system (CNS) projection neuron and optic nerve crush (ONC) injury model, to optimize a mature miRNA arm-specific quantification method for characterizing the developmental regulation of miR-1247-5p in RGCs, investigated whether injury affects its expression, and tested whether upregulating miR-1247-5p-mimic in RGCs promotes neuroprotection and axon regeneration. We found that, miR-1247-5p is developmentally-downregulated in RGCs, and is further downregulated after ONC. Importantly, RGC-specific upregulation of miR-1247-5p promoted neuroprotection and axon regeneration after injury in vivo. To gain insight into the underlying mechanisms, we analyzed by bulk-mRNA-seq embryonic and adult RGCs, along with adult RGCs transduced by miR-1247-5p-expressing viral vector, and identified developmentally-regulated cilial and mitochondrial biological processes, which were reinstated to their embryonic levels in adult RGCs by upregulation of miR-1247-5p. Since axon growth is also a developmentally-regulated process, in which mitochondrial dynamics play important roles, it is possible that miR-1247-5p promoted neuroprotection and axon regeneration through regulating mitochondrial functions.

Tissue-specific profiling of age-dependent miRNAomic changes in Caenorhabditis elegans

Ageing exhibits common and distinct features in various tissues, making it critical to decipher the tissue-specific ageing mechanisms. MiRNAs are essential regulators in ageing and are recently highlighted as a class of intercellular messengers. However, little is known about the tissue-specific transcriptomic changes of miRNAs during ageing. C. elegans is a well-established model organism in ageing research. Here, we profile the age-dependent miRNAomic changes in five isolated worm tissues. Besides the diverse ageing-regulated miRNA expression across tissues, we discover numerous miRNAs in the tissues without their transcription. We further profile miRNAs in the extracellular vesicles and find that worm miRNAs undergo inter-tissue trafficking via these vesicles in an age-dependent manner. Using these datasets, we uncover the interaction between body wall muscle-derived mir-1 and DAF-16/FOXO in the intestine, suggesting mir-1 as a messenger in inter-tissue signalling. Taken together, we systematically investigate worm miRNAs in the somatic tissues and extracellular vesicles during ageing, providing a valuable resource to study tissue-autonomous and nonautonomous functions of miRNAs in ageing.

Computational tools supporting known miRNA identification

The study of small RNAs is a field that is expanding quickly. Other functional short RNA molecules other than microRNAs, and gene expression regulators, have been found in animals and plants. MicroRNAs play a significant role in host-microbe interactions, and parasite microRNAs may affect the host's innate immunity. Furthermore, short RNAs are intriguing non-invasive biomarker possibilities because they can be found in physiological fluids. These trends suggest that for many researchers, quick and simple techniques for expression profiling and subsequent downstream analysis of miRNA-seq data are crucial. We selected sRNAtoolbox to make integrated sRNA research easier. Each tool can be used separately or to explore and analyze sRNAbench results in further depth. A special focus was placed on the tools' usability. We review available miRNA research tools to have an overview of the evaluation of the tools. Mainly we evaluate the tool sRNAtoolbox.

Circulating isomiRs May Be Superior Biomarkers Compared to Their Corresponding miRNAs: A Pilot Biomarker Study of Using isomiR-Ome to Detect Coronary Calcium-Based Cardiovascular Risk in Patients with NAFLD

Circulating miRNAs are increasingly being considered as biomarkers in various medical contexts, but the value of analyzing isomiRs (isoforms of canonical miRNA sequences) has not frequently been assessed. Here we hypothesize that an in-depth analysis of the full circulating miRNA landscape could identify specific isomiRs that are stronger biomarkers, compared to their corresponding miRNA, for identifying increased CV risk in patients with non-alcoholic fatty liver disease (NAFLD)-a clinical unmet need. Plasma miRNAs were sequenced with next-generation sequencing (NGS). Liver fat content was measured with magnetic-resonance spectrometry (MRS); CV risk was determined, beyond using traditional biomarkers, by a CT-based measurement of coronary artery calcium (CAC) score and the calculation of a CAC score-based CV-risk percentile (CAC-CV%). This pilot study included n = 13 patients, age > 45 years, with an MRS-measured liver fat content of ≥5% (wt/wt), and free of overt CVD. NGS identified 1103 miRNAs and 404,022 different isomiRs, of which 280 (25%) and 1418 (0.35%), respectively, passed an abundance threshold. Eighteen (sixteen/two) circulating miRNAs correlated positively/negatively, respectively, with CAC-CV%, nine of which also significantly discriminated between high/low CV risk through ROC-AUC analysis. IsomiR-ome analyses uncovered 67 isomiRs highly correlated (R ≥ 0.55) with CAC-CV%. Specific isomiRs of miRNAs 101-3p, 144-3p, 421, and 484 exhibited stronger associations with CAC-CV% compared to their corresponding miRNA. Additionally, while miRNAs 140-3p, 223-3p, 30e-5p, and 342-3p did not correlate with CAC-CV%, specific isomiRs with altered seed sequences exhibited a strong correlation with coronary atherosclerosis burden. Their predicted isomiRs-specific targets were uniquely enriched (compared to their canonical miRNA sequence) in CV Disease (CVD)-related pathways. Two of the isomiRs exhibited discriminative ROC-AUC, and another two showed a correlation with reverse cholesterol transport from cholesterol-loaded macrophages to ApoB-depleted plasma. In summary, we propose a pipeline for exploring circulating isomiR-ome as an approach to uncover novel and strong CVD biomarkers.

JAK1 Is a Novel Target of Tumor- and Invasion-Suppressive microRNA 494-5p in Colorectal Cancer

MiR-494-5p expression has been suggested to be associated with colorectal cancer (CRC) and its metastases in our previous studies. However, functional investigations on the molecule-mediating actions of this miR in CRC are lacking. In silico analysis in the present study revealed a putative binding sequence within the 3'UTR of JAK1. Overexpression of miR-494-5p in cultured CRC significantly reduced the luciferase activity of a reporter plasmid containing the wild-type JAK1-3'UTR, which was abolished by seed sequence mutation. Furthermore, the overexpression of miR-494-5p in CRC cell lines led to a significant reduction in JAK1 expression, proliferation, in vitro migration, and invasion. These effects were abolished by co-transfection with a specific double-stranded RNA that inhibits endogenous miR-494-5p. Moreover, IL-4-induced migration, invasion, and phosphorylation of JAK1, STAT6, and AKT proteins were reduced after an overexpression of this miR, suggesting that this miR affects one of the most essential pathways in CRC. A Kaplan-Meier plotter analysis revealed that patients with high JAK1 expression show reduced survival. Together, these data suggest that miR-494-5p physically inhibits the expression of JAK1 at the translational level as well as in migration and invasion, supporting the hypothesis of miR-494-5p as an early tumor suppressor and inhibitor of early steps of metastasis in CRC.

miR-7 is recruited to the high molecular weight RNA-induced silencing complex in CD8+ T cells upon activation and suppresses IL-2 signaling

Increasing evidence suggests mammalian Argonaute (Ago) proteins partition into distinct complexes within cells, but there is still little biochemical or functional understanding of the miRNAs differentially associated with these complexes. In naïve T cells, Ago2 is found almost exclusively in low molecular weight (LMW) complexes which are associated with miRNAs but not their target mRNAs. Upon T-cell activation, a proportion of these Ago2 complexes move into a newly formed high molecular weight (HMW) RNA-induced silencing complex (RISC), which is characterized by the presence of the GW182 protein that mediates translational repression. Here, we demonstrate distinct partitioning of miRNAs and isomiRs in LMW versus HMW RISCs upon antigen-mediated activation of CD8+ T cells. We identify miR-7 as highly enriched in HMW RISC and demonstrate that miR-7 inhibition leads to increased production of IL-2 and up-regulation of the IL-2 receptor, the transferrin receptor, CD71 and the amino acid transporter, CD98. Our data support a model where recruitment of miR-7 to HMW RISC restrains IL-2 signaling and the metabolic processes regulated by IL-2.

Interrogating the Role of miR-125b and Its 3'isomiRs in Protection against Hypoxia

MiR-125b has therapeutic potential in the amelioration of myocardial ischemic injury. MicroRNA isomiRs, with either 5' or 3' addition or deletion of nucleotide(s), have been reported from next-generation sequencing data (NGS). However, due to technical challenges, validation and functional studies of isomiRs are few. In this study, we discovered using NGS, four 3'isomiRs of miR-125b, i.e., addition of A (adenosine), along with deletions of A, AG (guanosine) and AGU (uridine) from rat and sheep heart. These findings were validated using RT-qPCR. Comprehensive functional studies were carried out in the H9C2 hypoxia model. After miR-125b, isomiRs of Plus A, Trim A, AG and AGU mimic transfection, the H9C2 cells were subjected to hypoxic challenge. As assessed using cell viability, apoptosis, CCK-8 and LDH release, miR-125b and isomiRs were all protective against hypoxia. However, Plus A and Trim A were more effective than miR-125b, whilst Trim AG and Trim AGU had far weaker effects than miR-125b. Interestingly, both the gene regulation profile and apoptotic gene validation indicated a major overlap among miR-125b, Plus A and Trim A, whilst Trims AG and AGU revealed a different profile compared to miR-125b. Conclusions: miR-125b and its 3' isomiRs are expressed stably in the heart. miR-125b and isomiRs with addition or deletion of A might function concurrently and concordantly under specific physiological and pathophysiological conditions. In-depth understanding of isomiRs' metabolism and function will contribute to better miRNA therapeutic drug design.

The Typical tRNA Co-Expresses Multiple 5' tRNA Halves Whose Sequences and Abundances Depend on Isodecoder and Isoacceptor and Change with Tissue Type, Cell Type, and Disease

Transfer RNA-derived fragments (tRFs) are noncoding RNAs that arise from either mature transfer RNAs (tRNAs) or their precursors. One important category of tRFs comprises the tRNA halves, which are generated through cleavage at the anticodon. A given tRNA typically gives rise to several co-expressed 5'-tRNA halves (5'-tRHs) that differ in the location of their 3' ends. These 5'-tRHs, even though distinct, have traditionally been treated as indistinguishable from one another due to their near-identical sequences and lengths. We focused on co-expressed 5'-tRHs that arise from the same tRNA and systematically examined their exact sequences and abundances across 10 different human tissues. To this end, we manually curated and analyzed several hundred human RNA-seq datasets from NCBI's Sequence Run Archive (SRA). We grouped datasets from the same tissue into their own collection and examined each group separately. We found that a given tRNA produces different groups of co-expressed 5'-tRHs in different tissues, different cell lines, and different diseases. Importantly, the co-expressed 5'-tRHs differ in their sequences, absolute abundances, and relative abundances, even among tRNAs with near-identical sequences from the same isodecoder or isoacceptor group. The findings suggest that co-expressed 5'-tRHs that are produced from the same tRNA or closely related tRNAs have distinct, context-dependent roles. Moreover, our analyses show that cell lines modeling the same tissue type and disease may not be interchangeable when it comes to experimenting with tRFs.

The Role of MicroRNAs in Liver Functioning: from Biogenesis to Therapeutic Approaches (Review)

Molecular diagnostics based on small non-coding RNA molecules (in particular microRNA) is a new direction in modern biomedicine and is considered a promising method for identification of a wide range of pathologies at an early stage, clinical phenotype assessment, as well as monitoring the course of the disease, evaluation of therapy efficacy and the risk of the disease recurrence. Currently, the role of microRNAs as the most important epigenetic regulator in cancer development has been proven within the studies of normal and pathogenic processes. However, currently, there are insignificant studies devoted to studying the role of microRNAs in functioning of other organs and tissues, as well as to development of possible therapeutic approaches based on microRNAs. A huge number of metabolic processes in the liver are controlled by microRNAs, which creates enormous potential for the use of microRNAs as a diagnostic marker and makes it a target for therapeutic intervention in metabolic, oncological, and even viral diseases of this organ. This review examines various aspects of biological functions of microRNAs in different types of liver cells. Both canonical and non-canonical pathways of biogenesis, epigenetic regulation mediated by microRNAs, as well as the microRNAs role in intercellular communication and the course of viral diseases are shown. The potential of microRNAs as a diagnostic marker for various liver pathologies is described, as well as therapeutic approaches and medicines based on microRNAs, which are approved for clinical use and currently being developed.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

Genetic Regulation of Human isomiR Biogenesis

MicroRNAs play a critical role in regulating gene expression post-transcriptionally. Variations in mature microRNA sequences, known as isomiRs, arise from imprecise cleavage and nucleotide substitution or addition. These isomiRs can target different mRNAs or compete with their canonical counterparts, thereby expanding the scope of miRNA post-transcriptional regulation. Our study investigated the relationship between cis-acting single-nucleotide polymorphisms (SNPs) in precursor miRNA regions and isomiR composition, represented by the ratio of a specific 5'-isomiR subtype to all isomiRs identified for a particular mature miRNA. Significant associations between 95 SNP-isomiR pairs were identified. Of note, rs6505162 was significantly associated with both the 5'-extension of hsa-miR-423-3p and the 5'-trimming of hsa-miR-423-5p. Comparison of breast cancer and normal samples revealed that the expression of both isomiRs was significantly higher in tumors than in normal tissues. This study sheds light on the genetic regulation of isomiR maturation and advances our understanding of post-transcriptional regulation by microRNAs.

RNA-RNA competitive interactions: a molecular civil war ruling cell physiology and diseases

The idea that proteins are the main determining factors in the functioning of cells and organisms, and their dysfunctions are the first cause of pathologies, has been predominant in biology and biomedicine until recently. This protein-centered view was too simplistic and failed to explain the physiological and pathological complexity of the cell. About 80% of the human genome is dynamically and pervasively transcribed, mostly as non-protein-coding RNAs (ncRNAs), which competitively interact with each other and with coding RNAs generating a complex RNA network regulating RNA processing, stability, and translation and, accordingly, fine-tuning the gene expression of the cells. Qualitative and quantitative dysregulations of RNA-RNA interaction networks are strongly involved in the onset and progression of many pathologies, including cancers and degenerative diseases. This review will summarize the RNA species involved in the competitive endogenous RNA network, their mechanisms of action, and involvement in pathological phenotypes. Moreover, it will give an overview of the most advanced experimental and computational methods to dissect and rebuild RNA networks.

Integrative genetic and genomic networks identify microRNA associated with COPD and ILD

Chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD) are clinically and molecularly heterogeneous diseases. We utilized clustering and integrative network analyses to elucidate roles for microRNAs (miRNAs) and miRNA isoforms (isomiRs) in COPD and ILD pathogenesis. Short RNA sequencing was performed on 351 lung tissue samples of COPD (n = 145), ILD (n = 144) and controls (n = 64). Five distinct subclusters of samples were identified including 1 COPD-predominant cluster and 2 ILD-predominant clusters which associated with different clinical measurements of disease severity. Utilizing 262 samples with gene expression and SNP microarrays, we built disease-specific genetic and expression networks to predict key miRNA regulators of gene expression. Members of miR-449/34 family, known to promote airway differentiation by repressing the Notch pathway, were among the top connected miRNAs in both COPD and ILD networks. Genes associated with miR-449/34 members in the disease networks were enriched among genes that increase in expression with airway differentiation at an air-liquid interface. A highly expressed isomiR containing a novel seed sequence was identified at the miR-34c-5p locus. 47% of the anticorrelated predicted targets for this isomiR were distinct from the canonical seed sequence for miR-34c-5p. Overexpression of the canonical miR-34c-5p and the miR-34c-5p isomiR with an alternative seed sequence down-regulated NOTCH1 and NOTCH4. However, only overexpression of the isomiR down-regulated genes involved in Ras signaling such as CRKL and GRB2. Overall, these findings elucidate molecular heterogeneity inherent across COPD and ILD patients and further suggest roles for miR-34c in regulating disease-associated gene-expression.

Phospho-RNA-Seq Highlights Specific Small RNA Profiles in Plasma Extracellular Vesicles

Small RNAs (sRNAs) are bioactive molecules that can be detected in biofluids, reflecting physiological and pathological states. In plasma, sRNAs are found within extracellular vesicles (EVs) and in extravesicular compartments, offering potential sources of highly sensitive biomarkers. Deep sequencing strategies to profile sRNAs favor the detection of microRNAs (miRNAs), the best-known class of sRNAs. Phospho-RNA-seq, through the enzymatic treatment of sRNAs with T4 polynucleotide kinase (T4-PNK), has been recently developed to increase the detection of thousands of previously inaccessible RNAs. In this study, we investigated the value of phospho-RNA-seq on both the EVs and extravesicular plasma subfractions. Phospho-RNA-seq increased the proportion of sRNAs used for alignment and highlighted the diversity of the sRNA transcriptome. Unsupervised clustering analysis using sRNA counts matrices correctly classified the EVs and extravesicular samples only in the T4-PNK treated samples, indicating that phospho-RNA-seq stresses the features of sRNAs in each plasma subfraction. Furthermore, T4-PNK treatment emphasized specific miRNA variants differing in the 5'-end (5'-isomiRs) and certain types of tRNA fragments in each plasma fraction. Phospho-RNA-seq increased the number of tissue-specific messenger RNA (mRNA) fragments in the EVs compared with the extravesicular fraction, suggesting that phospho-RNA-seq favors the discovery of tissue-specific sRNAs in EVs. Overall, the present data emphasizes the value of phospho-RNA-seq in uncovering RNA-based biomarkers in EVs.

Differentially Expressed Conserved Plant Vegetative Phase-Change-Related microRNAs in Mature and Rejuvenated Silver Birch In Vitro Propagated Tissues

In plants, phase change from the juvenile stage to maturity involves physiological and anatomical changes, which are initiated and controlled by evolutionary highly conserved microRNAs. This process is of particular significance for the in vitro propagation of woody plant species, as individuals or tissues that have undergone the transition to vegetative maturity are recalcitrant to propagation. Conserved miRNAs differentially expressed between juvenile (including rejuvenated) and mature silver birch tissues were identified using high-throughput sequencing of small RNA libraries. Expression of some miR156 isoforms was high in juvenile tissues and has been previously reported to regulate phase transitions in a range of species. Additional miRNAs, such as miR394 and miR396, that were previously reported to be highly expressed in juvenile woody plant tissues were also differentially expressed in this study. However, expression of miR172, previously reported to be highly expressed in mature tissues, was low in all sample types in this study. The obtained results will provide insight for further investigation of the molecular mechanisms regulating vegetative phase change in silver birch and other perennial woody plant species, by analysing a wider range of genotypes, tissue types and maturation stages. This knowledge can potentially assist in identification of rejuvenated material at an earlier stage than currently possible, increasing the efficiency of silver birch in vitro propagation.

Circulating Small Noncoding RNA Profiling as a Potential Biomarker of Atherosclerotic Plaque Composition in Type 1 Diabetes

OBJECTIVE

Cardiovascular disease (CVD) accounts for most deaths in patients with type 1 diabetes (T1D); however, the determinants of plaque composition are unknown. miRNAs regulate gene expression, participate in the development of atherosclerosis, and represent promising CVD biomarkers. This study analyzed the circulating miRNA expression profile in T1D with either carotid calcified (CCP) or fibrous plaque (CFP).

RESEARCH DESIGN AND METHODS

Circulating small noncoding RNAs were sequenced and quantified using next-generation sequencing and bioinformatic analysis in an exploratory set of 26 subjects with T1D with CCP and in 25 with CFP. Then, in a validation set of 40 subjects with CCP, 40 with CFP, and 24 control subjects with T1D, selected miRNA expression was measured by digital droplet PCR. Putative gene targets enriched for pathways implicated in atherosclerosis/vascular calcification/diabetes were analyzed. The patients' main clinical characteristics were also recorded.

RESULTS

miR-503-5p, let-7d-5p, miR-106b-3p, and miR-93-5p were significantly upregulated, while miR-10a-5p was downregulated in patients with CCP compared with CFP (all fold change >±1.5; P < 0.05). All candidate miRNAs showed a significant correlation with LDL-cholesterol, direct for the upregulated and inverse for the downregulated miRNA, in CCP. Many target genes of upregulated miRNAs in CCP participate in osteogenic differentiation, apoptosis, inflammation, cholesterol metabolism, and extracellular matrix organization.

CONCLUSIONS

These findings characterize miRNAs and their signature in the regulatory network of carotid plaque phenotype in T1D, providing new insights into plaque pathophysiology and possibly novel biomarkers of plaque composition.

Diversity and Differential Expression of MicroRNAs in the Human Skeletal Muscle with Distinct Fiber Type Composition

The ratio of fast- and slow-twitch fibers in human skeletal muscle is variable and largely determined by genetic factors. In this study, we investigated the contribution of microRNA (miRNA) in skeletal muscle fiber type composition. The study involved biopsy samples of the vastus lateralis muscle from 24 male participants with distinct fiber type ratios. The miRNA study included samples from five endurance athletes and five power athletes with the predominance of slow-twitch (61.6-72.8%) and fast-twitch (69.3-80.7%) fibers, respectively. Total and small RNA were extracted from tissue samples. Total RNA sequencing (N = 24) revealed 352 differentially expressed genes between the groups with the predominance of fast- and slow-twitch muscle fibers. Small RNA sequencing showed upregulation of miR-206, miR-501-3p and miR-185-5p, and downregulation of miR-499a-5p and miR-208-5p in the group of power athletes with fast-twitch fiber predominance. Two miRtronic miRNAs, miR-208b-3p and miR-499a-5p, had strong correlations in expression with their host genes (MYH7 and MYH7B, respectively). Correlations between the expression of miRNAs and their experimentally validated messenger RNA (mRNA) targets were calculated, and 11 miRNA-mRNA interactions with strong negative correlations were identified. Two of them belonged to miR-208b-3p and miR-499a-5p, indicating their regulatory links with the expression of CDKN1A and FOXO4, respectively.

One locus, several functional RNAs-emerging roles of the mechanisms responsible for the sequence variability of microRNAs

With the development of modern molecular genetics, the original "one gene-one enzyme" hypothesis has been outdated. For protein coding genes, the discovery of alternative splicing and RNA editing provided the biochemical background for the RNA repertoire of a single locus, which also serves as an important pillar for the enormous protein variability of the genomes. Non-protein coding RNA genes were also revealed to produce several RNA species with distinct functions. The loci of microRNAs (miRNAs), encoding for small endogenous regulatory RNAs, were also found to produce a population of small RNAs, rather than a single defined product. This review aims to present the mechanisms contributing to the astonishing variability of miRNAs revealed by the new sequencing technologies. One important source is the careful balance of arm selection, producing sequentially different 5p- or 3p-miRNAs from the same pre-miRNA, thereby broadening the number of regulated target RNAs and the phenotypic response. In addition, the formation of 5', 3' and polymorphic isomiRs, with variable end and internal sequences also leads to a higher number of targeted sequences, and increases the regulatory output. These miRNA maturation processes, together with other known mechanisms such as RNA editing, further increase the potential outcome of this small RNA pathway. By discussing the subtle mechanisms behind the sequence diversity of miRNAs, this review intends to reveal this engaging aspect of the inherited "RNA world", how it contributes to the almost infinite molecular variability among living organisms, and how this variability can be exploited to treat human diseases.

Differential co-expression network analysis with DCoNA reveals isomiR targeting aberrations in prostate cancer

MOTIVATION

One of the standard methods of high-throughput RNA sequencing analysis is differential expression. However, it does not detect changes in molecular regulation. In contrast to the standard differential expression analysis, differential co-expression one aims to detect pairs or clusters whose mutual expression changes between two conditions.

RESULTS

We developed Differential Co-expression Network Analysis (DCoNA)-an open-source statistical tool that allows one to identify pair interactions, which correlation significantly changes between two conditions. Comparing DCoNA with the state-of-the-art analog, we showed that DCoNA is a faster, more accurate and less memory-consuming tool. We applied DCoNA to prostate mRNA/miRNA-seq data collected from The Cancer Genome Atlas (TCGA) and compared predicted regulatory interactions of miRNA isoforms (isomiRs) and their target mRNAs between normal and cancer samples. As a result, almost all highly expressed isomiRs lost negative correlation with their targets in prostate cancer samples compared to ones without the pathology. One exception to this trend was the canonical isomiR of hsa-miR-93-5p acquiring cancer-specific targets. Further analysis showed that cancer aggressiveness simultaneously increased with the expression level of this isomiR in both TCGA primary tumor samples and 153 blood plasma samples of P. Hertsen Moscow Oncology Research Institute patients' cohort analyzed by miRNA microarrays.

AVAILABILITY AND IMPLEMENTATION

Source code and documentation of DCoNA are available at https://github.com/zhiyanov/DCoNA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The dual role of microRNA (miR)-20b in cancers: Friend or foe?

MicroRNAs, as non-coding transcripts, modulate gene expression through RNA silencing under normal physiological conditions. Their aberrant expression has strongly associated with tumorigenesis and cancer development. MiR-20b is one of the crucial miRNAs that regulate essential biological processes such as cell proliferation, apoptosis, autophagy, and migration. Deregulated levels of miR-20b contribute to the early- and advanced stages of cancer. On the other hand, investigations emphasize the tumor suppressor ability of miR-20b. High-throughput strategies are developed to identify miR-20b potential targets, providing the proper insight into its molecular mechanism of action. Moreover, accumulated results suggest that miR-20b exerts its effects through diverse signaling pathways, including PI3K/AKT/mTOR and ERK axes. Restoration of the altered expression levels of miR-20b induces cell apoptosis and reduces invasion and migration. Further, miR-20b can be used as a biomarker in cancer. The current comprehensive review could lead to a better understanding of the miR-20b in either tumorigenesis or tumor regression that may open new avenues for cancer treatment. Video Abstract.

Genetics and RNA Regulation of Uveal Melanoma

Uveal melanoma (UM) is the most common intraocular malignant tumor and the most frequent melanoma not affecting the skin. While the rate of UM occurrence is relatively low, about 50% of patients develop metastasis, primarily to the liver, with lethal outcome despite medical treatment. Notwithstanding that UM etiopathogenesis is still under investigation, a set of known mutations and chromosomal aberrations are associated with its pathogenesis and have a relevant prognostic value. The most frequently mutated genes are BAP1, EIF1AX, GNA11, GNAQ, and SF3B1, with mutually exclusive mutations occurring in GNAQ and GNA11, and almost mutually exclusive ones in BAP1 and SF3B1, and BAP1 and EIF1AX. Among chromosomal aberrations, monosomy of chromosome 3 is the most frequent, followed by gain of chromosome 8q, and full or partial loss of chromosomes 1 and 6. In addition, epigenetic mechanisms regulated by non-coding RNAs (ncRNA), namely microRNAs and long non-coding RNAs, have also been investigated. Several papers investigating the role of ncRNAs in UM have reported that their dysregulated expression affects cancer-related processes in both in vitro and in vivo models. This review will summarize current findings about genetic mutations, chromosomal aberrations, and ncRNA dysregulation establishing UM biology.

Alarmins and MicroRNAs, a New Axis in the Genesis of Respiratory Diseases: Possible Therapeutic Implications

It is well ascertained that airway inflammation has a key role in the genesis of numerous respiratory pathologies, including asthma, chronic obstructive pulmonary disease, and acute respiratory distress syndrome. Pulmonary tissue inflammation and anti-inflammatory responses implicate an intricate relationship between local and infiltrating immune cells and structural pulmonary cells. Alarmins are endogenic proteins discharged after cell injury in the extracellular microenvironment. The purpose of our review is to highlight the alterations in respiratory diseases involving some alarmins, such as high mobility group box 1 (HMGB1) and interleukin (IL)-33, and their inter-relationships and relationships with genetic non-coding material, such as microRNAs. The role played by these alarmins in some pathophysiological processes confirms the existence of an axis composed of HMGB1 and IL-33. These alarmins have been implicated in ferroptosis, the onset of type 2 inflammation and airway alterations. Moreover, both factors can act on non-coding genetic material capable of modifying respiratory function. Finally, we present an outline of alarmins and RNA-based therapeutics that have been proposed to treat respiratory pathologies.

Identification of extremely GC-rich micro RNAs for RT-qPCR data normalization in human plasma

We aimed at extending the repertoire of high-quality miRNA normalizers for reverse transcription-quantitative PCR (RT-qPCR) of human plasma with special emphasis on the extremely guanine-cytosine-rich portion of the miRNome. For high-throughput selection of stable candidates, microarray technology was preferred over small-RNA sequencing (sRNA-seq) since the latter underrepresented miRNAs with a guanine-cytosine (GC) content of at least 75% (p = 0.0002, n = 2). miRNA abundances measured on the microarray were ranked for consistency and uniformity using nine normalization approaches. The eleven most stable sequences included miRNAs of moderate, but also extreme GC content (45%-65%: miR-320d, miR-425-5p, miR-185-5p, miR-486-5p; 80%-95%: miR-1915-3p, miR-3656-5p, miR-3665-5p, miR-3960-5p, miR-4488-5p, miR-4497 and miR-4787-5p). In contrast, the seven extremely GC-rich miRNAs were not found in the two plasma miRNomes screened by sRNA-seq. Stem-loop RT-qPCR was employed for stability verification in 32 plasma samples of healthy male Caucasians (age range: 18-55 years). In general, inter-individual variance of miRNA abundance was low or very low as indicated by coefficient of variation (CV) values of 0.6%-8.2%. miR-3665 and miR-1915-3p outperformed in this analysis (CVs: 0.6 and 2.4%, respectively). The eight most stable sequences included four extremely GC-rich miRNAs (miR-1915-3p, miR-3665, miR-4787-5p and miR-4497). The best-performing duo normalization factor (NF) for the condition of human plasma, miR-320d and miR-4787-5p, also included a GC-extreme miRNA. In summary, the identification of extremely guanine-cytosine-rich plasma normalizers will help to increase accuracy of PCR-based miRNA quantification, thus raise the potential that miRNAs become markers for psychological stress reactions or early and precise diagnosis of clinical phenotypes. The novel miRNAs might also be useful for orthologous contexts considering their conservation in related animal genomes.

Induction of Somatic Embryogenesis in Plants: Different Players and Focus on WUSCHEL and WUS-RELATED HOMEOBOX (WOX) Transcription Factors

In plants, other cells can express totipotency in addition to the zygote, thus resulting in embryo differentiation; this appears evident in apomictic and epiphyllous plants. According to Haberlandt's theory, all plant cells can regenerate a complete plant if the nucleus and the membrane system are intact. In fact, under in vitro conditions, ectopic embryos and adventitious shoots can develop from many organs of the mature plant body. We are beginning to understand how determination processes are regulated and how cell specialization occurs. However, we still need to unravel the mechanisms whereby a cell interprets its position, decides its fate, and communicates it to others. The induction of somatic embryogenesis might be based on a plant growth regulator signal (auxin) to determine an appropriate cellular environment and other factors, including stress and ectopic expression of embryo or meristem identity transcription factors (TFs). Still, we are far from having a complete view of the regulatory genes, their target genes, and their action hierarchy. As in animals, epigenetic reprogramming also plays an essential role in re-establishing the competence of differentiated cells to undergo somatic embryogenesis. Herein, we describe the functions of WUSCHEL-RELATED HOMEOBOX (WOX) transcription factors in regulating the differentiation-dedifferentiation cell process and in the developmental phase of in vitro regenerated adventitious structures.

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.