Qianlan Yao

Subpathway-GM: identification of metabolic subpathways via joint power of interesting genes and metabolites and their topologies within pathways

Various ‘omics’ technologies, including microarrays and gas chromatography mass spectrometry, can be used to identify hundreds of interesting genes, proteins and metabolites, such as differential genes, proteins and metabolites associated with diseases. Identifying metabolic pathways has become an invaluable aid to understanding the genes and metabolites associated with studying conditions. However, the classical methods used to identify pathways fail to accurately consider joint power of interesting gene/metabolite and the key regions impacted by them within metabolic pathways. In this study, we propose a powerful analytical method referred to as Subpathway-GM for the identification of metabolic subpathways. This provides a more accurate level of pathway analysis by integrating information from genes and metabolites, and their positions and cascade regions within the given pathway. We analyzed two colorectal cancer and one metastatic prostate cancer data sets and demonstrated that Subpathway-GM was able to identify disease-relevant subpathways whose corresponding entire pathways might be ignored using classical entire pathway identification methods. Further analysis indicated that the power of a joint genes/metabolites and subpathway strategy based on their topologies may play a key role in reliably recalling disease-relevant subpathways and finding novel subpathways.

The Implications of Relationships between Human Diseases and Metabolic Subpathways

One of the challenging problems in the etiology of diseases is to explore the relationships between initiation and progression of diseases and abnormalities in local regions of metabolic pathways. To gain insight into such relationships, we applied the “k-clique” subpathway identification method to all disease-related gene sets. For each disease, the disease risk regions of metabolic pathways were then identified and considered as subpathways associated with the disease. We finally built a disease-metabolic subpathway network (DMSPN). Through analyses based on network biology, we found that a few subpathways, such as that of cytochrome P450, were highly connected with many diseases, and most belonged to fundamental metabolisms, suggesting that abnormalities of fundamental metabolic processes tend to cause more types of diseases. According to the categories of diseases and subpathways, we tested the clustering phenomenon of diseases and metabolic subpathways in the DMSPN. The results showed that both disease nodes and subpathway nodes displayed slight clustering phenomenon. We also tested correlations between network topology and genes within disease-related metabolic subpathways, and found that within a disease-related subpathway in the DMSPN, the ratio of disease genes and the ratio of tissue-specific genes significantly increased as the number of diseases caused by the subpathway increased. Surprisingly, the ratio of essential genes significantly decreased and the ratio of housekeeping genes remained relatively unchanged. Furthermore, the coexpression levels between disease genes and other types of genes were calculated for each subpathway in the DMSPN. The results indicated that those genes intensely influenced by disease genes, including essential genes and tissue-specific genes, might be significantly associated with the disease diversity of subpathways, suggesting that different kinds of genes within a disease-related subpathway may play significantly differential roles on the diversity of diseases caused by the corresponding subpathway.

Topologically inferring risk-active pathways toward precise cancer classification by directed random walk

MOTIVATION The accurate prediction of disease status is a central challenge in clinical cancer research. Microarray-based gene biomarkers have been identified to predict outcome and outperform traditional clinical parameters. However, the robustness of the individual gene biomarkers is questioned because of their little reproducibility between different cohorts of patients. Substantial progress in treatment requires advances in methods to identify robust biomarkers. Several methods incorporating pathway information have been proposed to identify robust pathway markers and build classifiers at the level of functional categories rather than of individual genes. However, current methods consider the pathways as simple gene sets but ignore the pathway topological information, which is essential to infer a more robust pathway activity. RESULTS Here, we propose a directed random walk (DRW)-based method to infer the pathway activity. DRW evaluates the topological importance of each gene by capturing the structure information embedded in the directed pathway network. The strategy of weighting genes by their topological importance greatly improved the reproducibility of pathway activities. Experiments on 18 cancer datasets showed that the proposed method yielded a more accurate and robust overall performance compared with several existing gene-based and pathway-based classification methods. The resulting risk-active pathways are more reliable in guiding therapeutic selection and the development of pathway-specific therapeutic strategies. AVAILABILITY DRW is freely available at http://210.46.85.180:8080/DRWPClass/

Characterizing the Network of Drugs and Their Affected Metabolic Subpathways

