Inferring causal molecular networks: empirical assessment through a community-based effort.

It remains unclear whether causal, rather than merely correlational, relationships in molecular networks can be inferred in complex biological settings. Here we describe the HPN-DREAM network inference challenge, which focused on learning causal influences in signaling networks. We used phosphoprotein data from cancer cell lines as well as in silico data from a nonlinear dynamical model. Using the phosphoprotein data, we scored more than 2,000 networks submitted by challenge participants. The networks spanned 32 biological contexts and were scored in terms of causal validity with respect to unseen interventional data. A number of approaches were effective, and incorporating known biology was generally advantageous. Additional sub-challenges considered time-course prediction and visualization. Our results suggest that learning causal relationships may be feasible in complex settings such as disease states. Furthermore, our scoring approach provides a practical way to empirically assess inferred molecular networks in a causal sense.

Multi-omic analyses of changes in the tumor microenvironment of pancreatic adenocarcinoma following neoadjuvant treatment with anti-PD-1 therapy.

Successful pancreatic ductal adenocarcinoma (PDAC) immunotherapy necessitates optimization and maintenance of activated effector T cells (Teff). We prospectively collected and applied multi-omic analyses to paired pre- and post-treatment PDAC specimens collected in a platform neoadjuvant study of granulocyte-macrophage colony-stimulating factor-secreting allogeneic PDAC vaccine (GVAX) vaccine ± nivolumab (anti-programmed cell death protein 1 [PD-1]) to uncover sensitivity and resistance mechanisms. We show that GVAX-induced tertiary lymphoid aggregates become immune-regulatory sites in response to GVAX + nivolumab. Higher densities of tumor-associated neutrophils (TANs) following GVAX + nivolumab portend poorer overall survival (OS). Increased T cells expressing CD137 associated with cytotoxic Teff signatures and correlated with increased OS. Bulk and single-cell RNA sequencing found that nivolumab alters CD4+ T cell chemotaxis signaling in association with CD11b+ neutrophil degranulation, and CD8+ T cell expression of CD137 was required for optimal T cell activation. These findings provide insights into PD-1-regulated immune pathways in PDAC that should inform more effective therapeutic combinations that include TAN regulators and T cell activators.

Cytokines secreted by stromal cells in TNBC microenvironment as potential targets for cancer therapy.

In triple-negative breast cancer (TNBC), the lack of therapeutic markers and effective targeted therapies result in an incurable metastatic disease associated with a poor prognosis. Crosstalks within the tumor microenvironment (TME), including those between cancer and stromal cells, affect the tumor heterogeneity, growth, and metastasis. Previously, we have demonstrated that IL-6, IL-8, and CCL5 play a significant role in TNBC growth and metastasis. In this study, we performed a systematic analysis of cytokine factors secreted from four stromal components (fibroblasts, macrophages, lymphatic endothelial cells, and blood microvascular endothelial cells) induced by four TNBC cell types. Through bioinformatic analysis, we selected putative candidates of secreted factors from stromal cells, which are involved in EMT activity, cell proliferation, metabolism, and matrisome pathways. Among the candidates, LCN2, GM-CSF, CST3, IL-6, IL-8, and CHI3L1 are ranked highly. Significantly, Lipocalin-2 (LCN2) is upregulated in the crosstalk of stromal cells and four different TNBC cells. We validated the increase of LCN2 secreted from four stromal cells induced by TNBC cells. Using a specific LCN2 antibody, we observed the inhibition of TNBC cell growth and migration. Taken together, these results propose secreted factors as molecular targets to treat TNBC progression via crosstalk with stromal components.

Correction: Novel Insight into Mutational Landscape of Head and Neck Squamous Cell Carcinoma.

[This corrects the article DOI: 10.1371/journal.pone.0093102.].

Epigenetic regulation of gene expression in cancer: techniques, resources and analysis.

