The DIANA-mirExTra web server: from gene expression data to microRNA function

PVALB Was Identified as an Independent Prognostic Factor for HCC Closely Related to Immunity, and Its Absence Accelerates Tumor Progression by Regulating NK Cell Infiltration

Purpose

Hepatocellular carcinoma is the most common primary liver cancer, with poor prognosis. Complex immune microenvironment of the liver is linked to the development of HCC. PVALB is a calcium-binding protein which has been described as a cancer suppressor gene in thyroid cancer and glioma. Nevertheless, the role of PVALB in HCC is unknown.

Materials and Methods

We obtained data from TCGA and GSE54236 datasets. MCP-counter, WGCNA and LASSO model were applied to identify PVALB. With UALCAN, MethSurv, and other websites, we probed the expression, methylation and survival of PVALB. LinkedOmics and GSEA were adopted for functional analysis, while TIMER, TISIDB, Kaplan-Meier plotter, TIDE databases were utilized to evaluate the relevance of PVALB to the tumor immune microenvironment and predict immunotherapy efficacy. TargetScan, DIANA, LncRNASNP2 databases and relevant experiments were employed to construct ceRNA network. Finally, molecular docking and drug sensitivity of PVALB were characterized by GeneMANIA, CTD, and so on.

Results

PVALB was recognized as a gene associated with HCC and NK cell. Its expression was down-regulated in HCC tissue, which lead to adverse prognosis. Besides, the hypomethylation of PVALB was related to its reduced expression. Notably, PVALB was tightly linked to immune, and its reduced expression attenuated the anticancer effect of NK cells via the Fas/FasL pathway, leading to a adverse outcome. The lnc-YY1AP1-3/hsa-miR-6735-5p/PVALB axis may regulate the PVALB expression. Finally, we found immunotherapy might be a viable treatment option.

Conclusion

In a word, PVALB is a prognostic indicator, whose low expression facilitates HCC progression by impacting NK cell infiltration.

Fetal sex modulates placental microRNA expression, potential microRNA-mRNA interactions, and levels of amino acid transporter expression and substrates: INFAT study subpopulation analysis of n-3 LCPUFA intervention during pregnancy and associations with offspring body composition

Background

Previously, we revealed sexually dimorphic mRNA expression and responsiveness to maternal dietary supplementation with n-3 long-chain polyunsaturated fatty acids (LCPUFA) in placentas from a defined INFAT study subpopulation. Here, we extended these analyses and explored the respective placental microRNA expression, putative microRNA-mRNA interactions, and downstream target processes as well as their associations with INFAT offspring body composition.

Results

We performed explorative placental microRNA profiling, predicted microRNA-mRNA interactions by bioinformatics, validated placental target microRNAs and their putative targets by RT-qPCR and western blotting, and measured amino acid levels in maternal and offspring cord blood plasma and placenta. microRNA, mRNA, protein, and amino acid levels were associated with each other and with offspring body composition from birth to 5 years of age. Forty-six differentially regulated microRNAs were found. Validations identified differential expression for microRNA-99a (miR-99a) and its predicted target genes mTOR, SLC7A5, encoding L-type amino acid transporter 1 (LAT1), and SLC6A6, encoding taurine transporter (TauT), and their prevailing significant sexually dimorphic regulation. Target mRNA levels were mostly higher in placentas from control male than from female offspring, whereas respective n-3 LCPUFA responsive target upregulation was predominantly found in female placentas, explaining the rather balanced expression levels between the sexes present only in the intervention group. LAT1 and TauT substrates tryptophan and taurine, respectively, were significantly altered in both maternal plasma at 32 weeks’ gestation and cord plasma following intervention, but not in the placenta. Several significant associations were observed for miR-99a, mTOR mRNA, SLC7A5 mRNA, and taurine and tryptophan in maternal and cord plasma with offspring body composition at birth, 1 year, 3 and 5 years of age.

Conclusions

Our data suggest that the analyzed targets may be part of a sexually dimorphic molecular regulatory network in the placenta, possibly modulating gene expression per se and/or counteracting n-3 LCPUFA responsive changes, and thereby stabilizing respective placental and fetal amino acid levels. Our data propose placental miR-99, SLC7A5 mRNA, and taurine and tryptophan levels in maternal and fetal plasma as potentially predictive biomarkers for offspring body composition.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00345-x.

Regulation of Pancreatic Fibrosis by Acinar Cell-Derived Exosomal miR-130a-3p via Targeting of Stellate Cell PPAR-γ

Introduction

As endogenous miRNA carriers, exosomes play a role in the pathophysiological processes of various diseases. However, their functions and regulation mechanisms in pancreatic fibrosis remain unclear.

Methods

In this study, an RNA microarray was used to detect differentially expressed exosomal miR-130a-3p in AR42J cells before and after taurolithocholate (TLC) treatment. mRNA-seq was used to screen differentially expressed genes before and after pancreatic stellate cell (PSC) activation. We used the STRING database to construct a protein-protein interaction (PPI) network for differentially expressed genes, used CytoNCA to analyze the centrality of the PPI network, and identified 10 essential proteins in the biological network. Then, the TargetScan and miRanda databases were used to predict the target genes of miR-130a-3p. The intersections of the target genes and the mRNAs encoding the 10 essential proteins were identified to construct miR-130a-3p/peroxisome proliferator-activated receptor gamma (PPAR-γ) pairs. Fluorescence labeling of exosomes and dynamic tracing showed that exosomes can fuse with the cell membranes of PSCs and transport miR-130a-3p into PSCs. A luciferase reporter gene assay was used to confirm that miR-130a-3p can bind to PPAR-γ to inhibit PPAR-γ expression. In vitro and in vivo functional experiments were performed for gain-of-function studies and loss-of-function studies, respectively.

Results

The studies showed that acinar cell-derived exosomal miR-130a-3p promotes PSC activation and collagen formation through targeting of stellate cellular PPAR-γ. Knockdown of miR-130a-3p significantly improved pancreatic fibrosis. Notably, miR-130a-3p knockdown reduced serum levels of hyaluronic acid (HA) and β-amylase and increased the C-peptide level to protect endocrine and exocrine pancreatic functions and the function of endothelial cells.

Conclusion

This study revealed that the exosomal miR-130a-3p/PPAR-γ axis participates in PSC activation and the mechanism of chronic pancreatitis (CP) with fibrosis, thus providing a potential new target for the treatment of chronic pancreatic fibrosis.

Transcriptional Regulation of Cancer Immune Checkpoints: Emerging Strategies for Immunotherapy

