Identification of key genes and therapeutic drugs for cocaine addiction using integrated bioinformatics analysis.

Introduction: Cocaine is a highly addictive drug that is abused due to its excitatory effect on the central nervous system. It is critical to reveal the mechanisms of cocaine addiction and identify key genes that play an important role in addiction. Methods: In this study, we proposed a centrality algorithm integration strategy to identify key genes in a protein-protein interaction (PPI) network constructed by deferential genes from cocaine addiction-related datasets. In order to investigate potential therapeutic drugs for cocaine addiction, a network of targeted relationships between nervous system drugs and key genes was established. Results: Four key genes (JUN, FOS, EGR1, and IL6) were identified and well validated using CTD database correlation analysis, text mining, independent dataset analysis, and enrichment analysis methods, and they might serve as biomarkers of cocaine addiction. A total of seventeen drugs have been identified from the network of targeted relationships between nervous system drugs and key genes, of which five (disulfiram, cannabidiol, dextroamphetamine, diazepam, and melatonin) have been shown in the literature to play a role in the treatment of cocaine addiction. Discussion: This study identified key genes and potential therapeutic drugs for cocaine addiction, which provided new ideas for the research of the mechanism of cocaine addiction.

Guiding Drug Repositioning for Cancers Based on Drug Similarity Networks.

Drug repositioning aims to discover novel clinical benefits of existing drugs, is an effective way to develop drugs for complex diseases such as cancer and may facilitate the process of traditional drug development. Meanwhile, network-based computational biology approaches, which allow the integration of information from different aspects to understand the relationships between biomolecules, has been successfully applied to drug repurposing. In this work, we developed a new strategy for network-based drug repositioning against cancer. Combining the mechanism of action and clinical efficacy of the drugs, a cancer-related drug similarity network was constructed, and the correlation score of each drug with a specific cancer was quantified. The top 5% of scoring drugs were reviewed for stability and druggable potential to identify potential repositionable drugs. Of the 11 potentially repurposable drugs for non-small cell lung cancer (NSCLC), 10 were confirmed by clinical trial articles and databases. The targets of these drugs were significantly enriched in cancer-related pathways and significantly associated with the prognosis of NSCLC. In light of the successful application of our approach to colorectal cancer as well, it provides an effective clue and valuable perspective for drug repurposing in cancer.

Leader gene identification for digestive system cancers based on human subcellular location and cancer-related characteristics in protein-protein interaction networks.

Stomach, liver, and colon cancers are the most common digestive system cancers leading to mortality. Cancer leader genes were identified in the current study as the genes that contribute to tumor initiation and could shed light on the molecular mechanisms in tumorigenesis. An integrated procedure was proposed to identify cancer leader genes based on subcellular location information and cancer-related characteristics considering the effects of nodes on their neighbors in human protein-protein interaction networks. A total of 69, 43, and 64 leader genes were identified for stomach, liver, and colon cancers, respectively. Furthermore, literature reviews and experimental data including protein expression levels and independent datasets from other databases all verified their association with corresponding cancer types. These final leader genes were expected to be used as diagnostic biomarkers and targets for new treatment strategies. The procedure for identifying cancer leader genes could be expanded to open up a window into the mechanisms, early diagnosis, and treatment of other cancer types.

Integrating Expression Data-Based Deep Neural Network Models with Biological Networks to Identify Regulatory Modules for Lung Adenocarcinoma.

Lung adenocarcinoma is the most common type of primary lung cancer, but the regulatory mechanisms during carcinogenesis remain unclear. The identification of regulatory modules for lung adenocarcinoma has become one of the hotspots of bioinformatics. In this paper, multiple deep neural network (DNN) models were constructed using the expression data to identify regulatory modules for lung adenocarcinoma in biological networks. First, the mRNAs, lncRNAs and miRNAs with significant differences in the expression levels between tumor and non-tumor tissues were obtained. MRNA DNN models were established and optimized to mine candidate mRNAs that significantly contributed to the DNN models and were in the center of an interaction network. Another DNN model was then constructed and potential ceRNAs were screened out based on the contribution of each RNA to the model. Finally, three modules comprised of miRNAs and their regulated mRNAs and lncRNAs with the same regulation direction were identified as regulatory modules that regulated the initiation of lung adenocarcinoma through ceRNAs relationships. They were validated by literature and functional enrichment analysis. The effectiveness of these regulatory modules was evaluated in an independent lung adenocarcinoma dataset. Regulatory modules for lung adenocarcinoma identified in this study provided a reference for regulatory mechanisms during carcinogenesis.