Cancer is a complex disease, driven by aberrant activity in numerous signaling pathways in even individual malignant cells. Epigenetic changes are critical mediators of these functional changes that drive and maintain the malignant phenotype. Changes in DNA methylation, histone acetylation and methylation, noncoding RNAs, posttranslational modifications are all epigenetic drivers in cancer, independent of changes in the DNA sequence. These epigenetic alterations were once thought to be crucial only for the malignant phenotype maintenance. Now, epigenetic alterations are also recognized as critical for disrupting essential pathways that protect the cells from uncontrolled growth, longer survival and establishment in distant sites from the original tissue. In this review, we focus on DNA methylation and chromatin structure in cancer. The precise functional role of these alterations is an area of active research using emerging high-throughput approaches and bioinformatics analysis tools. Therefore, this review also describes these high-throughput measurement technologies, public domain databases for high-throughput epigenetic data in tumors and model systems and bioinformatics algorithms for their analysis. Advances in bioinformatics data that combine these epigenetic data with genomics data are essential to infer the function of specific epigenetic alterations in cancer. These integrative algorithms are also a focus of this review. Future studies using these emerging technologies will elucidate how alterations in the cancer epigenome cooperate with genetic aberrations during tumor initiation and progression. This deeper understanding is essential to future studies with epigenetics biomarkers and precision medicine using emerging epigenetic therapies.

Functional characterization of alternatively spliced GSN in head and neck squamous cell carcinoma.

We have recently performed the characterization of alternative splicing events (ASEs) in head and neck squamous cell carcinoma, which allows dysregulation of protein expression common for cancer cells. Such analysis demonstrated a high ASE prevalence among tumor samples, including tumor-specific alternative splicing in the GSN gene.In vitro studies confirmed that overall expression of either ASE-GSN or wild-type GSN (WT-GSN) isoform inversely correlated with cell proliferation, whereas the high ratio of ASE-GSN to WT-GSN correlated with increased cellular invasion. Additionally, a change in expression of either isoform caused compensatory changes in expression of the other isoform. Our results suggest that the overall expression and the balance between GSN isoforms are mediating factors in proliferation, while increased overall expression of ASE-GSN is specific to cancer tissues. As a result, we propose ASE-GSN can serve not only as a biomarker of disease and disease progression, but also as a neoantigen for head and neck squamous cell carcinoma treatment, for which only a limited number of disease-specific targeted therapies currently exist.

Integrated time course omics analysis distinguishes immediate therapeutic response from acquired resistance.

BACKGROUND: Targeted therapies specifically act by blocking the activity of proteins that are encoded by genes critical for tumorigenesis. However, most cancers acquire resistance and long-term disease remission is rarely observed. Understanding the time course of molecular changes responsible for the development of acquired resistance could enable optimization of patients' treatment options. Clinically, acquired therapeutic resistance can only be studied at a single time point in resistant tumors. METHODS: To determine the dynamics of these molecular changes, we obtained high throughput omics data (RNA-sequencing and DNA methylation) weekly during the development of cetuximab resistance in a head and neck cancer in vitro model. The CoGAPS unsupervised algorithm was used to determine the dynamics of the molecular changes associated with resistance during the time course of resistance development. RESULTS: CoGAPS was used to quantify the evolving transcriptional and epigenetic changes. Applying a PatternMarker statistic to the results from CoGAPS enabled novel heatmap-based visualization of the dynamics in these time course omics data. We demonstrate that transcriptional changes result from immediate therapeutic response or resistance, whereas epigenetic alterations only occur with resistance. Integrated analysis demonstrates delayed onset of changes in DNA methylation relative to transcription, suggesting that resistance is stabilized epigenetically. CONCLUSIONS: Genes with epigenetic alterations associated with resistance that have concordant expression changes are hypothesized to stabilize the resistant phenotype. These genes include FGFR1, which was associated with EGFR inhibitors resistance previously. Thus, integrated omics analysis distinguishes the timing of molecular drivers of resistance. This understanding of the time course progression of molecular changes in acquired resistance is important for the development of alternative treatment strategies that would introduce appropriate selection of new drugs to treat cancer before the resistant phenotype develops.

