Identifying Significantly Perturbed Subnetworks in Cancer Using Multiple Protein-Protein Interaction Networks.

BACKGROUND: The identification of cancer driver genes and key molecular pathways has been the focus of large-scale cancer genome studies. Network-based methods detect significantly perturbed subnetworks as putative cancer pathways by incorporating genomics data with the topological information of PPI networks. However, commonly used PPI networks have distinct topological structures, making the results of the same method vary widely when applied to different networks. Furthermore, emerging context-specific PPI networks often have incomplete topological structures, which pose serious challenges for existing subnetwork detection algorithms. METHODS: In this paper, we propose a novel method, referred to as MultiFDRnet, to address the above issues. The basic idea is to model a set of PPI networks as a multiplex network to preserve the topological structure of individual networks, while introducing dependencies among them, and, then, to detect significantly perturbed subnetworks on the modeled multiplex network using all the structural information simultaneously. RESULTS: To illustrate the effectiveness of the proposed approach, an extensive benchmark analysis was conducted on both simulated and real cancer data. The experimental results showed that the proposed method is able to detect significantly perturbed subnetworks jointly supported by multiple PPI networks and to identify novel modular structures in context-specific PPI networks.

Association of SNPs in the PAI1 Gene with Disease Recurrence and Clinical Outcome in Bladder Cancer.

PURPOSE: Bladder cancer (BCa) is one of the most common cancer types worldwide and is characterized by a high rate of recurrence. In previous studies, we and others have described the functional influence of plasminogen activator inhibitor-1 (PAI1) in bladder cancer development. While polymorphisms in PAI1 have been associated with increased risk and worsened prognosis in some cancers, the mutational status of PAI1 in human bladder tumors has not been well defined. METHODS: In this study, we evaluated the mutational status of PAI1 in a series of independent cohorts, comprised of a total of 660 subjects. RESULTS: Sequencing analyses identified two clinically relevant 3' untranslated region (UTR) single nucleotide polymorphisms (SNPs) in PAI1 (rs7242; rs1050813). Somatic SNP rs7242 was present in human BCa cohorts (overall incidence of 72%; 62% in Caucasians and 72% in Asians). In contrast, the overall incidence of germline SNP rs1050813 was 18% (39% in Caucasians and 6% in Asians). Furthermore, Caucasian patients with at least one of the described SNPs had worse recurrence-free survival and overall survival (p = 0.03 and p = 0.03, respectively). In vitro functional studies demonstrated that SNP rs7242 increased the anti-apoptotic effect of PAI1, and SNP rs1050813 was related to a loss of contact inhibition associated with cellular proliferation when compared to wild type. CONCLUSION: Further investigation of the prevalence and potential downstream influence of these SNPs in bladder cancer is warranted.

Alignment-free comparison of metagenomics sequences via approximate string matching.

Summary: Quantifying pairwise sequence similarities is a key step in metagenomics studies. Alignment-free methods provide a computationally efficient alternative to alignment-based methods for large-scale sequence analysis. Several neural network-based methods have recently been developed for this purpose. However, existing methods do not perform well on sequences of varying lengths and are sensitive to the presence of insertions and deletions. In this article, we describe the development of a new method, referred to as AsMac that addresses the aforementioned issues. We proposed a novel neural network structure for approximate string matching for the extraction of pertinent information from biological sequences and developed an efficient gradient computation algorithm for training the constructed neural network. We performed a large-scale benchmark study using real-world data that demonstrated the effectiveness and potential utility of the proposed method. Availability and implementation: The open-source software for the proposed method and trained neural-network models for some commonly used metagenomics marker genes were developed and are freely available at www.acsu.buffalo.edu/~yijunsun/lab/AsMac.html. Supplementary information: Supplementary data are available at Bioinformatics online.

Computational approach to modeling microbiome landscapes associated with chronic human disease progression.

