miRSystem: An Integrated System for Characterizing Enriched Functions and Pathways of MicroRNA Targets

Background

Many prediction tools for microRNA (miRNA) targets have been developed, but inconsistent predictions were observed across multiple algorithms, which can make further analysis difficult. Moreover, the nomenclature of human miRNAs changes rapidly. To address these issues, we developed a web-based system, miRSystem, for converting queried miRNAs to the latest annotation and predicting the function of miRNA by integrating miRNA target gene prediction and function/pathway analyses.

Results

First, queried miRNA IDs were converted to the latest annotated version to prevent potential conflicts resulting from multiple aliases. Next, by combining seven algorithms and two validated databases, potential gene targets of miRNAs and their functions were predicted based on the consistency across independent algorithms and observed/expected ratios. Lastly, five pathway databases were included to characterize the enriched pathways of target genes through bootstrap approaches. Based on the enriched pathways of target genes, the functions of queried miRNAs could be predicted.

Conclusions

MiRSystem is a user-friendly tool for predicting the target genes and their associated pathways for many miRNAs simultaneously. The web server and the documentation are freely available at http://mirsystem.cgm.ntu.edu.tw/.     

A comparison of performance of plant miRNA target prediction tools and the characterization of features for genome-wide target prediction

Background

Deep-sequencing has enabled the identification of large numbers of miRNAs and siRNAs, making the high-throughput target identification a main limiting factor in defining their function. In plants, several tools have been developed to predict targets, majority of them being trained on Arabidopsis datasets. An extensive and systematic evaluation has not been made for their suitability for predicting targets in species other than Arabidopsis. Nor, these have not been evaluated for their suitability for high-throughput target prediction at genome level.

Results

We evaluated the performance of 11 computational tools in identifying genome-wide targets in Arabidopsis and other plants with procedures that optimized score-cutoffs for estimating targets. Targetfinder was most efficient [89% ‘precision’ (accuracy of prediction), 97% ‘recall’ (sensitivity)] in predicting ‘true-positive’ targets in Arabidopsis miRNA-mRNA interactions. In contrast, only 46% of true positive interactions from non-Arabidopsis species were detected, indicating low ‘recall’ values. Score optimizations increased the ‘recall’ to only 70% (corresponding ‘precision’: 65%) for datasets of true miRNA-mRNA interactions in species other than Arabidopsis. Combining the results of Targetfinder and psRNATarget delivers high true positive coverage, whereas the intersection of psRNATarget and Tapirhybrid outputs deliver highly ‘precise’ predictions. The large number of ‘false negative’ predictions delivered from non-Arabidopsis datasets by all the available tools indicate the diversity in miRNAs-mRNA interaction features between Arabidopsis and other species. A subset of miRNA-mRNA interactions differed significantly for features in seed regions as well as the total number of matches/mismatches.

Conclusion

Although, many plant miRNA target prediction tools may be optimized to predict targets with high specificity in Arabidopsis, such optimized thresholds may not be suitable for many targets in non-Arabidopsis species. More importantly, non-conventional features of miRNA-mRNA interaction may exist in plants indicating alternate mode of miRNA target recognition. Incorporation of these divergent features would enable next-generation of algorithms to better identify target interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-348) 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.     

MicroRNA expression profiling of the fifth-instar posterior silk gland of Bombyx mori

Background

The growth and development of the posterior silk gland and the biosynthesis of the silk core protein at the fifth larval instar stage of Bombyx mori are of paramount importance for silk production.

Results

Here, aided by next-generation sequencing and microarry assay, we profile 1,229 microRNAs (miRNAs), including 728 novel miRNAs and 110 miRNA/miRNA* duplexes, of the posterior silk gland at the fifth larval instar. Target gene prediction yields 14,222 unique target genes from 1,195 miRNAs. Functional categorization classifies the targets into complex pathways that include both cellular and metabolic processes, especially protein synthesis and processing.

Conclusion

The enrichment of target genes in the ribosome-related pathway indicates that miRNAs may directly regulate translation. Our findings pave a way for further functional elucidation of these miRNAs and their targets in silk production.

Electronic supplementary material

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

Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach

Background

While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia.

Methodology/Principal Findings

A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue.

Conclusions/Significance

This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.     

Downregulation of miR-31, miR-155, and miR-564 in chronic myeloid leukemia cells

Background/Aims

MicroRNAs (miRNAs) are short non-coding regulatory RNAs that control gene expression and play an important role in cancer development and progression. However, little is known about the role of miRNAs in chronic myeloid leukemia (CML). Our objective is to decipher a miRNA expression signature associated with CML and to determine potential target genes and signaling pathways affected by these signature miRNAs.

Results

Using miRNA microarrays and miRNA real-time PCR we characterized the miRNAs expression profile of CML cell lines and patients in reference to non-CML cell lines and healthy blood. Of all miRNAs tested, miR-31, miR-155, and miR-564 were down-regulated in CML cells. Down-regulation of these miRNAs was dependent on BCR-ABL activity. We next analyzed predicted targets and affected pathways of the deregulated miRNAs. As expected, in K562 cells, the expression of several of these targets was inverted to that of the miRNA putatively regulating them. Reassuringly, the analysis identified CML as the main disease associated with these miRNAs. MAPK, ErbB, mammalian target of rapamycin (mTOR) and vascular endothelial growth factor (VEGF) were the main molecular pathways related with these expression patterns. Utilizing Venn diagrams we found appreciable overlap between the CML-related miRNAs and the signaling pathways-related miRNAs.

Conclusions

The miRNAs identified in this study might offer a pivotal role in CML. Nevertheless, while these data point to a central disease, the precise molecular pathway/s targeted by these miRNAs is variable implying a high level of complexity of miRNA target selection and regulation. These deregulated miRNAs highlight new candidate gene targets allowing for a better understanding of the molecular mechanism underlying the development of CML, and propose possible new avenues for therapeutic treatment.     