A fundamental issue in biology and medicine is illustration of the overall drug impact which is always the consequence of changes in local regions of metabolic pathways (subpathways). To gain insights into the global relationship between drugs and their affected metabolic subpathways, we constructed a drug–metabolic subpathway network (DRSN). This network included 3925 significant drug–metabolic subpathway associations representing drug dual effects. Through analyses based on network biology, we found that if drugs were linked to the same subpathways in the DRSN, they tended to share the same indications and side effects. Furthermore, if drugs shared more subpathways, they tended to share more side effects. We then calculated the association score by integrating drug-affected subpathways and disease-related subpathways to quantify the extent of the associations between each drug class and disease class. The results showed some close drug–disease associations such as sex hormone drugs and cancer suggesting drug dual effects. Surprisingly, most drugs displayed close associations with their side effects rather than their indications. To further investigate the mechanism of drug dual effects, we classified all the subpathways in the DRSN into therapeutic and non-therapeutic subpathways representing drug therapeutic effects and side effects. Compared to drug side effects, the therapeutic effects tended to work through tissue-specific genes and these genes tend to be expressed in the adrenal gland, liver and kidney; while drug side effects always occurred in the liver, bone marrow and trachea. Taken together, the DRSN could provide great insights into understanding the global relationship between drugs and metabolic subpathways.

Identification of miRNA-Mediated Core Gene Module for Glioma Patient Prediction by Integrating High-Throughput miRNA, mRNA Expression and Pathway Structure

The prognosis of glioma patients is usually poor, especially in patients with glioblastoma (World Health Organization (WHO) grade IV). The regulatory functions of microRNA (miRNA) on genes have important implications in glioma cell survival. However, there are not many studies that have investigated glioma survival by integrating miRNAs and genes while also considering pathway structure. In this study, we performed sample-matched miRNA and mRNA expression profilings to systematically analyze glioma patient survival. During this analytical process, we developed pathway-based random walk to identify a glioma core miRNA-gene module, simultaneously considering pathway structure information and multi-level involvement of miRNAs and genes. The core miRNA-gene module we identified was comprised of four apparent sub-modules; all four sub-modules displayed a significant correlation with patient survival in the testing set (P-values≤0.001). Notably, one sub-module that consisted of 6 miRNAs and 26 genes also correlated with survival time in the high-grade subgroup (WHO grade III and IV), P-value = 0.0062. Furthermore, the 26-gene expression signature from this sub-module had robust predictive power in four independent, publicly available glioma datasets. Our findings suggested that the expression signatures, which were identified by integration of miRNA and gene level, were closely associated with overall survival among the glioma patients with various grades.

Dissection of miRNA-miRNA Interaction in Esophageal Squamous Cell Carcinoma

The relationships between miRNAs and their regulatory influences in esophageal carcinoma remain largely unknown. Accumulated evidence suggests that delineation of subpathways within an entire pathway can underlie complex diseases. To analyze the regulation of differentially expressed miRNAs in subpathways of esophageal squamous cell carcinoma (ESCC), we constructed bipartite miRNA and subpathway networks to determine miRNA regulatory influences on subpathways. The miRNA-subpathway network indicated that miRNAs regulate numerous subpathways. Two principal biological networks were derived from the miRNA-subpathway network by the hypergeometric test. This miRNA-miRNA network revealed the co-regulation of subpathways between the upregulated and downregulated miRNAs. Subpathway-subpathway networks characterized scale free, small world, and modular architecture. K-clique analysis revealed co-regulation of subpathways between certain downregulated and upregulated miRNAs. When ESCC patients were grouped according to their expression levels of paired upregulation of miR-31 and downregulation of miR-338-3p, survival time analysis revealed a significant difference based on miR-31-miR-338-3p interaction. These findings can facilitate the understanding of the biological meaning of miRNA-miRNA interactions with either the same or opposite expression trend.

Network based analyses of gene expression profile of LCN2 overexpression in esophageal squamous cell carcinoma

