Dynamic CD8+ T cell responses to cancer immunotherapy in human regional lymph nodes are disrupted in metastatic lymph nodes.

CD8+ T cell responses are critical for anti-tumor immunity. While extensively profiled in the tumor microenvironment, recent studies in mice identified responses in lymph nodes (LNs) as essential; however, the role of LNs in human cancer patients remains unknown. We examined CD8+ T cells in human head and neck squamous cell carcinomas, regional LNs, and blood using mass cytometry, single-cell genomics, and multiplexed ion beam imaging. We identified progenitor exhausted CD8+ T cells (Tpex) that were abundant in uninvolved LN and clonally related to terminally exhausted cells in the tumor. After anti-PD-L1 immunotherapy, Tpex in uninvolved LNs reduced in frequency but localized near dendritic cells and proliferating intermediate-exhausted CD8+ T cells (Tex-int), consistent with activation and differentiation. LN responses coincided with increased circulating Tex-int. In metastatic LNs, these response hallmarks were impaired, with immunosuppressive cellular niches. Our results identify important roles for LNs in anti-tumor immune responses in humans.

Discovering dominant tumor immune archetypes in a pan-cancer census.

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

The PhyloFacts FAT-CAT web server: ortholog identification and function prediction using fast approximate tree classification.

The PhyloFacts 'Fast Approximate Tree Classification' (FAT-CAT) web server provides a novel approach to ortholog identification using subtree hidden Markov model-based placement of protein sequences to phylogenomic orthology groups in the PhyloFacts database. Results on a data set of microbial, plant and animal proteins demonstrate FAT-CAT's high precision at separating orthologs and paralogs and robustness to promiscuous domains. We also present results documenting the precision of ortholog identification based on subtree hidden Markov model scoring. The FAT-CAT phylogenetic placement is used to derive a functional annotation for the query, including confidence scores and drill-down capabilities. PhyloFacts' broad taxonomic and functional coverage, with >7.3 M proteins from across the Tree of Life, enables FAT-CAT to predict orthologs and assign function for most sequence inputs. Four pipeline parameter presets are provided to handle different sequence types, including partial sequences and proteins containing promiscuous domains; users can also modify individual parameters. PhyloFacts trees matching the query can be viewed interactively online using the PhyloScope Javascript tree viewer and are hyperlinked to various external databases. The FAT-CAT web server is available at http://phylogenomics.berkeley.edu/phylofacts/fatcat/.

Decoding functional hematopoietic progenitor cells in the adult human lung.

The bone marrow is the main site of blood cell production in adults, however, rare pools of hematopoietic stem and progenitor cells with self-renewal and differentiation potential have been found in extramedullary organs. The lung is primarily known for its role in gas exchange but has recently been described as a site of blood production in mice. Here, we show that functional hematopoietic precursors reside in the extravascular spaces of the human lung, at a frequency similar to the bone marrow, and are capable of proliferation and engraftment. The organ-specific gene signature of pulmonary and medullary CD34+ hematopoietic progenitors indicates greater baseline activation of immune, megakaryocyte/platelet and erythroid-related pathways in lung progenitors. Spatial transcriptomics mapped blood progenitors in the lung to a vascular-rich alveolar interstitium niche. These results identify the lung as a pool for uniquely programmed blood stem and progenitor cells with the potential to support hematopoiesis in humans.

Tumor-associated macrophages restrict CD8+ T cell function through collagen deposition and metabolic reprogramming of the breast cancer microenvironment.

Tumor progression is accompanied by fibrosis, a condition of excessive extracellular matrix accumulation, which is associated with diminished antitumor immune infiltration. Here we demonstrate that tumor-associated macrophages (TAMs) respond to the stiffened fibrotic tumor microenvironment (TME) by initiating a collagen biosynthesis program directed by transforming growth factor-ß. A collateral effect of this programming is an untenable metabolic milieu for productive CD8+ T cell antitumor responses, as collagen-synthesizing macrophages consume environmental arginine, synthesize proline and secrete ornithine that compromises CD8+ T cell function in female breast cancer. Thus, a stiff and fibrotic TME may impede antitumor immunity not only by direct physical exclusion of CD8+ T cells but also through secondary effects of a mechano-metabolic programming of TAMs, which creates an inhospitable metabolic milieu for CD8+ T cells to respond to anticancer immunotherapies.

Defining and using immune archetypes to classify and treat cancer.

