MinePath: Mining for Phenotype Differential Sub-paths in Molecular Pathways

Artificial intelligence in pancreatic cancer: diagnosis, limitations, and the future prospects-a narrative review

PURPOSE

This review aims to explore the role of AI in the application of pancreatic cancer management and make recommendations to minimize the impact of the limitations to provide further benefits from AI use in the future.

METHODS

A comprehensive review of the literature was conducted using a combination of MeSH keywords, including "Artificial intelligence", "Pancreatic cancer", "Diagnosis", and "Limitations".

RESULTS

The beneficial implications of AI in the detection of biomarkers, diagnosis, and prognosis of pancreatic cancer have been explored. In addition, current drawbacks of AI use have been divided into subcategories encompassing statistical, training, and knowledge limitations; data handling, ethical and medicolegal aspects; and clinical integration and implementation.

CONCLUSION

Artificial intelligence (AI) refers to computational machine systems that accomplish a set of given tasks by imitating human intelligence in an exponential learning pattern. AI in gastrointestinal oncology has continued to provide significant advancements in the clinical, molecular, and radiological diagnosis and intervention techniques required to improve the prognosis of many gastrointestinal cancer types, particularly pancreatic cancer.

Inference of a Dynamic Aging-related Biological Subnetwork via Network Propagation

Gene expression (GE)data capture valuable condition-specific information ("condition" can mean a biological process, disease stage, age, patient, etc.)However, GE analyses ignore physical interactions between gene products, i.e., proteins. Because proteins function by interacting with each other, and because biological networks (BNs)capture these interactions, BN analyses are promising. However, current BN data fail to capture condition-specific information. Recently, GE and BN data have been integrated using network propagation (NP)to infer condition-specific BNs. However, existing NP-based studies result in a static condition-specific subnetwork, even though cellular processes are dynamic. A dynamic process of our interest is human aging. We use prominent existing NP methods in a new task of inferring a dynamic rather than static condition-specific (aging-related)subnetwork. Then, we study evolution of network structure with age - we identify proteins whose network positions significantly change with age and predict them as new aging-related candidates. We validate the predictions via e.g., functional enrichment analyses and literature search. Dynamic network inference via NP yields higher prediction quality than the only existing method for inferring a dynamic aging-related BN, which does not use NP. Our data and code are available at https://nd.edu/~cone/dynetinf.

Transcriptional effects of fingolimod treatment on peripheral T cells in relapsing remitting multiple sclerosis patients

Aim: To investigate the transcriptional changes induced by Fingolimod (FTY) in T cells of relapsing remitting multiple sclerosis patients. Patients & methods: Transcriptomic changes after 6 months of FTY therapy were evaluated on T cells from 24 relapsing remitting multiple sclerosis patients through RNA-sequencing, followed by technical validation and pathway analysis. Results: Among differentially expressed genes, CX3CR1 and CCR7 resulted strongly up- and down-regulated, respectively. Two relevant genes were validated with quantitative PCR and we largely confirmed findings from two previous microarray-based studies with similar design. Pathway analysis pointed to an involvement of processes related to immune function and cell migration. Conclusion: Our data support the evidence that FTY induces major transcriptional changes in genes involved in immune response and cell trafficking in T lymphocytes.

PATHOME-Drug: a subpathway-based polypharmacology drug-repositioning method

MOTIVATION

Drug repositioning reveals novel indications for existing drugs and in particular, diseases with no available drugs. Diverse computational drug repositioning methods have been proposed by measuring either drug-treated gene expression signatures or the proximity of drug targets and disease proteins found in prior networks. However, these methods do not explain which signaling subparts allow potential drugs to be selected, and do not consider polypharmacology, i.e. multiple targets of a known drug, in specific subparts.

RESULTS

Here, to address the limitations, we developed a subpathway-based polypharmacology drug repositioning method, PATHOME-Drug, based on drug-associated transcriptomes. Specifically, this tool locates subparts of signaling cascading related to phenotype changes (e.g. disease status changes), and identifies existing approved drugs such that their multiple targets are enriched in the subparts. We show that our method demonstrated better performance for detecting signaling context and specific drugs/compounds, compared to WebGestalt and clusterProfiler, for both real biological and simulated datasets. We believe that our tool can successfully address the current shortage of targeted therapy agents.

AVAILABILITY AND IMPLEMENTATION