A microbial community is a dynamic system undergoing constant change in response to internal and external stimuli. These changes can have significant implications for human health. However, due to the difficulty in obtaining longitudinal samples, the study of the dynamic relationship between the microbiome and human health remains a challenge. Here, we introduce a novel computational strategy that uses massive cross-sectional sample data to model microbiome landscapes associated with chronic disease development. The strategy is based on the rationale that each static sample provides a snapshot of the disease process, and if the number of samples is sufficiently large, the footprints of individual samples populate progression trajectories, which enables us to recover disease progression paths along a microbiome landscape by using computational approaches. To demonstrate the validity of the proposed strategy, we developed a bioinformatics pipeline and applied it to a gut microbiome dataset available from a Crohn's disease study. Our analysis resulted in one of the first working models of microbial progression for Crohn's disease. We performed a series of interrogations to validate the constructed model. Our analysis suggested that the model recapitulated the longitudinal progression of microbial dysbiosis during the known clinical trajectory of Crohn's disease. By overcoming restrictions associated with complex longitudinal sampling, the proposed strategy can provide valuable insights into the role of the microbiome in the pathogenesis of chronic disease and facilitate the shift of the field from descriptive research to mechanistic studies.

A Diagnostic Gene Expression Signature for Bladder Cancer Can Stratify Cases into Prescribed Molecular Subtypes and Predict Outcome.

Bladder cancer is a biologically heterogeneous disease with variable clinical presentations, outcomes and responses to therapy. Thus, the clinical utility of single biomarkers for the detection and prediction of biological behavior of bladder cancer is limited. We have previously identified and validated a bladder cancer diagnostic signature composed of 10 biomarkers, which has been incorporated into a multiplex immunoassay bladder cancer test, Oncuria™. In this study, we evaluate whether these 10 biomarkers can assist in the prediction of bladder cancer clinical outcomes. Tumor gene expression and patient survival data from bladder cancer cases from The Cancer Genome Atlas (TCGA) were analyzed. Alignment between the mRNA expression of 10 biomarkers and the TCGA 2017 subtype classification was assessed. Kaplan-Meier analysis of multiple gene expression datasets indicated that high expression of the combined 10 biomarkers correlated with a significant reduction in overall survival. The analysis of three independent, publicly available gene expression datasets confirmed that multiplex prognostic models outperformed single biomarkers. In total, 8 of the 10 biomarkers from the Oncuria™ test were significantly associated with either luminal or basal molecular subtypes, and thus, the test has the potential to assist in the prediction of clinical outcome.

PAI-1 is a potential transcriptional silencer that supports bladder cancer cell activity.

The extracellular activity of Plasminogen activator inhibitor-1 (PAI-1) is well described, acting as an inhibitor of tissue plasminogen activator and urokinase-type plasminogen activator, impacting fibrinolysis. Recent studies have revealed a pro-tumorigenic role of PAI-1 in human cancers, via the regulation of angiogenesis and tumor cell survival. In this study, immunohistochemical staining of 939 human bladder cancer specimens showed that PAI-1 expression levels correlated with tumor grade, tumor stage and overall survival. The typical subcellular localization of PAI-1 is cytoplasmic, but in approximately a quarter of the cases, PAI-1 was observed to be localized to both the tumor cell cytoplasm and the nucleus. To investigate the potential function of nuclear PAI-1 in tumor biology we applied chromatin immunoprecipitation (ChIP)-sequencing, gene expression profiling, and rapid immunoprecipitation mass spectrometry to a pair of bladder cancer cell lines. ChIP-sequencing revealed that PAI-1 can bind DNA at distal intergenic regions, suggesting a role as a transcriptional coregulator. The downregulation of PAI-1 in bladder cancer cell lines caused the upregulation of numerous genes, and the integration of ChIP-sequence and RNA-sequence data identified 57 candidate genes subject to PAI-1 regulation. Taken together, the data suggest that nuclear PAI-1 can influence gene expression programs and support malignancy.

Application of a multiplex urinalysis test for the prediction of intravesical BCG treatment response: A pilot study.