DNA methylation regulates TMEM16A/ANO1 expression through multiple CpG islands in head and neck squamous cell carcinoma.

ANO1 is a calcium-activated chloride channel that is frequently overexpressed in head and neck squamous cell carcinoma (HNSCC) and other cancers. While ANO1 expression negatively correlates with survival in several cancers, its epigenetic regulation is poorly understood. We analyzed HNSCC samples from TCGA and a separate dataset of HPV+ oropharyngeal squamous cell carcinoma (OPSCC) samples to identify differentially methylated regions. E6 and E7 transfected normal oral keratinocytes (NOK) were used to induce hypermethylation of the ANO1 promoter. We found three CpG islands that correlated with ANO1 expression, including two positively correlated with expression. Using two HNSCC datasets with differential expression of ANO1, we showed hypermethylation of positively correlated CpG islands potentiates ANO1 expression. E7 but not E6 transfection of NOK cells led to hypermethylation of a positively correlated CpG island without a change in ANO1 expression. ANO1 promoter methylation was also correlated with patient survival. Our results are the first to show the contribution of positively correlated CpG's for regulating gene expression in HNSCC. Hypermethylation of the ANO1 promoter was strongly correlated with but not sufficient to increase ANO1 expression, suggesting methylation of positively correlated CpG's likely serves as an adjunct to other mechanisms of ANO1 activation.

Integrated Analysis of Whole-Genome ChIP-Seq and RNA-Seq Data of Primary Head and Neck Tumor Samples Associates HPV Integration Sites with Open Chromatin Marks.

Chromatin alterations mediate mutations and gene expression changes in cancer. Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) has been utilized to study genome-wide chromatin structure in human cancer cell lines, yet numerous technical challenges limit comparable analyses in primary tumors. Here we have developed a new whole-genome analytic pipeline to optimize ChIP-Seq protocols on patient-derived xenografts from human papillomavirus-related (HPV+) head and neck squamous cell carcinoma (HNSCC) samples. We further associated chromatin aberrations with gene expression changes from a larger cohort of the tumor and normal samples with RNA-Seq data. We detect differential histone enrichment associated with tumor-specific gene expression variation, sites of HPV integration in the human genome, and HPV-associated histone enrichment sites upstream of cancer driver genes, which play central roles in cancer-associated pathways. These comprehensive analyses enable unprecedented characterization of the complex network of molecular changes resulting from chromatin alterations that drive HPV-related tumorigenesis. Cancer Res; 77(23); 6538-50. ©2017 AACR.

RNAKinetics: a web server that models secondary structure kinetics of an elongating RNA.

