Tools for Sequence-Based miRNA Target Prediction: What to Choose?

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.

Small Extracellular Vesicles From Infarcted and Failing Heart Accelerate Tumor Growth

BACKGROUND

Myocardial infarction (MI) and heart failure are associated with an increased incidence of cancer. However, the mechanism is complex and unclear. Here, we aimed to test our hypothesis that cardiac small extracellular vesicles (sEVs), particularly cardiac mesenchymal stromal cell-derived sEVs (cMSC-sEVs), contribute to the link between post-MI left ventricular dysfunction (LVD) and cancer.

METHODS

We purified and characterized sEVs from post-MI hearts and cultured cMSCs. Then, we analyzed cMSC-EV cargo and proneoplastic effects on several lines of cancer cells, macrophages, and endothelial cells. Next, we modeled heterotopic and orthotopic lung and breast cancer tumors in mice with post-MI LVD. We transferred cMSC-sEVs to assess sEV biodistribution and its effect on tumor growth. Finally, we tested the effects of sEV depletion and spironolactone treatment on cMSC-EV release and tumor growth.

RESULTS

Post-MI hearts, particularly cMSCs, produced more sEVs with proneoplastic cargo than nonfailing hearts did. Proteomic analysis revealed unique protein profiles and higher quantities of tumor-promoting cytokines, proteins, and microRNAs in cMSC-sEVs from post-MI hearts. The proneoplastic effects of cMSC-sEVs varied with different types of cancer, with lung and colon cancers being more affected than melanoma and breast cancer cell lines. Post-MI cMSC-sEVs also activated resting macrophages into proangiogenic and protumorigenic states in vitro. At 28-day follow-up, mice with post-MI LVD developed larger heterotopic and orthotopic lung tumors than did sham-MI mice. Adoptive transfer of cMSC-sEVs from post-MI hearts accelerated the growth of heterotopic and orthotopic lung tumors, and biodistribution analysis revealed accumulating cMSC-sEVs in tumor cells along with accelerated tumor cell proliferation. sEV depletion reduced the tumor-promoting effects of MI, and adoptive transfer of cMSC-sEVs from post-MI hearts partially restored these effects. Finally, spironolactone treatment reduced the number of cMSC-sEVs and suppressed tumor growth during post-MI LVD.

CONCLUSIONS

Cardiac sEVs, specifically cMSC-sEVs from post-MI hearts, carry multiple protumorigenic factors. Uptake of cMSC-sEVs by cancer cells accelerates tumor growth. Treatment with spironolactone significantly reduces accelerated tumor growth after MI. Our results provide new insight into the mechanism connecting post-MI LVD to cancer and propose a translational option to mitigate this deadly association.

circAMN1-Mediated Ferroptosis Regulates the Expulsion of Placenta in Trophoblast Cells

After delivery, the death of trophoblast cells can promote the expulsion of the placenta. Ferroptosis, an iron-dependent programmed cell death, is involved in mammalian development. Circular RNAs are associated with placental development; however, it is unclear whether circular RNAs regulate the expulsion of fetal membranes through ferroptosis. The gene expression profiles in the tail vein blood of Holstein cows with normal and retained placentas were investigated using RNA sequencing and a GSE214588 dataset. circAMN1 and SLC39A8 expression was significantly downregulated in the blood of cows with a retained placenta, whereas miR-205_R-1 expression was significantly upregulated. We validated erastin-induced ferroptosis in trophoblast cells. Transfection with si-circAMN1 and miR-205_R-1 mimic reduced intracellular total iron, Fe2+, and glutathione disulfide levels; increased intracellular glutathione levels and glutathione/glutathione disulfide; and enhanced cell viability in these cells. In contrast, transfection with pcDNA3.1 circAMN1 and an miR-205_R-1 inhibitor promoted ferroptosis. As an miR-205_R-1 sponge, circAMN1 regulated the expression of SLC39A8 to control erastin-induced ferroptosis and regulated the proliferation, invasion, and migration of trophoblast cells. Our findings provide a theoretical basis for studying the mechanism by which programmed cell death regulates fetal membrane expulsion and indicate its potential as a therapeutic target for placenta retention.

Integrating functional scoring and regulatory data to predict the effect of non-coding SNPs in a complex neurological disease