LCN2 (lipocalin 2) is a member of the lipocalin family of proteins that transport small, hydrophobic ligands. LCN2 is elevated in various cancers including esophageal squamous cell carcinoma (ESCC). In this study, LCN2 was overexpressed in the EC109 ESCC cell line and we applied integrated analyses of the gene expression data to identify protein-protein interactions (PPI) network to enhance our understanding of the role of LCN2 in ESCC. Through further mining of PPI sub-networks, hundreds of differentially expressed genes (DEGs) were identified to interact with thousands of other proteins. Subcellular localization analyses found the DEGs and their directly or indirectly interacting proteins distributed in multiple layers, which was applied to analyze the possible paths between two DEGs. Gene Ontology annotation generated a functional annotation map and found hundreds of significant terms, especially those associated with the known and potential roles of LCN2 protein. The algorithm of Random Walk with Restart was applied to prioritize the DEGs and identified several cancer-related DEGs ranked closest to LCN2 protein. These analyses based on PPI network have greatly expanded our understanding of the mRNA expression profile of LCN2 overexpresssion for future examination of the roles and mechanisms of LCN2.

Identifying disease related sub-pathways for analysis of genome-wide association studies

Most methods for genome-wide association studies (GWAS) focus on discovering a single genetic variant, but the pathogenesis of complex diseases is thought to arise from the joint effect of multiple genetic variants. Information about pathway structure, such as the interactions and distances between gene products within pathways, can help us learn more about the functions and joint effect of genes associated with disease risk. We developed a novel sub-pathway based approach to study the joint effect of multiple genetic variants that are modestly associated with disease. The approach prioritized sub-pathways based on the significance values of single nucleotide polymorphisms (SNPs) and the interactions and distances between gene products within pathways. We applied the method to seven complex diseases. The result showed that our method can efficiently identify statistically significant sub-pathways associated with the pathogenesis of complex diseases. The approach identified sub-pathways that may inform the interpretation of GWAS data.

A novel dysregulated pathway identification analysis based on global influence of within-pathway effects and crosstalk between pathways

Identifying dysregulated pathways from high-throughput experimental data in order to infer underlying biological insights is an important task. Current pathway-identification methods focus on single pathways in isolation; however, consideration of crosstalk between pathways could improve our understanding of alterations in biological states. We propose a novel method of pathway analysis based on global influence (PAGI) to identify dysregulated pathways, by considering both within-pathway effects and crosstalk between pathways. We constructed a global gene–gene network based on the relationships among genes extracted from a pathway database. We then evaluated the extent of differential expression for each gene, and mapped them to the global network. The random walk with restart algorithm was used to calculate the extent of genes affected by global influence. Finally, we used cumulative distribution functions to determine the significance values of the dysregulated pathways. We applied the PAGI method to five cancer microarray datasets, and compared our results with gene set enrichment analysis and five other methods. Based on these analyses, we demonstrated that PAGI can effectively identify dysregulated pathways associated with cancer, with strong reproducibility and robustness. We implemented PAGI using the freely available R-based and Web-based tools (http://bioinfo.hrbmu.edu.cn/PAGI).

Prioritizing Candidate Disease Metabolites Based on Global Functional Relationships between Metabolites in the Context of Metabolic Pathways

Identification of key metabolites for complex diseases is a challenging task in todays medicine and biology. A special disease is usually caused by the alteration of a series of functional related metabolites having a global influence on the metabolic network. Moreover, the metabolites in the same metabolic pathway are often associated with the same or similar disease. Based on these functional relationships between metabolites in the context of metabolic pathways, we here presented a pathway-based random walk method called PROFANCY for prioritization of candidate disease metabolites. Our strategy not only takes advantage of the global functional relationships between metabolites but also sufficiently exploits the functionally modular nature of metabolic networks. Our approach proved successful in prioritizing known metabolites for 71 diseases with an AUC value of 0.895. We also assessed the performance of PROFANCY on 16 disease classes and found that 4 classes achieved an AUC value over 0.95. To investigate the robustness of the PROFANCY, we repeated all the analyses in two metabolic networks and obtained similar results. Then we applied our approach to Alzheimers disease (AD) and found that a top ranked candidate was potentially related to AD but had not been reported previously. Furthermore, our method was applicable to prioritize the metabolites from metabolomic profiles of prostate cancer. The PROFANCY could identify prostate cancer related-metabolites that are supported by literatures but not considered to be significantly differential by traditional differential analysis. We also developed a freely accessible web-based and R-based tool at http://bioinfo.hrbmu.edu.cn/PROFANCY.