Integrated analysis of microRNA-target interactions with clinical outcomes for cancers

Background

Clinical statement alone is not enough to predict the progression of disease. Instead, the gene expression profiles have been widely used to forecast clinical outcomes. Many genes related to survival have been identified, and recently miRNA expression signatures predicting patient survival have been also investigated for several cancers. However, miRNAs and their target genes associated with clinical outcomes have remained largely unexplored.

Methods

Here, we demonstrate a survival analysis based on the regulatory relationships of miRNAs and their target genes. The patient survivals for the two major cancers, ovarian cancer and glioblastoma multiforme (GBM), are investigated through the integrated analysis of miRNA-mRNA interaction pairs.

Results

We found that there is a larger survival difference between two patient groups with an inversely correlated expression profile of miRNA and mRNA. It supports the idea that signatures of miRNAs and their targets related to cancer progression can be detected via this approach.

Conclusions

This integrated analysis can help to discover coordinated expression signatures of miRNAs and their target mRNAs that can be employed for therapeutics in human cancers.

     

Comprehensive annotation of microRNA expression profiles

Background

MicroRNAs (miRNAs) regulate many biological processes by post-translational gene silencing. Analysis of miRNA expression profiles is a reliable method for investigating particular biological processes due to the stability of miRNA and the development of advanced sequencing methods. However, this approach is limited by the broad specificity of miRNAs, which may target several mRNAs.

Result

In this study, we developed a method for comprehensive annotation of miRNA array or deep sequencing data for investigation of cellular biological effects. Using this method, the specific pathways and biological processes involved in Alzheimer’s disease were predicted with high correlation in four independent samples. Furthermore, this method was validated for evaluation of cadmium telluride (CdTe) nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways associated with CdTe nanoparticle exposure, which is consistent with previous studies.

Conclusions

Our findings contribute to the validation of miRNA microarray or deep sequencing results for early diagnosis of disease and evaluation of the biological safety of new materials and drugs.

     

SAGE detects microRNA precursors

Background

MicroRNAs (miRNAs) have been shown to play important roles in regulating gene expression. Since miRNAs are often evolutionarily conserved and their precursors can be folded into stem-loop hairpins, many miRNAs have been predicted. Yet experimental confirmation is difficult since miRNA expression is often specific to particular tissues and developmental stages.

Results

Analysis of 29 human and 230 mouse longSAGE libraries revealed the expression of 22 known and 10 predicted mammalian miRNAs. Most were detected in embryonic tissues. Four SAGE tags detected in human embryonic stem cells specifically match a cluster of four human miRNAs (mir-302a, b, c&d) known to be expressed in embryonic stem cells. LongSAGE data also suggest the existence of a mouse homolog of human and rat mir-493.

Conclusion

The observation that some orphan longSAGE tags uniquely match miRNA precursors provides information about the expression of some known and predicted miRNAs.

     

Upregulation of microRNA 344a-3p is involved in curcumin induced apoptosis in RT4 schwannoma cells

Background

Schwannoma arising from peripheral nervous sheaths is a benign tumor.

Methods

To evaluate cell cytotoxicity, (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) tetrazolium reduction and terminal deoxynucleotidyltransferase UTP nick-end labeling (TUNEL) assays were used. A microRNA (miRNA) array was used to identify the miRNAs involved in curcumin-induced apoptosis. To examine miRNA expression, quantitative RT-PCR was used.

Results

In this study, curcumin exerted cellular cytotoxicity against RT4 schwannoma cells, with an increase in TUNEL-positive cells. Curcumin also activated the expression of apoptotic proteins, such as polyADP ribose polymerase, caspase-3, and caspase-9. The miRNA array revealed that seven miRNAs (miRNA 350, miRNA 17-2-3p, let 7e-3p, miRNA1224, miRNA 466b-1-3p, miRNA 18a-5p, and miRNA 322-5p) were downregulated following treatment with both 10 and 20 μM curcumin in RT4 cells, while four miRNAs (miRNA122-5p, miRNA 3473, miRNA182, and miRNA344a-3p) were upregulated. Interestingly, transfection with a miRNA 344a-3p mimic downregulated the mRNA expression of Bcl2 and upregulated that of Bax, Curcumin treatment in RT 4 cells also reduced the mRNA expression of Bcl2 and enhanced expression of Bax, Overexpression of miRNA344a-3p mimic combined with curcumin treatment activated the expression of apoptotic proteins, including procaspase-9 and cleaved caspase-3 while inhibition of miRNA 344a-3p using miR344a-3p inhibitor repressed cleaved caspase-3 and -9 in curcumin treated RT-4 cells compared to control.

