Splitpea: quantifying protein interaction network rewiring changes due to alternative splicing in cancer.

Protein-protein interactions play an essential role in nearly all biological processes, and it has become increasingly clear that in order to better understand the fundamental processes that underlie disease, we must develop a strong understanding of both their context specificity (e.g., tissue-specificity) as well as their dynamic nature (e.g., how they respond to environmental changes). While network-based approaches have found much initial success in the application of protein-protein interactions (PPIs) towards systems-level explorations of biology, they often overlook the fact that large numbers of proteins undergo alternative splicing. Alternative splicing has not only been shown to diversify protein function through the generation of multiple protein isoforms, but also remodel PPIs and affect a wide range diseases, including cancer. Isoform-specific interactions are not well characterized, so we develop a computational approach that uses domain-domain interactions in concert with differential exon usage data from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression project (GTEx). Using this approach, we can characterize PPIs likely disrupted or possibly even increased due to splicing events for individual TCGA cancer patient samples relative to a matched GTEx normal tissue background.

An activation to memory differentiation trajectory of tumor-infiltrating lymphocytes informs metastatic melanoma outcomes.

There is a need for better classification and understanding of tumor-infiltrating lymphocytes (TILs). Here, we applied advanced functional genomics to interrogate 9,000 human tumors and multiple single-cell sequencing sets using benchmarked T cell states, comprehensive T cell differentiation trajectories, human and mouse vaccine responses, and other human TILs. Compared with other T cell states, enrichment of T memory/resident memory programs was observed across solid tumors. Trajectory analysis of single-cell melanoma CD8+ TILs also identified a high fraction of memory/resident memory-scoring TILs in anti-PD-1 responders, which expanded post therapy. In contrast, TILs scoring highly for early T cell activation, but not exhaustion, associated with non-response. Late/persistent, but not early activation signatures, prognosticate melanoma survival, and co-express with dendritic cell and IFN-γ response programs. These data identify an activation-like state associated to poor response and suggest successful memory conversion, above resuscitation of exhaustion, is an under-appreciated aspect of successful anti-tumoral immunity.

SynNotch-CAR T cells overcome challenges of specificity, heterogeneity, and persistence in treating glioblastoma.

Treatment of solid cancers with chimeric antigen receptor (CAR) T cells is plagued by the lack of ideal target antigens that are both absolutely tumor specific and homogeneously expressed. We show that multi-antigen prime-and-kill recognition circuits provide flexibility and precision to overcome these challenges in the context of glioblastoma. A synNotch receptor that recognizes a specific priming antigen, such as the heterogeneous but tumor-specific glioblastoma neoantigen epidermal growth factor receptor splice variant III (EGFRvIII) or the central nervous system (CNS) tissue-specific antigen myelin oligodendrocyte glycoprotein (MOG), can be used to locally induce expression of a CAR. This enables thorough but controlled tumor cell killing by targeting antigens that are homogeneous but not absolutely tumor specific. Moreover, synNotch-regulated CAR expression averts tonic signaling and exhaustion, maintaining a higher fraction of the T cells in a naïve/stem cell memory state. In immunodeficient mice bearing intracerebral patient-derived xenografts (PDXs) with heterogeneous expression of EGFRvIII, a single intravenous infusion of EGFRvIII synNotch-CAR T cells demonstrated higher antitumor efficacy and T cell durability than conventional constitutively expressed CAR T cells, without off-tumor killing. T cells transduced with a synNotch-CAR circuit primed by the CNS-specific antigen MOG also exhibited precise and potent control of intracerebral PDX without evidence of priming outside of the brain. In summary, by using circuits that integrate recognition of multiple imperfect but complementary antigens, we improve the specificity, completeness, and persistence of T cells directed against glioblastoma, providing a general recognition strategy applicable to other solid tumors.

Discriminatory Power of Combinatorial Antigen Recognition in Cancer T Cell Therapies.

Precise discrimination of tumor from normal tissues remains a major roadblock for therapeutic efficacy of chimeric antigen receptor (CAR) T cells. Here, we perform a comprehensive in silico screen to identify multi-antigen signatures that improve tumor discrimination by CAR T cells engineered to integrate multiple antigen inputs via Boolean logic, e.g., AND and NOT. We screen >2.5 million dual antigens and ∼60 million triple antigens across 33 tumor types and 34 normal tissues. We find that dual antigens significantly outperform the best single clinically investigated CAR targets and confirm key predictions experimentally. Further, we identify antigen triplets that are predicted to show close to ideal tumor-versus-normal tissue discrimination for several tumor types. This work demonstrates the potential of 2- to 3-antigen Boolean logic gates for improving tumor discrimination by CAR T cell therapies. Our predictions are available on an interactive web server resource (antigen.princeton.edu).

Serum cytokine levels in breast cancer patients during neoadjuvant treatment with bevacizumab.

