Inferring the perturbed microRNA regulatory networks from gene expression data using a network propagation based method

Background

MicroRNAs (miRNAs) are a class of endogenous small regulatory RNAs. Identifications of the dys-regulated or perturbed miRNAs and their key target genes are important for understanding the regulatory networks associated with the studied cellular processes. Several computational methods have been developed to infer the perturbed miRNA regulatory networks by integrating genome-wide gene expression data and sequence-based miRNA-target predictions. However, most of them only use the expression information of the miRNA direct targets, rarely considering the secondary effects of miRNA perturbation on the global gene regulatory networks.

Results

We proposed a network propagation based method to infer the perturbed miRNAs and their key target genes by integrating gene expressions and global gene regulatory network information. The method used random walk with restart in gene regulatory networks to model the network effects of the miRNA perturbation. Then, it evaluated the significance of the correlation between the network effects of the miRNA perturbation and the gene differential expression levels with a forward searching strategy. Results show that our method outperformed several compared methods in rediscovering the experimentally perturbed miRNAs in cancer cell lines. Then, we applied it on a gene expression dataset of colorectal cancer clinical patient samples and inferred the perturbed miRNA regulatory networks of colorectal cancer, including several known oncogenic or tumor-suppressive miRNAs, such as miR-17, miR-26 and miR-145.

Conclusions

Our network propagation based method takes advantage of the network effect of the miRNA perturbation on its target genes. It is a useful approach to infer the perturbed miRNAs and their key target genes associated with the studied biological processes using gene expression data.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-255) contains supplementary material, which is available to authorized users.     

MiR-107 and MiR-185 Can Induce Cell Cycle Arrest in Human Non Small Cell Lung Cancer Cell Lines

Background

MicroRNAs (miRNAs) are short single stranded noncoding RNAs that suppress gene expression through either translational repression or degradation of target mRNAs. The annealing between messenger RNAs and 5′ seed region of miRNAs is believed to be essential for the specific suppression of target gene expression. One miRNA can have several hundred different targets in a cell. Rapidly accumulating evidence suggests that many miRNAs are involved in cell cycle regulation and consequentially play critical roles in carcinogenesis.

Methodology/Principal Findings

Introduction of synthetic miR-107 or miR-185 suppressed growth of the human non-small cell lung cancer cell lines. Flow cytometry analysis revealed these miRNAs induce a G1 cell cycle arrest in H1299 cells and the suppression of cell cycle progression is stronger than that by Let-7 miRNA. By the gene expression analyses with oligonucleotide microarrays, we find hundreds of genes are affected by transfection of these miRNAs. Using miRNA-target prediction analyses and the array data, we listed up a set of likely targets of miR-107 and miR-185 for G1 cell cycle arrest and validate a subset of them using real-time RT-PCR and immunoblotting for CDK6.

Conclusions/Significance

We identified new cell cycle regulating miRNAs, miR-107 and miR-185, localized in frequently altered chromosomal regions in human lung cancers. Especially for miR-107, a large number of down-regulated genes are annotated with the gene ontology term ‘cell cycle’. Our results suggest that these miRNAs may contribute to regulate cell cycle in human malignant tumors.     

Integrated analyses of zebrafish miRNA and mRNA expression profiles identify miR-29b and miR-223 as potential regulators of optic nerve regeneration

Background

Unlike mammals, zebrafish have the ability to regenerate damaged parts of their central nervous system (CNS) and regain functionality of the affected area. A better understanding of the molecular mechanisms involved in zebrafish regeneration may therefore provide insight into how CNS repair might be induced in mammals. Although many studies have described differences in gene expression in zebrafish during CNS regeneration, the regulatory mechanisms underpinning the differential expression of these genes have not been examined.

Results

We used microarrays to analyse and integrate the mRNA and microRNA (miRNA) expression profiles of zebrafish retina after optic nerve crush to identify potential regulatory mechanisms that underpin central nerve regeneration. Bioinformatic analysis identified 3 miRNAs and 657 mRNAs that were differentially expressed after injury. We then combined inverse correlations between our miRNA expression and mRNA expression, and integrated these findings with target predictions from TargetScan Fish to identify putative miRNA-gene target pairs. We focused on two over-expressed miRNAs (miR-29b and miR-223), and functionally validated seven of their predicted gene targets using RT-qPCR and luciferase assays to confirm miRNA-mRNA binding. Gene ontology analysis placed the miRNA-regulated genes (eva1a, layna, nefmb, ina, si:ch211-51a6.2, smoc1, sb:cb252) in key biological processes that included cell survival/apoptosis, ECM-cytoskeleton signaling, and heparan sulfate proteoglycan binding,

