Mutual enrichment in aggregated ranked lists with applications to gene expression regulation

miRNA normalization enables joint analysis of several datasets to increase sensitivity and to reveal novel miRNAs differentially expressed in breast cancer

Different miRNA profiling protocols and technologies introduce differences in the resulting quantitative expression profiles. These include differences in the presence (and measurability) of certain miRNAs. We present and examine a method based on quantile normalization, Adjusted Quantile Normalization (AQuN), to combine miRNA expression data from multiple studies in breast cancer into a single joint dataset for integrative analysis. By pooling multiple datasets, we obtain increased statistical power, surfacing patterns that do not emerge as statistically significant when separately analyzing these datasets. To merge several datasets, as we do here, one needs to overcome both technical and batch differences between these datasets. We compare several approaches for merging and jointly analyzing miRNA datasets. We investigate the statistical confidence for known results and highlight potential new findings that resulted from the joint analysis using AQuN. In particular, we detect several miRNAs to be differentially expressed in estrogen receptor (ER) positive versus ER negative samples. In addition, we identify new potential biomarkers and therapeutic targets for both clinical groups. As a specific example, using the AQuN-derived dataset we detect hsa-miR-193b-5p to have a statistically significant over-expression in the ER positive group, a phenomenon that was not previously reported. Furthermore, as demonstrated by functional assays in breast cancer cell lines, overexpression of hsa-miR-193b-5p in breast cancer cell lines resulted in decreased cell viability in addition to inducing apoptosis. Together, these observations suggest a novel functional role for this miRNA in breast cancer. Packages implementing AQuN are provided for Python and Matlab: https://github.com/YakhiniGroup/PyAQN.

This work demonstrates a practical approach to the joint-analysis of multiple miRNA expression profiling datasets acquired with different measurement technologies. The use of different platforms in miRNA profiling can lead to major differences in results. In particular, some miRNA species are less amenable to detection and quantification by certain platforms or designs. Our approach, termed AQuN, combines quantile normalization with special attention to missing entities, to normalize miRNA expression across datasets, technologies, designs and platforms. As we show, our proposed approach uncovers patterns of interest that would not have emerged as statistically significant when analyzing the datasets individually or with other standard-practice normalization methods.

Amongst our findings, we noted a previously undocumented miRNA that is significantly over-expressed in samples from estrogen-receptor positive breast cancer patients as compared to samples from estrogen-receptor negative patients. We further investigated this miRNA, hsa-miR-193b-5p, and experimentally show, in cell lines, that its expression level impacts the viability of tumor cells. AQuN is available to the community in the form of Python and Matlab packages. The joint-processed data is also made available for further investigation.

CogNet: classification of gene expression data based on ranked active-subnetwork-oriented KEGG pathway enrichment analysis

Most of the traditional gene selection approaches are borrowed from other fields such as statistics and computer science, However, they do not prioritize biologically relevant genes since the ultimate goal is to determine features that optimize model performance metrics not to build a biologically meaningful model. Therefore, there is an imminent need for new computational tools that integrate the biological knowledge about the data in the process of gene selection and machine learning. Integrative gene selection enables incorporation of biological domain knowledge from external biological resources. In this study, we propose a new computational approach named CogNet that is an integrative gene selection tool that exploits biological knowledge for grouping the genes for the computational modeling tasks of ranking and classification. In CogNet, the pathfindR serves as the biological grouping tool to allow the main algorithm to rank active-subnetwork-oriented KEGG pathway enrichment analysis results to build a biologically relevant model. CogNet provides a list of significant KEGG pathways that can classify the data with a very high accuracy. The list also provides the genes belonging to these pathways that are differentially expressed that are used as features in the classification problem. The list facilitates deep analysis and better interpretability of the role of KEGG pathways in classification of the data thus better establishing the biological relevance of these differentially expressed genes. Even though the main aim of our study is not to improve the accuracy of any existing tool, the performance of the CogNet outperforms a similar approach called maTE while obtaining similar performance compared to other similar tools including SVM-RCE. CogNet was tested on 13 gene expression datasets concerning a variety of diseases.

