Identifying tightly regulated and variably expressed networks by Differential Rank Conservation (DIRAC)

Exploring casual effects and shared molecular mechanism between psoriasis and liver cancer through Mendelian randomization and comprehensive bioinformatic analyses

Psoriasis (Ps), a chronic inflammatory disease affecting approximately 2 % of the global population, has been associated with an increased risk of liver cancer in observational studies. However, their causal relationships as well as underlying shared molecular mechanisms between Ps and liver cancer remain unclear. Using bidirectional Mendelian randomization analysis, we revealed that a genetic predisposition to liver cancer increased the risk of Ps in European and East Asian populations but not the other way around. Moreover, we analyzed three transcriptomic datasets of patients with Ps and liver cancer from open-source databases. Differentially expressed genes (DEGs) and disease-specific gene co-expression module analyses revealed that cell-cycle dysregulation was the shared mechanism of Ps and liver cancer. Moreover, we identified a rank-conservative gene signature shared between these two diseases, which demonstrated significance in diagnostic and prognostic predictions. These findings provided valuable insights into the interconnections between Ps and liver cancer, which may be helpful to guide therapeutic management.

Lifespan-extending interventions induce consistent patterns of fatty acid oxidation in mouse livers

Aging manifests as progressive deteriorations in homeostasis, requiring systems-level perspectives to investigate the gradual molecular dysregulation of underlying biological processes. Here, we report systemic changes in the molecular regulation of biological processes under multiple lifespan-extending interventions. Differential Rank Conservation (DIRAC) analyses of mouse liver proteomics and transcriptomics data show that mechanistically distinct lifespan-extending interventions (acarbose, 17α-estradiol, rapamycin, and calorie restriction) generally tighten the regulation of biological modules. These tightening patterns are similar across the interventions, particularly in processes such as fatty acid oxidation, immune response, and stress response. Differences in DIRAC patterns between proteins and transcripts highlight specific modules which may be tightened via augmented cap-independent translation. Moreover, the systemic shifts in fatty acid metabolism are supported through integrated analysis of liver transcriptomics data with a mouse genome-scale metabolic model. Our findings highlight the power of systems-level approaches for identifying and characterizing the biological processes involved in aging and longevity.

An overview of metabolomic and proteomic profiling in bipolar disorder and its clinical value

INTRODUCTION

Bipolar disorder (BD) is a complex psychiatric disease characterized by alternating mood episodes. As for any other psychiatric illness, currently there is no biochemical test that is able to support diagnosis or therapeutic decisions for BD. In this context, the discovery and validation of biomarkers are interesting strategies that can be achieved through proteomics and metabolomics.

AREAS COVERED

In this descriptive review, a literature search including original articles and systematic reviews published in the last decade was performed with the objective to discuss the results of BD proteomic and metabolomic profiling analyses and indicate proteins and metabolites (or metabolic pathways) with potential clinical value.

EXPERT OPINION

A large number of proteins and metabolites have been reported as potential BD biomarkers; however, most studies do not reach biomarker validation stages. An effort from the scientific community should be directed toward the validation of biomarkers and the development of simplified bioanalytical techniques or protocols to determine them in biological samples, in order to translate proteomic and metabolomic findings into clinical routine assays.

The Coaching for Cognition in Alzheimer's (COCOA) trial: Study design

Comprehensive treatment of Alzheimer's disease (AD) requires not only pharmacologic treatment but also management of existing medical conditions and lifestyle modifications including diet, cognitive training, and exercise. We present the design and methodology for the Coaching for Cognition in Alzheimer's (COCOA) trial. AD and other dementias result from the interplay of multiple interacting dysfunctional biological systems. Monotherapies have had limited success. More interventional studies are needed to test the effectiveness of multimodal multi-domain therapies for dementia prevention and treatment. Multimodal therapies use multiple interventions to address multiple systemic causes and potentiators of cognitive decline and functional loss; they can be personalized, as different sets of etiologies and systems responsive to therapy may be present in different individuals. COCOA is designed to test the hypothesis that coached multimodal interventions beneficially alter the trajectory of cognitive decline for individuals on the spectrum of AD and related dementias (ADRD). COCOA is a two-arm prospective randomized controlled trial (RCT). COCOA collects psychometric, clinical, lifestyle, genomic, proteomic, metabolomic, and microbiome data at multiple timepoints across 2 years for each participant. These data enable systems biology analyses. One arm receives standard of care and generic healthy aging recommendations. The other arm receives standard of care and personalized data-driven remote coaching. The primary outcome measure is the Memory Performance Index (MPI), a measure of cognition. The MPI is a summary statistic of the MCI Screen (MCIS). Secondary outcome measures include the Functional Assessment Staging Test (FAST), a measure of function. COCOA began enrollment in January 2018. We hypothesize that multimodal interventions will ameliorate cognitive decline and that data-driven health coaching will increase compliance, assist in personalizing multimodal interventions, and improve outcomes for patients, particularly for those in the early stages of the AD spectrum.

Highlights

The Coaching for Cognition in Alzheimer's (COCOA) trial tests personalized multimodal lifestyle interventions for Alzheimer's disease and related dementias.Dense longitudinal molecular data will be useful for future studies.Increased use of Hill's criteria in analyses may advance knowledge generation.Remote coaching may be an effective intervention.Because lifestyle interventions are inexpensive, they may be particularly valuable in reducing global socioeconomic disparities in dementia care.

IndGOterm: a qualitative method for the identification of individually dysregulated GO terms in cancer

Individual pathway analysis can dissect heterogeneities among different cancer patients and provide efficient guidelines for individualized therapy. However, the existence of the batch effect brings extensive limitations for the application of many individual methods for pathway analysis. Previously, researchers proposed that methods based on within-sample relative expression ordering (REO) of the genes are notably insensitive to ‘batch effects’. In this article, we focus on the Gene Ontology (GO) database and propose an individual qualitative GO term analysis method (IndGOterm) based on the REO of genes. Compared with some current widely used single-sample enrichment analysis methods, such as ssGSEA and GSVA, IndGOterm has a predominance of ignoring the batch effects caused by diverse technologies. Through the survival and drug responses analysis, we found IndGOterm could capture more terms connected to cancer than other single-sample enrichment analysis methods. Furthermore, through the application of IndGOterm, we found some terms that present different dysregulation models that manifest heterogenetic in homologous patients. Collectively, these results attested that IndGOterm could capture useful information from patients and be a useful tool to reveal the intrinsic characteristic of cancer. An open-source R statistical analysis package ‘IndGOterm’ is available at https://github.com/robert19960424/IndGOterm.

Synergistic AHR Binding Pathway with EMT Effects on Serous Ovarian Tumors Recognized by Multidisciplinary Integrated Analysis

Epithelial ovarian cancers (EOCs) are fatal and obstinate among gynecological malignancies in advanced stage or relapsed status, with serous carcinomas accounting for the vast majority. Unlike EOCs, borderline ovarian tumors (BOTs), including serous BOTs, maintain a semimalignant appearance. Using gene ontology (GO)-based integrative analysis, we analyzed gene set databases of serous BOTs and serous ovarian carcinomas for dysregulated GO terms and pathways and identified multiple differentially expressed genes (DEGs) in various aspects. The SRC (SRC proto-oncogene, non-receptor tyrosine kinase) gene and dysfunctional aryl hydrocarbon receptor (AHR) binding pathway consistently influenced progression-free survival and overall survival, and immunohistochemical staining revealed elevated expression of related biomarkers (SRC, ARNT, and TBP) in serous BOT and ovarian carcinoma samples. Epithelial-mesenchymal transition (EMT) is important during tumorigenesis, and we confirmed the SNAI2 (Snail family transcriptional repressor 2, SLUG) gene showing significantly high performance by immunohistochemistry. During serous ovarian tumor formation, activated AHR in the cytoplasm could cooperate with SRC, enter cell nuclei, bind to AHR nuclear translocator (ARNT) together with TATA-Box Binding Protein (TBP), and act on DNA to initiate AHR-responsive genes to cause tumor or cancer initiation. Additionally, SNAI2 in the tumor microenvironment can facilitate EMT accompanied by tumorigenesis. Although it has not been possible to classify serous BOTs and serous ovarian carcinomas as the same EOC subtype, the key determinants of relevant DEGs (SRC, ARNT, TBP, and SNAI2) found here had a crucial role in the pathogenetic mechanism of both tumor types, implying gradual evolutionary tendencies from serous BOTs to ovarian carcinomas. In the future, targeted therapy could focus on these revealed targets together with precise detection to improve therapeutic effects and patient survival rates.

