Toward the human cellular microRNAome

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.

Bioinformatic Analysis of Autism-Related miRNAs and Their PoTential as Biomarkers for Autism Epigenetic Inheritance

BACKGROUND/OBJECTIVES

The dysregulation of miRNA expression in samples from autistic individuals indicates that they are involved in autism. The participation of miRNAs in paternal epigenetic inheritance has also been reported. This study used bioinformatics tools to analyze the literature and genetic databases to search for miRNAs associated with autism, aiming to explore their suitability to investigate paternal epigenetic inheritance.

METHODS

Autism-related miRNAs were searched in public databases using bioinformatic tools (miRNA-to-genes analysis). The genes targeted by these autism-related miRNAs, which are common to neurons, sperm, and PBMCs, were identified. Enrichment analyses were performed to identify the biological processes regulated by the candidate miRNAs. Autism-related miRNAs were also identified by an inverse analysis (genes-to-miRNA analysis), starting from autism-related genes.

RESULTS

In the miRNA-to-gene analysis, 416 miRNAs involved in autism were found, of which 77 were expressed in sperm, PBMCs, and neurons. From these, 18 were differentially expressed in the brain and in at least one peripheral sample (saliva or blood), suggesting that they might be suitable to be used in the investigation of autism biomarkers and inheritance. In the genes-to-miRNA analysis, 36 miRNAs were identified, from which 9 coincided with the results of direct analysis.

CONCLUSIONS

Although there is no consensus about miRNAs related to autism, there are candidate miRNAs that show clear potential to be explored as biomarkers. The coincidence in the expression of miRNAs in sperm, neurons, and PBMCs indicates that they are valuable biological samples to study the role of miRNAs in the paternal epigenetic inheritance of autism.

Investigation of the role of miRNA variants in neurodegenerative brain diseases

Introduction

miRNAs are small noncoding elements known to regulate different molecular processes, including developmental and executive functions in the brain. Dysregulation of miRNAs could contribute to brain neurodegeneration, as suggested by miRNA profiling studies of individuals suffering from neurodegenerative brain diseases (NBDs). Here, we report rare miRNA variants in patients with Alzheimer's dementia (AD) and frontotemporal dementia (FTD).

Methods

We initially used whole exome sequencing data in a subset of FTD patients (n = 209) from Flanders-Belgium. We then performed targeted resequencing of variant-harboring miRNAs in an additional subset of FTD patients (n = 126) and control individuals (n = 426). Lastly, we sequenced the MIR885 locus in a Flanders-Belgian AD cohort (n = 947) and a total number of n = 755 controls.

Results

WES identified rare seed variants in MIR656, MIR423, MIR122 and MIR885 in FTD patients. Most of these miRNAs bind to FTD-associated genes, implicated in different biological pathways. Additionally, some miRNA variants create novel binding sites for genes associated with FTD. Sequencing of the MIR885 locus in the AD cohort initially showed a significant enrichment of MIR885 variants in AD patients compared to controls (SKAT-O, p-value = 0.026). Genetic association was not maintained when we included sex and APOE status as covariates. Using the miRVaS prediction tool, variants rs897551430 and rs993255773 appeared to evoke significant structural changes in the primary miRNA. These variants are also predicted to strongly downregulate mature miR885 levels, in line with what is reported for MIR885 in the context of AD.

Discussion

Functional investigation of miRNAs/variants described in this study could propose novel miRNA-mediated molecular cascades in FTD and AD pathogenicity. Furthermore, we believe that the genetic evidence presented here suggests a role for MIR885 in molecular mechanisms involved in AD and warrants genetic follow-up in larger cohorts to explore this hypothesis.

Multi-Omics Study Reveals Nc886/vtRNA2-1 as a Positive Regulator of Prostate Cancer Cell Immunity

Noncoding RNA 886 has emerged as a pivotal regulatory RNA with distinct functions across tissues, acting as a regulator of protein activity by directly binding to protein partners. While it is well recognized as a tumor suppressor in prostate cancer, the underlying molecular mechanisms remain elusive. To discover the principal pathways regulated by nc886 in prostate cancer, we used a transcriptomic and proteomic approach analyzing malignant DU145, LNCaP, PC3, and benign RWPE-1 prostate cell line models transiently transfected with in vitro transcribed nc886 or antisense oligonucleotides. Multiomics revelead a significant enrichment of immune system-related pathways across the cell lines, including cytokines and interferon signaling. The interferon response provoked by nc886 was validated by functional assays. The invariability of PKR phosphorylation and NF-κB activity in the gain/loss of nc886 function experiments and the positive regulation of innate immunity suggest a PKR-independent mechanism of nc886 action. Accordingly, the GSEA of the PRAD-TCGA data set revealed immune stimulation as the nc886 most associated node also in the clinical setting. Our study showed that the reduction of nc886 leads to a blunted immune response in prostate cancer.

Beyond the Bubble: A Debate on microRNA Sorting Into Extracellular Vesicles

Over the past decade, a scientific field has been developed demonstrating microRNAs (miRNAs) to be actively sorted into extracellular vesicles via specific nucleotide motifs that interact with discrete RNA-binding proteins. These miRNAs are proposed to be transported into recipient cells in which they can regulate specific cellular pathways. This mechanism could have enormous potential in explaining how cells signal and regulate other cells nearby or at a distance. Tens of studies have built this theme of a regulated transport of miRNAs. However, some concerns exist about this field. Taken together, there are concerns of a lack of a consistent motif, RNA-binding protein, or preferential miRNA involved in this process. In this study, we provide an expert and extensive analysis of the field that makes the cases for and against an active sorting mechanism. We provide potential explanations on why there is a lack of agreement. Most importantly, we provide ideas on how to move this field forward with more rigor and reproducibility. It is hoped that by engaging in a scientific debate of the pros and cons of this field, more rigorous experiments can be performed to conclusively demonstrate this biological activity.

Mapping the microRNA landscape in the older adult brain and its genetic contribution to neuropsychiatric conditions

MicroRNAs (miRNAs) play a crucial role in regulating gene expression and influence many biological processes. Despite their importance, understanding of how genetic variation affects miRNA expression in the brain and how this relates to brain disorders remains limited. Here we investigated these questions by identifying microRNA expression quantitative trait loci (miR-QTLs), or genetic variants associated with brain miRNA levels, using genome-wide small RNA sequencing profiles from dorsolateral prefrontal cortex samples of 604 older adult donors of European ancestry. Here we show that nearly half (224 of 470) of the analyzed miRNAs have associated miR-QTLs, many of which fall in regulatory regions such as brain promoters and enhancers. We also demonstrate that intragenic miRNAs often have genetic regulation independent from their host genes. Furthermore, by integrating our findings with 16 genome-wide association studies of psychiatric and neurodegenerative disorders, we identified miRNAs that likely contribute to bipolar disorder, depression, schizophrenia and Parkinson's disease. These findings advance understanding of the genetic regulation of miRNAs and their role in brain health and disease.

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.

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.

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.

MiR-30a-5p isoform -1|1 promotes the progression of gastric cancer by inhibiting TMEM66 and reducing intratumoral cytotoxic T cells