Identification of cardiomyopathy-related core genes through human metabolic networks and expression data.

BACKGROUND: Cardiomyopathy is a complex type of myocardial disease, and its incidence has increased significantly in recent years. Dilated cardiomyopathy (DCM) and ischemic cardiomyopathy (ICM) are two common and indistinguishable types of cardiomyopathy. RESULTS: Here, a systematic multi-omics integration approach was proposed to identify cardiomyopathy-related core genes that could distinguish normal, DCM and ICM samples using cardiomyopathy expression profile data based on a human metabolic network. First, according to the differentially expressed genes between different states (DCM/ICM and normal, or DCM and ICM) of samples, three sets of initial modules were obtained from the human metabolic network. Two permutation tests were used to evaluate the significance of the Pearson correlation coefficient difference score of the initial modules, and three candidate modules were screened out. Then, a cardiomyopathy risk module that was significantly related to DCM and ICM was determined according to the significance of the module score based on Markov random field. Finally, based on the shortest path between cardiomyopathy known genes, 13 core genes related to cardiomyopathy were identified. These core genes were enriched in pathways and functions significantly related to cardiomyopathy and could distinguish between samples of different states. CONCLUSION: The identified core genes might serve as potential biomarkers of cardiomyopathy. This research will contribute to identifying potential biomarkers of cardiomyopathy and to distinguishing different types of cardiomyopathy.

Candidate gene prioritization for chronic obstructive pulmonary disease using expression information in protein-protein interaction networks.

BACKGROUND: Identifying or prioritizing genes for chronic obstructive pulmonary disease (COPD), one type of complex disease, is particularly important for its prevention and treatment. METHODS: In this paper, a novel method was proposed to Prioritize genes using Expression information in Protein-protein interaction networks with disease risks transferred between genes (abbreviated as PEP). A weighted COPD PPI network was constructed using expression information and then COPD candidate genes were prioritized based on their corresponding disease risk scores in descending order. RESULTS: Further analysis demonstrated that the PEP method was robust in prioritizing disease candidate genes, and superior to other existing prioritization methods exploiting either topological or functional information. Top-ranked COPD candidate genes and their significantly enriched functions were verified to be related to COPD. The top 200 candidate genes might be potential disease genes in the diagnosis and treatment of COPD. CONCLUSIONS: The proposed method could provide new insights to the research of prioritizing candidate genes of COPD or other complex diseases with expression information from sequencing or microarray data.

Identification of genes associated with cancer prognosis in glioma: In silico data analyses of Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA).

Glioma is one particular type of brain malignancy which is highly complex and usually has a poor prognosis. Despite the limited diagnostic level of glioma, the survival time of affected patients broadly varies. Here, we conducted a detailed analysis, regarding the differences in patient survival time, to discover potential survival-related genes in glioma as well as their putative regulatory mechanisms. To contextualize the acquisition of these potential prognosis markers in large populations, particularly in China, we combined CGGA and The Cancer Genome Atlas (TCGA) databases to properly identify genes that are significantly related to survival. Our workflow combined a series of analytical approaches, including differential analysis, survival time, co-expression, clinical correlation analysis, ROC curve evaluation and prediction ability. Our results indicate that the four particular genes - PLAT, IGFBP2, BCAT1, SERPINH1 could be used as independent prognostic marker genes. These genes have also shown good prognostic ability in distinct populations, reiterating the robustness and value of these prognostic markers.

Human Microbe-Disease Association Prediction by a Novel Double-Ended Random Walk with Restart.