The study of immune evasion has gained a well-deserved eminence in cancer research by successfully developing a new class of therapeutics, immune checkpoint inhibitors, such as pembrolizumab and nivolumab, anti-PD-1 antibodies. By aiming at the immune checkpoint blockade (ICB), these new therapeutics have advanced cancer treatment with notable increases in overall survival and tumor remission. However, recent reports reveal that 40-60% of patients fail to benefit from ICB therapy due to acquired resistance or tumor relapse. This resistance may stem from increased expression of co-inhibitory immune checkpoints or alterations in the tumor microenvironment that promotes immune suppression. Because these mechanisms are poorly elucidated, the transcription factors that regulate immune checkpoints, known as "master regulators", have garnered interest. These include AP-1, IRF-1, MYC, and STAT3, which are known to regulate PD/PD-L1 and CTLA-4. Identifying these and other potential master regulators as putative therapeutic targets or biomarkers can be facilitated by mining cancer literature, public datasets, and cancer genomics resources. In this review, we describe recent advances in master regulator identification and characterization of the mechanisms underlying immune checkpoints regulation, and discuss how these master regulators of immune checkpoint molecular expression can be targeted as a form of auxiliary therapeutic strategy to complement traditional immunotherapy.

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.

Human microRNA-30 inhibits influenza virus infection by suppressing the expression of SOCS1, SOCS3, and NEDD4

Influenza A virus (IAV) has evolved multiple mechanisms to compromise type I interferon (IFN) responses. The antiviral function of IFN is mainly exerted by activating the JAK/STAT signalling and subsequently inducing IFN‐stimulated gene (ISG) production. However, the mechanism by which IAV combat the type I IFN signalling pathway is not fully elucidated. In this study, we explored the roles of human microRNAs modulated by IAV infection in type I IFN responses. We demonstrated that microRNA‐30 (miR‐30) family members were downregulated by IAV infection. Our data showed that the forced expression of miR‐30 family members inhibited IAV proliferation, while miR‐30 family member inhibitors promoted IAV proliferation. Mechanistically, we found that miR‐30 family members targeted and reduced SOCS1 and SOCS3 expression, and thus relieved their inhibiting effects on IFN/JAK/STAT signalling pathway. In addition, miR‐30 family members inhibited the expression of NEDD4, a negative regulator of IFITM3, which is important for host defence against influenza viruses. Our findings suggest that IAV utilises a novel strategy to restrain host type I IFN‐mediated antiviral immune responses by decreasing the expression of miR‐30 family members, and add a new way to understand the mechanism of immune escape caused by influenza viruses.

miRmapper: A Tool for Interpretation of miRNA⁻mRNA Interaction Networks

It is estimated that 30% of all genes in the mammalian cells are regulated by microRNA (miRNAs). The most relevant miRNAs in a cellular context are not necessarily those with the greatest change in expression levels between healthy and diseased tissue. Differentially expressed (DE) miRNAs that modulate a large number of messenger RNA (mRNA) transcripts ultimately have a greater influence in determining phenotypic outcomes and are more important in a global biological context than miRNAs that modulate just a few mRNA transcripts. Here, we describe the development of a tool, “miRmapper”, which identifies the most dominant miRNAs in a miRNA–mRNA network and recognizes similarities between miRNAs based on commonly regulated mRNAs. Using a list of miRNA–target gene interactions and a list of DE transcripts, miRmapper provides several outputs: (1) an adjacency matrix that is used to calculate miRNA similarity utilizing the Jaccard distance; (2) a dendrogram and (3) an identity heatmap displaying miRNA clusters based on their effect on mRNA expression; (4) a miRNA impact table and (5) a barplot that provides a visual illustration of this impact. We tested this tool using nonmetastatic and metastatic bladder cancer cell lines and demonstrated that the most relevant miRNAs in a cellular context are not necessarily those with the greatest fold change. Additionally, by exploiting the Jaccard distance, we unraveled novel cooperative interactions between miRNAs from independent families in regulating common target mRNAs; i.e., five of the top 10 miRNAs act in synergy.

From Oxidative Stress Damage to Pathways, Networks, and Autophagy via MicroRNAs

Oxidative stress can alter the expression level of many microRNAs (miRNAs), but how these changes are integrated and related to oxidative stress responses is poorly understood. In this article, we addressed this question by using in silico tools. We reviewed the literature for miRNAs whose expression is altered upon oxidative stress damage and used them in combination with various databases and software to predict common gene targets of oxidative stress-modulated miRNAs and affected pathways. Furthermore, we identified miRNAs that simultaneously target the predicted oxidative stress-modulated miRNA gene targets. This generated a list of novel candidate miRNAs potentially involved in oxidative stress responses. By literature search and grouping of pathways and cellular responses, we could classify these candidate miRNAs and their targets into a larger scheme related to oxidative stress responses. To further exemplify the potential of our approach in free radical research, we used our explorative tools in combination with ingenuity pathway analysis to successfully identify new candidate miRNAs involved in the ubiquitination process, a master regulator of cellular responses to oxidative stress and proteostasis. Lastly, we demonstrate that our approach may also be useful to identify novel candidate connections between oxidative stress-related miRNAs and autophagy. In summary, our results indicate novel and important aspects with regard to the integrated biological roles of oxidative stress-modulated miRNAs and demonstrate how this type of in silico approach can be useful as a starting point to generate hypotheses and guide further research on the interrelation between miRNA-based gene regulation, oxidative stress signaling pathways, and autophagy.

Identification of a serum-induced transcriptional signature associated with metastatic cervical cancer

Objective

Tumor cells that escape local tissue control can convert inflammatory cells from tumor suppressors to tumor promoters. Moreover, soluble immune-modulating factors secreted from the tumor environment can be difficult to identify in patient serum due to their low abundance. We used an alternative strategy to infer a metastatic signature induced by sera of cervical cancer patients.

Methods

Sera from patients with local and metastatic cervical cancer were used to induce a disease-specific transcriptional signature in cultured, healthy peripheral blood mononuclear cells (PBMCs). An empirical Bayesian method, EBarrays, was used to identify differentially expressed (DE) genes with a target false discovery rate of <5%. Ingenuity Pathway Analysis (IPA) software was used to detect the top molecular and cellular functions associated with the DE genes. IPA and in silco analysis was used to pinpoint candidate upstream regulators, including cancer-related microRNAs (miRNAs).

Results

We identified enriched pathways in the metastatic cervical group related to immune surveillance functions, such as downregulation of engulfment, accumulation, and phagocytosis of hematopoietic cells. The predicted top upstream genes were IL-10 and immunoglobulins. In silco analysis identified miRNAs predicted to drive the transcriptional signature. Two of the 4 miRNAs (miR-23a-3p and miR-944) were validated in a cohort of women with local and metastatic cervical cancer.

Conclusions

This study supports the use of a cell-based assay that uses PBMC “reporters” to predict biologically relevant factors in patient serum. Further, disease-specific transcriptional signatures induced by patient sera have the potential to differentiate patients with local versus metastatic disease.

A Functional Genomics Approach to Henipavirus Research: The Role of Nuclear Proteins, MicroRNAs and Immune Regulators in Infection and Disease

Hendra and Nipah viruses (family Paramyxoviridae, genus Henipavirus) are zoonotic RNA viruses that cause lethal disease in humans and are designated as Biosafety Level 4 (BSL4) agents. Moreover, henipaviruses belong to the same group of viruses that cause disease more commonly in humans such as measles, mumps and respiratory syncytial virus. Due to the relatively recent emergence of the henipaviruses and the practical constraints of performing functional genomics studies at high levels of containment, our understanding of the henipavirus infection cycle is incomplete. In this chapter we describe recent loss-of-function (i.e. RNAi) functional genomics screens that shed light on the henipavirus-host interface at a genome-wide level. Further to this, we cross-reference RNAi results with studies probing host proteins targeted by henipavirus proteins, such as nuclear proteins and immune modulators. These functional genomics studies join a growing body of evidence demonstrating that nuclear and nucleolar host proteins play a crucial role in henipavirus infection. Furthermore these studies will underpin future efforts to define the role of nucleolar host-virus interactions in infection and disease.