Gastric cancer is histologically classified into the intestinal subtype, which forms tubular structures, and the aggressive diffuse subtype, characterized by rapid invasion and poor prognosis. The variety and quantity of miRNA isoforms between different histological subtypes of gastric cancer were unknown. Through systematic filtering, we found that more diverse miR-30a-5p isoforms was present in the diffuse subtype of gastric cancer, and was associated with patients' worse survival independent of tumor stage based on the TCGA miRNA-seq data. Among all nine isoforms of miR-30a-5p, miR-30a-5p -1|1 was more abundant than the archetype of miR-30a-5p. Higher expression of miR-30a-5p -1|1 was observed in patients with advanced tumor stage and poor survival. Furthermore, miR-30a-5p -1|1 could promote the metastasis of gastric cancer cells both in vitro and in vivo by down-regulating TMEM66. In clinical samples, decreased expression of TMEM66 was characteristic of gastric cancer, and the low level of TMEM66 correlated with deceased CD8 positive cells in the tumor microenvironment probably due to decreased cytokines production. In conclusion, the variety of miR-30a-5p isoforms correlates with worse survival in gastric cancer patients. Moreover, miR-30a-5p -1|1 could promote gastric cancer metastasis by inhibiting TMEM66 and the infiltration of intratumoral CD8 positive cells.

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.

Role of CircCHD2 in the pathogenesis of gestational diabetes mellitus by regulating autophagy via miR-33b-3p/ULK1 axis

Gestational diabetes mellitus (GDM) is a common pregnancy complication with a high incidence in women; however, its pathophysiology remains unknown. Our previous study suggested that the circCHD2/miR-33b-3p/ULK1 axis may be involved in GDM pathogenesis. However, the mechanism through which circCHD2 regulates GDM development requires further investigation. We found that high-glucose (HG, 25 mmol/L) significantly induced the expression of circCHD2, increased autophagy and apoptosis, and decreased cell viability in human placental trophoblast HTR-8/SVneo cells. In contrast, the downregulation of circCHD2 significantly attenuated the effects of HG on HTR-8/SVneo cells. MiR-33b-3p downregulated in the placenta of GDM patients was reduced by HG and detected as a target of circCHD2 using bioinformatics analysis, a dual-luciferase reporter assay, and qRT-PCR assay. Further studies showed that the inhibition of miR-33b-3p significantly blocked the effects of circCHD2 downregulation on cell viability, apoptosis, and autophagy in HG-treated HTR-8/SVneo cells. ULK1 is a target of miR-33b-3p, based on bioinformatics analysis, a dual-luciferase reporter assay, qRT-PCR assay, and Western blot analysis. Compared to miR-33b-3p, ULK1 is upregulated in the placenta of GDM patients. ULK1 overexpression notably blocked the effects of miR-33b-3p mimics on cell viability, apoptosis, and autophagy in HG-treated HTR-8/SVneo cells. These findings suggested that circCHD2 acts as an autophagy promoter via the miR-33b-3p/ULK1 axis to induce apoptosis in HTR-8/SVneo cells, suggesting that circCHD2 is a potential diagnostic and therapeutic target for GDM.

Visceral adipocyte size links obesity with dysmetabolism more than fibrosis, and both can be estimated by circulating miRNAs

OBJECTIVE

In obesity, adipocyte hypertrophy is detrimental to health, but its' interrelation with fibrosis in the visceral adipose tissue (VAT) depot remains unclear. Because VAT is less accessible via biopsy, biomarkers for VAT quality are needed. The authors hypothesized that VAT adipocyte size and fibrosis are interrelated and can be estimated by circulating microRNAs (circ-miRNAs), contributing to subphenotyping obesity.

METHODS

Adipocyte size and AT fibrosis were estimated in n = 43 participants (BMI ≥ 30 kg/m2 ). Circ-miRNAs were sequenced (Next Generation Sequencing).

RESULTS

Participants with above- versus below-median VAT adipocyte area exhibited metabolic dysfunction but lower total and pericellular fibrosis. VAT adipocyte size remained associated with metabolic dysfunction even when controlling for BMI or VAT fibrosis in the entire cohort, as in matched-pairs subanalyses. Next Generation Sequencing uncovered 22 and 6 circ-miRNAs associated with VAT adipocyte size and fibrosis, respectively, with miRNA-130b-3p common to both analyses. The combination of miRNA-130b-3p + miR-150-5p + high-density lipoprotein cholesterol discriminated among those with large versus small VAT adipocytes (receiver operating characteristic-area under the curve: 0.872 [95% CI: 0.747-0.996]), whereas miRNA-130b-3p + miRNA-15a-5p + high-density lipoprotein cholesterol discriminated among those with low and high fibrosis (receiver operating characteristic-area under the curve: 0.823 [95% CI: 0.676-0.97]).

CONCLUSIONS

These findings suggest that VAT adipocyte size and fibrosis are inversely correlated in obesity and can be estimated by distinct circ-miRNAs, providing a potential tool to subphenotype obesity via a liquid biopsy-like approach to assess VAT health in nonsurgical patients.

Associations Between Epigenetic Age Acceleration and microRNA Expression Among U.S. Firefighters

Epigenetic changes may be biomarkers of health. Epigenetic age acceleration (EAA), the discrepancy between epigenetic age measured via epigenetic clocks and chronological age, is associated with morbidity and mortality. However, the intersection of epigenetic clocks with microRNAs (miRNAs) and corresponding miRNA-based health implications have not been evaluated. We analyzed DNA methylation and miRNA profiles from blood sampled among 332 individuals enrolled across 2 U.S.-based firefighter occupational studies (2015-2018 and 2018-2020). We considered 7 measures of EAA in leukocytes (PhenoAge, GrimAge, Horvath, skin-blood, and Hannum epigenetic clocks, and extrinsic and intrinsic epigenetic age acceleration). We identified miRNAs associated with EAA using individual linear regression models, adjusted for sex, race/ethnicity, chronological age, and cell type estimates, and investigated downstream effects of associated miRNAs with miRNA enrichment analyses and genomic annotations. On average, participants were 38 years old, 88% male, and 75% non-Hispanic white. We identified 183 of 798 miRNAs associated with EAA (FDR q < 0.05); 126 with PhenoAge, 59 with GrimAge, 1 with Horvath, and 1 with the skin-blood clock. Among miRNAs associated with Horvath and GrimAge, there were 61 significantly enriched disease annotations including age-related metabolic and cardiovascular conditions and several cancers. Enriched pathways included those related to proteins and protein modification. We identified miRNAs associated with EAA of multiple epigenetic clocks. PhenoAge had more associations with individual miRNAs, but GrimAge and Horvath had greater implications for miRNA-associated pathways. Understanding the relationship between these epigenetic markers could contribute to our understanding of the molecular underpinnings of aging and aging-related diseases.

Targeting a cell-specific microRNA repressor of CXCR4 ameliorates atherosclerosis in mice

