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

Liquid Biopsy-Based Exo-oncomiRNAs Can Predict Prostate Cancer Aggressiveness

Liquid biopsy-based biomarkers, including microRNAs packaged within extracellular vesicles, are promising tools for patient management. The cytokine tumor necrosis factor-like weak inducer of apoptosis (TWEAK) is related to PCa progression and is found in the semen of patients with PCa. TWEAK can induce the transfer of exo-oncomiRNAs from tumor cells to body fluids, and this process might have utility in non-invasive PCa prognosis. We investigated TWEAK-regulated exo-microRNAs in semen and in post-digital rectal examination urine from patients with different degrees of PCa aggressiveness. We first identified 14 exo-oncomiRNAs regulated by TWEAK in PCa cells in vitro, and subsequently validated those using liquid biopsies from 97 patients with PCa. Exo-oncomiR-221-3p, -222-3p and -31-5p were significantly higher in the semen of high-risk patients than in low-risk peers, whereas exo-oncomiR-193-3p and -423-5p were significantly lower in paired samples of post-digital rectal examination urine. A panel of semen biomarkers comprising exo-oncomiR-221-3p, -222-3p and TWEAK was designed that could correctly classify 87.5% of patients with aggressive PCa, with 85.7% specificity and 76.9% sensitivity with an area under the curve of 0.857. We additionally found that TWEAK modulated two exo-oncomiR-221-3p targets, TCF12 and NLK. Overall, we show that liquid biopsy detection of TWEAK-regulated exo-oncomiRNAs can improve PCa prognosis prediction.

MiR-338-3p improved lung adenocarcinoma by AKAP12 suppression

INTRODUCTION

This study aimed to explore the biological functions of AKAP12 in lung adenocarcinoma and investigate the interaction between AKAP12 and miR-338-3p.

MATERIAL AND METHODS

Sixty-one differentially expressed genes in lung adenocarcinoma and adjacent normal tissues were first analyzed by TCGA. Immunohistochemistry and quantitative reverse transcription PCR (qRT-PCR) were further utilized to confirm aberrant AKAP12 expression in tumor tissues. The influences of AKAP12 on proliferation, invasion and migration, and apoptosis of lung adenocarcinoma were investigated by clone formation assay and MTT assay, transwell assay, and flow cytometry analysis respectively. TargetScan and miRanda databases predicted the binding sites of miRNAs on AKAP12 3'-UTR and structure changes were validated by RNA folding form. The target relationship between miR-338-3p and AKAP12 was confirmed by the dual-luciferase reporter system. Disease-free survival (DFS) and overall survival (OS) curves were generated with Kaplan-Meier plotter according to the TCGA data and the correlation among AKAP12 expression, miR-338-3p expression and prognosis was also analyzed.

RESULTS

AKAP12 was upregulated in lung adenocarcinoma tissues and cells (all p < 0.01), and negatively correlated with prognosis outcomes of patients (both p < 0.05). High expression of AKAP12 promoted proliferation, invasion and migration of cancer cells, and inhibited cell apoptosis (all p < 0.05). MiR-338-3p could directly bind to the 3'-UTR of AKAP12 and showed most significant suppression on AKAP12 expression among four predicted miRNAs (all p < 0.01). Additionally, miR-338-3p could suppress AKAP12 in lung adenocarcinoma, improving prognosis (all p < 0.05).

CONCLUSIONS

AKAP12 acts as a tumor promoter in lung adenocarcinoma development. Upregulation of MiR-338-3p could suppress AKAP12 expression in lung cancer cells and contribute to a better prognosis.

Analysis of microRNA Regulation in Single Cells

MicroRNAs (miRNAs) regulate gene expression by binding to mRNAs. Consequently, they reduce target gene expression levels and expression variability, also known as "noise." Single-cell RNA sequencing (scRNA-seq) technology has been used to study miRNA and mRNA expression in single cells, and has demonstrated its strength in quantifying cell-to-cell variation. Here we describe how to investigate miRNA regulation using data with both mRNA and miRNA expression in single cell format. We show that miRNAs reduce the expression levels and also expression noise of target genes in single cells. Finally, we also discuss potential improvements in experimental design and computational analysis of scRNA-seq in order to reduce or partition the technical noise.

MicroRNA and mRNA profiling of cerebral cortex in a transgenic mouse model of Alzheimer's disease by RNA sequencing

In a previous study, we found that long non-coding genes in Alzheimer’s disease (AD) are a result of endogenous gene disorders caused by the recruitment of microRNA (miRNA) and mRNA, and that miR-200a-3p and other representative miRNAs can mediate cognitive impairment and thus serve as new biomarkers for AD. In this study, we investigated the abnormal expression of miRNA and mRNA and the pathogenesis of AD at the epigenetic level. To this aim, we performed RNA sequencing and an integrative analysis of the cerebral cortex of the widely used amyloid precursor protein and presenilin-1 double transgenic mouse model of AD. Overall, 129 mRNAs and 68 miRNAs were aberrantly expressed. Among these, eight down-regulated miRNAs and seven up-regulated miRNAs appeared as promising noninvasive biomarkers and therapeutic targets. The main enriched signaling pathways involved mitogen-activated kinase protein, phosphatidylinositol 3-kinase-protein kinase B, mechanistic target of rapamycin kinase, forkhead box O, and autophagy. An miRNA-mRNA network between dysregulated miRNAs and corresponding target genes connected with AD progression was also constructed. These miRNAs and mRNAs are potential biomarkers and therapeutic targets for new treatment strategies, early diagnosis, and prevention of AD. The present results provide a novel perspective on the role of miRNAs and mRNAs in AD. This study was approved by the Experimental Animal Care and Use Committee of Institute of Medicinal Biotechnology of Beijing, China (approval No. IMB-201909-D6) on September 6, 2019.

Transcription Factors ZEB1 and CREB Promote the Transcription of Bovine ABHD5 Gene

Alpha/beta hydrolase domain 5 (ABHD5) plays a significant role in intracellular lipid metabolism, which is regulated by a complex network of transcription factors. The transcriptional regulation of the ABHD5 gene in cattle and other livestock, however, has not been previously investigated. Investigations in humans and animal models indicate that the transcription factors zinc finger E-box binding homeobox 1 (ZEB1) and cAMP-response element binding protein (CREB) may play important roles in the transcriptional regulation of ABHD5 in cattle. Our comparison of the sequence similarities in the transcription factor binding sites in Bos taurus, Bos indicus, Bos mutus, and Homo sapiens revealed high homology. Based on the data collected by the Cistrome Data Browser and its visualization window, we found that ZEB1 and CREB have significant ChIP-seq enrichments in the 5'-untranslated region (5' UTR) of the human ABHD5 gene. In bovine adipocytes, we detected ZEB1 and CREB binding sites in the ABHD5 gene. Mutations in the ZEB1 and CREB binding sites significantly reduced the promoter activity (p < 0.05 and p < 0.01, respectively). Moreover, electrophoretic mobility shift assays and chromatin immunoprecipitation (ChIP) assays demonstrated the binding of the transcription factors in vivo and in vitro, respectively. And overexpression or silencing the expression of the ZEB1 and CREB, respectively, resulted in significant changes to the ABHD5 promoter activity. Collectively, these results indicate that ZEB1 and CREB are important transcription factors that regulate ABHD5 gene expression in bovine adipocytes. They further our understanding of the transcriptional regulation and biological functions of the bovine ABHD5 gene.

Up-regulated spinal microRNAs induce aggregation of superoxide dismutase 1 protein in canine degenerative myelopathy

Canine degenerative myelopathy (DM) is a fatal progressive neurodegenerative disease. Mutations in the superoxide dismutase 1 (SOD1) gene have been shown to be the major risk factor for DM, and it is hypothesized that neural degeneration is caused by a "gain of toxic function" of mutant SOD1. In this study, the spinal cord microRNA (miRNA) profiles of DM-affected dogs were investigated to elucidate the pathomechanisms of DM. Quantification of 277 miRNAs identified three up-regulated miRNAs and 18 down-regulated miRNAs in the spinal cords of DM-affected dogs. Based on gene ontology analysis, the target cluster of up-regulated miRNAs was associated with protein expression or modification and cellular response, and that of down-regulated miRNAs was associated with tissue development. In these clusters, we focused on the mechanism of protein ubiquitination. Polyubiquitination assay demonstrated that canine SOD1 proteins were polyubiquitinated and degraded by proteasomes. Immunohistochemistry of the spinal cords of DM-affected dogs showed that mutant SOD1 aggregations were not ubiquitin immunopositive. Using cultured cells, co-transfection of canine SOD1 and up-regulated miRNA in DM-affected dogs demonstrated that miR-23a, miR-142 and miR-221 significantly increased the proportion of cells with mutant SOD1 aggregation. These results suggested that up-regulated miRNAs in the spinal cords of DM-affected dogs may inhibit ubiquitination of misfolded SOD1 protein and induce mutant SOD1 aggregations, leading to further progression of degenerative processes in the DM pathology.