miR-331-3p and Aurora Kinase inhibitor II co-treatment suppresses prostate cancer tumorigenesis and progression

RNA-based therapeutics could represent a new avenue of cancer treatment. miRNA 331-3p (miR-331-3p) is implicated in prostate cancer (PCa) as a putative tumor suppressor, but its functional activity and synergy with other anti-tumor agents is largely unknown. We found miR-331-3p expression in PCa tumors was significantly decreased compared to non-malignant matched tissue. Analysis of publicly available PCa gene expression data sets showed miR-331-3p expression negatively correlated with Gleason Score, tumor stage, lymph node involvement and PSA value, and was significantly down regulated in tumor tissue relative to normal prostate tissue. Overexpression of miR-331-3p reduced PCa cell growth, migration and colony formation, as well as xenograft tumor initiation, proliferation and survival of mice. Microarray analysis identified seven novel targets of miR-331-3p in PCa. The 3’-untranslated regions of PLCγ1 and RALA were confirmed as targets of miR-331-3p, with mutation analyses confirming RALA as a direct target. Expression of miR-331-3p or RALA siRNA in PCa cells reduced RALA expression, proliferation, migration and colony formation in vitro. RALA expression positively correlated with Gleason grade in two separate studies, as well as in a PCa tissue microarray. Co-treatment using siRALA with an Aurora Kinase inhibitor (AKi-II) decreased colony formation of PCa cells while the combination of AKi-II with miR-331-3p resulted in significant reduction of PCa cell proliferation in vitro and PCa xenograft growth in vivo. Thus, miR-331-3p directly targets the RALA pathway and the addition of the AKi-II has a synergistic effect on tumor growth inhibition, suggesting a potential role as combination therapy in PCa.

Functional Analysis of miRNAs Using the DIANA Tools Online Suite

microRNAs (miRNAs) are central regulators of gene expression. They are actively studied for their involvement in numerous physiological and pathological conditions but also as diagnostic biomarkers or promising therapeutic targets. The increased complexity of the miRNA interactomes hinders straightforward interpretation of miRNA expression differences between states and conditions. To this end, functional analysis web servers process and combine experimental and in silico data, enabling researchers to uncover targeted pathways and transcriptional mechanisms that are hidden within numerous interactions and vast expression datasets. DIANA-tools ( www.microrna.gr ) is a web server hosting state-of-the-art utilities and databases for miRNA functional investigation. Available utilities cover a wide scope of different needs and research scenarios, rendering DIANA website a one-stop-shop for miRNA analyses. The most commonly utilized databases and algorithms include DIANA-microT-CDS, DIANA-TarBase v7.0, DIANA-lncBase v2.0, DIANA-miRGen v3.0, DIANA-miRPath v3.0, and DIANA-mirExTra v2.0.In the presented protocol, we will utilize different online tools in order to explore miRNA functions and to identify probable targets of interest for downstream analyses and wet lab experiments. The combined use of different applications from the DIANA suite can shed light to numerous different aspects of miRNA regulation and regulatory function, without the necessity for extensive bioinformatics expertise or computational infrastructure.

Practical Bioinformatics Analysis of MiRNA Data Using Online Tools

The purpose of this chapter is to provide a starting point for the analysis of miRNA array data, using freely available online suites of tools. This chapter does not describe how to perform analysis of primary array data, rather how to use the top differentially regulated miRNA (returned from comparing one miRNA group to another) as the starting point for further practical analysis.Here we describe the methods and tools required to identify targets worthy of additional investigation, using the identified miRNA as a starting point. Importantly, this additional information (pathways targeted, gene expression, mRNA targets, miRNA families) can be used to positively inform any project.

DIANA-mirExTra v2.0: Uncovering microRNAs and transcription factors with crucial roles in NGS expression data