Conclusions

Our findings demonstrate that curcumin induces apoptosis in schwannoma cells via miRNA 344a-3p. Thus, curcumin may serve as a potent therapeutic agent for the treatment of schwannoma.

     

A novel two-layer SVM model in miRNA Drosha processing site detection

Background

MicroRNAs (miRNAs) are a large class of non-coding RNAs with important functions wide spread in animals, plants and viruses. Studies showed that an RNase III family member called Drosha recognizes most miRNAs, initiates their processing and determines the mature miRNAs. The Drosha processing sites identification will shed some light on both miRNA identification and understanding the mechanism of Drosha processing.

Methods

We developed a computational method for Drosha processing site predicting, named as DroshaPSP, which employs a two-layer mathematical model to integrate structure feature in the first layer and sequence features in the second layer. The performance of DroshaPSP was estimated by 5-fold cross-validation and measured by ACC (accuracy), Sn (sensitivity), Sp (specificity), P (precision) and MCC (Matthews correlation coefficient).

Results

The results of testing DroshaPSP on the miRNA data of Drosophila melanogaster indicated that the Sn, Sp, and MCC thereof reach to 0.86, 0.99 and 0.86 respectively.

Conclusions

We found the Shannon entropy, a chemical kinetics feature, is a significant feature in telling the true sites among the nearby sites and improving the performance.

     

Novel link prediction for large-scale miRNA-lncRNA interaction network in a bipartite graph

Background

Current knowledge and data on miRNA-lncRNA interactions is still limited and little effort has been made to predict target lncRNAs of miRNAs. Accumulating evidences suggest that the interaction patterns between lncRNAs and miRNAs are closely related to relative expression level, forming a titration mechanism. It could provide an effective approach for characteristic feature extraction. In addition, using the coding non-coding co-expression network and sequence data could also help to measure the similarities among miRNAs and lncRNAs. By mathematically analyzing these types of similarities, we come up with two findings that (i) lncRNAs/miRNAs tend to collaboratively interact with miRNAs/lncRNAs of similar expression profiles, and vice versa, and (ii) those miRNAs interacting with a cluster of common target genes tend to jointly target at the common lncRNAs.

Methods

In this work, we developed a novel group preference Bayesian collaborative filtering model called GBCF for picking up a top-k probability ranking list for an individual miRNA or lncRNA based on the known miRNA-lncRNA interaction network.

Results

To evaluate the effectiveness of GBCF, leave-one-out and k-fold cross validations as well as a series of comparison experiments were carried out. GBCF achieved the values of area under ROC curve of 0.9193, 0.8354+/??0.0079, 0.8615+/??0.0078, and 0.8928+/??0.0082 based on leave-one-out, 2-fold, 5-fold, and 10-fold cross validations respectively, demonstrating its reliability and robustness.

Conclusions

GBCF could be used to select potential lncRNA targets of specific miRNAs and offer great insights for further researches on ceRNA regulation network.

     

Profiling microRNAs in individuals at risk of progression to rheumatoid arthritis

Background

Individuals at risk of rheumatoid arthritis (RA) demonstrate systemic autoimmunity in the form of anti-citrullinated peptide antibodies (ACPA). MicroRNAs (miRNAs) are implicated in established RA. This study aimed to (1) compare miRNA expression between healthy individuals and those at risk of and those that develop RA, (2) evaluate the change in expression of miRNA from “at-risk” to early RA and (3) explore whether these miRNAs could inform a signature predictive of progression from “at-risk” to RA.

Methods

We performed global profiling of 754 miRNAs per patient on a matched serum sample cohort of 12 anti-cyclic citrullinated peptide (CCP)?+?“at-risk” individuals that progressed to RA. Each individual had a serum sample from baseline and at time of detection of synovitis, forming the matched element. Healthy controls were also studied. miRNAs with a fold difference/fold change of four in expression level met our primary criterion for selection as candidate miRNAs. Validation of the miRNAs of interest was conducted using custom miRNA array cards on matched samples (baseline and follow up) in 24 CCP+ individuals; 12 RA progressors and 12 RA non-progressors.