Tumours are surrounded by a host immune system that can suppress or promote tumour growth. The tumour microenvironment (TME) has often been framed as a singular entity, suggesting a single type of immune state that is defective and in need of therapeutic intervention. By contrast, the past few years have highlighted a plurality of immune states that can surround tumours. In this Perspective, we suggest that different TMEs have 'archetypal' qualities across all cancers - characteristic and repeating collections of cells and gene-expression profiles at the level of the bulk tumour. We discuss many studies that together support a view that tumours typically draw from a finite number (around 12) of 'dominant' immune archetypes. In considering the likely evolutionary origin and roles of these archetypes, their associated TMEs can be predicted to have specific vulnerabilities that can be leveraged as targets for cancer treatment with expected and addressable adverse effects for patients.

Transcriptional space-time mapping identifies concerted immune and stromal cell patterns and gene programs in wound healing and cancer.

Tissue repair responses in metazoans are highly coordinated by different cell types over space and time. However, comprehensive single-cell-based characterization covering this coordination is lacking. Here, we captured transcriptional states of single cells over space and time during skin wound closure, revealing choreographed gene-expression profiles. We identified shared space-time patterns of cellular and gene program enrichment, which we call multicellular "movements" spanning multiple cell types. We validated some of the discovered space-time movements using large-volume imaging of cleared wounds and demonstrated the value of this analysis to predict "sender" and "receiver" gene programs in macrophages and fibroblasts. Finally, we tested the hypothesis that tumors are like "wounds that never heal" and found conserved wound healing movements in mouse melanoma and colorectal tumor models, as well as human tumor samples, revealing fundamental multicellular units of tissue biology for integrative studies.

Critical role of CD206+ macrophages in organizing anti-tumor immunity.

Tumor-associated macrophages (TAMs) are frequently and simplistically categorized as immunosuppressive, and one molecule prominently used to highlight their so-called 'M2' state is the surface protein CD206. However, direct evidence of the impact of macrophages remains impaired by the lack of sufficiently penetrant and specific tools to manipulate them in vivo. We thus made a novel conditional CD206 knock-in mouse to specifically visualize and/or deplete these TAMs. Early depletion of CD206+ macrophages and monocytes (here, 'MonoMacs') strikingly led to an indirect loss of a key anti-tumor network of NK cells, conventional type I dendritic cells (cDC1) and CD8 T cells. Among myeloid cells, we found that the CD206+ TAMs are the primary producers of CXCL9, the well-established chemoattractant for CXCR3-expressing NK and CD8 T cells. In contrast, a population of stress-responsive TAMs ("Hypoxic" or Spp1+) and immature monocytes, which remain following depletion, expressed vastly diminished levels of CXCL9. We confirmed that the missing NK and CD8 T cells are the primary producers of the cDC1-attracting chemokine Xcl1 and cDC1 growth factor Flt3l. Consistent with the loss of this critical network, CD206+ TAM depletion decreased tumor control in mice. Likewise, in humans, the CD206+ MonoMac signature correlated robustly with stimulatory cDC1 signature genes. Together, these findings negate the classification of CD206+ macrophages as immunosuppressive and instead illuminate the role of this majority of TAMs in organizing a critical tumor-reactive archetype of immunity.

Tumor cell heterogeneity drives spatial organization of the intratumoral immune response in squamous cell skin carcinoma.

Intratumoral heterogeneity (ITH)-defined as genetic and cellular diversity within a tumor-is linked to failure of immunotherapy and an inferior anti-tumor immune response. The underlying mechanism of this association is unknown. To address this question, we modeled heterogeneous tumors comprised of a pro-inflammatory ("hot") and an immunosuppressive ("cold") tumor population, labeled with YFP and RFP tags respectively to enable precise spatial tracking. The resulting mixed-population tumors exhibited distinct regions comprised of YFP+ (hot) cells, RFP+ (cold) cells, or a mixture. We found that tumor regions occupied by hot tumor cells (YFP+) harbored more total T cells and a higher frequency of Th1 cells and IFNγ+ CD8 T cells compared to regions occupied by cold tumor cells (RFP+), whereas immunosuppressive macrophages showed the opposite spatial pattern. We identified the chemokine CX3CL1, produced at higher levels by our cold tumors, as a mediator of intratumoral macrophage accumulation, particularly immunosuppressive CD206Hi macrophages. Furthermore, we examined the response of heterogeneous tumors to a therapeutic combination of PD-1 blockade and CD40 agonist on a region-by-region basis. While the combination successfully increases Th1 abundance in "cold" tumor regions, it fails to bring overall T cell activity to the same level as seen in "hot" regions. The presence of the "cold" cells thus ultimately leads to a failure of the therapy to induce tumor rejection. Collectively, our results demonstrate that the organization of heterogeneous tumor cells has a profound impact on directing the spatial organization and function of tumor-infiltrating immune cells as well as on responses to immunotherapy.

Multimodal identification of rare potent effector CD8 T cells in solid tumors.

