Non-linear tests for identifying differentially expressed genes or genetic networks

Computational Methods to Investigate the Impact of miRNAs on Pathways

Pathway analysis is a wide class of methods allowing to determine the alteration of functional processes in complex diseases. However, biological pathways are still partial, and knowledge coming from posttranscriptional regulators has started to be considered in a systematic way only recently. Here we will give a global and updated view of the main pathway and subpathway analysis methodologies, focusing on the improvements obtained through the recent introduction of microRNAs as regulatory elements in these frameworks.

Post-transcriptional knowledge in pathway analysis increases the accuracy of phenotypes classification

Motivation

Prediction of phenotypes from high-dimensional data is a crucial task in precision biology and medicine. Many technologies employ genomic biomarkers to characterize phenotypes. However, such elements are not sufficient to explain the underlying biology. To improve this, pathway analysis techniques have been proposed. Nevertheless, such methods have shown lack of accuracy in phenotypes classification.

Results

Here we propose a novel methodology called MITHrIL (Mirna enrIched paTHway Impact anaLysis) for the analysis of signaling pathways, which extends the work of Tarca et al., 2009. MITHrIL augments pathways with missing regulatory elements, such as microRNAs, and their interactions with genes. The method takes as input the expression values of genes and/or microRNAs and returns a list of pathways sorted according to their degree of deregulation, together with the corresponding statistical significance (p-values). Our analysis shows that MITHrIL outperforms its competitors even in the worst case. In addition, our method is able to correctly classify sets of tumor samples drawn from TCGA.

Availability

MITHrIL is freely available at the following URL: http://alpha.dmi.unict.it/mithril/

Detecting Disease Specific Pathway Substructures through an Integrated Systems Biology Approach

In the era of network medicine, pathway analysis methods play a central role in the prediction of phenotype from high throughput experiments. In this paper, we present a network-based systems biology approach capable of extracting disease-perturbed subpathways within pathway networks in connection with expression data taken from The Cancer Genome Atlas (TCGA). Our system extends pathways with missing regulatory elements, such as microRNAs, and their interactions with genes. The framework enables the extraction, visualization, and analysis of statistically significant disease-specific subpathways through an easy to use web interface. Our analysis shows that the methodology is able to fill the gap in current techniques, allowing a more comprehensive analysis of the phenomena underlying disease states.

A decision analysis model for KEGG pathway analysis

Background

The knowledge base-driven pathway analysis is becoming the first choice for many investigators, in that it not only can reduce the complexity of functional analysis by grouping thousands of genes into just several hundred pathways, but also can increase the explanatory power for the experiment by identifying active pathways in different conditions. However, current approaches are designed to analyze a biological system assuming that each pathway is independent of the other pathways.

Results

A decision analysis model is developed in this article that accounts for dependence among pathways in time-course experiments and multiple treatments experiments. This model introduces a decision coefficient—a designed index, to identify the most relevant pathways in a given experiment by taking into account not only the direct determination factor of each Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway itself, but also the indirect determination factors from its related pathways. Meanwhile, the direct and indirect determination factors of each pathway are employed to demonstrate the regulation mechanisms among KEGG pathways, and the sign of decision coefficient can be used to preliminarily estimate the impact direction of each KEGG pathway. The simulation study of decision analysis demonstrated the application of decision analysis model for KEGG pathway analysis.

Conclusions

A microarray dataset from bovine mammary tissue over entire lactation cycle was used to further illustrate our strategy. The results showed that the decision analysis model can provide the promising and more biologically meaningful results. Therefore, the decision analysis model is an initial attempt of optimizing pathway analysis methodology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1285-1) contains supplementary material, which is available to authorized users.

A de novo transcriptome analysis shows that modulation of the JAK-STAT signaling pathway by salmonid alphavirus subtype 3 favors virus replication in macrophage/dendritic-like TO-cells

BACKGROUND

The Janus kinase (Jak) and signaling transducer activator of transcription (Stat) pathway mediates the signaling of genes required for cellular development and homeostasis. To elucidate the effect of type I IFN on the Jak/stat pathway in salmonid alphavirus subtype 3 (SAV3) infected macrophage/dendritic like TO-cells derived from Atlantic salmon (Salmo salar L) headkidney leukocytes, we used a differential transcriptome analysis by RNA-seq and the Kyoto encyclopedia of genes and genomes (KEGGs) pathway analysis to generate a repertoire of de novo assembled genes from type I IFN treated and non-treated TO-cells infected with SAV3.

RESULTS