maTE: discovering expressed interactions between microRNAs and their targets

MOTIVATION

Disease is often manifested via changes in transcript and protein abundance. MicroRNAs (miRNAs) are instrumental in regulating protein abundance and may measurably influence transcript levels. miRNAs often target more than one mRNA (for humans, the average is three), and mRNAs are often targeted by more than one miRNA (for the genes considered in this study, the average is also three). Therefore, it is difficult to determine the miRNAs that may cause the observed differential gene expression. We present a novel approach, maTE, which is based on machine learning, that integrates information about miRNA target genes with gene expression data. maTE depends on the availability of a sufficient amount of patient and control samples. The samples are used to train classifiers to accurately classify the samples on a per miRNA basis. Multiple high scoring miRNAs are used to build a final classifier to improve separation.

RESULTS

The aim of the study is to find a set of miRNAs causing the regulation of their target genes that best explains the difference between groups (e.g. cancer versus control). maTE provides a list of significant groups of genes where each group is targeted by a specific miRNA. For the datasets used in this study, maTE generally achieves an accuracy well above 80%. Also, the results show that when the accuracy is much lower (e.g. ∼50%), the set of miRNAs provided is likely not causative of the difference in expression. This new approach of integrating miRNA regulation with expression data yields powerful results and is independent of external labels and training data. Thereby, this approach allows new avenues for exploring miRNA regulation and may enable the development of miRNA-based biomarkers and drugs.

AVAILABILITY AND IMPLEMENTATION

The KNIME workflow, implementing maTE, is available at Bioinformatics online.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular harvesting with electroporation for tissue profiling

Recent developments in personalized medicine are based on molecular measurement steps that guide personally adjusted medical decisions. A central approach to molecular profiling consists of measuring DNA, RNA, and/or proteins in tissue samples, most notably in and around tumors. This measurement yields molecular biomarkers that are potentially predictive of response and of tumor type. Current methods in cancer therapy mostly use tissue biopsy as the starting point of molecular profiling. Tissue biopsies involve a physical resection of a small tissue sample, leading to localized tissue injury, bleeding, inflammation and stress, as well as to an increased risk of metastasis. Here we developed a technology for harvesting biomolecules from tissues using electroporation. We show that tissue electroporation, achieved using a combination of high-voltage short pulses, 50 pulses 500 V cm-1, 30 µs, 1 Hz, with low-voltage long pulses 50 pulses 50 V cm-1, 10 ms, delivered at 1 Hz, allows for tissue-specific extraction of RNA and proteins. We specifically tested RNA and protein extraction from excised kidney and liver samples and from excised HepG2 tumors in mice. Further in vivo development of extraction methods based on electroporation can drive novel approaches to the molecular profiling of tumors and of tumor environment and to related diagnosis practices.

A combined sequence and structure based method for discovering enriched motifs in RNA from in vivo binding data

RNA binding proteins (RBPs) play an important role in regulating many processes in the cell. RBPs often recognize their RNA targets in a specific manner. In addition to the RNA primary sequence, the structure of the RNA has been shown to play a central role in RNA recognition by RBPs. In recent years, many experimental approaches, both in vitro and in vivo, were developed and employed to identify and characterize RBP targets and extract their binding specificities. In vivo binding techniques, such as CrossLinking and ImmunoPrecipitation (CLIP)-based methods, enable the characterization of protein binding sites on RNA targets. However, these methods do not provide information regarding the structural preferences of the protein. While methods to obtain the structure of RNA are available, inferring both the sequence and the structure preferences of RBPs remains a challenge. Here we present SMARTIV, a novel computational tool for discovering combined sequence and structure binding motifs from in vivo RNA binding data relying on the sequences of the target sites, the ranking of their binding scores and their predicted secondary structure. The combined motifs are provided in a unified representation that is informative and easy for visual perception. We tested the method on CLIP-seq data from different platforms for a variety of RBPs. Overall, we show that our results are highly consistent with known binding motifs of RBPs, offering additional information on their structural preferences.