Differential expression analysis (DEA) is one of the main instruments utilized for revealing molecular mechanisms in pathological and physiological conditions. DIANA-mirExTra v2.0 (http://www.microrna.gr/mirextrav2) performs a combined DEA of mRNAs and microRNAs (miRNAs) to uncover miRNAs and transcription factors (TFs) playing important regulatory roles between two investigated states. The web server uses as input miRNA/RNA-Seq read count data sets that can be uploaded for analysis. Users can combine their data with 350 small-RNA-Seq and 65 RNA-Seq in-house analyzed libraries which are provided by DIANA-mirExTra v2.0.

The web server utilizes miRNA:mRNA, TF:mRNA and TF:miRNA interactions derived from extensive experimental data sets. More than 450 000 miRNA interactions and 2 000 000 TF binding sites from specific or high-throughput techniques have been incorporated, while accurate miRNA TSS annotation is obtained from microTSS experimental/in silico framework. These comprehensive data sets enable users to perform analyses based solely on experimentally supported information and to uncover central regulators within sequencing data: miRNAs controlling mRNAs and TFs regulating mRNA or miRNA expression. The server also supports predicted miRNA:gene interactions from DIANA-microT-CDS for 4 species (human, mouse, nematode and fruit fly). DIANA-mirExTra v2.0 has an intuitive user interface and is freely available to all users without any login requirement.

Bioinformatic tools for microRNA dissection

Recently, microRNAs (miRNAs) have emerged as important elements of gene regulatory networks. MiRNAs are endogenous single-stranded non-coding RNAs (~22-nt long) that regulate gene expression at the post-transcriptional level. Through pairing with mRNA, miRNAs can down-regulate gene expression by inhibiting translation or stimulating mRNA degradation. In some cases they can also up-regulate the expression of a target gene. MiRNAs influence a variety of cellular pathways that range from development to carcinogenesis. The involvement of miRNAs in several human diseases, particularly cancer, makes them potential diagnostic and prognostic biomarkers. Recent technological advances, especially high-throughput sequencing, have led to an exponential growth in the generation of miRNA-related data. A number of bioinformatic tools and databases have been devised to manage this growing body of data. We analyze 129 miRNA tools that are being used in diverse areas of miRNA research, to assist investigators in choosing the most appropriate tools for their needs.

Identification of Molecular Markers of Delayed Graft Function Based on the Regulation of Biological Ageing

Introduction

Delayed graft function is a prevalent clinical problem in renal transplantation for which there is no objective system to predict occurrence in advance. It can result in a significant increase in the necessity for hospitalisation post-transplant and is a significant risk factor for other post-transplant complications.

Methodology

The importance of microRNAs (miRNAs), a specific subclass of small RNA, have been clearly demonstrated to influence many pathways in health and disease. To investigate the influence of miRNAs on renal allograft performance post-transplant, the expression of a panel of miRNAs in pre-transplant renal biopsies was measured using qPCR. Expression was then related to clinical parameters and outcomes in two independent renal transplant cohorts.

Results

Here we demonstrate, in two independent cohorts of pre-implantation human renal allograft biopsies, that a novel pre-transplant renal performance scoring system (GRPSS), can determine the occurrence of DGF with a high sensitivity (>90%) and specificity (>60%) for donor allografts pre-transplant, using just three senescence associated microRNAs combined with donor age and type of organ donation.

Conclusion

These results demonstrate a relationship between pre-transplant microRNA expression levels, cellular biological ageing pathways and clinical outcomes for renal transplantation. They provide for a simple, rapid quantitative molecular pre-transplant assay to determine post-transplant allograft function and scope for future intervention. Furthermore, these results demonstrate the involvement of senescence pathways in ischaemic injury during the organ transplantation process and an indication of accelerated bio-ageing as a consequence of both warm and cold ischaemia.

Subclinical Detection of Diabetic Cardiomyopathy with MicroRNAs: Challenges and Perspectives

The prevalence of cardiac diabetic diseases has been increased around the world, being the most common cause of death and disability among diabetic patients. In particular, diabetic cardiomyopathy is characterized with a diastolic dysfunction and cardiac remodelling without signs of hypertension and coronary artery diseases. In an early stage, it is an asymptomatic disease; however, clinical studies demonstrate that diabetic myocardia are more vulnerable to injury derived by acute myocardial infarct and are the worst prognosis for rehabilitation. Currently, biochemical and imaging diagnostic methods are unable to detect subclinical manifestation of the disease (prior to diastolic dysfunction). In this review, we elaborately discuss the current scientific evidences to propose circulating microRNAs as promising biomarkers for early detection of diabetic cardiomyopathy and, then, to identify patients at high risk of diabetic cardiomyopathy development. Moreover, here we summarise the research strategies to identify miRNAs as potential biomarkers, present limitations, challenges, and future perspectives.

Human MicroRNA miR-532-5p Exhibits Antiviral Activity against West Nile Virus via Suppression of Host Genes SESTD1 and TAB3 Required for Virus Replication

West Nile virus (WNV) is a mosquito-transmitted flavivirus that naturally circulates between mosquitos and birds but can also infect humans, causing severe neurological disease. The early host response to WNV infection in vertebrates primarily relies on the type I interferon pathway; however, recent studies suggest that microRNAs (miRNAs) may also play a notable role. In this study, we assessed the role of host miRNAs in response to WNV infection in human cells. We employed small RNA sequencing (RNA-seq) analysis to determine changes in the expression of host miRNAs in HEK293 cells infected with an Australian strain of WNV, Kunjin (WNVKUN), and identified a number of host miRNAs differentially expressed in response to infection. Three of these miRNAs were confirmed to be significantly upregulated in infected cells by quantitative reverse transcription (qRT)-PCR and Northern blot analyses, and one of them, miR-532-5p, exhibited a significant antiviral effect against WNVKUN infection. We have demonstrated that miR-532-5p targets and downregulates expression of the host genes SESTD1 and TAB3 in human cells. Small interfering RNA (siRNA) depletion studies showed that both SESTD1 and TAB3 were required for efficient WNVKUN replication. We also demonstrated upregulation of mir-532-5p expression and a corresponding decrease in the expression of its targets, SESTD1 and TAB3, in the brains of WNVKUN -infected mice. Our results show that upregulation of miR-532-5p and subsequent suppression of the SESTD1 and TAB3 genes represent a host antiviral response aimed at limiting WNVKUN infection and highlight the important role of miRNAs in controlling RNA virus infections in mammalian hosts.

IMPORTANCE

West Nile virus (WNV) is a significant viral pathogen that poses a considerable threat to human health across the globe. There is no specific treatment or licensed vaccine available for WNV, and deeper insight into how the virus interacts with the host is required to facilitate their development. In this study, we addressed the role of host microRNAs (miRNAs) in antiviral response to WNV in human cells. We identified miR-532-5p as a novel antiviral miRNA and showed that it is upregulated in response to WNV infection and suppresses the expression of the host genes TAB3 and SESTD1 required for WNV replication. Our results show that upregulation of miR-532-5p and subsequent suppression of the SESTD1 and TAB3 genes represent an antiviral response aimed at limiting WNV infection and highlight the important role of miRNAs in controlling virus infections in mammalian hosts.

Regulation signature of miR-143 and miR-26 in porcine Salmonella infection identified by binding site enrichment analysis

Salmonella infects many vertebrate species, and pigs colonized with Salmonella are typically Salmonella carriers. Transcriptomic analysis of the response to Salmonella infection in whole blood has been reported for the pig. The objective of this study is to identify the important miRNAs involved in Salmonella infection using binding site enrichment analysis. We predicted porcine microRNA (miRNA) binding sites in the 3' UTR of protein-coding genes for all miRNA families. Based on those predictions, we analyzed miRNA-binding sites for mRNAs expressed in peripheral blood to investigate the functional importance of miRNAs in Salmonella infection in pig. Enrichment analysis revealed that binding sites of five miRNAs (including miR-143, -9839, -26, -2483, and -4335) were significantly over represented for the differentially expressed gene sets. Real-time PCR results indicated that selected members of this miRNA group (miR-143, -26, and -4335) were differentially expressed in whole blood after Salmonella inoculation. The luciferase reporter assay showed that ATP6V1A and IL13RA1 were targets of miR-143 and that miR-26 regulates BINP3L and ARL6IP6. The results strongly suggest that miR-143 and miR-26 play important regulatory roles in the development of Salmonella infection in pig.

MicroRNA-363 targets myosin 1B to reduce cellular migration in head and neck cancer

Background

Squamous cell carcinoma of the head and neck (SCCHN) remains a prevalent and devastating disease. Recently, there has been an increase in SCCHN cases that are associated with high-risk human papillomavirus (HPV) infection. The clinical characteristics of HPV-positive and HPV-negative SCCHN are known to be different but their molecular features are only recently beginning to emerge. MicroRNAs (miRNAs, miRs) are small, non-coding RNAs that are likely to play significant roles in cancer initiation and progression where they may act as oncogenes or tumor suppressors. Previous studies in our laboratory showed that miR-363 is overexpressed in HPV-positive compared to HPV-negative SCCHN cell lines, and the HPV type 16-E6 oncoprotein upregulates miR-363 in SCCHN cell lines. However, the functional role of miR-363 in SCCHN in the context of HPV infection remains to be elucidated.

Methods

We analyzed miR-363 levels in SCCHN tumors with known HPV-status from The Cancer Genome Atlas (TCGA) and an independent cohort from our institution. Cell migration studies were conducted following the overexpression of miR-363 in HPV-negative cell lines. Bioinformatic tools and a luciferase reporter assay were utilized to confirm that miR-363 targets the 3’-UTR of myosin 1B (MYO1B). MYO1B mRNA and protein expression levels were evaluated following miR-363 overexpression in HPV-negative SCCHN cell lines. Small interfering RNA (siRNA) knockdown of MYO1B was performed to assess the phenotypic implication of reduced MYO1B expression in SCCHN cell lines.

Results

MiR-363 was found to be overexpressed in HPV-16-positive compared to the HPV-negative SCCHN tumors. Luciferase reporter assays performed in HPV-negative JHU028 cells confirmed that miR-363 targets one of its two potential binding sites in the 3’UTR of MYO1B. MYO1B mRNA and protein levels were reduced upon miR-363 overexpression in four HPV-negative SCCHN cell lines. Increased miR-363 expression or siRNA knockdown of MYO1B expression reduced Transwell migration of SCCHN cell lines, indicating that the miR-363-induced migration attenuation of SCCHN cells may act through MYO1B downregulation.

Conclusions

These findings demonstrate that the overexpression of miR-363 reduces cellular migration in head and neck cancer and reveal the biological relationship between miR-363, myosin 1b, and HPV-positive SCCHN.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1888-3) contains supplementary material, which is available to authorized users.

