Identification of microRNA-regulated gene networks by expression analysis of target genes

The Relevance of Reperfusion Stroke Therapy for miR-9-3p and miR-9-5p Expression in Acute Stroke-A Preliminary Study

Reperfusion stroke therapy is a modern treatment that involves thrombolysis and the mechanical removal of thrombus from the extracranial and/or cerebral arteries, thereby increasing penumbra reperfusion. After reperfusion therapy, 46% of patients are able to live independently 3 months after stroke onset. MicroRNAs (miRNAs) are essential regulators in the development of cerebral ischemia/reperfusion injury and the efficacy of the applied treatment. The first aim of this study was to examine the change in serum miRNA levels via next-generation sequencing (NGS) 10 days after the onset of acute stroke and reperfusion treatment. Next, the predictive values of the bioinformatics analysis of miRNA gene targets for the assessment of brain ischemic response to reperfusion treatment were explored. Human serum samples were collected from patients on days 1 and 10 after stroke onset and reperfusion treatment. The samples were subjected to NGS and then validated using qRT-PCR. Differentially expressed miRNAs (DEmiRNAs) were used for enrichment analysis. Hsa-miR-9-3p and hsa-miR-9-5p expression were downregulated on day 10 compared to reperfusion treatment on day 1 after stroke. The functional analysis of miRNA target genes revealed a strong association between the identified miRNA and stroke-related biological processes related to neuroregeneration signaling pathways. Hsa-miR-9-3p and hsa-miR-9-5p are potential candidates for the further exploration of reperfusion treatment efficacy in stroke patients.

Bioinformatic Analysis of miR-200b/429 and Hub Gene Network in Cervical Cancer

The miR-200b/429 located at 1p36 is a highly conserved miRNA cluster emerging as a critical regulator of cervical cancer. Using publicly available miRNA expression data from TCGA and GEO followed by independent validation, we aimed to identify the association between miR-200b/429 expression and cervical cancer. miR-200b/429 cluster was significantly overexpressed in cancer samples compared to normal samples. miR-200b/429 expression did not correlate with patient survival; however, its overexpression correlated with histological type. Protein-protein interaction analysis of 90 target genes of miR-200b/429 identified EZH2, FLT1, IGF2, IRS1, JUN, KDR, SOX2, MYB, ZEB1, and TIMP2 as the top ten hub genes. PI3K-AKT and MAPK signaling pathways emerged as major target pathways of miR-200b/429 and their hub genes. Kaplan-Meier survival analysis showed the expression of seven miR-200b/429 target genes (EZH2, FLT1, IGF2, IRS1, JUN, SOX2, and TIMP2) to influence the overall survival of patients. The miR-200a-3p and miR-200b-5p could help predict cervical cancer with metastatic potential. The cancer hallmark enrichment analysis showed hub genes to promote growth, sustained proliferation, resistance to apoptosis, induction of angiogenesis, activation of invasion, and metastasis, enabling replicative immortality, evading immune destruction, and tumor-promoting inflammation. The drug-gene interaction analysis identified 182 potential drugs to interact with 27 target genes of miR-200b/429 with paclitaxel, doxorubicin, dabrafenib, bortezomib, docetaxel, ABT-199, eribulin, vorinostat, etoposide, and mitoxantrone emerging as the top ten best candidate drugs. Taken together, miR-200b/429 and associated hub genes can be helpful for prognostic application and clinical management of cervical cancer.

Specific upregulation of extracellular miR-6238 in particulate matter-induced acute lung injury and its immunomodulation

Acute lung injury (ALI) and acute respiratory distress syndrome (ARDS) are life-threatening diseases characterized by a severe inflammatory response and the destruction of alveolar epithelium and endothelium. ALI/ARDS is caused by pathogens and toxic environmental stimuli, such as particulate matter (PM). However, the general symptoms of ALI/ARDS are similar, and determining the cause of lung injury is often challenging. In this study, we investigated whether there is a critical miRNA that characterizes PM-induced ALI. We found that the expression of miR-6238 is specifically upregulated in lung tissue and lung-derived extracellular vesicles (EVs) in response to PM exposure. Notably, bacterial endotoxin (Lipopolysaccharide; LPS or peptidoglycan; PTG) does not induce the expression of miR-6238 in the lung. Instead, the expression of miR-155 is dramatically increased in LPS-induced ALI. We further demonstrated that human lung epithelial cells and macrophages predominantly produce miR-6238 and miR-155, respectively. Mechanistically, EV-miR-6238 is effectively internalized into alveolar macrophages (AMs) and regulates inflammatory responses in vivo. CXCL3 is a main target of miR-6238 in AMs and modulates neutrophil infiltration into the lung alveoli. Collectively, our findings suggest that miR-6238 is a novel regulator of pulmonary inflammation and a putative biomarker that distinguishes PM-induced ALI from endotoxin (LPS/PTG)-mediated ALI.

Chromosome 1p status in neuroblastoma correlates with higher expression levels of miRNAs targeting neuronal differentiation pathway

Neuroblastoma (NB) is characterized by acquired segmental and numerical chromosome aberrations. Although deletions of distal 1p and 11q are frequent alterations, no candidate tumor suppressor gene residing in these chromosomal sites could be identified so far. In the present study, we detected the genomic imbalances of six neuroblastoma cell lines using the multiplex ligation-dependent probe amplification (MLPA) technique and the microRNA (miRNA) expression profiles of the cell lines by a microarray study. According to MLPA results, we aimed to assess the miRNA expression profiles of the cell lines harboring 11q and 1p deletions. The cell lines with 1p deletions revealed statistically significant higher levels of expression for 29 miRNAs in contrast to the cell lines without 1p deletion in microarray study. We also performed GO enrichment analysis for predicted targets of the differentially expressed miRNAs. According to GO enrichment analysis, miRNAs that showed the high change in expression was associated with neuronal differentiation. We showed that hsa-miR-494, hsa-miR-495, and hsa-miR-543 target most of mRNAs in neuronal differentiation pathway. Although limited to the cell lines, our results highly suggest that NBs with different segmental chromosome abnormalities may have different dysregulated miRNA expression signatures that target the genes involved in neuronal differentiation.

Identification of non-coding RNA related prognosis biomarkers based on ceRNA network in thyroid cancer

Introduction: Thyroid cancer (THCA) has become a serious malignant tumor worldwide. Identification of non-coding RNA related regulators is very necessary to improve the knowledge of THCA treatment. The aim of this study was to identify novel therapeutic targets and prognosis biomarkers for predicting pathological characteristics and subsequently treating THCA. Methods: We investigated the alterations of miRNAs, mRNAs and lncRNAs in THCA. Functional enrichment and clustering analysis were conducted for these aberrantly expressed RNAs. Multiple interaction networks among miRNAs, mRNAs and lncRNAs were constructed and the functional modules associated with THCA patients' prognosis were identified. Furthermore, we evaluated the prognostic roles of the important miRNAs, mRNAs and lncRNAs in THCA and investigated the regulatory potential of non-coding RNAs on immune cell infiltration. Results: We firstly identified that miR-4709-3p and miR-146b-3p could significantly classify patients into high/low risk groups, which may be potential prognosis biomarkers of THCA. Secondly, we constructed a THCA-related miRNA-mRNA network, which displayed small world network topological characters. Two THCA-related functional modules were identified from the miRNA-mRNA network by MCODE. Results showed that two modules could implicate in known cancer pathways, such as apoptosis and focal adhesion. Thirdly, a THCA-related miRNA-lncRNA network was constructed. A subnetwork of miRNA-lncRNA network showed strong prognosis effect in THCA. Fourthly, we constructed a THCA-related mRNA-lncRNA network and detected several typical lncRNA-miRNA-mRNA crosstalk, such as AC068138, BCL2, miR-21 and miR-146b, which had good prognosis effect in THCA. Immune infiltration results showed that lncRNAs LA16c-329F2, RP11-395N3, RP11-423H2, RP11-399B17 and RP11-1036E20 were high related to neutrophil and dendritic cell infiltration. Discussion: Non-coding RNA-mediated gene regulatory network has the strong regulatory potential in pathological processes of THCA. All these results could help us uncover the non-coding RNA-mediated regulatory mechanism in THCA.

The Role of Cluster C19MC in Pre-Eclampsia Development