Conclusion

Our results suggest a key role for miR-29b and miR-223 in zebrafish regeneration. The identification of miRNA regulation in a zebrafish injury model provides a framework for future studies in which to investigate not only the cellular processes required for CNS regeneration, but also how these mechanisms might be regulated to promote successful repair and return of function in the injured mammalian brain.

Electronic supplementary material

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

miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors

Introduction

Few studies have performed expression profiling of both miRNA and mRNA from the same primary breast carcinomas. In this study we present and analyze data derived from expression profiling of 799 miRNAs in 101 primary human breast tumors, along with genome-wide mRNA profiles and extensive clinical information.

Methods

We investigate the relationship between these molecular components, in terms of their correlation with each other and with clinical characteristics. We use a systems biology approach to examine the correlative relationship between miRNA and mRNAs using statistical enrichment methods.

Results

We identify statistical significant differential expression of miRNAs between molecular intrinsic subtypes, and between samples with different levels of proliferation. Specifically, we point to miRNAs significantly associated with TP53 and ER status. We also show that several cellular processes, such as proliferation, cell adhesion and immune response, are strongly associated with certain miRNAs. We validate the role of miRNAs in regulating proliferation using high-throughput lysate-microarrays on cell lines and point to potential drivers of this process.

Conclusion

This study provides a comprehensive dataset as well as methods and system-level results that jointly form a basis for further work on understanding the role of miRNA in primary breast cancer.     

Expression and localization of microRNAs in perinatal rat pancreas: role of miR-21 in regulation of cholesterol metabolism

Objective

To investigate the expression of pancreatic microRNAs (miRNAs) during the period of perinatal beta-cell expansion and maturation in rats, determine the localization of these miRNAs and perform a pathway analysis with predicted target mRNAs expressed in perinatal pancreas.

Research Design and Methods

RNA was extracted from whole pancreas at embryonic day 20 (E20), on the day of birth (P0) and two days after birth (P2) and hybridized to miRNA microarrays. Differentially expressed miRNAs were verified by northern blotting and their pancreatic localization determined by in situ hybridization. Pathway analysis was done using regulated sets of mRNAs predicted as targets of the miRNAs. Possible target genes were tested using reporter-gene analysis in INS-1E cells.

Results

Nine miRNAs were differentially expressed perinatally, seven were confirmed to be regulated at the level of the mature miRNA. The localization studies showed endocrine localization of six of these miRNAs (miR-21, -23a, -29a, -125b-5p, -376b-3p and -451), and all were expressed in exocrine cells at one time point at least. Pathways involving metabolic processes, terpenoid and sterol metabolism were selectively affected by concomitant regulation by miRNAs and mRNAs, and Srebf1 was validated as a target of miR-21.

Conclusions

The findings suggest that miRNAs are involved in the functional maturation of pancreatic exocrine and endocrine tissue following birth. Pathway analysis of target genes identify changes in sterol metabolism around birth as being selectively affected by differential miRNA expression during this period.     

Integrated miRNA-mRNA Analysis Revealing the Potential Roles of miRNAs in Chordomas

Introduction

Emerging evidence suggests that microRNAs (miRNAs) are crucially involved in tumorigenesis and that paired expression profiles of miRNAs and mRNAs can be used to identify functional miRNA-target relationships with high precision. However, no studies have applied integrated analysis to miRNA and mRNA profiles in chordomas. The purpose of this study was to provide insights into the pathogenesis of chordomas by using this integrated analysis method.

Methods

Differentially expressed miRNAs and mRNAs of chordomas (n = 3) and notochord tissues (n = 3) were analyzed by using microarrays with hierarchical clustering analysis. Subsequently, the target genes of the differentially expressed miRNAs were predicted and overlapped with the differentially expressed mRNAs. Then, GO and pathway analyses were performed for the intersecting genes.

Results

The microarray analysis indicated that 33 miRNAs and 2,791 mRNAs were significantly dysregulated between the two groups. Among the 2,791 mRNAs, 911 overlapped with putative miRNA target genes. A pathway analysis showed that the MAPK pathway was consistently enriched in the chordoma tissue and that miR-149-3p, miR-663a, miR-1908, miR-2861 and miR-3185 likely play important roles in the regulation of MAPK pathways. Furthermore, the Notch signaling pathway and the loss of the calcification or ossification capacity of the notochord may also be involved in chordoma pathogenesis.

