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.

Joint embedding of biological networks for cross-species functional alignment.

MOTIVATION: Model organisms are widely used to better understand the molecular causes of human disease. While sequence similarity greatly aids this cross-species transfer, sequence similarity does not imply functional similarity, and thus, several current approaches incorporate protein-protein interactions to help map findings between species. Existing transfer methods either formulate the alignment problem as a matching problem which pits network features against known orthology, or more recently, as a joint embedding problem. RESULTS: We propose a novel state-of-the-art joint embedding solution: Embeddings to Network Alignment (ETNA). ETNA generates individual network embeddings based on network topological structure and then uses a Natural Language Processing-inspired cross-training approach to align the two embeddings using sequence-based orthologs. The final embedding preserves both within and between species gene functional relationships, and we demonstrate that it captures both pairwise and group functional relevance. In addition, ETNA's embeddings can be used to transfer genetic interactions across species and identify phenotypic alignments, laying the groundwork for potential opportunities for drug repurposing and translational studies. AVAILABILITY AND IMPLEMENTATION: https://github.com/ylaboratory/ETNA.

Patterns and rates of exonic de novo mutations in autism spectrum disorders.

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.

Network2Canvas: network visualization on a canvas with enrichment analysis.

MOTIVATION: Networks are vital to computational systems biology research, but visualizing them is a challenge. For networks larger than ∼100 nodes and ∼200 links, ball-and-stick diagrams fail to convey much information. To address this, we developed Network2Canvas (N2C), a web application that provides an alternative way to view networks. N2C visualizes networks by placing nodes on a square toroidal canvas. The network nodes are clustered on the canvas using simulated annealing to maximize local connections where a node's brightness is made proportional to its local fitness. The interactive canvas is implemented in HyperText Markup Language (HTML)5 with the JavaScript library Data-Driven Documents (D3). We applied N2C to visualize 30 canvases made from human and mouse gene-set libraries and 6 canvases made from the Food and Drug Administration (FDA)-approved drug-set libraries. Given lists of genes or drugs, enriched terms are highlighted on the canvases, and their degree of clustering is computed. Because N2C produces visual patterns of enriched terms on canvases, a trained eye can detect signatures instantly. In summary, N2C provides a new flexible method to visualize large networks and can be used to perform and visualize gene-set and drug-set enrichment analyses. AVAILABILITY: N2C is freely available at http://www.maayanlab.net/N2C and is open source. CONTACT: avi.maayan@mssm.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Disentangling associations between complex traits and cell types with seismic.

Integrating single-cell RNA sequencing (scRNA-seq) with Genome-Wide Association Studies (GWAS) can help reveal GWAS-associated cell types, furthering our understanding of the cell-type-specific biological processes underlying complex traits and disease. However, current methods have technical limitations that hinder them from making systematic, scalable, interpretable disease-cell-type associations. In order to rapidly and accurately pinpoint associations, we develop a novel framework, seismic, which characterizes cell types using a new specificity score. We compare seismic with alternative methods across over 1,000 cell type characterizations at different granularities and 28 traits, demonstrating that seismic both corroborates findings and identifies trait-relevant cell groups which are not apparent through other methodologies. Furthermore, as part of the seismic framework, the specific genes driving cell type-trait associations can easily be accessed and analyzed, enabling further biological insights. The advantages of seismic are particularly salient in neurodegenerative diseases such as Parkinson's and Alzheimer's, where disease pathology has not only cell-specific manifestations, but also brain region-specific differences. Interestingly, a case study of Alzheimer's disease reveals the importance of considering GWAS endpoints, as studies relying on clinical diagnoses consistently identify microglial associations, while GWAS with a tau biomarker endpoint reveals neuronal associations. In general, seismic is a computationally efficient, powerful, and interpretable approach for identifying associations between complex traits and cell type-specific expression.

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.

FNV: light-weight flash-based network and pathway viewer.

MOTIVATION: Network diagrams are commonly used to visualize biochemical pathways by displaying the relationships between genes, proteins, mRNAs, microRNAs, metabolites, regulatory DNA elements, diseases, viruses and drugs. While there are several currently available web-based pathway viewers, there is still room for improvement. To this end, we have developed a flash-based network viewer (FNV) for the visualization of small to moderately sized biological networks and pathways. SUMMARY: Written in Adobe ActionScript 3.0, the viewer accepts simple Extensible Markup Language (XML) formatted input files to display pathways in vector graphics on any web-page providing flexible layout options, interactivity with the user through tool tips, hyperlinks and the ability to rearrange nodes on the screen. FNV was utilized as a component in several web-based systems, namely Genes2Networks, Lists2Networks, KEA, ChEA and PathwayGenerator. In addition, FVN can be used to embed pathways inside pdf files for the communication of pathways in soft publication materials. AVAILABILITY: FNV is available for use and download along with the supporting documentation and sample networks at http://www.maayanlab.net/FNV. CONTACT: avi.maayan@mssm.edu.

Recovering protein-protein and domain-domain interactions from aggregation of IP-MS proteomics of coregulator complexes.