The Fate of miRNA* Strand through Evolutionary Analysis: Implication for Degradation As Merely Carrier Strand or Potential Regulatory Molecule?

Background

During typical microRNA (miRNA) biogenesis, one strand of a ∼22 nt RNA duplex is preferentially selected for entry into a silencing complex, whereas the other strand, known as the passenger strand or miRNA* strand, is degraded. Recently, some miRNA* sequences were reported as guide miRNAs with abundant expression. Here, we intended to discover evolutionary implication of the fate of miRNA* strand by analyzing miRNA/miRNA* sequences across vertebrates.

Principal Findings

Mature miRNAs based on gene families were well conserved especially for their seed sequences across vertebrates, while their passenger strands always showed various divergence patterns. The divergence mainly resulted from divergence of different animal species, homologous miRNA genes and multicopy miRNA hairpin precursors. Some miRNA* sequences were phylogenetically conserved in seed and anchor sequences similar to mature miRNAs, while others revealed high levels of nucleotide divergence despite some of their partners were highly conserved. Most of those miRNA precursors that could generate abundant miRNAs from both strands always were well conserved in sequences of miR-#-5p and miR-#-3p, especially for their seed sequences.

Conclusions

The final fate of miRNA* strand, either degraded as merely carrier strand or expressed abundantly as potential functional guide miRNA, may be destined across evolution. Well-conserved miRNA* strands, particularly conservation in seed sequences, maybe afford potential opportunities for contributing to regulation network. The study will broaden our understanding of potential functional miRNA* species.     

Diversity,expression and mRNA targeting abilities of Argonaute-targeting miRNAs among selected vascular plants

Background

Micro (mi)RNAs are important regulators of plant development. Across plant lineages, Dicer-like 1 (DCL1) proteins process long ds-like structures to produce micro (mi) RNA duplexes in a stepwise manner. These miRNAs are incorporated into Argonaute (AGO) proteins and influence expression of RNAs that have sequence complementarity with miRNAs. Expression levels of AGOs are greatly regulated by plants in order to minimize unwarranted perturbations using miRNAs to target mRNAs coding for AGOs. AGOs may also have high promoter specificity-sometimes expression of AGO can be limited to just a few cells in a plant. Viral pathogens utilize various means to counter antiviral roles of AGOs including hijacking the host encoded miRNAs to target AGOs. Two host encoded miRNAs namely miR168 and miR403 that target AGOs have been described in the model plant Arabidopsis and such a mechanism is thought to be well conserved across plants because AGO sequences are well conserved.

Results

We show that the interaction between AGO mRNAs and miRNAs is species-specific due to the diversity in sequences of two miRNAs that target AGOs, sequence diversity among corresponding target regions in AGO mRNAs and variable expression levels of these miRNAs among vascular plants. We used miRNA sequences from 68 plant species representing 31 plant families for this analysis. Sequences of miR168 and miR403 are not conserved among plant lineages, but surprisingly they differ drastically in their sequence diversity and expression levels even among closely related plants. Variation in miR168 expression among plants correlates well with secondary structures/length of loop sequences of their precursors.

Conclusions

Our data indicates a complex AGO targeting interaction among plant lineages due to miRNA sequence diversity and sequences of miRNA targeting regions among AGO mRNAs, thus leading to the assumption that the perturbations by viruses that use host miRNAs to target antiviral AGOs can only be species-specific. We also show that rapid evolution and likely loss of expression of miR168 isoforms in tobacco is related to the insertion of MITE-like transposons between miRNA and miRNA* sequences, a possible mechanism showing how miRNAs are lost in few plant lineages even though other close relatives have abundantly expressing miRNAs.

Electronic supplementary material

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

miR-29b,miR-205 and miR-221 Enhance Chemosensitivity to Gemcitabine in HuH28 Human Cholangiocarcinoma Cells

Background and Aims

Cholangiocarcinoma (CCA) is highly resistant to chemotherapy, including gemcitabine (Gem) treatment. MicroRNAs (miRNAs) are endogenous, non-coding, short RNAs that can regulate multiple genes expression. Some miRNAs play important roles in the chemosensitivity of tumors. Here, we examined the relationship between miRNA expression and the sensitivity of CCA cells to Gem.

Methods

Microarray analysis was used to determine the miRNA expression profiles of two CCA cell lines, HuH28 and HuCCT1. To determine the effect of candidate miRNAs on Gem sensitivity, expression of each candidate miRNA was modified via either transfection of a miRNA mimic or transfection of an anti-oligonucleotide. Ontology-based programs were used to identify potential target genes of candidate miRNAs that were confirmed to affect the Gem sensitivity of CCA cells.

Results

HuCCT1 cells were more sensitive to Gem than were HuH28 cells, and 18 miRNAs were differentially expressed whose ratios over ± 2log2 between HuH28 and HuCCT1. Among these 18 miRNAs, ectopic overexpression of each of three downregulated miRNAs in HuH28 (miR-29b, miR-205, miR-221) restored Gem sensitivity to HuH28. Suppression of one upregulated miRNA in HuH28, miR-125a-5p, inhibited HuH28 cell proliferation independently to Gem treatment. Selective siRNA-mediated downregulation of either of two software-predicted targets, PIK3R1 (target of miR-29b and miR-221) or MMP-2 (target of miR-29b), also conferred Gem sensitivity to HuH28.

Conclusions

miRNA expression profiling was used to identify key miRNAs that regulate Gem sensitivity in CCA cells, and software that predicts miRNA targets was used to identify promising target genes for anti-tumor therapies.