A high concentration of circulating vascular endothelial growth factor (VEGF) in cancer patients is associated with an aggressive tumor phenotype. Here, serum levels of 27 cytokines and blood cell counts were assessed in breast cancer patients receiving neoadjuvant chemotherapy with or without bevacizumab (Bev) in a randomized cohort of 132 patients with non-metastatic HER2-negative tumors. Cytokine levels were determined prior to treatment and at various time-points. The cytotoxic chemotherapy regimen of fluorouracil, epirubicin, and cyclophosphamide (FEC) had a profound impact on both circulating white blood cells and circulating cytokine levels. At the end of FEC treatment, the global decrease in cytokine levels correlated with the drop in white blood cell counts and was significantly greater in the patients of the Bev arm for cytokines, such as VEGF-A, IL-12, IP-10 and IL-10. Among patients who received Bev, those with pathological complete response (pCR) exhibited significantly lower levels of VEGF-A, IFN-γ, TNF-α and IL-4 than patients without pCR. This effect was not observed in the chemotherapy-only arm. Certain circulating cytokine profiles were found to correlate with different immune cell types at the tumor site. For the Bev arm patients, the serum cytokine levels correlated with higher levels of cytotoxic T cells at the end of the therapy regimen, which was indicative of treatment response. The higher response rate for Bev-treated patients and stronger correlations between serum cytokine levels and infiltrating CD8T cells merits further investigation.

Data-driven analysis of immune infiltrate in a large cohort of breast cancer and its association with disease progression, ER activity, and genomic complexity.

The tumor microenvironment is now widely recognized for its role in tumor progression, treatment response, and clinical outcome. The intratumoral immunological landscape, in particular, has been shown to exert both pro-tumorigenic and anti-tumorigenic effects. Identifying immunologically active or silent tumors may be an important indication for administration of therapy, and detecting early infiltration patterns may uncover factors that contribute to early risk. Thus far, direct detailed studies of the cell composition of tumor infiltration have been limited; with some studies giving approximate quantifications using immunohistochemistry and other small studies obtaining detailed measurements by isolating cells from excised tumors and sorting them using flow cytometry. Herein we utilize a machine learning based approach to identify lymphocyte markers with which we can quantify the presence of B cells, cytotoxic T-lymphocytes, T-helper 1, and T-helper 2 cells in any gene expression data set and apply it to studies of breast tissue. By leveraging over 2,100 samples from existing large scale studies, we are able to find an inherent cell heterogeneity in clinically characterized immune infiltrates, a strong link between estrogen receptor activity and infiltration in normal and tumor tissues, changes with genomic complexity, and identify characteristic differences in lymphocyte expression among molecular groupings. With our extendable methodology for capturing cell type specific signal we systematically studied immune infiltration in breast cancer, finding an inverse correlation between beneficial lymphocyte infiltration and estrogen receptor activity in normal breast tissue and reduced infiltration in estrogen receptor negative tumors with high genomic complexity.

Bioinformatics Approaches to Profile the Tumor Microenvironment for Immunotherapeutic Discovery.

In the microenvironment of a malignancy, tumor cells do not exist in isolation, but rather in a diverse ecosystem consisting not only of heterogeneous tumor-cell clones, but also normal cell types such as fibroblasts, vasculature, and an extensive pool of immune cells at numerous possible stages of activation and differentiation. This results in a complex interplay of diverse cellular signaling systems, where the immune cell component is now established to influence cancer progression and therapeutic response. It is experimentally difficult and laborious to comprehensively and systematically profile these distinct cell types from heterogeneous tumor samples in order to capitalize on potential therapeutic and biomarker discoveries. One emerging solution to address this challenge is to computationally extract cell-type specific information directly from bulk tumors. Such in silico approaches are advantageous because they can capture both the cell-type specific profiles and the tissue systems level of cell-cell interactions. Accurately and comprehensively predicting these patterns in tumors is an important challenge to overcome, not least given the success of immunotherapeutic drug treatment of several human cancers. This is especially challenging for subsets of closely related immune cell phenotypes with relatively small gene expression differences, which have critical functional distinctions. Here, we outline the existing and emerging novel bioinformatics strategies that can be used to profile the tumor immune landscape.

IFNγ-Dependent Tissue-Immune Homeostasis Is Co-opted in the Tumor Microenvironment.

Homeostatic programs balance immune protection and self-tolerance. Such mechanisms likely impact autoimmunity and tumor formation, respectively. How homeostasis is maintained and impacts tumor surveillance is unknown. Here, we find that different immune mononuclear phagocytes share a conserved steady-state program during differentiation and entry into healthy tissue. IFNγ is necessary and sufficient to induce this program, revealing a key instructive role. Remarkably, homeostatic and IFNγ-dependent programs enrich across primary human tumors, including melanoma, and stratify survival. Single-cell RNA sequencing (RNA-seq) reveals enrichment of homeostatic modules in monocytes and DCs from human metastatic melanoma. Suppressor-of-cytokine-2 (SOCS2) protein, a conserved program transcript, is expressed by mononuclear phagocytes infiltrating primary melanoma and is induced by IFNγ. SOCS2 limits adaptive anti-tumoral immunity and DC-based priming of T cells in vivo, indicating a critical regulatory role. These findings link immune homeostasis to key determinants of anti-tumoral immunity and escape, revealing co-opting of tissue-specific immune development in the tumor microenvironment.

