A dictionary on microRNAs and their putative target pathways

While in the last decade mRNA expression profiling was among the most popular research areas, over the past years the study of non-coding RNAs, especially microRNAs (miRNAs), has gained increasing interest. For almost 900 known human miRNAs hundreds of pretended targets are known. However, there is only limited knowledge about putative systemic effects of changes in the expression of miRNAs and their regulatory influence. We determined for each known miRNA the biochemical pathways in the KEGG and TRANSPATH database and the Gene Ontology categories that are enriched with respect to its target genes. We refer to these pathways and categories as target pathways of the corresponding miRNA. Investigating target pathways of miRNAs we found a strong relation to disease-related regulatory pathways, including mitogen-activated protein kinase (MAPK) signaling cascade, Transforming growth factor (TGF)-beta signaling pathway or the p53 network. Performing a sophisticated analysis of differentially expressed genes of 13 cancer data sets extracted from gene expression omnibus (GEO) showed that targets of specific miRNAs were significantly deregulated in these sets. The respective miRNA target analysis is also a novel part of our gene set analysis pipeline GeneTrail. Our study represents a comprehensive theoretical analysis of the relationship between miRNAs and their predicted target pathways. Our target pathways analysis provides a ‘miRNA-target pathway’ dictionary, which enables researchers to identify target pathways of differentially regulated miRNAs.     

胃癌多药耐药相关microRNA的筛选和鉴定

目的:研究胃癌多药耐药相关microRNA并对其进行鉴定、靶基因预测和预测靶基因的生物信息学分析。方法:运用microRNA芯片对胃癌多药耐药细胞SGC7901/ADR和其亲本细胞SGC7901进行microRNA表达谱分析;采用实时定量PCR的方法对差异表达的miRNA进行验证;再运用生物信息学方法对差异表达的miRNA进行靶基因预测;再对预测的靶基因进行GO和KEGG通路分析。结果:与SGC7901相比SGC7901/ADR表达上调超过2倍的miRNA有6个,表达下调超过2倍的有11个。实时定量PCR对共同差异表达的microRNA进行验证显示与芯片结果的一致性。对这17个差异表达的miRNA进行靶基因预测,再对预测得到的靶基因进行GO和KEGG通路分析显示预测的靶基因参与了肿瘤相关通路、MAPK通路、Focal Adhesion通路等。结论:我们初步筛选得到了胃癌多药耐药相关miRNA并对其进行了生物信息学分析,为进一步地探索miRNA在胃癌多药耐药中的作用及其分子机制奠定了基础。     

Bioinformatic analyses of microRNA-targeted genes and microarray-identified genes correlated with Barrett's esophagus

Barrett's esophagus (BE) is defined as a metaplasia condition in the distal esophagus, in which the native squamous epithelium lining is replaced by a columnar epithelium with or without intestinal metaplasia. It is commonly accepted that BE is a precancerous lesion for esophageal adenocarcinoma. The aim of this study was to investigate the aberrant microRNAs (miRNAs) and differentially expressed genes (DEGs) associated with BE based on online microarray datasets. One miRNA and five gene expression profiling datasets were retrieved from the Gene Expression Omnibus Database. Aberrant microRNAs and DEGs were obtained using R/Bioconductor statistical analysis language and software. 23 dysregulated miRNAs and 632 DEGs demonstrating consistent expression tendencies in the five gene microarrays were identified in BE. Moreover, 1962 target genes of aberrant miRNAs were predicted using three bioinformatic tools, namely TargetScan, RNA22-HSA and miRDB. Ultimately, 93 target DEGs were obtained, after which functional annotation was performed on DAVID Bioinformatics Resources. Among Gene Ontology (GO) biological processes, digestive tract development and epithelial cell differentiation have demonstrated significant associations with BE pathogenesis. In addition, analysis of the KEGG pathways has revealed associations with cancer. To enable further study, one miRNA-target DEGs regulatory network was constructed using Cytoscape. 6 target DEGs demonstrated higher-degree distributions in the network, and ROC analysis indicated that FNDC3B may be the best potential biomarker for BE diagnosis. The data presented herein may provide new perspectives for exploring BE pathogenesis and may offer hits with regard to potential biomarkers in BE diagnosis, prediction and therapeutic evaluation.     

Investigation of key microRNAs associated with hepatocellular carcinoma using small RNA-seq data