The CXC chemokine receptor 4 (CXCR4) in endothelial cells (ECs) and vascular smooth muscle cells (VSMCs) is crucial for vascular integrity. The atheroprotective functions of CXCR4 in vascular cells may be counteracted by atherogenic functions in other nonvascular cell types. Thus, strategies for cell-specifically augmenting CXCR4 function in vascular cells are crucial if this receptor is to be useful as a therapeutic target in treating atherosclerosis and other vascular disorders. Here, we identified miR-206-3p as a vascular-specific CXCR4 repressor and exploited a target-site blocker (CXCR4-TSB) that disrupted the interaction of miR-206-3p with CXCR4 in vitro and in vivo. In vitro, CXCR4-TSB enhanced CXCR4 expression in human and murine ECs and VSMCs to modulate cell viability, proliferation, and migration. Systemic administration of CXCR4-TSB in Apoe-deficient mice enhanced Cxcr4 expression in ECs and VSMCs in the walls of blood vessels, reduced vascular permeability and monocyte adhesion to endothelium, and attenuated the development of diet-induced atherosclerosis. CXCR4-TSB also increased CXCR4 expression in B cells, corroborating its atheroprotective role in this cell type. Analyses of human atherosclerotic plaque specimens revealed a decrease in CXCR4 and an increase in miR-206-3p expression in advanced compared with early lesions, supporting a role for the miR-206-3p-CXCR4 interaction in human disease. Disrupting the miR-206-3p-CXCR4 interaction in a cell-specific manner with target-site blockers is a potential therapeutic approach that could be used to treat atherosclerosis and other vascular diseases.

Modulation of Influenza A virus NS1 expression reveals prioritization of host response antagonism at single-cell resolution

Influenza A virus (IAV) is an important human respiratory pathogen that causes significant seasonal epidemics and potential devastating pandemics. As part of its life cycle, IAV encodes the multifunctional protein NS1, that, among many roles, prevents immune detection and limits interferon (IFN) production. As distinct host immune pathways exert different selective pressures against IAV, as replication progresses, we expect a prioritization in the host immune antagonism by NS1. In this work, we profiled bulk transcriptomic differences in a primary bronchial epithelial cell model facing IAV infections at distinct NS1 levels. We further demonstrated that, at single cell level, the intracellular amount of NS1 in-part shapes the heterogeneity of the host response. We found that modulation of NS1 levels reveal a ranking in its inhibitory roles: modest NS1 expression is sufficient to inhibit immune detection, and thus the expression of pro-inflammatory cytokines (including IFNs), but higher levels are required to inhibit IFN signaling and ISG expression. Lastly, inhibition of chaperones related to the unfolded protein response requires the highest amount of NS1, often associated with later stages of viral replication. This work demystifies some of the multiple functions ascribed to IAV NS1, highlighting the prioritization of NS1 in antagonizing the different pathways involved in the host response to IAV infection.

AP003352.1/miR-141-3p axis enhances the proliferation of osteosarcoma by LPAR3

Osteosarcoma (OS) is a highly malignant tumor with a poor prognosis and a growing incidence. LncRNAs and microRNAs control the occurrence and development process of osteosarcoma through ceRNA patterns. The LPAR3 gene is important in cancer cell proliferation, apoptosis and disease development. However, the regulatory mechanism of the ceRNA network through which LPAR3 participates in osteosarcoma has not been clarified. Herein, our study demonstrated that the AP003352.1/miR-141-3p axis drives LPAR3 expression to induce the malignant progression of osteosarcoma. First, the expression of LPAR3 is regulated by the changes in AP003352.1 and miR-141-3p. Similar to the ceRNA of miR-141-3p, AP003352.1 regulates the expression of LPAR3 through this mechanism. In addition, the regulation of AP003352.1 in malignant osteosarcoma progression depends to a certain degree on miR-141-3p. Importantly, the AP003352.1/miR-141-3p/LPAR3 axis can better serve as a multi-gene diagnostic marker for osteosarcoma. In conclusion, our research reveals a new ceRNA regulatory network, which provides a novel potential target for the diagnosis and treatment of osteosarcoma.

Omics Overview of the SPARC Gene in Mesothelioma

The SPARC gene plays multiple roles in extracellular matrix synthesis and cell shaping, associated with tumor cell migration, invasion, and metastasis. The SPARC gene is also involved in the epithelial-mesenchymal transition (EMT) process, which is a critical phenomenon leading to a more aggressive cancer cell phenotype. SPARC gene overexpression has shown to be associated with poor survival in the mesothelioma (MESO) cohort from the TCGA database, indicating that this gene may be a powerful prognostic factor in MESO. Its overexpression is correlated with the immunosuppressive tumor microenvironment. Here, we summarize the omics advances of the SPARC gene, including the summary of SPARC gene expression associated with prognosis in pancancer and MESO, the immunosuppressive microenvironment, and cancer cell stemness. In addition, SPARC might be targeted by microRNAs. Notably, despite the controversial functions on angiogenesis, SPARC may directly or indirectly contribute to tumor angiogenesis in MESO. In conclusion, SPARC is involved in tumor invasion, metastasis, immunosuppression, cancer cell stemness, and tumor angiogenesis, eventually impacting patient survival. Strategies targeting this gene may provide novel therapeutic approaches to the treatment of MESO.

Expression Microdissection for the Analysis of miRNA in a Single-Cell Type

Cell-specific microRNA (miRNA) expression estimates are important in characterizing the localization of miRNA signaling within tissues. Much of these data are obtained from cultured cells, a process known to significantly alter miRNA expression levels. Thus, our knowledge of in vivo cell miRNA expression estimates is poor. We previously demonstrated expression microdissection-miRNA-sequencing (xMD-miRNA-seq) to acquire in vivo estimates, directly from formalin-fixed tissues, albeit with a limited yield. In this study, we optimized each step of the xMD process, including tissue retrieval, tissue transfer, film preparation, and RNA isolation, to increase RNA yields and ultimately show strong enrichment for in vivo miRNA expression by qPCR array. These method improvements, such as the development of a noncrosslinked ethylene vinyl acetate membrane, resulted in a 23- to 45-fold increase in miRNA yield, depending on the cell type. By qPCR, miR-200a increased by 14-fold in xMD-derived small intestine epithelial cells, with a concurrent 336-fold reduction in miR-143 relative to the matched nondissected duodenal tissue. xMD is now an optimized method to obtain robust in vivo miRNA expression estimates from cells. xMD will allow formalin-fixed tissues from surgical pathology archives to make theragnostic biomarker discoveries.

Mouse microglia express unique miRNA-mRNA networks to facilitate age-specific functions in the developing central nervous system

Microglia regulate multiple processes in the central nervous system, exhibiting a considerable level of cellular plasticity which is facilitated by an equally dynamic transcriptional environment. While many gene networks that regulate microglial functions have been characterised, the influence of epigenetic regulators such as small non-coding microRNAs (miRNAs) is less well defined. We have sequenced the miRNAome and mRNAome of mouse microglia during brain development and adult homeostasis, identifying unique profiles of known and novel miRNAs. Microglia express both a consistently enriched miRNA signature as well as temporally distinctive subsets of miRNAs. We generated robust miRNA-mRNA networks related to fundamental developmental processes, in addition to networks associated with immune function and dysregulated disease states. There was no apparent influence of sex on miRNA expression. This study reveals a unique developmental trajectory of miRNA expression in microglia during critical stages of CNS development, establishing miRNAs as important modulators of microglial phenotype.