Conclusion

This study provides an integrated dataset of the miRNA and mRNA profiles in chordomas, and the results demonstrate that not only the MAPK pathway and its related miRNAs but also the Notch pathway may be involved in chordoma development. The occurrence of chordoma may be associated with dysfunctional calcification or ossification of the notochord.     

Integrated miRNA and mRNA Expression Profiling in Inflamed Colon of Patients with Ulcerative Colitis

Background

Ulcerative colitis (UC) is associated with differential colonic expression of genes involved in immune response (e.g. IL8) and barrier integrity (e.g. cadherins). MicroRNAs (miRNAs) are regulators of gene expression and are involved in various immune-related diseases. In this study, we investigated (1) if miRNA expression in UC mucosa is altered and (2) if any of these changes correlate with mucosal mRNA expression. Integration of mRNA and miRNA expression profiling may allow the identification of functional links between dysregulated miRNAs and their target mRNA.

Methodology

Colonic mucosal biopsies were obtained from 17 UC (10 active and 7 inactive) patients and 10 normal controls. Total RNA was used to analyze miRNA and mRNA expression via Affymetrix miRNA 2.0 and Affymetrix Human Gene 1.0ST arrays, respectively. Both miRNA and gene expression profiles were integrated by correlation analysis to identify dysregulated miRNAs with their corresponding predicted target mRNA. Microarray data were validated with qRT-PCR. Regulation of IL8 and CDH11 expression by hsa-miR-200c-3p was determined by luciferase reporter assays.

Results

When comparing active UC patients vs. controls, 51 miRNAs and 1543 gene probe sets gave significantly different signals. In contrast, in inactive UC vs. controls, no significant miRNA expression differences were found while 155 gene probe sets had significantly different signals. We then identified potential target genes of the significantly dysregulated miRNAs and genes in active UC vs. controls and found a highly significant inverse correlation between hsa-miR-200c-3p and IL8, an inflammatory marker, and between hsa-miR-200c-3p and CDH11, a gene related to intestinal epithelial barrier function. We could demonstrate that hsa-miR-200c-3p directly regulates IL8 and CDH11 expression.

Conclusion

Differential expression of immune- and barrier-related genes in inflamed UC mucosa may be influenced by altered expression of miRNAs. Integrated analysis of miRNA and mRNA expression profiles revealed hsa-miR-200c-3p for use of miRNA mimics as therapeutics.     

Comparative mRNA and microRNA expression profiling of three genitourinary cancers reveals common hallmarks and cancer-specific molecular events

Background

Genome-wide gene expression profile using deep sequencing technologies can drive the discovery of cancer biomarkers and therapeutic targets. Such efforts are often limited to profiling the expression signature of either mRNA or microRNA (miRNA) in a single type of cancer.

Methodology

Here we provided an integrated analysis of the genome-wide mRNA and miRNA expression profiles of three different genitourinary cancers: carcinomas of the bladder, kidney and testis.

Principal Findings

Our results highlight the general or cancer-specific roles of several genes and miRNAs that may serve as candidate oncogenes or suppressors of tumor development. Further comparative analyses at the systems level revealed that significant aberrations of the cell adhesion process, p53 signaling, calcium signaling, the ECM-receptor and cell cycle pathways, the DNA repair and replication processes and the immune and inflammatory response processes were the common hallmarks of human cancers. Gene sets showing testicular cancer-specific deregulation patterns were mainly implicated in processes related to male reproductive function, and general disruptions of multiple metabolic pathways and processes related to cell migration were the characteristic molecular events for renal and bladder cancer, respectively. Furthermore, we also demonstrated that tumors with the same histological origins and genes with similar functions tended to group together in a clustering analysis. By assessing the correlation between the expression of each miRNA and its targets, we determined that deregulation of ‘key’ miRNAs may result in the global aberration of one or more pathways or processes as a whole.

Conclusions

This systematic analysis deciphered the molecular phenotypes of three genitourinary cancers and investigated their variations at the miRNA level simultaneously. Our results provided a valuable source for future studies and highlighted some promising genes, miRNAs, pathways and processes that may be useful for diagnostic or therapeutic applications.