Integrative analysis of transcriptomes highlights potential functions of transfer-RNA-derived small RNAs in experimental intracerebral hemorrhage

Transfer-RNA-derived small RNAs (tsRNAs) are a novel class of short non-coding RNAs, that possess regulatory functions. However, their biological roles in hemorrhagic stroke are not understood. In this study, by RNA sequencing, we investigated the tsRNA expression profiles of intracerebral hemorrhagic rat brains in the chronic phase. A total of 331 tsRNAs were identified (308 in sham and 309 in intracerebral hemorrhage). Among them, the validation revealed that 7 tsRNAs (1 up-regulated and 6 down-regulated) were significantly changed. Subsequently, we predicted the target mRNAs of the 7 tsRNAs. Through integrative analysis, the predicted targets were validated by mRNA microarray data. Moreover, we confirmed the functions of tsRNAs targeting mRNAs in vitro. Furthermore, using bioinformatics tools and databases, we developed a tsRNA-mRNA-pathway interaction network to visualize their potential functions. Bioinformatics analyses and confirmatory experiments indicated that the altered genes were mainly enriched in several signaling pathways. These pathways were interrelated with intracerebral hemorrhage, such as response to oxidative stress, endocytosis, and regulation of G protein-coupled receptor signaling pathway. In summary, this study systematically revealed the profiles of tsRNAs after an experimental intracerebral hemorrhage. These results may provide novel therapeutic targets following a hemorrhagic stroke in the chronic phase.

Circ_0075829 facilitates the progression of pancreatic carcinoma by sponging miR-1287-5p and activating LAMTOR3 signalling

Pancreatic cancer (PC) is a leading cause of cancer‐related mortality globally. Though increasing evidence has demonstrated that circular RNAs (circRNAs) are linked to the development and progression of cancers, the biological functions of circRNAs in PC remain largely unexplored so far. Based on previous studies, Hsc_circ_0075829 (circ_0075829) was screened out and then further identified in PC clinical specimens and cell lines by real‐time PCR. After the stability tests, a series of in vitro and in vivo functional experiments were performed to investigate the role of circ_0075829 in PC development. Furthermore, fluorescent in situ hybridization (FISH), bioinformatics tools, dual‐luciferase assays and rescue experiments were conducted to clarify the regulatory mechanisms of circ_0075829 in SW1990 and BxPC‐3 cells. Compared with paracancerous tissues, the expression of circ_0075829 was increased in PC tissues, which was positively correlated with the clinical features of PC. Knockdown of circ_0075829 significantly suppressed the proliferative, migratory and invasive rates of SW1990 and BxPC‐3 cells both in vitro and in vivo. Bioinformatics analysis and dual‐luciferase reporter gene assay indicated that circ_0075829 could bind to miR‐1287‐5p. Mechanism research and rescue experiments demonstrated that circ_0075829 could regulate the LAMTOR3/p‐ERK signalling pathway via sponging miR‐1287‐5p in PC cell lines. Our data reveal that the circ_0075829 could facilitate the proliferation and metastasis of PC through circ_0075829/miR‐1287‐5p/LAMTOR3 axis.

The Role of Circular RNAs in Pancreatic Ductal Adenocarcinoma and Biliary-Tract Cancers

Simple Summary

Pancreatic and biliary tract cancers often present with non-specific symptoms, resulting in diagnosis at a late stage. This may be too late for curative surgery. Earlier detection and characterisation may guide treatment options and increase survival. Natural “circles” of RNA (circRNAs) are shown to regulate cancer-related genes, and act as cancer “biomarkers”. Recent research has shown that circRNAs are both abundant and stable, both of which are desirable characteristics for clinically useful biomarkers. In this systematic review, we describe the roles of circRNAs in pancreatic and biliary tract cancers, summarise the current published research and explore their utility as a biomarker. A total of 32 articles were included: 22 considering Pancreatic Cancer, 7 for Bile Duct Cancer and 3 for Gallbladder Cancer. CircRNA proved an exciting prospect as a biomarker for these cancers and future work should continue to develop and expand this field of research.

Abstract

Pancreatic Ductal Adenocarcinoma (PDAC) and biliary-tract cancers (BTC) often present at a late stage, and consequently patients have poor survival-outcomes. Circular RNAs (circRNAs) are non-coding RNA molecules whose role in tumourigenesis has recently been realised. They are stable, conserved and abundant, with tissue-specific expression profiles. Therefore, significant interest has arisen in their use as potential biomarkers for PDAC and BTC. High-throughput methods and more advanced bioinformatic techniques have enabled better profiling and progressed our understanding of how circRNAs may function in the competing endogenous RNA (ceRNA) network to influence the transcriptome in these cancers. Therefore, the aim of this systematic review was to describe the roles of circRNAs in PDAC and BTC, their potential as biomarkers, and their function in the wider ceRNA network in regulating microRNAs and the transcriptome. Medline, Embase, Scopus and PubMed were systematically reviewed to identify all the studies addressing circRNAs in PDAC and BTC. A total of 32 articles were included: 22 considering PDAC, 7 for Cholangiocarcinoma (CCA) and 3 for Gallbladder Cancer (GBC). There were no studies investigating Ampullary Cancer. Dysregulated circRNA expression was associated with features of malignancy in vitro, in vivo, and ex vivo. Overall, there have been very few PDAC and BTC tissues profiled for circRNA signatures. Therefore, whilst the current studies have demonstrated some of their functions in these cancers, further work is required to elucidate their potential role as cancer biomarkers in tissue, biofluids and biopsies.

miR-421 promotes the viability of A549 lung cancer cells by targeting forkhead box O1

MicroRNA (miR)-421 has been reported to serve various important roles in numerous types of cancer, including neuroblastoma and gastric cancer. However, to the best of our knowledge, few reports have determined the role of miR-421 in lung cancer. The aim of the current study was to analyze the expression levels of miR-421 in A549 lung cancer cells, to determine the target gene of miR-421, and to investigate the function and mechanism of miR-421 in cellular cytotoxicity. miR-421 expression levels were analyzed in A549 lung cancer cells using reverse transcription-quantitative PCR, a MTT assay was performed to determine the effect of miR-421 on A549 cell cytotoxicity and the protein expression levels of forkhead box O1 (FOXO1) were determined via western blotting. The target gene of miR-421 was predicted and verified using TargetScan and a dual-luciferase reporter assay, respectively. The results revealed that miR-421 expression levels were significantly upregulated in A549 lung cancer cell lines compared with the normal cells (P<0.01). Additionally, it was discovered that miR-421 promoted A549 cell viability (P<0.01) compared with A549 transfected with negative control. miR-421 was also identified to bind to the 3′-untranslated region of FOXO1. In A549 cells transfected with miR-421-mimics, the expression levels of phosphorylated (p)-AKT, p-glycogen synthase kinase-3β, p-retinoblastoma and cyclin D1 were significantly upregulated (P<0.01), whereas the expression levels of FOXO1 and p21 were significantly downregulated (P<0.01) compared with the control group. In conclusion, the results of the present study suggested that miR-421 may promote the viability of A549 lung cancer cells by targeting FOXO1 and modulating cell cycle, indicating that targeting miR-421 and FOXO1 may represent future therapeutic strategies for the treatment of patients with lung cancer.

TMEM63C, a Potential Novel Target for Albuminuria Development, Is Regulated by MicroRNA-564 and Transforming Growth Factor beta in Human Renal Cells

INTRODUCTION

Transmembrane protein (TMEM) 63C is a member of the TMEM gene family and was recently linked to glomerular filtration barrier function and albuminuria. Its molecular function and expression regulation are largely unknown.

OBJECTIVE

In this study, we set out to characterize the regulating impact of microRNAs (miRNAs) such as miRNA-564 (miR-564) on TMEM63C expression in renal cells. Also, we examined the influence of transforming growth factor beta (TGF-ß) on TMEM63C expression and the potential impact of TMEM63C inhibition on epithelial-mesenchymal transition (EMT) in renal cells and on cell viability in human embryonic kidney 293 cells (HEK 293).

METHODS

Expression analyses were done using real-time PCR and Western blot. Dual luciferase assay was performed to determine the miRNA-mediated expression control. Cell viability was assessed via trypan blue exclusion staining.

RESULTS AND CONCLUSIONS