The web-service is available at http://statgen.snu.ac.kr/software/pathome. The source codes and data are available at https://github.com/labnams/pathome-drug.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent advances in artificial intelligence for pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) remains the most lethal type of cancer. The 5-year survival rate for patients with early-stage diagnosis can be as high as 20%, suggesting that early diagnosis plays a pivotal role in the prognostic improvement of PDAC cases. In the medical field, the broad availability of biomedical data has led to the advent of the "big data" era. To overcome this deadly disease, how to fully exploit big data is a new challenge in the era of precision medicine. Artificial intelligence (AI) is the ability of a machine to learn and display intelligence to solve problems. AI can help to transform big data into clinically actionable insights more efficiently, reduce inevitable errors to improve diagnostic accuracy, and make real-time predictions. AI-based omics analyses will become the next alterative approach to overcome this poor-prognostic disease by discovering biomarkers for early detection, providing molecular/genomic subtyping, offering treatment guidance, and predicting recurrence and survival. Advances in AI may therefore improve PDAC survival outcomes in the near future. The present review mainly focuses on recent advances of AI in PDAC for clinicians. We believe that breakthroughs will soon emerge to fight this deadly disease using AI-navigated precision medicine.

Transcriptomic Analysis of Peripheral Monocytes upon Fingolimod Treatment in Relapsing Remitting Multiple Sclerosis Patients

Fingolimod (FTY), a second-line oral drug approved for relapsing remitting Multiple Sclerosis (RRMS) acts in preventing lymphocyte migration outside lymph nodes; moreover, several lines of evidence suggest that it also inhibits myeloid cell activation. In this study, we investigated the transcriptional changes induced by FTY in monocytes in order to better elucidate its mechanism of action. CD14⁺ monocytes were collected from 24 RRMS patients sampled at baseline and after 6 months of treatment and RNA profiles were obtained through next-generation sequencing. We conducted pathway and sub-paths analysis, followed by centrality analysis of cell-specific interactomes on differentially expressed genes (DEGs). We investigated also the predictive role of baseline monocyte transcription profile in influencing the response to FTY therapy. We observed a marked down-regulation effect (60 down-regulated vs. 0 up-regulated genes). Most of the down-regulated DEGs resulted related with monocyte activation and migration like IL7R, CCR7 and the Wnt signaling mediators LEF1 and TCF7. The involvement of Wnt signaling was also confirmed by subpaths analyses. Furthermore, pathway and network analyses showed an involvement of processes related to immune function and cell migration. Baseline transcriptional profile of the HLA class II gene HLA-DQA1 and HLA-DPA1 were associated with evidence of disease activity after 2 years of treatment. Our data support the evidence that FTY induces major transcriptional changes in monocytes, mainly regarding genes involved in cell trafficking and immune cell activation. The baseline transcriptional levels of genes associated with antigen presenting function were associated with disease activity after 2 years of FTY treatment.

Genome-scale mechanistic modeling of signaling pathways made easy: A bioconductor/cytoscape/web server framework for the analysis of omic data

Genome-scale mechanistic models of pathways are gaining importance for genomic data interpretation because they provide a natural link between genotype measurements (transcriptomics or genomics data) and the phenotype of the cell (its functional behavior). Moreover, mechanistic models can be used to predict the potential effect of interventions, including drug inhibitions. Here, we present the implementation of a mechanistic model of cell signaling for the interpretation of transcriptomic data as an R/Bioconductor package, a Cytoscape plugin and a web tool with enhanced functionality which includes building interpretable predictors, estimation of the effect of perturbations and assessment of the effect of mutations in complex scenarios.

Deep learning models in genomics; are we there yet?

With the evolution of biotechnology and the introduction of the high throughput sequencing, researchers have the ability to produce and analyze vast amounts of genomics data. Since genomics produce big data, most of the bioinformatics algorithms are based on machine learning methodologies, and lately deep learning, to identify patterns, make predictions and model the progression or treatment of a disease. Advances in deep learning created an unprecedented momentum in biomedical informatics and have given rise to new bioinformatics and computational biology research areas. It is evident that deep learning models can provide higher accuracies in specific tasks of genomics than the state of the art methodologies. Given the growing trend on the application of deep learning architectures in genomics research, in this mini review we outline the most prominent models, we highlight possible pitfalls and discuss future directions. We foresee deep learning accelerating changes in the area of genomics, especially for multi-scale and multimodal data analysis for precision medicine.