Comparing Host Module Activation Patterns and Temporal Dynamics in Infection by Influenza H1N1 Viruses

Influenza is a serious global health threat that shows varying pathogenicity among different virus strains. Understanding similarities and differences among activated functional pathways in the host responses can help elucidate therapeutic targets responsible for pathogenesis. To compare the types and timing of functional modules activated in host cells by four influenza viruses of varying pathogenicity, we developed a new DYNAmic MOdule (DYNAMO) method that addresses the need to compare functional module utilization over time. This integrative approach overlays whole genome time series expression data onto an immune-specific functional network, and extracts conserved modules exhibiting either different temporal patterns or overall transcriptional activity. We identified a common core response to influenza virus infection that is temporally shifted for different viruses. We also identified differentially regulated functional modules that reveal unique elements of responses to different virus strains. Our work highlights the usefulness of combining time series gene expression data with a functional interaction map to capture temporal dynamics of the same cellular pathways under different conditions. Our results help elucidate conservation of the immune response both globally and at a granular level, and provide mechanistic insight into the differences in the host response to infection by influenza strains of varying pathogenicity.

Dysregulated Immunological Functionome and Dysfunctional Metabolic Pathway Recognized for the Pathogenesis of Borderline Ovarian Tumors by Integrative Polygenic Analytics

The pathogenesis and molecular mechanisms of ovarian low malignant potential (LMP) tumors or borderline ovarian tumors (BOTs) have not been fully elucidated to date. Surgery remains the cornerstone of treatment for this disease, and diagnosis is mainly made by histopathology to date. However, there is no integrated analysis investigating the tumorigenesis of BOTs with open experimental data. Therefore, we first utilized a functionome-based speculative model from the aggregated obtainable datasets to explore the expression profiling data among all BOTs and two major subtypes of BOTs, serous BOTs (SBOTs) and mucinous BOTs (MBOTs), by analyzing the functional regularity patterns and clustering the separate gene sets. We next prospected and assembled the association between these targeted biomolecular functions and their related genes. Our research found that BOTs can be accurately recognized by gene expression profiles by means of integrative polygenic analytics among all BOTs, SBOTs, and MBOTs; the results exhibited the top 41 common dysregulated biomolecular functions, which were sorted into four major categories: immune and inflammatory response-related functions, cell membrane- and transporter-related functions, cell cycle- and signaling-related functions, and cell metabolism-related functions, which were the key elements involved in its pathogenesis. In contrast to previous research, we identified 19 representative genes from the above classified categories (IL6, CCR2 for immune and inflammatory response-related functions; IFNG, ATP1B1, GAS6, and PSEN1 for cell membrane- and transporter-related functions; CTNNB1, GATA3, and IL1B for cell cycle- and signaling-related functions; and AKT1, SIRT1, IL4, PDGFB, MAPK3, SRC, TWIST1, TGFB1, ADIPOQ, and PPARGC1A for cell metabolism-related functions) that were relevant in the cause and development of BOTs. We also noticed that a dysfunctional pathway of galactose catabolism had taken place among all BOTs, SBOTs, and MBOTs from the analyzed gene set databases of canonical pathways. With the help of immunostaining, we verified significantly higher performance of interleukin 6 (IL6) and galactose-1-phosphate uridylyltransferase (GALT) among BOTs than the controls. In conclusion, a bioinformatic platform of gene-set integrative molecular functionomes and biophysiological pathways was constructed in this study to interpret the complicated pathogenic pathways of BOTs, and these important findings demonstrated the dysregulated immunological functionome and dysfunctional metabolic pathway as potential roles during the tumorigenesis of BOTs and may be helpful for the diagnosis and therapy of BOTs in the future.

Metabolic-Pathway-Based Subtyping of Triple-Negative Breast Cancer Reveals Potential Therapeutic Targets

Triple-negative breast cancer (TNBC) remains an unmet medical challenge. We investigated metabolic dysregulation in TNBCs by using our multi-omics database (n = 465, the largest to date). TNBC samples were classified into three heterogeneous metabolic-pathway-based subtypes (MPSs) with distinct metabolic features: MPS1, the lipogenic subtype with upregulated lipid metabolism; MPS2, the glycolytic subtype with upregulated carbohydrate and nucleotide metabolism; and MPS3, the mixed subtype with partial pathway dysregulation. These subtypes were validated by metabolomic profiling of 72 samples. These three subtypes had distinct prognoses, molecular subtype distributions, and genomic alterations. Moreover, MPS1 TNBCs were more sensitive to metabolic inhibitors targeting fatty acid synthesis, whereas MPS2 TNBCs showed higher sensitivity to inhibitors targeting glycolysis. Importantly, inhibition of lactate dehydrogenase could enhance tumor response to anti-PD-1 immunotherapy in MPS2 TNBCs. Collectively, our analysis demonstrated the metabolic heterogeneity of TNBCs and enabled the development of personalized therapies targeting unique tumor metabolic profiles.

The Potential Role of Complement System in the Progression of Ovarian Clear Cell Carcinoma Inferred from the Gene Ontology-Based Immunofunctionome Analysis

Ovarian clear cell carcinoma (OCCC) is the second most common epithelial ovarian carcinoma (EOC). It is refractory to chemotherapy with a worse prognosis after the preliminary optimal debulking operation, such that the treatment of OCCC remains a challenge. OCCC is believed to evolve from endometriosis, a chronic immune/inflammation-related disease, so that immunotherapy may be a potential alternative treatment. Here, gene set-based analysis was used to investigate the immunofunctionomes of OCCC in early and advanced stages. Quantified biological functions defined by 5917 Gene Ontology (GO) terms downloaded from the Gene Expression Omnibus (GEO) database were used. DNA microarray gene expression profiles were used to convert 85 OCCCs and 136 normal controls into to the functionome. Relevant offspring were as extracted and the immunofunctionomes were rebuilt at different stages by machine learning. Several dysregulated pathogenic functions were found to coexist in the immunopathogenesis of early and advanced OCCC, wherein the complement-activation-alternative-pathway may be the headmost dysfunctional immunological pathway in duality for carcinogenesis at all OCCC stages. Several immunological genes involved in the complement system had dual influences on patients' survival, and immunohistochemistrical analysis implied the higher expression of C3a receptor (C3aR) and C5a receptor (C5aR) levels in OCCC than in controls.

Improvement in the Current Therapies for Hepatocellular Carcinoma Using a Systems Medicine Approach

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related death primarily due to the lack of effective targeted therapies. Despite the distinct morphological and phenotypic patterns of HCC, treatment strategies are restricted to relatively homogeneous therapies, including multitargeted tyrosine kinase inhibitors and immune checkpoint inhibitors. Therefore, more effective therapy options are needed to target dysregulated metabolic and molecular pathways in HCC. Integrative genomic profiling of HCC patients provides insight into the most frequently mutated genes and molecular targets, including telomerase reverse transcriptase, the TP53 gene, and the Wnt/β-catenin signaling pathway oncogene (CTNNB1). Moreover, emerging techniques, such as genome-scale metabolic models may elucidate the underlying cancer-specific metabolism, which allows for the discovery of potential drug targets and identification of biomarkers. De novo lipogenesis has been revealed as consistently upregulated since it is required for cell proliferation in all HCC patients. The metabolic network-driven stratification of HCC patients in terms of redox responses, utilization of metabolites, and subtype-specific pathways may have clinical implications to drive the development of personalized medicine. In this review, the current and emerging therapeutic targets in light of molecular approaches and metabolic network-based strategies are summarized, prompting effective treatment of HCC patients.

The recent progress and therapy in endometriosis-associated ovarian cancer