Pre-eclampsia is a placenta-related complication occurring in 2-10% of all pregnancies. miRNAs are a group of non-coding RNAs regulating gene expression. There is evidence that C19MC miRNAs are involved in the development of the placenta. Deregulation of chromosome 19 microRNA cluster (C19MC) miRNAs expression leads to impaired cell differentiation, abnormal trophoblast invasion and pathological angiogenesis, which can lead to the development of pre-eclampsia. Information was obtained through a review of articles available in PubMed Medline. Articles on the role of the C19MC miRNA in the development of pre-eclampsia published in 2009-2022 were analyzed. This review article summarizes the current data on the role of the C19MC miRNA in the development of pre-eclampsia. They indicate a significant increase in the expression of most C19MC miRNAs in placental tissue and a high level of circulating fractions in serum and plasma, both in the first and/or third trimester in women with PE. Only for miR-525-5p, low levels of plasma expression were noted in the first trimester, and in the placenta in the third trimester. The search for molecular factors indicating the development of pre-eclampsia before the onset of clinical symptoms seems to be a promising diagnostic route. Identifying women at risk of developing pre-eclampsia at the pre-symptomatic stage would avoid serious complications in both mothers and fetuses. We believe that miRNAs belonging to cluster C19MC could be promising biomarkers of pre-eclampsia development.

Identification of miRNA-mRNA-TFs regulatory network and crucial pathways involved in asthma through advanced systems biology approaches

Asthma is a life-threatening and chronic inflammatory lung disease that is posing a true global health challenge. The genetic basis of the disease is fairly well examined. However, the molecular crosstalk between microRNAs (miRNAs), target genes, and transcription factors (TFs) networks and their contribution to disease pathogenesis and progression is not well explored. Therefore, this study was aimed at dissecting the molecular network between mRNAs, miRNAs, and TFs using robust computational biology approaches. The transcriptomic data of bronchial epithelial cells of severe asthma patients and healthy controls was studied by different systems biology approaches like differentially expressed gene detection, functional enrichment, miRNA-target gene pairing, and mRNA-miRNA-TF molecular networking. We detected the differential expression of 1703 (673 up-and 1030 down-regulated) genes and 71 (41 up-and 30 down-regulated) miRNAs in the bronchial epithelial cells of asthma patients. The DEGs were found to be enriched in key pathways like IL-17 signaling (KEGG: 04657), Th1 and Th2 cell differentiation (KEGG: 04658), and the Th17 cell differentiation (KEGG: 04659) (p-values = 0.001). The results from miRNAs-target gene pairs-transcription factors (TFs) have detected the key roles of 3 miRs (miR-181a-2-3p; miR-203a-3p; miR-335-5p), 6 TFs (TFAM, FOXO1, GFI1, IRF2, SOX9, and HLF) and 32 miRNA target genes in eliciting autoimmune reactions in bronchial epithelial cells of the respiratory tract. Through systemic implementation of comprehensive system biology tools, this study has identified key miRNAs, TFs, and miRNA target gene pairs as potential tissue-based asthma biomarkers.

Impact of smoking-induced dysregulated human miRNAs in chronic disease development and their potential use in prognostic and therapeutic purposes

MicroRNAs (miRNAs) are evolutionary conserved small noncoding RNA molecules with a significant ability to regulate gene expression at the posttranscriptional level either through translation repression or messenger RNA degradation. miRNAs are differentially expressed in various pathophysiological conditions, affecting the course of the disease by modulating several critical target genes. As the persistence of irreversible molecular changes caused by cigarette smoking is central to the pathogenesis of various chronic diseases, several studies have shown its direct correlation with the dysregulation of different miRNAs, affecting numerous essential biological processes. This review provides an insight into the current status of smoking-induced miRNAs dysregulation in chronic diseases such as COPD, atherosclerosis, pulmonary hypertension, and different cancers and explores the diagnostic/prognostic potential of miRNA-based biomarkers and their efficacy as therapeutic targets.

Upregulation of BMI1-suppressor miRNAs (miR-200c, miR-203) during terminal differentiation of colon epithelial cells

BACKGROUND

MicroRNAs (miRNAs) are key regulators of stem cell functions, including self-renewal and differentiation. In this study, we aimed to identify miRNAs that are upregulated during terminal differentiation in the human colon epithelium, and elucidate their role in the mechanistic control of stem cell properties.

METHODS

"Bottom-of-the-crypt" (EPCAM⁺/CD44⁺/CD66a^low) and "top-of-the-crypt" (EPCAM⁺/CD44^neg/CD66a^high) epithelial cells from 8 primary colon specimens (6 human, 2 murine) were purified by flow cytometry and analyzed for differential expression of 335 miRNAs. The miRNAs displaying the highest upregulation in "top-of-the-crypt" (terminally differentiated) epithelial cells were tested for positive correlation and association with survival outcomes in a colon cancer RNA-seq database (n = 439 patients). The two miRNAs with the strongest "top-of-the-crypt" expression profile were evaluated for capacity to downregulate self-renewal effectors and inhibit in vitro proliferation of colon cancer cells, in vitro organoid formation by normal colon epithelial cells and in vivo tumorigenicity by patient-derived xenografts (PDX).

RESULTS

Six miRNAs (miR-200a, miR-200b, miR-200c, miR-203, miR-210, miR-345) were upregulated in "top-of-the-crypt" cells and positively correlated in expression among colon carcinomas. Overexpression of the three miRNAs with the highest inter-correlation coefficients (miR-200a, miR-200b, miR-200c) associated with improved survival. The top two over-expressed miRNAs (miR-200c, miR-203) cooperated synergistically in suppressing expression of BMI1, a key regulator of self-renewal in stem cell populations, and in inhibiting proliferation, organoid-formation and tumorigenicity of colon epithelial cells.

CONCLUSION

In the colon epithelium, terminal differentiation associates with the coordinated upregulation of miR-200c and miR-203, which cooperate to suppress BMI1 and disable the expansion capacity of epithelial cells.

Integrated computational analysis reveals HOX genes cluster as oncogenic drivers in head and neck squamous cell carcinoma

Alterations in homeobox (HOX) gene expression are involved in the progression of several cancer types including head and neck squamous cell carcinoma (HNSCC). However, regulation of the entire HOX cluster in the pathophysiology of HNSCC is still elusive. By using different comprehensive databases, we have identified the significance of differentially expressed HOX genes (DEHGs) in stage stratification and HPV status in the cancer genome atlas (TCGA)-HNSCC datasets. The genetic and epigenetic alterations, druggable genes, their associated functional pathways and their possible association with cancer hallmarks were identified. We have performed extensive analysis to identify the target genes of DEHGs driving HNSCC. The differentially expressed HOX cluster-embedded microRNAs (DEHMs) in HNSCC and their association with HOX-target genes were evaluated to construct a regulatory network of the HOX cluster in HNSCC. Our analysis identified sixteen DEHGs in HNSCC and determined their importance in stage stratification and HPV infection. We found a total of 55 HNSCC driver genes that were identified as targets of DEHGs. The involvement of DEHGs and their targets in cancer-associated signaling mechanisms have confirmed their role in pathophysiology. Further, we found that their oncogenic nature could be targeted by using the novel and approved anti-neoplastic drugs in HNSCC. Construction of the regulatory network depicted the interaction between DEHGs, DEHMs and their targets genes in HNSCC. Hence, aberrantly expressed HOX cluster genes function in a coordinated manner to drive HNSCC. It could provide a broad perspective to carry out the experimental investigation, to understand the underlying oncogenic mechanism and allow the discovery of new clinical biomarkers for HNSCC.

A ceRNA Network Composed of Survival-Related lncRNAs, miRNAs, and mRNAs in Clear Cell Renal Carcinoma

Clear cell renal carcinoma (ccRCC) is one of the most common renal carcinomas worldwide, which has worse prognosis compared with other subtypes of tumors. We propose a potential RNA regulatory mechanism associated with ccRCC progression. Accordingly, we screened out clinical factors and the expression of RNAs and miRNAs of ccRCC from the TCGA database. 9 lncRNAs (FGF12-AS2, WT1-AS, TRIM36-IT1, AC009093.1, LINC00443, TCL6, COL18A1-AS1, AC110619.1, HOTTIP), 2 miRNAs (mir-155 and mir-21), and 3 mRNAs (COL4A4, ERMP1, PRELID2) were selected from differential expression RNAs and built predictive survival models. The survival models performed very well in predicting prognosis and were found to be highly correlated with tumor stage. In addition, the survival-related lncRNA-miRNA-mRNA (ceRNA) network was constructed by 18 RNAs including 12 mRNAs, 2 miRNAs, and 4 lncRNAs. It is found that the “ECM-receptor interaction,” “Pathways in cancer,” and “Chemokine signaling pathway” as the main pathways in KEGG pathway analysis. Overall, we established predictive survival model and ceRNA network based on multivariate Cox regression analysis. It may open a new approach and potential biomarkers for clinical prognosis and treatment of ccRCC patients.

miRCoop: Identifying Cooperating miRNAs via Kernel Based Interaction Tests

