Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Spatial analysis of human lung cancer reveals organized immune hubs enriched for stem-like CD8 T cells and associated with immunotherapy response

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially-localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens, and found that they were associated with beneficial responses to PD-1-blockade. Immunity hubs were enriched for many interferon-stimulated genes, T cells in multiple differentiation states, and CXCL9/10/11 + macrophages that preferentially interact with CD8 T cells. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcomes, distinct from mature tertiary lymphoid structures, and enriched for stem-like TCF7+PD-1+ CD8 T cells and activated CCR7 + LAMP3 + dendritic cells, as well as chemokines that organize these cells. These results elucidate the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Peripheral immune cell reactivity and neural response to reward in patients with depression and anhedonia

Increased levels of peripheral cytokines have been previously associated with depression in preclinical and clinical research. Although the precise nature of peripheral immune dysfunction in depression remains unclear, evidence from animal studies points towards a dysregulated response of peripheral leukocytes as a risk factor for stress susceptibility. This study examined dynamic release of inflammatory blood factors from peripheral blood mononuclear cells (PBMC) in depressed patients and associations with neural and behavioral measures of reward processing. Thirty unmedicated patients meeting criteria for unipolar depressive disorder and 21 healthy control volunteers were enrolled. PBMCs were isolated from whole blood and stimulated ex vivo with lipopolysaccharide (LPS). Olink multiplex assay was used to analyze a large panel of inflammatory proteins. Participants completed functional magnetic resonance imaging with an incentive flanker task to probe neural responses to reward anticipation, as well as clinical measures of anhedonia and pleasure including the Temporal Experience of Pleasure Scale (TEPS) and the Snaith-Hamilton Pleasure Scale (SHAPS). LPS stimulation revealed larger increases in immune factors in depressed compared to healthy subjects using an aggregate immune score (t₄₉= 2.83, p = 0.007). Higher peripheral immune score was associated with reduced neural responses to reward anticipation within the ventral striatum (VS) (r = -0.39, p = 0.01), and with reduced anticipation of pleasure as measured with the TEPS anticipatory sub-score (r = -0.318, p = 0.023). Our study provides new evidence suggesting that dynamic hyper-reactivity of peripheral leukocytes in depressed patients is associated with blunted activation of the brain reward system and lower subjective anticipation of pleasure.

Systems Analysis Implicates WAVE2 Complex in the Pathogenesis of Developmental Left-Sided Obstructive Heart Defects

•

Combining CHD phenotype–driven gene set enrichment and CRISPR knockdown screening in zebrafish is an effective approach to identifying novel CHD genes.

•

Mutations affecting genes coding for the WAVE2 protein complex and small GTPase-mediated signaling are associated with LVOTO lesions.

•

WAVE2 complex genes brk1, nckap1, and wasf2 and regulators of small GTPase signaling cul3a and racgap1 are critical to zebrafish heart development.

Genetic variants are the primary driver of congenital heart disease (CHD) pathogenesis. However, our ability to identify causative variants is limited. To identify causal CHD genes that are associated with specific molecular functions, the study used prior knowledge to filter de novo variants from 2,881 probands with sporadic severe CHD. This approach enabled the authors to identify an association between left ventricular outflow tract obstruction lesions and genes associated with the WAVE2 complex and regulation of small GTPase-mediated signal transduction. Using CRISPR zebrafish knockdowns, the study confirmed that WAVE2 complex proteins brk1, nckap1, and wasf2 and the regulators of small GTPase signaling cul3a and racgap1 are critical to cardiac development.

Length of gestation and birth weight are associated with indices of combined kidney biomarkers in early childhood

Infants born prematurely or with low birth weights are more susceptible to kidney dysfunction throughout their lives. Multiple proteins measured in urine are noninvasive biomarkers of subclinical kidney damage, but few studies have examined the joint effects of multiple biomarkers. We conducted an exploratory study of 103 children in the Programing Research in Obesity, Growth, Environment, and Social Stressors (PROGRESS) longitudinal birth cohort, and measured nine proteins selected a priori in banked spot urine samples collected at ages 4-6. The goal of our study was to explore the combined effects of kidney damage biomarkers previously associated with birth outcomes. To do this, we generated kidney biomarker indices using weighted quantile sum regression and assessed associations with length of gestation or birth weight. A decile increase in each kidney biomarker index was associated with 2-day shorter gestations (β = -2.0, 95% CI: -3.2, -0.9) and 59-gram lower birth weights (β = -58.5, 95% CI: -98.3, -18.7), respectively. Weights highlighting the contributions showed neutrophil gelatinase-associated lipocalin (NGAL) (60%) and osteopontin (19%) contributed most to the index derived for gestational age. NGAL (66%) and beta-2-microglobulin (10%) contributed most to the index derived for birth weight. Joint analyses of multiple kidney biomarkers can provide integrated measures of kidney dysfunction and improved statistical assessments compared to biomarkers assessed individually. Additionally, shorter gestations and lower birth weights may contribute to subclinical kidney damage measurable in childhood.