Endometriosis-associated ovarian cancers (EAOCs) including endometrioid and clear cell ovarian carcinoma are subgroups of epithelial ovarian carcinomas (EOCs), which is generally acknowledged as the most lethal gynecological malignancy. Endometriosis (ES), a common clinical disease among women, presents with clinical symptoms of pelvic pain, infertility, or adnexal masses with the formation of endometrioma. It has long been considered to be a potential risk factor for developing EOCs, mainly of endometrioid and clear cell subtypes. Here, we compiled data from previous researches on deregulated molecular functions among ES and EOCs using gene set-based integrative analysis to decipher molecular and genetic relationships between ovarian ES and EOCs, especially EAOCs. We conclude that epidermal growth factor receptor (ERBB) and Phosphoinositide 3-kinases (PI3K)-related pathways are important in the carcinogenesis of type I EOCs, including clear cell, endometrioid, and mucinous ovarian carcinoma. Dysfunctional molecular pathways, such as deregulated oxidoreductase activity, metabolism, hormone activity, inflammatory response, innate immune response, and cell-cell signaling, played key roles in the malignant transformation of EAOCs. Nine genes related to inflammasome complex and inflammasome-related pathway were identified, indicating the importance of inflammation/immunity in EAOC transformation. We also collect progressive treatments of EAOC focused on targeted therapies and immunotherapy so far. This summarized information can contribute toward effective detection and treatment of EAOCs in the future.

Differential Variation Analysis Enables Detection of Tumor Heterogeneity Using Single-Cell RNA-Sequencing Data

Tumor heterogeneity provides a complex challenge to cancer treatment and is a critical component of therapeutic response, disease recurrence, and patient survival. Single-cell RNA-sequencing (scRNA-seq) technologies have revealed the prevalence of intratumor and intertumor heterogeneity. Computational techniques are essential to quantify the differences in variation of these profiles between distinct cell types, tumor subtypes, and patients to fully characterize intratumor and intertumor molecular heterogeneity. In this study, we adapted our algorithm for pathway dysregulation, Expression Variation Analysis (EVA), to perform multivariate statistical analyses of differential variation of expression in gene sets for scRNA-seq. EVA has high sensitivity and specificity to detect pathways with true differential heterogeneity in simulated data. EVA was applied to several public domain scRNA-seq tumor datasets to quantify the landscape of tumor heterogeneity in several key applications in cancer genomics such as immunogenicity, metastasis, and cancer subtypes. Immune pathway heterogeneity of hematopoietic cell populations in breast tumors corresponded to the amount of diversity present in the T-cell repertoire of each individual. Cells from head and neck squamous cell carcinoma (HNSCC) primary tumors had significantly more heterogeneity across pathways than cells from metastases, consistent with a model of clonal outgrowth. Moreover, there were dramatic differences in pathway dysregulation across HNSCC basal primary tumors. Within the basal primary tumors, there was increased immune dysregulation in individuals with a high proportion of fibroblasts present in the tumor microenvironment. These results demonstrate the broad utility of EVA to quantify intertumor and intratumor heterogeneity from scRNA-seq data without reliance on low-dimensional visualization. SIGNIFICANCE: This study presents a robust statistical algorithm for evaluating gene expression heterogeneity within pathways or gene sets in single-cell RNA-seq data.

Dysregulated immunological and metabolic functions discovered by a polygenic integrative analysis for PCOS

RESEARCH QUESTION

Polycystic ovary syndrome (PCOS) is a complex disease and its pathophysiology is still unclear. This polygenic study may provide some clues.

DESIGN

A polygenic, functionome-based study with the ovarian gene expression profiles downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database, including 48 PCOS and 181 normal control samples. These profiles were converted to the gene set regularity (GSR) indices, which were computed by the modified differential rank conversion algorithm and were defined by the gene ontology terms.

RESULTS

Machine learning could accurately recognize the patterns of functional regularities between PCOS and normal controls. The significantly aberrant functions in PCOS included transporter activity, catalytic activity, the receptor signalling pathway via signal transducer and activator of transcription (STAT), the cellular metabolic process, and immune response.

CONCLUSION

This study provided a comprehensive view of the dysregulated functions and information for further studies on the management of PCOS.

Splice Expression Variation Analysis (SEVA) for inter-tumor heterogeneity of gene isoform usage in cancer

Motivation

Current bioinformatics methods to detect changes in gene isoform usage in distinct phenotypes compare the relative expected isoform usage in phenotypes. These statistics model differences in isoform usage in normal tissues, which have stable regulation of gene splicing. Pathological conditions, such as cancer, can have broken regulation of splicing that increases the heterogeneity of the expression of splice variants. Inferring events with such differential heterogeneity in gene isoform usage requires new statistical approaches.

Results

We introduce Splice Expression Variability Analysis (SEVA) to model increased heterogeneity of splice variant usage between conditions (e.g. tumor and normal samples). SEVA uses a rank-based multivariate statistic that compares the variability of junction expression profiles within one condition to the variability within another. Simulated data show that SEVA is unique in modeling heterogeneity of gene isoform usage, and benchmark SEVA's performance against EBSeq, DiffSplice and rMATS that model differential isoform usage instead of heterogeneity. We confirm the accuracy of SEVA in identifying known splice variants in head and neck cancer and perform cross-study validation of novel splice variants. A novel comparison of splice variant heterogeneity between subtypes of head and neck cancer demonstrated unanticipated similarity between the heterogeneity of gene isoform usage in HPV-positive and HPV-negative subtypes and anticipated increased heterogeneity among HPV-negative samples with mutations in genes that regulate the splice variant machinery. These results show that SEVA accurately models differential heterogeneity of gene isoform usage from RNA-seq data.

Availability and implementation

SEVA is implemented in the R/Bioconductor package GSReg.

Contact

bahman@jhu.edu or favorov@sensi.org or ejfertig@jhmi.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comparative transcriptomic analysis of human placentae at term and preterm delivery

Preterm birth affects 1 out of every 10 infants in the United States, resulting in substantial neonatal morbidity and mortality. Currently, there are few predictive markers and few treatment options to prevent preterm birth. A healthy, functioning placenta is essential to positive pregnancy outcomes. Previous studies have suggested that placental pathology may play a role in preterm birth etiology. Therefore, we tested the hypothesis that preterm placentae may exhibit unique transcriptomic signatures compared to term samples reflective of their abnormal biology leading to this adverse outcome. We aggregated publicly available placental villous microarray data to generate a preterm and term sample dataset (n = 133, 55 preterm placentae and 78 normal term placentae). We identified differentially expressed genes using the linear regression for microarray (LIMMA) package and identified perturbations in known biological networks using Differential Rank Conservation (DIRAC). We identified 129 significantly differentially expressed genes between term and preterm placenta with 96 genes upregulated and 33 genes downregulated (P-value <0.05). Significant changes in gene expression in molecular networks related to Tumor Protein 53 and phosphatidylinositol signaling were identified using DIRAC. We have aggregated a uniformly normalized transcriptomic dataset and have identified novel and established genes and pathways associated with developmental regulation of the placenta and potential preterm birth pathology. These analyses provide a community resource to integrate with other high-dimensional datasets for additional insights in normal placental development and its disruption.

Discovery of therapeutic agents for prostate cancer using genome-scale metabolic modeling and drug repositioning

BACKGROUND

Genome-scale metabolic models (GEMs) offer insights into cancer metabolism and have been used to identify potential biomarkers and drug targets. Drug repositioning is a time- and cost-effective method of drug discovery that can be applied together with GEMs for effective cancer treatment.

METHODS

In this study, we reconstruct a prostate cancer (PRAD)-specific GEM for exploring prostate cancer metabolism and also repurposing new therapeutic agents that can be used in development of effective cancer treatment. We integrate global gene expression profiling of cell lines with >1000 different drugs through the use of prostate cancer GEM and predict possible drug-gene interactions.

FINDINGS

We identify the key reactions with altered fluxes based on the gene expression changes and predict the potential drug effect in prostate cancer treatment. We find that sulfamethoxypyridazine, azlocillin, hydroflumethiazide, and ifenprodil can be repurposed for the treatment of prostate cancer based on an in silico cell viability assay. Finally, we validate the effect of ifenprodil using an in vitro cell assay and show its inhibitory effect on a prostate cancer cell line.

INTERPRETATION

Our approach demonstate how GEMs can be used to predict therapeutic agents for cancer treatment based on drug repositioning. Besides, it paved a way and shed a light on the applicability of computational models to real-world biomedical or pharmaceutical problems.

Precision Medicine in Pancreatic Disease-Knowledge Gaps and Research Opportunities: Summary of a National Institute of Diabetes and Digestive and Kidney Diseases Workshop