MiR-564 reduced TMEM63C expression in HEK 293 and human podocytes (hPC). The treatment of renal cells with TGF-ß led to an increased expression of TMEM63C. Moreover, a reduced TMEM63C expression was associated with a changed ratio of EMT marker proteins such as α-smooth muscle actin versus E-cadherin in HEK 293 and decreased nephrin expression in hPC. In addition, cell viability was reduced upon inhibition of TMEM63C expression in HEK 293. This study demonstrates first mechanisms involved in TMEM63C expression regulation and a link to EMT in renal cells.

Transboundary Pathogenic microRNA Analysis Framework for Crop Fungi Driven by Biological Big Data and Artificial Intelligence Model

Plant fungal diseases have been affecting the world's agricultural production and economic levels for a long time, such as rice blast, gray tomato mold, potato late blight etc. Recent studies have shown that fungal pathogens transmit microRNA as an effector to host plants for infection. However, bioassay-based verification analysis is time-consuming and challenging, and it is difficult to analyze from a global perspective. With the accumulation of fungal and plant-related data, data analysis methods can be used to analyze pathogenic fungal microRNA further. Based on the microRNA expression data of fungal pathogens infecting plants before and after, this paper discusses the selection strategy of sample data, the extraction strategy of pathogenic fungal microRNA, the prediction strategy of a fungal pathogenic microRNA target gene, the bicluster-based fungal pathogenic microRNA functional analysis strategy and experimental verification methods. A general analysis pipeline based on machine learning and bicluster-based function module was proposed for plant-fungal pathogenic microRNA.The pipeline proposed in this paper is applied to the infection process of Magnaporthe oryzae and the infection process of potato late blight. It has been verified to prove the feasibility of the pipeline. It can be extended to other relevant crop pathogen research, providing a new idea for fungal research on plant diseases. It can be used as a reference for understanding the interaction between fungi and plants.

Ischemic dilated cardiomyopathy pathophysiology through microRNA-16-5p

INTRODUCTION AND OBJECTIVES

The expression levels of microRNA-16-5p (miR-16) are upregulated in ischemic cardiomyopathy and in animal models of ischemic dilated cardiomyopathy (iDCM), inducing myocardial apoptosis. We investigated the role of miR-16 in the adaptive cellular response associated with endoplasmic reticulum (ER) stress and autophagy in the apoptotic iDCM environment.

METHODS

We quantified the miR-16 plasma levels of 168 participants-76 controls, 60 iDCM patients, and 32 familial DCM patients with the pathogenic variant of BAG3-by quantitative real-time polymerase chain reaction and correlated the levels with patient variables. The effects of intracellular miR-16 overexpression were analyzed in a human cardiac cell line. Apoptosis and cell viability were measured, as well as the levels of markers associated with ER stress, cardiac injury, and autophagy.

RESULTS

Plasma miR-16 levels were upregulated in iDCM patients (P=.039). A multivariate logistic regression model determined the association of miR-16 with iDCM clinical variables (P <.001). In vitro, miR-16 overexpression increased apoptosis (P=.02) and reduced cell viability (P=.008). Furthermore, it induced proapoptotic components of ER stress, based on upregulation of the PERK/CHOP pathway. However, we observed augmentation of autophagic flux (P <.001) without lysosomal blockade by miR-16 as a possible cytoprotective mechanism.

CONCLUSIONS

MiR-16 is specifically associated with iDCM. In an ischemic setting, miR-16 activates ER stress and promotes inflammation followed by autophagy in human cardiac cells. Thus, autophagy may be an attempt to maintain cellular homeostasis in response to misfolded/aggregated proteins related to ER stress, prior to apoptosis.

Circulating microRNAs in Breast Milk and Their Potential Impact on the Infant

MicroRNAs (MiRNAs) are small RNA molecules that can exert regulatory functions in gene expression. MiRNAs have been identified in diverse tissues and biological fluids, both in the context of health and disease. Breastfeeding has been widely recognized for its superior nutritional benefits; however, a number of bioactive compounds have been found to transcend these well-documented nutritional contributions. Breast milk was identified as a rich source of miRNAs. There has been increasing interest about their potential ability to transfer to the offspring as well as what their specific involvement is within the benefits of breast milk in the infant. In comparison to breast milk, formula milk lacks many of the benefits of breastfeeding, which is thought to be a result of the absence of some of these bioactive compounds. In recent years, the miRNA profile of breast milk has been widely studied, along with the possible transfer mechanisms throughout the infant's digestive tract and the role of miRNA-modulated genes and their potential protective and regulatory functions. Nonetheless, to date, the current evidence is not consistent, as many methodological limitations have been identified; hence, discrepancies exits about the biological functions of miRNAs. Further research is needed to provide thorough knowledge in this field.

The Hsa_circ_0091579/miR-940/TACR1 Axis Regulates the Development of Hepatocellular Carcinoma

Purpose

Circular RNAs (circRNAs) play important roles in hepatocellular carcinoma (HCC) development. The circRNA hsa_circ_0091579 (circ_0091579) is dysregulated in HCC, while the mechanism of circ_0091579 in HCC development is largely unknown.

Patients and Methods

Thirty paired cancer and adjacent normal tissues were harvested from HCC patients. SNU-387 and Huh7 cells were cultured in this study. circ_0091579, microRNA-940 (miR-940) and tachykinin-1 receptor (TACR1) abundances were measured via quantitative reverse transcription-polymerase chain reaction or Western blot. Cell viability, migration, invasion, colony ability, cell cycle distribution and apoptosis were assessed via 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide, transwell assay, colony formation assay and flow cytometry. The interaction among circ_0091579, miR-940 and TACR1 was tested via dual-luciferase reporter analysis. The anti-HCC role of circ_0091579 knockdown in vivo was investigated using xenograft model.

Results

circ_0091579 expression was enhanced in HCC tissue samples and cells. circ_0091579 silence inhibited cell viability, migration, invasion and colony formation, induced cell cycle arrest at G0/G1 phase, and promoted apoptosis in HCC cells. miR-940 was targeted via circ_0091579 and miR-940 knockdown reversed the suppressive effect of circ_0091579 silence on HCC development. miR-940 targeted TACR1 to repress HCC development. circ_0091579 could regulate TACR1 expression by mediating miR-940. Down-regulation of circ_0091579 decreased xenograft tumor growth.

Conclusion

Knockdown of circ_0091579 repressed HCC development by mediating miR-940/TACR1 axis, indicating a new pathogenesis of HCC.

Robust partitioning of microRNA targets from downstream regulatory changes

The biological impact of microRNAs (miRNAs) is determined by their targets, and robustly identifying direct miRNA targets remains challenging. Existing methods suffer from high false-positive rates and are unable to effectively differentiate direct miRNA targets from downstream regulatory changes. Here, we present an experimental and computational framework to deconvolute post-transcriptional and transcriptional changes using a combination of RNA-seq and PRO-seq. This novel approach allows us to systematically profile the regulatory impact of a miRNA. We refer to this approach as CARP: Combined Analysis of RNA-seq and PRO-seq. We apply CARP to multiple miRNAs and show that it robustly distinguishes direct targets from downstream changes, while greatly reducing false positives. We validate our approach using Argonaute eCLIP-seq and ribosome profiling, demonstrating that CARP defines a comprehensive repertoire of targets. Using this approach, we identify miRNA-specific activity of target sites within the open reading frame. Additionally, we show that CARP facilitates the dissection of complex changes in gene regulatory networks triggered by miRNAs and identification of transcription factors that mediate downstream regulatory changes. Given the robustness of the approach, CARP would be particularly suitable for dissecting miRNA regulatory networks in vivo.

Popular Computational Tools Used for miRNA Prediction and Their Future Development Prospects

MicroRNAs (miRNAs) are 19-24 nucleotide (nt)-long noncoding, single-stranded RNA molecules that play significant roles in regulating the gene expression, growth, and development of plants and animals. From the year that miRNAs were first discovered until the beginning of the twenty-first century, researchers used experimental methods such as cloning and sequencing to identify new miRNAs and their roles in the posttranscriptional regulation of protein synthesis. Later, in the early 2000s, informatics approaches to the discovery of new miRNAs began to be implemented. With increasing knowledge about miRNA, more efficient algorithms have been developed for computational miRNA prediction. The miRNA research community, hoping for greater coverage and faster results, has shifted from cumbersome and expensive traditional experimental approaches to computational approaches. These computational methods started with homology-based comparisons of known miRNAs with orthologs in the genomes of other species; this method could identify a known miRNA in new species. Second-generation sequencing and next-generation sequencing of mRNA at different developmental stages and in specific tissues, in combination with a better search and alignment algorithm, have accelerated the process of predicting novel miRNAs in a particular species. Using the accumulated annotated miRNA sequence information, researchers have been able to design ab initio algorithms for miRNA prediction independent of genome sequence knowledge. Here, the methods recently used for miRNA computational prediction are summarized and classified into the following four categories: homology-based, target-based, scoring-based, and machine-learning-based approaches. Finally, the future developmental directions of miRNA prediction methods are discussed.