BACKGROUND: Intravesical Bacillus Calmette-Guerin (BCG), a live attenuated tuberculosis vaccine that acts as a non-specific immune system stimulant, is the most effective adjuvant treatment for patients with intermediate or high-risk non-muscle-invasive bladder cancer (NMIBC). However, to date, there are no reliable tests that are predictive of BCG treatment response. In this study, we evaluated the performance of OncuriaTM, a bladder cancer detection test, to predict response to intravesical BCG. METHODS: OncuriaTM data was evaluated in voided urine samples obtained from a prospectively collected cohort of 64 subjects with intermediate or high risk NMIBC prior to treatment with intravesical BCG. The OncuriaTM test, which measures 10 cancer-associated biomarkers was performed in an independent clinical laboratory. The ability of the test to identify those patients in whom BCG is ineffective against tumor recurrence was tested. Predictive models were derived using supervised learning and cross-validation analyses. Model performance was assessed using ROC curves. RESULTS: Pre-treatment urinary concentrations of MMP9, VEGFA, CA9, SDC1, PAI1, APOE, A1AT, ANG and MMP10 were increased in patients who developed disease recurrence. A combinatorial predictive model of treatment outcome achieved an AUROC 0.89 [95% CI: 0.80-0.99], outperforming any single biomarker, with a test sensitivity of 81.8% and a specificity of 84.9%. Hazard ratio analysis revealed that patients with higher urinary levels of ANG, CA9 and MMP10 had a significantly higher risk of disease recurrence. CONCLUSIONS: Monitoring the urinary levels of a cancer-associated biomarker panel enabled the discrimination of patients who did not respond to intravesical BCG therapy. With further study, the multiplex OncuriaTM test may be applicable for the clinical evaluation of bladder cancer patients considering intravesical BCG treatment.

Influencing Factors on the Oncuria™ Urinalysis Assay: An Experimental Model.

BACKGROUND: The Oncuria™ urine test for the detection of bladder cancer measures a multiplex protein signature. In this study, we investigated the influence of urinary cellularity, protein, and hematuria on the performance of the Oncuria™ test in an ex vivo experimental model. MATERIALS AND METHODS: Pooled urine from healthy subjects was spiked with cultured benign (UROtsa) or malignant cells (T24), cellular proteins, or whole blood. The resulting samples were analyzed using the Oncuria™ test following the manufacturer's instructions. RESULTS: Urine samples obtained from healthy subjects were negative for bladder cancer by Oncuria™ test criteria. The majority of the manipulated conditions did not result in a false-positive test. The addition of whole blood (high concentration) did result in a false-positive result, but this was abrogated by sample centrifugation prior to analysis. The addition of cellular proteins (high concentration) resulted in a positive Oncuria™ test, and this was unaffected by pre-analysis sample centrifugation. CONCLUSIONS: The Oncuria™ multiplex test performed well in the ex vivo experimental model and shows promise for clinical application. The identification of patients who require additional clinical evaluation could reduce the need to subject patients who do not have bladder cancer to frequent, uncomfortable and expensive cystoscopic examinations, thus benefiting both patients and the healthcare system.

Diagnostic performance of Oncuria™, a urinalysis test for bladder cancer.

BACKGROUND: Due to insufficient accuracy, urine-based assays currently have a limited role in the management of patients with bladder cancer. The identification of multiplex molecular signatures associated with disease has the potential to address this deficiency and to assist with accurate, non-invasive diagnosis and monitoring. METHODS: To evaluate the performance of Oncuria™, a multiplex immunoassay for bladder detection in voided urine samples. The test was evaluated in a multi-institutional cohort of 362 prospectively collected subjects presenting for bladder cancer evaluation. The parallel measurement of 10 biomarkers (A1AT, APOE, ANG, CA9, IL8, MMP9, MMP10, PAI1, SDC1 and VEGFA) was performed in an independent clinical laboratory. The ability of the test to identify patients harboring bladder cancer was assessed. Bladder cancer status was confirmed by cystoscopy and tissue biopsy. The association of biomarkers and demographic factors was evaluated using linear discriminant analysis (LDA) and predictive models were derived using supervised learning and cross-validation analyses. Diagnostic performance was assessed using ROC curves. RESULTS: The combination of the 10 biomarkers provided an AUROC 0.93 [95% CI 0.87-0.98], outperforming any single biomarker. The addition of demographic data (age, sex, and race) into a hybrid signature improved the diagnostic performance AUROC 0.95 [95% CI 0.90-1.00]. The hybrid signature achieved an overall sensitivity of 0.93, specificity of 0.93, PPV of 0.65 and NPV of 0.99 for bladder cancer classification. Sensitivity values of the diagnostic panel for high-grade bladder cancer, low-grade bladder cancer, MIBC and NMIBC were 0.94, 0.89, 0.97 and 0.93, respectively. CONCLUSIONS: Urinary levels of a biomarker panel enabled the accurate discrimination of bladder cancer patients and controls. The multiplex Oncuria™ test can achieve the efficient and accurate detection and monitoring of bladder cancer in a non-invasive patient setting.