A workshop on research gaps and opportunities for Precision Medicine in Pancreatic Disease was sponsored by the National Institute of Diabetes and Digestive Kidney Diseases on July 24, 2019, in Pittsburgh. The workshop included an overview lecture on precision medicine in cancer and 4 sessions: (1) general considerations for the application of bioinformatics and artificial intelligence; (2) omics, the combination of risk factors and biomarkers; (3) precision imaging; and (4) gaps, barriers, and needs to move from precision to personalized medicine for pancreatic disease. Current precision medicine approaches and tools were reviewed, and participants identified knowledge gaps and research needs that hinder bringing precision medicine to pancreatic diseases. Most critical were (a) multicenter efforts to collect large-scale patient data sets from multiple data streams in the context of environmental and social factors; (b) new information systems that can collect, annotate, and quantify data to inform disease mechanisms; (c) novel prospective clinical trial designs to test and improve therapies; and (d) a framework for measuring and assessing the value of proposed approaches to the health care system. With these advances, precision medicine can identify patients early in the course of their pancreatic disease and prevent progression to chronic or fatal illness.

Key Immunological Functions Involved in the Progression of Epithelial Ovarian Serous Carcinoma Discovered by the Gene Ontology-Based Immunofunctionome Analysis

Serous carcinoma (SC) is the most common and lethal subtype of epithelial ovarian carcinoma; immunotherapy is a potential treatment for SC, however, the global immunological functions of SC as well as their change during the progression of SC have not been investigated in detail till now. We conducted a genome-wide integrative analysis to investigate the immunofunctionomes of SC at four tumor stages by quantifying the immunological functions defined by the Gene Ontology gene sets. DNA microarray gene expression profiles of 1100 SCs and 136 normal ovarian tissue controls were downloaded from the Gene Expression Omnibus database and converted to the functionome. Then the immunofunctionomes were reconstructed by extracting the offspring from the functionome for the four SC staging groups. The key immunological functions extracted from immunofunctionomes with a series of filters revealed that the immunopathy of SC consisted of a group of deregulated functions with the core members including B cell activation and differentiation, regulation of leukocyte chemotaxis/cellular extravasation, antigen receptor mediated signaling pathway, T helper mediated immunity and macrophage activation; and the auxiliary elements included leukocyte mediated immunity, regulation of inflammatory response, T cell differentiation, mononuclear cell migration, megakaryocyte differentiation, complement activation and cytokine production. These deregulated immunological functions reveal the candidates to target in the immunotherapy.

TLR5 decoy receptor as a novel anti-amyloid therapeutic for Alzheimer's disease

There is considerable interest in harnessing innate immunity to treat Alzheimer's disease (AD). Here, we explore whether a decoy receptor strategy using the ectodomain of select TLRs has therapeutic potential in AD. AAV-mediated expression of human TLR5 ectodomain (sTLR5) alone or fused to human IgG4 Fc (sTLR5Fc) results in robust attenuation of amyloid β (Aβ) accumulation in a mouse model of Alzheimer-type Aβ pathology. sTLR5Fc binds to oligomeric and fibrillar Aβ with high affinity, forms complexes with Aβ, and blocks Aβ toxicity. Oligomeric and fibrillar Aβ modulates flagellin-mediated activation of human TLR5 but does not, by itself, activate TLR5 signaling. Genetic analysis shows that rare protein coding variants in human TLR5 may be associated with a reduced risk of AD. Further, transcriptome analysis shows altered TLR gene expression in human AD. Collectively, our data suggest that TLR5 decoy receptor-based biologics represent a novel and safe Aβ-selective class of biotherapy in AD.

Digitizing omics profiles by divergence from a baseline

Technological advances enable increasingly comprehensive profiling of the molecular landscapes of cells, and these data can inform the personalized treatment of complex diseases. Two major obstacles are the complexity of these data and the high degree of person-to-person heterogeneity. We develop a highly simplified, personalized data representation by comparing the profile of an individual to the range of landscapes in a baseline population, thereby mimicking basic clinical diagnostic testing for departures of selected variables from normal levels. Moreover, our method can be applied to any data modality and at any level of granularity, from single features to any subset of features treated as a single entity, for example the gene expression levels in a pathway. Experiments involve both healthy human tissues and various cancer subtypes.

Data collected from omics technologies have revealed pervasive heterogeneity and stochasticity of molecular states within and between phenotypes. A prominent example of such heterogeneity occurs between genome-wide mRNA, microRNA, and methylation profiles from one individual tumor to another, even within a cancer subtype. However, current methods in bioinformatics, such as detecting differentially expressed genes or CpG sites, are population-based and therefore do not effectively model intersample diversity. Here we introduce a unified theory to quantify sample-level heterogeneity that is applicable to a single omics profile. Specifically, we simplify an omics profile to a digital representation based on the omics profiles from a set of samples from a reference or baseline population (e.g., normal tissues). The state of any subprofile (e.g., expression vector for a subset of genes) is said to be “divergent” if it lies outside the estimated support of the baseline distribution and is consequently interpreted as “dysregulated” relative to that baseline. We focus on two cases: single features (e.g., individual genes) and distinguished subsets (e.g., regulatory pathways). Notably, since the divergence analysis is at the individual sample level, dysregulation can be analyzed probabilistically; for example, one can estimate the probability that a gene or pathway is divergent in some population. Finally, the reduction in complexity facilitates a more “personalized” and biologically interpretable analysis of variation, as illustrated by experiments involving tissue characterization, disease detection and progression, and disease–pathway associations.

Comparative analysis of gene expression in maternal peripheral blood and monocytes during spontaneous preterm labor

BACKGROUND

Preterm birth is the leading cause of newborn death worldwide, and is associated with significant cognitive and physiological challenges in later life. There is a pressing need to define the mechanisms that initiate spontaneous preterm labor, and for development of novel clinical biomarkers to identify high-risk pregnancies. Most preterm birth studies utilize fetal tissues, and there is limited understanding of the transcriptional changes that occur in mothers undergoing spontaneous preterm labor. Earlier work revealed that a specific population of maternal peripheral leukocytes (macrophages/monocytes) play an active role in the initiation of labor. Thus, we hypothesized that there are dynamic gene expression changes in maternal blood leukocytes during preterm labor.

OBJECTIVE

Using next-generation sequencing we aim to characterize the transcriptome in whole blood leukocytes and peripheral monocytes of women undergoing spontaneous preterm labor compared to healthy pregnant women who subsequently delivered at full term.

STUDY DESIGN

RNA sequencing was performed in both whole blood and peripheral monocytes from women who underwent preterm labor (24-34 weeks of gestation, N = 20) matched for gestational age to healthy pregnant controls (N = 30). All participants were a part of the Ontario Birth Study cohort (Toronto, Ontario, Canada).

RESULTS

We identified significant differences in expression of 262 genes in peripheral monocytes and 184 genes in whole blood of women who were in active spontaneous preterm labor compared to pregnant women of the same gestational age not undergoing labor, with 43 of these genes differentially expressed in both whole blood and peripheral monocytes. ADAMTS2 expression was significantly increased in women actively undergoing spontaneous preterm labor, which we validated through digital droplet reverse transcriptase polymerase chain reaction. Intriguingly, we have also identified a number of gene sets including signaling by stem cell factor-KIT, nucleotide metabolism, and trans-Golgi network vesicle budding, which exhibited changes in relative gene expression that was predictive of preterm labor status in both maternal whole blood and peripheral monocytes.

CONCLUSION

This study is the first to investigate changes in both whole blood leukocytes and peripheral monocytes of women actively undergoing spontaneous preterm labor through robust transcript measurements from RNA sequencing. Our unique study design overcame confounding based on gestational age by collecting blood samples from women matched by gestational age, allowing us to study transcriptomic changes directly related to the active preterm parturition. We performed RNA profiling using whole genome sequencing, which is highly sensitive and allowed us to identify subtle changes in specific genes. ADAMTS2 expression emerged as a marker of prematurity within peripheral blood leukocytes, an accessible tissue that plays a functional role in signaling during the onset of labor. We identified changes in relative gene expression in a number of gene sets related to signaling in monocytes and whole blood of women undergoing spontaneous preterm labor compared to controls. These genes and pathways may help identify potential targets for the development of novel drugs for preterm birth prevention.

Integrating the dysregulated inflammasome-based molecular functionome in the malignant transformation of endometriosis-associated ovarian carcinoma

