Motivation: The emergence of network medicine not only offers more opportunities for better and more complete understanding of the molecular complexities of diseases, but also serves as a promising tool for identifying new drug targets and establishing new relationships among diseases that enable drug repositioning. Computational approaches for drug repositioning by integrating information from multiple sources and multiple levels have the potential to provide great insights to the complex relationships among drugs, targets, disease genes and diseases at a system level. Results: In this article, we have proposed a computational framework based on a heterogeneous network model and applied the approach on drug repositioning by using existing omics data about diseases, drugs and drug targets. The novelty of the framework lies in the fact that the strength between a disease–drug pair is calculated through an iterative algorithm on the heterogeneous graph that also incorporates drug-target information. Comprehensive experimental results show that the proposed approach significantly outperforms several recent approaches. Case studies further illustrate its practical usefulness . Availability and implementation:http://cbc.case.edu Contact:jingli@cwru.edu Supplementary information: Supplementary data are available at Bioinformatics online.

To identify factors that regulate gut microbiota density and the impact of varied microbiota density on health, we assayed this fundamental ecosystem property in fecal samples across mammals, human disease, and therapeutic interventions. Physiologic features of the host (carrying capacity) and the fitness of the gut microbiota shape microbiota density. Therapeutic manipulation of microbiota density in mice altered host metabolic and immune homeostasis. In humans, gut microbiota density was reduced in Crohn’s disease, ulcerative colitis, and ileal pouch-anal anastomosis. The gut microbiota in recurrent Clostridium difficile infection had lower density and reduced fitness that were restored by fecal microbiota transplantation. Understanding the interplay between microbiota and disease in terms of microbiota density, host carrying capacity, and microbiota fitness provide new insights into microbiome structure and microbiome targeted therapeutics.

A key issue in drug development is to understand the hidden relationships among drugs and targets. Computational methods for novel drug target predictions can greatly reduce time and costs compared with experimental methods. In this paper, we propose a network based computational approach for novel drug and target association predictions. More specifically, a heterogeneous drug-target graph, which incorporates known drug-target interactions as well as drug-drug and target-target similarities, is first constructed. Based on this graph, a novel graph-based inference method is introduced. Compared with two state-of-the-art methods, large-scale cross-validation results indicate that the proposed method can greatly improve novel target predictions.

Background: Patient stratification based on molecular subtypes is an important strategy for cancer precision medicine. Deriving clinically informative cancer molecular subtypes from transcriptomic data generated on whole tumor tissue samples is a non-trivial task, especially given the various non-cancer cellular elements intertwined with cancer cells in the tumor microenvironment. Methods: We developed a computational deconvolution method, DeClust, that stratifies patients into subtypes based on cancer cell-intrinsic signals identified by distinguishing cancer-type-specific signals from non-cancer signals in bulk tumor transcriptomic data. DeClust differs from most existing methods by directly incorporating molecular subtyping of solid tumors into the deconvolution process and outputting molecular subtype-specific tumor reference profiles for the cohort rather than individual tumor profiles. In addition, DeClust does not require reference expression profiles or signature matrices as inputs and estimates cancer-type-specific microenvironment signals from bulk tumor transcriptomic data. Results: DeClust was evaluated on both simulated data and 13 solid tumor datasets from The Cancer Genome Atlas (TCGA). DeClust performed among the best, relative to existing methods, for estimation of cellular composition. Compared to molecular subtypes reported by TCGA or other similar approaches, the subtypes generated by DeClust had higher correlations with cancer-intrinsic genomic alterations (e.g., somatic mutations and copy number variations) and lower correlations with tumor purity. While DeClust-identified subtypes were not more significantly associated with survival in general, DeClust identified a poor prognosis subtype of clear cell renal cancer, papillary renal cancer, and lung adenocarcinoma, all of which were characterized by CDKN2A deletions. As a reference profile-free deconvolution method, the tumor-type-specific stromal profiles and cancer cell-intrinsic subtypes generated by DeClust were supported by single-cell RNA sequencing data. Conclusions: DeClust is a useful tool for cancer cell-intrinsic molecular subtyping of solid tumors. DeClust subtypes, together with the tumor-type-specific stromal profiles generated by this pan-cancer study, may lead to mechanistic and clinical insights across multiple tumor types.

Chronic hepatitis B virus (HBV) infection leads to liver fibrosis, which is a major risk factor in hepatocellular carcinoma (HCC) and an independent risk factor of recurrence after HCC tumor resection. The HBV genome can be inserted into the human genome, and chronic inflammation may trigger somatic mutations. However, how HBV integration and other genomic changes contribute to the risk of tumor recurrence with regards to the different degree of liver fibrosis is not clearly understood. We sequenced mRNAs of 21 pairs of tumor and distant non-neoplastic liver tissues of HBV-HCC patients and performed comprehensive genomic analyses of our RNAseq data and public available HBV-HCC sequencing data. We developed a robust pipeline for sensitively identifying HBV integration sites based on sequencing data. Simulations showed that our method outperformed existing methods. Applying it to our data, 374 and 106 HBV host genes were identified in non-neoplastic liver and tumor tissues, respectively. When applying it to other RNA sequencing datasets, consistently more HBV integrations were identified in non-neoplastic liver than in tumor tissues. HBV host genes identified in non-neoplastic liver samples significantly overlapped with known tumor suppressor genes. More significant enrichment of tumor suppressor genes was observed among HBV host genes identified from patients with tumor recurrence, indicating the potential risk of tumor recurrence driven by HBV integration in non-neoplastic liver tissues. We also compared SNPs of each sample with SNPs in a cancer census database and inferred samples’ pathogenic SNP loads. Pathogenic SNP loads in non-neoplastic liver tissues were consistently higher than those in normal liver tissues. Additionally, HBV host genes identified in non-neoplastic liver tissues significantly overlapped with pathogenic somatic mutations, suggesting that HBV integration and somatic mutations targeting the same set of genes are important to tumorigenesis. HBV integrations and pathogenic mutations showed distinct patterns between low and high liver fibrosis patients with regards to tumor recurrence. The results suggest that HBV integrations and pathogenic SNPs in non-neoplastic tissues are important for tumorigenesis and different recurrence risk models are needed for patients with low and high degrees of liver fibrosis.

Identifying disease genes is crucial to the understanding of disease pathogenesis, and to the improvement of disease diagnosis and treatment. In recent years, many researchers have proposed approaches to prioritize candidate genes by considering the relationship of candidate genes and existing known disease genes, reflected in other data sources. In this paper, we propose an expandable framework for gene prioritization that can integrate multiple heterogeneous data sources by taking advantage of a unified graphic representation. Gene-gene relationships and gene-disease relationships are then defined based on the overall topology of each network using a diffusion kernel measure. These relationship measures are in turn normalized to derive an overall measure across all networks, which is utilized to rank all candidate genes. Based on the informativeness of available data sources with respect to each specific disease, we also propose an adaptive threshold score to select a small subset of candidate genes for further validation studies. We performed large scale cross-validation analysis on 110 disease families using three data sources. Results have shown that our approach consistently outperforms other two state of the art programs. A case study using Parkinson disease (PD) has identified four candidate genes (UBB, SEPT5, GPR37 and TH) that ranked higher than our adaptive threshold, all of which are involved in the PD pathway. In particular, a very recent study has observed a deletion of TH in a patient with PD, which supports the importance of the TH gene in PD pathogenesis. A web tool has been implemented to assist scientists in their genetic studies.