An efficient and effective method to identify significantly perturbed subnetworks in cancer.

The identification of key functional biological networks from high-dimensional genomics data is pivotal for cancer research. Here, we introduce FDRnet, a method for the detection of molecular subnetworks in cancer, which addresses several challenges in pathway analysis. FDRnet detects key subnetworks by solving a mixed-integer linear programming problem, using a given upper bound of false discovery rate (FDR) as a budget constraint, and minimizing a conductance score to find dense subgraphs around seed genes. A large-scale benchmark study was performed on both simulation and cancer genomics data. FDRnet outperformed other methods in the ability to detect functionally homogeneous subnetworks in a scale-free biological network, to control FDRs of the genes in detected subnetworks, to improve computational efficiency and to integrate multi-omics data. By overcoming the limitations of existing approaches, FDRnet can facilitate the detection of key functional pathways in cancer and other genetic diseases.

Deep-learning approach to identifying cancer subtypes using high-dimensional genomic data.

MOTIVATION: Cancer subtype classification has the potential to significantly improve disease prognosis and develop individualized patient management. Existing methods are limited by their ability to handle extremely high-dimensional data and by the influence of misleading, irrelevant factors, resulting in ambiguous and overlapping subtypes. RESULTS: To address the above issues, we proposed a novel approach to disentangling and eliminating irrelevant factors by leveraging the power of deep learning. Specifically, we designed a deep-learning framework, referred to as DeepType, that performs joint supervised classification, unsupervised clustering and dimensionality reduction to learn cancer-relevant data representation with cluster structure. We applied DeepType to the METABRIC breast cancer dataset and compared its performance to state-of-the-art methods. DeepType significantly outperformed the existing methods, identifying more robust subtypes while using fewer genes. The new approach provides a framework for the derivation of more accurate and robust molecular cancer subtypes by using increasingly complex, multi-source data. AVAILABILITY AND IMPLEMENTATION: An open-source software package for the proposed method is freely available at http://www.acsu.buffalo.edu/~yijunsun/lab/DeepType.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Molecular Profiles of Matched Primary and Metastatic Tumor Samples Support a Linear Evolutionary Model of Breast Cancer.

The interpretation of accumulating genomic data with respect to tumor evolution and cancer progression requires integrated models. We developed a computational approach that enables the construction of disease progression models using static sample data. Application to breast cancer data revealed a linear, branching evolutionary model with two distinct trajectories for malignant progression. Here, we used the progression model as a foundation to investigate the relationships between matched primary and metastasis breast tumor samples. Mapping paired data onto the model confirmed that molecular breast cancer subtypes can shift during progression and supported directional tumor evolution through luminal subtypes to increasingly malignant states. Cancer progression modeling through the analysis of available static samples represents a promising breakthrough. Further refinement of a roadmap of breast cancer progression will facilitate the development of improved cancer diagnostics, prognostics, and targeted therapeutics. SIGNIFICANCE: Analysis of matched primary and metastatic tumor samples supports a unidirectional, linear cancer evolution process and sheds light on longstanding issues regarding the origins of molecular subtypes and their progression relationships.

Monoclonal Antibody against CXCL1 (HL2401) as a Novel Agent in Suppressing IL6 Expression and Tumoral Growth.