Most SNPs associated with complex diseases seem to lie in non-coding regions of the genome; however, their contribution to gene expression and disease phenotype remains poorly understood. Here, we established a workflow to provide assistance in prioritising the functional relevance of non-coding SNPs of candidate genes as susceptibility loci in polygenic neurological disorders. To illustrate the applicability of our workflow, we considered the multifactorial disorder migraine as a model to follow our step-by-step approach. We annotated the overlap of selected SNPs with regulatory elements and assessed their potential impact on gene expression based on publicly available prediction algorithms and functional genomics information. Some migraine risk loci have been hypothesised to reside in non-coding regions and to be implicated in the neurotransmission pathway. In this study, we used a set of 22 non-coding SNPs from neurotransmission and synaptic machinery-related genes previously suggested to be involved in migraine susceptibility based on our candidate gene association studies. After prioritising these SNPs, we focused on non-reported ones that demonstrated high regulatory potential: (1) VAMP2_rs1150 (3' UTR) was predicted as a target of hsa-mir-5010-3p miRNA, possibly disrupting its own gene expression; (2) STX1A_rs6951030 (proximal enhancer) may affect the binding affinity of zinc-finger transcription factors (namely ZNF423) and disturb TBL2 gene expression; and (3) SNAP25_rs2327264 (distal enhancer) expected to be in a binding site of ONECUT2 transcription factor. This study demonstrated the applicability of our practical workflow to facilitate the prioritisation of potentially relevant non-coding SNPs and predict their functional impact in multifactorial neurological diseases.

Whole-transcriptome sequencing revealed the role of noncoding RNAs in susceptibility and resistance of Pekin ducks to DHAV-3

As the most prevalent pathogen of duck viral hepatitis (DVH), duck hepatitis A virus genotype 3 (DHAV-3) has caused huge economic losses to the duck industry in China. Herein, we obtained whole-transcriptome sequencing data of susceptible (S) and resistant (R) Pekin duckling samples at 0 h, 12 h, and 24 h after DHAV-3 infection. We found that DHAV-3 infection induces 5,396 differentially expressed genes (DEGs), 85 differentially expressed miRNAs (DEMs), and 727 differentially expressed lncRNAs (DELs) at 24 hpi in S vs. R ducks, those upregulated genes were enriched in inflammation and cell communications pathways and downregulated genes were related to metabolic processes. Upregulated genes showed high connectivity with the miR-33, miR-193, and miR-11591, and downregulated genes were mainly regulated by miR-2954, miR-125, and miR-146b. With the construction of lncRNA-miRNA-mRNA axis, we further identified a few aberrantly expressed lncRNAs (e.g., MSTRG.36194.1, MSTRG.50601.1, MSTRG.34328.7, and MSTRG.29445.1) that regulate expression of hub genes (e.g., THBD, CLIC2, IL8, ACOX2, GPHN, SMLR1, and HAO1) by sponging those highly connected miRNAs. Altogether, our findings defined a dual role of ncRNAs in immune and metabolic regulation during DHAV-3 infection, suggesting potential new targets for treating DHAV-3 infected ducks.

MicroRNA-329-3p inhibits the Wnt/β-catenin pathway and proliferation of osteosarcoma cells by targeting transcription factor 7-like 1

An important factor in the emergence and progression of osteosarcoma (OS) is the dysregulated expression of microRNAs (miRNAs). Transcription factor 7-like 1 (TCF7L1), a member of the T cell factor/lymphoid enhancer factor (TCF/LEF) transcription factor family, interacts with the Wnt signaling pathway regulator β-catenin and acts as a DNA-specific binding protein. This study sought to elucidate the impact of the interaction between miR-329-3p and TCF7L1 on the growth and apoptosis of OS and analyze the regulatory expression relationship between miRNA and mRNA in osteosarcoma cells using a variety of approaches. MiR329-3p was significantly downregulated, while TCF7L1 was considerably up-regulated in all examined OS cell lines. Additionally, a clinical comparison study was performed using the TCGA database. Subsequently, the regulatory relationship between miR-329-3p and TCF7L1 on the proliferation and apoptosis of OS cells was verified through in vitro and in vivo experiments. When miR-329-3p was transfected into the OS cell line, the expression of TCF7L1 decreased, the proliferation of OS cells was inhibited, the cytoskeleton disintegrated, and the nucleus condensed to form apoptotic bodies. The expression of proteins that indicate apoptosis increased simultaneously. The cell cycle was arrested in the G0/G1 phase, and the G1/S transition was blocked. The introduction of miR-329-3p also inhibited downstream Cyclin D1 of the Wnt pathway. Xenograft experiments indicated that the overexpression of miR-329-3p significantly inhibited the growth of OS xenografts in nude mice, and the expression of TCF7L1 and c-Myc in tumor tissues decreased. MiR-329-3p was significantly reduced in OS cells and played a suppressive role in tumorigenesis and proliferation by targeting TCF7L1 both in vitro and in vivo. Osteosarcoma cell cycle arrest and pathway inhibition were observed upon the regulation of TCF7L1 by miR-329-3p. Summarizing these results, it can be inferred that miR-329-3p exerts anticancer effects in osteosarcoma by inhibiting TCF7L1.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

miR-212/132 attenuates OVA-induced airway inflammation by inhibiting mast cells activation through MRGPRX2 and ASAP1

Allergic asthma is a chronic inflammatory disease of airways involving complex mechanisms, including MAS-related GPR family member X2 (MRGPRX2) and its orthologue MRGPRB2 on mast cells (MCs). Although miRNAs have been previously shown to related to allergic asthma, the role of miR-212/132 in this process has not been studied. In this study, the predicted pairing of miRNAs and MRGPRX2 (MRGPRB2) mRNAs was carried out by online databases and the function was verify using in vivo and in vitro experiments. Database prediction showed that miR-212/132 interact with MRGPRX2 and MRGPRB2. miR-212/132 mimics alleviated MRGPRB2 mRNA expression as well as pathology changes in lungs and AHR of mice with airway inflammation in vivo. The expression level of MRGPRB2 in the mice lungs after inhaled OVA was also decreased by miR-212/132 mimics. Meanwhile, miR-212/132 inhibited MCs degranulation and cytokines release triggered by C48/80 in vitro. Further, ASAP1 (ARF GTPase-Activating Protein 1) was selected from the junction related pathways using RNAseq and KEGG enrichment. ASAP1 mRNA level was upregulated in airway inflammation and MCs activation and decreased by miR-212/132 mimics. miR-212/132 attenuated OVA-induced airway inflammation by inhibiting MCs activation through MRGPRX2 and ASAP1.

Differential Expression of miRNAs, lncRNAs, and circRNAs between Ovaries and Testes in Common Carp (Cyprinus carpio)

Female common carp grow faster than male individuals, implying that rearing females could be more profitable in aquaculture. Non-coding RNAs (ncRNAs) serve as versatile regulators with multiple functions in diverse biological processes. However, the roles of ncRNAs in the sex differentiation of common carp are less studied. In this study, we investigated the expression profiles of ncRNAs, including miRNAs, lncRNAs, and circRNAs, in the gonads to comprehend the roles of ncRNAs in sex differentiation in common carp. A substantial number of differentially expressed (DE) ncRNAs in ovaries and testes were identified. Some miRNAs, notably miR-205, miR-214, and miR-460-5p, might modulate hormone synthesis and thus maintain sex. A novel miRNA, novel_158, was predicted to suppress the expression of foxl3. DE lncRNAs were associated with oocyte meiosis, GnRH signaling pathways, and steroid biosynthesis, while DE circRNA target genes were enriched in the WNT signaling pathway and MAPK signaling pathway. We also analyzed ncRNA-mRNA interactions to shed light on the crosstalk between competing endogenous RNAs (ceRNAs), which is the critical mechanism by which lncRNAs and circRNAs function. Some lncRNAs and circRNAs may be able to competitively bind novel_313, a new miRNA, and thus regulate hsd17β3. Our research will provide a valuable resource for understanding the genetic basis of gonadal differentiation and development in common carp.

Significant role of some miRNAs as biomarkers for the degree of obesity

BACKGROUND

Obesity is one of the most serious problems over the world. MicroRNAs have developed as main mediators of metabolic processes, playing significant roles in physiological processes. Thus, the present study aimed to evaluate the expressions of (miR-15a, miR-Let7, miR-344, and miR-365) and its relationship with the different classes in obese patients.

METHODS

A total of 125 individuals were enrolled in the study and classified into four groups: healthy non-obese controls (n = 50), obese class I (n = 24), obese class II (n = 17), and obese class III (n = 34) concerning body mass index (BMI < 30 kg/m2, BMI 30-34.9 kg/m2, BMI 35-39.9 kg/m2 and BMI ≥ 40 kg/m2, respectively). BMI and the biochemical measurements (fasting glucose, total cholesterol, triglycerides, HDL and LDL, urea, creatinine, AST, and ALT) were determined. The expressions of (miR-15a, miR-Let7, miR-344, and miR-365) were detected through quantitative real-time PCR (RT-qPCR).

RESULTS

There was a significant difference between different obese classes and controls (P < 0.05) concerning (BMI, TC, TG, HDL, and LDL). In contrast, fasting glucose, kidney, and liver functions had no significant difference. Our data revealed that the expression of miR-15a and miR-365 were significantly associated with different obese classes. But the circulating miR-Let7 and miR-344 were not significantly related to obesity in different classes.

CONCLUSION

Our study indicated that miR-15a and miR-365 might consider as biomarkers for the obesity development into different obese classes. Thus, the relationship between regulatory microRNAs and disease has been the object of intense investigation.

Exploring the Relationship between Fusion Genes and MicroRNAs in Cancer

Fusion genes are key cancer driver genes that can be used as potential drug targets in precision therapies, and they can also serve as accurate diagnostic and prognostic biomarkers. The fusion genes can cause microRNA (miRNA/miR) aberrations in many types of cancer. Nevertheless, whether fusion genes incite miRNA aberrations as one of their many critical oncogenic functionalities for driving carcinogenesis needs further investigation. Recent discoveries of miRNA genes that are present within the regions of genomic rearrangements that initiate fusion gene-based intronic miRNA dysregulation have brought the fusion genes into the limelight and revealed their unexplored potential in the field of cancer biology. Fusion gene-based 'promoter-switch' event aberrantly activate the miRNA-related upstream regulatory signals, while fusion-based coding region alterations disrupt the original miRNA coding loci. Fusion genes can potentially regulate the miRNA aberrations regardless of the protein-coding capability of the resultant fusion transcript. Studies on out-of-frame fusion and nonrecurrent fusion genes that cause miRNA dysregulation have attracted the attention of researchers on fusion genes from an oncological perspective and therefore could have potential implications in cancer therapies. This review will provide insights into the role of fusion genes and miRNAs, and their possible interrelationships in cancer.

The landscape of microRNA interaction annotation: analysis of three rare disorders as a case study

In recent years, a huge amount of data on ncRNA interactions has been described in scientific papers and databases. Although considerable effort has been made to annotate the available knowledge in public repositories, there are still significant discrepancies in how different resources capture and interpret data on ncRNA functional and physical associations. In the present paper, we present a collection of microRNA-mRNA interactions annotated from the scientific literature following recognized standard criteria and focused on microRNAs, which regulate genes associated with rare diseases as a case study. The list of protein-coding genes with a known role in specific rare diseases was retrieved from the Genome England PanelApp, and associated microRNA-mRNA interactions were annotated in the IntAct database and compared with other datasets. RNAcentral identifiers were used for unambiguous, stable identification of ncRNAs. The information about the interaction was enhanced by a detailed description of the cell types and experimental conditions, providing a computer-interpretable summary of the published data, integrated with the huge amount of protein interactions already gathered in the database. Furthermore, for each interaction, the binding sites of the microRNA are precisely mapped on a well-defined mRNA transcript of the target gene. This information is crucial to conceive and design optimal microRNA mimics or inhibitors to interfere in vivo with a deregulated process. As these approaches become more feasible, high-quality, reliable networks of microRNA interactions are needed to help, for instance, in the selection of the best target to be inhibited and to predict potential secondary off-target effects. Database URL https://www.ebi.ac.uk/intact.

Turning off a few overexpressed genes in prostate cancer with microRNAs using a 7mer-seed match model

PURPOSE

This research aims to identify the miRNAs that could target the genes overexpressed in prostate cancer so that miRNA-based therapeutics could be developed.

METHODS

A 7mer-m8 model of microRNA targeting was utilized in order to analyse the relationship between microRNAs and overexpressed genes. The efficiency of miRNA binding was investigated using various parameters namely free energy (AMFE), GC and GC3 content, translation efficiency, cosine similarity metric, mRNA stability, free energy of RNA duplex, and base compositional difference. BLAST2GO software was used to elucidate the functional roles of the genes overexpressed in prostate cancer.

RESULTS

The current research reveals that the coding sequences of the genes were found targeted with multiple miRNAs. For instance, the HPN gene was targeted by the microRNA miR-4279 at two distinct sites i.e. 263-278 and 746-761 in the coding sequence. In the present study, it was observed that the target region of the genes exhibited a comparatively high GC and GC3 contents in comparison to the flanking regions. A low translational rate and weak relationship between RSCU and tRNA were obtained which may be due to the absence of optimal codons.

CONCLUSION

In this study, we have uncovered the human miRNAs that have potential for binding to the coding sequences of 14 most overexpressed genes in prostate cancer and thereby could silence those genes.

Alternative splicing impacts microRNA regulation within coding regions

MicroRNAs (miRNAs) are small non-coding RNA molecules that bind to target sites in different gene regions and regulate post-transcriptional gene expression. Approximately 95% of human multi-exon genes can be spliced alternatively, which enables the production of functionally diverse transcripts and proteins from a single gene. Through alternative splicing, transcripts might lose the exon with the miRNA target site and become unresponsive to miRNA regulation. To check this hypothesis, we studied the role of miRNA target sites in both coding and non-coding regions using six cancer data sets from The Cancer Genome Atlas (TCGA) and Parkinson's disease data from PPMI. First, we predicted miRNA target sites on mRNAs from their sequence using TarPmiR. To check whether alternative splicing interferes with this regulation, we trained linear regression models to predict miRNA expression from transcript expression. Using nested models, we compared the predictive power of transcripts with miRNA target sites in the coding regions to that of transcripts without target sites. Models containing transcripts with target sites perform significantly better. We conclude that alternative splicing does interfere with miRNA regulation by skipping exons with miRNA target sites within the coding region.

Bta-miR-484 Targets SFRP1 and Affects Preadipocytes Proliferation, Differentiation, and Apoptosis

MicroRNAs (miRNAs) are essential regulators of numerous biological processes in animals, including adipogenesis. Despite the abundance of miRNAs associated with adipogenesis, their exact mechanisms of action remain largely unknown. Our study highlights the role of bta-miR-484 as a major regulator of adipocyte proliferation, apoptosis, and differentiation. Here, we demonstrated that the expression of bta-miR-484 initially increased during adipogenesis before decreasing. Overexpression of bta-miR-484 in adipocytes ultimately inhibited cell proliferation and differentiation, reduced the number of EdU fluorescence-stained cells, increased the number of G1 phase cells, reduced the number of G2 and S phase cells, and downregulated the expression of proliferation markers (CDK2 and PCNA) and differentiation markers (CEBPA, FABP4, and LPL). Additionally, overexpression of bta-miR-484 promoted the expression of apoptosis-related genes (Caspase 3, Caspase 9, and BAX), and increased the number of apoptotic cells observed via flow cytometry. In contrast, bta-miR-484 inhibition in adipocytes yielded opposite effects to those observed during bta-miR-484 overexpression. Moreover, luciferase reporter assays confirmed SFRP1 as a target gene of bta-miR-484, and revealed that bta-miR-484 downregulates SFRP1 mRNA expression. These findings offer compelling evidence that bta-miR-484 targets SFRP1, inhibits proliferation and differentiation, and promotes apoptosis. Therefore, these results offer novel insights into the bta-miR-484 regulation of adipocyte growth and development.

Anti-glioma effect of ginseng-derived exosomes-like nanoparticles by active blood-brain-barrier penetration and tumor microenvironment modulation

Inhibition of tumor growth and normalization of immune responses in the tumor microenvironment (TME) are critical issues for improving cancer therapy. However, in the treatment of glioma, effective nanomedicine has limited access to the brain because of the blood-brain barrier (BBB). Previously, we demonstrated nano-sized ginseng-derived exosome-like nanoparticles (GENs) consisting of phospholipids including various bioactive components, and evaluated anti-tumor immune responses in T cells and Tregs to inhibit tumor progression. It was found that the enhanced targeting ability of GENs to the BBB and glioma induced a significant therapeutic effect and exhibited strong efficacy in recruiting M1 macrophage expression in the TME. GENs were demonstrated to be successful candidates in glioma therapeutics both in vitro and in vivo, suggesting excellent potential for inhibiting glioma progression and regulating tumor-associated macrophages (TAMs).

Noncoding RNA circuitry in melanoma onset, plasticity, and therapeutic response

Melanoma, the cancer of the melanocyte, is the deadliest form of skin cancer with an aggressive nature, propensity to metastasize and tendency to resist therapeutic intervention. Studies have identified that the re-emergence of developmental pathways in melanoma contributes to melanoma onset, plasticity, and therapeutic response. Notably, it is well known that noncoding RNAs play a critical role in the development and stress response of tissues. In this review, we focus on the noncoding RNAs, including microRNAs, long non-coding RNAs, circular RNAs, and other small RNAs, for their functions in developmental mechanisms and plasticity, which drive onset, progression, therapeutic response and resistance in melanoma. Going forward, elucidation of noncoding RNA-mediated mechanisms may provide insights that accelerate development of novel melanoma therapies.

MicroRNA-574-5p targeting HOXC6 expression inhibits the hepatocyte lipid uptake to alleviate non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is one of the main causes of liver disease that has reached its last stage. Over the past few years, evidence for miRNAs' centrality in NAFLD pathogenesis has accumulated. According to some studies, miR-574-5p plays a role in lipid metabolism. However, research on the relationship between miR-574-5p and NAFLD is lacking. For in vivo experiments, we induced the NAFLD mice model with a high-fat diet (HFD). AgomiR-574-5p was injected intravenously into HFD-fed mice for eight weeks, and qPCR was used to identify the expression of miR-574-5p in the serum. In in vitro experiments, The treatment of L-O2 cells with a miR-574-5p mimic resulted in a significant reduction in lipid deposition, suggesting that miR-574-5p can inhibit lipid accumulation and lipid formation induced by OA. The dual-luciferase reporter gene assay revealed that miR-574-5p targets the 3' UTR region of HOXC6 directly. We discovered that OA-induced lipid accumulation in hepatocytes might be mediated through the miR-574-5p-HOXC6 signaling axis. Additional research is required in order to determine the specific mechanism by which HOXC6 downstream pathways are involved in the miR-574-5p induced lipid uptake.

Computational approaches for circRNAs prediction and in silico characterization

Circular RNAs (circRNAs) are single-stranded and covalently closed non-coding RNA molecules originated from RNA splicing. Their functions include regulatory potential over other RNA species, such as microRNAs, messenger RNAs and RNA binding proteins. For circRNA identification, several algorithms are available and can be classified in two major types: pseudo-reference-based and split-alignment-based approaches. In general, the data generated from circRNA transcriptome initiatives is deposited on public specific databases, which provide a large amount of information on different species and functional annotations. In this review, we describe the main computational resources for the identification and characterization of circRNAs, covering the algorithms and predictive tools to evaluate its potential role in a particular transcriptomics project, including the public repositories containing relevant data and information for circRNAs, recapitulating their characteristics, reliability and amount of data reported.

New Insight into Mechanisms of Cardiovascular Diseases: An Integrative Analysis Approach to Identify TheranoMiRNAs

MiRNAs regulate both physiological and pathological heart functions. Altered expression of miRNAs is associated with cardiovascular diseases (CVDs), making miRNAs attractive therapeutic strategies for the diagnosis and treatment of heart diseases. A recent publication defined, for the first time, the term theranoMiRNA, meaning the miRNAs that may be used both for diagnosis and treatment. The use of in silico tools may be considered fundamental for these purposes, clarifying several molecular aspects, suggesting future directions for in vivo studies. This study aims to explore different bioinformatic tools in order to clarify miRNA interactions with candidate genes, demonstrating the need to use a computational approach when establishing the most probable associations between miRNAs and target genes. This study focused on the functions of miR-133a-3p, miR-21-5p, miR-499a-5p, miR-1-3p, and miR-126-3p, providing an up-to-date overview, and suggests future lines of research in the identification of theranoMiRNAs related to CVDs. Based on the results of the present study, we elucidated the molecular mechanisms that could be linked between miRNAs and CVDs, confirming that these miRNAs play an active role in the genesis and development of heart damage. Given that CVDs are the leading cause of death in the world, the identification of theranoMiRNAs is crucial, hence the need for a definition of in vivo studies in order to obtain further evidence in this challenging field of research.

The non-coding genome in Autism Spectrum Disorders

Autism Spectrum Disorders (ASD) are a group of neurodevelopmental disorders (NDDs) characterized by difficulties in social interaction and communication, repetitive behavior, and restricted interests. While ASD have been proven to have a strong genetic component, current research largely focuses on coding regions of the genome. However, non-coding DNA, which makes up for ∼99% of the human genome, has recently been recognized as an important contributor to the high heritability of ASD, and novel sequencing technologies have been a milestone in opening up new directions for the study of the gene regulatory networks embedded within the non-coding regions. Here, we summarize current progress on the contribution of non-coding alterations to the pathogenesis of ASD and provide an overview of existing methods allowing for the study of their functional relevance, discussing potential ways of unraveling ASD's "missing heritability".

Integration of miRNA's Theranostic Potential with Nanotechnology: Promises and Challenges for Parkinson's Disease Therapeutics

Despite the wide research going on in Parkinson's disease (PD), the burden of PD still remains high and continues to increase. The current drugs available for the treatment of PD are only aimed at symptomatic control. Hence, research is mainly focused on identifying the novel therapeutic targets that can be effectively targeted in order to slow down or culminate the disease progression. Recently the role of microRNAs (miRNAs) in the regulation of various pathological mechanisms of PD has been thoroughly explored and many of them were found to be dysregulated in the biological samples of PD patients. These miRNAs can be used as diagnostic markers and novel therapeutic options to manage PD. The delivery of miRNAs to the target site in brain is a challenging job owing to their nature of degradability by endonucleases as well as poor blood brain barrier (BBB) permeability. Nanoparticles appear to be the best solution to effectively encase the miRNA in their core as well as cross the BBB to deliver them into brain. Functionalisation of these nanoparticles further enhances the site-specific delivery.

Erythrocyte miRNA-92a-3p interactions with PfEMP1 as determinants of clinical malaria

Based on the recently added high throughput analysis data on small noncoding RNAs in modulating disease pathophysiology of malaria, we performed an integrative computational analysis for exploring the role of human-host erythrocytic microRNAs (miRNAs) and their influence on parasite survival and host homeostasis. An in silico analysis was performed on transcriptomic datasets accessed from PlasmoDB and Gene Expression Omnibus (GEO) repositories analyzed using miRanda, miRTarBase, mirDIP, and miRDB to identify the candidate miRNAs that were further subjected to network analysis using MCODE and DAVID. This was followed by immune infiltration analysis and screening for RNA degradation mechanisms. Seven erythrocytic miRNAs, miR-451a, miR-92a-3p, miR-16-5p, miR-142-3p, miR-15b-5p, miR-19b-3p, and miR-223-3p showed favourable interactions with parasite genes expressed during blood stage infection. The miR-92a-3p that targeted the virulence gene PfEMP1 showed drastic reduction during infection. Performing pathway analysis for the human-host gene targets for the miRNA identified TOB1, TOB2, CNOT4, and XRN1 genes that are associated to RNA degradation processes, with the exoribonuclease XRN1, highly enriched in the malarial samples. On evaluating the role of exoribonucleases in miRNA degradation further, the pattern of Plasmodium falciparum_XRN1 showed increased levels during infection thus suggesting a defensive role for parasite survival. This study identifies miR-92a-3p, a member of C13orf25/ miR-17-92 cluster, as a novel miRNA inhibitor of the crucial parasite genes responsible for symptomatic malaria. Evidence for a plausible link to chromosome 13q31.3 loci controlling the epigenetic disease regulation is also suggested.

Insights into Online microRNA Bioinformatics Tools

MicroRNAs (miRNAs) are members of the small non-coding RNA family regulating gene expression at the post-transcriptional level. MiRNAs have been found to have critical roles in various biological and pathological processes. Research in this field has significantly progressed, with increased recognition of the importance of miRNA regulation. As a result of the vast data and information available regarding miRNAs, numerous online tools have emerged to address various biological questions related to their function and influence across essential cellular processes. This review includes a brief introduction to available resources for an investigation covering aspects such as miRNA sequences, target prediction/validation, miRNAs associated with disease, pathway analysis and genetic variants within miRNAs.

Association of OPRM1, MIR23B, and MIR107 genetic variability with acute pain, chronic pain and adverse effects after postoperative tramadol and paracetamol treatment in breast cancer

BACKGROUND

Tramadol is an opioid analgesic often used for pain management after breast cancer surgery. Its analgesic activity is due to the activation of the μ-opioid receptor, encoded by the OPRM1 gene. This study investigated the association of genetic variability in OPRM1 and its regulatory miRNA genes with outcomes of tramadol/paracetamol treatment after breast cancer surgery with axillary lymphadenectomy.

PATIENTS AND METHODS

The study included 113 breast cancer patients after breast cancer surgery with axillary lymphadenectomy treated with either 75/650 mg or 37.5/325 mg of tramadol with paracetamol for pain relief within the randomized clinical trial KCT 04/2015-DORETAonko/si at the Institute of Oncology Ljubljana. All patients were genotyped for OPRM1 rs1799971 and rs677830, MIR23B rs1011784, and MIR107 rs2296616 using competitive allele-specific PCR. The association of genetic factors with acute and chronic pain as well as adverse effects of tramadol treatment was evaluated using logistic regression, Fisher's exact test, and Mann-Whitney test.

RESULTS

The investigated OPRM1 related polymorphisms were not associated with acute pain assessed with the VAS scale within four weeks after surgery (all P > 0.05). Carriers of at least one polymorphic OPRM1 rs1799971 allele had a higher risk of constipation in the first four weeks after surgery compared to non-carriers (OR = 4.5, 95% CI = 1.6-12.64, P = 0.004). Carriers of at least one polymorphic OPRM1 rs677830 allele had a higher risk of constipation after third week of tramadol treatment (OR = 3.11, 95% CI = 1.08-8.89, P = 0.035). Furthermore, carriers of two polymorphic MIR23B rs1011784 alleles had a higher risk of nausea after 28 days of tramadol treatment (OR = 7.35, 95% CI = 1.27-42.6, P = 0.026), while heterozygotes for MIR107 rs2296616 allele had a lower risk of nausea after 21 days of tramadol treatment (OR = 0.21, 95% CI = 0.05-0.87, P = 0.031). In carriers of two polymorphic MIR107 rs2296616 alleles, chronic pain was significantly more common than in carriers of two wild-type alleles (P = 0.004). Carriers of at least one polymorphic MIR23B rs1011784 allele experienced more neuropathic pain after adjustment for tramadol dose (OR = 2.85, 95% CI = 1.07-7.59, P = 0.036), while carriers of at least one polymorphic OPRM1 rs677830 allele experienced less neuropathic pain compared to carriers of two wild-type alleles (OR = 0.38, 95% CI = 0.15-0.99, P = 0.047).

CONCLUSIONS

Genetic variability of OPRM1 and genes coding for miRNAs that could affect OPRM1 expression may be associated with adverse effects of tramadol/paracetamol treatment as well as with chronic and neuropathic pain after breast cancer surgery with axillary lymphadenectomy.

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

RNA-seq analysis reveals the critical role of the novel lncRNA BIANCR in intramuscular adipogenesis through the ERK1/2 signaling pathway

BACKGROUND

Long non-coding RNAs (lncRNAs) regulate numerous biological processes, including adipogenesis. Research on adipogenesis will assist in the treatment of human metabolic diseases and improve meat quality in livestock, such as the content of intramuscular fat (IMF). However, the significance of lncRNAs in intramuscular adipogenesis remains unclear. This research aimed to reveal the lncRNAs transcriptomic profiles in the process of bovine intramuscular adipogenesis and to identify the lncRNAs involved in the adipogenesis of bovine intramuscular adipocytes.

RESULTS

In this research, a landscape of lncRNAs was identified with RNA-seq in bovine intramuscular adipocytes at four adipogenesis stages (0 d, 3 d, 6 d, and 9 d after differentiation). A total of 7035 lncRNAs were detected, including 3396 novel lncRNAs. Based on the results of differential analysis, co-expression analysis, and functional prediction, we focused on the bovine intramuscular adipogenesis-associated long non-coding RNA (BIANCR), a novel lncRNA that may have an important regulatory function. The knockdown of BIANCR inhibited proliferation and promoted apoptosis of intramuscular preadipocytes. Moreover, BIANCR knockdown inhibited intramuscular adipogenesis by regulating the ERK1/2 signaling pathway.

CONCLUSION

This study obtained the landscape of lncRNAs during adipogenesis in bovine intramuscular adipocytes. BIANCR plays a crucial role in adipogenesis through the ERK1/2 signaling pathway. The results are noteworthy for improving beef meat quality, molecular breeding, and metabolic disease research.

Molecular Mediators and Regulators of Retinal Angiogenesis

BACKGROUND

Retinal neovascularization is the major cause of vision loss that affects both adults and young children including premature babies. It has been a major pathology in several retinal diseases like age-related macular degeneration (AMD), diabetic retinopathy (DR) and retinopathy of prematurity (ROP). Current treatment modalities such as anti-VEGF therapy, laser are not suitable for every patient and response to these therapies is highly variable. Thus, there is a need to investigate newer therapeutic targets for DR, ROP and AMD, based on a clear understanding of disease pathology and regulatory mechanisms involved.

METHOD

Appropriate articles published till February 2021 were extracted from PUBMED using keywords like ocular angiogenesis, DR, ROP, AMD, miRNA, mRNA, and cirMiRNA and containvaluable information regarding the involvement of miRNA in causing neovascularization. After compiling the list of miRNA regulating mRNA expression in angiogenesis and neovascularaization, their interactions were studied using online available tool MIENTURNET (http://userver.bio.uniroma1.it/apps/mienturnet/). The pathways involved in these processes were also predicted using the same tool.

RESULTS

Most of the studies have explored potential targets like HIF1-α, PDGF, TGFβ, FGF, etc., for their involvement in pathological angiogenesis in different retinal diseases. The regulatory role of microRNA (miRNA) has also been explored in various retinal ocular pathologies. This review highlights regulatory mechanism of cellular and circulatory miRNAs and their interactions with the genes involved in retinal neovascularization. The role of long noncoding RNA (ncRNA) in the regulation of genes involved in different pathways is also noteworthy and discussed in this review.

CONCLUSION

This review highlights the potential regulatory mechanism/pathways involved in retinal neovascularization and its implications in retinal diseases and for identifying new drug targets.

HDAC Inhibitor SAHA Alleviates Pyroptosis by up-regulating miR-340 to Inhibit NEK7 Signaling in Subarachnoid Hemorrhage

BACKGROUND

Subarachnoid hemorrhage (SAH) is a cerebral hemorrhagic disease with a high disability and fatality rate. Cell pyroptosis is involved in the brain injury following SAH. Here, we explored the effect of HDAC inhibitor SAHA against cell pyroptosis after SAH.

METHODS

The rat SAH model was established by endovascular perforation and the rat microglia were treated with 25 μm oxyhemoglobin (OxyHb) for 24 h to mimic SAH model in vitro. Neurological score and brain edema were assessed in rat SAH model. TUNEL staining detected apoptosis. qRT-PCR and western blotting were employed to detect expression levels of miR-340, NEK7 and inflammatory cytokines. ELISA assay determined the secretion of IL-1β and IL-18 in rat serum and cell supernatant. A lactate dehydrogenase (LDH) kit measured the LDH activity in rat primary microglia. Microglia pyroptosis was detected by flow cytometry. RIP and dual luciferase reporter assays confirmed the binding relationship between miR-340 and NEK7.

RESULTS

SAHA alleviated neurological dysfunction, inflammatory injury and microglia pyroptosis in SAH rats. SAHA suppressed LDH release, inflammatory factor expression and pyroptosis in microglia treated with OxyHb. Meanwhile, SAHA increased miR-340 expression and inhibited NEK7 level in vivo and in vitro SAH models. Further, miR-340 directly targeted NEK7 to inhibit the NLRP3 signaling pathway. Knockdown of miR-340 or overexpression of NEK7 reversed the suppressive effects of SAHA on microglia inflammation and pyroptosis. Additionally, knockdown of NEK7 impaired microglia inflammation and pyroptosis induced by miR-340 inhibitor.

CONCLUSION

HDAC inhibitor SAHA ameliorates microglia pyroptosis in SAH through triggering miR-340 expression to suppress NEK7 signaling. This novel mechanism provides promise for SAHA in SAH treatment.

Tumour Necrosis Factor-Alpha (TNF-α)-Induced Metastatic Phenotype in Colorectal Cancer Epithelial Cells: Mechanistic Support for the Role of MicroRNA-21

Colorectal cancer is driven by genetic and epigenetic changes in cells to confer phenotypes that promote metastatic transformation and development. Tumour necrosis factor-alpha (TNF-α), a pro-inflammatory mediator, regulates cellular communication within the tumour microenvironment and is associated with the progression of the metastatic phenotype. Oncogenic miR-21 has been shown to be overexpressed in most solid tumours, including colorectal cancer, and is known to target proteins involved in metastatic transformation. In this study, we investigated the relationship between TNF-α and miR-21 regulation in colorectal cancer epithelial cells (SW480 and HCT116). We observed that TNF-α, at concentrations reported to be present in serum and tumour tissue from colorectal cancer patients, upregulated miR-21 expression in both cell lines. TNF-α treatment also promoted cell migration, downregulation of the expression of E-cadherin, a marker of epithelial to mesenchymal transition, and anti-apoptotic BCL-2 (a validated target for miR-21). Knockdown of miR-21 had the opposite effect on each of these TNF-a induced phenotypic changes. Additionally, in the SW480 cell line, although TNF-α treatment selectively induced expression of a marker of metastatic progression VEGF-A, it failed to affect MMP2 expression or invasion activity. Our data indicate that exposing colorectal cancer epithelial cells to TNF-α, at concentrations occurring in the serum and tumour microenvironment of colorectal cancer patients, upregulated miR-21 expression and promoted the metastatic phenotype.

RNA silencing proteins and small RNAs in oomycete plant pathogens and biocontrol agents

Introduction

Oomycetes cause several damaging diseases of plants and animals, and some species also act as biocontrol agents on insects, fungi, and other oomycetes. RNA silencing is increasingly being shown to play a role in the pathogenicity of Phytophthora species, either through trans-boundary movement of small RNAs (sRNAs) or through expression regulation of infection promoting effectors.

Methods

To gain a wider understanding of RNA silencing in oomycete species with more diverse hosts, we mined genome assemblies for Dicer-like (DCL), Argonaute (AGO), and RNA dependent RNA polymerase (RDRP) proteins from Phytophthora plurivora, Ph. cactorum, Ph. colocasiae, Pythium oligandrum, Py. periplocum, and Lagenidium giganteum. Moreover, we sequenced small RNAs from the mycelium stage in each of these species.

Results and discussion

Each of the species possessed a single DCL protein, but they differed in the number and sequence of AGOs and RDRPs. SRNAs of 21nt, 25nt, and 26nt were prevalent in all oomycetes analyzed, but the relative abundance and 5' base preference of these classes differed markedly between genera. Most sRNAs mapped to transposons and other repeats, signifying that the major role for RNA silencing in oomycetes is to limit the expansion of these elements. We also found that sRNAs may act to regulate the expression of duplicated genes. Other sRNAs mapped to several gene families, and this number was higher in Pythium spp., suggesting a role of RNA silencing in regulating gene expression. Genes for most effector classes were the source of sRNAs of variable size, but some gene families showed a preference for specific classes of sRNAs, such as 25/26 nt sRNAs targeting RxLR effector genes in Phytophthora species. Novel miRNA-like RNAs (milRNAs) were discovered in all species, and two were predicted to target transcripts for RxLR effectors in Ph. plurivora and Ph. cactorum, indicating a putative role in regulating infection. Moreover, milRNAs from the biocontrol Pythium species had matches in the predicted transcriptome of Phytophthora infestans and Botrytis cinerea, and L. giganteum milRNAs matched candidate genes in the mosquito Aedes aegypti. This suggests that trans-boundary RNA silencing may have a role in the biocontrol action of these oomycetes.

Innovative approaches in transforming microRNAs into therapeutic tools

MicroRNA (miRNA) is regarded as a prominent genetic regulator, as it can fine-tune an entire biological pathway by targeting multiple target genes. This characteristic makes miRNAs promising therapeutic tools to reinstate cell functions that are disrupted as a consequence of diseases. Currently, miRNA replacement by miRNA mimics and miRNA inhibition by anti-miRNA oligonucleotides are the main approaches to utilizing miRNA molecules for therapeutic purposes. Nevertheless, miRNA-based therapeutics are hampered by major issues such as off-target effects, immunogenicity, and uncertain delivery platforms. Over the past few decades, several innovative approaches have been established to minimize off-target effects, reduce immunostimulation, and provide efficient transfer to the target cells in which these molecules exert their function. Recent achievements have led to the testing of miRNA-based drugs in clinical trials, and these molecules may become next-generation therapeutics for medical intervention. Despite the achievement of exciting milestones, the dosage of miRNA administration remains unclear, and ways to address this issue are proposed. Elucidating the current status of the main factors of therapeutic miRNA would allow further developments and innovations to achieve safe therapeutic tools. This article is categorized under: RNA in Disease and Development > RNA in Disease Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action.

MicroRNA-Target Identification: A Combinatorial In Silico Approach

Contemporary computational target prediction tools with their distinctive properties and stringency have been playing a vital role in pursuing putative targets for a solitary miRNA or a subcategory of miRNAs. These tools utilize a defined set of probabilistic algorithms, machine learning techniques, and information of experimentally validated miRNA targets to provide the best selection. However, there are numerous false-positive predictions, and a method to choose an archetypal approach and put the data provided by the prediction tools into context is still lacking. Moreover, sensitivity, specificity, and overall efficiency of a single tool have not yet been achieved. Therefore, a systematic combination of selective online tools combining elementary and advanced factors of miRNA target identification might reinforce the current target prediction regime. The focus of this study was to build a comprehensive workflow by combining six available online tools to facilitate the current understanding of miRNA-target prediction.

Resistance to white spot syndrome virus in the European shore crab is associated with suppressed virion trafficking and heightened immune responses

Introduction

All decapod crustaceans are considered potentially susceptible to White Spot Syndrome Virus (WSSV) infection, but the degree of White Spot Disease (WSD) susceptibility varies widely between species. The European shore crab Carcinus maenas can be infected with the virus for long periods of time without signs of disease. Given the high mortality rate of susceptible species, the differential susceptibility of these resistant hosts offers an opportunity to investigate mechanisms of disease resistance.

Methods

Here, the temporal transcriptional responses (mRNA and miRNA) of C. maenas following WSSV injection were analysed and compared to a previously published dataset for the highly WSSV susceptible Penaeus vannamei to identify key genes, processes and pathways contributing to increased WSD resistance.

Results

We show that, in contrast to P. vannamei, the transcriptional response during the first 2 days following WSSV injection in C. maenas is limited. During the later time points (7 days onwards), two groups of crabs were identified, a recalcitrant group where no replication of the virus occurred, and a group where significant viral replication occurred, with the transcriptional profiles of the latter group resembling those of WSSV-susceptible species. We identify key differences in the molecular responses of these groups to WSSV injection.

Discussion

We propose that increased WSD resistance in C. maenas may result from impaired WSSV endocytosis due to the inhibition of internal vesicle budding by dynamin-1, and a delay in movement to the nucleus caused by the downregulation of cytoskeletal transcripts required for WSSV cytoskeleton docking, during early stages of the infection. This response allows resistant hosts greater time to fine-tune immune responses associated with miRNA expression, apoptosis and the melanisation cascade to defend against, and clear, invading WSSV. These findings suggest that the initial stages of infection are key to resistance to WSSV in the crab and highlight possible pathways that could be targeted in farmed crustacean to enhance resistance to WSD.

Hsa_circ_0015278 Regulates FLT3-ITD AML Progression via Ferroptosis-Related Genes

AML with the FLT3-ITD mutation seriously threatens human health. The mechanism by which circRNAs regulate the pathogenesis of FLT3-ITD mutant-type AML through ferroptosis-related genes (FerRGs) remains unclear. Differentially expressed circRNAs and mRNAs were identified from multiple integrated data sources. The target miRNAs and mRNAs of the circRNAs were predicted using various databases. The PPI network, ceRNA regulatory network, GO, and KEGG enrichment analyses were performed. The "survival" and the "pROC" R packages were used for K-M and ROC analysis, respectively. GSEA, immune infiltration analysis, and clinical subgroup analysis were performed. Finally, circRNAs were validated by Sanger sequencing and qRT-PCR. In our study, 77 DECircs-1 and 690 DECircs-2 were identified. Subsequently, 11 co-up-regulated DECircs were obtained by intersecting DECircs-1 and DECircs-2. The target miRNAs of the circRNAs were screened by CircInteractome, circbank, and circAtlas. Utilizing TargetScan, ENCORI, and miRWalk, the target mRNAs of the miRNAs were uncovered. Ultimately, 73 FerRGs were obtained, and the ceRNA regulatory network was constructed. Furthermore, MAPK3 and CD44 were significantly associated with prognosis. qRT-PCR results confirmed that has_circ_0015278 was significantly overexpressed in FLT3-ITD mutant-type AML. In summary, we constructed the hsa_circ_0015278/miRNAs/FerRGs signaling axis, which provides new insight into the pathogenesis and therapeutic targets of AML with FLT3-ITD mutation.

Identification of serum MiRNAs as candidate biomarkers for non-small cell lung cancer diagnosis

BACKGROUND

Lung cancer is one of the most common solid tumors worldwide and the leading cause of cancer-associated death. Non-small cell lung cancer (NSCLC) is accounts for approximately 85% of all the lung cancers and lung squamous carcinoma (SCC) and adenocarcinoma (ADC) are the main subtypes of NSCLC. Early diagnose using serum biomarkers could improve the overall survival of patients. In this study, we aimed to identify miRNAs from serum with clinical utility in the diagnosis of NSCLC.

METHODS

Ten patients with SCC, ten patients with ADC and five noncancerous individuals were enrolled in the screening cohort. miRNA expression levels in serum were measured by microarray analysis. Candidate miRNAs were validated by real-time quantitative polymerase chain reaction analysis in a validation cohort of 78 NSCLC patients and 44 noncancerous individuals. Receiver operating characteristic curves were used to assess the diagnostic performance of serum miRNAs for NSCLC. Logistic regression was used to evaluate the diagnostic value of the combination of markers.

RESULTS

Six candidate miRNAs were differentially expressed between NSCLC patients and noncancerous individuals in the screening set (fold change > 2, p < 0.05). Among them, expression levels of miR-3149 and miR-4769.3p were confirmed to be significantly increased in tumor serum in the validation set. The area under the curve values of miR-3149 and miR-4769.3p in distinguishing NSCLC patients from noncancerous controls were 0.830 and 0.735, respectively. When combined with tumor markers CEA and Cyfra21-1, the joint diagnostic model increased the area under the curve to 0.898.

CONCLUSION

Serum miRNAs miR-3149 and miR-4769.3p were up-regulated in NSCLC and may be potential biomarkers for early diagnosis of lung cancer.

Tumor Suppressor miRNA-503 Inhibits Cell Invasion in Head and Neck Cancer through the Wnt Signaling Pathway via the WNT3A/MMP Molecular Axis

Head and neck cancer (HNC) is the fifth most common cancer worldwide, and its incidence and death rates have been consistently high throughout the past decades. MicroRNAs (miRNAs) have recently gained significant attention because of their role in the regulation of a variety of biological processes via post-transcriptional silencing mechanisms. Previously, we determined a specific profile of miRNAs associated with HNC using a miRNA microarray analysis. Of the 23 miRNAs with highly altered expression in HNC cells, miR-503 was the most significantly downregulated miRNA. In this study, we confirmed that miR-503 acts as a tumor suppressor, as our results showed decreased levels of miR-503 in cancer cells and patients with HNC. We further characterized the role of miR-503 in the malignant functions of HNC. Although there was a minimal effect on cell growth, miR-503 was found to inhibit cellular invasion significantly. Algorithm-based studies identified multiple potential target genes and pathways associated with oncogenic mechanisms. The candidate target gene, WNT3A, was confirmed to be downregulated by miR-503 at both the mRNA and protein levels and validated by a reporter assay. Furthermore, miR-503 modulated multiple invasion-associated genes, including matrix metalloproteinases (MMPs), through the Wnt downstream signaling pathway. Overall, this study demonstrates that miR-503 suppresses HNC malignancy by inhibiting cell invasion through the Wnt signaling pathway via the WNT3A/MMP molecular axis. The modulation of miR-503 may be a novel therapeutic approach to intervene in cancer invasion.

5’isomiR-183-5p|+2 elicits tumor suppressor activity in a negative feedback loop with E2F1

Background

MicroRNAs (miRNAs) and isomiRs play important roles in tumorigenesis as essential regulators of gene expression. 5’isomiRs exhibit a shifted seed sequence compared to the canonical miRNA, resulting in different target spectra and thereby extending the phenotypic impact of the respective common pre-miRNA. However, for most miRNAs, expression and function of 5’isomiRs have not been studied in detail yet. Therefore, this study aims to investigate the functions of miRNAs and their 5’isomiRs.

Methods

The expression of 5’isomiRs was assessed in The Cancer Genome Atlas (TCGA) breast cancer patient dataset. Phenotypic effects of miR-183 overexpression in triple-negative breast cancer (TNBC) cell lines were investigated in vitro and in vivo by quantifying migration, proliferation, tumor growth and metastasis. Direct targeting of E2F1 by miR-183-5p|+2 was validated with a 3’UTR luciferase assay and linked to the phenotypes of isomiR overexpression.

Results

TCGA breast cancer patient data indicated that three variants of miR-183-5p are highly expressed and upregulated, namely miR-183-5p|0, miR-183-5p|+1 and miR-183-5p|+2. However, TNBC cell lines displayed reduced proliferation and invasion upon overexpression of pre-miR-183. While invasion was reduced individually by all three isomiRs, proliferation and cell cycle progression were specifically inhibited by overexpression of miR-183-5p|+2. Proteomic analysis revealed reduced expression of E2F target genes upon overexpression of this isomiR, which could be attributed to direct targeting of E2F1, specifically by miR-183-5p|+2. Knockdown of E2F1 partially phenocopied the effect of miR-183-5p|+2 overexpression on cell proliferation and cell cycle. Gene set enrichment analysis of TCGA and METABRIC patient data indicated that the activity of E2F strongly correlated with the expression of miR-183-5p, suggesting transcriptional regulation of the miRNA by a factor of the E2F family. Indeed, in vitro, expression of miR-183-5p was regulated by E2F1. Hence, miR-183-5p|+2 directly targeting E2F1 appears to be part of a negative feedback loop potentially fine-tuning its activity.

Conclusions

This study demonstrates that 5’isomiRs originating from the same arm of the same pre-miRNA (i.e. pre-miR-183-5p) may exhibit different functions and thereby collectively contribute to the same phenotype. Here, one of three isomiRs was shown to counteract expression of the pre-miRNA by negatively regulating a transcriptional activator (i.e. E2F1). We speculate that this might be part of a regulatory mechanism to prevent uncontrolled cell proliferation, which is disabled during cancer progression.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02380-8.

Microarray Analysis Reveals Changes in tRNA-Derived Small RNAs (tsRNAs) Expression in Mice with Septic Cardiomyopathy

Background: tRNA-derived small RNAs (tsRNAs) as a novel non-coding RNA have been studied in many cardiovascular diseases, but the relationship between tsRNAs and septic cardiomyopathy has not been investigated. We sought to analyze changes of the expression profile of tsRNAs in septic cardiomyopathy and reveal an important role for tsRNAs. Methods: We constructed a sepsis model by cecal ligation and puncture (CLP) in mice, and microarray analysis was used to find differentially expressed tsRNAs. Quantitative real-time PCR was used to verify the expression of tsRNAs and the interference effect of angiogenin (ANG), a key nuclease producing tsRNAs. Bioinformatics analysis was used to predict target genes and functions. CCK-8 and LDH release assays were used to detect cell viability and cell death. Results: A total of 158 tsRNAs were screened, of which 101 were up-regulated and 57 were down-regulated. A total of 8 tsRNAs were verified by qPCR, which was consistent with microarray results. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses suggest that these tsRNAs may be associated with the Wnt signaling pathway and participate in cellular process. The expression of tsRNAs decreased after the interference of the key nuclease ANG, while CCK-8 suggested a corresponding decrease in cell viability and an increase in the release of LDH (cell death), indicating that tsRNAs can protect cardiomyocytes during the development of septic cardiomyopathy, reduced cardiomyocyte death. Conclusions: A total of 158 tsRNAs changed significantly in septic cardiomyopathy, and these tsRNAs may play a protective role in the development of septic cardiomyopathy.

miRNome-transcriptome analysis unveils the key regulatory pathways involved in the tumorigenesis of tongue squamous cell carcinoma

Tongue squamous cell carcinoma (TSCC) is considered the most common malignant tumor among the oral squamous cell carcinomas with a poor prognosis. Understanding the underlying molecular mechanisms that underpin TSCC and its treatments is the focus of the research. Deregulated expression of microRNAs (miRNAs) has recently been implicated in various biological processes linked to cancer. Therefore, in this study, we attempted to investigate miRNAs and their targets expressed in TSCC, which could be involved in its oncogenesis. We performed next-generation sequencing of small RNAs and transcriptomes in H357 TSCC cell line and human oral keratinocytes as a control to find miRNAs and mRNAs that are differentially expressed (DE), which were then supplemented with additional expression datasets from databases, yielding 269 DE miRNAs and 2094 DE genes. The target prediction followed by pathway and disease function analysis revealed that the DE targets were significantly associated with the key processes and pathways, such as apoptosis, epithelial-mesenchymal transition, endocytosis and vascular endothelial growth factor signaling pathways. Furthermore, the top 12 DE targets were chosen based on their involvement in more than one cancer-related pathway, of which 6 genes are targeted by miR-128-3p. Real-time quantitative PCR validation of this miRNA and its targets in H357 and SCC9 TSCC cells confirmed their possible targeting from their reciprocal expression, with MAP2K7 being a critical target that might be involved in oncogenesis and progression of TSCC by acting as a tumor suppressor. Further research is underway to understand how miR-128-3p regulates oncogenesis in TSCC via MAP2K7 and associated pathways.

A Computationally Constructed lncRNA-Associated Competing Triplet Network in Clear Cell Renal Cell Carcinoma

Long noncoding RNAs (lncRNAs) are revealed to be involved in the tumorigenesis and progression of human malignancies mediated by microRNA (miRNA) via the competing endogenous RNA (ceRNA) mechanism, a newly proposed “RNA language.” However, the lncRNA-associated competing triplet (lncACT) network among ceRNA transcripts in clear cell renal cell carcinoma (ccRCC) is currently lacking. We carried out differential expression analysis to identify aberrantly expressed lncRNAs, miRNAs, and mRNAs by analyzing the RNA-seq data of 420 ccRCC tissues and 71 noncancerous kidney tissues obtained from The Cancer Genome Atlas (TCGA). Then, a ccRCC-specific ceRNA network was built using computational algorithms, including miRcode, TargetScan, miRanda, and miRTarBase. In total, 1491 dysregulated lncRNAs were found between normal renal tissues and ccRCC (fold $\text{change}>4$ and false discovery $\text{rate}<0.01$ ). A ceRNA network that comprised of 46 DElncRNAs, 11 DEmiRNAs, and 55 DEmRNAs was established by integrating the lncRNA/miRNA and miRNA/mRNA interactions into lncACTs. Several lncRNAs were identified to be significantly associated with clinical features of ccRCC patients. Notably, four key lncRNAs (TCL6, HOTTIP, HULC, and PCGEM1) were tightly correlated with both patients’ clinical characteristics and overall survival (log-rank $P<0.05$ ), indicating their potential important roles in ccRCC. HOTTIP may be a potential prognostic and therapeutic molecular marker for ccRCC patients. Collectively, our results provide a comprehensive view of the lncRNA-associated ceRNA regulatory network for a better understanding of the mechanisms and prognosis biomarkers for ccRCC.

Differentially expressed tRNA-derived fragments in bovine fetuses with assisted reproduction induced congenital overgrowth syndrome

Background: As couples struggle with infertility and livestock producers wish to rapidly improve genetic merit in their herd, assisted reproductive technologies (ART) have become increasingly popular in human medicine as well as the livestock industry. Utilizing ART can cause an increased risk of congenital overgrowth syndromes, such as Large Offspring Syndrome (LOS) in ruminants. A dysregulation of transcripts has been observed in bovine fetuses with LOS, which is suggested to be a cause of the phenotype. Our recent study identified variations in tRNA expression in LOS individuals, leading us to hypothesize that variations in tRNA expression can influence the availability of their processed regulatory products, tRNA-derived fragments (tRFs). Due to their resemblance in size to microRNAs, studies suggest that tRFs target mRNA transcripts and regulate gene expression. Thus, we have sequenced small RNA isolated from skeletal muscle and liver of day 105 bovine fetuses to elucidate the mechanisms contributing to LOS. Moreover, we have utilized our previously generated tRNA sequencing data to analyze the contribution of tRNA availability to tRF abundance. Results: 22,289 and 7,737 unique tRFs were predicted in the liver and muscle tissue respectively. The greatest number of reads originated from 5' tRFs in muscle and 5' halves in liver. In addition, mitochondrial (MT) and nuclear derived tRF expression was tissue-specific with most MT-tRFs and nuclear tRFs derived from LysUUU and iMetCAU in muscle, and AsnGUU and GlyGCC in liver. Despite variation in tRF abundance within treatment groups, we identified differentially expressed (DE) tRFs across Control-AI, ART-Normal, and ART-LOS groups with the most DE tRFs between ART-Normal and ART-LOS groups. Many DE tRFs target transcripts enriched in pathways related to growth and development in the muscle and tumor development in the liver. Finally, we found positive correlation coefficients between tRNA availability and tRF expression in muscle (R = 0.47) and liver (0.6). Conclusion: Our results highlight the dysregulation of tRF expression and its regulatory roles in LOS. These tRFs were found to target both imprinted and non-imprinted genes in muscle as well as genes linked to tumor development in the liver. Furthermore, we found that tRNA transcription is a highly modulated event that plays a part in the biogenesis of tRFs. This study is the first to investigate the relationship between tRNA and tRF expression in combination with ART-induced LOS.

miR-146a and miR-200b alter cognition by targeting NMDA receptor subunits

MicroRNAs fine-tune gene regulation and can be targeted for therapeutic purposes. We investigated the physiological roles of miR-146a and miR-200b that are differentially expressed in neurological disorders such as Alzheimer's disease and schizophrenia, particularly in learning and memory mechanisms. Using bioinformatics tools and luciferase assay, we show interaction of these miRNAs with transcripts of N-methyl-D-aspartate receptor (NMDAR) subunits Grin2A and Grin2B. Overexpression of these miRNAs in primary hippocampal neurons caused downregulation of GluN2B and GluN2A proteins. Stereotactic injections of these miRNAs into rat hippocampus caused cognitive deficits in multiple behavioral tests with decreased protein levels of GluN1, GluN2A, GluN2B, AMPAR subunit GluR1, and Neuregulin 1. In pharmacologically treated rat models [MK-801 treated and methylazoxymethanol acetate (MAM) treated], we found upregulated levels of these miRNAs, implying their involvement in downregulating NMDAR subunits in these models. These results suggest the importance of miR-146a-5p and miR-200b-3p in hippocampus-dependent learning and memory.

Role of epigallocatechin-3- gallate in the regulation of known and novel microRNAs in breast carcinoma cells

Breast cancer comprises 30% of all cancer cases among the world’s women population. MicroRNAs are small, endogenous, non-coding RNAs that regulate cell proliferating and apoptotic pathways by modulating expressions of related genes. Phytochemicals like epigallocatechin-3-gallate (EGCG) are known to have a chemotherapeutic effect on cancer often through the regulation of microRNAs. The aim is to find out the key known and novel miRNAs, which are controlled by EGCG in breast cancer cell line MDA-MB-231. Next-generation sequencing (NGS) revealed 1,258 known and 330 novel miRNAs from untreated and 83 μM EGCG (IC50 value of EGCG) treated cells. EGCG modulated 873 known and 47 novel miRNAs in the control vs. treated sample. The hypothesis of EGCG being a great modulator of miRNAs that significantly control important cancer-causing pathways has been established by analyzing with Kyoto Encyclopedia of Genes and Genomes (KEGG) and Protein Analysis Through Evolutionary Relationships (PANTHER) database. Validation of known and novel miRNA expression differences in untreated vs. treated cells was done using qPCR. From this study, a few notable miRNAs were distinguished that can be used as diagnostics as well as prognostic markers for breast cancer.

Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.