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.     

Distribution of miRNA genes in the pig genome

Background

Recent completion of swine genome may simplify the production of swine as a large biomedical model. Here we studied sequence and location of known swine miRNA genes, key regulators of protein-coding genes at the level of RNA, and compared them to human and mouse data to prioritize future molecular studies.

Results

Distribution of miRNA genes in pig genome shows no particular relation to different genomic features including protein coding genes - proportions of miRNA genes in intergenic regions, introns and exons roughly agree with the size of these regions in the pig genome. Our analyses indicate that host genes harbouring intragenic miRNAs are longer from other protein-coding genes, however, no important GO enrichment was found. Swine mature miRNAs show high sequence similarity to their human and mouse orthologues. Location of miRNA genes relative to protein-coding genes is also similar among studied species, however, there are differences in the precise position in particular intergenic regions and within particular hosts. The most prominent difference between pig and human miRNAs is a large group of pig-specific sequences (53% of swine miRNAs). We found no evidence that this group of evolutionary new pig miRNAs is different from old miRNAs genes with respect to genomic location except that they are less likely to be clustered.

Conclusions

There are differences in precise location of orthologues miRNA genes in particular intergenic regions and within particular hosts, and their meaning for coexpression with protein-coding genes deserves experimental studies. Functional studies of a large group of pig-specific sequences in future may reveal limits of the pig as a model organism to study human gene expression.

Electronic supplementary material

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

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.     

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing

Background

MiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.

Results

We identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.

Conclusions

We provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N₂-fixing symbiotic nodules.

Electronic supplementary material

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

Microarray-Based Approach Identifies Differentially Expressed MicroRNAs in Porcine Sexually Immature and Mature Testes

Background

MicroRNAs (miRNAs) are short non-coding RNA molecules which are proved to be involved in mammalian spermatogenesis. Their expression and function in the porcine germ cells are not fully understood.

Methodology

We employed a miRNA microarray containing 1260 unique miRNA probes to evaluate the miRNA expression patterns between sexually immature (60-day) and mature (180-day) pig testes. One hundred and twenty nine miRNAs representing 164 reporter miRNAs were expressed differently (p<0.1). Fifty one miRNAs were significantly up-regulated and 78 miRNAs were down-regulated in mature testes. Nine of these differentially expressed miRNAs were validated using quantitative RT-PCR assay. Totally 15919 putative miRNA-target sites were detected by using RNA22 method to align 445 NCBI pig cDNA sequences with these 129 differentially expressed miRNAs, and seven putative target genes involved in spermatogenesis including DAZL, RNF4 gene were simply confirmed by quantitative RT-PCR.

Conclusions

Overall, the results of this study indicated specific miRNAs expression in porcine testes and suggested that miRNAs had a role in regulating spermatogenesis.     

miR-503 represses human cell proliferation and directly targets the oncogene DDHD2 by non-canonical target pairing

Background

The pathways regulating the transition of mammalian cells from quiescence to proliferation are mediated by multiple miRNAs. Despite significant improvements in our understanding of miRNA targeting, the majority of miRNA regulatory networks are still largely unknown and require experimental validation.

Results

Here we identified miR-503, miR-103, and miR-494 as negative regulators of proliferation in primary human cells. We experimentally determined their genome wide target profiles using RNA-induced silencing complex (RISC) immunoprecipitations and gene expression profiling. Analysis of the genome wide target profiles revealed evidence of extensive regulation of gene expression through non-canonical target pairing by miR-503. We identified the proto-oncogene DDHD2 as a target of miR-503 that requires pairing outside of the canonical 5′ seed region of miR-503, representing a novel mode of miRNA-target pairing. Further bioinformatics analysis implicated miR-503 and DDHD2 in breast cancer tumorigenesis.

Conclusions

Our results provide an extensive genome wide set of targets for miR-503, miR-103, and miR-494, and suggest that miR-503 may act as a tumor suppressor in breast cancer by its direct non-canonical targeting of DDHD2.

Electronic supplementary material

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

Genome-wide identification of vegetative phase transition-associated microRNAs and target predictions using degradome sequencing in Malus hupehensis

Background

A long juvenile period between germination and flowering is a common characteristic among fruit trees, including Malus hupehensis (Pamp.) Rehd., which is an apple rootstock widely used in China. microRNAs (miRNAs) play an important role in the regulation of phase transition and reproductive growth processes.

Results

M. hupehensis RNA libraries, one adult and one juvenile phase, were constructed using tree leaves and underwent high-throughput sequencing. We identified 42 known miRNA families and 172 novel miRNAs. We also identified 127 targets for 25 known miRNA families and 168 targets for 35 unique novel miRNAs using degradome sequencing. The identified miRNA targets were categorized into 58 biological processes, and the 123 targets of known miRNAs were associated with phase transition processes. The KEGG analysis revealed that these targets were involved in starch and sucrose metabolism, and plant hormone signal transduction. Expression profiling of miRNAs and their targets indicated multiple regulatory functions in the phase transition. The higher expression level of mdm-miR156 and lower expression level of mdm-miR172 in the juvenile phase leaves implied that these two small miRNAs regulated the phase transition. mdm-miR160 and miRNA393, which regulate genes involved in auxin signal transduction, could also be involved in controlling this process. The identification of known and novel miRNAs and their targets provides new information on this regulatory process in M. hupehensis, which will contribute to the understanding of miRNA functions during growth, phase transition and reproduction in woody fruit trees.

Conclusions

The combination of sRNA and degradome sequencing can be used to better illustrate the profiling of hormone-regulated miRNAs and miRNA targets involving complex regulatory networks, which will contribute to the understanding of miRNA functions during growth, phase transition and reproductive growth in perennial woody fruit trees.