The coexistence of endometriosis (ES) with ovarian clear cell carcinoma (CCC) or endometrioid carcinoma (EC) suggested that malignant transformation of ES leads to endometriosis associated ovarian carcinoma (EAOC). However, there is still lack of an integrating data analysis of the accumulated experimental data to provide the evidence supporting the hypothesis of EAOC transformation. Herein we used a function-based analytic model with the publicly available microarray datasets to investigate the expression profiling between ES, CCC, and EC. We analyzed the functional regularity pattern of the three type of samples and hierarchically clustered the gene sets to identify key mechanisms regulating the malignant transformation of EAOC. We identified a list of 18 genes (NLRP3, AIM2, PYCARD, NAIP, Caspase-4, Caspase-7, Caspase-8, TLR1, TLR7, TOLLIP, NFKBIA, TNF, TNFAIP3, INFGR2, P2RX7, IL-1B, IL1RL1, IL-18) closely related to inflammasome complex, indicating an important role of inflammation/immunity in EAOC transformation. We next explore the association between these target genes and patient survival using Gene Expression Omnibus (GEO), and found significant correlation between the expression levels of the target genes and the progression-free survival. Interestingly, high expression levels of AIM2 and NLRP3, initiating proteins of inflammasomes, were significantly correlated with poor progression-free survival. Immunohistochemistry staining confirmed a correlation between high AIM2 and high Ki-67 in clinical EAOC samples, supporting its role in disease progression. Collectively, we established a bioinformatic platform of gene-set integrative molecular functionome to dissect the pathogenic pathways of EAOC, and demonstrated a key role of dysregulated inflammasome in modulating the malignant transformation of EAOC.

Gene set-based analysis of mucinous ovarian carcinoma

OBJECTIVE

Mucinous ovarian carcinoma (MOC) is an uncommon subtype of epithelial ovarian cancers, and the pathogenesis is still poorly understood because of its rarity. We conducted a gene set-based analysis to investigate the pathogenesis of MOC by integrating microarray gene expression datasets based on the regularity of functions defined by gene ontology or canonical pathway databases.

MATERIALS AND METHODS

Forty-five pairs of MOC and normal ovarian tissue sample gene expression profiles were downloaded from the National Center for Biotechnology Information Gene Expression Omnibus database. The gene expression profiles were converted to the gene set regularity indexes by measuring the change of gene expression ordering in a gene set. Then the pathogenesis of MOC was investigated with the differences of function regularity with the gene set regularity indexes between the MOC and normal control samples.

RESULTS

The informativeness of the gene set regularity indexes was sufficient for machine learning to accurately recognize and classify the functional regulation patterns with an accuracy of 99.44%. The statistical analysis revealed that the GTPase regulators and receptor tyrosine kinase erbB-2 (ERBB2) were the most important aberrations; the exploratory factor analysis revealed phosphoinositide 3-kinase-activating kinase, G-protein coupled receptor pathway, oxidoreductase activity, immune response, peptidase activity, regulation of translation, and transport and channel activity were also involved in the pathogenesis of MOC.

CONCLUSION

Investigating the pathogenesis of MOC with the functionome provided a comprehensive view of the deregulated functions of this disease. In addition to GTPase regulators and ERBB2, a plenty of deregulated functions such as phosphoinositide 3-kinase, G-protein coupled receptor pathway, and immune response also participated in the interaction network of MOC pathogenesis.

Discovering the Deregulated Molecular Functions Involved in Malignant Transformation of Endometriosis to Endometriosis-Associated Ovarian Carcinoma Using a Data-Driven, Function-Based Analysis

The clinical characteristics of clear cell carcinoma (CCC) and endometrioid carcinoma EC) are concomitant with endometriosis (ES), which leads to the postulation of malignant transformation of ES to endometriosis-associated ovarian carcinoma (EAOC). Different deregulated functional areas were proposed accounting for the pathogenesis of EAOC transformation, and there is still a lack of a data-driven analysis with the accumulated experimental data in publicly-available databases to incorporate the deregulated functions involved in the malignant transformation of EOAC. We used the microarray gene expression datasets of ES, CCC and EC downloaded from the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) database. Then, we investigated the pathogenesis of EAOC by a data-driven, function-based analytic model with the quantified molecular functions defined by 1454 Gene Ontology (GO) term gene sets. This model converts the gene expression profiles to the functionome consisting of 1454 quantified GO functions, and then, the key functions involving the malignant transformation of EOAC can be extracted by a series of filters. Our results demonstrate that the deregulated oxidoreductase activity, metabolism, hormone activity, inflammatory response, innate immune response and cell-cell signaling play the key roles in the malignant transformation of EAOC. These results provide the evidence supporting the specific molecular pathways involved in the malignant transformation of EAOC.

GRAPE: a pathway template method to characterize tissue-specific functionality from gene expression profiles

BACKGROUND

Personalizing treatment regimes based on gene expression profiles of individual tumors will facilitate management of cancer. Although many methods have been developed to identify pathways perturbed in tumors, the results are often not generalizable across independent datasets due to the presence of platform/batch effects. There is a need to develop methods that are robust to platform/batch effects and able to identify perturbed pathways in individual samples.

RESULTS

We present Gene-Ranking Analysis of Pathway Expression (GRAPE) as a novel method to identify abnormal pathways in individual samples that is robust to platform/batch effects in gene expression profiles generated by multiple platforms. GRAPE first defines a template consisting of an ordered set of pathway genes to characterize the normative state of a pathway based on the relative rankings of gene expression levels across a set of reference samples. This template can be used to assess whether a sample conforms to or deviates from the typical behavior of the reference samples for this pathway. We demonstrate that GRAPE performs well versus existing methods in classifying tissue types within a single dataset, and that GRAPE achieves superior robustness and generalizability across different datasets. A powerful feature of GRAPE is the ability to represent individual gene expression profiles as a vector of pathways scores. We present applications to the analyses of breast cancer subtypes and different colonic diseases. We perform survival analysis of several TCGA subtypes and find that GRAPE pathway scores perform well in comparison to other methods.

CONCLUSIONS

GRAPE templates offer a novel approach for summarizing the behavior of gene-sets across a collection of gene expression profiles. These templates offer superior robustness across distinct experimental batches compared to existing methods. GRAPE pathway scores enable identification of abnormal gene-set behavior in individual samples using a non-competitive approach that is fundamentally distinct from popular enrichment-based methods. GRAPE may be an appropriate tool for researchers seeking to identify individual samples displaying abnormal gene-set behavior as well as to explore differences in the consensus gene-set behavior of groups of samples. GRAPE is available in R for download at https://CRAN.R-project.org/package=GRAPE .

Identifying disease-associated pathways in one-phenotype data based on reversal gene expression orderings

Due to the invasiveness nature of tissue biopsy, it is common that investigators cannot collect sufficient normal controls for comparison with diseased samples. We developed a pathway enrichment tool, DRFunc, to detect significantly disease-disrupted pathways by incorporating normal controls from other experiments. The method was validated using both microarray and RNA-seq expression data for different cancers. The high concordant differentially ranked (DR) gene pairs were identified between cases and controls from different independent datasets. The DR gene pairs were used in the DRFunc algorithm to detect significantly disrupted pathways in one-phenotype expression data by combing controls from other studies. The DRFunc algorithm was exemplified by the detection of significant pathways in glioblastoma samples. The algorithm can also be used to detect altered pathways in the datasets with weak expression signals, as shown by the analysis on the expression data of chemotherapy-treated breast cancer samples.

Bioinformatics Meets Biomedicine: OncoFinder, a Quantitative Approach for Interrogating Molecular Pathways Using Gene Expression Data

We propose a biomathematical approach termed OncoFinder (OF) that enables performing both quantitative and qualitative analyses of the intracellular molecular pathway activation. OF utilizes an algorithm that distinguishes the activator/repressor role of every gene product in a pathway. This method is applicable for the analysis of any physiological, stress, malignancy, and other conditions at the molecular level. OF showed a strong potential to neutralize background-caused differences between experimental gene expression data obtained using NGS, microarray and modern proteomics techniques. Importantly, in most cases, pathway activation signatures were better markers of cancer progression compared to the individual gene products. OF also enables correlating pathway activation with the success of anticancer therapy for individual patients. We further expanded this approach to analyze impact of micro RNAs (miRs) on the regulation of cellular interactome. Many alternative sources provide information about miRs and their targets. However, instruments elucidating higher level impact of the established total miR profiles are still largely missing. A variant of OncoFinder termed MiRImpact enables linking miR expression data with its estimated outcome on the regulation of molecular processes, such as signaling, metabolic, cytoskeleton, and DNA repair pathways. MiRImpact was used to establish cancer-specific and cytomegaloviral infection-linked interactomic signatures for hundreds of molecular pathways. Interestingly, the impact of miRs appeared orthogonal to pathway regulation at the mRNA level, which stresses the importance of combining all available levels of gene regulation to build a more objective molecular model of cell.