Computational Methods and Software Tools for Functional Analysis of miRNA Data

miRNAs are important regulators of gene expression that play a key role in many biological processes. High-throughput techniques allow researchers to discover and characterize large sets of miRNAs, and enrichment analysis tools are becoming increasingly important in decoding which miRNAs are implicated in biological processes. Enrichment analysis of miRNA targets is the standard technique for functional analysis, but this approach carries limitations and bias; alternatives are currently being proposed, based on direct and curated annotations. In this review, we describe the two workflows of miRNAs enrichment analysis, based on target gene or miRNA annotations, highlighting statistical tests, software tools, up-to-date databases, and functional annotations resources in the study of metazoan miRNAs.

Dietary impact of a plant-derived microRNA on the gut microbiome

Background:

Global estimations of 4 billion people living on plant-based diets signify tremendous diversity in plant consumption and their assorted miRNAs, which presents a challenging model to experimentally address how plant-based miRNAs impact the microbiome. Here we establish baseline gut microbiome composition for a mouse model deficient in the specific mammalian miR-146a shown to alter gut microbiomes. We then asses the effect on the gut microbiome when miR-146a-deficient mice are fed a transgenic plant-based diet expressing the murine-derived miR-146a. Mice deficient in miR-146a were maintained either on a baseline diet until 7 weeks of age (day 0) and then fed either vector or miR-146a-expressing plant-based diets for 21 days. The gut microbiomes of mice were examined by comparing the V4 region of 16S rRNA gene sequences of DNA isolated from fecal samples at days 0 (baseline diet) and 21 (vector or miR-146a expressing plant-based diets).

Results:

Beta-diversity analysis demonstrated that the transition from baseline chow to a plant-based diet resulted in significant longitudinal shifts in microbial community structure attributable to increased fiber intake. Bipartite network analysis suggests that miR-146a-deficient mice fed a plant diet rich in miR-146a have a microbiome population modestly different than mice fed an isogenic control plant diet deficient in miR-146a.

Conclusion:

A mouse diet composed of a transgenic plant expressing a mouse miR-146a may fine tune microbial communities but does not appear to have global effects on microbiome structure and composition.

Key MicroRNA's and Their Targetome in Adrenocortical Cancer

Adrenocortical Carcinoma (ACC) is a rare but aggressive malignancy with poor prognosis and limited response to available systemic therapies. Although complete surgical resection gives the best chance for long-term survival, ACC has a two-year recurrence rate of 50%, which poses a therapeutic challenge. High throughput analyses focused on characterizing the molecular signature of ACC have revealed specific micro-RNAs (miRNAs) that are associated with aggressive tumor phenotypes. MiRNAs are small non-coding RNA molecules that regulate gene expression by inhibiting mRNA translation or degrading mRNA transcripts and have been generally implicated in carcinogenesis. This review summarizes the current insights into dysregulated miRNAs in ACC tumorigenesis, their known functions, and specific targetomes. In addition, we explore the possibility of particular miRNAs to be exploited as clinical biomarkers in ACC and as potential therapeutics.

Individualized multi-omic pathway deviation scores using multiple factor analysis

Malignant progression of normal tissue is typically driven by complex networks of somatic changes, including genetic mutations, copy number aberrations, epigenetic changes, and transcriptional reprogramming. To delineate aberrant multi-omic tumor features that correlate with clinical outcomes, we present a novel pathway-centric tool based on the multiple factor analysis framework called padma. Using a multi-omic consensus representation, padma quantifies and characterizes individualized pathway-specific multi-omic deviations and their underlying drivers, with respect to the sampled population. We demonstrate the utility of padma to correlate patient outcomes with complex genetic, epigenetic, and transcriptomic perturbations in clinically actionable pathways in breast and lung cancer.

miR-125 family regulates XIRP1 and FIH in response to myocardial infarction

MicroRNAs (miRNAs) are powerful regulators of protein expression. Many play important roles in cardiac development and disease. While several miRNAs and targets have been well characterized, the abundance of miRNAs and the numerous potential targets for each suggest that the vast majority of these interactions have yet to be described. The goal of this study was to characterize miRNA expression in the mouse heart after coronary artery ligation (LIG) and identify novel mRNA targets altered during the initial response to ischemic stress. We performed small RNA sequencing (RNA-Seq) of ischemic heart tissue 1 day and 3 days after ligation and identified 182 differentially expressed miRNAs. We then selected relevant mRNA targets from all potential targets by correlating miRNA and mRNA expression from a corresponding RNA-Seq data set. From this analysis we chose to focus, as proof of principle, on two miRNAs from the miR-125 family, miR-125a and miR-351, and two of their potential mRNA targets, Xin actin-binding repeat-containing protein 1 (XIRP1) and factor inhibiting hypoxia-inducible factor (FIH). We found miR-125a to be less abundant and XIRP1 more abundant after ligation. In contrast, the related murine miRNA miR-351 was substantially upregulated in response to ischemic injury, and FIH expression correspondingly decreased. Luciferase reporter assays confirmed direct interactions between these miRNAs and targets. In summary, we utilized a correlative analysis strategy combining miRNA and mRNA expression data to identify functional miRNA-mRNA relationships in the heart after ligation. These findings provide insight into the response to ischemic injury and suggest future therapeutic targets.

Construction and Analysis of Competing Endogenous RNA Networks for Breast Cancer Based on TCGA Dataset

Background

Long noncoding RNAs (lncRNAs) act as competing endogenous RNAs for microRNAs in cancer metastasis. However, the roles of lncRNA-mediated competing endogenous RNA (ceRNA) networks for breast cancer (BC) are still unclear. Material and Methods. The expression profiles of mRNAs, lncRNAs, and miRNAs with BC were extracted from The Cancer Genome Atlas database. Weighted gene coexpression network analysis was conducted to extract differentially expressed mRNAs (DEmRNAs) that might be core genes. Through miRWalk, TargetScan, and miRDB to predict the target genes, an abnormal lncRNA-miRNA-mRNA ceRNA network with BC was constructed. The survival possibilities of mRNAs, miRNAs, and lncRNAs for patients with BC were determined by Kaplan-Meier survival curves and Oncomine.

Results

We identified 2134 DEmRNAs, 1059 differentially expressed lncRNAs (DElncRNAs), and 86 differentially expressed miRNAs (DEmiRNAs). We then compose a ceRNA network for BC, including 72 DElncRNAs, 8 DEmiRNAs, and 12 DEmRNAs. After verification, 2 lncRNAs (LINC00466, LINC00460), 1 miRNA (Hsa-mir-204), and 5 mRNAs (TGFBR2, CDH2, CHRDL1, FGF2, and CHL1) were meaningful as prognostic biomarkers for BC patients. In the ceRNA network, we found that three axes were present in 10 RNAs related to the prognosis of BC, namely, LINC00466-Hsa-mir-204-TGFBR2, LINC00466-Hsa-mir-204-CDH2, and LINC00466-Hsa-mir-204-CHRDL1.

Conclusion

This study highlighted lncRNA-miRNA-mRNA ceRNA related to the pathogenesis of BC, which might be used for latent diagnostic biomarkers and therapeutic targets for BC.

Stabilization of HIF-1α in Human Retinal Endothelial Cells Modulates Expression of miRNAs and Proangiogenic Growth Factors

Retinal hypoxia is one of the causative factors of diabetic retinopathy and is also one of the triggers of VEGF release. We hypothesized that specific dysregulated miRNAs in diabetic retinopathy could be linked to hypoxia-induced damage in human retinal endothelial cells (HRECs). We investigated in HRECs the effects of chemical (CoCl₂) hypoxia on the expression of HIF-1α, VEGF, PlGF, and of a focused set of miRNAs. We found that miR-20a-5p, miR-20b-5p, miR-27a-3p, miR-27b-3p, miR-206-3p, miR-381-3p correlated also with expression of TGFβ signaling pathway genes in HRECs, challenged with chemical hypoxic stimuli. In conclusion, our data suggest that retinal angiogenesis would be promoted, at least under HIF-1α activation, by upregulation of PlGF and other factors such as miRNAs, VEGFA, and TGFβ1.

Circulating microRNA trafficking and regulation: computational principles and practice