Single-cell immune landscape of human atherosclerotic plaques

Atherosclerosis is driven by multifaceted contributions of the immune system within the circulation and at vascular focal sites. However, specific characteristics of dysregulated immune cells within atherosclerotic lesions that lead to clinical events such as ischemic stroke or myocardial infarction are poorly understood. Here, using single-cell proteomic and transcriptomic analyses, we uncovered distinct features of both T cells and macrophages in carotid artery plaques of patients with clinically symptomatic disease (recent stroke or transient ischemic attack) compared to asymptomatic disease (no recent stroke). Plaques from symptomatic patients were characterized by a distinct subset of CD4+ T cells and by T cells that were activated and differentiated. Moreover, some T cell subsets in these plaques presented markers of T cell exhaustion. Additionally, macrophages from these plaques contained alternatively activated phenotypes, including subsets associated with plaque vulnerability. In plaques from asymptomatic patients, T cells and macrophages were activated and displayed evidence of interleukin-1β signaling. The identification of specific features of innate and adaptive immune cells in plaques that are associated with cerebrovascular events may enable the design of more precisely tailored cardiovascular immunotherapies.

Clustergrammer, a web-based heatmap visualization and analysis tool for high-dimensional biological data

Most tools developed to visualize hierarchically clustered heatmaps generate static images. Clustergrammer is a web-based visualization tool with interactive features such as: zooming, panning, filtering, reordering, sharing, performing enrichment analysis, and providing dynamic gene annotations. Clustergrammer can be used to generate shareable interactive visualizations by uploading a data table to a web-site, or by embedding Clustergrammer in Jupyter Notebooks. The Clustergrammer core libraries can also be used as a toolkit by developers to generate visualizations within their own applications. Clustergrammer is demonstrated using gene expression data from the cancer cell line encyclopedia (CCLE), original post-translational modification data collected from lung cancer cells lines by a mass spectrometry approach, and original cytometry by time of flight (CyTOF) single-cell proteomics data from blood. Clustergrammer enables producing interactive web based visualizations for the analysis of diverse biological data.

Differential cytokine contributions of perivascular haematopoietic stem cell niches

Arterioles and sinusoids of the bone marrow (BM) are accompanied by stromal cells that express nerve/glial antigen 2 (NG2) and leptin receptor (LepR), and constitute specialized niches that regulate quiescence and proliferation of haematopoietic stem cells (HSCs). However, how niche cells differentially regulate HSC functions remains unknown. Here, we show that the effects of cytokines regulating HSC functions are dependent on the producing cell sources. Deletion of chemokine C-X-C motif ligand 12 (Cxcl12) or stem cell factor (Scf) from all perivascular cells marked by nestin-GFP dramatically depleted BM HSCs. Selective Cxcl12 deletion from arteriolar NG2+ cells, but not from sinusoidal LepR+ cells, caused HSC reductions and altered HSC localization in BM. By contrast, deletion of Scf in LepR+ cells, but not NG2+ cells, led to reductions in BM HSC numbers. These results uncover distinct contributions of cytokines derived from perivascular cells in separate vascular niches to HSC maintenance.

Extraction and analysis of signatures from the Gene Expression Omnibus by the crowd

Gene expression data are accumulating exponentially in public repositories. Reanalysis and integration of themed collections from these studies may provide new insights, but requires further human curation. Here we report a crowdsourcing project to annotate and reanalyse a large number of gene expression profiles from Gene Expression Omnibus (GEO). Through a massive open online course on Coursera, over 70 participants from over 25 countries identify and annotate 2,460 single-gene perturbation signatures, 839 disease versus normal signatures, and 906 drug perturbation signatures. All these signatures are unique and are manually validated for quality. Global analysis of these signatures confirms known associations and identifies novel associations between genes, diseases and drugs. The manually curated signatures are used as a training set to develop classifiers for extracting similar signatures from the entire GEO repository. We develop a web portal to serve these signatures for query, download and visualization.

L1000CDS2: LINCS L1000 characteristic direction signatures search engine