To identify key microRNAs (miRNAs) associated with hepatocellular carcinoma (HCC) using small RNA-seq data. Small RNA-seq data for two HCC samples and two normal samples were downloaded from NCBI Gene Expression Omnibus. MiRNAs were identified through database search. Differentially expressed miRNAs were screened out with t test and their target genes were retrieved. Functional enrichment analysis was performed to uncover their biological functions. Regulatory networks and core metabolic networks were also constructed to present the global patterns. In addition, new miRNAs and their target genes were predicted. A total of 59 differentially expressed miRNAs were obtained, 12 up-regulated and 47 down-regulated. A total of 3,306 target genes were retrieved for eight miRNAs. Pathway enrichment analysis for the target genes showed that “pathways in cancer” and “MAPK signaling pathway” were significantly over-represented. Functional enrichment analysis found that “biological regulation” and “macromolecule modification” were significantly related to the target genes. Two regulatory networks were constructed for up- and down-regulated differentially expressed miRNAs with information from Ingenuity Pathway Analysis database. Two metabolic networks were also established based upon “pathways in cancer” and “MAPK signaling pathway”, consisting of miRNAs, target genes, compounds and others genes. Moreover, a number of new miRNAs and relevant target genes were predicted. Our study discloses a number of miRNAs as well as genes which may be involved in the development of HCC and these findings are beneficial in guiding future researches.     

Prognostic value of the tumor-specific ceRNA network in epithelial ovarian cancer

Epithelial ovarian cancer is one of the leading causes of cancer-related death worldwide. Growing evidence indicates that multiple complex altered pathways play important regulatory roles in the development and progression of a variety of cancers, including epithelial ovarian cancer. However, the underlying mechanisms remain unclear. First, we identified differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in epithelial ovarian cancer by comparing the expression profiles between epithelial ovarian cancer samples and normal tissue samples in different GEO datasets. Then, GO- and KEGG-pathway-enrichment analyses were applied to investigate the primary functions of the overlapped differentially expressed mRNAs. Moreover, the primary enriched genes were used to construct the signal-network with Cytoscape software. In addition, we integrated the relationship among lncRNAs-miRNAs-mRNAs to create a competing endogenous RNA network. Finally, mRNAs that were associated with patient prognosis in epithelial ovarian cancer were selected using univariate Cox regression analysis. A total of 2,225 mRNAs, 336 lncRNAs, and 14 miRNAs were shown to be differentially expressed in epithelial ovarian cancer compared with normal tissues. The dysregulated mRNAs were primarily enriched in cell division and signal transduction, according to Gene Ontology, whereas, according to KEGG, they were primarily enriched in metabolic pathways and pathways in cancer. A total of 10 mRNAs were associated with patient prognosis in ovarian cancer. This study identifies a novel lncRNA–miRNA–mRNA network, which may suggest potential molecular mechanisms underlying the development of epithelial ovarian cancer, providing new insights for survival prediction and interventional strategies for epithelial ovarian cancer.     

Retracted: Deep integrative analysis of microRNA-mRNA regulatory networks for biomarker and target discovery in chondrosarcoma

Chondrosarcoma (CHS) is a common malignant bone sarcoma and its occurrence increases with age. microRNAs (miRNAs) are a class of noncoding RNAs that participate in various biological processes and disease pathogenesis by targeting functional messenger RNA (mRNA). However, the modulation of miRNAs in CHS remains largely unknown. In this study, we performed integrative analysis to explore the expression profiles of miRNAs and mRNAs, together with their interaction networks in human CHS tissues and cell lines by RNA-seq (miRNA and mRNA). A total of 133 and 796 differentially expressed miRNAs and mRNAs were identified (|Fold change| ≥ 2 and P-value ≤ 0.5). miRNA-mRNA regulatory interactions between 55 miRNAs and 242 mRNAs were screened by the Pearson correlation analysis and target prediction. mRNAs in the network were enriched to 145 Gene Ontology terms and 35 Kyoto Encyclopedia of Genes and Genomes pathways. Specifically, some key regulators (hub-miRNAs) in the network (miR-622, miR-4539, miR-145, miR-25, and miR-96) were suggested to play important regulatory roles in the pathogenesis of CHS. In addition, functional experiments validated that miR-622 regulated CHS cell proliferation by targeting bone morphogenetic protein 1 (BMP1).     

miRNAs may play a major role in the control of gene expression in key pathobiological processes in Chagas disease cardiomyopathy

Chronic Chagas disease cardiomyopathy (CCC), an especially aggressive inflammatory dilated cardiomyopathy caused by lifelong infection with the protozoan Trypanosoma cruzi, is a major cause of cardiomyopathy in Latin America. Although chronic myocarditis may play a major pathogenetic role, little is known about the molecular mechanisms responsible for its severity. The aim of this study is to study the genes and microRNAs expression in tissues and their connections in regards to the pathobiological processes. To do so, we integrated for the first time global microRNA and mRNA expression profiling from myocardial tissue of CCC patients employing pathways and network analyses. We observed an enrichment in biological processes and pathways associated with the immune response and metabolism. IFNγ, TNF and NFkB were the top upstream regulators. The intersections between differentially expressed microRNAs and differentially expressed target mRNAs showed an enrichment in biological processes such as Inflammation, inflammation, Th1/IFN-γ-inducible genes, fibrosis, hypertrophy, and mitochondrial/oxidative stress/antioxidant response. MicroRNAs also played a role in the regulation of gene expression involved in the key cardiomyopathy-related processes fibrosis, hypertrophy, myocarditis and arrhythmia. Significantly, a discrete number of differentially expressed microRNAs targeted a high number of differentially expressed mRNAs (>20) in multiple processes. Our results suggest that miRNAs orchestrate expression of multiple genes in the major pathophysiological processes in CCC heart tissue. This may have a bearing on pathogenesis, biomarkers and therapy.     