MiRImpact, a new bioinformatic method using complete microRNA expression profiles to assess their overall influence on the activity of intracellular molecular pathways

MicroRNAs (miRs) are short noncoding RNA molecules that regulate expression of target mRNAs. Many published sources provide information about miRs and their targets. However, bioinformatic tools elucidating higher level impact of the established total miR profiles, are still largely missing. Recently, we developed a method termed OncoFinder enabling quantification of the activities of intracellular molecular pathways basing on gene expression data. Here we propose a new technique, MiRImpact, which enables to link miR expression data with its estimated outcome on the regulation of molecular pathways, like signaling, metabolic, cytoskeleton rearrangement, and DNA repair pathways. MiRImpact uses OncoFinder rationale for pathway activity calculations, with the major distinctions that (i) it deals with the concentrations of miRs - known regulators of gene products participating in molecular pathways, and (ii) miRs are considered as negative regulators of target molecules, if other is not specified. MiRImpact operates with 2 types of databases: for molecular targets of miRs and for gene products participating in molecular pathways. We applied MiRImpact to compare regulation of human bladder cancer-specific signaling pathways at the levels of mRNA and miR expression. We took 2 most complete alternative databases of experimentally validated miR targets – miRTarBase and DianaTarBase, and an OncoFinder database featuring 2725 gene products and 271 signaling pathways. We showed that the impact of miRs is orthogonal to pathway regulation at the mRNA level, which stresses the importance of studying posttranscriptional regulation of gene expression. We also report characteristic set of miR and mRNA regulation features linked with bladder cancer.

Molecular pathway activation features linked with transition from normal skin to primary and metastatic melanomas in human

Melanoma is the most aggressive and dangerous type of skin cancer, but its molecular mechanisms remain largely unclear. For transcriptomic data of 478 primary and metastatic melanoma, nevi and normal skin samples, we performed high-throughput analysis of intracellular molecular networks including 592 signaling and metabolic pathways. We showed that at the molecular pathway level, the formation of nevi largely resembles transition from normal skin to primary melanoma. Using a combination of bioinformatic machine learning algorithms, we identified 44 characteristic signaling and metabolic pathways connected with the formation of nevi, development of primary melanoma, and its metastases. We created a model describing formation and progression of melanoma at the level of molecular pathway activation. We discovered six novel associations between activation of metabolic molecular pathways and progression of melanoma: for allopregnanolone biosynthesis, L-carnitine biosynthesis, zymosterol biosynthesis (inhibited in melanoma), fructose 2, 6-bisphosphate synthesis and dephosphorylation, resolvin D biosynthesis (activated in melanoma), D-myo-inositol hexakisphosphate biosynthesis (activated in primary, inhibited in metastatic melanoma). Finally, we discovered fourteen tightly coordinated functional clusters of molecular pathways. This study helps to decode molecular mechanisms underlying the development of melanoma.

Genome scale metabolic modeling of cancer

Cancer cells reprogram metabolism to support rapid proliferation and survival. Energy metabolism is particularly important for growth and genes encoding enzymes involved in energy metabolism are frequently altered in cancer cells. A genome scale metabolic model (GEM) is a mathematical formalization of metabolism which allows simulation and hypotheses testing of metabolic strategies. It has successfully been applied to many microorganisms and is now used to study cancer metabolism. Generic models of human metabolism have been reconstructed based on the existence of metabolic genes in the human genome. Cancer specific models of metabolism have also been generated by reducing the number of reactions in the generic model based on high throughput expression data, e.g. transcriptomics and proteomics. Targets for drugs and bio markers for diagnostics have been identified using these models. They have also been used as scaffolds for analysis of high throughput data to allow mechanistic interpretation of changes in expression. Finally, GEMs allow quantitative flux predictions using flux balance analysis (FBA). Here we critically review the requirements for successful FBA simulations of cancer cells and discuss the symmetry between the methods used for modeling of microbial and cancer metabolism. GEMs have great potential for translational research on cancer and will therefore become of increasing importance in the future.

A method for predicting target drug efficiency in cancer based on the analysis of signaling pathway activation

A new generation of anticancer therapeutics called target drugs has quickly developed in the 21st century. These drugs are tailored to inhibit cancer cell growth, proliferation, and viability by specific interactions with one or a few target proteins. However, despite formally known molecular targets for every "target" drug, patient response to treatment remains largely individual and unpredictable. Choosing the most effective personalized treatment remains a major challenge in oncology and is still largely trial and error. Here we present a novel approach for predicting target drug efficacy based on the gene expression signature of the individual tumor sample(s). The enclosed bioinformatic algorithm detects activation of intracellular regulatory pathways in the tumor in comparison to the corresponding normal tissues. According to the nature of the molecular targets of a drug, it predicts whether the drug can prevent cancer growth and survival in each individual case by blocking the abnormally activated tumor-promoting pathways or by reinforcing internal tumor suppressor cascades. To validate the method, we compared the distribution of predicted drug efficacy scores for five drugs (Sorafenib, Bevacizumab, Cetuximab, Sorafenib, Imatinib, Sunitinib) and seven cancer types (Clear Cell Renal Cell Carcinoma, Colon cancer, Lung adenocarcinoma, non-Hodgkin Lymphoma, Thyroid cancer and Sarcoma) with the available clinical trials data for the respective cancer types and drugs. The percent of responders to a drug treatment correlated significantly (Pearson's correlation 0.77 p = 0.023) with the percent of tumors showing high drug scores calculated with the current algorithm.

Gene Set-Based Integrative Analysis Revealing Two Distinct Functional Regulation Patterns in Four Common Subtypes of Epithelial Ovarian Cancer

Clear cell (CCC), endometrioid (EC), mucinous (MC) and high-grade serous carcinoma (SC) are the four most common subtypes of epithelial ovarian carcinoma (EOC). The widely accepted dualistic model of ovarian carcinogenesis divided EOCs into type I and II categories based on the molecular features. However, this hypothesis has not been experimentally demonstrated. We carried out a gene set-based analysis by integrating the microarray gene expression profiles downloaded from the publicly available databases. These quantified biological functions of EOCs were defined by 1454 Gene Ontology (GO) term and 674 Reactome pathway gene sets. The pathogenesis of the four EOC subtypes was investigated by hierarchical clustering and exploratory factor analysis. The patterns of functional regulation among the four subtypes containing 1316 cases could be accurately classified by machine learning. The results revealed that the ERBB and PI3K-related pathways played important roles in the carcinogenesis of CCC, EC and MC; while deregulation of cell cycle was more predominant in SC. The study revealed that two different functional regulation patterns exist among the four EOC subtypes, which were compatible with the type I and II classifications proposed by the dualistic model of ovarian carcinogenesis.

Gene Set-Based Functionome Analysis of Pathogenesis in Epithelial Ovarian Serous Carcinoma and the Molecular Features in Different FIGO Stages

Serous carcinoma (SC) is the most common subtype of epithelial ovarian carcinoma and is divided into four stages by the Federation of Gynecologists and Obstetrics (FIGO) staging system. Currently, the molecular functions and biological processes of SC at different FIGO stages have not been quantified. Here, we conducted a whole-genome integrative analysis to investigate the functions of SC at different stages. The function, as defined by the GO term or canonical pathway gene set, was quantified by measuring the changes in the gene expressional order between cancerous and normal control states. The quantified function, i.e., the gene set regularity (GSR) index, was utilized to investigate the pathogenesis and functional regulation of SC at different FIGO stages. We showed that the informativeness of the GSR indices was sufficient for accurate pattern recognition and classification for machine learning. The function regularity presented by the GSR indices showed stepwise deterioration during SC progression from FIGO stage I to stage IV. The pathogenesis of SC was centered on cell cycle deregulation and accompanied with multiple functional aberrations as well as their interactions.

Stratification of Hepatocellular Carcinoma Patients Based on Acetate Utilization

Hepatocellular carcinoma (HCC) is a deadly form of liver cancer that is increasingly prevalent. We analyzed global gene expression profiling of 361 HCC tumors and 49 adjacent noncancerous liver samples by means of combinatorial network-based analysis. We investigated the correlation between transcriptome and proteome of HCC and reconstructed a functional genome-scale metabolic model (GEM) for HCC. We identified fundamental metabolic processes required for cell proliferation using the network centric view provided by the GEM. Our analysis revealed tight regulation of fatty acid biosynthesis (FAB) and highly significant deregulation of fatty acid oxidation in HCC. We predicted mitochondrial acetate as an emerging substrate for FAB through upregulation of mitochondrial acetyl-CoA synthetase (ACSS1) in HCC. We analyzed heterogeneous expression of ACSS1 and ACSS2 between HCC patients stratified by high and low ACSS1 and ACSS2 expression and revealed that ACSS1 is associated with tumor growth and malignancy under hypoxic conditions in human HCC.

Subtype prediction in pediatric acute myeloid leukemia: classification using differential network rank conservation revisited

Background

One of the most important application spectrums of transcriptomic data is cancer phenotype classification. Many characteristics of transcriptomic data, such as redundant features and technical artifacts, make over-fitting commonplace. Promising classification results often fail to generalize across datasets with different sources, platforms, or preprocessing. Recently a novel differential network rank conservation (DIRAC) algorithm to characterize cancer phenotypes using transcriptomic data. DIRAC is a member of a family of algorithms that have shown useful for disease classification based on the relative expression of genes. Combining the robustness of this family’s simple decision rules with known biological relationships, this systems approach identifies interpretable, yet highly discriminate networks. While DIRAC has been briefly employed for several classification problems in the original paper, the potentials of DIRAC in cancer phenotype classification, and especially robustness against artifacts in transcriptomic data have not been fully characterized yet.

Results

In this study we thoroughly investigate the potentials of DIRAC by applying it to multiple datasets, and examine the variations in classification performances when datasets are (i) treated and untreated for batch effect; (ii) preprocessed with different techniques. We also propose the first DIRAC-based classifier to integrate multiple networks. We show that the DIRAC-based classifier is very robust in the examined scenarios. To our surprise, the trained DIRAC-based classifier even translated well to a dataset with different biological characteristics in the presence of substantial batch effects that, as shown here, plagued the standard expression value based classifier. In addition, the DIRAC-based classifier, because of the integrated biological information, also suggests pathways to target in specific subtypes, which may enhance the establishment of personalized therapy in diseases such as pediatric AML. In order to better comprehend the prediction power of the DIRAC-based classifier in general, we also performed classifications using publicly available datasets from breast and lung cancer. Furthermore, multiple well-known classification algorithms were utilized to create an ideal test bed for comparing the DIRAC-based classifier with the standard gene expression value based classifier. We observed that the DIRAC-based classifier greatly outperforms its rival.

Conclusions

Based on our experiments with multiple datasets, we propose that DIRAC is a promising solution to the lack of generalizability in classification efforts that uses transcriptomic data. We believe that superior performances presented in this study may motivate other to initiate a new aline of research to explore the untapped power of DIRAC in a broad range of cancer types.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0737-3) contains supplementary material, which is available to authorized users.