The library of integrated network-based cellular signatures (LINCS) L1000 data set currently comprises of over a million gene expression profiles of chemically perturbed human cell lines. Through unique several intrinsic and extrinsic benchmarking schemes, we demonstrate that processing the L1000 data with the characteristic direction (CD) method significantly improves signal to noise compared with the MODZ method currently used to compute L1000 signatures. The CD processed L1000 signatures are served through a state-of-the-art web-based search engine application called L1000CDS². The L1000CDS² search engine provides prioritization of thousands of small-molecule signatures, and their pairwise combinations, predicted to either mimic or reverse an input gene expression signature using two methods. The L1000CDS² search engine also predicts drug targets for all the small molecules profiled by the L1000 assay that we processed. Targets are predicted by computing the cosine similarity between the L1000 small-molecule signatures and a large collection of signatures extracted from the gene expression omnibus (GEO) for single-gene perturbations in mammalian cells. We applied L1000CDS² to prioritize small molecules that are predicted to reverse expression in 670 disease signatures also extracted from GEO, and prioritized small molecules that can mimic expression of 22 endogenous ligand signatures profiled by the L1000 assay. As a case study, to further demonstrate the utility of L1000CDS², we collected expression signatures from human cells infected with Ebola virus at 30, 60 and 120 min. Querying these signatures with L1000CDS² we identified kenpaullone, a GSK3B/CDK2 inhibitor that we show, in subsequent experiments, has a dose-dependent efficacy in inhibiting Ebola infection in vitro without causing cellular toxicity in human cell lines. In summary, the L1000CDS² tool can be applied in many biological and biomedical settings, while improving the extraction of knowledge from the LINCS L1000 resource.

Enrichr: a comprehensive gene set enrichment analysis web server 2016 update

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

Enrichr: a comprehensive gene set enrichment analysis web server 2016 update

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

Genomics, epigenomics, transcriptomics, proteomics and metabolomics efforts rapidly generate a plethora of data on the activity and levels of biomolecules within mammalian cells. At the same time, curation projects that organize knowledge from the biomedical literature into online databases are expanding. Hence, there is a wealth of information about genes, proteins and their associations, with an urgent need for data integration to achieve better knowledge extraction and data reuse. For this purpose, we developed the Harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins from over 70 major online resources. We extracted, abstracted and organized data into ∼72 million functional associations between genes/proteins and their attributes. Such attributes could be physical relationships with other biomolecules, expression in cell lines and tissues, genetic associations with knockout mouse or human phenotypes, or changes in expression after drug treatment. We stored these associations in a relational database along with rich metadata for the genes/proteins, their attributes and the original resources. The freely available Harmonizome web portal provides a graphical user interface, a web service and a mobile app for querying, browsing and downloading all of the collected data. To demonstrate the utility of the Harmonizome, we computed and visualized gene–gene and attribute–attribute similarity networks, and through unsupervised clustering, identified many unexpected relationships by combining pairs of datasets such as the association between kinase perturbations and disease signatures. We also applied supervised machine learning methods to predict novel substrates for kinases, endogenous ligands for G-protein coupled receptors, mouse phenotypes for knockout genes, and classified unannotated transmembrane proteins for likelihood of being ion channels. The Harmonizome is a comprehensive resource of knowledge about genes and proteins, and as such, it enables researchers to discover novel relationships between biological entities, as well as form novel data-driven hypotheses for experimental validation.

Database URL: http://amp.pharm.mssm.edu/Harmonizome.

Enrichr: a comprehensive gene set enrichment analysis web server 2016 update

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

Drug/Cell-line Browser: interactive canvas visualization of cancer drug/cell-line viability assay datasets

SUMMARY

Recently, several high profile studies collected cell viability data from panels of cancer cell lines treated with many drugs applied at different concentrations. Such drug sensitivity data for cancer cell lines provide suggestive treatments for different types and subtypes of cancer. Visualization of these datasets can reveal patterns that may not be obvious by examining the data without such efforts. Here we introduce Drug/Cell-line Browser (DCB), an online interactive HTML5 data visualization tool for interacting with three of the recently published datasets of cancer cell lines/drug-viability studies. DCB uses clustering and canvas visualization of the drugs and the cell lines, as well as a bar graph that summarizes drug effectiveness for the tissue of origin or the cancer subtypes for single or multiple drugs. DCB can help in understanding drug response patterns and prioritizing drug/cancer cell line interactions by tissue of origin or cancer subtype.

AVAILABILITY AND IMPLEMENTATION

DCB is an open source Web-based tool that is freely available at: http://www.maayanlab.net/LINCS/DCB CONTACT: avi.maayan@mssm.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

LINCS Canvas Browser: interactive web app to query, browse and interrogate LINCS L1000 gene expression signatures

For the Library of Integrated Network-based Cellular Signatures (LINCS) project many gene expression signatures using the L1000 technology have been produced. The L1000 technology is a cost-effective method to profile gene expression in large scale. LINCS Canvas Browser (LCB) is an interactive HTML5 web-based software application that facilitates querying, browsing and interrogating many of the currently available LINCS L1000 data. LCB implements two compacted layered canvases, one to visualize clustered L1000 expression data, and the other to display enrichment analysis results using 30 different gene set libraries. Clicking on an experimental condition highlights gene-sets enriched for the differentially expressed genes from the selected experiment. A search interface allows users to input gene lists and query them against over 100 000 conditions to find the top matching experiments. The tool integrates many resources for an unprecedented potential for new discoveries in systems biology and systems pharmacology. The LCB application is available at http://www.maayanlab.net/LINCS/LCB. Customized versions will be made part of the http://lincscloud.org and http://lincs.hms.harvard.edu websites.