Although miRNAs can cause widespread changes in expression programs, single miRNAs typically induce mild repression on their targets. Cooperativity among miRNAs is reported as one strategy to overcome this constraint. Expanding the catalog of synergistic miRNAs is critical for understanding gene regulation and for developing miRNA-based therapeutics. In this study, we develop miRCoop to identify synergistic miRNA pairs that have weak or no repression on the target mRNA individually, but when act together, induce strong repression. miRCoop uses kernel-based statistical interaction tests, together with miRNA and mRNA target information. We apply our approach to patient data of two different cancer types. In kidney cancer, we identify 66 putative triplets. For 64 of these triplets, there is at least one common transcription factor that potentially regulates all participating RNAs of the triplet, supporting a functional association among them. Furthermore, we find that identified triplets are enriched for certain biological processes that are relevant to kidney cancer. Some of the synergistic miRNAs are very closely encoded in the genome, hinting a functional association among them. In applying the method on tumor data with the primary liver site, we find 3105 potential triplet interactions. We believe miRCoop can aid our understanding of the complex regulatory interactions in different health and disease states of the cell and can help in designing miRNA-based therapies. Matlab code for the methodology is provided in https://github.com/guldenolgun/miRCoop.

Changes in adipose tissue microRNA expression across the menstrual cycle in regularly menstruating females: a pilot study

Cyclical changes in hormone profiles across the menstrual cycle are associated with alterations in metabolic control. MicroRNAs (miRNAs) contribute to regulating metabolic control, including adipose tissue metabolism. How fluctuations in hormonal profiles across the menstrual cycle affect adipose tissue miRNA expression remains unknown. Eleven healthy, regularly menstruating females underwent four sampling visits across their menstrual cycle. Subcutaneous abdominal adipose tissue and venous blood samples were collected at each sampling visit. Luteinizing hormone (LH) tests, calendar counting, and serum hormone concentrations were used to determine menstrual cycle phases: early-follicular (EF), late-follicular (LF), postovulatory (PO), and midluteal (ML). Serum follicle-stimulating hormone, LH, estrogen, progesterone, and testosterone were determined using multiplex magnetic bead panels and enzyme-linked immunosorbent assays. Global adipose tissue miRNA expression levels were determined via microarray in a subset of participants (n = 8) and 17 candidate miRNAs were validated by RT-qPCR in the whole cohort (n = 11). Global analysis of adipose tissue miRNA expression identified 33 miRNAs significantly altered across the menstrual cycle; however, no significant differences remained after correcting for multiple testing (P > 0.05). RT-qPCR analysis of candidate miRNAs revealed miR-497-5p expression was significantly altered across the menstrual cycle ([Formula: see text] = 0.18, P = 0.03); however, post hoc tests did not reveal any significant differences between menstrual cycle phases (P > 0.05). miR-30c-5p was associated with testosterone concentration (R² = 0.13, P = 0.033). These pilot data indicate differences in adipose tissue miRNAs in healthy women across the menstrual cycle and a weak association with ovarian hormones. Further research in larger sample sizes is required to confirm regulation of miRNA expression across the menstrual cycle.

MicroRNA Targeting

MicroRNAs (miRNAs) are small noncoding elements that play essential roles in the posttranscriptional regulation of biochemical processes. miRNAs recognize and target multiple mRNAs; therefore, investigating miRNA dysregulation is an indispensable strategy to understand pathological conditions and to design innovative drugs. Targeting miRNAs in diseases improve outcomes of several therapeutic strategies thus, this present study highlights miRNA targeting methods through experimental assays and bioinformatics tools. The first part of this review focuses on experimental miRNA targeting approaches for elucidating key biochemical pathways. A growing body of evidence about the miRNA world reveals the fact that it is not possible to uncover these molecules' structural and functional characteristics related to the biological processes with a deterministic approach. Instead, a systemic point of view is needed to truly understand the facts behind the natural complexity of interactions and regulations that miRNA regulations present. This task heavily depends both on computational and experimental capabilities. Fortunately, several miRNA bioinformatics tools catering to nonexperts are available as complementary wet-lab approaches. For this purpose, this work provides recent research and information about computational tools for miRNA targeting research.

The Performance Comparison of Gene Co-expression Networks of Breast and Prostate Cancer using Different Selection Criteria

Gene co-expression networks (GCN) present undirected relations between genes to understand molecular structures behind the diseases, including cancer. The utilization of various biological datasets and gene network inference (GNI) algorithms can reveal meaningful gene-gene interactions of GCNs. This study applies three GNI algorithms on mRNA gene expression, RNA-Seq, and miRNA-target genes datasets to infer GCNs of breast and prostate cancers. To evaluate the performance of the GCNs, we utilize overlap analysis via literature data, topological assessment, and Gene Ontology-based biological assessment. The results emphasize how the selection of biological datasets and GNI algorithms affect the performance results on different evaluation criteria. GCNs on microarray gene expression data slightly outperform in overlap analysis. Also, GCNs on RNA-Seq and gene expression datasets follow scale-free topology. The biological assessment results are close to each other on all biological datasets. C3NET algorithm-based GCNs did not contain any biological assessment modules; therefore, it is not optimal for biological assessment. GNI algorithms' selection did not change the overlap analysis and topological assessment results. Our primary objective is to compare the performance results of biological datasets and GNI algorithms based on different evaluation criteria. For this purpose, we developed the GNIAP R package that enables users to select different GNI algorithms to infer GCNs. The GNIAP R package also provides literature-based overlap analysis, and topological and biological analyses on GCNs. Users can access the GNIAP R package via https://github.com/ozgurcingiz/GNIAP .

Emerging Evidence of the Functional Impact of the miR379/miR656 Cluster (C14MC) in Breast Cancer

Many microRNAs exist in clusters that share comparable sequence homology and may target genes in a common pathway. The miR-379/miR-656 (C14MC) cluster is imprinted in the DLK1-Dio3 region of 14q32.3 and contains 42 miRNAs. It plays a functional role in numerous biological pathways including vascular remodeling and early development. With many C14MC miRNAs highlighted as potential tumor suppressors in a variety of cancers, the role of this cluster in breast cancer (BC) has garnered increased attention in recent years. This review focuses on C14MC in BC, providing an overview of the constituent miRNAs and addressing each in terms of functional impact, potential target genes/pathways, and, where relevant, biomarker capacity. Studies have revealed the regulation of key factors in disease progression and metastasis including tyrosine kinase pathways and factors critical to epithelial–mesenchymal transition (EMT). This has potentially important clinical implications, with EMT playing a critical role in BC metastasis and tyrosine kinase inhibitors (TKIs) in widespread use for the treatment of BC. While the majority of studies have reported tumor-suppressing roles for these miRNAs, some have highlighted their potential as oncomiRs. Understanding the collective contribution of miRNAs within C14MC to BC may support improved understanding of disease etiology and present novel approaches to targeted therapy.

Long noncoding RNAs in cancer metastasis

Metastasis is a major contributor to cancer-associated deaths. It is characterized by a multistep process that occurs through the acquisition of molecular and phenotypic changes enabling cancer cells from a primary tumour to disseminate and colonize at distant organ sites. Over the past decade, the discovery and characterization of long noncoding RNAs (lncRNAs) have revealed the diversity of their regulatory roles, including key contributions throughout the metastatic cascade. Here, we review how lncRNAs promote metastasis by functioning in discrete pro-metastatic steps including the epithelial-mesenchymal transition, invasion and migration and organotrophic colonization, and by influencing the metastatic tumour microenvironment, often by interacting within ribonucleoprotein complexes or directly with other nucleic acid entities. We discuss well-characterized lncRNAs with in vivo phenotypes and highlight mechanistic commonalities such as convergence with the TGFβ-ZEB1/ZEB2 axis or the nuclear factor-κB pathway, in addition to lncRNAs with controversial mechanisms and the influence of methodologies on mechanistic interpretation. Furthermore, some lncRNAs can help identify tumours with increased metastatic risk and spur novel therapeutic strategies, with several lncRNAs having shown potential as novel targets for antisense oligonucleotide therapy in animal models. In addition to well-characterized examples of lncRNAs functioning in metastasis, we discuss controversies and ongoing challenges in lncRNA biology. Finally, we present areas for future study for this rapidly evolving field.

Self-Organization and Information Processing: From Basic Enzymatic Activities to Complex Adaptive Cellular Behavior