Reassessment of miRNA variant (isomiRs) composition by small RNA sequencing

IsomiRs, sequence variants of mature microRNAs, are usually detected and quantified using high-throughput sequencing. Many examples of their biological relevance have been reported, but sequencing artifacts identified as artificial variants might bias biological inference and therefore need to be ideally avoided. We conducted a comprehensive evaluation of 10 different small RNA sequencing protocols, exploring both a theoretically isomiR-free pool of synthetic miRNAs and HEK293T cells. We calculated that, with the exception of two protocols, less than 5% of miRNA reads can be attributed to library preparation artifacts. Randomized-end adapter protocols showed superior accuracy, with 40% of true biological isomiRs. Nevertheless, we demonstrate concordance across protocols for selected miRNAs in non-templated uridyl additions. Notably, NTA-U calling and isomiR target prediction can be inaccurate when using protocols with poor single-nucleotide resolution. Our results highlight the relevance of protocol choice for biological isomiRs detection and annotation, which has key potential implications for biomedical applications.

Direct observation of the evolution of cell-type-specific microRNA expression signatures supports the hematopoietic origin model of endothelial cells

The evolution of specialized cell-types is a long-standing interest of biologists, but given the deep time-scales very difficult to reconstruct or observe. microRNAs have been linked to the evolution of cellular complexity and may inform on specialization. The endothelium is a vertebrate-specific specialization of the circulatory system that enabled a critical new level of vasoregulation. The evolutionary origin of these endothelial cells is unclear. We hypothesized that Mir-126, an endothelial cell-specific microRNA may be informative. We here reconstruct the evolutionary history of Mir-126. Mir-126 likely appeared in the last common ancestor of vertebrates and tunicates, which was a species without an endothelium, within an intron of the evolutionary much older EGF Like Domain Multiple (Egfl) locus. Mir-126 has a complex evolutionary history due to duplications and losses of both the host gene and the microRNA. Taking advantage of the strong evolutionary conservation of the microRNA among Olfactores, and using RNA in situ hybridization, we localized Mir-126 in the tunicate Ciona robusta. We found exclusive expression of the mature Mir-126 in granular amebocytes, supporting a long-proposed scenario that endothelial cells arose from hemoblasts, a type of proto-endothelial amoebocyte found throughout invertebrates. This observed change of expression of Mir-126 from proto-endothelial amoebocytes in the tunicate to endothelial cells in vertebrates is the first direct observation of the evolution of a cell-type in relation to microRNA expression indicating that microRNAs can be a prerequisite of cell-type evolution.

miR-126 promotes M1 to M2 macrophage phenotype switching via VEGFA and KLF4

Background

Macrophage polarization and microRNA play crucial roles in the development of atherosclerosis (AS). The M1 macrophage phenotype contributes to the formation of plaques, while the M2 macrophage phenotype resolves inflammation and promotes tissue repair. MiR-126 has been found to play a role in regulating macrophage polarization in the context of AS. However, the exact mechanism of miR-126 requires further research.

Methods

The foam cell model was established by stimulating THP-1 with oxidized low-density lipoprotein (ox-LDL). We transfected foam cells with miR-126 mimic and its negative control. The transfection of miR-126 was implemented by riboFECT CP transfection kit. The levels of miR-126 and M1/M2 associated genes in foam cells were quantified using reverse transcription-quantitative PCR (RT-qPCR). Additionally, the expressions of CD86+ and CD206+ cells in foam cells were determined by flow cytometry. Western blotting and RT-qPCR were used to determine the protein and mRNA levels of the vascular endothelial growth factor A (VEGFA) and the transcriptional regulator Krüppel-like factor 4 (KLF4), respectively. Additionally, we detected endothelial cell migration after co-culturing endothelial cells and macrophages. MG-132 was used to indirectly activate the expression of VEGFA, and the expression of KLF4 was also evaluated.

Results

The activation of apoptosis and production of foam cells were boosted by the addition of ox-LDL. We transfected foam cells with miR-126 mimic and its negative control and observed that miR-126 greatly suppressed foam cell development and inhibited phagocytosis. Moreover, it caused pro-inflammatory M1 macrophages to switch to the anti-inflammatory M2 phenotype. This was reflected by the increase in anti-inflammatory gene expression and the decrease in pro-inflammatory gene expression. Additionally, miR-126 dramatically decreased the expressions of VEGFA and KLF4. The protein-protein interaction network analysis showed a significantly high correlation between miR-126, VEGFA, and KLF4. MiR-126 may also promote EC migration by activating macrophage PPAR γ expression and effectively suppressing macrophage inflammation. MG-132 indirectly activated the expression of VEGFA, and the expression of KLF4 also significantly increased, which indicates a direct or indirect relationship between VEGFA and KLF4.

Conclusion

Our study shows that miR-126 can reverse ox-LDL-mediated phagocytosis and apoptosis in macrophages. Consequently, the potential role of miR-126 was manifested in regulating macrophage function and promoting vascular endothelial migration.

Dialog beyond the Grave: Necrosis in the Tumor Microenvironment and Its Contribution to Tumor Growth

Damage-associated molecular patterns (DAMPs) are endogenous molecules released from the necrotic cells dying after exposure to various stressors. After binding to their receptors, they can stimulate various signaling pathways in target cells. DAMPs are especially abundant in the microenvironment of malignant tumors and are suspected to influence the behavior of malignant and stromal cells in multiple ways often resulting in promotion of cell proliferation, migration, invasion, and metastasis, as well as increased immune evasion. This review will start with a reminder of the main features of cell necrosis, which will be compared to other forms of cell death. Then we will summarize the various methods used to assess tumor necrosis in clinical practice including medical imaging, histopathological examination, and/or biological assays. We will also consider the importance of necrosis as a prognostic factor. Then the focus will be on the DAMPs and their role in the tumor microenvironment (TME). We will address not only their interactions with the malignant cells, frequently leading to cancer progression, but also with the immune cells and their contribution to immunosuppression. Finally, we will emphasize the role of DAMPs released by necrotic cells in the activation of Toll-like receptors (TLRs) and the possible contributions of TLRs to tumor development. This last point is very important for the future of cancer therapeutics since there are attempts to use TLR artificial ligands for cancer therapeutics.

Updated review of advances in microRNAs and complex diseases: experimental results, databases, webservers and data fusion

MicroRNAs (miRNAs) are gene regulators involved in the pathogenesis of complex diseases such as cancers, and thus serve as potential diagnostic markers and therapeutic targets. The prerequisite for designing effective miRNA therapies is accurate discovery of miRNA-disease associations (MDAs), which has attracted substantial research interests during the last 15 years, as reflected by more than 55 000 related entries available on PubMed. Abundant experimental data gathered from the wealth of literature could effectively support the development of computational models for predicting novel associations. In 2017, Chen et al. published the first-ever comprehensive review on MDA prediction, presenting various relevant databases, 20 representative computational models, and suggestions for building more powerful ones. In the current review, as the continuation of the previous study, we revisit miRNA biogenesis, detection techniques and functions; summarize recent experimental findings related to common miRNA-associated diseases; introduce recent updates of miRNA-relevant databases and novel database releases since 2017, present mainstream webservers and new webserver releases since 2017 and finally elaborate on how fusion of diverse data sources has contributed to accurate MDA prediction.