Rapid advances in genomics discovery tools and a growing realization of microRNA's implication in intercellular communication have led to a proliferation of studies of circulating microRNA sorting and regulation across cells and different species. Although sometimes, reaching controversial scientific discoveries and conclusions, these studies have yielded new insights in the functional roles of circulating microRNA and a plethora of analytical methods and tools. Here, we consider this body of work in light of key computational principles underpinning discovery of circulating microRNAs in terms of their sorting and targeting, with the goal of providing practical guidance for applications that is focused on the design and analysis of circulating microRNAs and their context-dependent regulation. We survey a broad range of informatics methods and tools that are available to the researcher, discuss their key features, applications and various unsolved problems and close this review with prospects and broader implication of this field.

Autologous Exosome Transfer: A New Personalised Treatment Concept to Prevent Colitis in a Murine Model

BACKGROUND AND AIMS

Epigenetic information delivered by intestinal exosomes can be useful for diagnosing intestinal diseases, such as ulcerative colitis, but the therapeutic effects of intestinal exosomes have not been fully exploited. We herein developed an autologous exosome therapy that could treat intestinal disease without any risk of inducing a systemic immunological reaction.

METHODS

Intestinal exosomes were isolated and purified from faeces by our newly developed multi-step sucrose gradient ultracentrifugation method. Lipopolysaccharide [LPS]-activated macrophages were employed to test the in vitro anti-inflammatory ability of intestinal exosomes. To evaluate the in vivo anti-inflammatory activity of our system, we gavaged dextran sulphate sodium [DSS]-induced colitic mice with their own healing phase intestinal exosomes.

RESULTS

Mouse intestinal exosomes are round extracellular vesicles with a hydrodynamic diameter of ~140 [±20] nm and a surface charge of ~-12 [±3] mV. Among the exosomes obtained at four different stages of DSS-induced ulcerative colitis [1, before treatment; 2, DSS-treated; 3, healing phase; and 4, back to normal], the healing phase exosomes showed the best in vitro anti-inflammatory effects and promotion of wound healing. Moreover, oral co-administration of autologous healing phase exosomes with DSS was found to significantly reduce the risk of a second round of DSS-induced ulcerative colitis in mice.

CONCLUSIONS

Intestinal exosomes obtained during the healing phase that follows induced intestinal inflammation could strongly promote wound healing in the host. Oral administration of autologous exosomes from the healing phase could be a safe and effective approach for treating the ulcerative colitis of a given patient in the context of personalised medicine.

Loss of miR-101-3p Promotes Transmigration of Metastatic Breast Cancer Cells through the Brain Endothelium by Inducing COX-2/MMP1 Signaling

Brain metastases represent one of the incurable end stages in breast cancer (BC). Developing effective or preventive treatments is hampered by a lack of knowledge on the molecular mechanisms driving brain metastasis. Transmigration of BC cells through the brain endothelium is a key event in the pathogenesis of brain metastasis. In this study, we identified miR-101-3p as a critical micro-RNA able to reduce transmigration of BC cells through the brain endothelium. Our results revealed that miR-101-3p expression is downregulated in brain metastatic BC cells compared to less invasive variants, and varies inversely compared to the brain metastatic propensity of BC cells. Using a loss-and-gain of function approach, we found that miR-101-3p downregulation increased transmigration of BC cells through the brain endothelium in vitro by inducing COX-2 expression in cancer cells, whereas ectopic restoration of miR-101-3p exerted a metastasis-reducing effect. In regulatory experiments, we found that miR-101-3p mediated its effect by modulating COX-2-MMP1 signaling capable of degrading the inter-endothelial junctions (claudin-5 and VE-cadherin), key components of the brain endothelium. These findings suggest that miR-101-3p plays a critical role in the transmigration of breast cancer cells through the brain endothelium by modulating the COX-2-MMP1 signaling and thus may serve as a therapeutic target that can be exploited to prevent or suppress brain metastasis in human breast cancer.

Coding and Non-coding RNAs: Molecular Basis of Forest-Insect Outbreaks

Insect population dynamics are closely related to 'human' ecological and economic environments, and a central focus of research is outbreaks. However, the lack of molecular-based investigations restricts our understanding of the intrinsic mechanisms responsible for insect outbreaks. In this context, the moth Dendrolimus punctatus Walker can serve as an ideal model species for insect population dynamics research because it undergoes periodic outbreaks. Here, high-throughput whole-transcriptome sequencing was performed using D. punctatus, sampled during latent and outbreak periods, to systemically explore the molecular basis of insect outbreaks and to identify the involved non-coding RNA (ncRNA) regulators, namely microRNAs, long non-coding RNAs, and circular RNAs. Differentially expressed mRNAs of D. punctatus from different outbreak periods were involved in developmental, reproductive, immune, and chemosensory processes; results that were consistent with the physiological differences in D. punctatus during differing outbreak periods. Targets analysis of the non-coding RNAs indicated that long non-coding RNAs could be the primary ncRNA regulators of D. punctatus outbreaks, while circular RNAs mainly regulated synapses and cell junctions. The target genes of differentially expressed microRNAs mainly regulated the metabolic and reproductive pathways during the D. punctatus outbreaks. Developmental, multi-organismal, and reproductive processes, as well as biological adhesion, characterized the competing endogenous RNA network. Chemosensory and immune genes closely related to the outbreak of D. punctatus were further analyzed in detail: from their ncRNA regulators' analysis, we deduce that both lncRNA and miRNA may play significant roles. This is the first report to examine the molecular basis of coding and non-coding RNAs' roles in insect outbreaks. The results provide potential biomarkers for control targets in forest insect management, as well as fresh insights into underlying outbreak-related mechanisms, which could be used for improving insect control strategies in the future.

Reveal the Regulation Patterns of Prognosis-Related miRNAs and lncRNAs Across Solid Tumors in the Cancer Genome Atlas

BACKGROUND

The dysregulation of non-coding RNAs (ncRNAs) such as miRNAs and lncRNAs are associated with the pathogenesis and progression in multiple cancers including solid tumors. Comprehensive investigations of prognosis-related ncRNA markers could promote the development of therapeutic strategies for solid tumors, but rarely reported.

METHODS

By taking advantage of The Cancer Genome Atlas (TCGA), pan-cancer prognosis analysis (PCPA) models were firstly constructed based on miRNA and lncRNA expression profiles of 8,450 samples in 19 solid tumors. Further, the co-occurrence and exclusivity among ncRNA markers were systematically analyzed for different cancers.

RESULTS

In identified ncRNA makers, 71% of the miRNA markers were shared in multiple cancers, whereas 96% of the lncRNA markers were cancer-specific. Moreover, to analyze the regulation patterns of prognosis-related ncRNAs at the pan-cancer level, miRNA markers were further annotated into eight carcinogenic pathways. Results represented that approximately 86% of these miRNA markers could regulate the PI3K-Akt signaling pathway, while only 48% for the Notch signaling pathway. Finally, among 126 common genes that participated in eight carcinogenic pathways, BCL2, CSNK2A1, EGFR, PDGFRA, and VEGFA were proposed as potential drug targets for multiple cancers.

CONCLUSION

The prognosis analysis and regulation characteristics of ncRNAs presented in this study may help to facilitate the discovery of anti-cancer drugs for multiple solid tumors.

Identification of lncRNA and mRNA Biomarkers in Osteoarthritic Degenerative Meniscus by Weighted Gene Coexpression Network and Competing Endogenous RNA Network Analysis

BACKGROUND

Long noncoding RNAs (lncRNAs) play a crucial role in varieties of biological processes. This study is aimed at investigating meniscal degeneration-specific lncRNAs and mRNAs and their related networks in knee osteoarthritis (KOA).

METHODS

The dataset GSE98918, which included 24 meniscus samples and related clinical data, was downloaded from the Gene Expression Omnibus database. The differentially expressed lncRNAs and mRNAs in the meniscus between KOA and control groups were identified. Based on the enriched differentially expressed lncRNAs and mRNAs, we constructed the coexpression network using WGCNA (weighted correlation network analysis) and identified the critical module related to KOA. For mRNAs in the key module, gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were carried out using the DAVID database. A competing endogenous RNA network (ceRNA) based on the screened mRNAs, lncRNAs, and related miRNAs was constructed to reveal presumptive biomarkers further. Finally, the hub lncRNAs and mRNAs were screened, and the diagnostic value was evaluated using a receiver operating characteristic (ROC) curve. Hub mRNAs were validated using the dataset GSE113825.

RESULTS