Microorganisms in the human body play a vital role in metabolism, immune defense, nutrient absorption, cancer control, and prevention of pathogen colonization. More and more biological and clinical studies have shown that the imbalance of microbial communities is closely related to the occurrence and development of various complex human diseases. Finding potential microbial-disease associations is critical for understanding the pathology of a few diseases and thus further improving disease diagnosis and prognosis. In this study, we proposed a novel computational model to predict disease-associated microbes. Specifically, we first constructed a heterogeneous interconnection network based on known microbe-disease associations deposited in a few databases, the similarity between diseases, and the similarity between microorganisms. We then predicted novel microbe-disease associations by a new method called the double-ended restart random walk model (DRWHMDA) implemented on the interconnection network. In addition, we performed case studies of colon cancer and asthma for further evaluation. The results indicate that 10 and 9 of the top 10 microorganisms predicted to be associated with colorectal cancer and asthma were validated by relevant literatures, respectively. Our method is expected to be effective in identifying disease-related microorganisms and will help to reveal the relationship between microorganisms and complex human diseases.

Network-based integration method for potential breast cancer gene identification.

Breast cancer is the most common female death-causing cancer worldwide. A network-based integration method was proposed to identify potential breast cancer genes. First, genes were prioritized using a gene prioritization algorithm by the strategy of disease risks transferred between genes in a network with weighted vertexes and edges. Our prioritization algorithm was effectives and robust for top-ranked seed gene number and higher area under the curve values compared to ToppGene and ToppNet. Then, 20 potential breast cancer genes were identified as common genes of the top 50 candidate genes for their robustness in multiple prioritizations. These genes could accurately classify tumor and normal samples of all and paired sample sets and three independent datasets. Of potential breast cancer genes, 18 were verified by literature and 2 were novel genes that need further study. This study would contribute to the understanding of the genetic architecture for the diagnosis and treatment of breast cancer.

A susceptibility biomarker identification strategy based on significantly differentially expressed ceRNA triplets for ischemic cardiomyopathy.

Ischemic cardiomyopathy (ICM) is a common human heart disease that causes death. No effective biomarkers for ICM could be found in existing databases, which is detrimental to the in-depth study of this disease. In the present study, ICM susceptibility biomarkers were identified using a proposed strategy based on RNA-Seq and miRNA-Seq data of ICM and normal samples. Significantly differentially expressed competing endogenous RNA (ceRNA) triplets were constructed using permutation tests and differentially expressed mRNAs, miRNAs and lncRNAs. Candidate ICM susceptible genes were screened out as differentially expressed genes in significantly differentially expressed ceRNA triplets enriched in ICM-related functional classes. Finally, eight ICM susceptibility genes and their significantly correlated lncRNAs with high classification accuracy were identified as ICM susceptibility biomarkers. These biomarkers would contribute to the diagnosis and treatment of ICM. The proposed strategy could be extended to other complex diseases without disease biomarkers in public databases.

Identification of breast cancer risk modules via an integrated strategy.

Breast cancer is one of the most common malignant cancers among females worldwide. This complex disease is not caused by a single gene, but resulted from multi-gene interactions, which could be represented by biological networks. Network modules are composed of genes with significant similarities in terms of expression, function and disease association. Therefore, the identification of disease risk modules could contribute to understanding the molecular mechanisms underlying breast cancer. In this paper, an integrated disease risk module identification strategy was proposed according to a multi-objective programming model for two similarity criteria as well as significance of permutation tests in Markov random field module score, function consistency score and Pearson correlation coefficient difference score. Three breast cancer risk modules were identified from a breast cancer-related interaction network. Genes in these risk modules were confirmed to play critical roles in breast cancer by literature review. These risk modules were enriched in breast cancer-related pathways or functions and could distinguish between breast tumor and normal samples with high accuracy for not only the microarray dataset used for breast cancer risk module identification, but also another two independent datasets. Our integrated strategy could be extended to other complex diseases to identify their risk modules and reveal their pathogenesis.

Association of PD-1 polymorphisms with the risk and prognosis of lung adenocarcinoma in the northeastern Chinese Han population.