miRSCAPE - inferring miRNA expression from scRNA-seq data

Our understanding of miRNA activity at cellular resolution is thwarted by the inability of standard scRNA-seq protocols to capture miRNAs. We introduce a novel tool, miRSCAPE, to infer miRNA expression in a sample from its RNA-seq profile. We establish miRSCAPE's accuracy in 10 tumor and normal cohorts demonstrating its superiority over alternatives. miRSCAPE accurately infers cell type-specific miRNA activities (predicted versus observed fold-difference correlation ∼0.81) in two independent scRNA-seq datasets. We apply miRSCAPE to infer miRNA activities in scRNA clusters in pancreatic and lung adenocarcinomas, as well as in 56 cell types in the human cell landscape (HCL). In pancreatic and breast cancer scRNA-seq data, miRSCAPE recapitulates miRNAs associated with stemness and epithelial-mesenchymal transition (EMT) cell states, respectively. Overall, miRSCAPE recapitulates and refines miRNA biology at cellular resolution. miRSCAPE is freely available and is easily applicable to scRNA-seq data to infer miRNA activities at cellular resolution.

TENT2, TUT4, and TUT7 selectively regulate miRNA sequence and abundance

TENTs generate miRNA isoforms by 3' tailing. However, little is known about how tailing regulates miRNA function. Here, we generate isogenic HEK293T cell lines in which TENT2, TUT4 and TUT7 are knocked out individually or in combination. Together with rescue experiments, we characterize TENT-specific effects by deep sequencing, Northern blot and in vitro assays. We find that 3' tailing is not random but highly specific. In addition to its known adenylation, TENT2 contributes to guanylation and uridylation on mature miRNAs. TUT4 uridylates most miRNAs whereas TUT7 is dispensable. Removing adenylation has a marginal impact on miRNA levels. By contrast, abolishing uridylation leads to dysregulation of a set of miRNAs. Besides let-7, miR-181b and miR-222 are negatively regulated by TUT4/7 via distinct mechanisms while the miR-888 cluster is upregulated specifically by TUT7. Our results uncover the selective actions of TENTs in generating 3' isomiRs and pave the way to investigate their functions.

A curated human cellular microRNAome based on 196 primary cell types

BACKGROUND

An incomplete picture of the expression distribution of microRNAs (miRNAs) across human cell types has long hindered our understanding of this important regulatory class of RNA. With the continued increase in available public small RNA sequencing datasets, there is an opportunity to more fully understand the general distribution of miRNAs at the cell level.

RESULTS

From the NCBI Sequence Read Archive, we obtained 6,054 human primary cell datasets and processed 4,184 of them through the miRge3.0 small RNA sequencing alignment software. This dataset was curated down, through shared miRNA expression patterns, to 2,077 samples from 196 unique cell types derived from 175 separate studies. Of 2,731 putative miRNAs listed in miRBase (v22.1), 2,452 (89.8%) were detected. Among reasonably expressed miRNAs, 108 were designated as cell specific/near specific, 59 as infrequent, 52 as frequent, 54 as near ubiquitous, and 50 as ubiquitous. The complexity of cellular microRNA expression estimates recapitulates tissue expression patterns and informs on the miRNA composition of plasma.

CONCLUSIONS