Lymphocyte Invasion in IC10/Basal-Like Breast Tumors Is Associated with Wild-Type TP53.

UNLABELLED: Lymphocytic infiltration is associated with better prognosis in several epithelial malignancies including breast cancer. The tumor suppressor TP53 is mutated in approximately 30% of breast adenocarcinomas, with varying frequency across molecular subtypes. In this study of 1,420 breast tumors, we tested for interaction between TP53 mutation status and tumor subtype determined by PAM50 and integrative cluster analysis. In integrative cluster 10 (IC10)/basal-like breast cancer, we identify an association between lymphocytic infiltration, determined by an expression score, and retention of wild-type TP53. The expression-derived score agreed with the degree of lymphocytic infiltration assessed by pathologic review, and application of the Nanodissect algorithm was suggestive of this infiltration being primarily of cytotoxic T lymphocytes (CTL). Elevated expression of this CTL signature was associated with longer survival in IC10/Basal-like tumors. These findings identify a new link between the TP53 pathway and the adaptive immune response in estrogen receptor (ER)-negative breast tumors, suggesting a connection between TP53 inactivation and failure of tumor immunosurveillance. IMPLICATIONS: The association of lymphocytic invasion of ER-negative breast tumors with the retention of wild-type TP53 implies a novel protective connection between TP53 function and tumor immunosurveillance.

Sets2Networks: network inference from repeated observations of sets.

BACKGROUND: The skeleton of complex systems can be represented as networks where vertices represent entities, and edges represent the relations between these entities. Often it is impossible, or expensive, to determine the network structure by experimental validation of the binary interactions between every vertex pair. It is usually more practical to infer the network from surrogate observations. Network inference is the process by which an underlying network of relations between entities is determined from indirect evidence. While many algorithms have been developed to infer networks from quantitative data, less attention has been paid to methods which infer networks from repeated co-occurrence of entities in related sets. This type of data is ubiquitous in the field of systems biology and in other areas of complex systems research. Hence, such methods would be of great utility and value. RESULTS: Here we present a general method for network inference from repeated observations of sets of related entities. Given experimental observations of such sets, we infer the underlying network connecting these entities by generating an ensemble of networks consistent with the data. The frequency of occurrence of a given link throughout this ensemble is interpreted as the probability that the link is present in the underlying real network conditioned on the data. Exponential random graphs are used to generate and sample the ensemble of consistent networks, and we take an algorithmic approach to numerically execute the inference method. The effectiveness of the method is demonstrated on synthetic data before employing this inference approach to problems in systems biology and systems pharmacology, as well as to construct a co-authorship collaboration network. We predict direct protein-protein interactions from high-throughput mass-spectrometry proteomics, integrate data from Chip-seq and loss-of-function/gain-of-function followed by expression data to infer a network of associations between pluripotency regulators, extract a network that connects 53 cancer drugs to each other and to 34 severe adverse events by mining the FDA's Adverse Events Reporting Systems (AERS), and construct a co-authorship network that connects Mount Sinai School of Medicine investigators. The predicted networks and online software to create networks from entity-set libraries are provided online at http://www.maayanlab.net/S2N. CONCLUSIONS: The network inference method presented here can be applied to resolve different types of networks in current systems biology and systems pharmacology as well as in other fields of research.

Genes2FANs: connecting genes through functional association networks.

BACKGROUND: Protein-protein, cell signaling, metabolic, and transcriptional interaction networks are useful for identifying connections between lists of experimentally identified genes/proteins. However, besides physical or co-expression interactions there are many ways in which pairs of genes, or their protein products, can be associated. By systematically incorporating knowledge on shared properties of genes from diverse sources to build functional association networks (FANs), researchers may be able to identify additional functional interactions between groups of genes that are not readily apparent. RESULTS: Genes2FANs is a web based tool and a database that utilizes 14 carefully constructed FANs and a large-scale protein-protein interaction (PPI) network to build subnetworks that connect lists of human and mouse genes. The FANs are created from mammalian gene set libraries where mouse genes are converted to their human orthologs. The tool takes as input a list of human or mouse Entrez gene symbols to produce a subnetwork and a ranked list of intermediate genes that are used to connect the query input list. In addition, users can enter any PubMed search term and then the system automatically converts the returned results to gene lists using GeneRIF. This gene list is then used as input to generate a subnetwork from the user's PubMed query. As a case study, we applied Genes2FANs to connect disease genes from 90 well-studied disorders. We find an inverse correlation between the counts of links connecting disease genes through PPI and links connecting diseases genes through FANs, separating diseases into two categories. CONCLUSIONS: Genes2FANs is a useful tool for interpreting the relationships between gene/protein lists in the context of their various functions and networks. Combining functional association interactions with physical PPIs can be useful for revealing new biology and help form hypotheses for further experimentation. Our finding that disease genes in many cancers are mostly connected through PPIs whereas other complex diseases, such as autism and type-2 diabetes, are mostly connected through FANs without PPIs, can guide better strategies for disease gene discovery. Genes2FANs is available at: http://actin.pharm.mssm.edu/genes2FANs.