Electronic supplementary material

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

An integrative approach to identify hexaploid wheat miRNAome associated with development and tolerance to abiotic stress

Background

Wheat is a major staple crop with broad adaptability to a wide range of environmental conditions. This adaptability involves several stress and developmentally responsive genes, in which microRNAs (miRNAs) have emerged as important regulatory factors. However, the currently used approaches to identify miRNAs in this polyploid complex system focus on conserved and highly expressed miRNAs avoiding regularly those that are often lineage-specific, condition-specific, or appeared recently in evolution. In addition, many environmental and biological factors affecting miRNA expression were not yet considered, resulting still in an incomplete repertoire of wheat miRNAs.

Results

We developed a conservation-independent technique based on an integrative approach that combines machine learning, bioinformatic tools, biological insights of known miRNA expression profiles and universal criteria of plant miRNAs to identify miRNAs with more confidence. The developed pipeline can potentially identify novel wheat miRNAs that share features common to several species or that are species specific or clade specific. It allowed the discovery of 199 miRNA candidates associated with different abiotic stresses and development stages. We also highlight from the raw data 267 miRNAs conserved with 43 miRBase families. The predicted miRNAs are highly associated with abiotic stress responses, tolerance and development. GO enrichment analysis showed that they may play biological and physiological roles associated with cold, salt and aluminum (Al) through auxin signaling pathways, regulation of gene expression, ubiquitination, transport, carbohydrates, gibberellins, lipid, glutathione and secondary metabolism, photosynthesis, as well as floral transition and flowering.

Conclusion

This approach provides a broad repertoire of hexaploid wheat miRNAs associated with abiotic stress responses, tolerance and development. These valuable resources of expressed wheat miRNAs will help in elucidating the regulatory mechanisms involved in freezing and Al responses and tolerance mechanisms as well as for development and flowering. In the long term, it may help in breeding stress tolerant plants.

Electronic supplementary material

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

Role of MicroRNA Modulation in the Interferon-α/Ribavirin Suppression of HIV-1 In Vivo

Background

Interferon-α (IFN-α) treatment suppresses HIV-1 viremia and reduces the size of the HIV-1 latent reservoir. Therefore, investigation of the molecular and immunologic effects of IFN-α may provide insights that contribute to the development of novel prophylactic, therapeutic and curative strategies for HIV-1 infection. In this study, we hypothesized that microRNAs (miRNAs) contribute to the IFN-α-mediated suppression of HIV-1. To inform the development of novel miRNA-based antiretroviral strategies, we investigated the effects of exogenous IFN-α treatment on global miRNA expression profile, HIV-1 viremia, and potential regulatory networks between miRNAs and cell-intrinsic anti-HIV-1 host factors in vivo.

Methods

Global miRNA expression was examined in longitudinal PBMC samples obtained from seven HIV/HCV-coinfected, antiretroviral therapy-naïve individuals before, during, and after pegylated interferon-α/ribavirin therapy (IFN-α/RBV). We implemented novel hybrid computational-empirical approaches to characterize regulatory networks between miRNAs and anti-HIV-1 host restriction factors.

Results

miR-422a was the only miRNA significantly modulated by IFN-α/RBV in vivo (p<0.0001, paired t test; FDR<0.037). Our interactome mapping revealed extensive regulatory involvement of miR-422a in p53-dependent apoptotic and pyroptotic pathways. Based on sequence homology and inverse expression relationships, 29 unique miRNAs may regulate anti-HIV-1 restriction factor expression in vivo.

Conclusions

The specific reduction of miR-422a is associated with exogenous IFN-α treatment, and likely contributes to the IFN-α suppression of HIV-1 through the enhancement of anti-HIV-1 restriction factor expression and regulation of genes involved in programmed cell death. Moreover, our regulatory network analysis presents additional candidate miRNAs that may be targeted to enhance anti-HIV-1 restriction factor expression in vivo.     

Layered Signaling Regulatory Networks Analysis of Gene Expression Involved in Malignant Tumorigenesis of Non-Resolving Ulcerative Colitis via Integration of Cross-Study Microarray Profiles

Background

Ulcerative colitis (UC) was the most frequently diagnosed inflammatory bowel disease (IBD) and closely linked to colorectal carcinogenesis. By far, the underlying mechanisms associated with the disease are still unclear. With the increasing accumulation of microarray gene expression profiles, it is profitable to gain a systematic perspective based on gene regulatory networks to better elucidate the roles of genes associated with disorders. However, a major challenge for microarray data analysis is the integration of multiple-studies generated by different groups.

Methodology/Principal Findings

In this study, firstly, we modeled a signaling regulatory network associated with colorectal cancer (CRC) initiation via integration of cross-study microarray expression data sets using Empirical Bayes (EB) algorithm. Secondly, a manually curated human cancer signaling map was established via comprehensive retrieval of the publicly available repositories. Finally, the co-differently-expressed genes were manually curated to portray the layered signaling regulatory networks.

Results

Overall, the remodeled signaling regulatory networks were separated into four major layers including extracellular, membrane, cytoplasm and nucleus, which led to the identification of five core biological processes and four signaling pathways associated with colorectal carcinogenesis. As a result, our biological interpretation highlighted the importance of EGF/EGFR signaling pathway, EPO signaling pathway, T cell signal transduction and members of the BCR signaling pathway, which were responsible for the malignant transition of CRC from the benign UC to the aggressive one.

Conclusions

The present study illustrated a standardized normalization approach for cross-study microarray expression data sets. Our model for signaling networks construction was based on the experimentally-supported interaction and microarray co-expression modeling. Pathway-based signaling regulatory networks analysis sketched a directive insight into colorectal carcinogenesis, which was of significant importance to monitor disease progression and improve therapeutic interventions.