MicroRNA-132 enhances transition from inflammation to proliferation during wound healing

Wound healing is a complex process that is characterized by an initial inflammatory phase followed by a proliferative phase. This transition is a critical regulatory point; however, the factors that mediate this process are not fully understood. Here, we evaluated microRNAs (miRs) in skin wound healing and characterized the dynamic change of the miRNome in human skin wounds. miR-132 was highly upregulated during the inflammatory phase of wound repair, predominantly expressed in epidermal keratinocytes, and peaked in the subsequent proliferative phase. TGF-β1 and TGF-β2 induced miR-132 expression in keratinocytes, and transcriptome analysis of these cells revealed that miR-132 regulates a large number of immune response- and cell cycle-related genes. In keratinocytes, miR-132 decreased the production of chemokines and the capability to attract leukocytes by suppressing the NF-κB pathway. Conversely, miR-132 increased activity of the STAT3 and ERK pathways, thereby promoting keratinocyte growth. Silencing of the miR-132 target heparin-binding EGF-like growth factor (HB-EGF) phenocopied miR-132 overexpression in keratinocytes. Using mouse and human ex vivo wound models, we found that miR-132 blockade delayed healing, which was accompanied by severe inflammation and deficient keratinocyte proliferation. Together, our results indicate that miR-132 is a critical regulator of skin wound healing that facilitates the transition from the inflammatory to the proliferative phase.

Tissue elasticity regulated tumor gene expression: implication for diagnostic biomarkers of primitive neuroectodermal tumor

BACKGROUND

The tumor microenvironment consists of both physical and chemical factors. Tissue elasticity is one physical factor contributing to the microenvironment of tumor cells. To test the importance of tissue elasticity in cell culture, primitive neuroectodermal tumor (PNET) stem cells were cultured on soft polyacrylamide (PAA) hydrogel plates that mimics the elasticity of brain tissue compared with PNET on standard polystyrene (PS) plates. We report the molecular profiles of PNET grown on either PAA or PS.

METHODOLOGY/PRINCIPAL FINDINGS

A whole-genome microarray profile of transcriptional expression between the two culture conditions was performed as a way to probe effects of substrate on cell behavior in culture. The results showed more genes downregulated on PAA compared to PS. This led us to propose microRNA (miRNA) silencing as a potential mechanism for downregulation. Bioinformatic analysis predicted a greater number of miRNA binding sites from the 3' UTR of downregulated genes and identified as specific miRNA binding sites that were enriched when cells were grown on PAA-this supports the hypothesis that tissue elasticity plays a role in influencing miRNA expression. Thus, Dicer was examined to determine if miRNA processing was affected by tissue elasticity. Dicer genes were downregulated on PAA and had multiple predicted miRNA binding sites in its 3' UTR that matched the miRNA binding sites found enriched on PAA. Many differentially regulated genes were found to be present on PS but downregulated on PAA were mapped onto intron sequences. This suggests expression of alternative polyadenylation sites within intron regions that provide alternative 3' UTRs and alternative miRNA binding sites. This results in tissue specific transcriptional downregulation of mRNA in humans by miRNA. We propose a mechanism, driven by the physical characteristics of the microenvironment by which downregulation of genes occur. We found that tissue elasticity-mediated cytokines (TGFβ2 and TNFα) signaling affect expression of ECM proteins.

CONCLUSIONS

Our results suggest that tissue elasticity plays important roles in miRNA expression, which, in turn, regulate tumor growth or tumorigenicity.

APP intracellular domain acts as a transcriptional regulator of miR-663 suppressing neuronal differentiation

Amyloid precursor protein (APP) is best known for its involvement in the pathogenesis of Alzheimer's disease. We have previously demonstrated that APP intracellular domain (AICD) regulates neurogenesis; however, the mechanisms underlying AICD-mediated regulation of neuronal differentiation are not yet fully characterized. Using genome-wide chromatin immunoprecipitation approaches, we found that AICD is specifically recruited to the regulatory regions of several microRNA genes, and acts as a transcriptional regulator for miR-663, miR-3648 and miR-3687 in human neural stem cells. Functional assays show that AICD negatively modulates neuronal differentiation through miR-663, a primate-specific microRNA. Microarray data further demonstrate that miR-663 suppresses the expression of multiple genes implicated in neurogenesis, including FBXL18 and CDK6. Our results indicate that AICD has a novel role in suppression of neuronal differentiation via transcriptional regulation of miR-663 in human neural stem cells.

Sequence-dependent off-target inhibition of TLR7/8 sensing by synthetic microRNA inhibitors

Anti-microRNA (miRNA) oligonucleotides (AMOs) with 2'-O-Methyl (2'OMe) residues are commonly used to study miRNA function and can achieve high potency, with low cytotoxicity. Not withstanding this, we demonstrate the sequence-dependent capacity of 2'OMe AMOs to inhibit Toll-like receptor (TLR) 7 and 8 sensing of immunostimulatory RNA, independent of their miRNA-targeting function. Through a screen of 29 AMOs targeting common miRNAs, we found a subset of sequences highly inhibitory to TLR7 sensing in mouse macrophages. Interspecies conservation of this inhibitory activity was confirmed on TLR7/8 activity in human peripheral blood mononuclear cells. Significantly, we identified a core motif governing the inhibitory activity of these AMOs, which is present in more than 50 AMOs targeted to human miRNAs in miRBaseV20. DNA/locked nucleic acids (LNA) AMOs synthesized with a phosphorothioate backbone also inhibited TLR7 sensing in a sequence-dependent manner, demonstrating that the off-target effects of AMOs are not restricted to 2'OMe modification. Taken together, our work establishes the potential for off-target effects of AMOs on TLR7/8 function, which should be taken into account in their therapeutic development and in vivo application.

miR-24 limits aortic vascular inflammation and murine abdominal aneurysm development