This study represents the most complete reference, to date, of miRNA expression patterns by primary cell type. The data are available through the human cellular microRNAome track at the UCSC Genome Browser (https://genome.ucsc.edu/cgi-bin/hgHubConnect) and an R/Bioconductor package (https://bioconductor.org/packages/microRNAome/).

FBP1 /miR-24-1/enhancer axis activation blocks renal cell carcinoma progression via Warburg effect

Warburg effect is a pivotal hallmark of cancers and appears prevalently in renal cell carcinoma (RCC). FBP1 plays a negative role in Warburg effect as a rate-limiting enzyme in gluconeogenesis, yet its mechanism in RCC remains to be further characterized. Herein, we revealed that FBP1 was downregulated in RCC tissue samples and was related to the poor survival rate of RCC. Strikingly, miR-24-1 whose DNA locus is overlapped with enhancer region chr9:95084940-95087024 was closely linked with the depletion of FBP1 in RCC. Of note, miRNAs like miR-24-1 whose DNA loci are enriched with H3K27ac and H3K4me1 modifications are belonging to nuclear activating miRNAs (NamiRNAs), which surprisingly upregulate target genes in RCC through enhancer beyond the conventional role of repressing target gene expression. Moreover, miR-24-1 reactivated the expression of FBP1 to suppress Warburg effect in RCC cells, and subsequently inhibited proliferation and metastasis of RCC cells. In mechanism, the activating role of miR-24-1 was dependent on enhancer integrity by dual luciferase reporter assay and CRISPR/Cas9 system. Ultimately, animal assay in vivo validated the suppressive function of FBP1 on 786-O and ACHN cells. Collectively, the current study highlighted that activation of FBP1 by enhancer-overlapped miR-24-1 is capable of contributing to Warburg effect repression through which RCC progression is robustly blocked, providing an alternative mechanism for RCC development and as well implying a potential clue for RCC treatment strategy.

hnRNPC induces isoform shifts in miR-21-5p leading to cancer development

MicroRNA (miRNA) processing is a critical step in mature miRNA production. Its dysregulation leads to an increase in miRNA isoforms with heterogenous 5'-ends (isomiRs), which can recognize distinct target sites because of their shifted seed sequence. Although some miRNA genes display productive expression of their 5'-isomiRs in cancers, how their production is controlled and how 5'-isomiRs affect tumor progression have yet to be explored. In this study, based on integrative analyses of high-throughput sequencing data produced by our group and publicly available data, we demonstrate that primary miR-21 (pri-miR-21) is processed into the cancer-specific isomiR isomiR-21-5p | ±1, which suppresses growth hormone receptor (GHR) in liver cancer. Treatment with antagomirs against isomiR-21-5p | ±1 inhibited the in vitro tumorigenesis of liver cancer cells and allowed the recovery of GHR, whereas the introduction of isomiR-21-5p | ±1 mimics attenuated these effects. These effects were validated in a mouse model of spontaneous liver cancer. Heterogeneous nuclear ribonucleoprotein C and U2 small nuclear RNA auxiliary factor 2 were predicted to bind upstream of pre-miR-21 via a poly-(U) motif and influence Drosha processing to induce the production of isomiR-21-5p | ±1. Our findings suggest an oncogenic function for the non-canonical isomiR-21-5p | ±1 in liver cancer, and its production was shown to be regulated by hnRNPC.

A plant immune protein enables broad antitumor response by rescuing microRNA deficiency

Cancer cells are featured with uncontrollable activation of cell cycle, and microRNA deficiency drives tumorigenesis. The RNA-dependent RNA polymerase (RDR) is essential for small-RNA-mediated immune response in plants but is absent in vertebrates. Here, we show that ectopic expression of plant RDR1 can generally inhibit cancer cell proliferation. In many human primary tumors, abnormal microRNA isoforms with 1-nt-shorter 3' ends are widely accumulated. RDR1 with nucleotidyltransferase activity can recognize and modify the problematic AGO2-free microRNA duplexes with mononucleotides to restore their 2 nt overhang structure, which eventually rescues AGO2-loading efficiency and elevates global miRNA expression to inhibit cancer cell-cycle specifically. The broad antitumor effects of RDR1, which can be delivered by an adeno-associated virus, are visualized in multiple xenograft tumor models in vivo. Altogether, we reveal the widespread accumulation of aberrant microRNA isoforms in tumors and develop a plant RDR1-mediated antitumor stratagem by editing and repairing defective microRNAs.

miRNA/mRNA co-profiling identifies the miR-200 family as a central regulator of SMC quiescence

miRNAs are versatile regulators of smooth muscle cell (SMC) fate and behavior in vascular development and disease. Targeted loss-of-function studies have established the relevance of specific miRNAs in controlling SMC differentiation or mediating phenotypic modulation. Our goal was to characterize SMC miRNAome and its contribution to transcriptome changes during phenotypic modulation. Small RNA sequencing revealed that dedifferentiation led to the differential expression of over 50 miRNAs in cultured SMC. miRNA/mRNA comparison predicted that over a third of SMC transcript expression was regulated by differentially expressed miRNAs. Our screen identified the miR-200 cluster as highly downregulated during dedifferentiation. miR-200 maintains SMC quiescence and represses proliferation, migration, and neointima formation, in part by targeting Quaking, a central SMC phenotypic switching mediator. Our study unraveled the substantial contribution of miRNAs in regulating the SMC transcriptome and identified the miR-200 cluster as a pro-quiescence mechanism and a potential inhibitor of vascular restenosis.

Are miRNAs Dynamic Biomarkers in Keratoconus? A Review of the Literature

Aim: A review of miRNA (microRNA) profiling studies in keratoconus. Methods: Literature search strategy—PubMed central database, using miRNA or microRNA and keratoconus as keywords. Results: Eleven experimental or clinical studies on humans regarding miRNA and keratoconus, published in English between 2009 and 2020 were retrieved. Conclusion: The publications regarding the role of miRNAs in keratoconus are scarce and diverse but provide some valuable information about potential new mechanisms of keratoconus development and progression. The cornea expresses almost 300 different miRNAs, 18 of which are specific, and miR-184 is by far the most abundant, with expression restricted to central basal and suprabasal epithelial cells. Mutations in the seed region of MIR184 were proved to be rare and nonspecific in patients with isolated keratoconus. Overall, in keratoconus, a total of 29 miRNAs were upregulated, and 11 were downregulated. It appeared that miR-143-3p, miR-182-5p, and miR-92a-3p were highly expressed, while the miRNAs connected to cell–cell junction, cell division, and motor activity were downregulated. In less advanced forms, altered expression of four miRNAs—miR-151a-3p, miR-194-5p, miR-195-5p, miR-185-5p—was proved in the cone epithelium; in contrast, in advanced keratoconus, the expression of miR-151a-3p and miR-194-5p remained altered, changes in the expression of miR-195 and miR-185 were not reported, and the expression of miR-138-5p, miR-146b-5p, miR-28-5p, and miR-181a-2-3p was also altered in the corneal epithelium. Keratoconus is a dynamic process of corneal stromal thinning that might result from a dynamic miRNA expression in the corneal epithelium exposed to environmental and behavioral factors causing repetitive traumas. Further experimental studies are needed to prove this hypothesis.

MicroRNA Regulation of Bone Marrow Mesenchymal Stem Cell Chondrogenesis: Toward Articular Cartilage

The production of a clinically useful engineered cartilage is an outstanding and unmet clinical need. High-throughput RNA sequencing provides a means of characterizing the molecular phenotype of populations of cells and can be leveraged to better understand differences among source cells, derivative engineered tissues, and target phenotypes. In this study, small RNA sequencing is utilized to comprehensively characterize the microRNA transcriptomes (miRNomes) of native human neonatal articular cartilage and human bone marrow-derived mesenchymal stem cells (hBM-MSCs) differentiating into cartilage organoids, contrasting the microRNA regulation of engineered cartilage with that of a promising target phenotype. Five dominant microRNAs are upregulated during cartilage organoid differentiation and disproportionately regulate transcription factors: miR-148a-3p, miR-140-3p, miR-27b-3p, miR-140-5p, and miR-181a-5p. Two microRNAs that dominate the miRNomes of hBM-MSCs, miR-21-5p and miR-143-3p, persist throughout the differentiation process and may limit the ability of these cells to differentiate into an engineered cartilage resembling target native articular cartilage. By using predictive bioinformatics tools and antagomir inhibition, these persistent microRNAs are shown to destabilize the mRNA of genes with known or potential roles in cartilage biology including FGF18, TGFBR2, TET1, STOX2, ARAP2, N4BP2L1, LHX9, NFIA, and RPS6KA5. These results shed light on the extent to which only a few microRNAs contribute to the complex regulatory environment of hBM-MSCs for engineered tissues. Impact statement MicroRNAs are emerging as important controlling elements in the differentiation of human bone marrow-derived mesenchymal stem cells (hBM-MSCs). By using a robust bioinformatic approach and further validation in vitro, here we provide a comprehensive characterization of the microRNA transcriptomes (miRNomes) of a commonly studied and clinically promising source of multipotent cells (hBM-MSCs), a gold standard model of in vitro chondrogenesis (hBM-MSC-derived cartilage organoids), and an attractive in vivo target phenotype for clinically useful engineered cartilage (neonatal articular cartilage). These analyses highlighted a specific set of microRNAs involved in the chondrogenic program that could be manipulated to acquire a more robust articular cartilage-like phenotype. This characterization provides researchers in the cartilage tissue engineering field a useful atlas with which to contextualize microRNA involvement in complex differentiation pathways.

Flexible pri-miRNA structures enable tunable production of 5’ isomiRs

The Drosha cleavage of a pri-miRNA defines mature microRNA sequence. Drosha cleavage at alternative positions generates 5’ isoforms (isomiRs) which have distinctive functions. To understand how pri-miRNA structures influence Drosha cleavage, we performed a systematic analysis of the maturation of endogenous pri-miRNAs and their variants both in vitro and in vivo. We show that in addition to previously known features, the overall structural flexibility of pri-miRNA impact Drosha cleavage fidelity. Internal loops and nearby G · U wobble pairs on the pri-miRNA stem induce the use of non-canonical cleavage sites by Drosha, resulting in 5’ isomiR production. By analysing patient data deposited in the Cancer Genome Atlas, we provide evidence that alternative Drosha cleavage of pri-miRNAs is a tunable process that responds to the level of pri-miRNA-associated RNA-binding proteins. Together, our findings reveal that Drosha cleavage fidelity can be modulated by altering pri-miRNA structure, a potential mechanism underlying 5’ isomiR biogenesis in tumours.

Graphical Abstract

Blockade of the pro-fibrotic reaction mediated by the miR-143/-145 cluster enhances the responses to targeted therapy in melanoma

Lineage dedifferentiation toward a mesenchymal‐like state displaying myofibroblast and fibrotic features is a common mechanism of adaptive and acquired resistance to targeted therapy in melanoma. Here, we show that the anti‐fibrotic drug nintedanib is active to normalize the fibrous ECM network, enhance the efficacy of MAPK‐targeted therapy, and delay tumor relapse in a preclinical model of melanoma. Acquisition of this resistant phenotype and its reversion by nintedanib pointed to miR‐143/‐145 pro‐fibrotic cluster as a driver of this mesenchymal‐like phenotype. Upregulation of the miR‐143/‐145 cluster under BRAFi/MAPKi therapy was observed in melanoma cells in vitro and in vivo and was associated with an invasive/undifferentiated profile. The 2 mature miRNAs generated from this cluster, miR‐143‐3p and miR‐145‐5p, collaborated to mediate transition toward a drug‐resistant undifferentiated mesenchymal‐like state by targeting Fascin actin‐bundling protein 1 (FSCN1), modulating the dynamic crosstalk between the actin cytoskeleton and the ECM through the regulation of focal adhesion dynamics and mechanotransduction pathways. Our study brings insights into a novel miRNA‐mediated regulatory network that contributes to non‐genetic adaptive drug resistance and provides proof of principle that preventing MAPKi‐induced pro‐fibrotic stromal response is a viable therapeutic opportunity for patients on targeted therapy.

A comprehensive framework for analysis of microRNA sequencing data in metastatic colorectal cancer

Although microRNAs (miRNAs) contribute to all hallmarks of cancer, miRNA dysregulation in metastasis remains poorly understood. The aim of this work was to reliably identify miRNAs associated with metastatic progression of colorectal cancer (CRC) using novel and previously published next-generation sequencing (NGS) datasets generated from 268 samples of primary (pCRC) and metastatic CRC (mCRC; liver, lung and peritoneal metastases) and tumor adjacent tissues. Differential expression analysis was performed using a meticulous bioinformatics pipeline, including only bona fide miRNAs, and utilizing miRNA-tailored quality control and processing. Five miRNAs were identified as up-regulated at multiple metastatic sites Mir-210_3p, Mir-191_5p, Mir-8-P1b_3p [mir-141–3p], Mir-1307_5p and Mir-155_5p. Several have previously been implicated in metastasis through involvement in epithelial-to-mesenchymal transition and hypoxia, while other identified miRNAs represent novel findings. The use of a publicly available pipeline facilitates reproducibility and allows new datasets to be added as they become available. The set of miRNAs identified here provides a reliable starting-point for further research into the role of miRNAs in metastatic progression.

The Role of miRNAs as New Molecular Biomarkers for Dating the Age of Wound Production: A Systematic Review

BACKGROUND

The timing of wounds production is a significant issue in forensic pathology. Although various methods have been evaluated, obtaining an accurate dating of lesions is still a challenge. The pathologist uses many parameters to value wound age, such as histological and immunohistochemical. In recent years, there have been many studies regarding the use of miRNAs in wound-age estimation; indeed, miRNAs have multiple potential uses in forensic pathology.

SCOPE

This review aims to verify the efficacy and feasibility of miRNAs as a tool for determining the timing of lesions.

MATERIALS AND METHODS

The authors conducted the systematic review according to the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines. PubMed was used as a search engine to find articles published between January, 1st 2016 and October, 1st 2021, to evaluate the current state of the art regarding wound-age estimation.

RESULTS

A total of 256 articles were collected; after screening according to PRISMA guidelines, the systematic review included 8 articles. The studies included in this review were all Original articles evaluating the use of biomarkers for wound-age determination.

DISCUSSION AND CONCLUSION

The literature review showed that analysis of miRNA is an innovative field of study with significant potentiality in forensic pathology. There are few studies, and almost all of them are at an early stage. The challenge is to understand how to standardize the samples' selection to obtain reliable experimental data. This observation represents a necessary prerequisite to planning further clinical trials.

OUP accepted manuscript

Alternative polyadenylation increases transcript diversities at the 3’ end, regulating biological processes including cell differentiation, embryonic development and cancer progression. Here, we present a Bayesian method SCAPE, which enables de novo identification and quantification of polyadenylation (pA) sites at single-cell level by utilizing insert size information. We demonstrated its accuracy and robustness and identified 31 558 sites from 36 mouse organs, 43.8% (13 807) of which were novel. We illustrated that APA isoforms were associated with miRNAs binding and regulated in tissue-, cell type-and tumor-specific manners where no difference was found at gene expression level, providing an extra layer of information for cell clustering. Furthermore, we found genome-wide dynamic changes of APA usage during erythropoiesis and induced pluripotent stem cell (iPSC) differentiation, suggesting APA contributes to the functional flexibility and diversity of single cells. We expect SCAPE to aid the analyses of cellular dynamics and diversities in health and disease.

Pre-Analytical Modification of Serum miRNAs: Diagnostic Reliability of Serum miRNAs in Hemolytic Diseases

Circulating microRNAs (miRNAs) are useful biomarkers of hemolysis. Since blood cells are the main origins of circulating miRNAs, we evaluated blood cell-related pre-analytical modification of the miRNA signatures during blood drawing and serum processing. The levels of miRNA before and after ex vivo blood drawing were analyzed with the reverse transcriptase-based polymerase chain reaction method. Furthermore, the changes of miRNA signatures caused by different time-lag between blood drawing and serum preparation by 24 h were evaluated. Finally, we compared the miRNA levels between leftover samples and samples of hemolytic diseases. Blood drawing procedure induced increments of red blood cell (RBC)-related miRNAs (miR-451a, miR-486) about 2-fold. One hour standing of blood samples before serum separation induced almost the same increases in RBC-related miRNAs. To test the clinical usefulness of miR-451a as a biomarker of hemolytic diseases, we analyzed miRNAs of samples from 10 normal subjects, 30 leftover samples in the clinical laboratory, and 20 samples from patients with hemolytic diseases. Serum miR-451a significantly increased in patients with hemolytic anemia more than the levels of pre-analytical modification. In conclusion, the pre-analytical modification of serum miRNAs did not disturb the usefulness of RBC-derived miRNAs as biomarkers of hemolytic diseases.

Sarcopenia: Etiology, Nutritional Approaches, and miRNAs

Sarcopenia, an age-related decline in skeletal muscle mass and function, dramatically affects the quality of life. Although there is a consensus that sarcopenia is a multifactorial syndrome, the etiology and underlying mechanisms are not yet delineated. Moreover, research about nutritional interventions to prevent the development of sarcopenia is mainly focused on the amount and quality of protein intake. The impact of several nutrition strategies that consider timing of food intake, anti-inflammatory nutrients, metabolic control, and the role of mitochondrial function on the progression of sarcopenia is not fully understood. This narrative review summarizes the metabolic background of this phenomenon and proposes an integral nutritional approach (including dietary supplements such as creatine monohydrate) to target potential molecular pathways that may affect reduce or ameliorate the adverse effects of sarcopenia. Lastly, miRNAs, in particular those produced by skeletal muscle (MyomiR), might represent a valid tool to evaluate sarcopenia progression as a potential rapid and early biomarker for diagnosis and characterization.

New technologies for improved relevance in miRNA research

After a number of years of research in the field of miRNA, the robustness and biological relevance of many published articles is increasingly being questioned. We propose the use of new RNA-seq approaches, genome editing technologies, and updated public databases to improve the quality, reliability, and relevance of published data.

miRge3.0: a comprehensive microRNA and tRF sequencing analysis pipeline

MicroRNAs and tRFs are classes of small non-coding RNAs, known for their roles in translational regulation of genes. Advances in next-generation sequencing (NGS) have enabled high-throughput small RNA-seq studies, which require robust alignment pipelines. Our laboratory previously developed miRge and miRge2.0, as flexible tools to process sequencing data for annotation of miRNAs and other small-RNA species and further predict novel miRNAs using a support vector machine approach. Although miRge2.0 is a leading analysis tool in terms of speed with unique quantifying and annotation features, it has a few limitations. We present miRge3.0 that provides additional features along with compatibility to newer versions of Cutadapt and Python. The revisions of the tool include the ability to process Unique Molecular Identifiers (UMIs) to account for PCR duplicates while quantifying miRNAs in the datasets, correct erroneous single base substitutions in miRNAs with miREC and an accurate mirGFF3 formatted isomiR tool. miRge3.0 also has speed improvements benchmarked to miRge2.0, Chimira and sRNAbench. Finally, miRge3.0 output integrates into other packages for a streamlined analysis process and provides a cross-platform Graphical User Interface (GUI). In conclusion miRge3.0 is our third generation small RNA-seq aligner with improvements in speed, versatility and functionality over earlier iterations.

Circulating miRNAs as Novel Non-Invasive Biomarkers to Aid the Early Diagnosis of Suspicious Breast Lesions for Which Biopsy Is Recommended

Simple Summary

In population-based screens, tissue biopsy remains the standard practice for women with imaging that suggests breast cancer. We examined circulating microRNAs as minimally invasive diagnostic biomarkers to discriminate malignant from benign breast lesions. A retrospective cohort of plasma samples divided into training and testing sets and a prospective cohort of women with suspicious imaging findings who underwent tissue biopsy were investigated through a global microRNA profile by OpenArray. Seven signatures, involving 5 specific miRNAs (miR-625, miR-423-5p, miR-370-3p, miR-181c, and miR-301b), were identified and validated in the testing set. Among the 7 signatures, the discriminatory performances of 5 of them were confirmed in the prospective cohort.

Abstract

In population-based screens, tissue biopsy remains the standard practice for women with imaging that suggests breast cancer. We examined circulating microRNAs as minimally invasive diagnostic biomarkers to discriminate malignant from benign breast lesions. miRNAs were analyzed by OpenArray in a retrospective cohort of plasma samples including 100 patients with malignant (T), 89 benign disease (B), and 99 healthy donors (HD) divided into training and testing sets and a prospective cohort (BABE) of 289 women with suspicious imaging findings who underwent tissue biopsy. miRNAs associated with disease status were identified by univariate analysis and then combined into signatures by multivariate logistic regression models. By combining 16 miRNAs differentially expressed in the T vs. HD comparison, 26 signatures were also able to significantly discriminate T from B disease. Seven of them, involving 5 specific miRNAs (miR-625, miR-423-5p, miR-370-3p, miR-181c, and miR-301b), were statistically validated in the testing set. Among the 7 signatures, the discriminatory performances of 5 were confirmed in the prospective BABE Cohort. This study identified 5 circulating miRNAs that, properly combined, distinguish malignant from benign breast disease in women with a high likelihood of malignancy.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

Profiling of Primary and Mature miRNA Expression in Atherosclerosis-Associated Cell Types

[Figure: see text].

COPD-associated miR-145-5p is downregulated in early-decline FEV1 trajectories in childhood asthma

BACKGROUND

Many microRNAs (miRNAs) have been associated with asthma and chronic obstructive pulmonary disease (COPD). Longitudinal lung function growth trajectories of children with asthma-normal growth, reduced growth (RG), early decline (ED), and RG with an ED (RGED)-have been observed, with RG and RGED associated with adverse outcomes, including COPD.

OBJECTIVE

Our aim was to determine whether circulating miRNAs from an early age in children with asthma would be prognostic of reduced lung function growth patterns over the next 16 years.

METHODS

We performed small RNA sequencing on sera from 492 children aged 5 to 12 years with mild-to-moderate asthma from the CAMP clinical trial, who were subsequently followed for 12 to 16 years. miRNAs were assessed for differential expression between previously assigned lung function growth patterns.

RESULTS

We had 448 samples and 259 miRNAs for differential analysis. In a comparison of the normal and the most severe group (ie, normal growth compared with RGED), we found 1 strongly dysregulated miRNA, hsa-miR-145-5p (P < 8.01E-05). This miR was downregulated in both ED groups (ie, ED and RGED). We verified that miR-145-5p was strongly associated with airway smooth muscle cell growth in vitro.

CONCLUSION

Our results showed that miR-145-5p is associated with the ED patterns of lung function growth leading to COPD in children with asthma and additionally increases airway smooth muscle cell proliferation. This represents a significant extension of our understanding of the role of miR-145-5p in COPD and suggests that reduced expression of miR-145-5p is a risk factor for ED of long-term lung function.

Crosstalk between microRNA expression and DNA methylation drives the hormone-dependent phenotype of breast cancer

BACKGROUND

Abnormal DNA methylation is observed as an early event in breast carcinogenesis. However, how such alterations arise is still poorly understood. microRNAs (miRNAs) regulate gene expression at the post-transcriptional level and play key roles in various biological processes. Here, we integrate miRNA expression and DNA methylation at CpGs to study how miRNAs may affect the breast cancer methylome and how DNA methylation may regulate miRNA expression.

METHODS

miRNA expression and DNA methylation data from two breast cancer cohorts, Oslo2 (n = 297) and The Cancer Genome Atlas (n = 439), were integrated through a correlation approach that we term miRNA-methylation Quantitative Trait Loci (mimQTL) analysis. Hierarchical clustering was used to identify clusters of miRNAs and CpGs that were further characterized through analysis of mRNA/protein expression, clinicopathological features, in silico deconvolution, chromatin state and accessibility, transcription factor binding, and long-range interaction data.

RESULTS

Clustering of the significant mimQTLs identified distinct groups of miRNAs and CpGs that reflect important biological processes associated with breast cancer pathogenesis. Notably, two major miRNA clusters were related to immune or fibroblast infiltration, hence identifying miRNAs associated with cells of the tumor microenvironment, while another large cluster was related to estrogen receptor (ER) signaling. Studying the chromatin landscape surrounding CpGs associated with the estrogen signaling cluster, we found that miRNAs from this cluster are likely to be regulated through DNA methylation of enhancers bound by FOXA1, GATA2, and ER-alpha. Further, at the hub of the estrogen cluster, we identified hsa-miR-29c-5p as negatively correlated with the mRNA and protein expression of DNA methyltransferase DNMT3A, a key enzyme regulating DNA methylation. We found deregulation of hsa-miR-29c-5p already present in pre-invasive breast lesions and postulate that hsa-miR-29c-5p may trigger early event abnormal DNA methylation in ER-positive breast cancer.

CONCLUSIONS

We describe how miRNA expression and DNA methylation interact and associate with distinct breast cancer phenotypes.