Rationale: The expression of the chemokine (C-X-C motif) ligand 1 (CXCL1), an inflammatory protein, has been reported to be up-regulated in many human cancers. The mechanisms through which aberrant cellular CXCL1 levels promote specific steps in tumor growth and progression are unknown. Methods: We described the anticancer effects and mechanism of action of HL2401, a monoclonal antibody directed at CXCL1 with in vitro and in vivo data on bladder and prostate cancers. Results: HL2401 inhibited proliferation and invasion of bladder and prostate cells along with disrupting endothelial sprouting in vitro. Furthermore, novel mechanistic investigations revealed that CXCL1 expression stimulated interleukin 6 (IL6) expression and repressed tissue inhibitor of metalloproteinase 4 (TIMP4). Systemic administration of HL2401 in mice bearing bladder and prostate xenograft tumors retarded tumor growth through the inhibition of cellular proliferation and angiogenesis along with an induction of apoptosis. Our findings reveal a previously undocumented relationship between CXCL1, IL6 and TIMP4 in solid tumor biology. Principal conclusions: Taken together, our results argue that CXCL1 plays an important role in sustaining the growth of bladder and prostate tumors via up-regulation of IL6 and down-regulation of TIMP4. Targeting these critical interactions with a CXCL1 monoclonal antibody offers a novel strategy to therapeutically manage bladder and prostate cancers.

Meta-analysis of a 10-plex urine-based biomarker assay for the detection of bladder cancer.

A 10-plex urine-based bladder cancer (BCa) diagnostic signature has the potential to non-invasively predict the presence of BCa in at-risk patients, as reported in various case-control studies. The present meta-analysis was performed to re-evaluate and demonstrate the robustness and consistency of the diagnostic utility of the 10-plex urine-based diagnostic assay. We re-analyzed primary data collected in five previously published case-control studies on the 10-plex diagnostic assay. Studies reported the sensitivity and specificity of ten urinary protein biomarkers for the detection of BCa, including interleukin 8, matrix metalloproteinases 9 and 10, angiogenin, apolipoprotein E, syndecan 1, alpha-1 antitrypsin, plasminogen activator inhibitor-1, carbonic anhydrase 9, and vascular endothelial growth factor A. Data were extracted and reviewed independently by two investigators. Log odds ratios (ORs) were calculated to determine how strongly the 10-plex biomarker panel and individual biomarkers are associated with the presence of BCa. Data pooled from 1,173 patients were analyzed. The log OR for each biomarker was improved by 1.5 or greater with smaller 95% CI in our meta-analysis of the overall cohort compared with each analysis of an individual cohort. The combination of the ten biomarkers showed a higher log OR (log OR: 3.46, 95% CI: 2.60-4.31) than did any single biomarker irrespective of histological grade or disease stage of tumors. We concluded that the 10-plex BCa-associated diagnostic signature demonstrated a higher potential to identify BCa when compared to any single biomarker. Our results justify further advancement of the 10-plex protein-based diagnostic signature toward clinical application.

Genetic deletion of sphingosine kinase 1 suppresses mouse breast tumor development in an HER2 transgenic model.

Aberrant sphingolipid metabolism has been reported to promote breast cancer progression. Sphingosine kinase 1 (SphK1) is a key metabolic enzyme for the formation of pro-survival S1P from pro-apoptotic ceramide. The role of SphK1 in breast cancer has been well studied in estrogen receptor (ER)-positive breast cancer; however, its role in human epidermal growth factor 2 (HER2)-positive breast cancer remains unclear. Here, we show that genetic deletion of SphK1 significantly reduced mammary tumor development with reduced tumor incidence and multiplicity in the MMTV-neu transgenic mouse model. Gene expression analysis revealed significant reduction of claudin-2 (CLDN2) expression in tumors from SphK1 deficient mice, suggesting that CLDN2 may mediate SphK1's function. It is remarkable that SphK1 deficiency in HER2-positive breast cancer model inhibited tumor formation by the different mechanism from ER-positive breast cancer. In vitro experiments demonstrated that overexpression of SphK1 in ER-/PR-/HER2+ human breast cancer cells enhanced cell proliferation, colony formation, migration and invasion. Furthermore, immunostaining of SphK1 and CLDN2 in HER2-positive human breast tumors revealed a correlation in high-grade disease. Taken together, these findings suggest that SphK1 may play a pivotal role in HER2-positive breast carcinogenesis. Targeting SphK1 may represent a novel approach for HER2-positive breast cancer chemoprevention and/or treatment.

Prognostic cancer gene signatures share common regulatory motifs.