We screened 208 significantly differentially expressed lncRNAs and mRNAs in menisci between the KOA and non-KOA samples, which were enriched in sixteen modules using WGCNA, especially the green module. Coexpression network based on the enriched differentially expressed lncRNAs and mRNAs in the green module uncovered 5 lncRNAs and 56 mRNAs. The lncRNA-miRNA-mRNA ceRNA network revealed that lnc-HLA-DQA1-5, lnc-RP11-127H5.1.1-1, lnc-RTN2-1, IGFBP4 (insulin-like growth factor binding protein 4), and KLF2 (Kruppel-like factor 2) were significantly correlated with the meniscus degeneration of KOA. ROC curve analysis revealed that these hub lncRNAs and mRNAs showed excellent diagnostic value for KOA.

CONCLUSIONS

These hub lncRNAs and mRNAs were potential prognostic biomarkers for the meniscus degeneration of KOA. Further studies are required to validate these new biomarkers and better understand the pathological process of the meniscus degeneration of KOA.

Prediction of miRNA targets by learning from interaction sequences

MicroRNAs (miRNAs) are involved in a diverse variety of biological processes through regulating the expression of target genes in the post-transcriptional level. So, it is of great importance to discover the targets of miRNAs in biological research. But, due to the short length of miRNAs and limited sequence complementarity to their gene targets in animals, it is challenging to develop algorithms to predict the targets of miRNA accurately. Here we developed a new miRNA target prediction algorithm using a multilayer convolutional neural network. Our model learned automatically the interaction patterns of the experiment-validated miRNA:target-site chimeras from the raw sequence, avoiding hand-craft selection of features by domain experts. The performance on test dataset is inspiring, indicating great generalization ability of our model. Moreover, considering the stability of miRNA:target-site duplexes, our method also showed good performance to predict the target transcripts of miRNAs.

Human Follicle-Stimulating Hormone ß Subunit Expression Depends on FOXL2 and SMAD4

Follicle-stimulating hormone (FSH), an essential regulator of mammalian fertility, is synthesized by pituitary gonadotrope cells in response to activins. In mice, activins signal via SMAD3, SMAD4, and FOXL2 to regulate transcription of the FSHβ subunit (Fshb) gene. Gonadotrope-specific deletion of Foxl2, alone or in combination with Smad4, renders mice FSH-deficient. Whether human FSHB expression is similarly regulated is not known. Here, we used a combination of transgenic and conditional knockout mouse strains to assess the roles of activins, FOXL2, and SMAD4 in regulation of the human FSHB gene. First, we cultured pituitaries from mice harboring a human FSHB transgene (hFSHB mice) and measured both murine Fshb and human FSHB messenger ribonucleic acid (mRNA) expression in response to exogenous activins or two antagonists of endogenous activin-like signaling (follistatin-288 and SB431542). Both murine Fshb and human FSHB expression were stimulated by activins and reduced by the inhibitors. Next, we analyzed human FSHB expression in hFSHB mice carrying floxed Foxl2 and Smad4 alleles. Cre-mediated ablation of FOXL2 and SMAD4 strongly reduced basal and activin-stimulated murine Fshb and human FSHB expression in cultured pituitaries. Finally, the hFSHB transgene was previously shown to rescue FSH production and fertility in Fshb knockout mice. However, gonadotrope-specific Foxl2/Smad4 knockout females carrying the hFSHB transgene have significantly reduced murine Fshb and human FSHB pituitary mRNA levels and are hypogonadal. Collectively, these data suggest that similar to Fshb regulation in mice, FOXL2 and SMAD4 play essential roles in human FSHB expression.

LncRNA LINC00483 promotes gastric cancer development through regulating MAPK1 expression by sponging miR-490-3p

Background

Previous studies have shown that long noncoding RNA (lncRNA) LINC00483 was aberrantly expressed in human cancers, including gastric cancer. However, the regulatory mechanism of this lncRNA in gastric cancer remains largely unknown. The present study aimed to investigate the effect of LINC00483 on gastric cancer development and explore the potential regulatory network of LINC00483/microRNA (miR)-490-3p/mitogen-activated protein kinase 1 (MAPK1).

Methods

Thirty patients with gastric cancer were recruited for tissues collection. The expression levels of LINC00483, miR-490-3p and MAPK1 were detected by quantitative real-time polymerase chain reaction or western blot. Cell viability, apoptosis, migration and invasion were determined by MTT, flow cytometry, transwell assays and western blot, respectively. The target association between miR-490-3p and LINC00483 or MAPK1 was confirmed by luciferase reporter assay. Xenograft model was established to assess the function of LINC00483 in vivo.

Results

LINC00483 and MAPK1 levels were increased in gastric cancer tissues and cells. Knockdown of LINC00483 or MAPK1 inhibited cells viability, migration and invasion but promoted apoptosis in gastric cancer cells. Moreover, MAPK1 overexpression attenuated the effect of LINC00483 knockdown on gastric cancer development. LINC00483 could increase MAPK1 expression by competitively sponging miR-490-3p. miR-490-3p overexpression suppressed gastric cancer development, which was abated by introduction of LINC00483. Besides, inhibition of LINC00483 decreased xenograft tumor growth by regulating miR-490-3p/MAPK1 axis.

Conclusion

Knockdown of LINC00483 inhibited gastric cancer development in vitro and in vivo by increasing miR-490-3p and decreasing MAPK1, elucidating a novel mechanism for understanding the development of gastric cancer.

Weighted gene correlation network analysis reveals novel biomarkers associated with mesenchymal stromal cell differentiation in early phase

Osteoporosis is a major public health problem that is associated with high morbidity and mortality, and its prevalence is increasing as the world’s population ages. Therefore, understanding the molecular basis of the disease is becoming a high priority. In this regard, studies have shown that an imbalance in adipogenic and osteogenic differentiation of bone marrow mesenchymal stem cells (MSCs) is associated with osteoporosis. In this study, we conducted a Weighted Gene Co-Expression Network Analysis to identify gene modules associated with the differentiation of bone marrow MSCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genome enrichment analysis showed that the most significant module, the brown module, was enriched with genes involved in cell cycle regulation, which is in line with the initial results published using these data. In addition, the Cytoscape platform was used to identify important hub genes and lncRNAs correlated with the gene modules. Furthermore, differential gene expression analysis identified 157 and 40 genes that were upregulated and downregulated, respectively, after 3 h of MSCs differentiation. Interestingly, regulatory network analysis, and comparison of the differentially expressed genes with those in the brown module identified potential novel biomarker genes, including two transcription factors (ZNF740, FOS) and two hub genes (FOXQ1, SGK1), which were further validated for differential expression in another data set of differentiation of MSCs. Finally, Gene Set Enrichment Analysis suggested that the two most important candidate hub genes are involved in regulatory pathways, such as the JAK-STAT and RAS signaling pathways. In summary, we have revealed new molecular mechanisms of MSCs differentiation and identified novel genes that could be used as potential therapeutic targets for the treatment of osteoporosis.

A review of databases predicting the effects of SNPs in miRNA genes or miRNA-binding sites

Modern precision medicine comprises the knowledge and understanding of individual differences in the genomic sequence of patients to provide tailor-made treatments. Regularly, such variants are considered in coding regions only, and their effects are predicted based on their impact on the amino acid sequence of expressed proteins. However, assessing the effects of variants in noncoding elements, in particular microRNAs (miRNAs) and their binding sites, is important as well, as a single miRNA can influence the expression patterns of many genes at the same time. To analyze the effects of variants in miRNAs and their target sites, several databases storing variant impact predictions have been published. In this review, we will compare the core functionalities and features of these databases and discuss the importance of up-to-date data resources in the context of web applications. Finally, we will outline some recommendations for future developments in the field.

Comparison of Adult Testis and Ovary MicroRNA Expression Profiles in Reeves' Pond Turtles (Mauremys reevesii) With Temperature-Dependent Sex Determination

Some differentially expressed genes (DEGs) that encode key enzymes involved in steroidogenic biosynthesis (CYP19A1) and key molecules related to gonadal functions (DMRT1, SOX9, AMH, FOXL2, WNT4, RSPO2, and GDF9) have been identified in adult gonadal RNA-seq studies of Reeves' pond turtle (Mauremys reevesii) with temperature-dependent sex determination (TSD). Gonadal functional maintenance and gametogenesis comprises a highly regulated and coordinated biological process, and increasing evidence indicates that microRNAs (miRNAs) may be involved in this dynamic program. However, it is not clear how the regulatory network comprising miRNAs changes the expression levels of these genes. In this study, miRNA sequencing of adult testis and ovary tissues from M. reevesii detected 25 known and 379 novel miRNAs, where 60 miRNAs were differentially expressed in the testis and ovary. A total of 1,477 target genes based on the differentially expressed miRNAs were predicted, where 221 target genes also exhibited differential expression. To verify the accuracy of the sequencing data, 10 differentially expressed miRNAs were validated by quantitative reverse transcription real-time PCR, and were found to be consistent with the transcriptome sequencing results. Moreover, several miRNA/target gene pairs, i.e., mre-let-7a-5p/mre-let-7e-5p and CYP19A1, mre-miR-200a-3p and DMRT1, mre-miR-101-3p and SOX9, and mre-miR-138-5p and AMH were identified. To explore the regulatory role of miRNAs, we conducted target gene enrichment analysis of the miRNAs and 221 target genes in the regulatory network. The signaling pathways related to gonadal functional maintenance and gametogenesis based on the DEGs and target genes were then compared. Our findings provide crucial information to facilitate further research into the regulatory mechanisms involving miRNAs in turtle species with TSD.