BACKGROUND: Lung cancer is a leading cause of death from cancer worldwide, especially non-small cell lung cancer (NSCLC). The marker of progression in lung adenocarcinoma, the main type of NSCLC, has been rarely studied. Programmed death 1 (PD-1) is an effective drug target for the treatment of NSCLC. Our study aimed to examine the PD-1 role in the disease process. The study of the effect of polymorphisms on the progression of lung adenocarcinoma in the Han population of Northeast China may provide a valuable reference for the research and application of these drugs. METHODS: Chi-square test, Wilcoxon rank sum test, and classification efficiency assessment were used to test SNPs of PD-1 in 287 patients and combined with clinical information. RESULTS: We successfully identified biomarkers (rs2227981, rs2227982, and rs3608432) that could distinguish between lung adenocarcinoma patients of early stages and late stages. Multiple clinical indicators showed significant differences among different SNPs and cancer stages. Furthermore, this gene was confirmed to effectively distinguish the stages of lung adenocarcinoma with RNA-seq data in TCGA. CONCLUSIONS: Out study indicated that the PD-1 gene and the SNPs on it could be used as markers for distinguishing lung adenocarcinoma staging in the Northeast Han population. Our investigation into the link between PD-1 polymorphisms and lung adenocarcinoma would help to provide guidance for the treatment and prognosis of lung adenocarcinoma.

Candidate gene prioritization for non-communicable diseases based on functional information: Case studies.

Candidate gene prioritization for complex non-communicable diseases is essential to understanding the mechanism and developing better means for diagnosing and treating these diseases. Many methods have been developed to prioritize candidate genes in protein-protein interaction (PPI) networks. Integrating functional information/similarity into disease-related PPI networks could improve the performance of prioritization. In this study, a candidate gene prioritization method was proposed for non-communicable diseases considering disease risks transferred between genes in weighted disease PPI networks with weights for nodes and edges based on functional information. Here, three types of non-communicable diseases with pathobiological similarity, Type 2 diabetes (T2D), coronary artery disease (CAD) and dilated cardiomyopathy (DCM), were used as case studies. Literature review and pathway enrichment analysis of top-ranked genes demonstrated the effectiveness of our method. Better performance was achieved after comparing our method with other existing methods. Pathobiological similarity among these three diseases was further investigated for common top-ranked genes to reveal their pathogenesis.

Prioritizing complex disease risk genes by integrating multiple data.

Complex diseases, such as obesity, type II diabetes and chronic obstructive pulmonary disease (COPD) as metabolic disorder-related diseases are major concern for worldwide public health in the 21st century. The identification of these disease risk genes has attracted increasing interest in computational systems biology. In this paper, a novel method was proposed to prioritize disease risk genes (PDRG) by integrating functional annotations, protein interactions and gene expression information to assess similarity between genes in a disease-related metabolic network. The gene prioritization method was successfully carried out for obesity and COPD, the effectiveness of which was superior to those of ToppGene and ToppNet in both literature validation and recall rate by LOOCV. Our method could be applied broadly to other metabolism-related diseases, helping to prioritize novel disease risk genes, and could shed light on diagnosis and effective therapies.

The Pathogenesis of Atherosclerosis Based on Human Signaling Networks and Stem Cell Expression Data.

Atherosclerosis is a common and complex disease, whose morbidity increased significantly. Here, an integrated approach was proposed to elucidate systematically the pathogenesis of atherosclerosis from a systems biology point of view. Two weighted human signaling networks were constructed based on atherosclerosis related gene expression data of stem cells. Then, 37 candidate Atherosclerosis-risk Modules were detected using four kinds of permutation tests. Five Atherosclerosis-risk Modules (three Absent Modules and two Emerging Modules) enriched in functions significantly associated with disease genes were identified and verified to be associated with the maintenance of normal biological process and the pathogenesis and development of atherosclerosis. Especially for Atherosclerosis-risk Emerging Module P96, it could distinguish between normal and disease samples by Supporting Vector Machine with the average expression value of the module as classification feature. These identified modules and their genes may act as potential atherosclerosis biomarkers. Our study would shed light on the signal transduction of atherosclerosis, and provide new insights to its pathogenesis from the perspective of stem cells.