Identification and treatment of abdominal aortic aneurysm (AAA) remain among the most prominent challenges in vascular medicine. MicroRNAs (miRNAs) are crucial regulators of cardiovascular pathology and represent intriguing targets to limit AAA expansion. Here we show, by using two established murine models of AAA disease along with human aortic tissue and plasma analysis, that miR-24 is a key regulator of vascular inflammation and AAA pathology. In vivo and in vitro studies reveal chitinase 3-like 1 (Chi3l1) to be a major target and effector under the control of miR-24, regulating cytokine synthesis in macrophages as well as their survival, promoting aortic smooth muscle cell migration and cytokine production, and stimulating adhesion molecule expression in vascular endothelial cells. We further show that modulation of miR-24 alters AAA progression in animal models, and that miR-24 and CHI3L1 represent novel plasma biomarkers of AAA disease progression in humans.

Abdominal aortic aneurysm (AAA) is a potentially fatal and often asymptomatic disease whose causes remain unclear. Here the authors show that a microRNA, miR-24, and its target, the glycoprotein chitinase 3-like 1, represent key regulators of AAA development.

miRror2.0: a platform for assessing the joint action of microRNAs in cell regulation

microRNAs (miRNAs) are short, noncoding RNAs that negatively regulate the levels of mRNA post-transcriptionally. Recent experiments revealed thousands of mRNA-miRNA pairs in which multiple miRNAs may bind the same transcript. These results raised the notion of miRNAs teamwork for a wide range of cellular context. miRror2.0 utilizes the miRNA-target predictions from over a dozen programs and resources and unifies them under a common statistical basis. The platform, called miRror2.0, considers the combinatorial regulation by miRNAs in different tissues, cell lines and under a broad range of conditions. A flexible setting permits the selection of the preferred combination of miRNA-target prediction resources as well as the statistical parameters for the analysis. miRror2.0 covers six major model organisms including human and mouse. Importantly, the system is capable of analyzing hundreds of genes that were subjected to miRNAs' regulation. Activating miRror2.0 by introducing thousands of genes from miRNA overexpression experiments successfully identified the objective miRNAs. The output from miRror2.0 is a list of genes that is optimally regulated by a defined set of miRNAs. A symmetric application of miRror2.0 starts with a set of miRNAs, and the system then seeks the preferred set of genes that are regulated by that miRNA composition. The results from miRror2.0 are empowered by an iterative procedure called PSI-miRror. PSI-miRror tests the robustness of miRror2.0 prediction. It allows a refinement of the initial list of genes in view of the miRNAs that optimally regulate this list. We present miRror2.0 as a valuable resource for supporting cellular experimentalists that seek recovery of combinatorial regulation by miRNAs from noisy experimental data. miRror2.0 is available at http://www.mirrorsuite.cs.huji.ac.il .

MicroRNA target identification: lessons from hypoxamiRs

SIGNIFICANCE

MicroRNAs (miRNAs) are small noncoding RNAs that have emerged as key regulators of many physiological and pathological processes, including those relevant to hypoxia such as cancer, neurological dysfunctions, myocardial infarction, and lung diseases.

RECENT ADVANCES

During the last 5 years, miRNAs have been shown to play a role in the regulation of the cellular response to hypoxia. The identification of several bona fide targets of these hypoxamiRs has underlined their pleiotropic functions and the complexity of the molecular rules directing miRNA::target transcript pairing.

CRITICAL ISSUES

This review outlines the main in silico and experimental approaches used to identify the targetome of hypoxamiRs and presents new recent relevant methodologies for future studies.

FUTURE DIRECTIONS

Since hypoxia plays key roles in many pathophysiological conditions, the precise characterization of regulatory hypoxamiRs networks will be instrumental both at a fundamental level and for their future potential therapeutic applications.

Inferring condition-specific miRNA activity from matched miRNA and mRNA expression data

MOTIVATION

MicroRNAs (miRNAs) play crucial roles in complex cellular networks by binding to the messenger RNAs (mRNAs) of protein coding genes. It has been found that miRNA regulation is often condition-specific. A number of computational approaches have been developed to identify miRNA activity specific to a condition of interest using gene expression data. However, most of the methods only use the data in a single condition, and thus, the activity discovered may not be unique to the condition of interest. Additionally, these methods are based on statistical associations between the gene expression levels of miRNAs and mRNAs, so they may not be able to reveal real gene regulatory relationships, which are causal relationships.

RESULTS

We propose a novel method to infer condition-specific miRNA activity by considering (i) the difference between the regulatory behavior that an miRNA has in the condition of interest and its behavior in the other conditions; (ii) the causal semantics of miRNA-mRNA relationships. The method is applied to the epithelial-mesenchymal transition (EMT) and multi-class cancer (MCC) datasets. The validation by the results of transfection experiments shows that our approach is effective in discovering significant miRNA-mRNA interactions. Functional and pathway analysis and literature validation indicate that the identified active miRNAs are closely associated with the specific biological processes, diseases and pathways. More detailed analysis of the activity of the active miRNAs implies that some active miRNAs show different regulation types in different conditions, but some have the same regulation types and their activity only differs in different conditions in the strengths of regulation.

AVAILABILITY AND IMPLEMENTATION

The R and Matlab scripts are in the Supplementary materials.

Investigating microRNA-target interaction-supported tissues in human cancer tissues based on miRNA and target gene expression profiling

Recent studies have revealed that a small non-coding RNA, microRNA (miRNA) down-regulates its mRNA targets. This effect is regarded as an important role in various biological processes. Many studies have been devoted to predicting miRNA-target interactions. These studies indicate that the interactions may only be functional in some specific tissues, which depend on the characteristics of an miRNA. No systematic methods have been established in the literature to investigate the correlation between miRNA-target interactions and tissue specificity through microarray data. In this study, we propose a method to investigate miRNA-target interaction-supported tissues, which is based on experimentally validated miRNA-target interactions. The tissue specificity results by our method are in accordance with the experimental results in the literature.

AVAILABILITY AND IMPLEMENTATION

Our analysis results are available at http://tsmti.mbc.nctu.edu.tw/ and http://www.stat.nctu.edu.tw/hwang/tsmti.html.

Computational prediction of microRNA genes

The computational identification of novel microRNA (miRNA) genes is a challenging task in bioinformatics. Massive amounts of data describing unknown functional RNA transcripts have to be analyzed for putative miRNA candidates with automated computational pipelines. Beyond those miRNAs that meet the classical definition, high-throughput sequencing techniques have revealed additional miRNA-like molecules that are derived by alternative biogenesis pathways. Exhaustive bioinformatics analyses on such data involve statistical issues as well as precise sequence and structure inspection not only of the functional mature part but also of the whole precursor sequence of the putative miRNA. Apart from a considerable amount of species-specific miRNAs, the majority of all those genes are conserved at least among closely related organisms. Some miRNAs, however, can be traced back to very early points in the evolution of eukaryotic species. Thus, the investigation of the conservation of newly found miRNA candidates comprises an important step in the computational annotation of miRNAs.Topics covered in this chapter include a review on the obvious problem of miRNA annotation and family definition, recommended pipelines of computational miRNA annotation or detection, and an overview of current computer tools for the prediction of miRNAs and their limitations. The chapter closes discussing how those bioinformatic approaches address the problem of faithful miRNA prediction and correct annotation.

MicroRNAs: The Missing Link in the Biology of Graft-Versus-Host Disease?