Results

We report on the first study to use matched serum samples and a comprehensive miRNA array approach to identify in particular, three miRNAs (miR-22, miR-486-3p, and miR-382) associated with progression from systemic autoimmunity to RA inflammation. MiR-22 demonstrated significant fold difference between progressors and non-progressors indicating a potential biomarker role for at-risk individuals.

Conclusions

This first study using a cohort with matched serum samples provides important mechanistic insights in the transition from systemic autoimmunity to inflammatory disease for future investigation, and with further evaluation, might also serve as a predictive biomarker.

     

Estimating survival time of patients with glioblastoma multiforme and characterization of the identified microRNA signatures

Background

Though glioblastoma multiforme (GBM) is the most frequently occurring brain malignancy in adults, clinical treatment still faces challenges due to poor prognoses and tumor relapses. Recently, microRNAs (miRNAs) have been extensively used with the aim of developing accurate molecular therapies, because of their emerging role in the regulation of cancer-related genes. This work aims to identify the miRNA signatures related to survival of GBM patients for developing molecular therapies.

Results

This work proposes a support vector regression (SVR)-based estimator, called SVR-GBM, to estimate the survival time in patients with GBM using their miRNA expression profiles. SVR-GBM identified 24 out of 470 miRNAs that were significantly associated with survival of GBM patients. SVR-GBM had a mean absolute error of 0.63 years and a correlation coefficient of 0.76 between the real and predicted survival time. The 10 top-ranked miRNAs according to prediction contribution are as follows: hsa-miR-222, hsa-miR-345, hsa-miR-587, hsa-miR-526a, hsa-miR-335, hsa-miR-122, hsa-miR-24, hsa-miR-433, hsa-miR-574 and hsa-miR-320. Biological analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway on the identified miRNAs revealed their influence in GBM cancer.

Conclusion

The proposed SVR-GBM using an optimal feature selection algorithm and an optimized SVR to identify the 24 miRNA signatures associated with survival of GBM patients. These miRNA signatures are helpful to uncover the individual role of miRNAs in GBM prognosis and develop miRNA-based therapies.

     

Identification of glioblastoma gene prognosis modules based on weighted gene co-expression network analysis

Background

Glioblastoma multiforme, the most prevalent and aggressive brain tumour, has a poor prognosis. The molecular mechanisms underlying gliomagenesis remain poorly understood. Therefore, molecular research, including various markers, is necessary to understand the occurrence and development of glioma.

Method

Weighted gene co-expression network analysis (WGCNA) was performed to construct a gene co-expression network in TCGA glioblastoma samples. Gene ontology (GO) and pathway-enrichment analysis were used to identify significance of gene modules. Cox proportional hazards regression model was used to predict outcome of glioblastoma patients.

Results

We performed weighted gene co-expression network analysis (WGCNA) and identified a gene module (yellow module) related to the survival time of TCGA glioblastoma samples. Then, 228 hub genes were calculated based on gene significance (GS) and module significance (MS). Four genes (OSMR + SOX21?+?MED10?+?PTPRN) were selected to construct a Cox proportional hazards regression model with high accuracy (AUC?=?0.905). The prognostic value of the Cox proportional hazards regression model was also confirmed in GSE16011 dataset (GBM: n?=?156).

Conclusion

We developed a promising mRNA signature for estimating overall survival in glioblastoma patients.

     

Bioavailability of transgenic microRNAs in genetically modified plants

Background

Transgenic expression of small RNAs is a prevalent approach in agrobiotechnology for the global enhancement of plant foods. Meanwhile, emerging studies have, on the one hand, emphasized the potential of transgenic microRNAs (miRNAs) as novel dietary therapeutics and, on the other, suggested potential food safety issues if harmful miRNAs are absorbed and bioactive. For these reasons, it is necessary to evaluate the bioavailability of transgenic miRNAs in genetically modified crops.

Results