Concurrent SAV3 infection with type I IFN treatment of TO-cells suppressed SAV3 structural protein (SP) expression by 2log10 at 2 days post infection compared to SAV3 infection without IFN treatment which paved way to evaluating the impact of type I IFN on expression of Jak/stat pathway genes in SAV3 infected TO-cells. In the absence of type I IFN treatment, SAV3 downregulated several Jak/stat pathway genes that included type I and II receptor genes, Jak2, tyrosine kinase 2 (Tyk2), Stat3 and Stat5 pointing to possible failure to activate the Jak/stat signaling pathway and inhibition of signal transducers caused by SAV3 infection. Although the suppressor of cytokine signaling (SOCS) genes 1 and 3 were upregulated in the IFN treated cells, only SOCS3 was downregulated in the SAV3 infected cells which points to inhibition of SOCS3 by SAV3 infection in TO-cells.

CONCLUSION

Data presented in this study shows that SAV3 infection downregulates several genes of the Jak/stat pathway, which could be an immune evasion strategy, used to block the transcription of antiviral genes that would interfere with SAV3 replication in TO-cells. Overall, we have shown that combining de novo assembly with pathway based transcriptome analyses provides a contextual approach to understanding the molecular networks of genes that form the Jak/stat pathway in TO-cells infected by SAV3.

A flexible nonparametric approach to find candidate genes associated with disease in microarray experiments

Very often biologists are interested to know the biological function of a particular gene. Its true biological function may depend on other genes. Finding other genes in the same biological pathway of that gene may enhance further understanding of its biological function. Therefore, we are interested in finding other candidate genes whose expression values are highly correlated with that of a "seed" gene. The "seed" gene, which is known and associated with a disease, is used as a reference to extract candidate genes from microarray experiments and enriched pathways. We propose a nonparametric procedure for selecting the candidate genes. The expression levels for these candidate genes are correlated with that of a "seed" gene in microarray experiments. The proposed test statistic compares two Area Under Receiver Operating Characteristic Curves (AUC) for gene pairs, taking implicit correlation between two AUCs into account. The performance of our method is compared to the other well-known methods through the use of simulation and real data analysis.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Gene set analysis for self-contained tests: complex null and specific alternative hypotheses

MOTIVATION

The analysis of differentially expressed gene sets became a routine in the analyses of gene expression data. There is a multitude of tests available, ranging from aggregation tests that summarize gene-level statistics for a gene set to true multivariate tests, accounting for intergene correlations. Most of them detect complex departures from the null hypothesis but when the null hypothesis is rejected, the specific alternative leading to the rejection is not easily identifiable.

RESULTS

In this article we compare the power and Type I error rates of minimum-spanning tree (MST)-based non-parametric multivariate tests with several multivariate and aggregation tests, which are frequently used for pathway analyses. In our simulation study, we demonstrate that MST-based tests have power that is for many settings comparable with the power of conventional approaches, but outperform them in specific regions of the parameter space corresponding to biologically relevant configurations. Further, we find for simulated and for gene expression data that MST-based tests discriminate well against shift and scale alternatives. As a general result, we suggest a two-step practical analysis strategy that may increase the interpretability of experimental data: first, apply the most powerful multivariate test to find the subset of pathways for which the null hypothesis is rejected and second, apply MST-based tests to these pathways to select those that support specific alternative hypotheses.

CONTACT

gvglazko@uams.edu or yrahmatallah@uams.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Assessment method for a power analysis to identify differentially expressed pathways

Gene expression data can provide a very rich source of information for elucidating the biological function on the pathway level if the experimental design considers the needs of the statistical analysis methods. The purpose of this paper is to provide a comparative analysis of statistical methods for detecting the differentially expression of pathways (DEP). In contrast to many other studies conducted so far, we use three novel simulation types, producing a more realistic correlation structure than previous simulation methods. This includes also the generation of surrogate data from two large-scale microarray experiments from prostate cancer and ALL. As a result from our comprehensive analysis of parameter configurations, we find that each method should only be applied if certain conditions of the data from a pathway are met. Further, we provide method-specific estimates for the optimal sample size for microarray experiments aiming to identify DEP in order to avoid an underpowered design. Our study highlights the sensitivity of the studied methods on the parameters of the system.

Ten years of pathway analysis: current approaches and outstanding challenges

Pathway analysis has become the first choice for gaining insight into the underlying biology of differentially expressed genes and proteins, as it reduces complexity and has increased explanatory power. We discuss the evolution of knowledge base–driven pathway analysis over its first decade, distinctly divided into three generations. We also discuss the limitations that are specific to each generation, and how they are addressed by successive generations of methods. We identify a number of annotation challenges that must be addressed to enable development of the next generation of pathway analysis methods. Furthermore, we identify a number of methodological challenges that the next generation of methods must tackle to take advantage of the technological advances in genomics and proteomics in order to improve specificity, sensitivity, and relevance of pathway analysis.

GENETIC NETWORK ANALYSIS BY QUASI-BAYESIAN METHOD