One of the main aims of current biology is to understand the origin of the molecular organization that underlies the complex dynamic architecture of cellular life. Here, we present an overview of the main sources of biomolecular order and complexity spanning from the most elementary levels of molecular activity to the emergence of cellular systemic behaviors. First, we have addressed the dissipative self-organization, the principal source of molecular order in the cell. Intensive studies over the last four decades have demonstrated that self-organization is central to understand enzyme activity under cellular conditions, functional coordination between enzymatic reactions, the emergence of dissipative metabolic networks (DMN), and molecular rhythms. The second fundamental source of order is molecular information processing. Studies on effective connectivity based on transfer entropy (TE) have made possible the quantification in bits of biomolecular information flows in DMN. This information processing enables efficient self-regulatory control of metabolism. As a consequence of both main sources of order, systemic functional structures emerge in the cell; in fact, quantitative analyses with DMN have revealed that the basic units of life display a global enzymatic structure that seems to be an essential characteristic of the systemic functional metabolism. This global metabolic structure has been verified experimentally in both prokaryotic and eukaryotic cells. Here, we also discuss how the study of systemic DMN, using Artificial Intelligence and advanced tools of Statistic Mechanics, has shown the emergence of Hopfield-like dynamics characterized by exhibiting associative memory. We have recently confirmed this thesis by testing associative conditioning behavior in individual amoeba cells. In these Pavlovian-like experiments, several hundreds of cells could learn new systemic migratory behaviors and remember them over long periods relative to their cell cycle, forgetting them later. Such associative process seems to correspond to an epigenetic memory. The cellular capacity of learning new adaptive systemic behaviors represents a fundamental evolutionary mechanism for cell adaptation.

MicroRNAs in ascending thoracic aortic aneurysms

Thoracic Aortic Aneurysm (TAA) is characterized by the dilation of the aorta and is fatal if not diagnosed and treated appropriately. The underlying genetic mechanisms have not been completely delineated, so better knowledge of the physiopathology of TAAs is needed to improve detection and therapy. MicroRNAs (miRNAs) regulate gene expression post-transcriptionally and are known to be involved in cardiovascular diseases (CVDs). The current study aimed to identify miRNAs that can be used as possible biomarkers for the early diagnosis of patients with ascending TAAs (ATAAs). MiRNA expression was profiled by NanoString nCounter technology using 12 samples including tissue and pre- and post-surgical plasma from ATAA patients. Four miRNAs were selected and further validated by real time polymerase chain reaction (RT-PCR) in 22 plasma samples from which three miRNAs (hsa-miR140-5p, hsa-miR-191-5p and hsa-miR-214-3p) showed significant expression level differences between the two types of plasma samples. Further analyses of the corresponding predicted target genes by these miRNAs, revealed two genes (Myotubularin-related protein 4 (MTMR4) and Phosphatase 1 catalytic subunit β (PPP1CB)) whose expression was inversely correlated with the expression of their respective miRNAs. Overall, in this pilot study, we identified three miRNAs that might serve as potential biomarkers and therapeutic targets in ATAA.

miRGTF-net: Integrative miRNA-gene-TF network analysis reveals key drivers of breast cancer recurrence