Integrated miRNA/mRNA Counter-Expression Analysis Highlights Oxidative Stress-Related Genes CCR7 and FOXO1 as Blood Markers of Coronary Arterial Disease

Our interest in the mechanisms of atherosclerosis progression (ATHp) has led to the recent identification of 13 miRNAs and 1285 mRNAs whose expression was altered during ATHp. Here, we deepen the functional relationship among these 13 miRNAs and genes associated to oxidative stress, a crucial step in the onset and progression of vascular disease. We first compiled a list of genes associated to the response to oxidative stress (Oxstress genes) by performing a reverse Gene Ontology analysis (rGO, from the GO terms to the genes) with the GO terms GO0006979, GO1902882, GO1902883 and GO1902884, which included a total of 417 unique Oxstress genes. Next, we identified 108 putative targets of the 13 miRNAs among these unique Oxstress genes, which were validated by an integrated miRNA/mRNA counter-expression analysis with the 1285 mRNAs that yielded 14 genes, Map2k1, Mapk1, Mapk9, Dapk1, Atp2a2, Gata4, Fos, Egfr, Foxo1, Ccr7, Vkorc1l1, Rnf7, Kcnh3, and Mgat3. GO enrichment analysis and a protein–protein-interaction network analysis (PPI) identified most of the validated Oxstress transcripts as components of signaling pathways, highlighting a role for MAP signaling in ATHp. Lastly, expression of these Oxstress transcripts was measured in PBMCs from patients suffering severe coronary artery disease, a serious consequence of ATHp. This allowed the identification of FOXO1 and CCR7 as blood markers downregulated in CAD. These results are discussed in the context of the interaction of the Oxstress transcripts with the ATHp-associated miRNAs.

Prediction of the miRNA interactome - Established methods and upcoming perspectives

MicroRNAs (miRNAs) are well-studied small noncoding RNAs involved in post-transcriptional gene regulation in a wide range of organisms, including mammals. Their function is mediated by base pairing with their target RNAs. Although many features required for miRNA-mediated repression have been described, the identification of functional interactions is still challenging. In the last two decades, numerous Machine Learning (ML) models have been developed to predict their putative targets. In this review, we summarize the biological knowledge and the experimental data used to develop these ML models. Recently, Deep Neural Network-based models have also emerged in miRNA interaction modeling. We thus outline established and emerging models to give a perspective on the future developments needed to improve the identification of genes directly regulated by miRNAs.

Retinal biomarkers and pharmacological targets for Hermansky-Pudlak syndrome 7

Deletion of dystrobrevin binding protein 1 has been linked to Hermansky-Pudlak syndrome type 7 (HPS-7), a rare disease characterized by oculocutaneous albinism and retinal dysfunction. We studied dysbindin-1 null mutant mice (Dys^−/−) to shed light on retinal neurodevelopment defects in HPS-7. We analyzed the expression of a focused set of miRNAs in retina of wild type (WT), Dys^+/− and Dys^−/− mice. We also investigated the retinal function of these mice through electroretinography (ERG). We found that miR-101-3p, miR-137, miR-186-5p, miR-326, miR-382-5p and miR-876-5p were up-regulated in Dys^−/−mice retina. Dys^−/− mice showed significant increased b-wave in ERG, compared to WT mice. Bioinformatic analysis highlighted that dysregulated miRNAs target synaptic plasticity and dopaminergic signaling pathways, affecting retinal functions of Dys^−/− mice. Overall, the data indicate potential mechanisms in retinal neurodevelopment of Dys^−/− mice, which may have translational significance in HSP-7 patients, both in terms of diagnostic/prognostic biomarkers and novel pharmacological targets.

Identification of microRNAs and their target genes in the placenta as biomarkers of inflammation

OBJECTIVE

Recently, microRNA (miRNA) has been identified both as a powerful regulator involved in various biological processes through the regulation of numerous genes and as an effective biomarker for the prediction and diagnosis of various disease states. The objective of this study was to identify and validate miRNAs and their target genes involved in inflammation in placental tissue.

METHODS

Microarrays were utilized to obtain miRNA and gene expression profiles from placentas with or without inflammation obtained from nine normal pregnant women and 10 preterm labor patients. Quantitative real-time polymerase chain reaction and Western blots were performed to validate the miRNAs and differentially-expressed genes in the placentas with inflammation. Correlations between miRNA and target gene expression were confirmed by luciferase assays in HTR-8/SVneo cells.

RESULTS

We identified and validated miRNAs and their target genes that were differentially expressed in placentas with inflammation. We also demonstrated that several miRNAs (miR-371a-5p, miR-3065-3p, miR-519b-3p, and miR-373-3p) directly targeted their target genes (LEF1, LOX, ITGB4, and CD44). However, some miRNAs and their direct target genes showed no correlation in tissue samples. Interestingly, miR-373-3p and miR-3065-3p were markedly regulated by lipopolysaccharide (LPS) treatment, although the expression of their direct targets CD44 and LOX was not altered by LPS treatment.

CONCLUSION

These results provide candidate miRNAs and their target genes that could be used as placental biomarkers of inflammation. These candidates may be useful for further miRNA-based biomarker development.

Single nucleotide alterations in MicroRNAs and human cancer-A not fully explored field

MicroRNAs are ~20 nt long small noncoding RNAs that are processed from stem-looped precursors and function mainly as posttranscriptional regulators of protein coding genes through binding to 3'-untranslated regions of messenger RNAs to inhibit the translation or cause RNA degradation. It is predicted microRNAs could regulate up to half of all human genes and are proved to play important roles in human diseases including cancer. They bind to target mRNAs based on complementary binding which is dominated by the so-called "seed" region which are the 5' 2-8 bases of the microRNA. Due to the small size in nature, even a single nucleotide variation in the precursor region especially those located in the seed regions could show big influence. Here, I summarized and reviewed the current knowledge of these single nucleotide alterations in microRNAs in human cancer including (i) common SNPs in the precursor region, (ii) isomiRs, (iii) somatic mutations of microRNAs. Briefly, this is an underexploited field and clearly, warrants further studies to reveal their biological and clinical significances. I believe they will be key to advancing personalized medicine.

The expression of small RNAs in exosomes of follicular fluid altered in human polycystic ovarian syndrome

Polycystic ovary syndrome (PCOS) can cause reproductive disorders that may affect oocyte quality from punctured follicles in human follicular fluid (HFF). The non-coding RNA family includes micro RNA (miRNA), piwi-interacting RNA (piRNA) and transfer RNA (tRNA); these non-coding RNA transcripts play diverse functions and are implicated in a variety of diseases and health conditions, including infertility. In this study, to explore the role of HFF exosomes in PCOS, we extracted and sequenced RNA from HFF exosomes of PCOS patients and compared the analysis results with those of non-PCOS control group. The HFF exosomes were successfully isolated and characterized in a variety of ways. The sequencing results of the HFF exosomal RNA showed that about 6.6% of valid reads in the PCOS group and 8.6% in the non-PCOS group were successfully mapped to the human RNA database. Using a hierarchical clustering method, we found there were ten small RNA sequences whose expression was significantly different between the PCOS and non-PCOS groups. We chose six of them to predict target genes of interest for further GO analysis, and pathway analysis showed that the target genes are mainly involved in biosynthesis of amino acids, glycine, serine and glycosaminoglycan, as well as threonine metabolism. Therefore, the small RNA sequences contained in HFF EXs may play a key role in the mechanism that drives PCOS pathogenesis, and thereby can act as molecular biomarkers for PCOS diagnosis in the future.

Restoration of microRNA-197 expression suppresses oncogenicity in fibrosarcoma through negative regulation of RAN

