This is a PDF file of a peer-reviewed paper that has been accepted for publication. Although unedited, the content has been subjected to preliminary formatting. Nature is providing this early version of the typeset paper as a service to our authors and readers. The text and figures will undergo copyediting and a proof review before the paper is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers apply.

Mass cytometry allows high-resolution dissection of the cellular composition of the immune system. However, the high-dimensionality, large size, and non-linear structure of the data poses considerable challenges for the data analysis. In particular, dimensionality reduction-based techniques like t-SNE offer single-cell resolution but are limited in the number of cells that can be analyzed. Here we introduce Hierarchical Stochastic Neighbor Embedding (HSNE) for the analysis of mass cytometry data sets. HSNE constructs a hierarchy of non-linear similarities that can be interactively explored with a stepwise increase in detail up to the single-cell level. We apply HSNE to a study on gastrointestinal disorders and three other available mass cytometry data sets. We find that HSNE efficiently replicates previous observations and identifies rare cell populations that were previously missed due to downsampling. Thus, HSNE removes the scalability limit of conventional t-SNE analysis, a feature that makes it highly suitable for the analysis of massive high-dimensional data sets.

Purpose: Human papillomavirus (HPV)-associated oropharyngeal squamous cell cancer (OPSCC) has a much better prognosis than HPV-negative OPSCC, and this is linked to dense tumor immune infiltration. As the viral antigens may trigger potent immunity, we studied the relationship between the presence of intratumoral HPV-specific T-cell responses, the immune contexture in the tumor microenvironment, and clinical outcome.Experimental Design: To this purpose, an in-depth analysis of tumor-infiltrating immune cells in a prospective cohort of 97 patients with HPV16-positive and HPV16-negative OPSCC was performed using functional T-cell assays, mass cytometry (CyTOF), flow cytometry, and fluorescent immunostaining of tumor tissues. Key findings were validated in a cohort of 75 patients with HPV16-positive OPSCC present in the publicly available The Cancer Genome Atlas database.Results: In 64% of the HPV16-positive tumors, type I HPV16-specific T cells were present. Their presence was not only strongly related to a better overall survival, a smaller tumor size, and less lymph node metastases but also to a type I-oriented tumor microenvironment, including high numbers of activated CD161 þ T cells, CD103 þ tissue-resident T cells, dendritic cells (DC), and DC-like macrophages. Conclusions:The viral antigens trigger a tumor-specific T-cell response that shapes a favorable immune contexture for the response to standard therapy. Hence, reinforcement of HPV16-specific T-cell reactivity is expected to boost this process. Clin Cancer Res; 24(3); 634-47. Ó2017 AACR.

SUMMARY Tissue-resident lymphocytes play a key role in immune surveillance, but it remains unclear how these inherently stable cell populations respond to chronic inflammation. In the setting of celiac disease (CeD), where exposure to dietary antigen can be controlled, gluten-induced inflammation triggered a profound depletion of naturally occurring Vγ4+/Vδ1+ intraepithelial lymphocytes (IELs) with innate cytolytic properties and specificity for the butyrophilin-like (BTNL) molecules BTNL3/BTNL8. Creation of a new niche with reduced expression of BTNL8 and loss of Vγ4+/Vδ1+ IELs was accompanied by the expansion of gluten-sensitive, interferon-γ-producing Vδ1+ IELs bearing T cell receptors (TCRs) with a shared non-germ-line-encoded motif that failed to recognize BTNL3/BTNL8. Exclusion of dietary gluten restored BTNL8 expression but was insufficient to reconstitute the physiological Vγ4+/Vδ1+ subset among TCRγδ+ IELs. Collectively, these data show that chronic inflammation permanently reconfigures the tissue-resident TCRγδ+ IEL compartment in CeD.