Analysis of regulatory networks is a powerful framework for identification and quantification of intracellular interactions. We introduce miRGTF-net, a novel tool for construction of miRNA-gene-TF networks. We consider multiple transcriptional and post-transcriptional interaction types, including regulation of gene and miRNA expression by transcription factors, gene silencing by miRNAs, and co-expression of host genes with their intronic miRNAs. The underlying algorithm uses information on experimentally validated interactions as well as integrative miRNA/mRNA expression profiles in a given set of samples. The latter ensures simultaneous tissue-specificity and biological validity of interactions. We applied miRGTF-net to paired miRNA/mRNA-sequencing data of breast cancer samples from The Cancer Genome Atlas (TCGA). Together with topological analysis of the constructed network we showed that considered players can form reliable prognostic gene signatures for ER-positive breast cancer. A number of signatures demonstrated remarkably high accuracy on transcriptomic data obtained by both microarrays and RNA sequencing from several independent patient cohorts. Furthermore, an essential part of prognostic genes were identified as direct targets of transcription factor E2F1. The putative interplay between estrogen receptor alpha and E2F1 was suggested as a potential recurrence factor in patients treated with tamoxifen. Source codes of miRGTF-net are available at GitHub (https://github.com/s-a-nersisyan/miRGTF-net).

HOTAIR inhibits the proliferation of glioblastoma cells by targeting miR-219

BACKGROUND

To investigate the influences of HOX transcript antisense ribonucleic acid (HOTAIR) on the proliferation and apoptosis of glioblastoma cells by targeting micro RNA (miR)-219.

OBJECTIVE

With glioblastoma cell line U87 as the object, the changes in expression levels of HOTAIR and miR-219 in each group were detected via quantitative real-time polymerase chain reaction (qRT-PCR) after HOTAIR in U87 cell lines was knocked down using small interfering RNA (siRNA). Then the cell proliferation in each group was determined using Cell Counting Kit-8 (CCK-8) and colony formation assays. Flow cytometry was applied to detect the cell apoptosis, and Western blotting assay was adopted to measure the changes in protein levels of Cyclin D1 and Bcl-2-associated X protein (Bax). After miR-219 knockdown with siRNA, the changes in expression levels of HOTAIR and miR-219 in each group were examined through qRT-PCR, and the cell proliferation was tested by CCK-8 assay.

RESULTS

After interference inHOTAIR using siRNA, compared with those in control group, the RNA expression level of HOTAIR was decreased remarkably (p< 0.05), the RNA expression level of miR-219 was increased notably (p< 0.05), the cell proliferation rate was inhibited evidently (p< 0.05), the apoptosis rate was enhanced obviously (p< 0.05), the protein expression level of Cyclin D1 declined markedly (p< 0.05), and the protein expression level of Bax rose distinctly (p< 0.05) in HOTAIR-siRNA group. After miR-219 knockdown with siRNA, the cell proliferation rate was raised remarkably (p< 0.05), but there was no significant change in the RNA expression level of HOTAIR (p> 0.05).

CONCLUSION

HOTAIR can repress the proliferation and promote the apoptosis of glioblastoma cells by targeting miR-219.

Integrated bioinformatic analysis of miR-15a/16-1 cluster network in cervical cancer

The miR-15a/16-1 cluster is abnormally expressed in cervical cancer (CC) tissues and plays a vital role in cervical carcinogenesis. We aimed to evaluate the miR-15a/16-1 expression in healthy and cancerous cervical tissues, identify the associated networks, and to test its prognostic significance. miR-15a/16-1-MC expressions were analyzed in TCGA-CESC datasets by UALCAN, GEPIA2, and Datasetviewer. miR-15a/16-1 validated targets were extracted from mirTarBase and in silico functional analysis of the target genes were performed using WebGestalt. The interaction networks were constructed by the miRNet, STRING, and NetworkAnalyst tools. The prognostic significance and metastatic potential of the target genes were predicted using UALCAN and HCMDB. The FDA approved drugs to target miR-15a/16-1 and target gene network in CC were performed using DGIdb, STITCH and PanDrugs. TCGA-CESC and GEO data analysis suggested significant overexpression of miR-15a/16-1 in CC samples. The Kaplan-Meier survival analysis showed that miR-15a and its four target genes (BCL2, CCNE1, NUP50, and RBPJ) influence the overall survival of CC patients. Among the 66 differentially expressed target genes, 12 of them are linked to head, neck, or lung metastasis. Functional enrichment analysis predicted the association of this cluster with p53 signaling, human papillomavirus infection, PI3-AKT signaling pathway, and pathways in cancer. Drug-gene interaction analysis showed 52 potential FDA approved drugs to interact with the miR-15a/16-1 target genes. Nine of the 52 drugs are currently used as a chemotherapeutic agent for the treatment of CC patients. The present study shows that miR-15a/16-1 expression can be used as a clinical marker and target for therapy in CC.

Identification of miR-802-5p and its involvement in type 2 diabetes mellitus

MicroRNAs (miRNA) are recently discovered endogenous, small noncoding RNAs (of 22 nucleotides) that play pivotal roles in gene regulation. They are involved in post-transcriptional control of gene expression. miRNAs are emerging as important regulators of cell proliferation, development, cancer formation, stress responses, cell death and physiological conditions. Increasing evidence has demonstrated the human miRNAs bind to their target mRNA sequences with perfect or near-perfect sequence complementarily. This provides a powerful strategy for discovering potential type 2 diabetes mellitus (T2DM) targets and gives the probability to exploit them for diagnostic and therapeutic causes. About 6% of the world population is affected by T2DM, and it is recognized as a global epidemic by the World Health Organization. At present there is no valid biomarker to control or manage T2DM. Therefore, the present study applied a mature sequence of miRNAs from publicly accessible databases to identify the miRNA from T2DM expressed sequence tags, and the results are detailed and discussed below.

miR-631 Inhibits Intrahepatic Metastasis of Hepatocellular Carcinoma by Targeting PTPRE

MicroRNAs (miRNAs) have been reported to play critical roles in the pathological development of hepatocellular carcinoma (HCC), one of the most common cancers in the world. Our study aims to explore the expression, function and mechanism of miR-631 in HCC. Our findings are that expression of miR-631 is significantly down-regulated in HCC tissue compared with that in adjacent non-cancerous tissue, and low expression of miR-631 in HCC tissue is associated with cirrhosis, multiple tumors, incomplete tumor encapsulation, poor tumor differentiation, and high TNM stage. Our test results showed that miR-631 could inhibit migration, invasion, epithelial–mesenchymal transition (EMT) and intrahepatic metastasis of HCC. Receptor-type protein tyrosine phosphatase epsilon (PTPRE) as a downstream target of miR-631 could promote migration, invasion and EMT of HCC cells. Besides, the expression of PTPRE had a negative correlation with the expression of miR-631 both in vivo and in vitro, and increasing expression of PTPRE could reverse inhibitory effects of miR-631 in HCC cells. In sum, our study first demonstrated that miR-631 targeted PTPRE to inhibit intrahepatic metastasis in HCC. We gain insights from these findings into the mechanism of miRNAs regulation in HCC metastasis and further introduce a novel therapeutic target for HCC treatment.

A compound downregulation of SRRM2 and miR-27a-3p with upregulation of miR-27b-3p in PBMCs of Parkinson's patients is associated with the early stage onset of disease

Parkinson’s disease (PD) is diagnosed when motor symptoms emerges, which almost 70% of dopamine neurons are lost. Therefore, early diagnosis of PD is crucial to prevent the progress of disease. Blood-based biomarkers, which are minimally invasive, potentially used for diagnosis of PD, including miRNAs. The aim of this study was to assess whether SRRM2 and miR-27a/b-3p could act as early diagnostic biomarkers for PD. Total RNAs from PBMCs of 30 PD’s patients and 14 healthy age and gender matched subjects was extracted. The expression levels of respective genes were assessed. Data were presented applying a two-tailed unpaired t-test and one-way ANOVA. We observed significant down-regulation of SRRM2 (p = 0.0002) and miR-27a-3p (p = 0.0001), and up-regulation of miR-27b-3p (p = 0.02) in PBMCs of Parkinson's patients. Down-regulation of miR-27a-3p is associated with increasing disease severity, whereas the up-regulation of miR-27b-3p was observed mostly at HY-1 and disease duration between 3–5 years. There was a negative correlation between SRRM2 and miR-27b-3p expressions, and miR-27a-3p positively was correlated with miR-27b-3p. Based on functional enrichment analysis, SRRM2 and miR-27a/b-3p acted on common functional pathways. miR-27a/b-3p could potentially predict the progression and severity of PD. Although both miRs had no similarity on expression, a positive correlation between both miRs was identified, supporting their potential role as biomarkers in clinical PD stages. Of note that SRRM2 and miR-27a-3p were able to distinguish PD patients from healthy individuals. Functional analysis of the similarity between genes associated with SRRM2 and miR-27a/b-3p indicates common functional pathways and their dysfunction correlates with molecular etiopathology mechanisms of PD onset.

Circulating miRNAs in HER2-Positive and Triple Negative Breast Cancers: Potential Biomarkers and Therapeutic Targets

Breast cancer is one of the most prevalent diseases among women worldwide and is highly associated with cancer-related mortality. Of the four major molecular subtypes, HER2-positive and triple-negative breast cancer (TNBC) comprise more than 30% of all breast cancers. While the HER2-positive subtype lacks estrogen and progesterone receptors and overexpresses HER2, the TNBC subtype lacks estrogen, progesterone and HER2 receptors. Although advances in molecular biology and genetics have substantially ameliorated breast cancer disease management, targeted therapies for the treatment of estrogen-receptor negative breast cancer patients are still restricted, particularly for TNBC. On the other hand, it has been demonstrated that microRNAs, miRNAs or small non-coding RNAs that regulate gene expression are involved in diverse biological processes, including carcinogenesis. Moreover, circulating miRNAs in serum/plasma are among the most promising diagnostic/therapeutic tools as they are stable and relatively easy to quantify. Various circulating miRNAs have been identified in several human cancers including specific breast cancer subtypes. This review aims to discuss the role of circulating miRNAs as potential diagnostic and prognostic biomarkers as well as therapeutic targets for estrogen-receptor negative breast cancers, HER2+ and triple negative.

miR760 regulates ATXN1 levels via interaction with its 5' untranslated region

Identifying modifiers of dosage-sensitive genes involved in neurodegenerative disorders is imperative to discover novel genetic risk factors and potential therapeutic entry points. In this study, we focus on Ataxin-1 (ATXN1), a dosage-sensitive gene involved in the neurodegenerative disease spinocerebellar ataxia type 1 (SCA1). While the precise maintenance of ATXN1 levels is essential to prevent disease, the mechanisms that regulate ATXN1 expression remain largely unknown. We demonstrate that ATXN1's unusually long 5' untranslated region (5' UTR) negatively regulates its expression via posttranscriptional mechanisms. Based on recent reports that microRNAs (miRNAs) can interact with both 3' and 5' UTRs to regulate their target genes, we identify miR760 as a negative regulator that binds to a conserved site in ATXN1's 5' UTR to induce RNA degradation and translational inhibition. We found that delivery of Adeno-associated virus (AAV)-expressing miR760 in the cerebellum reduces ATXN1 levels in vivo and mitigates motor coordination deficits in a mouse model of SCA1. These findings provide new insights into the regulation of ATXN1 levels, present additional evidence for miRNA-mediated gene regulation via 5' UTR binding, and raise the possibility that noncoding mutations in the ATXN1 locus may act as risk factors for yet to be discovered progressive ataxias.

The Similar Effects of miR-512-3p and miR-519a-2-5p on the Promotion of Hepatocellular Carcinoma: Different Tunes Sung With Equal Skill

Although the therapeutic methods of hepatocellular carcinoma (HCC) have made great advances, the current situation is that HCC is the common malignancy. Our previous bioinformatic study presented that two members of C19MC (mir-512-1 and mir-519a-2) acted as crucial roles in the HCC progression. In this study, we first demonstrated that the miR-512-3p and miR-519a-2-5p, which were spliced from the mir-512-1 and mir-519a-2, were the functional mature miRNAs. Meanwhile, both miR-512-3p and miR-519a-2-5p were significantly upregulated in human HCC samples and HCC cell lines. The miR-512-3p and miR-519a-2-5p promoted the proliferation, invasion, and metastasis in vitro and in vivo. Moreover, the two miRNAs co-targeted the downstream tumor suppressors MAP3K2 and MAP2K4 and subsequently achieved the HCC progression. In the clinical cohort, high expression of miR-512-3p and miR-519a-2-5p acted as two risk factors for HCC recurrence and distinguished patients with poor tumor-free survival after radical resection. The integration of the two miRNAs into the AJCC staging system significantly improved the accuracy for the prediction of HCC recurrence. Our study suggests that miR-512-3p and miR-519a-2-5p have similar effects on the promotion of HCC progression. They can be robust markers for the prediction of HCC recurrence and therapy targets.

Dynamical gene regulatory networks are tuned by transcriptional autoregulation with microRNA feedback

Concepts from dynamical systems theory, including multi-stability, oscillations, robustness and stochasticity, are critical for understanding gene regulation during cell fate decisions, inflammation and stem cell heterogeneity. However, the prevalence of the structures within gene networks that drive these dynamical behaviours, such as autoregulation or feedback by microRNAs, is unknown. We integrate transcription factor binding site (TFBS) and microRNA target data to generate a gene interaction network across 28 human tissues. This network was analysed for motifs capable of driving dynamical gene expression, including oscillations. Identified autoregulatory motifs involve 56% of transcription factors (TFs) studied. TFs that autoregulate have more interactions with microRNAs than non-autoregulatory genes and 89% of autoregulatory TFs were found in dual feedback motifs with a microRNA. Both autoregulatory and dual feedback motifs were enriched in the network. TFs that autoregulate were highly conserved between tissues. Dual feedback motifs with microRNAs were also conserved between tissues, but less so, and TFs regulate different combinations of microRNAs in a tissue-dependent manner. The study of these motifs highlights ever more genes that have complex regulatory dynamics. These data provide a resource for the identification of TFs which regulate the dynamical properties of human gene expression.

RPmirDIP: Reciprocal Perspective improves miRNA targeting prediction

MicroRNAs (miRNAs) are short, non-coding RNAs that interact with messenger RNA (mRNA) to accomplish critical cellular activities such as the regulation of gene expression. Several machine learning methods have been developed to improve classification accuracy and reduce validation costs by predicting which miRNA will target which gene. Application of these predictors to large numbers of unique miRNA–gene pairs has resulted in datasets comprising tens of millions of scored interactions; the largest among these is mirDIP. We here demonstrate that miRNA target prediction can be significantly improved (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p < 0.001$$\end{document}p<0.001) through the application of the Reciprocal Perspective (RP) method, a cascaded, semi-supervised machine learning method originally developed for protein-protein interaction prediction. The RP method, aptly named RPmirDIP, augments the original mirDIP prediction scores by leveraging local thresholds from the two complimentary views available to each miRNA–gene pair, rather than apply a traditional global decision threshold. Application of this novel RPmirDIP predictor promises to help identify new, unexpected miRNA–gene interactions. A dataset of RPmirDIP-scored interactions are made available to the scientific community at cu-bic.ca/RPmirDIP and 10.5683/SP2/LD8JKJ.

MMiRNA-Viewer2, a bioinformatics tool for visualizing functional annotation for MiRNA and MRNA pairs in a network

Background

Although there are many studies on the characteristics of miRNA-mRNA interactions using miRNA and mRNA sequencing data, the complexity of the change of the correlation coefficients and expression values of the miRNA-mRNA pairs between tumor and normal samples is still not resolved, and this hinders the potential clinical applications. There is an urgent need to develop innovative methodologies and tools that can characterize and visualize functional consequences of cancer risk gene and miRNA pairs while analyzing the tumor and normal samples simultaneously.

Results

We developed an innovative bioinformatics tool for visualizing functional annotation of miRNA-mRNA pairs in a network, known as MMiRNA-Viewer². The tool takes mRNA and miRNA interaction pairs and visualizes mRNA and miRNA regulation network. Moreover, our MMiRNA-Viewer² web server integrates and displays the mRNA and miRNA gene annotation information, signaling cascade pathways and direct cancer association between miRNAs and mRNAs. Functional annotation and gene regulatory information can be directly retrieved from our web server, which can help users quickly identify significant interaction sub-network and report possible disease or cancer association. The tool can identify pivotal miRNAs or mRNAs that contribute to the complexity of cancer, while engaging modern next-generation sequencing technology to analyze the tumor and normal samples concurrently. We compared our tools with other visualization tools.

Conclusion

Our MMiRNA-Viewer² serves as a multitasking platform in which users can identify significant interaction clusters and retrieve functional and cancer-associated information for miRNA-mRNA pairs between tumor and normal samples. Our tool is applicable across a range of diseases and cancers and has advantages over existing tools.

Expression profiling and analysis of some miRNAs in subcutaneous white adipose tissue during development of obesity

Background

MicroRNAs are emerging as new mediators in the regulation of adipocyte physiology and have been approved to play a role in obesity. Despite several studies have focused on microRNA expression profiles and functions in different metabolic tissues, little is known about their response to nutritional interventions in white adipose tissue during obesity stages, and whether they differ in this response to weight-reduction strategy is poorly understood. Our objectives were to study the dysregulation of some miRNAs in subcutaneous inguinal white adipose tissue during weight change, expansion/reduction; in response to both a high-fat diet and switching to a normal diet feeding, and to evaluate them as potential biomarkers and therapeutic targets for early obesity management

Method

A hundred 6-week-old male Wister rats were randomly divided into a normal diet group (N.D), a high-fat diet group (H.F.D), and a switched to a normal diet group (H.F.D/N.D). At the beginning and at intervals 2 weeks, serum lipid, hormone levels, total body fat mass, and inguinal subcutaneous white adipose tissue mass (WAT) measurements were recorded using dual-energy X-ray absorptiometry (DEXA). The expression levels of microRNAs were evaluated using real-time PCR.

Results

Significant alterations were observed in serum glucose, lipid profile, and adipokine hormones during the early stages of obesity development. Alteration in rno-mir 30a-5p, rno-mir 133a-5p, and rno-mir 107-5p expression levels were observed at more than one time point. While rno-let-7a-5p, rno-mir 193a-5p, and rno-mir125a-5p were downregulated and rno-mir130a-5p was upregulated at all time points within 2 to 4 weeks in response to H.F.D feeding for 10 weeks. The impact of switching to normal diet has a reversed effect on lipid profile, adipokine hormone levels, and some miRNAs. The bioinformatics results have identified a novel and important pathway related to inflammatory signalling.

Conclusion

Our research demonstrated significant alterations in some adipocyte-expressed miRNAs after a short time of high caloric diet consumption. This provides further evidence of the significant role of nutrition as an epigenetic factor in regulation of lipid and glucose metabolism genes by modulating of related key miRNAs. Therefore, we suggest that miRNAs could be used as biomarkers for adiposity during diet-induced obesity. Perhaps limitation in calories intake is a way to manipulate obesity and associated metabolic disorders. Further studies are needed to fully elucidate the role of microRNAs in the development of obesity

Unravelling the role of long non-coding RNA - LINC01087 in breast cancer

Apoptosis is a 'programmed fate' of all cells participating in diverse physiological and pathological conditions. The role of critical regulators and their involvement in this complex multi-stage process of apoptosis weaved around non-coding RNAs (ncRNAs) is poorly deciphered in breast carcinoma (BC). Aberrant expression patterns of the ncRNAs and their interacting partners, either ncRNAs or coding RNAs or proteins at any point along these pathways, may lead to the malignant transformation of the affected cells, tumour metastasis and resistance to anticancer drugs. Longest non-coding type of ncRNAs (lncRNAs) have been considered as critical factors for the development and progression of breast cancer. The aim of our study was to identify set of novel lncRNAs interacting with microRNAs (miRNAs) or proteins that were significantly dysregulated in breast cancer using RNA-Sequencing (RNA-Seq) technique in different samples acting as oncogenic drivers contributing to cancerous phenotype involved in post-transcriptional processing of RNAs. Four lncRNAs; LINC01087, lnc-CLSTN2-1:1, lnc-c7orf65-3:3 and LINC01559:2 were selected for further analysis. Gene expression analysis of over-expressed LINC01087 in vitro reduced both cell viability and apoptosis. We integrated miRNA and mRNA (hsa-miR-548 and AKT1) expression profiles with curated regulations with lncRNA (LINC01087) which has not been previously associated with any breast cancer type, using different computational tools. The network (lncRNA→ miRNA→ mRNA) is promising for the identification of carcinoma associated genes and apoptosis signaling path highlighting the potential roles of LINC01087, hsa-miR548n, AKT1 gene which may play crucial role in proliferation.

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

MOTIVATION

Since the initial discovery of microRNAs as post-transcriptional, regulatory key players in the 1990s, a total number of $2656$ mature microRNAs have been publicly described for Homo sapiens. As discovery of new miRNAs is still on-going, target identification remains to be an essential and challenging step preceding functional annotation analysis. One key challenge for researchers seems to be the selection of the most appropriate tool out of the larger multiverse of published solutions for a given research study set-up.

RESULTS

In this review we collectively describe the field of in silico target prediction in the course of time and point out long withstanding principles as well as recent developments. By compiling a catalog of characteristics about the 98 prediction methods and identifying common and exclusive traits, we signpost a simplified mechanism to address the problem of application selection. Going further we devised interpretation strategies for common types of output as generated by frequently used computational methods. To this end, our work specifically aims to make prospective users aware of common mistakes and practical questions that arise during the application of target prediction tools.

AVAILABILITY

An interactive implementation of our recommendations including materials shown in the manuscript is freely available at https://www.ccb.uni-saarland.de/mtguide.

microRNAs Regulating Human and Mouse Naïve Pluripotency

microRNAs are ~22bp nucleotide non-coding RNAs that play important roles in the post-transcriptional regulation of gene expression. Many studies have established that microRNAs are important for cell fate choices, including the naïve to primed pluripotency state transitions, and their intermediate state, the developmentally suspended diapause state in early development. However, the full extent of microRNAs associated with these stage transitions in human and mouse remain under-explored. By meta-analysis of microRNA-seq, RNA-seq, and metabolomics datasets from human and mouse, we found a set of microRNAs, and importantly, their experimentally validated target genes that show consistent changes in naïve to primed transitions (microRNA up, target genes down, or vice versa). The targets of these microRNAs regulate developmental pathways (e.g., the Hedgehog-pathway), primary cilium, and remodeling of metabolic processes (oxidative phosphorylation, fatty acid metabolism, and amino acid transport) during the transition. Importantly, we identified 115 microRNAs that significantly change in the same direction in naïve to primed transitions in both human and mouse, many of which are novel candidate regulators of pluripotency. Furthermore, we identified 38 microRNAs and 274 target genes that may be involved in diapause, where embryonic development is temporarily suspended prior to implantation to uterus. The upregulated target genes suggest that microRNAs activate stress response in the diapause stage. In conclusion, we provide a comprehensive resource of microRNAs and their target genes involved in naïve to primed transition and in the paused intermediate, the embryonic diapause stage.

Biclustering analysis of transcriptome big data identifies condition-specific microRNA targets

We present a novel approach to identify human microRNA (miRNA) regulatory modules (mRNA targets and relevant cell conditions) by biclustering a large collection of mRNA fold-change data for sequence-specific targets. Bicluster targets were assessed using validated messenger RNA (mRNA) targets and exhibited on an average 17.0% (median 19.4%) improved gain in certainty (sensitivity + specificity). The net gain was further increased up to 32.0% (median 33.4%) by incorporating functional networks of targets. We analyzed cancer-specific biclusters and found that the PI3K/Akt signaling pathway is strongly enriched with targets of a few miRNAs in breast cancer and diffuse large B-cell lymphoma. Indeed, five independent prognostic miRNAs were identified, and repression of bicluster targets and pathway activity by miR-29 was experimentally validated. In total, 29 898 biclusters for 459 human miRNAs were collected in the BiMIR database where biclusters are searchable for miRNAs, tissues, diseases, keywords and target genes.

Cotargeting among microRNAs in the brain

MicroRNAs (miRNAs) play roles in diverse developmental and disease processes. Distinct miRNAs have hundreds to thousands of conserved mRNA binding sites but typically direct only modest repression via single sites. Cotargeting of individual mRNAs by different miRNAs could potentially achieve stronger and more complex patterns of repression. By comparing target sets of different miRNAs, we identified hundreds of pairs of miRNAs that share more mRNA targets than expected (often by twofold or more) relative to stringent controls. Genetic perturbations revealed a functional overlap in neuronal differentiation for the cotargeting pair miR-138/miR-137. Clustering of all cotargeting pairs revealed a group of nine predominantly brain-enriched miRNAs that share many targets. In reporter assays, subsets of these miRNAs together repressed gene expression by five- to 10-fold, often showing cooperative repression. Together, our results uncover an unexpected pattern in which combinations of miRNAs collaborate to robustly repress cotargets, and suggest important developmental roles for cotargeting.

miR-190-5p in human diseases

miRNAs, a major class of small noncoding RNAs approximately 18-25 nucleotides in length, function by repressing the expression of target genes through binding to complementary sequences in the 3'-UTRs of target genes. Emerging evidence has highlighted their important roles in numerous diseases, including human cancers. Recently, miR-190 has been shown to be dysregulated in various types of human cancers that participates in cancer-related biological processes, including proliferation, apoptosis, metastasis, drug resistance, by regulating associated target genes, and to predict cancer diagnosis and prognosis. In this review, we summarized the roles of miR-190-5p in human diseases, especially in human cancers. Then we classified its target genes in tumorigenesis and progression, which might provide evidence for cancer diagnosis and prognosis, promising tools for cancer treatment, or leads for further investigation.

miRNA-dependent regulation of STIM1 expression in breast cancer

Store-operated Ca²⁺ entry (SOCE) has been shown to be important for breast cancer metastasis in xenograft mouse models. The ER Ca²⁺ sensor STIM1 and Orai plasma membrane Ca²⁺ channels molecularly mediate SOCE. Here we investigate the role of the microRNA machinery in regulating STIM1 expression. We show that STIM1 expression is regulated post-transcriptionally by the miRNA machinery and identify miR-223 and miR-150 as regulators of STIM1 expression in the luminal non-aggressive MCF7 breast cancer cell line. In contrast, STIM1 expression in the more aggressive basal triple-negative MDA-MB-231 cell line is not significantly modulated by a single miRNA species but is rather upregulated due to inhibition of the miRNA machinery through downregulation of Ago2. Consistently, overexpression of Ago2 results in decreased STIM1 protein levels in MDA-MB-231 cells. Clinically, STIM1 and Ago2 expression levels do not correlate with breast cancer progression, however in the basal subtype high STIM1 expression is associated with poorer survival. Our findings show that STIM1 expression is differentially regulated by the miRNA machinery in different cell types and argue for a role for this regulation in breast cancer.

Gene Regulatory Networks in Peripheral Mononuclear Cells Reveals Critical Regulatory Modules and Regulators of Multiple Sclerosis

Multiple sclerosis (MS) is a complex, demyelinating disease with the involvement of autoimmunity and neurodegeneration. Increasing efforts have been made towards identifying the diagnostic markers to differentiate the classes of MS from other similar neurological conditions. Using a systems biology approach, we constructed four types of gene regulatory networks (GRNs) involved in peripheral blood mononuclear cells (PBMCs). The regulatory strength of each GRN across primary progressive MS (PPMS), relapsing-remitting MS (RRMS), secondary progressive MS (SPMS), and control were evaluated by an integrity algorithm. Among the constructed GRNs (referred as TF_gene_miRNA), POU3F2_CDK6_hsa-miR-590-3p, MEIS1_CASC3_hsa-miR-1261, STAT3_OGG1_hsa-miR-298, and TCF4_FMR1_hsa-miR-301b were top-ranked and differentially regulated in all classes of MS compared to control. These GRNs showed potential involvement in regulating various molecular pathways such as interleukin, integrin, glypican, sphingosine phosphate, androgen, and Wnt signaling pathways. For validation, the qPCR analysis of the GRN components (TFs, gene, and miRNAs) in PBMCs of healthy controls (n = 30), RRMS (n = 14), PPMS (n = 13) and SPMS (n = 12) were carried out. Real-time expression analysis of GRNs showed a similar regulatory pattern as derived from our systems biology approach. Also, our study provided several novel GRNs that regulate unique and common molecular mechanisms between MS conditions. Hence, these regulatory components of GRNs will help to understand the disease mechanism across MS classes and further insight may though light towards diagnosis.

Cluster miRNAs and cancer: Diagnostic, prognostic and therapeutic opportunities

MiRNAs are class of noncoding RNA important for gene expression regulation in many plants, animals and viruses. MiRNA clusters contain a set of two or more miRNA encoding genes, transcribed together as polycistronic miRNAs. Currently, there are approximately 159 miRNA clusters reported in the human genome consisting of miRNAs ranging from two or more miRNA genes. A large proportion of clustered miRNAs resides in and around the fragile sites or cancer associated genomic hotspots and plays an important role in carcinogenesis. Altered expression of miRNA cluster can be pro-tumorigenic or anti-tumorigenic and can be targeted for clinical management of cancer. Over the past few years, manipulation of miRNA clusters expression is attempted for experimental purpose as well as for diagnostic, prognostic and therapeutic applications in cancer. Re-expression of miRNAs by epigenetic therapy, genome editing such as clustered regulatory interspaced short palindromic repeats (CRISPR) and miRNA mowers showed promising results in cancer therapy. In this review, we focused on the potential of miRNA clusters as a biomarker for diagnosis, prognosis, targeted therapy as well as strategies for modulating their expression in a therapeutic context. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Processing > Processing of Small RNAs RNA in Disease and Development > RNA in Disease Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Integrated analysis of lncRNA-miRNA-mRNA ceRNA network in squamous cell carcinoma of tongue

Background

Numerous studies have highlighted that long non-coding RNAs (lncRNAs) can bind to microRNA (miRNA) sites as competing endogenous RNAs (ceRNAs), thereby affecting and regulating the expression of mRNAs and target genes. These lncRNA-associated ceRNAs have been theorized to play a significant role in cancer initiation and progression. However, the roles and functions of the lncRNA-miRNA-mRNA ceRNA network in squamous cell carcinoma of the tongue (SCCT) are still unclear.

Methods

The miRNA, mRNA and lncRNA expression profiles from 138 patients with SCCT were downloaded from The Cancer Genome Atlas database. We identified the differential expression of miRNAs, mRNAs, and lncRNAs using the limma package of R software. We used the clusterProfiler package for GO and KEGG pathway annotations. The survival package was used to estimate survival analysis according to the Kaplan-Meier curve. Finally, the GDCRNATools package was used to construct the lncRNA-miRNA-mRNA ceRNA network.

Results

In total, 1943 SCCT-specific mRNAs, 107 lncRNAs and 100 miRNAs were explored. Ten mRNAs (CSRP2, CKS2, ADGRG6, MB21D1, GMNN, RIPOR3, RAD51, PCLAF, ORC1, NAGS), 9 lncRNAs (LINC02560, HOXC13 − AS, FOXD2 − AS1, AC105277.1, AC099850.3, STARD4 − AS1, SLC16A1 − AS1, MIR503HG, MIR100HG) and 8 miRNAs (miR − 654, miR − 503, miR − 450a, miR − 379, miR − 369, miR − 190a, miR − 101, and let−7c) were found to be significantly associated with overall survival (log-rank p < 0.05). Based on the analysis of the lncRNA-miRNA-mRNA ceRNA network, one differentially expressed (DE) lncRNA, five DEmiRNAs, and three DEmRNAs were demonstrated to be related to the pathogenesis of SCCT.

Conclusions

In this study, we described the gene regulation by the lncRNA-miRNA-mRNA ceRNA network in the progression of SCCT. We propose a new lncRNA-associated ceRNA that could help in the diagnosis and treatment of SCCT.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5983-8) contains supplementary material, which is available to authorized users.

mirDIP 4.1-integrative database of human microRNA target predictions

MicroRNAs are important regulators of gene expression, achieved by binding to the gene to be regulated. Even with modern high-throughput technologies, it is laborious and expensive to detect all possible microRNA targets. For this reason, several computational microRNA-target prediction tools have been developed, each with its own strengths and limitations. Integration of different tools has been a successful approach to minimize the shortcomings of individual databases. Here, we present mirDIP v4.1, providing nearly 152 million human microRNA-target predictions, which were collected across 30 different resources. We also introduce an integrative score, which was statistically inferred from the obtained predictions, and was assigned to each unique microRNA-target interaction to provide a unified measure of confidence. We demonstrate that integrating predictions across multiple resources does not cumulate prediction bias toward biological processes or pathways. mirDIP v4.1 is freely available at http://ophid.utoronto.ca/mirDIP/.

miR-181a/b downregulation exerts a protective action on mitochondrial disease models

Mitochondrial diseases (MDs) are a heterogeneous group of devastating and often fatal disorders due to defective oxidative phosphorylation. Despite the recent advances in mitochondrial medicine, effective therapies are still not available for these conditions. Here, we demonstrate that the microRNAs miR-181a and miR-181b (miR-181a/b) regulate key genes involved in mitochondrial biogenesis and function and that downregulation of these miRNAs enhances mitochondrial turnover in the retina through the coordinated activation of mitochondrial biogenesis and mitophagy. We thus tested the effect of miR-181a/b inactivation in different animal models of MDs, such as microphthalmia with linear skin lesions and Leber's hereditary optic neuropathy. We found that miR-181a/b downregulation strongly protects retinal neurons from cell death and significantly ameliorates the disease phenotype in all tested models. Altogether, our results demonstrate that miR-181a/b regulate mitochondrial homeostasis and that these miRNAs may be effective gene-independent therapeutic targets for MDs characterized by neuronal degeneration.

PITX2 enhances progression of lung adenocarcinoma by transcriptionally regulating WNT3A and activating Wnt/β-catenin signaling pathway

Background

The homeodomain transcription factor, PITX2 is associated with tumorigenesis of multiple cancers. In this research, we aimed to study the expression, function and mechanism of PITX2 in lung adenocarcinoma (LUAD).

Methods

The TCGA dataset was used to analyze the expression and clinical significance of PITX2 in LUAD. The expression of PITX2 in tumor samples and LUAD cell lines was examined by quantitative real-time PCR (qRT-PCR) and western blotting. Small interfering RNAs (siRNAs) were constructed to knockdown PITX2 and to determine the physiological function of PITX2 in vitro. Xenograft model was used to confirm the role of PITX2 in vivo.

Results

PITX2 was overexpressed in LUAD and patients with high level of PITX2 had a worse overall survival and an advanced clinical stage. Knockdown of PITX2 inhibited cell proliferation, migration and invasion of LUAD cells. Further study revealed that the oncogenic role of PITX2 was dependent on activating Wnt/β-catenin signaling pathway, especially by transcriptionally regulating the Wnt gene family member, WNT3A. Lastly, we identified miR-140-5p as a negative mediator of PITX2 by binding its 3′UTR and ectopic expression of miR-140-5p inhibited progression of LUAD cells via suppressing the expression of PITX2.

Conclusions

Up-regulation of PITX2 acts as an oncogene in LUAD by activating Wnt/β-catenin signaling pathway, suggesting that PITX2 may serve as a novel diagnostic and prognostic biomarker in LUAD.

Electronic supplementary material

The online version of this article (10.1186/s12935-019-0800-7) contains supplementary material, which is available to authorized users.

Estivation-responsive microRNAs in a hypometabolic terrestrial snail

When faced with extreme environmental conditions, the milk snail (Otala lactea) enters a state of dormancy known as estivation. This is characterized by a strong reduction in metabolic rate to <30% of normal resting rate that is facilitated by various behavioural, physiological, and molecular mechanisms. Herein, we investigated the regulation of microRNA in the induction of estivation. Changes in the expression levels of 75 highly conserved microRNAs were analysed in snail foot muscle, of which 26 were significantly upregulated during estivation compared with controls. These estivation-responsive microRNAs were linked to cell functions that are crucial for long-term survival in a hypometabolic state including anti-apoptosis, cell-cycle arrest, and maintenance of muscle functionality. Several of the microRNA responses by snail foot muscle also characterize hypometabolism in other species and support the existence of a conserved suite of miRNA responses that regulate environmental stress responsive metabolic rate depression across phylogeny.

Circulating miRNAs as a marker of metastatic disease and prognostic factor in metastatic breast cancer

Background

Circulating miRNAs (miRs) are increasingly recognized as potential biomarkers in cancer. We aimed to evaluate the differential expression of miR-23b and miR-190 which are involved in tumor dormancy, miR-21 involved in metastasis and miR-200b and miR-200c involved in epithelial-mesenchymal transition (EMT) and metastasis, in the plasma of patients with early and metastatic breast cancer (MBC). We also aimed to identify associations of the expression levels with patient and disease characteristics and outcomes in metastatic patients treated with first-line chemotherapy.

Results

miR-21 (p < 0.001), miR-23b (p = 0.033), miR-200b (p < 0.001) and miR-200c (p < 0.001) expression was higher in metastatic compared to early breast cancer. ROC curve analysis showed that miR-21 (AUC = 0.722; p < 0.001) and miR-200b (AUC = 0.720; p < 0.001) distinguished with high accuracy among the two disease states, whereas the combination of miR-21, miR-190, miR-200b and miR-200c, further improved accuracy (AUC = 0.797; p < 0.001). High miR-200b expression independently predicted for shorter OS (p = 0.026) in MBC. High expression of both miR23b and miR-190 emerged as a strong independent factor associated with shorter PFS (p = 0.001) in de novo metastatic patients and high miR-200b independently predicted for decreased OS in the HER2-negative subgroup (p = 0.007).

Materials and Methods

Blood samples were obtained from patients with early (n = 133) and MBC (n = 110) before adjuvant or first-line chemotherapy, respectively. Plasma miRNA expression levels were assessed by RT-qPCR and were classified as high or low according to the median values.

Conclusions

Our results are in support of the concept that circulating miRNAs represent a tool with significant diagnostic and prognostic implications in breast cancer.

miRNA-27a-3p and miRNA-222-3p as Novel Modulators of Phosphodiesterase 3a (PDE3A) in Cerebral Microvascular Endothelial Cells

Endothelial dysfunction is a key element in cerebral small vessel disease (CSVD), which may cause stroke and cognitive decline. Cyclic nucleotide signaling modulates endothelial function. The cyclic adenosine monophosphate-degrading enzyme phosphodiesterase 3 (PDE3) is an important treatment target which may be modulated by microRNAs (miRNAs) important for regulating gene expression. We aimed to identify PDE3-targeting miRNAs to highlight potential therapeutic targets for endothelial dysfunction and CSVD. PDE3-targeting miRNAs were identified by in silico analysis (TargetScan, miRWalk, miRanda, and RNA22). The identified miRNAs were ranked on the basis of TargetScan context scores and their expression (log2 read counts) in a human brain endothelial cell line (hCMEC/D3) described recently. miRNAs were subjected to co-expression meta-analysis (CoMeTa) to create miRNA clusters. The pathways targeted by the miRNAs were assigned functional annotations via the KEGG pathway and COOL. hCMEC/D3 cells were transfected with miRNA mimics miR-27a-3p and miR-222-3p, and the effect on PDE3A protein expression was analyzed by Western blotting. Only PDE3A is expressed in hCMEC/D3 cells. The in silico prediction identified 67 PDE3A-related miRNAs, of which 49 were expressed in hCMEC/D3 cells. Further analysis of the top two miRNA clusters (miR-221/miR-222 and miR-27a/miR-27b/miR-128) indicated a potential link to pathways relevant to cerebral and vascular integrity and repair. hCMEC/D3 cells transfected with miR-27a-3p and miR-222-3p mimics had reduced relative expression of PDE3A protein. PDE3A-related miRNAs miR-221/miR-222 and miR-27a/miR-27b/miR-128 are potentially linked to pathways essential for immune regulation as well as cerebral and vascular integrity/function. Furthermore, relative PDE3A protein expression was reduced by miR27a-3p and miR-222-3p.