As a pilot study, two transgenic Arabidopsis lines ectopically expressing unique miRNAs were compared and contrasted with the plant bioavailable small RNA MIR2911 for digestive stability and serum bioavailability. The expression levels of these transgenic miRNAs in Arabidopsis were found to be comparable to that of MIR2911 in fresh tissues. Assays of digestive stability in vitro and in vivo suggested the transgenic miRNAs and MIR2911 had comparable resistance to degradation. Healthy mice consuming diets rich in Arabidopsis lines expressing these miRNAs displayed MIR2911 in the bloodstream but no detectable levels of the transgenic miRNAs.

Conclusions

These preliminary results imply digestive stability and high expression levels of miRNAs in plants do not readily equate to bioavailability. This initial work suggests novel engineering strategies be employed to enhance miRNA bioavailability when attempting to use transgenic foods as a delivery platform.

     

Grouping miRNAs of similar functions via weighted information content of gene ontology

Background

Regulation mechanisms between miRNAs and genes are complicated. To accomplish a biological function, a miRNA may regulate multiple target genes, and similarly a target gene may be regulated by multiple miRNAs. Wet-lab knowledge of co-regulating miRNAs is limited. This work introduces a computational method to group miRNAs of similar functions to identify co-regulating miRNAsfrom a similarity matrix of miRNAs.

Results

We define a novel information content of gene ontology (GO) to measure similarity between two sets of GO graphs corresponding to the two sets of target genes of two miRNAs. This between-graph similarity is then transferred as a functional similarity between the two miRNAs. Our definition of the information content is based on the size of a GO term’s descendants, but adjusted by a weight derived from its depth level and the GO relationships at its path to the root node or to the most informative common ancestor (MICA). Further, a self-tuning technique and the eigenvalues of the normalized Laplacian matrix are applied to determine the optimal parameters for the spectral clustering of the similarity matrix of the miRNAs.

Conclusions

Experimental results demonstrate that our method has better clustering performance than the existing edge-based, node-based or hybrid methods. Our method has also demonstrated a novel usefulness for the function annotation of new miRNAs, as reported in the detailed case studies.

     

FMSM: a novel computational model for predicting potential miRNA biomarkers for various human diseases

Background

MicroRNA (miRNA) plays a key role in regulation mechanism of human biological processes, including the development of disease and disorder. It is necessary to identify potential miRNA biomarkers for various human diseases. Computational prediction model is expected to accelerate the process of identification.

Results

Considering the limitations of previously proposed models, we present a novel computational model called FMSM. It infers latent miRNA biomarkers involved in the mechanism of various diseases based on the known miRNA-disease association network, miRNA expression similarity, disease semantic similarity and Gaussian interaction profile kernel similarity. FMSM achieves reliable prediction performance in 5-fold and leave-one-out cross validations with area under ROC curve (AUC) values of 0.9629+/??0.0127 and 0.9433, respectively, which outperforms the state-of-the-art competitors and classical algorithms. In addition, 19 of top 25 predicted miRNAs have been validated to have associations with Colonic Neoplasms in case study.

Conclusions

A factored miRNA similarity based model and miRNA expression similarity substantially contribute to the well-performing prediction. The list of the predicted most latent miRNA biomarkers of various human diseases is publicized. It is anticipated that FMSM could serve as a useful tool guiding the future experimental validation for those promising miRNA biomarker candidates.

     

Epithelial-specific histone modification of the <Emphasis Type="Italic">miR-96/182</Emphasis> locus targeting <Emphasis Type="Italic">AMAP1</Emphasis> mRNA predisposes p53 to suppress cell invasion in epithelial cells

Background

TP53 mutations in cancer cells often evoke cell invasiveness, whereas fibroblasts show invasiveness in the presence of intact TP53. AMAP1 (also called DDEF1 or ASAP1) is a downstream effector of ARF6 and is essential for the ARF6-driven cell-invasive phenotype. We found that AMAP1 levels are under the control of p53 (TP53 gene product) in epithelial cells but not in fibroblasts, and here addressed that molecular basis of the epithelial-specific function of p53 in suppressing invasiveness via targeting AMAP1.

Methods

Using MDA-MB-231 cells expressing wild-type and p53 mutants, we identified miRNAs in which their expression is controlled by normal-p53. Among them, we identified miRNAs that target AMAP1 mRNA, and analyzed their expression levels and epigenetic statuses in epithelial cells and nonepithelial cells.

Results