EpCAM Knockdown Alters MicroRNA Expression in Retinoblastoma- Functional Implication of EpCAM Regulated MiRNA in Tumor Progression

The co-ordinated regulation of oncogenes along with miRNAs play crucial role in carcinogenesis. In retinoblastoma (RB), several miRNAs are known to be differentially expressed. Epithelial cell adhesion molecule (EpCAM) gene is involved in many epithelial cancers including, retinoblastoma (RB) tumorigenesis. EpCAM silencing effectively reduces the oncogenic miR-17-92 cluster. In order to investigate whether EpCAM has wider effect as an inducer or silencer of miRNAs, we performed a global microRNA expression profile in EpCAM siRNA knockdown Y79 cells. MicroRNA profiling in EpCAM silenced Y79 cells showed seventy-three significantly up regulated and thirty-six down regulated miRNAs. A subset of these miRNAs was also validated in tumors. Functional studies on Y79 and WERI-Rb-1 cells transfected with antagomirs against two miRNAs of miR-181c and miR-130b showed striking changes in tumor cell properties in RB cells. Treatment with anti-miR-181c and miR-130b showed significant decrease in cell viability and cell invasion. Increase in caspase-3 level was noticed in antagomir transfected cell lines indicating the induction of apoptosis. Possible genes altered by EpCAM influenced microRNAs were predicted by bioinformatic tools. Many of these belong to pathways implicated in cancer. The study shows significant influence of EpCAM on global microRNA expression. EpCAM regulated miR-181c and miR-130b may play significant roles in RB progression. EpCAM based targeted therapies may reduce carcinogenesis through several miRNAs and target genes.     

miR-425-5p as an exosomal biomarker for metastatic prostate cancer

Prostate cancer-related deaths are mostly caused by metastasis, which indicates the importance of identifying clinical prognostic biomarkers. In this study, we evaluated the expression profile of exosomal microRNAs (miRNAs) derived from metastatic prostate cancer (mPCa) cell lines (LNCaP and PC-3). miRNA signatures in exosomes and cells were evaluated by miRNA microarray analysis. Fourteen miRNAs were identified as candidates for specific noninvasive biomarkers. The expression of five miRNAs was validated using RT-qPCR, which confirmed that miR-205-5p, miR-148a-3p, miR-125b-5p, miR-183-5p, and miR-425-5p were differentially expressed in mPCa exosomes. Bioinformatic analyses showed that miR-425-5p was associated with residual tumor, pathologic T and N stages, and TP53 status in PCa samples. Gene ontology analysis of negatively correlated and predicted targeted genes showed enrichment of genes related to bone development pathways. The LinkedOmics database indicated that the potential target HSPB8 has a significant negative correlation with miR-425-5p. In conclusion, this study identified a panel of exosomal miRNAs with potential value as prognostic biomarkers for prostate cancer.     

应用高通量测序技术检测与奶牛乳产量和乳质量调控相关的miRNA

本实验将中国荷斯坦牛泌乳期高乳品质奶牛（H）和泌乳期低乳品质奶牛（L）乳腺组织作为实验对象,利用高通量测序技术进行了miRNA测序,与miRNA数据库比对,获得已知miRNA,整合miREvo和mirDeep2这两个miRNA预测软件,进行新miRNA分析,通过差异表达分析筛选组间差异miRNAs,获得56个差异表达miRNA(P <0.05,FDRq <0.05)并对差异表达miRNA进行靶基因预测;利用DAVID对靶基因进行GO(Gene Ontology)和信号通路富集分析。经过对靶基因筛选,发现了4个已报道与乳蛋白、乳脂紧密相关的功能基因：CSN3、SCD、LALBA和DGAT2。靶基因聚集的生物学功能多数参与了蛋白质和脂肪代谢,乳腺发育和分化,以及免疫功能。靶基因主要富集在MAPK 信号通路、甘油磷酸脂质代谢、缺氧诱导因子1和磷脂酰肌醇3激酶 蛋白激酶B信号转导通路。结果显示,靶基因主要富集在糖类代谢、脂肪代谢、蛋白质代谢、细胞凋亡以及免疫相关通路。     

Novel transcriptome data analysis implicates circulating microRNAs in epigenetic inheritance in mammals