Genetic network analysis provides an important statistical strategy for the study of gene–gene interactions. Although existing methods work well in practice, several opportunities for improvement remain. For example, the regulation coefficients of some of the existing methods are not easy to solve, nor are the solutions they provide unique. Also, as genetic network analysis are typically applied to small datasets with large number of parameters, having prior knowledge about the parameters is valuable and should be incorporated into the analysis. The uniqueness of the parameter estimate and computational simplicity are also desirable in practice. To address these problems, we considered a quasi-Bayesian method for the analysis of gene regulatory networks by a multivariate linear model in which the data distribution is a quasi-likelihood, and the inference is Bayesian. This method incorporates prior information on the regulatory relationships; the set of regulation coefficients has a unique closed-form solution, and is very simple to compute. The model is evaluated by simulation and illustrated using a real dataset. This method is simple to use, permits information updating, is flexible to incorporate desired features, and has closed-form solution. Simulation studies show that the model fits the data quite well.

Unite and conquer: univariate and multivariate approaches for finding differentially expressed gene sets

MOTIVATION

Recently, many univariate and several multivariate approaches have been suggested for testing differential expression of gene sets between different phenotypes. However, despite a wealth of literature studying their performance on simulated and real biological data, still there is a need to quantify their relative performance when they are testing different null hypotheses.

RESULTS

In this article, we compare the performance of univariate and multivariate tests on both simulated and biological data. In the simulation study we demonstrate that high correlations equally affect the power of both, univariate as well as multivariate tests. In addition, for most of them the power is similarly affected by the dimensionality of the gene set and by the percentage of genes in the set, for which expression is changing between two phenotypes. The application of different test statistics to biological data reveals that three statistics (sum of squared t-tests, Hotelling's T(2), N-statistic), testing different null hypotheses, find some common but also some complementing differentially expressed gene sets under specific settings. This demonstrates that due to complementing null hypotheses each test projects on different aspects of the data and for the analysis of biological data it is beneficial to use all three tests simultaneously instead of focusing exclusively on just one.

Identification of differential gene pathways with principal component analysis

MOTIVATION

Development of high-throughput technology makes it possible to measure expressions of thousands of genes simultaneously. Genes have the inherent pathway structure, where pathways are composed of multiple genes with coordinated biological functions. It is of great interest to identify differential gene pathways that are associated with the variations of phenotypes.

RESULTS

We propose the following approach for detecting differential gene pathways. First, we construct gene pathways using databases such as KEGG or GO. Second, for each pathway, we extract a small number of representative features, which are linear combinations of gene expressions and/or their transformations. Specifically, we propose using (i) principal components (PCs) of gene expression sets, (ii) PCs of expanded gene expression sets and (iii) expanded sets of PCs of gene expressions, as the representative features. Third, we identify differential gene pathways as those with representative features significantly associated with the variations of phenotypes, particularly disease clinical outcomes, in regression models. The false discovery rate approach is used to adjust for multiple comparisons. Analysis of three gene expression datasets suggests that (i) the proposed approach can effectively identify differential gene pathways; (ii) PCs that explain only a small amount of variations of gene expressions may bear significant associations between gene pathways and phenotypes; (iii) including second-order terms of gene expressions may lead to identification of new differential gene pathways; (iv) the proposed approach is relatively insensitive to additional noises; and (v) the proposed approach can identify gene pathways missed by alternative approaches.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Construction of a reference gene association network from multiple profiling data: application to data analysis

MOTIVATION

Gene expression profiling is an important tool for gaining insight into biology. Novel strategies are required to analyze the growing archives of microarray data and extract useful information from them. One area of interest is in the construction of gene association networks from collections of profiling data. Various approaches have been proposed to construct gene networks using profiling data, and these networks have been used in functional inference as well as in data visualization. Here, we investigated a non-parametric approach to translate profiling data into a gene network. We explored the characteristics and utility of the resulting network and investigated the use of network information in analysis of variance models and hypothesis testing.

RESULTS

Our work is composed of two parts: gene network construction and partitioning and hypothesis testing using sub-networks as groups. In the first part, multiple independently collected microarray datasets from the Gene Expression Omnibus data repository were analyzed to identify probe pairs that are positively co-regulated across the samples. A co-expression network was constructed based on a reciprocal ranking criteria and a false discovery rate analysis. We named this network Reference Gene Association (RGA) network. Then, the network was partitioned into densely connected sub-networks of probes using a multilevel graph partitioning algorithm. In the second part, we proposed a new, MANOVA-based approach that can take individual probe expression values as input and perform hypothesis testing at the sub-network level. We applied this MANOVA methodology to two published studies and our analysis indicated that the methodology is both effective and sensitive for identifying transcriptional sub-networks or pathways that are perturbed across treatments.