Coregulator proteins (CoRegs) are part of multi-protein complexes that transiently assemble with transcription factors and chromatin modifiers to regulate gene expression. In this study we analyzed data from 3,290 immuno-precipitations (IP) followed by mass spectrometry (MS) applied to human cell lines aimed at identifying CoRegs complexes. Using the semi-quantitative spectral counts, we scored binary protein-protein and domain-domain associations with several equations. Unlike previous applications, our methods scored prey-prey protein-protein interactions regardless of the baits used. We also predicted domain-domain interactions underlying predicted protein-protein interactions. The quality of predicted protein-protein and domain-domain interactions was evaluated using known binary interactions from the literature, whereas one protein-protein interaction, between STRN and CTTNBP2NL, was validated experimentally; and one domain-domain interaction, between the HEAT domain of PPP2R1A and the Pkinase domain of STK25, was validated using molecular docking simulations. The scoring schemes presented here recovered known, and predicted many new, complexes, protein-protein, and domain-domain interactions. The networks that resulted from the predictions are provided as a web-based interactive application at http://maayanlab.net/HT-IP-MS-2-PPI-DDI/.

Current progress and open challenges for applying deep learning across the biosciences.

Deep Learning (DL) has recently enabled unprecedented advances in one of the grand challenges in computational biology: the half-century-old problem of protein structure prediction. In this paper we discuss recent advances, limitations, and future perspectives of DL on five broad areas: protein structure prediction, protein function prediction, genome engineering, systems biology and data integration, and phylogenetic inference. We discuss each application area and cover the main bottlenecks of DL approaches, such as training data, problem scope, and the ability to leverage existing DL architectures in new contexts. To conclude, we provide a summary of the subject-specific and general challenges for DL across the biosciences.

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.

JCoDA: a tool for detecting evolutionary selection.

BACKGROUND: The incorporation of annotated sequence information from multiple related species in commonly used databases (Ensembl, Flybase, Saccharomyces Genome Database, Wormbase, etc.) has increased dramatically over the last few years. This influx of information has provided a considerable amount of raw material for evaluation of evolutionary relationships. To aid in the process, we have developed JCoDA (Java Codon Delimited Alignment) as a simple-to-use visualization tool for the detection of site specific and regional positive/negative evolutionary selection amongst homologous coding sequences. RESULTS: JCoDA accepts user-inputted unaligned or pre-aligned coding sequences, performs a codon-delimited alignment using ClustalW, and determines the dN/dS calculations using PAML (Phylogenetic Analysis Using Maximum Likelihood, yn00 and codeml) in order to identify regions and sites under evolutionary selection. The JCoDA package includes a graphical interface for Phylip (Phylogeny Inference Package) to generate phylogenetic trees, manages formatting of all required file types, and streamlines passage of information between underlying programs. The raw data are output to user configurable graphs with sliding window options for straightforward visualization of pairwise or gene family comparisons. Additionally, codon-delimited alignments are output in a variety of common formats and all dN/dS calculations can be output in comma-separated value (CSV) format for downstream analysis. To illustrate the types of analyses that are facilitated by JCoDA, we have taken advantage of the well studied sex determination pathway in nematodes as well as the extensive sequence information available to identify genes under positive selection, examples of regional positive selection, and differences in selection based on the role of genes in the sex determination pathway. CONCLUSIONS: JCoDA is a configurable, open source, user-friendly visualization tool for performing evolutionary analysis on homologous coding sequences. JCoDA can be used to rapidly screen for genes and regions of genes under selection using PAML. It can be freely downloaded at http://www.tcnj.edu/~nayaklab/jcoda.

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).

Selective Neuronal Vulnerability in Alzheimer's Disease: A Network-Based Analysis.

A major obstacle to treating Alzheimer's disease (AD) is our lack of understanding of the molecular mechanisms underlying selective neuronal vulnerability, a key characteristic of the disease. Here, we present a framework integrating high-quality neuron-type-specific molecular profiles across the lifetime of the healthy mouse, which we generated using bacTRAP, with postmortem human functional genomics and quantitative genetics data. We demonstrate human-mouse conservation of cellular taxonomy at the molecular level for neurons vulnerable and resistant in AD, identify specific genes and pathways associated with AD neuropathology, and pinpoint a specific functional gene module underlying selective vulnerability, enriched in processes associated with axonal remodeling, and affected by amyloid accumulation and aging. We have made all cell-type-specific profiles and functional networks available at http://alz.princeton.edu. Overall, our study provides a molecular framework for understanding the complex interplay between Aß, aging, and neurodegeneration within the most vulnerable neurons in AD.

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.

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.

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.

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.

ESCAPE: database for integrating high-content published data collected from human and mouse embryonic stem cells.

High content studies that profile mouse and human embryonic stem cells (m/hESCs) using various genome-wide technologies such as transcriptomics and proteomics are constantly being published. However, efforts to integrate such data to obtain a global view of the molecular circuitry in m/hESCs are lagging behind. Here, we present an m/hESC-centered database called Embryonic Stem Cell Atlas from Pluripotency Evidence integrating data from many recent diverse high-throughput studies including chromatin immunoprecipitation followed by deep sequencing, genome-wide inhibitory RNA screens, gene expression microarrays or RNA-seq after knockdown (KD) or overexpression of critical factors, immunoprecipitation followed by mass spectrometry proteomics and phosphoproteomics. The database provides web-based interactive search and visualization tools that can be used to build subnetworks and to identify known and novel regulatory interactions across various regulatory layers. The web-interface also includes tools to predict the effects of combinatorial KDs by additive effects controlled by sliders, or through simulation software implemented in MATLAB. Overall, the Embryonic Stem Cell Atlas from Pluripotency Evidence database is a comprehensive resource for the stem cell systems biology community. Database URL: http://www.maayanlab.net/ESCAPE