iCOSSY: An Online Tool for Context-Specific Subnetwork Discovery from Gene Expression Data

Pathway analyses help reveal underlying molecular mechanisms of complex biological phenotypes. Biologists tend to perform multiple pathway analyses on the same dataset, as there is no single answer. It is often inefficient for them to implement and/or install all the algorithms by themselves. Online tools can help the community in this regard. Here we present an online gene expression analytical tool called iCOSSY which implements a novel pathway-based COntext-specific Subnetwork discoverY (COSSY) algorithm. iCOSSY also includes a few modifications of COSSY to increase its reliability and interpretability. Users can upload their gene expression datasets, and discover important subnetworks of closely interacting molecules to differentiate between two phenotypes (context). They can also interactively visualize the resulting subnetworks. iCOSSY is a web server that finds subnetworks that are differentially expressed in two phenotypes. Users can visualize the subnetworks to understand the biology of the difference.

An argument for mechanism-based statistical inference in cancer

Cancer is perhaps the prototypical systems disease, and as such has been the focus of extensive study in quantitative systems biology. However, translating these programs into personalized clinical care remains elusive and incomplete. In this perspective, we argue that realizing this agenda—in particular, predicting disease phenotypes, progression and treatment response for individuals—requires going well beyond standard computational and bioinformatics tools and algorithms. It entails designing global mathematical models over network-scale configurations of genomic states and molecular concentrations, and learning the model parameters from limited available samples of high-dimensional and integrative omics data. As such, any plausible design should accommodate: biological mechanism, necessary for both feasible learning and interpretable decision making; stochasticity, to deal with uncertainty and observed variation at many scales; and a capacity for statistical inference at the patient level. This program, which requires a close, sustained collaboration between mathematicians and biologists, is illustrated in several contexts, including learning biomarkers, metabolism, cell signaling, network inference and tumorigenesis.

Cladograms with Path to Event (ClaPTE): a novel algorithm to detect associations between genotypes or phenotypes using phylogenies

BACKGROUND

Associations between genotype and phenotype provide insight into the evolution of pathogenesis, drug resistance, and the spread of pathogens between hosts. However, common ancestry can lead to apparent associations between biologically unrelated features. The novel method Cladograms with Path to Event (ClaPTE) detects associations between character-pairs (either a pair of mutations or a mutation paired with a phenotype) while adjusting for common ancestry, using phylogenetic trees.

METHODS

ClaPTE tests for character-pairs changing close together on the phylogenetic tree, consistent with an associated character-pair. ClaPTE is compared to three existing methods (independent contrasts, mixed model, and likelihood ratio) to detect character-pair associations adjusted for common ancestry. Comparisons utilize simulations on gene trees for: HIV Env, HIV promoter, and bacterial DnaJ and GuaB; and case studies for Oseltamavir resistance in Influenza, and for DnaJ and GuaB. Simulated data include both true-positive/associated character-pairs, and true-negative/not-associated character-pairs, used to assess type I (frequency of p-values in true-negatives) and type II (sensitivity to true-positives) error control.

RESULTS AND CONCLUSIONS

ClaPTE has competitive sensitivity and better type I error control than existing methods. In the Influenza/Oseltamavir case study, ClaPTE reports no new permissive mutations but detects associations between adjacent (in primary sequence) amino acid positions which other methods miss. In the DnaJ and GuaB case study, ClaPTE reports more frequent associations between positions both from the same protein family than between positions from different families, in contrast to other methods. In both case studies, the results from ClaPTE are biologically plausible.

IL-10 alters immunoproteostasis in APP mice, increasing plaque burden and worsening cognitive behavior

Anti-inflammatory strategies are proposed to have beneficial effects in Alzheimer's disease. To explore how anti-inflammatory cytokine signaling affects Aβ pathology, we investigated the effects of adeno-associated virus (AAV2/1)-mediated expression of Interleukin (IL)-10 in the brains of APP transgenic mouse models. IL-10 expression resulted in increased Aβ accumulation and impaired memory in APP mice. A focused transcriptome analysis revealed changes consistent with enhanced IL-10 signaling and increased ApoE expression in IL-10-expressing APP mice. ApoE protein was selectively increased in the plaque-associated insoluble cellular fraction, likely because of direct interaction with aggregated Aβ in the IL-10-expressing APP mice. Ex vivo studies also show that IL-10 and ApoE can individually impair glial Aβ phagocytosis. Our observations that IL-10 has an unexpected negative effect on Aβ proteostasis and cognition in APP mouse models demonstrate the complex interplay between innate immunity and proteostasis in neurodegenerative diseases, an interaction we call immunoproteostasis.

Exploiting dependencies of pairwise comparison outcomes to predict patterns of gene response

BACKGROUND

The analysis of gene expression has played an important role in medical and bioinformatics research. Although it is known that a large number of samples is needed to determine the patterns of gene expression accurately, practical designs of gene expression studies occasionally have insufficient numbers of samples, making it difficult to ascertain true response patterns of variantly expressed genes.

RESULTS

We describe an approach to cope with the challenge of predicting true orders of gene response to treatments. We show that true patterns of gene response must be orderable sets. In experiments with few samples, we modify the conventional pairwise comparison tests and increase the significance level α intelligently to deduce orderable patterns, which are most likely true orders of gene response. Additionally, motivated by the fact that a gene can be involved in multiple biological functions, our method further resamples experimental replicates and predicts multiple response patterns for each gene.

CONCLUSIONS

This method can be useful in designing cost-effective experiments with small sample sizes. Patterns of highly-variantly expressed genes can be predicted by varying α intelligently. Furthermore, clusters are labeled meaningfully with patterns that describe precisely how genes in such clusters respond to treatments.