Graft-versus-host disease (GVHD) is still the major complication of allogeneic hematopoietic stem cell transplantation. Despite extensive studies in understanding the pathophysiology of GVHD, its pathogenesis remains unclear. Recently, important functions of microRNAs have been demonstrated in various autoimmune diseases and cancers such as psoriasis and lymphoma. This review highlights the need to investigate the role of microRNAs in GVHD and hypothesizes that microRNAs may be one of the missing links in our understanding of GVHD, with the potential for novel therapeutics.

BayMiR: inferring evidence for endogenous miRNA-induced gene repression from mRNA expression profiles

Background

Popular miRNA target prediction techniques use sequence features to determine the functional miRNA target sites. These techniques commonly ignore the cellular conditions in which miRNAs interact with their targets in vivo. Gene expression data are rich resources that can complement sequence features to take into account the context dependency of miRNAs.

Results

We introduce BayMiR, a new computational method, that predicts the functionality of potential miRNA target sites using the activity level of the miRNAs inferred from genome-wide mRNA expression profiles. We also found that mRNA expression variation can be used as another predictor of functional miRNA targets. We benchmarked BayMiR, the expression variation, Cometa, and the TargetScan “context scores” on two tasks: predicting independently validated miRNA targets and predicting the decrease in mRNA abundance in miRNA overexpression assays. BayMiR performed better than all other methods in both benchmarks and, surprisingly, the variation index performed better than Cometa and some individual determinants of the TargetScan context scores. Furthermore, BayMiR predicted miRNA target sets are more consistently annotated with GO and KEGG terms than similar sized random subsets of genes with conserved miRNA seed regions. BayMiR gives higher scores to target sites residing near the poly(A) tail which strongly favors mRNA degradation using poly(A) shortening. Our work also suggests that modeling multiplicative interactions among miRNAs is important to predict endogenous mRNA targets.

Conclusions

We develop a new computational method for predicting the target mRNAs of miRNAs. BayMiR applies a large number of mRNA expression profiles and successfully identifies the mRNA targets and miRNA activities without using miRNA expression data. The BayMiR package is publicly available and can be readily applied to any mRNA expression data sets.

ProteoMirExpress: inferring microRNA and protein-centered regulatory networks from high-throughput proteomic and mRNA expression data

MicroRNAs (miRNAs) regulate gene expression through translational repression and RNA degradation. Recently developed high-throughput proteomic methods measure gene expression changes at protein level and therefore can reveal the direct effects of miRNAs' translational repression. Here, we present a web server, ProteoMirExpress, that integrates proteomic and mRNA expression data together to infer miRNA-centered regulatory networks. With both types of high-throughput data from the users, ProteoMirExpress is able to discover not only miRNA targets that have decreased mRNA, but also subgroups of targets with suppressed proteins whose mRNAs are not significantly changed or with decreased mRNA whose proteins are not significantly changed, which are usually ignored by most current methods. Furthermore, both direct and indirect targets of miRNAs can be detected. Therefore, ProteoMirExpress provides more comprehensive miRNA-centered regulatory networks. We used several published data to assess the quality of our inferred networks and prove the value of our server. ProteoMirExpress is available online, with free access to academic users.

Genome-wide microRNA and messenger RNA profiling in rodent liver development implicates mir302b and mir20a in repressing transforming growth factor-beta signaling

MicroRNAs (miRNAs) are recently discovered small RNA molecules that regulate developmental processes, such as proliferation, differentiation, and apoptosis; however, the identity of miRNAs and their functions during liver development are largely unknown. Here we investigated the miRNA and gene expression profiles for embryonic day (E)8.5 endoderm, E14.5 Dlk1(+) liver cells (hepatoblasts), and adult liver by employing Illumina sequencing. We found that miRNAs were abundantly expressed at all three stages. Using K-means clustering analysis, 13 miRNA clusters with distinct temporal expression patterns were identified. mir302b, an endoderm-enriched miRNA, was identified as an miRNA whose predicted targets are expressed highly in E14.5 hepatoblasts but low in the endoderm. We validated the expression of mir302b in the endoderm by whole-mount in situ hybridization. Interestingly, mir20a, the most highly expressed miRNA in the endoderm library, was also predicted to regulate some of the same targets as mir302b. We found that through targeting Tgfbr2, mir302b and mir20a are able to regulate transforming growth factor beta (TGFβ) signal transduction. Moreover, mir302b can repress liver markers in an embryonic stem cell differentiation model. Collectively, we uncovered dynamic patterns of individual miRNAs during liver development, as well as miRNA networks that could be essential for the specification and differentiation of liver progenitors. (HEPATOLOGY 2013).

A new algorithm for integrated analysis of miRNA-mRNA interactions based on individual classification reveals insights into bladder cancer

Background

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression. It has been proposed that miRNAs play an important role in cancer development and progression. Their ability to affect multiple gene pathways by targeting various mRNAs makes them an interesting class of regulators.

Methodology/Principal Findings

We have developed an algorithm, Classification based Analysis of Paired Expression data of RNA (CAPE RNA), which is capable of identifying altered miRNA-mRNA regulation between tissues samples that assigns interaction states to each sample without preexisting stratification of groups. The distribution of the assigned interaction states compared to given experimental groups is used to assess the quality of a predicted interaction. We demonstrate the applicability of our approach by analyzing urothelial carcinoma and normal bladder tissue samples derived from 24 patients. Using our approach, normal and tumor tissue samples as well as different stages of tumor progression were successfully stratified. Also, our results suggest interesting differentially regulated miRNA-mRNA interactions associated with bladder tumor progression.

Conclusions/Significance

The need for tools that allow an integrative analysis of microRNA and mRNA expression data has been addressed. With this study, we provide an algorithm that emphasizes on the distribution of samples to rank differentially regulated miRNA-mRNA interactions. This is a new point of view compared to current approaches. From bootstrapping analysis, our ranking yields features that build strong classifiers. Further analysis reveals genes identified as differentially regulated by miRNAs to be enriched in cancer pathways, thus suggesting biologically interesting interactions.

The not-so-neutral role of microRNAs in neutrophil biology

The role of microRNAs (miRNAs) as fine-tuners of gene expression is now well established in most aspects of cellular biology. Critically, it is becoming apparent that characterization of miRNA regulation could further the understanding of elusive cellular processes. Here, I briefly review the current literature assessing the role of miRNAs in the modulation of neutrophil biology and discuss how the definition of such miRNA regulation could help in the better understanding of neutrophil function.

Working together: combinatorial regulation by microRNAs