Two gene set variation indexes as potential diagnostic tool for sepsis

Accurate diagnosis of sepsis remains challenging, new markers or combinations of markers are urgently needed. In the present study, we screened differentially expressed genes (DEGs) between sepsis and non-sepsis blood samples across three previously published gene expression data sets. Common upregulated and downregulated DEGs were ranked according to their average functional similarity. The ten genes (OLFM4, ORM1, CEP55, S100A12, S100P, LRG1, CEACAM8, MS4A4A, PLSCR1, and IL1R2) with the largest average functional similarity among the common upregulated genes and another ten genes (THEMIS, IL2RB, CD2, IL7R, CD3E, KLRB1, PVRIG, CCRR3, TGFBR3, and PLEKHA1) with the largest average functional similarity among the common downregulated genes were separately identified as the upregulated crucial gene set and the downregulated crucial gene set. Gene set variation analysis (GSVA) was used to obtain the GSVA index of each sample against the two crucial gene sets. Both the two crucial GSVA indexes may be robust markers for sepsis with high area under ROC curve. The diagnostic utility of the upregulated GSVA index was validated in another independent data set. Functional analyses revealed several sepsis-related pathways. In conclusion, we proposed two sepsis-related gene sets across multiple data sets and created two GSVA indexes with promising diagnostic value.

Patient empowerment for cancer patients through a novel ICT infrastructure

As a result of recent advances in cancer research and "precision medicine" approaches, i.e. the idea of treating each patient with the right drug at the right time, more and more cancer patients are being cured, or might have to cope with a life with cancer. For many people, cancer survival today means living with a complex and chronic condition. Surviving and living with or beyond cancer requires the long-term management of the disease, leading to a significant need for active rehabilitation of the patients. In this paper, we present a novel methodology employed in the iManageCancer project for cancer patient empowerment in which personal health systems, serious games, psychoemotional monitoring and other novel decision-support tools are combined into an integrated patient empowerment platform. We present in detail the ICT infrastructure developed and our evaluation with the involvement of cancer patients on two sites, a large-scale pilot for adults and a small-scale test for children. The evaluation showed mixed evidences on the improvement of patient empowerment, while ability to cope with cancer, including improvement in mood and resilience to cancer, increased for the participants of the adults' pilot.

A comparison of mechanistic signaling pathway activity analysis methods

Understanding the aspects of cell functionality that account for disease mechanisms or drug modes of action is a main challenge for precision medicine. Classical gene-based approaches ignore the modular nature of most human traits, whereas conventional pathway enrichment approaches produce only illustrative results of limited practical utility. Recently, a family of new methods has emerged that change the focus from the whole pathways to the definition of elementary subpathways within them that have any mechanistic significance and to the study of their activities. Thus, mechanistic pathway activity (MPA) methods constitute a new paradigm that allows recoding poorly informative genomic measurements into cell activity quantitative values and relate them to phenotypes. Here we provide a review on the MPA methods available and explain their contribution to systems medicine approaches for addressing challenges in the diagnostic and treatment of complex diseases.

Identification of Cancer Dysfunctional Subpathways by Integrating DNA Methylation, Copy Number Variation, and Gene-Expression Data