Learning dysregulated pathways in cancers from differential variability analysis

Analysis of gene sets can implicate activity in signaling pathways that is responsible for cancer initiation and progression, but is not discernible from the analysis of individual genes. Multiple methods and software packages have been developed to infer pathway activity from expression measurements for set of genes targeted by that pathway. Broadly, three major methodologies have been proposed: over-representation, enrichment, and differential variability. Both over-representation and enrichment analyses are effective techniques to infer differentially regulated pathways from gene sets with relatively consistent differentially expressed (DE) genes. Specifically, these algorithms aggregate statistics from each gene in the pathway. However, they overlook multivariate patterns related to gene interactions and variations in expression. Therefore, the analysis of differential variability of multigene expression patterns can be essential to pathway inference in cancers. The corresponding methodologies and software packages for such multivariate variability analysis of pathways are reviewed here. We also introduce a new, computationally efficient algorithm, expression variation analysis (EVA), which has been implemented along with a previously proposed algorithm, Differential Rank Conservation (DIRAC), in an open source R package, gene set regulation (GSReg). EVA inferred similar pathways as DIRAC at reduced computational costs. Moreover, EVA also inferred different dysregulated pathways than those identified by enrichment analysis.

Edge biomarkers for classification and prediction of phenotypes

In general, a disease manifests not from malfunction of individual molecules but from failure of the relevant system or network, which can be considered as a set of interactions or edges among molecules. Thus, instead of individual molecules, networks or edges are stable forms to reliably characterize complex diseases. This paper reviews both traditional node biomarkers and edge biomarkers, which have been newly proposed. These biomarkers are classified in terms of their contained information. In particular, we show that edge and network biomarkers provide novel ways of stably and reliably diagnosing the disease state of a sample. First, we categorize the biomarkers based on the information used in the learning and prediction steps. We then briefly introduce conventional node biomarkers, or molecular biomarkers without network information, and their computational approaches. The main focus of this paper is edge and network biomarkers, which exploit network information to improve the accuracy of diagnosis and prognosis. Moreover, by extracting both network and dynamic information from the data, we can develop dynamical network and edge biomarkers. These biomarkers not only diagnose the immediate pre-disease state but also detect the critical molecules or networks by which the biological system progresses from the healthy to the disease state. The identified critical molecules can be used as drug targets, and the critical state indicates the critical point of disease control. The paper also discusses representative biomarker-based methods.

Identification of phenotype deterministic genes using systemic analysis of transcriptional response

Systemic identification of deterministic genes for different phenotypes is a primary application of high-throughput expression profiles. However, gene expression differences cannot be used when the differences between groups are not significant. Therefore, novel methods incorporating features other than expression differences are required. We developed a promising method using transcriptional response as an operational feature, which is quantified as the correlation between expression levels of pathway genes and target genes of the pathway. We applied this method to identify causative genes associated with chemo-sensitivity to tamoxifen and epirubicin. Genes whose transcriptional response was dysregulated only in the drug-resistant patient group were chosen for in vitro validation in human breast cancer cells. Finally, we discovered two genes responsible for tamoxifen sensitivity and three genes associated with epirubicin sensitivity. The method we propose here can be widely applied to identify deterministic genes for different phenotypes with only minor differences in gene expression levels.

Automatic context-specific subnetwork discovery from large interaction networks

Genes act in concert via specific networks to drive various biological processes, including progression of diseases such as cancer. Under different phenotypes, different subsets of the gene members of a network participate in a biological process. Single gene analyses are less effective in identifying such core gene members (subnetworks) within a gene set/network, as compared to gene set/network-based analyses. Hence, it is useful to identify a discriminative classifier by focusing on the subnetworks that correspond to different phenotypes. Here we present a novel algorithm to automatically discover the important subnetworks of closely interacting molecules to differentiate between two phenotypes (context) using gene expression profiles. We name it COSSY (COntext-Specific Subnetwork discoverY). It is a non-greedy algorithm and thus unlikely to have local optima problems. COSSY works for any interaction network regardless of the network topology. One added benefit of COSSY is that it can also be used as a highly accurate classification platform which can produce a set of interpretable features.

Identification of direction in gene networks from expression and methylation

BACKGROUND

Reverse-engineering gene regulatory networks from expression data is difficult, especially without temporal measurements or interventional experiments. In particular, the causal direction of an edge is generally not statistically identifiable, i.e., cannot be inferred as a statistical parameter, even from an unlimited amount of non-time series observational mRNA expression data. Some additional evidence is required and high-throughput methylation data can viewed as a natural multifactorial gene perturbation experiment.

RESULTS

We introduce IDEM (Identifying Direction from Expression and Methylation), a method for identifying the causal direction of edges by combining DNA methylation and mRNA transcription data. We describe the circumstances under which edge directions become identifiable and experiments with both real and synthetic data demonstrate that the accuracy of IDEM for inferring both edge placement and edge direction in gene regulatory networks is significantly improved relative to other methods.

CONCLUSION

Reverse-engineering directed gene regulatory networks from static observational data becomes feasible by exploiting the context provided by high-throughput DNA methylation data.An implementation of the algorithm described is available at http://code.google.com/p/idem/.

Transcriptional analysis of aggressiveness and heterogeneity across grades of astrocytomas

Astrocytoma is the most common glioma, accounting for half of all primary brain and spinal cord tumors. Late detection and the aggressive nature of high-grade astrocytomas contribute to high mortality rates. Though many studies identify candidate biomarkers using high-throughput transcriptomic profiling to stratify grades and subtypes, few have resulted in clinically actionable results. This shortcoming can be attributed, in part, to pronounced lab effects that reduce signature robustness and varied individual gene expression among patients with the same tumor. We addressed these issues by uniformly preprocessing publicly available transcriptomic data, comprising 306 tumor samples from three astrocytoma grades (Grade 2, 3, and 4) and 30 non-tumor samples (normal brain as control tissues). Utilizing Differential Rank Conservation (DIRAC), a network-based classification approach, we examined the global and individual patterns of network regulation across tumor grades. Additionally, we applied gene-based approaches to identify genes whose expression changed consistently with increasing tumor grade and evaluated their robustness across multiple studies using statistical sampling. Applying DIRAC, we observed a global trend of greater network dysregulation with increasing tumor aggressiveness. Individual networks displaying greater differences in regulation between adjacent grades play well-known roles in calcium/PKC, EGF, and transcription signaling. Interestingly, many of the 90 individual genes found to monotonically increase or decrease with astrocytoma grade are implicated in cancer-affected processes such as calcium signaling, mitochondrial metabolism, and apoptosis. The fact that specific genes monotonically increase or decrease with increasing astrocytoma grade may reflect shared oncogenic mechanisms among phenotypically similar tumors. This work presents statistically significant results that enable better characterization of different human astrocytoma grades and hopefully can contribute towards improvements in diagnosis and therapy choices. Our results also identify a number of testable hypotheses relating to astrocytoma etiology that may prove helpful in developing much-needed biomarkers for earlier disease detection.

Pathway analyses and understanding disease associations

High throughput technologies have been applied to investigate the underlying mechanisms of complex diseases, identify disease-associations and help to improve treatment. However it is challenging to derive biological insight from conventional single gene based analysis of "omics" data from high throughput experiments due to sample and patient heterogeneity. To address these challenges, many novel pathway and network based approaches were developed to integrate various "omics" data, such as gene expression, copy number alteration, Genome Wide Association Studies, and interaction data. This review will cover recent methodological developments in pathway analysis for the detection of dysregulated interactions and disease-associated subnetworks, prioritization of candidate disease genes, and disease classifications. For each application, we will also discuss the associated challenges and potential future directions.

Integrative analysis of cancer-related signaling pathways

Identification and classification of cancer types and subtypes is a major issue in current cancer research. Whole genome expression profiling of cancer tissues is often the basis for such subtype classifications of tumors and different signatures for individual cancer types have been described. However, the search for best performing discriminatory gene-expression signatures covering more than one cancer type remains a relevant topic in cancer research as such a signature would help understanding the common changes in signaling networks in these disease types. In this work, we explore the idea of a top down approach for sample stratification based on a module-based network of cancer relevant signaling pathways. For assembly of this network, we consider several of the most established cancer pathways. We evaluate our sample stratification approach using expression data of human breast and ovarian cancer signatures. We show that our approach performs equally well to previously reported methods besides providing the advantage to classify different cancer types. Furthermore, it allows to identify common changes in network module activity of those cancer samples.