The RNAKinetics server (http://www.ig-msk.ru/RNA/kinetics) is a web interface for the newly developed RNAKinetics software. The software models the dynamics of RNA secondary structure by the means of kinetic analysis of folding transitions of a growing RNA molecule. The result of the modeling is a kinetic ensemble, i.e. a collection of RNA structures that are endowed with probabilities, which depend on time. This approach gives comprehensive probabilistic description of RNA folding pathways, revealing important kinetic details that are not captured by the traditional structure prediction methods. The access to the RNAKinetics server is free.

Novel insight into mutational landscape of head and neck squamous cell carcinoma.

Development of head and neck squamous cell carcinoma (HNSCC) is characterized by accumulation of mutations in several oncogenes and tumor suppressor genes. We have formerly described the mutation pattern of HNSCC and described NOTCH signaling pathway alterations. Given the complexity of the HNSCC, here we extend the previous study to understand the overall HNSCC mutation context and to discover additional genetic alterations. We performed high depth targeted exon sequencing of 51 highly actionable cancer-related genes with a high frequency of mutation across many cancer types, including head and neck. DNA from primary tumor tissues and matched normal tissues was analyzed for 37 HNSCC patients. We identified 26 non-synonymous or stop-gained mutations targeting 11 of 51 selected genes. These genes were mutated in 17 out of 37 (46%) studied HNSCC patients. Smokers harbored 3.2-fold more mutations than non-smokers. Importantly, TP53 was mutated in 30%, NOTCH1 in 8% and FGFR3 in 5% of HNSCC. HPV negative patients harbored 4-fold more TP53 mutations than HPV positive patients. These data confirm prior reports of the HNSCC mutational profile. Additionally, we detected mutations in two new genes, CEBPA and FES, which have not been previously reported in HNSCC. These data extend the spectrum of HNSCC mutations and define novel mutation targets in HNSCC carcinogenesis, especially for smokers and HNSCC without HPV infection.

Preferential activation of the hedgehog pathway by epigenetic modulations in HPV negative HNSCC identified with meta-pathway analysis.

Head and neck squamous cell carcinoma (HNSCC) is largely divided into two groups based on their etiology, human papillomavirus (HPV)-positive and -negative. Global DNA methylation changes are known to drive oncogene and tumor suppressor expression in primary HNSCC of both types. However, significant heterogeneity in DNA methylation within the groups results in different transcriptional profiles and clinical outcomes. We applied a meta-pathway analysis to link gene expression changes to DNA methylation in distinguishing HNSCC subtypes. This approach isolated specific epigenetic changes controlling expression in HPV- HNSCC that distinguish it from HPV+ HNSCC. Analysis of genes identified Hedgehog pathway activation specific to HPV- HNSCC. We confirmed that GLI1, the primary Hedgehog target, showed higher expression in tumors compared to normal samples with HPV- tumors having the highest GLI1 expression, suggesting that increased expression of GLI1 is a potential driver in HPV- HNSCC. Our algorithm for integration of DNA methylation and gene expression can infer biologically significant molecular pathways that may be exploited as therapeutics targets. Our results suggest that therapeutics targeting the Hedgehog pathway may be of benefit in HPV- HNSCC. Similar integrative analysis of high-throughput coupled DNA methylation and expression datasets may yield novel insights into deregulated pathways in other cancers.

Hybrid Modeling of Cell Signaling and Transcriptional Reprogramming and Its Application in C. elegans Development.

Modeling of signal driven transcriptional reprogramming is critical for understanding of organism development, human disease, and cell biology. Many current modeling techniques discount key features of the biological sub-systems when modeling multiscale, organism-level processes. We present a mechanistic hybrid model, GESSA, which integrates a novel pooled probabilistic Boolean network model of cell signaling and a stochastic simulation of transcription and translation responding to a diffusion model of extracellular signals. We apply the model to simulate the well studied cell fate decision process of the vulval precursor cells (VPCs) in C. elegans, using experimentally derived rate constants wherever possible and shared parameters to avoid overfitting. We demonstrate that GESSA recovers (1) the effects of varying scaffold protein concentration on signal strength, (2) amplification of signals in expression, (3) the relative external ligand concentration in a known geometry, and (4) feedback in biochemical networks. We demonstrate that setting model parameters based on wild-type and LIN-12 loss-of-function mutants in C. elegans leads to correct prediction of a wide variety of mutants including partial penetrance of phenotypes. Moreover, the model is relatively insensitive to parameters, retaining the wild-type phenotype for a wide range of cell signaling rate parameters.

An algorithm for identification of regulatory signals in unaligned DNA sequences, its testing and parallel implementation.

We describe an algorithm (IRSA) for identification of common regulatory signals in samples of unaligned DNA sequences. The algorithm was tested on randomly generated sequences of fixed length with implanted signal of length 15 with 4 mutations, and on natural upstream regions of bacterial genes regulated by PurR, ArgR and CRP. Then it was applied to upstream regions of orthologous genes from Escherichia coli and related genomes. Some new palindromic binding and direct repeats signals were identified. Finally we present a parallel version suitable for computers supporting the MPI protocol. This implementation is not strictly bounded by the number of available processors. The computation speed linearly depends on the number of processors.