MicroRNAs (miRNAs) negatively regulate gene expression level of mRNA post-transcriptionally. Deep sequencing and large-scale screening methods have yielded about 1,500 miRNA sequences in human. Each miRNA contains a seed sequence that is required, but not sufficient, for the correct matching with its targets. Recent technological advances make it possible to capture the miRNAs with their cognate mRNAs at the RISC complex. These experiments have revealed thousands of validated mRNA-miRNA pairing events. In the context of human stem cells, 90% of the identified transcripts appear to be paired with at least two different miRNAs.In this chapter, we present a comprehensive outline for a combinatorial regulation mode by miRNAs. Initially, we summarize the computational and experimental evidence that support a combined effect of multiple miRNAs. Then, we describe miRror2.0, a platform specifically convened to consider the likelihood of miRNAs cooperativity in view of the targets, tissues and cell lines. We show that results from miRror2.0 can be further refined by an iterative procedure, calls Psi-miRror that gauges the robustness of the regulation. We illustrate the combinatorial regulation projected onto graphs of human pathways and show that these pathways are amenable to disruption by a small set of miRNAs. Finally, we propose that miRNA combinatorial regulation is an attractive regulatory strategy not only at the level of single target, but also at the level of pathways and cellular homeostasis. The joint operation of miRNAs is a powerful means to overcome the low specificity inherent in each individual miRNA.

miRNA target enrichment analysis reveals directly active miRNAs in health and disease

microRNAs (miRNAs) are short non-coding regulatory RNA molecules. The activity of a miRNA in a biological process can often be reflected in the expression program that characterizes the outcome of the activity. We introduce a computational approach that infers such activity from high-throughput data using a novel statistical methodology, called minimum-mHG (mmHG), that examines mutual enrichment in two ranked lists. Based on this methodology, we provide a user-friendly web application that supports the statistical assessment of miRNA target enrichment analysis (miTEA) in the top of a ranked list of genes or proteins. Using miTEA, we analyze several target prediction tools by examining performance on public miRNA constitutive expression data. We also apply miTEA to analyze several integrative biology data sets, including a novel matched miRNA/mRNA data set covering nine human tissue types. Our novel findings include proposed direct activity of miR-519 in placenta, a direct activity of the oncogenic miR-15 in different healthy tissue types and a direct activity of the poorly characterized miR-768 in both healthy tissue types and cancer cell lines. The miTEA web application is available at http://cbl-gorilla.cs.technion.ac.il/miTEA/.

"Seed-Milarity" confers to hsa-miR-210 and hsa-miR-147b similar functional activity

Specificity of interaction between a microRNA (miRNA) and its targets crucially depends on the seed region located in its 5′-end. It is often implicitly considered that two miRNAs sharing the same biological activity should display similarity beyond the strict six nucleotide region that forms the seed, in order to form specific complexes with the same mRNA targets. We have found that expression of hsa-miR-147b and hsa-miR-210, though triggered by different stimuli (i.e. lipopolysaccharides and hypoxia, respectively), induce very similar cellular effects in term of proliferation, migration and apoptosis. Hsa-miR-147b only shares a “minimal” 6-nucleotides seed sequence with hsa-miR-210, but is identical with hsa-miR-147a over 20 nucleotides, except for one base located in the seed region. Phenotypic changes induced after heterologous expression of miR-147a strikingly differ from those induced by miR-147b or miR-210. In particular, miR-147a behaves as a potent inhibitor of cell proliferation and migration. These data fit well with the gene expression profiles observed for miR-147b and miR-210, which are very similar, and the gene expression profile of miR-147a, which is distinct from the two others. Bioinformatics analysis of all human miRNA sequences indicates multiple cases of miRNAs from distinct families exhibiting the same kind of similarity that would need to be further characterized in terms of putative functional redundancy. Besides, it implies that functional impact of some miRNAs can be masked by robust expression of miRNAs belonging to distinct families.

MicroRNA-125b down-regulates matrix metallopeptidase 13 and inhibits cutaneous squamous cell carcinoma cell proliferation, migration, and invasion

Cutaneous squamous cell carcinoma (cSCC) is the second most common human cancer. Although dysregulation of microRNAs (miRNAs) is known to be involved in a variety of cancers, the role of miRNAs in cSCC is unclear. In this study, we aimed to identify tumor suppressive and oncogenic miRNAs involved in the pathogenesis of cSCC. MiRNA expression profiles in healthy skins (n = 4) and cSCCs (n = 4) were analyzed using MicroRNA Low Density Array. MiR-125b expression was analyzed by quantitative real-time PCR and in situ hybridization in skin biopsies from 40 healthy donors, 13 actinic keratosis, and 74 cSCC patients. The effect of miR-125b was analyzed in wound closure, colony formation, migration, and invasion assays in two cSCC cell lines, UT-SCC-7 and A431. The genes regulated by miR-125b in cSCC were identified by microarray analysis and its direct target was validated by luciferase reporter assay. Comparing cSCC with healthy skin, we identified four up-regulated miRNAs (miR-31, miR-135b, miR-21, and miR-223) and 54 down-regulated miRNAs, including miR-125b, whose function was further examined. We found that miR-125b suppressed proliferation, colony formation, migratory, and invasive capacity of cSCC cells. Matrix metallopeptidase 13 (MMP13) was identified as a direct target suppressed by miR-125b, and there was an inverse relationship between the expression of miR-125b and MMP13 in cSCC. Knockdown of MMP13 expression phenocopied the effects of miR-125b overexpression. These findings provide a novel molecular mechanism by which MMP13 is up-regulated in cSCCs and indicate that miR-125b plays a tumor suppressive role in cSCC.

Induction of hepatocellular carcinoma by in vivo gene targeting

The distinct phenotypic and prognostic subclasses of human hepatocellular carcinoma (HCC) are difficult to reproduce in animal experiments. Here we have used in vivo gene targeting to insert an enhancer-promoter element at an imprinted chromosome 12 locus in mice, thereby converting ∼1 in 20,000 normal hepatocytes into a focus of HCC with a single genetic modification. A 300-kb chromosomal domain containing multiple mRNAs, snoRNAs, and microRNAs was activated surrounding the integration site. An identical domain was activated at the syntenic locus in a specific molecular subclass of spontaneous human HCCs with a similar histological phenotype, which was associated with partial loss of DNA methylation. These findings demonstrate the accuracy of in vivo gene targeting in modeling human cancer and suggest future applications in studying various tumors in diverse animal species. In addition, similar insertion events produced by randomly integrating vectors could be a concern for liver-directed human gene therapy.

MAGIA²: from miRNA and genes expression data integrative analysis to microRNA-transcription factor mixed regulatory circuits (2012 update)

MAGIA² (http://gencomp.bio.unipd.it/magia2) is an update, extension and evolution of the MAGIA web tool. It is dedicated to the integrated analysis of in silico target prediction, microRNA (miRNA) and gene expression data for the reconstruction of post-transcriptional regulatory networks. miRNAs are fundamental post-transcriptional regulators of several key biological and pathological processes. As miRNAs act prevalently through target degradation, their expression profiles are expected to be inversely correlated to those of the target genes. Low specificity of target prediction algorithms makes integration approaches an interesting solution for target prediction refinement. MAGIA² performs this integrative approach supporting different association measures, multiple organisms and almost all target predictions algorithms. Nevertheless, miRNAs activity should be viewed as part of a more complex scenario where regulatory elements and their interactors generate a highly connected network and where gene expression profiles are the result of different levels of regulation. The updated MAGIA² tries to dissect this complexity by reconstructing mixed regulatory circuits involving either miRNA or transcription factor (TF) as regulators. Two types of circuits are identified: (i) a TF that regulates both a miRNA and its target and (ii) a miRNA that regulates both a TF and its target.