We found that normal-p53 suppresses AMAP1 mRNA in cancer cells and normal epithelial cells, and that more than 30 miRNAs are induced by normal-p53. Among them, miR-96 and miR-182 were found to target the 3′-untranslated region of AMAP1 mRNA. Fibroblasts did not express these miRNAs at detectable levels. The ENCODE dataset demonstrated that the promoter region of the miR-183-96-182 cistron is enriched with H3K27 acetylation in epithelial cells, whereas this locus is enriched with H3K27 trimethylation in fibroblasts and other non-epithelial cells. miRNAs, such as miR-423, which are under the control of p53 but not associated with AMAP1 mRNA, demonstrated similar histone modifications at their gene loci in epithelial cells and fibroblasts, and were expressed in these cells.

Conclusion

Histone modifications of certain miRNA loci, such as the miR-183-96-182 cistron, are different between epithelial cells and non-epithelial cells. Such epithelial-specific miRNA regulation appears to provide the molecular basis for the epithelial-specific function of p53 in suppressing ARF6-driven invasiveness.

     

Microrna expression signatures predict patient progression and disease outcome in pediatric embryonal central nervous system neoplasms

Background

Although, substantial experimental evidence related to diagnosis and treatment of pediatric central nervous system (CNS) neoplasms have been demonstrated, the understanding of the etiology and pathogenesis of the disease remains scarce. Recent microRNA (miRNA)-based research reveals the involvement of miRNAs in various aspects of CNS development and proposes that they might compose key molecules underlying oncogenesis. The current study evaluated miRNA differential expression detected between pediatric embryonal brain tumors and normal controls to characterize candidate biomarkers related to diagnosis, prognosis and therapy.

Methods

Overall, 19 embryonal brain tumors; 15 Medulloblastomas (MBs) and 4 Atypical Teratoid/Rabdoid Tumors (AT/RTs) were studied. As controls, 13 samples were used; The First-Choice Human Brain Reference RNA and 12 samples from deceased children who underwent autopsy and were not present with any brain malignancy. RNA extraction was carried out using the Trizol method, whilst miRNA extraction was performed with the mirVANA miRNA isolation kit. The experimental approach included miRNA microarrays covering 1211 miRNAs. Quantitative Real-Time Polymerase Chain Reaction was performed to validate the expression profiles of miR-34a and miR-601 in all 32 samples initially screened with miRNA microarrays and in an additional independent cohort of 30 patients (21MBs and 9 AT/RTs). Moreover, meta-analyses was performed in total 27 embryonal tumor samples; 19 MBs, 8 ATRTs and 121 control samples. Twelve germinomas were also used as an independent validation cohort. All deregulated miRNAs were correlated to patients’ clinical characteristics and pathological measures.

Results

In several cases, there was a positive correlation between individual miRNA expression levels and laboratory or clinical characteristics. Based on that, miR-601 could serve as a putative tumor suppressor gene, whilst miR-34a as an oncogene. In general, miR-34a demonstrated oncogenic roles in all pediatric embryonal CNS neoplasms studied.

Conclusions

Deeper understanding of the aberrant miRNA expression in pediatric embryonal brain tumors might aid in the development of tumor-specific miRNA signatures, which could potentially afford promising biomarkers related to diagnosis, prognosis and patient targeted therapy.

     

Identification of cancer subtypes from single-cell RNA-seq data using a consensus clustering method

Background

Human cancers are complex ecosystems composed of cells with distinct molecular signatures. Such intratumoral heterogeneity poses a major challenge to cancer diagnosis and treatment. Recent advancements of single-cell techniques such as scRNA-seq have brought unprecedented insights into cellular heterogeneity. Subsequently, a challenging computational problem is to cluster high dimensional noisy datasets with substantially fewer cells than the number of genes.

Methods

In this paper, we introduced a consensus clustering framework conCluster, for cancer subtype identification from single-cell RNA-seq data. Using an ensemble strategy, conCluster fuses multiple basic partitions to consensus clusters.

Results

Applied to real cancer scRNA-seq datasets, conCluster can more accurately detect cancer subtypes than the widely used scRNA-seq clustering methods. Further, we conducted co-expression network analysis for the identified melanoma subtypes.

Conclusions

Our analysis demonstrates that these subtypes exhibit distinct gene co-expression networks and significant gene sets with different functional enrichment.