A rank-based statistical test for measuring synergistic effects between two gene sets

Integrated systems biology analysis of KSHV latent infection reveals viral induction and reliance on peroxisome mediated lipid metabolism

Kaposi’s Sarcoma associated Herpesvirus (KSHV), an oncogenic, human gamma-herpesvirus, is the etiological agent of Kaposi’s Sarcoma the most common tumor of AIDS patients world-wide. KSHV is predominantly latent in the main KS tumor cell, the spindle cell, a cell of endothelial origin. KSHV modulates numerous host cell-signaling pathways to activate endothelial cells including major metabolic pathways involved in lipid metabolism. To identify the underlying cellular mechanisms of KSHV alteration of host signaling and endothelial cell activation, we identified changes in the host proteome, phosphoproteome and transcriptome landscape following KSHV infection of endothelial cells. A Steiner forest algorithm was used to integrate the global data sets and, together with transcriptome based predicted transcription factor activity, cellular networks altered by latent KSHV were predicted. Several interesting pathways were identified, including peroxisome biogenesis. To validate the predictions, we showed that KSHV latent infection increases the number of peroxisomes per cell. Additionally, proteins involved in peroxisomal lipid metabolism of very long chain fatty acids, including ABCD3 and ACOX1, are required for the survival of latently infected cells. In summary, novel cellular pathways altered during herpesvirus latency that could not be predicted by a single systems biology platform, were identified by integrated proteomics and transcriptomics data analysis and when correlated with our metabolomics data revealed that peroxisome lipid metabolism is essential for KSHV latent infection of endothelial cells.

Kaposi’s Sarcoma herpesvirus (KSHV) is the etiologic agent of Kaposi’s Sarcoma, the most common tumor of AIDS patients. KSHV modulates host cell signaling and metabolism to maintain a life-long latent infection. To unravel the underlying cellular mechanisms modulated by KSHV, we used multiple global systems biology platforms to identify and integrate changes in both cellular protein expression and transcription following KSHV infection of endothelial cells, the relevant cell type for KS tumors. The analysis identified several interesting pathways including peroxisome biogenesis. Peroxisomes are small cytoplasmic organelles involved in redox reactions and lipid metabolism. KSHV latent infection increases the number of peroxisomes per cell and proteins involved in peroxisomal lipid metabolism are required for the survival of latently infected cells. In summary, through integration of multiple global systems biology analyses we were able to identify novel pathways that could not be predicted by one platform alone and found that lipid metabolism in a small cytoplasmic organelle is necessary for the survival of latent infection with a herpesvirus.

Discovering Pair-wise Synergies in Microarray Data

Informative gene selection can have important implications for the improvement of cancer diagnosis and the identification of new drug targets. Individual-gene-ranking methods ignore interactions between genes. Furthermore, popular pair-wise gene evaluation methods, e.g. TSP and TSG, are helpless for discovering pair-wise interactions. Several efforts to discover pair-wise synergy have been made based on the information approach, such as EMBP and FeatKNN. However, the methods which are employed to estimate mutual information, e.g. binarization, histogram-based and KNN estimators, depend on known data or domain characteristics. Recently, Reshef et al. proposed a novel maximal information coefficient (MIC) measure to capture a wide range of associations between two variables that has the property of generality. An extension from MIC(X; Y) to MIC(X₁; X₂; Y) is therefore desired. We developed an approximation algorithm for estimating MIC(X₁; X₂; Y) where Y is a discrete variable. MIC(X₁; X₂; Y) is employed to detect pair-wise synergy in simulation and cancer microarray data. The results indicate that MIC(X₁; X₂; Y) also has the property of generality. It can discover synergic genes that are undetectable by reference feature selection methods such as MIC(X; Y) and TSG. Synergic genes can distinguish different phenotypes. Finally, the biological relevance of these synergic genes is validated with GO annotation and OUgene database.

Statistical significance of combinatorial regulations

More than three transcription factors often work together to enable cells to respond to various signals. The detection of combinatorial regulation by multiple transcription factors, however, is not only computationally nontrivial but also extremely unlikely because of multiple testing correction. The exponential growth in the number of tests forces us to set a strict limit on the maximum arity. Here, we propose an efficient branch-and-bound algorithm called the "limitless arity multiple-testing procedure" (LAMP) to count the exact number of testable combinations and calibrate the Bonferroni factor to the smallest possible value. LAMP lists significant combinations without any limit, whereas the family-wise error rate is rigorously controlled under the threshold. In the human breast cancer transcriptome, LAMP discovered statistically significant combinations of as many as eight binding motifs. This method may contribute to uncover pathways regulated in a coordinated fashion and find hidden associations in heterogeneous data.