MicroRNAs (miRNAs) act as crucial regulators of biological pathways/processes by reinforcing transcriptional programs and moderating transcripts. Emerging evidences have shown the involvement of dysregulated miRNAs in pathophysiology of human diseases including several cancer types. Recently, miR-197-3p has been reported to play different roles in different cancers; however, its role in fibrosarcoma, a highly aggressive and malignant soft tissue sarcoma originated from the mesenchymal tissues, has not yet been studied. Therefore, this study aims to investigate the possible regulatory roles of miR-197-3p in the oncogenicity of fibrosarcoma. For this, we initially performed qRT-PCR of miR-197-3p, which we found to be downregulated in HT1080 human fibrosarcoma cells compared with IMR90-tert normal fibroblast cells. Subsequently, we performed gain-of-function study by employing several methods such as MTT assay, clonogenic assay, wound healing, flow cytometry cell cycle analysis, and acridine orange staining after transfecting HT1080 cells with miR-197-3p mimic. From these assays, we observed that miR-197-3p significantly inhibits viability, colony forming, and migration ability as well as triggers G2/M phase cell cycle arrest and autophagy in fibrosarcoma cells. To understand the mechanism through which miRNA performs these functions, we predicted its targets using TargetScan and performed pathway enrichment analysis after screening them by their expression in fibrosarcoma. Among the enriched targets, we found RAN (ras-related nuclear protein) to be a crucial target through which miR-197-3p represses tumorigenesis by binding to its 3´ UTR, validated by luciferase reporter assay. The tumor suppressive role of the miRNA was further confirmed by transfecting its mimic in RAN-overexpressed cells which showed significant attenuation in tumorigenic effect of RAN in fibrosarcoma as seen in different assays. Taken together, our study unveiled that miR-197-3p acts as an oncosuppressor in fibrosarcoma through G2/M phase arrest and induction of autophagy, and raises the possibility to act as a novel therapeutic intervention for the malignancy.

Downregulation of circulating miR 802-5p and miR 194-5p and upregulation of brain MEF2C along breast cancer brain metastasization

Breast cancer brain metastases (BCBMs) have been underinvestigated despite their high incidence and poor outcome. MicroRNAs (miRNAs), and particularly circulating miRNAs, regulate multiple cellular functions, and their deregulation has been reported in different types of cancer and metastasis. However, their signature in plasma along brain metastasis development and their relevant targets remain undetermined. Here, we used a mouse model of BCBM and next‐generation sequencing (NGS) to establish the alterations in circulating miRNAs during brain metastasis formation and development. We further performed bioinformatics analysis to identify their targets with relevance in the metastatic process. We additionally analyzed human resected brain metastasis samples of breast cancer patients for target expression validation. Breast cancer cells were injected in the carotid artery of mice to preferentially induce metastasis in the brain, and samples were collected at different timepoints (5 h, 3, 7, and 10 days) to follow metastasis development in the brain and in peripheral organs. Metastases were detected from 7 days onwards, mainly in the brain. NGS revealed a deregulation of circulating miRNA profile during BCBM progression, rising from 18% at 3 days to 30% at 10 days following malignant cells’ injection. Work was focused on those altered prior to metastasis detection, among which were miR‐802‐5p and miR‐194‐5p, whose downregulation was validated by qPCR. Using targetscan and diana tools, the transcription factor myocyte enhancer factor 2C (MEF2C) was identified as a target for both miRNAs, and its expression was increasingly observed in malignant cells along brain metastasis development. Its upregulation was also observed in peritumoral astrocytes pointing to a role of MEF2C in the crosstalk between tumor cells and astrocytes. MEF2C expression was also observed in human BCBM, validating the observation in mouse. Collectively, downregulation of circulating miR‐802‐5p and miR‐194‐5p appears as a precocious event in BCBM and MEF2C emerges as a new player in brain metastasis development.

Unveiling ncRNA regulatory axes in atherosclerosis progression

Completion of the human genome sequencing project highlighted the richness of the cellular RNA world, and opened the door to the discovery of a plethora of short and long non-coding RNAs (the dark transcriptome) with regulatory or structural potential, which shifted the balance of pathological gene alterations from coding to non-coding RNAs. Thus, disease risk assessment currently has to also evaluate the expression of new RNAs such as small micro RNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), competing endogenous RNAs (ceRNAs), retrogressed elements, 3'UTRs of mRNAs, etc. We are interested in the pathogenic mechanisms of atherosclerosis (ATH) progression in patients suffering Chronic Kidney Disease, and in this review, we will focus in the role of the dark transcriptome (non-coding RNAs) in ATH progression. We will focus in miRNAs and in the formation of regulatory axes or networks with their mRNA targets and with the lncRNAs that function as miRNA sponges or competitive inhibitors of miRNA activity. In this sense, we will pay special attention to retrogressed genomic elements, such as processed pseudogenes and Alu repeated elements, that have been recently seen to also function as miRNA sponges, as well as to the use or miRNA derivatives in gene silencing, anti-ATH therapies. Along the review, we will discuss technical developments associated to research in lncRNAs, from sequencing technologies to databases, repositories and algorithms to predict miRNA targets, as well as new approaches to miRNA function, such as integrative or enrichment analysis and their potential to unveil RNA regulatory networks.

The analysis of a ceRNA network and the correlation between lncRNA, miRNA, and mRNA in bladder cancer

Background

To explore the correlation between the lncRNA-miRNA-mRNA and ceRNA network through the differential expression analysis of lncRNAs, miRNAs and mRNAs in bladder cancer based on The Cancer Genome Atlas (TCGA) database combined with Gene Ontology (GO) and Kyoto Encyclopedia of Genes Genomes (KEGG) enrichment analysis.

Methods

Firstly, the expression profile data and corresponding clinical data of RNAs in bladder cancer were searched and downloaded from TCGA database, and aberrantly expressed long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) were screened and found by using TCGA database. The relationship between lncRNA-miRNA-mRNA was established by comparing these lncRNAs, miRNAs, and mRNAs, while the ceRNA network was constructed. Combined with the analysis of the GO annotation and KEGG pathway, the effects of lncRNA-miRNA-mRNA interaction on the development of bladder cancer were explored.

Results

A total of 1,742 differentially expressed lncRNA, 511 differentially expressed miRNAs, and 4,373 differentially expressed mRNAs were identified, and 328 lncRNAs, 73 miRNAs, and 677 mRNAs were screened by survival analysis. With the lncRNA-miRNA-mRNA correlation analysis, a ceRNA network consisting of 45 lncRNAs, 14 miRNAs, and 29 mRNAs was successfully constructed. The GO annotation and functional enrichment of target gene mRNAs in the network are mainly concentrated in the signal pathways and include fatty acid biosynthesis, gap junction, insulin signaling pathway, and the MAPK signaling pathway biological processes such as positive regulation of cellular process and system development.

Conclusions

We successfully identified the target gene correlating lncRNA, miRNA, and mRNA, and constructed a ceRNA network. Our findings can provide a potential target for the study of the occurrence, development, diagnosis, treatment, and prognosis of bladder cancer.

A novel computational approach to the silencing of Sugarcane Bacilliform Guadeloupe A Virus determines potential host-derived MicroRNAs in sugarcane (Saccharum officinarum L.)

Sugarcane Bacilliform Guadeloupe A Virus (SCBGAV, genus Badnavirus, family Caulimoviridae) is an emerging, deleterious pathogen of sugarcane which presents a substantial barrier to producing high sugarcane earnings. Sugarcane bacilliform viruses (SCBVs) are one of the main species that infect sugarcane. During the last 30 years, significant genetic changes in SCBV strains have been observed with a high risk of disease incidence associated with crop damage. SCBV infection may lead to significant losses in biomass production in susceptible sugarcane cultivars. The circular, double-stranded (ds) DNA genome of SCBGAV (7.4 Kb) is composed of three open reading frames (ORFs) on the positive strand that replicate by a reverse transcriptase. SCBGAV can infect sugarcane in a semipersistent manner via the insect vectors sugarcane mealybug species. In the current study, we used miRNA target prediction algorithms to identify and comprehensively analyze the genome-wide sugarcane (Saccharum officinarum L.)-encoded microRNA (miRNA) targets against the SCBGAV. Mature miRNA target sequences were retrieved from the miRBase (miRNA database) and were further analyzed for hybridization to the SCBGAV genome. Multiple computational approaches—including miRNA-target seed pairing, multiple target positions, minimum free energy, target site accessibility, maximum complementarity, pattern recognition and minimum folding energy for attachments—were considered by all algorithms. Among them, sof-miR396 was identified as the top effective candidate, capable of targeting the vital ORF3 of the SCBGAV genome. miRanda, RNA22 and RNAhybrid algorithms predicted hybridization of sof-miR396 at common locus position 3394. The predicted sugarcane miRNAs against viral mRNA targets possess antiviral activities, leading to translational inhibition by mRNA cleavage. Interaction network of sugarcane-encoded miRNAs with SCBGAV genes, created using Circos, allow analyze new targets. The finding of the present study acts as a first step towards the creation of SCBGAV-resistant sugarcane through the expression of the identified miRNAs.