Scientists have discovered various prognostic gene signatures (GSs) in different cancer types. Surprisingly, although different GSs from the same cancer type can be used to measure similar biological characteristics, often rarely is there a gene shared by different GSs. To explain such a paradox, we hypothesized that GSs from the same cancer type may be regulated by common regulatory motifs. To test this hypothesis, we carried out a comprehensive motif analysis on the prognostic GSs from five cancer types. We demonstrated that GSs from individual cancer type as well as across cancer types share regulatory motifs. We also observed that transcription factors that likely bind to these shared motifs have prognostic functions in cancers. Moreover, 75% of the predicted cofactors of these transcription factors may have cancer-related functions and some cofactors even have prognostic functions. In addition, there exist common microRNAs that regulate different GSs from individual cancer types and across cancer types, several of which are prognostic biomarkers for the corresponding cancer types. Our study suggested the existence of common regulatory mechanisms shared by GSs from individual cancer types and across cancer types, which shed light on the discovery of new prognostic GSs in cancers and the understanding of the regulatory mechanisms of cancers.

Tolfenamic acid-induced alterations in genes and pathways in pancreatic cancer cells.

Non-steroidal anti-inflammatory drugs (NSAIDs) are being tested extensively for their role in the treatment and prevention of several cancers. Typically NSAIDs exhibit anti-tumor activities via modulation of cyclooxygenase (COX)-dependent mechanisms, however, an anti-cancer NSAID tolfenamic acid (TA) is believed to work through COX-independent pathways. Results from our laboratory and others have demonstrated the anti-cancer activity of TA in various cancer models including pancreatic cancer. TA has been shown to modulate certain cellular processes including, apoptosis, reactive oxygen species and signaling. In this study, molecular profiling was performed to precisely understand the mode of action of TA. Three pancreatic cancer cell lines, L3.6pl, MIA PaCa-2, and Panc1 were treated with TA (50 µM for 48 h) and the changes in gene expression was evaluated using the Affymetrix GeneChip Human Gene ST Array platform. Microarray results were further validated using quantitative PCR for seven genes altered by TA treatment in all three cell lines. Functional analysis of differentially expressed genes (2 fold increase or decrease, p < 0.05) using Ingenuity Pathway Analysis software, revealed that TA treatment predominantly affected the genes involved in cell cycle, cell growth and proliferation, and cell death and survival. Promoter analysis of the differentially expressed genes revealed that they are enriched for Sp1 binding sites, suggesting that Sp1 could be a major contributor in mediating the effect of TA. The gene expression studies identified new targets involved in TA's mode of action, while supporting the hypothesis about the association of Sp1 in TA mediated effects in pancreatic cancer.

Computational approach for deriving cancer progression roadmaps from static sample data.

As with any biological process, cancer development is inherently dynamic. While major efforts continue to catalog the genomic events associated with human cancer, it remains difficult to interpret and extrapolate the accumulating data to provide insights into the dynamic aspects of the disease. Here, we present a computational strategy that enables the construction of a cancer progression model using static tumor sample data. The developed approach overcame many technical limitations of existing methods. Application of the approach to breast cancer data revealed a linear, branching model with two distinct trajectories for malignant progression. The validity of the constructed model was demonstrated in 27 independent breast cancer data sets, and through visualization of the data in the context of disease progression we were able to identify a number of potentially key molecular events in the advance of breast cancer to malignancy.

A microRNA biomarker panel for the non-invasive detection of bladder cancer.

The development of accurate, non-invasive urinary assays for bladder cancer would greatly facilitate the detection and management of a disease that has a high rate of recurrence and progression. In this study, we employed a discovery and validation strategy to identify microRNA signatures that can perform as a non-invasive bladder cancer diagnostic assay. Expression profiling of 754 human microRNAs (TaqMan low density arrays) was performed on naturally voided urine samples from a cohort of 85 subjects of known bladder disease status (27 with active BCa). A panel of 46 microRNAs significantly associated with bladder cancer were subsequently monitored in an independent cohort of 121 subjects (61 with active BCa) using quantitative real-time PCR (RT-PCR). Multivariable modeling identified a 25-target diagnostic signature that predicted the presence of BCa with an estimated sensitivity of 87% at a specificity of 100% (AUC 0.982). With additional validation, the monitoring of a urinary microRNA biomarker panel could facilitate the non-invasive evaluation of patients under investigation for BCa.