Regulatory Network of PD1 Signaling Is Associated with Prognosis in Glioblastoma Multiforme

Adjustment of spurious correlations in co-expression measurements from RNA-Sequencing data

MOTIVATION

Gene co-expression measurements are widely used in computational biology to identify coordinated expression patterns across a group of samples. Coordinated expression of genes may indicate that they are controlled by the same transcriptional regulatory program, or involved in common biological processes. Gene co-expression is generally estimated from RNA-Sequencing data, which are commonly normalized to remove technical variability. Here, we demonstrate that certain normalization methods, in particular quantile-based methods, can introduce false-positive associations between genes. These false-positive associations can consequently hamper downstream co-expression network analysis. Quantile-based normalization can, however, be extremely powerful. In particular, when preprocessing large-scale heterogeneous data, quantile-based normalization methods such as smooth quantile normalization can be applied to remove technical variability while maintaining global differences in expression for samples with different biological attributes.

RESULTS

We developed SNAIL (Smooth-quantile Normalization Adaptation for the Inference of co-expression Links), a normalization method based on smooth quantile normalization specifically designed for modeling of co-expression measurements. We show that SNAIL avoids formation of false-positive associations in co-expression as well as in downstream network analyses. Using SNAIL, one can avoid arbitrary gene filtering and retain associations to genes that only express in small subgroups of samples. This highlights the method's potential future impact on network modeling and other association-based approaches in large-scale heterogeneous data.

AVAILABILITY AND IMPLEMENTATION

The implementation of the SNAIL algorithm and code to reproduce the analyses described in this work can be found in the GitHub repository https://github.com/kuijjerlab/PySNAIL.

Heterogeneity in the gene regulatory landscape of leiomyosarcoma

Characterizing inter-tumor heterogeneity is crucial for selecting suitable cancer therapy, as the presence of diverse molecular subgroups of patients can be associated with disease outcome or response to treatment. While cancer subtypes are often characterized by differences in gene expression, the mechanisms driving these differences are generally unknown. We set out to model the regulatory mechanisms driving sarcoma heterogeneity based on patient-specific, genome-wide gene regulatory networks. We developed a new computational framework, PORCUPINE, which combines knowledge on biological pathways with permutation-based network analysis to identify pathways that exhibit significant regulatory heterogeneity across a patient population. We applied PORCUPINE to patient-specific leiomyosarcoma networks modeled on data from The Cancer Genome Atlas and validated our results in an independent dataset from the German Cancer Research Center. PORCUPINE identified 37 heterogeneously regulated pathways, including pathways representing potential targets for treatment of subgroups of leiomyosarcoma patients, such as FGFR and CTLA4 inhibitory signaling. We validated the detected regulatory heterogeneity through analysis of networks and chromatin states in leiomyosarcoma cell lines. We showed that the heterogeneity identified with PORCUPINE is not associated with methylation profiles or clinical features, thereby suggesting an independent mechanism of patient heterogeneity driven by the complex landscape of gene regulatory interactions.

Glioblastoma survival is associated with distinct proteomic alteration signatures post chemoirradiation in a large-scale proteomic panel

Background

Glioblastomas (GBM) are rapidly progressive, nearly uniformly fatal brain tumors. Proteomic analysis represents an opportunity for noninvasive GBM classification and biological understanding of treatment response.

Purpose

We analyzed differential proteomic expression pre vs. post completion of concurrent chemoirradiation (CRT) in patient serum samples to explore proteomic alterations and classify GBM by integrating clinical and proteomic parameters.

Materials and methods

82 patients with GBM were clinically annotated and serum samples obtained pre- and post-CRT. Serum samples were then screened using the aptamer-based SOMAScan® proteomic assay. Significant traits from uni- and multivariate Cox models for overall survival (OS) were designated independent prognostic factors and principal component analysis (PCA) was carried out. Differential expression of protein signals was calculated using paired t-tests, with KOBAS used to identify associated KEGG pathways. GSEA pre-ranked analysis was employed on the overall list of differentially expressed proteins (DEPs) against the MSigDB Hallmark, GO Biological Process, and Reactome databases with weighted gene correlation network analysis (WGCNA) and Enrichr used to validate pathway hits internally.

Results

