Complexity of the microRNA repertoire revealed by next-generation sequencing

MicroRNAs (miRNAs) have been implicated to play key roles in normal physiological functions, and altered expression of specific miRNAs has been associated with a number of diseases. It is of great interest to understand their roles and a prerequisite for such study is the ability to comprehensively and accurately assess the levels of the entire repertoire of miRNAs in a given sample. It has been shown that some miRNAs frequently have sequence variations termed isomirs. To better understand the extent of miRNA sequence heterogeneity and its potential implications for miRNA function and measurement, we conducted a comprehensive survey of miRNA sequence variations from human and mouse samples using next generation sequencing platforms. Our results suggest that the process of generating this isomir spectrum might not be random and that heterogeneity at the ends of miRNA affects the consistency and accuracy of miRNA level measurement. In addition, we have constructed a database from our sequencing data that catalogs the entire repertoire of miRNA sequences (http://galas.systemsbiology.net/cgi-bin/isomir/find.pl). This enables users to determine the most abundant sequence and the degree of heterogeneity for each individual miRNA species. This information will be useful both to better understand the functions of isomirs and to improve probe or primer design for miRNA detection and measurement.     

miTALOS v2: Analyzing Tissue Specific microRNA Function

MicroRNAs are involved in almost all biological processes and have emerged as regulators of signaling pathways. We show that miRNA target genes and pathway genes are not uniformly expressed across human tissues. To capture tissue specific effects, we developed a novel methodology for tissue specific pathway analysis of miRNAs. We incorporated the most recent and highest quality miRNA targeting data (TargetScan and StarBase), RNA-seq based gene expression data (EBI Expression Atlas) and multiple new pathway data sources to increase the biological relevance of the predicted miRNA-pathway associations. We identified new potential roles of miR-199a-3p, miR-199b-3p and the miR-200 family in hepatocellular carcinoma, involving the regulation of metastasis through MAPK and Wnt signaling. Also, an association of miR-571 and Notch signaling in liver fibrosis was proposed. To facilitate data update and future extensions of our tool, we developed a flexible database backend using the graph database neo4j. The new backend as well as the novel methodology were included in the updated miTALOS v2, a tool that provides insights into tissue specific miRNA regulation of biological pathways. miTALOS v2 is available at http://mips.helmholtz-muenchen.de/mitalos.     

WebPropagate: A Web Server for Network Propagation

Network propagation is a powerful tool for genetic analysis which is widely used to identify genes and genetic modules that underlie a process of interest. Here we provide a graphical, web-based platform (http://anat.cs.tau.ac.il/WebPropagate/) in which researchers can easily apply variants of this method to data sets of interest using up-to-date networks of protein–protein interactions in several organisms.     

RNAthermsw: Direct Temperature Simulations for Predicting the Location of RNA Thermometers

The mechanism of RNA thermometers is a subject of growing interest. Also known as RNA thermosensors, these temperature-sensitive segments of the mRNA regulate gene expression by changing their secondary structure in response to temperature fluctuations. The detection of RNA thermometers in various genes of interest is valuable as it could lead to the discovery of new thermometers participating in fundamental processes such as preferential translation during heat-shock. RNAthermsw is a user-friendly webserver for predicting the location of RNA thermometers using direct temperature simulations. It operates by analyzing dotted figures generated as a result of a moving window that performs successive energy minimization folding predictions. Inputs include the RNA sequence, window size, and desired temperature change. RNAthermsw can be freely accessed at http://www.cs.bgu.ac.il/~rnathemsw/RNAthemsw/ (with the slash sign at the end). The website contains a help page with explanations regarding the exact usage.     

MicroRNA-218 Is Deleted and Downregulated in Lung Squamous Cell Carcinoma

MicroRNAs (miRNAs) are a family of small, non-coding RNA species functioning as negative regulators of multiple target genes including tumour suppressor genes and oncogenes. Many miRNA gene loci are located within cancer-associated genomic regions. To identify potential new amplified oncogenic and/or deleted tumour suppressing miRNAs in lung cancer, we inferred miRNA gene dosage from high dimensional arrayCGH data. From miRBase v9.0 (http://microrna.sanger.ac.uk), 474 human miRNA genes were physically mapped to regions of chromosomal loss or gain identified from a high-resolution genome-wide arrayCGH study of 132 primary non-small cell lung cancers (NSCLCs) (a training set of 60 squamous cell carcinomas and 72 adenocarcinomas). MiRNAs were selected as candidates if their immediately flanking probes or host gene were deleted or amplified in at least 25% of primary tumours using both Analysis of Copy Errors algorithm and fold change (≥±1.2) analyses. Using these criteria, 97 miRNAs mapped to regions of aberrant copy number. Analysis of three independent published lung cancer arrayCGH datasets confirmed that 22 of these miRNA loci showed directionally concordant copy number variation. MiR-218, encoded on 4p15.31 and 5q35.1 within two host genes (SLIT2 and SLIT3), in a region of copy number loss, was selected as a priority candidate for follow-up as it is reported as underexpressed in lung cancer. We confirmed decreased expression of mature miR-218 and its host genes by qRT-PCR in 39 NSCLCs relative to normal lung tissue. This downregulation of miR-218 was found to be associated with a history of cigarette smoking, but not human papilloma virus. Thus, we show for the first time that putative lung cancer-associated miRNAs can be identified from genome-wide arrayCGH datasets using a bioinformatics mapping approach, and report that miR-218 is a strong candidate tumour suppressing miRNA potentially involved in lung cancer.     

Identifying Cancer Specific Functionally Relevant miRNAs from Gene Expression and miRNA-to-Gene Networks Using Regularized Regression

Identifying microRNA signatures for the different types and subtypes of cancer can result in improved detection, characterization and understanding of cancer and move us towards more personalized treatment strategies. However, using microRNA''s differential expression (tumour versus normal) to determine these signatures may lead to inaccurate predictions and low interpretability because of the noisy nature of miRNA expression data. We present a method for the selection of biologically active microRNAs using gene expression data and microRNA-to-gene interaction network. Our method is based on a linear regression with an elastic net regularization. Our simulations show that, with our method, the active miRNAs can be detected with high accuracy and our approach is robust to high levels of noise and missing information. Furthermore, our results on real datasets for glioblastoma and prostate cancer are confirmed by microRNA expression measurements. Our method leads to the selection of potentially functionally important microRNAs. The associations of some of our identified miRNAs with cancer mechanisms are already confirmed in other studies (hypoxia related hsa-mir-210 and apoptosis-related hsa-mir-296-5p). We have also identified additional miRNAs that were not previously studied in the context of cancer but are coherently predicted as active by our method and may warrant further investigation. The code is available in Matlab and R and can be downloaded on http://www.cs.toronto.edu/goldenberg/Anna_Goldenberg/Current_Research.html.     

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/.     

Regulation of Pattern Formation and Gene Amplification During Drosophila Oogenesis by the miR-318 microRNA

Pattern formation during epithelial development requires the coordination of multiple signaling pathways. Here, we investigate the functions of an ovary-enriched miRNA, miR-318, in epithelial development during Drosophila oogenesis. mir-318 maternal loss-of-function mutants were female-sterile and laid eggs with abnormal morphology. Removal of mir-318 disrupted the dorsal–anterior follicle cell patterning, resulting in abnormal dorsal appendages. mir-318 mutant females also produced thin and fragile eggshells due to impaired chorion gene amplification. We provide evidence that the ecdysone signaling pathway activates expression of miR-318 and that miR-318 cooperates with Tramtrack69 to control the switch from endocycling to chorion gene amplification during differentiation of the follicular epithelium. The multiple functions of miR-318 in oogenesis illustrate the importance of miRNAs in maintaining cell fate and in promoting the developmental transition in the female follicular epithelium.     

Hedgehog Signaling Strength Is Orchestrated by the mir-310 Cluster of MicroRNAs in Response to Diet

Since the discovery of microRNAs (miRNAs) only two decades ago, they have emerged as an essential component of the gene regulatory machinery. miRNAs have seemingly paradoxical features: a single miRNA is able to simultaneously target hundreds of genes, while its presence is mostly dispensable for animal viability under normal conditions. It is known that miRNAs act as stress response factors; however, it remains challenging to determine their relevant targets and the conditions under which they function. To address this challenge, we propose a new workflow for miRNA function analysis, by which we found that the evolutionarily young miRNA family, the mir-310s (mir-310/mir-311/mir-312/mir-313), are important regulators of Drosophila metabolic status. mir-310s-deficient animals have an abnormal diet-dependent expression profile for numerous diet-sensitive components, accumulate fats, and show various physiological defects. We found that the mir-310s simultaneously repress the production of several regulatory factors (Rab23, DHR96, and Ttk) of the evolutionarily conserved Hedgehog (Hh) pathway to sharpen dietary response. As the mir-310s expression is highly dynamic and nutrition sensitive, this signal relay model helps to explain the molecular mechanism governing quick and robust Hh signaling responses to nutritional changes. Additionally, we discovered a new component of the Hh signaling pathway in Drosophila, Rab23, which cell autonomously regulates Hh ligand trafficking in the germline stem cell niche. How organisms adjust to dietary fluctuations to sustain healthy homeostasis is an intriguing research topic. These data are the first to report that miRNAs can act as executives that transduce nutritional signals to an essential signaling pathway. This suggests miRNAs as plausible therapeutic agents that can be used in combination with low calorie and cholesterol diets to manage quick and precise tissue-specific responses to nutritional changes.     

TargetCompare: A web interface to compare simultaneous miRNAs targets

MicroRNAs (miRNAs) are small non-coding nucleotide sequences between 17 and 25 nucleotides in length that primarily function in the regulation of gene expression. A since miRNA has thousand of predict targets in a complex, regulatory cell signaling network. Therefore, it is of interest to study multiple target genes simultaneously. Hence, we describe a web tool (developed using Java programming language and MySQL database server) to analyse multiple targets of pre-selected miRNAs. We cross validated the tool in eight most highly expressed miRNAs in the antrum region of stomach. This helped to identify 43 potential genes that are target of at least six of the referred miRNAs. The developed tool aims to reduce the randomness and increase the chance of selecting strong candidate target genes and miRNAs responsible for playing important roles in the studied tissue.

Availability

http://lghm.ufpa.br/targetcompare     

PROmiRNA: a new miRNA promoter recognition method uncovers the complex regulation of intronic miRNAs

The regulation of intragenic miRNAs by their own intronic promoters is one of the open problems of miRNA biogenesis. Here, we describe PROmiRNA, a new approach for miRNA promoter annotation based on a semi-supervised statistical model trained on deepCAGE data and sequence features. We validate our results with existing annotation, PolII occupancy data and read coverage from RNA-seq data. Compared to previous methods PROmiRNA increases the detection rate of intronic promoters by 30%, allowing us to perform a large-scale analysis of their genomic features, as well as elucidate their contribution to tissue-specific regulation. PROmiRNA can be downloaded from http://promirna.molgen.mpg.de.