The fetus is thought to be protected from exposure to foreign antigens, yet CD45RO+ T cells reside in the fetal intestine. Here we combined functional assays with mass cytometry, single-cell RNA-sequencing and high-throughput T cell antigen receptor (TCR) sequencing to characterize the CD4+ T cell compartment in the human fetal intestine. We identified 22 CD4+ T cell clusters, including naive-like, regulatory-like and memory-like subpopulations, which were confirmed and further characterized at the transcriptional level. Memory-like CD4+ T cells had high expression of Ki-67, indicative of cell division, and CD5, a surrogate marker of TCR avidity, and produced the cytokines IFN-γ and IL-2. Pathway analysis revealed a differentiation trajectory associated with cellular activation and proinflammatory effector functions, and TCR repertoire analysis indicated clonal expansions, distinct repertoire characteristics and interconnections between subpopulations of memory-like CD4+ T cells. Imaging-mass cytometry indicated that memory-like CD4+ T cells colocalized with antigen-presenting cells. Collectively, these results provide evidence for the generation of memory-like CD4+ T cells in the human fetal intestine that is consistent with exposure to foreign antigens.

Inflammatory intestinal diseases are characterized by abnormal immune responses and affect distinct locations of the gastrointestinal tract. Although the role of several immune subsets in driving intestinal pathology has been studied, a system-wide approach that simultaneously interrogates all major lineages on a single-cell basis is lacking. We used high-dimensional mass cytometry to generate a system-wide view of the human mucosal immune system in health and disease. We distinguished 142 immune subsets and through computational applications found distinct immune subsets in peripheral blood mononuclear cells and intestinal biopsies that distinguished patients from controls. In addition, mucosal lymphoid malignancies were readily detected as well as precursors from which these likely derived. These findings indicate that an integrated high-dimensional analysis of the entire immune system can identify immune subsets associated with the pathogenesis of complex intestinal disorders. This might have implications for diagnostic procedures, immune-monitoring, and treatment of intestinal diseases and mucosal malignancies.

Figure 1: Cytosplore. Screenshot of our system with four widgets (adaptive settings, overview (a), embedding (b) and heatmap (c)), representing the workflow. Views can be rearranged or additional views of these types added. AbstractTo understand how the immune system works, one needs to have a clear picture of its cellular compositon and the cells' corresponding properties and functionality. Mass cytometry is a novel technique to determine the properties of single-cells with unprecedented detail. This amount of detail allows for much finer differentiation but also comes at the cost of more complex analysis. In this work, we present Cytosplore, implementing an interactive workflow to analyze mass cytometry data in an integrated system, providing multiple linked views, showing different levels of detail and enabling the rapid definition of known and unknown cell types. Cytosplore handles millions of cells, each represented as a high-dimensional data point, facilitates hypothesis generation and confirmation, and provides a significant speed up of the current workflow. We show the effectiveness of Cytosplore in a case study evaluation.

ObjectiveA comprehensive understanding of anticancer immune responses is paramount for the optimal application and development of cancer immunotherapies. We unravelled local and systemic immune profiles in patients with colorectal cancer (CRC) by high-dimensional analysis to provide an unbiased characterisation of the immune contexture of CRC.DesignThirty-six immune cell markers were simultaneously assessed at the single-cell level by mass cytometry in 35 CRC tissues, 26 tumour-associated lymph nodes, 17 colorectal healthy mucosa and 19 peripheral blood samples from 31 patients with CRC. Additionally, functional, transcriptional and spatial analyses of tumour-infiltrating lymphocytes were performed by flow cytometry, single-cell RNA-sequencing and multispectral immunofluorescence.ResultsWe discovered that a previously unappreciated innate lymphocyte population (Lin–CD7+CD127–CD56+CD45RO+) was enriched in CRC tissues and displayed cytotoxic activity. This subset demonstrated a tissue-resident (CD103+CD69+) phenotype and was most abundant in immunogenic mismatch repair (MMR)-deficient CRCs. Their presence in tumours was correlated with the infiltration of tumour-resident cytotoxic, helper and γδ T cells with highly similar activated (HLA-DR+CD38+PD-1+) phenotypes. Remarkably, activated γδ T cells were almost exclusively found in MMR-deficient cancers. Non-activated counterparts of tumour-resident cytotoxic and γδ T cells were present in CRC and healthy mucosa tissues, but not in lymph nodes, with the exception of tumour-positive lymph nodes.ConclusionThis work provides a blueprint for the understanding of the heterogeneous and intricate immune landscape of CRC, including the identification of previously unappreciated immune cell subsets. The concomitant presence of tumour-resident innate and adaptive immune cell populations suggests a multitargeted exploitation of their antitumour properties in a therapeutic setting.