3 clinical clusters of patients with differential survival were identified. 389 significantly DEPs pre vs. post-treatment were identified, including 284 upregulated and 105 downregulated, representing several pathways relevant to cancer metabolism and progression. The lowest survival group (median OS 13.2 months) was associated with DEPs affiliated with proliferative pathways and exhibiting distinct oppositional response including with respect to radiation therapy related pathways, as compared to better-performing groups (intermediate, median OS 22.4 months; highest, median OS 28.7 months). Opposite signaling patterns across multiple analyses in several pathways (notably fatty acid metabolism, NOTCH, TNFα via NF-κB, Myc target V1 signaling, UV response, unfolded protein response, peroxisome, and interferon response) were distinct between clinical survival groups and supported by WGCNA. 23 proteins were statistically signficant for OS with 5 (NETO2, CST7, SEMA6D, CBLN4, NPS) supported by KM.

Conclusion

Distinct proteomic alterations with hallmarks of cancer, including progression, resistance, stemness, and invasion, were identified in serum samples obtained from GBM patients pre vs. post CRT and corresponded with clinical survival. The proteome can potentially be employed for glioma classification and biological interrogation of cancer pathways.

The Network Zoo: a multilingual package for the inference and analysis of gene regulatory networks

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods.

Network-based approaches for modeling disease regulation and progression

Molecular interaction networks lay the foundation for studying how biological functions are controlled by the complex interplay of genes and proteins. Investigating perturbed processes using biological networks has been instrumental in uncovering mechanisms that underlie complex disease phenotypes. Rapid advances in omics technologies have prompted the generation of high-throughput datasets, enabling large-scale, network-based analyses. Consequently, various modeling techniques, including network enrichment, differential network extraction, and network inference, have proven to be useful for gaining new mechanistic insights. We provide an overview of recent network-based methods and their core ideas to facilitate the discovery of disease modules or candidate mechanisms. Knowledge generated from these computational efforts will benefit biomedical research, especially drug development and precision medicine. We further discuss current challenges and provide perspectives in the field, highlighting the need for more integrative and dynamic network approaches to model disease development and progression.

Latest updates on cellular and molecular biomarkers of gliomas

Gliomas are the most common central nervous system malignancies, compromising almost 80% of all brain tumors and is associated with significant mortality. The classification of gliomas has shifted from basic histological perspective to one that is based on molecular biomarkers. Treatment of this type of tumors consists currently of surgery, chemotherapy and radiation therapy. During the past years, there was a limited development of effective glioma diagnostics and therapeutics due to multiple factors including the presence of blood-brain barrier and the heterogeneity of this type of tumors. Currently, it is necessary to highlight the advantage of molecular diagnosis of gliomas to develop patient targeted therapies based on multiple oncogenic pathway. In this review, we will evaluate the development of cellular and molecular biomarkers for the diagnosis of gliomas and the impact of these diagnostic tools for better tailored and targeted therapies.

Sex differences in gene regulatory networks during mid-gestational brain development

Neurodevelopmental disorders differ considerably between males and females, and fetal brain development is one of the most critical periods to determine risk for these disorders. Transcriptomic studies comparing male and female fetal brain have demonstrated that the highest difference in gene expression occurs in sex chromosomes, but several autossomal genes also demonstrate a slight difference that has not been yet explored. In order to investigate biological pathways underlying fetal brain sex differences, we applied medicine network principles using integrative methods such as co-expression networks (CEMiTool) and regulatory networks (netZoo). The pattern of gene expression from genes in the same pathway tend to reflect biologically relevant phenomena. In this study, network analysis of fetal brain expression reveals regulatory differences between males and females. Integrating two different bioinformatics tools, our results suggest that biological processes such as cell cycle, cell differentiation, energy metabolism and extracellular matrix organization are consistently sex-biased. MSET analysis demonstrates that these differences are relevant to neurodevelopmental disorders, including autism.

gpuZoo: Cost-effective estimation of gene regulatory networks using the Graphics Processing Unit

Gene regulatory network inference allows for the modeling of genome-scale regulatory processes that are altered during development, in disease, and in response to perturbations. Our group has developed a collection of tools to model various regulatory processes, including transcriptional (PANDA, SPIDER) and post-transcriptional (PUMA) gene regulation, as well as gene regulation in individual samples (LIONESS). These methods work by postulating a network structure and then optimizing that structure to be consistent with multiple lines of biological evidence through repeated operations on data matrices. Although our methods are widely used, the corresponding computational complexity, and the associated costs and run times, do limit some applications. To improve the cost/time performance of these algorithms, we developed gpuZoo which implements GPU-accelerated calculations, dramatically improving the performance of these algorithms. The runtime of the gpuZoo implementation in MATLAB and Python is up to 61 times faster and 28 times less expensive than multi-core CPU implementation of the same methods. gpuZoo is available in MATLAB through the netZooM package https://github.com/netZoo/netZooM and in Python through the netZooPy package https://github.com/netZoo/netZooPy.