Experimental evidence supports a role of mobile small non-coding RNAs in mediating soma to germline hereditary information transfer in epigenetic inheritance in plants and worms. Similar evidence in mammals has not been reported so far. In this bioinformatic analysis, differentially expressed microRNAs (miRNAs) or mRNAs reported previously in genome level expression profiling studies related to or relevant in epigenetic inheritance in mammals were examined for circulating miRNA association. The reported sets of differentially expressed miRNAs or miRNAs that are known to target the reported sets of differentially expressed genes, in that order, showed enrichment of circulating miRNAs across environmental factors, tissues, life cycle stages, generations, genders and species. Circulating miRNAs commonly representing the expression profiles enriched various epigenetic processes. These results provide bioinformatic evidence for a role of circulating miRNAs in epigenetic inheritance in mammals.     

Identification and characterization of microRNAs and their target genes in Brassica oleracea

The microRNAs are a new class of small non-coding endogenous RNAs with lengths of approximately ~21 nt. MicroRNAs perform their biological function via the degradation of the target mRNAs or by inhibiting protein translation. Until recently, only limited numbers of miRNAs were identified in Brassica oleracea, a vegetable widely cultivated around the world. In present study, 193 potential miRNA candidates were identified from 17 expressed sequence tag (ESTs) and 152 genome survey sequences (GSSs) in B. oleracea. These miRNA candidates were classified into 70 families using a well-defined comparative genome-based computational analysis. Most miRNAs belong to the miRNA169, miR5021, miR156 and miR158 families. Of these, 36 miRNA families are firstly found in Brassica species. Around 1393 B. oleracea genes were predicted as candidate targets of 175 miRNAs. The mutual relationship between miRNAs and the candidate target genes was verified by checking differentially expression levels using quantitative real-time polymerase chain reaction (qRT-PCR) and 5' RLM-RACE analyses. These target genes participate in multiple biological and metabolic processes, including signal transduction, stress response, and plant development. Gene Ontology analysis shows that the 818, 514, and 265 target genes are involved in molecular functions, biological processes, and cellular component respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis suggests that these miRNAs might regulate 186 metabolic pathways, including those of lipid, energy, starch and sucrose, fatty acid and nitrogen.     

De-Regulated MicroRNAs in Pediatric Cancer Stem Cells Target Pathways Involved in Cell Proliferation,Cell Cycle and Development

Background

microRNAs (miRNAs) have been implicated in the control of many biological processes and their deregulation has been associated with many cancers. In recent years, the cancer stem cell (CSC) concept has been applied to many cancers including pediatric. We hypothesized that a common signature of deregulated miRNAs in the CSCs fraction may explain the disrupted signaling pathways in CSCs.

Methodology/Results

Using a high throughput qPCR approach we identified 26 CSC associated differentially expressed miRNAs (DEmiRs). Using BCmicrO algorithm 865 potential CSC associated DEmiR targets were obtained. These potential targets were subjected to KEGG, Biocarta and Gene Ontology pathway and biological processes analysis. Four annotated pathways were enriched: cell cycle, cell proliferation, p53 and TGF-beta/BMP. Knocking down hsa-miR-21-5p, hsa-miR-181c-5p and hsa-miR-135b-5p using antisense oligonucleotides and small interfering RNA in cell lines led to the depletion of the CSC fraction and impairment of sphere formation (CSC surrogate assays).

Conclusion

Our findings indicated that CSC associated DEmiRs and the putative pathways they regulate may have potential therapeutic applications in pediatric cancers.     

Mirna expression profiles identify drivers in colorectal and pancreatic cancers

Background and Aim

Altered expression of microRNAs (miRNAs) hallmarks many cancer types. The study of the associations of miRNA expression profile and cancer phenotype could help identify the links between deregulation of miRNA expression and oncogenic pathways.

Methods

Expression profiling of 866 human miRNAs in 19 colorectal and 17 pancreatic cancers and in matched adjacent normal tissues was investigated. Classical paired t-test and random forest analyses were applied to identify miRNAs associated with tissue-specific tumors. Network analysis based on a computational approach to mine associations between cancer types and miRNAs was performed.

Results

The merge between the two statistical methods used to intersect the miRNAs differentially expressed in colon and pancreatic cancers allowed the identification of cancer-specific miRNA alterations. By miRNA-network analysis, tissue-specific patterns of miRNA deregulation were traced: the driving miRNAs were miR-195, miR-1280, miR-140-3p and miR-1246 in colorectal tumors, and miR-103, miR-23a and miR-15b in pancreatic cancers.

Conclusion

MiRNA expression profiles may identify cancer-specific signatures and potentially useful biomarkers for the diagnosis of tissue specific cancers. miRNA-network analysis help identify altered miRNA regulatory networks that could play a role in tumor pathogenesis.