Prioritizing chronic obstructive pulmonary disease (COPD) candidate genes in COPD-related networks.

Chronic obstructive pulmonary disease (COPD) is a multi-factor disease, which could be caused by many factors, including disturbances of metabolism and protein-protein interactions (PPIs). In this paper, a weighted COPD-related metabolic network and a weighted COPD-related PPI network were constructed base on COPD disease genes and functional information. Candidate genes in these weighted COPD-related networks were prioritized by making use of a gene prioritization method, respectively. Literature review and functional enrichment analysis of the top 100 genes in these two networks suggested the correlation of COPD and these genes. The performance of our gene prioritization method was superior to that of ToppGene and ToppNet for genes from the COPD-related metabolic network or the COPD-related PPI network after assessing using leave-one-out cross-validation, literature validation and functional enrichment analysis. The top-ranked genes prioritized from COPD-related metabolic and PPI networks could promote the better understanding about the molecular mechanism of this disease from different perspectives. The top 100 genes in COPD-related metabolic network or COPD-related PPI network might be potential markers for the diagnosis and treatment of COPD.

Chronic obstructive pulmonary disease candidate gene prioritization based on metabolic networks and functional information.

Chronic obstructive pulmonary disease (COPD) is a multi-factor disease, in which metabolic disturbances played important roles. In this paper, functional information was integrated into a COPD-related metabolic network to assess similarity between genes. Then a gene prioritization method was applied to the COPD-related metabolic network to prioritize COPD candidate genes. The gene prioritization method was superior to ToppGene and ToppNet in both literature validation and functional enrichment analysis. Top-ranked genes prioritized from the metabolic perspective with functional information could promote the better understanding about the molecular mechanism of this disease. Top 100 genes might be potential markers for diagnostic and effective therapies.

Identification of susceptible genes for complex chronic diseases based on disease risk functional SNPs and interaction networks.

Complex chronic diseases are caused by the effects of genetic and environmental factors. Single nucleotide polymorphisms (SNPs), one common type of genetic variations, played vital roles in diseases. We hypothesized that disease risk functional SNPs in coding regions and protein interaction network modules were more likely to contribute to the identification of disease susceptible genes for complex chronic diseases. This could help to further reveal the pathogenesis of complex chronic diseases. Disease risk SNPs were first recognized from public SNP data for coronary heart disease (CHD), hypertension (HT) and type 2 diabetes (T2D). SNPs in coding regions that were classified into nonsense and missense by integrating several SNP functional annotation databases were treated as functional SNPs. Then, regions significantly associated with each disease were screened using random permutations for disease risk functional SNPs. Corresponding to these regions, 155, 169 and 173 potential disease susceptible genes were identified for CHD, HT and T2D, respectively. A disease-related gene product interaction network in environmental context was constructed for interacting gene products of both disease genes and potential disease susceptible genes for these diseases. After functional enrichment analysis for disease associated modules, 5 CHD susceptible genes, 7 HT susceptible genes and 3 T2D susceptible genes were finally identified, some of which had pleiotropic effects. Most of these genes were verified to be related to these diseases in literature. This was similar for disease genes identified from another method proposed by Lee et al. from a different aspect. This research could provide novel perspectives for diagnosis and treatment of complex chronic diseases and susceptible genes identification for other diseases.

Identification of breast cancer prognostic modules based on weighted protein-protein interaction networks.

Breast cancer is one of the leading causes of mortality in females. A number of prognostic markers have been identified, including single genes, multi-gene signatures and network modules; however, the robustness of these prognostic markers is insufficient. Thus, the present study proposed a more robust method to identify breast cancer prognostic modules based on weighted protein-protein interaction networks, by integrating four sets of disease-associated expression profiles. Three identified prognostic modules were closely associated with prognosis-associated functions and survival time, as determined by Cox regression and Kaplan-Meier survival analyses. The robustness of these modules was verified with an independent profile from another platform. Genes from these modules may be useful as breast cancer prognostic markers. The prognostic modules could be used to determine the prognoses of patients with breast cancer and characterize patient recovery.