A subpathway is defined as the local region of a biological pathway with specific biological functions. With the generation of large-scale sequencing data, there are more opportunities to study the molecular mechanisms of cancer development. It is necessary to investigate the potential impact of DNA methylation, copy number variation (CNV), and gene-expression changes in the molecular states of oncogenic dysfunctional subpathways. We propose a novel method, Identification of Cancer Dysfunctional Subpathways (ICDS), by integrating multi-omics data and pathway topological information to identify dysfunctional subpathways. We first calculated gene-risk scores by integrating the three following types of data: DNA methylation, CNV, and gene expression. Second, we performed a greedy search algorithm to identify the key dysfunctional subpathways within pathways for which the discriminative scores were locally maximal. Finally, a permutation test was used to calculate the statistical significance level for these key dysfunctional subpathways. We validated the effectiveness of ICDS in identifying dysregulated subpathways using datasets from liver hepatocellular carcinoma (LIHC), head-neck squamous cell carcinoma (HNSC), cervical squamous cell carcinoma, and endocervical adenocarcinoma. We further compared ICDS with methods that performed the same subpathway identification algorithm but only considered DNA methylation, CNV, or gene expression (defined as ICDS_M, ICDS_CNV, or ICDS_G, respectively). With these analyses, we confirmed that ICDS better identified cancer-associated subpathways than the three other methods, which only considered one type of data. Our ICDS method has been implemented as a freely available R-based tool (https://cran.r-project.org/web/packages/ICDS).

Identifying and exploiting gene-pathway interactions from RNA-seq data for binary phenotype

Background

RNA sequencing (RNA-seq) technology has identified multiple differentially expressed (DE) genes associated to complex disease, however, these genes only explain a modest part of variance. Omnigenic model assumes that disease may be driven by genes with indirect relevance to disease and be propagated by functional pathways. Here, we focus on identifying the interactions between the external genes and functional pathways, referring to gene-pathway interactions (GPIs). Specifically, relying on the relationship between the garrote kernel machine (GKM) and variance component test and permutations for the empirical distributions of score statistics, we propose an efficient analysis procedure as Permutation based gEne-pAthway interaction identification in binary phenotype (PEA).

Results

Various simulations show that PEA has well-calibrated type I error rates and higher power than the traditional likelihood ratio test (LRT). In addition, we perform the gene set enrichment algorithms and PEA to identifying the GPIs from a pan-cancer data (GES68086). These GPIs and genes possibly further illustrate the potential etiology of cancers, most of which are identified and some external genes and significant pathways are consistent with previous studies.

Conclusions

PEA is an efficient tool for identifying the GPIs from RNA-seq data. It can be further extended to identify the interactions between one variable and one functional set of other omics data for binary phenotypes.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0739-7) contains supplementary material, which is available to authorized users.

minepath.org: a free interactive pathway analysis web server

Minepath

( www.minepath.org ) is a web-based platform that elaborates on, and radically extends the identification of differentially expressed sub-paths in molecular pathways. Besides the network topology, the underlying MinePath algorithmic processes exploit exact gene-gene molecular relationships (e.g. activation, inhibition) and are able to identify differentially expressed pathway parts. Each pathway is decomposed into all its constituent sub-paths, which in turn are matched with corresponding gene expression profiles. The highly ranked, and phenotype inclined sub-paths are kept. Apart from the pathway analysis algorithm, the fundamental innovation of the MinePath web-server concerns its advanced visualization and interactive capabilities. To our knowledge, this is the first pathway analysis server that introduces and offers visualization of the underlying and active pathway regulatory mechanisms instead of genes. Other features include live interaction, immediate visualization of functional sub-paths per phenotype and dynamic linked annotations for the engaged genes and molecular relations. The user can download not only the results but also the corresponding web viewer framework of the performed analysis. This feature provides the flexibility to immediately publish results without publishing source/expression data, and get all the functionality of a web based pathway analysis viewer.

MIDAS: Mining differentially activated subpaths of KEGG pathways from multi-class RNA-seq data

Pathway based analysis of high throughput transcriptome data is a widely used approach to investigate biological mechanisms. Since a pathway consists of multiple functions, the recent approach is to determine condition specific sub-pathways or subpaths. However, there are several challenges. First, few existing methods utilize explicit gene expression information from RNA-seq. More importantly, subpath activity is usually an average of statistical scores, e.g., correlations, of edges in a candidate subpath, which fails to reflect gene expression quantity information. In addition, none of existing methods can handle multiple phenotypes. To address these technical problems, we designed and implemented an algorithm, MIDAS, that determines condition specific subpaths, each of which has different activities across multiple phenotypes. MIDAS utilizes gene expression quantity information fully and the network centrality information to determine condition specific subpaths. To test performance of our tool, we used TCGA breast cancer RNA-seq gene expression profiles with five molecular subtypes. 36 differentially activate subpaths were determined. The utility of our method, MIDAS, was demonstrated in four ways. All 36 subpaths are well supported by the literature information. Subsequently, we showed that these subpaths had a good discriminant power for five cancer subtype classification and also had a prognostic power in terms of survival analysis. Finally, in a performance comparison of MIDAS to a recent subpath prediction method, PATHOME, our method identified more subpaths and much more genes that are well supported by the literature information.

AVAILABILITY

http://biohealth.snu.ac.kr/software/MIDAS/.