Antitumor immunity is driven by CD8 T cells, yet we lack signatures for the exceptional effectors in tumors, amongst the vast majority of CD8 T cells undergoing exhaustion. By leveraging the measurement of a canonical T cell activation protein (CD69) together with its RNA (Cd69), we found a larger classifier for TCR stimulation-driven effector states in vitro and in vivo. This revealed exceptional 'star' effectors-highly functional cells distinguished amidst progenitor and terminally exhausted cells. Although rare in growing mouse and human tumors, they are prominent in mice during T cell-mediated tumor clearance, where they engage with tumor antigen and are superior in tumor cell killing. Employing multimodal CITE-Seq allowed de novo identification of similar rare effectors amidst T cell populations in human cancer. The identification of rare and exceptional immune states provides rational avenues for enhancement of antitumor immunity.

Cyclone: an accessible pipeline to analyze, evaluate and optimize multiparametric cytometry data.

In the past decade, high-dimensional single cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation which are computationally intense and difficult to evaluate and optimize. Here we present Cyclone, an analysis pipeline integrating dimensionality reduction, clustering, evaluation and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification, but also enables the unsupervised identification of lymphocytes and mononuclear phagocytes subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on variety of cytometry datasets which will further power immunology research and provide a scaffold for biological discovery.

Cyclone: an accessible pipeline to analyze, evaluate, and optimize multiparametric cytometry data.

In the past decade, high-dimensional single-cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation, which are computationally intense and difficult to evaluate and optimize. Here, we present Cytometry Clustering Optimization and Evaluation (Cyclone), an analysis pipeline integrating dimensionality reduction, clustering, evaluation, and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full-spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification but also enables the unsupervised identification of lymphocytes and mononuclear phagocyte subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on a variety of cytometry datasets, which will further power immunology research and provide a scaffold for biological discovery.

SATCHMO-JS: a webserver for simultaneous protein multiple sequence alignment and phylogenetic tree construction.

We present the jump-start simultaneous alignment and tree construction using hidden Markov models (SATCHMO-JS) web server for simultaneous estimation of protein multiple sequence alignments (MSAs) and phylogenetic trees. The server takes as input a set of sequences in FASTA format, and outputs a phylogenetic tree and MSA; these can be viewed online or downloaded from the website. SATCHMO-JS is an extension of the SATCHMO algorithm, and employs a divide-and-conquer strategy to jump-start SATCHMO at a higher point in the phylogenetic tree, reducing the computational complexity of the progressive all-versus-all HMM-HMM scoring and alignment. Results on a benchmark dataset of 983 structurally aligned pairs from the PREFAB benchmark dataset show that SATCHMO-JS provides a statistically significant improvement in alignment accuracy over MUSCLE, Multiple Alignment using Fast Fourier Transform (MAFFT), ClustalW and the original SATCHMO algorithm. The SATCHMO-JS webserver is available at http://phylogenomics.berkeley.edu/satchmo-js. The datasets used in these experiments are available for download at http://phylogenomics.berkeley.edu/satchmo-js/supplementary/.

Spatiotemporal co-dependency between macrophages and exhausted CD8+ T cells in cancer.

T cell exhaustion is a major impediment to antitumor immunity. However, it remains elusive how other immune cells in the tumor microenvironment (TME) contribute to this dysfunctional state. Here, we show that the biology of tumor-associated macrophages (TAMs) and exhausted T cells (Tex) in the TME is extensively linked. We demonstrate that in vivo depletion of TAMs reduces exhaustion programs in tumor-infiltrating CD8+ T cells and reinvigorates their effector potential. Reciprocally, transcriptional and epigenetic profiling reveals that Tex express factors that actively recruit monocytes to the TME and shape their differentiation. Using lattice light sheet microscopy, we show that TAM and CD8+ T cells engage in unique, long-lasting, antigen-specific synaptic interactions that fail to activate T cells but prime them for exhaustion, which is then accelerated in hypoxic conditions. Spatially resolved sequencing supports a spatiotemporal self-enforcing positive feedback circuit that is aligned to protect rather than destroy a tumor.

Holistic Characterization of Tumor Monocyte-to-Macrophage Differentiation Integrates Distinct Immune Phenotypes in Kidney Cancer.

The tumor immune microenvironment (TIME) is commonly infiltrated by diverse collections of myeloid cells. Yet, the complexity of myeloid-cell identity and plasticity has challenged efforts to define bona fide populations and determine their connections to T-cell function and their relationship to patient outcome. Here, we have leveraged single-cell RNA-sequencing analysis of several mouse and human tumors and found that monocyte-macrophage diversity is characterized by a combination of conserved lineage states as well as transcriptional programs accessed along the differentiation trajectory. We also found in mouse models that tumor monocyte-to-macrophage progression was profoundly tied to regulatory T cell (Treg) abundance. In human kidney cancer, heterogeneity in macrophage accumulation and myeloid composition corresponded to variance in, not only Treg density, but also the quality of infiltrating CD8+ T cells. In this way, holistic analysis of monocyte-to-macrophage differentiation creates a framework for critically different immune states.

Berkeley PHOG: PhyloFacts orthology group prediction web server.

Ortholog detection is essential in functional annotation of genomes, with applications to phylogenetic tree construction, prediction of protein-protein interaction and other bioinformatics tasks. We present here the PHOG web server employing a novel algorithm to identify orthologs based on phylogenetic analysis. Results on a benchmark dataset from the TreeFam-A manually curated orthology database show that PHOG provides a combination of high recall and precision competitive with both InParanoid and OrthoMCL, and allows users to target different taxonomic distances and precision levels through the use of tree-distance thresholds. For instance, OrthoMCL-DB achieved 76% recall and 66% precision on this dataset; at a slightly higher precision (68%) PHOG achieves 10% higher recall (86%). InParanoid achieved 87% recall at 24% precision on this dataset, while a PHOG variant designed for high recall achieves 88% recall at 61% precision, increasing precision by 37% over InParanoid. PHOG is based on pre-computed trees in the PhyloFacts resource, and contains over 366 K orthology groups with a minimum of three species. Predicted orthologs are linked to GO annotations, pathway information and biological literature. The PHOG web server is available at http://phylofacts.berkeley.edu/orthologs/.

Active surveillance characterizes human intratumoral T cell exhaustion.

Intratumoral T cells that might otherwise control tumors are often identified in an "exhausted" state, defined by specific epigenetic modifications and upregulation of genes such as CD38, cytotoxic T-lymphocyte-associated protein 4 (CTLA4), and programmed cell death 1 (PD1). Although the term might imply inactivity, there has been little study of this state at the phenotypic level in tumors to understand the extent of their incapacitation. Starting with the observation that T cells move more quickly through mouse tumors the longer they reside there and progress toward exhaustion, we developed a nonstimulatory, live-biopsy method for the real-time study of T cell behavior within individual patient tumors. Using 2-photon microscopy, we studied native CD8+ T cell interaction with antigen-presenting cells (APCs) and cancer cells in different microniches of human tumors and found that T cell speed was variable by region and by patient and was inversely correlated with local tumor density. Across a range of tumor types, we found a strong relationship between CD8+ T cell motility and the exhausted T cell state that corresponded with our observations made in mouse models in which exhausted T cells moved faster. Our study demonstrates T cell dynamic states in individual human tumors and supports the existence of an active program in "exhausted" T cells that extends beyond incapacitating them.

SCENITH: A Flow Cytometry-Based Method to Functionally Profile Energy Metabolism with Single-Cell Resolution.

Energetic metabolism reprogramming is critical for cancer and immune responses. Current methods to functionally profile the global metabolic capacities and dependencies of cells are performed in bulk. We designed a simple method for complex metabolic profiling called SCENITH, for single-cell energetic metabolism by profiling translation inhibition. SCENITH allows for the study of metabolic responses in multiple cell types in parallel by flow cytometry. SCENITH is designed to perform metabolic studies ex vivo, particularly for rare cells in whole blood samples, avoiding metabolic biases introduced by culture media. We analyzed myeloid cells in solid tumors from patients and identified variable metabolic profiles, in ways that are not linked to their lineage or their activation phenotype. SCENITH's ability to reveal global metabolic functions and determine complex and linked immune-phenotypes in rare cell subpopulations will contribute to the information needed for evaluating therapeutic responses or patient stratification.

A natural killer-dendritic cell axis defines checkpoint therapy-responsive tumor microenvironments.

Intratumoral stimulatory dendritic cells (SDCs) play an important role in stimulating cytotoxic T cells and driving immune responses against cancer. Understanding the mechanisms that regulate their abundance in the tumor microenvironment (TME) could unveil new therapeutic opportunities. We find that in human melanoma, SDC abundance is associated with intratumoral expression of the gene encoding the cytokine FLT3LG. FLT3LG is predominantly produced by lymphocytes, notably natural killer (NK) cells in mouse and human tumors. NK cells stably form conjugates with SDCs in the mouse TME, and genetic and cellular ablation of NK cells in mice demonstrates their importance in positively regulating SDC abundance in tumor through production of FLT3L. Although anti-PD-1 'checkpoint' immunotherapy for cancer largely targets T cells, we find that NK cell frequency correlates with protective SDCs in human cancers, with patient responsiveness to anti-PD-1 immunotherapy, and with increased overall survival. Our studies reveal that innate immune SDCs and NK cells cluster together as an excellent prognostic tool for T cell-directed immunotherapy and that these innate cells are necessary for enhanced T cell tumor responses, suggesting this axis as a target for new therapies.