IMMUNOLOGY. An interactive reference framework for modeling a dynamic immune system

SUMOylation inhibitors activate anti-tumor immunity by reshaping the immune microenvironment in a preclinical model of hepatocellular carcinoma

PURPOSE

High levels of heterogeneity and immunosuppression characterize the HCC immune microenvironment (TME). Unfortunately, the majority of hepatocellular carcinoma (HCC) patients do not benefit from immune checkpoint inhibitors (ICIs) therapy. New small molecule therapies for the treatment of HCC are the goal of our research.

METHODS

SUMOylation inhibitors (TAK-981 and ML-792) were evaluated for the treatment of preclinical mouse HCC models (including subcutaneous and orthotopic HCC models). We profile immune cell subsets from tumor samples after SUMOylation inhibitors treatment using single-cell RNA sequencing (scRNA-seq), mass cytometry (CyTOF), flow cytometry, and multiple immunofluorescences (mIF).

RESULTS

We discover that SUMOylation is higher in HCC patient samples compared to normal liver tissue. TAK-981 and ML-792 decrease SUMOylation at nanomolar levels in HCC cells and also successfully reduced the tumor burden. Analysis combining scRNA-seq and CyTOF demonstrate that treatment with SUMOylation inhibitors reduces the exhausted CD8+T (Tex) cells while enhancing the cytotoxic NK cells, M1 macrophages and cytotoxic T lymphocytes (CTL) in preclinical mouse HCC model. Furthermore, SUMOylation inhibitors have the potential to activate innate immune signals from CD8+T, NK and macrophages while promoting TNFα and IL-17 secretion. Most notably, SUMOylation inhibitors can directly alter the TME by adjusting the abundance of intestinal microbiota, thereby restoring anti-tumor immunity in HCC models.

CONCLUSIONS

This preclinical study suggests that SUMO signaling inhibitors may be beneficial for the treatment of HCC.

Cell-type specific molecular expression levels by restricted-dimensional cytometry

BACKGROUND

Cytometric analysis has been commonly used to delineate distinct cell subpopulations among peripheral blood mononuclear cells by the differential expression of surface receptors. This capability has reached its apogee with high-dimensional approaches such as mass cytometry and spectral cytometry that include simultaneous assessment of 20-50 analytes. Unfortunately, this approach also engenders significant complexity with analytical and interpretational pitfalls.

METHODS

Here, we demonstrate a complementary approach with restricted-dimensionality to assess cell-type specific intracellular molecular expression levels at exceptional levels of precision. The expression of five analytes was individually assessed in four mononuclear cell-types from peripheral blood.

RESULTS

Distinctions in expression levels were seen between cell-types and between samples from different donor groups. Mononuclear cell-type specific molecular expression levels distinguished pregnant from nonpregnant women and G-CSF-treated from untreated persons. Additionally, the precision of our analysis was sufficient to quantify a novel relationship between two molecules-Rel A and translocator protein-by correlational analysis.

CONCLUSIONS

Restricted-dimensional cytometry can provide a complementary approach to define characteristics of cell-type specific intracellular protein and phosphoantigen expression in mononuclear cells.

The tumor microenvironment of benign and malignant salivary gland tumors

BACKGROUND

Treatment of salivary gland tumors (SGTs) remains challenging. Little is known about the immune landscape of SGTs. We aimed to characterize the tumor microenvironment in benign and malignant SGTs.

METHODS

Eleven benign and nine malignant tumors were collected from patients undergoing curative intent surgery. Specimens were analyzed using mass cytometry by time-of-flight. Immune cell populations were manually gated, and T cells were clustered using the FlowSOM algorithm. Population frequencies were compared between high-grade and low-grade malignancies, corrected for multiple hypothesis testing.

RESULTS

There were trends towards increased CD4+ and CD8+ T cells among malignant tumors. High-grade malignancies exhibited trends towards higher frequencies of CD8+ PD-1+ CD39+ CD103+ exhausted T cells, CD4+ FoxP3+ TCF-1+ CD127- Tregs, and CD69+ CD25- CD4+ T cells compared to low-grade malignancies.

CONCLUSION

SGTs exhibit significant immunologic diversity. High-grade malignancies tended to have greater infiltration of exhausted CD8+ T cells and Tregs, which may guide future studies for immunotherapy strategies.

Neoadjuvant CD40 Agonism Remodels the Tumor Immune Microenvironment in Locally Advanced Esophageal/Gastroesophageal Junction Cancer

Sotigalimab is an agonistic anti-CD40 mAb that can modulate antitumor immune responses. In a phase II clinical trial of sotigalimab combined with neoadjuvant chemoradiation (CRT) in locally advanced esophageal/gastroesophageal junction (E/GEJ) cancer with the primary outcome of efficacy as measured by pathologic complete response (pCR) rate, the combination induced pCR in 38% of treated patients. We investigated the mechanism of action of sotigalimab in samples obtained from this clinical trial. Tumor biopsies and peripheral blood samples were collected at baseline, following an initial dose of sotigalimab, and at the time of surgery after CRT completion from six patients. High dimensional single-cell techniques were used, including combined single-cell RNA-sequencing and proteomics (CITEseq) and multiplexed ion beam imaging, to analyze immune responses. Sotigalimab dramatically remodeled the immune compartment in the periphery and within the tumor microenvironment (TME), increasing expression of molecules related to antigen processing and presentation and altering metabolic pathways in myeloid cells. Concomitant with these changes in myeloid cells, sotigalimab treatment primed new T cell clonotypes and increased the density and activation of T cells with enhanced cytotoxic function. Sotigalimab treatment also induced a decrease in the frequency of Tregs in the TME. These findings indicate that a single dose of sotigalimab leads to enhanced antigen presentation that can activate T cells and induce new T cell clones. This restructuring of the TME provides elements which are critical to the development of effective antitumor immune responses and improved clinical outcomes.

A novel mass cytometry protocol optimized for immunophenotyping of low-frequency antigen-specific T cells

Understanding antigen-specific T-cell responses, for example, following virus infections or allergen exposure, is of high relevance for the development of vaccines and therapeutics. We aimed on optimizing immunophenotyping of T cells after antigen stimulation by improving staining procedures for flow and mass cytometry. Our method can be used for primary cells of both mouse and human origin for the detection of low-frequency T-cell response using a dual-barcoding system for individual samples and conditions. First, live-cell barcoding was performed using anti-CD45 antibodies prior to an in vitro T-cell stimulation assay. Second, to discriminate between stimulation conditions and prevent cell loss, sample barcoding was combined with a commercial barcoding solution. This dual-barcoding approach is cell sparing and, therefore, particularly relevant for samples with low cell numbers. To further reduce cell loss and to increase debarcoding efficiency of multiplexed samples, we combined our dual-barcoding approach with a new centrifugation-free washing system by laminar flow (Curiox™). Finally, to demonstrate the benefits of our established protocol, we assayed virus-specific T-cell response in SARS-CoV-2-vaccinated and SARS-CoV-2-infected patients and compared with healthy non-exposed individuals by a high-parameter CyTOF analysis. We could reveal a heterogeneity of phenotypes among responding CD4, CD8, and gd-T cells following antigen-specific stimulations. Our protocol allows to assay antigen-specific responses of minute populations of T cells to virus-derived peptides, allergens, or other antigens from the same donor sample, in order to investigate qualitative and quantitative differences.

Epigenetic reprogramming shapes the cellular landscape of schwannoma

Mechanisms specifying cancer cell states and response to therapy are incompletely understood. Here we show epigenetic reprogramming shapes the cellular landscape of schwannomas, the most common tumors of the peripheral nervous system. We find schwannomas are comprised of 2 molecular groups that are distinguished by activation of neural crest or nerve injury pathways that specify tumor cell states and the architecture of the tumor immune microenvironment. Moreover, we find radiotherapy is sufficient for interconversion of neural crest schwannomas to immune-enriched schwannomas through epigenetic and metabolic reprogramming. To define mechanisms underlying schwannoma groups, we develop a technique for simultaneous interrogation of chromatin accessibility and gene expression coupled with genetic and therapeutic perturbations in single-nuclei. Our results elucidate a framework for understanding epigenetic drivers of tumor evolution and establish a paradigm of epigenetic and metabolic reprograming of cancer cells that shapes the immune microenvironment in response to radiotherapy.

A mass cytometry approach to track the evolution of T cell responses during infection and immunotherapy by paired T cell receptor repertoire and T cell differentiation state analysis

T cell receptor (TCR) recognition followed by clonal expansion is a fundamental feature of adaptive immune responses. Here, we developed a mass cytometric (CyTOF) approach combining antibodies specific for different TCR Vα- and Vβ-chains with antibodies against T cell activation and differentiation proteins to identify antigen-specific expansions of T cell subsets and assess aspects of cellular function. This strategy allowed for the identification of expansions of specific Vβ and Vα chain expressing CD8+ and CD4+ T cells with varying differentiation states in response to Listeria monocytogenes, tumors, and respiratory influenza infection. Expanded Vβ chain expressing T cells could be directly linked to the recognition of specific antigens from Listeria, tumor cells, or influenza. In the setting of influenza infection, we showed that the common therapeutic approaches of intramuscular vaccination or convalescent serum transfer altered the clonal diversity and differentiation state of responding T cells. Thus, we present a new method to monitor broad changes in TCR specificity paired with T cell differentiation during adaptive immune responses.

Single cell proteomics characterization of bone marrow hematopoiesis with distinct Ras pathway lesions

Normal hematopoiesis requires constant prolific production of different blood cell lineages by multipotent hematopoietic stem cells (HSC). Stem- and progenitor- cells need to balance dormancy with proliferation. How genetic alterations impact frequency, lineage potential, and metabolism of HSC is largely unknown. Here, we compared induced expression of KRAS G12D or RasGRP1 to normal hematopoiesis. At low-resolution, both Ras pathway lesions result in skewing towards myeloid lineages. Single-cell resolution CyTOF proteomics unmasked an expansion of HSC- and progenitor- compartments for RasGRP1, contrasted by a depletion for KRAS G12D . SCENITH™ quantitates protein synthesis with single-cell precision and corroborated that immature cells display low metabolic SCENITH™ rates. Both RasGRP1 and KRAS G12D elevated mean SCENITH™ signals in immature cells. However, RasGRP1-overexpressing stem cells retain a metabolically quiescent cell-fraction, whereas this fraction diminishes for KRAS G12D . Our temporal single cell proteomics and metabolomics datasets provide a resource of mechanistic insights into altered hematopoiesis at single cell resolution.

Trellis tree-based analysis reveals stromal regulation of patient-derived organoid drug responses

Patient-derived organoids (PDOs) can model personalized therapy responses; however, current screening technologies cannot reveal drug response mechanisms or how tumor microenvironment cells alter therapeutic performance. To address this, we developed a highly multiplexed mass cytometry platform to measure post-translational modification (PTM) signaling, DNA damage, cell-cycle activity, and apoptosis in >2,500 colorectal cancer (CRC) PDOs and cancer-associated fibroblasts (CAFs) in response to clinical therapies at single-cell resolution. To compare patient- and microenvironment-specific drug responses in thousands of single-cell datasets, we developed "Trellis"-a highly scalable, tree-based treatment effect analysis method. Trellis single-cell screening revealed that on-target cell-cycle blockage and DNA-damage drug effects are common, even in chemorefractory PDOs. However, drug-induced apoptosis is rarer, patient-specific, and aligns with cancer cell PTM signaling. We find that CAFs can regulate PDO plasticity-shifting proliferative colonic stem cells (proCSCs) to slow-cycling revival colonic stem cells (revCSCs) to protect cancer cells from chemotherapy.

An oncogenic phenoscape of colonic stem cell polarization

Cancer cells are regulated by oncogenic mutations and microenvironmental signals, yet these processes are often studied separately. To functionally map how cell-intrinsic and cell-extrinsic cues co-regulate cell fate, we performed a systematic single-cell analysis of 1,107 colonic organoid cultures regulated by (1) colorectal cancer (CRC) oncogenic mutations, (2) microenvironmental fibroblasts and macrophages, (3) stromal ligands, and (4) signaling inhibitors. Multiplexed single-cell analysis revealed a stepwise epithelial differentiation phenoscape dictated by combinations of oncogenes and stromal ligands, spanning from fibroblast-induced Clusterin (CLU)+ revival colonic stem cells (revCSCs) to oncogene-driven LRIG1+ hyper-proliferative CSCs (proCSCs). The transition from revCSCs to proCSCs is regulated by decreasing WNT3A and TGF-β-driven YAP signaling and increasing KRASG12D or stromal EGF/Epiregulin-activated MAPK/PI3K flux. We find that APC loss and KRASG12D collaboratively limit access to revCSCs and disrupt stromal-epithelial communication-trapping epithelia in the proCSC fate. These results reveal that oncogenic mutations dominate homeostatic differentiation by obstructing cell-extrinsic regulation of cell-fate plasticity.

Simultaneous selection of nanobodies for accessible epitopes on immune cells in the tumor microenvironment

In the rapidly advancing field of synthetic biology, there exists a critical need for technology to discover targeting moieties for therapeutic biologics. Here we present INSPIRE-seq, an approach that utilizes a nanobody library and next-generation sequencing to identify nanobodies selected for complex environments. INSPIRE-seq enables the parallel enrichment of immune cell-binding nanobodies that penetrate the tumor microenvironment. Clone enrichment and specificity vary across immune cell subtypes in the tumor, lymph node, and spleen. INSPIRE-seq identifies a dendritic cell binding clone that binds PHB2. Single-cell RNA sequencing reveals a connection with cDC1s, and immunofluorescence confirms nanobody-PHB2 colocalization along cell membranes. Structural modeling and docking studies assist binding predictions and will guide nanobody selection. In this work, we demonstrate that INSPIRE-seq offers an unbiased approach to examine complex microenvironments and assist in the development of nanobodies, which could serve as active drugs, modified to become drugs, or used as targeting moieties.

Mass cytometry as a tool in target validation and drug discovery

Mass cytometry provides highly multiparametric data at a single cell level, coupling the specificity and sensitivity of time-of-flight mass spectrometry with the single-cell throughput of flow cytometry. It offers great value in interrogating the potentially heterogenous impact that a drug may have on a biological system, allowing an investigator to capture not just changes in cell behavior, but how these changes may differ between cell subtypes. In this chapter, we review the technical details of the platform as well as its limitations, before describing our approach to planning and running a mass cytometry experiment. A series of method modules, spanning the staining process through to data cleaning, are described that are then combined to create three separate experiments. The first experiment illustrates a core process in mass cytometry: the validation and titration of a metal-conjugated antibody reporter. The second experiment explores the impact of a kinase inhibitor on cell cycle and apoptosis pathways of a human myeloma cell line. And the third experiment exploits the multiparametric capability of mass cytometry, by exploring the differential expression changes in a transcription factor upon drug treatment across the cellular compartments of a peripheral blood mononuclear cell sample.

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.

Evolving Approach to Clinical Cytometry for Immunodeficiencies and Other Immune Disorders

Primary immunodeficiencies were initially identified on the basis of recurrent, severe or unusual infections. Subsequently, it was noted that these diseases can also manifest with autoimmunity, autoinflammation, allergy, lymphoproliferation and malignancy, hence a conceptual change and their renaming as inborn errors of immunity. Ongoing advances in flow cytometry provide the opportunity to expand or modify the utility and scope of existing laboratory tests in this field to mirror this conceptual change. Here we have used the B cell subset, variably known as CD21low B cells, age-associated B cells and T-bet+ B cells, as an example to demonstrate this possibility.

Computational Methods for Single-Cell Proteomics

Advances in single-cell proteomics technologies have resulted in high-dimensional datasets comprising millions of cells that are capable of answering key questions about biology and disease. The advent of these technologies has prompted the development of computational tools to process and visualize the complex data. In this review, we outline the steps of single-cell and spatial proteomics analysis pipelines. In addition to describing available methods, we highlight benchmarking studies that have identified advantages and pitfalls of the currently available computational toolkits. As these technologies continue to advance, robust analysis tools should be developed in tandem to take full advantage of the potential biological insights provided by these data.

Transitional dendritic cells are distinct from conventional DC2 precursors and mediate proinflammatory antiviral responses

High-dimensional approaches have revealed heterogeneity amongst dendritic cells (DCs), including a population of transitional DCs (tDCs) in mice and humans. However, the origin and relationship of tDCs to other DC subsets has been unclear. Here we show that tDCs are distinct from other well-characterized DCs and conventional DC precursors (pre-cDCs). We demonstrate that tDCs originate from bone marrow progenitors shared with plasmacytoid DCs (pDCs). In the periphery, tDCs contribute to the pool of ESAM+ type 2 DCs (DC2s), and these DC2s have pDC-related developmental features. Different from pre-cDCs, tDCs have less turnover, capture antigen, respond to stimuli and activate antigen-specific naïve T cells, all characteristics of differentiated DCs. Different from pDCs, viral sensing by tDCs results in IL-1β secretion and fatal immune pathology in a murine coronavirus model. Our findings suggest that tDCs are a distinct pDC-related subset with a DC2 differentiation potential and unique proinflammatory function during viral infections.

Single-Cell Proteomics with Spatial Attributes: Tools and Techniques

Now-a-days, the single-cell proteomics (SCP) concept is attracting interest, especially in clinical research, because it can identify the proteomic signature specific to diseased cells. This information is very essential when dealing with the progression of certain diseases, such as cancer, diabetes, Alzheimer's, etc. One of the major drawbacks of conventional destructive proteomics is that it gives an average idea about the protein expression profile in the disease condition. During the extraction of the protein from a biopsy or blood sample, proteins may come from both diseased cells and adjacent normal cells or any other cells from the disease environment. Again, SCP along with spatial attributes is utilized to learn about the heterogeneous function of a single protein. Before performing SCP, it is necessary to isolate single cells. This can be done by various techniques, including fluorescence-activated cell sorting (FACS), magnetic-activated cell sorting (MACS), laser capture microdissection (LCM), microfluidics, manual cell picking/micromanipulation, etc. Among the different approaches for proteomics, mass spectrometry-based proteomics tools are widely used for their high resolution as well as sensitivity. This Review mainly focuses on the mass spectrometry-based approaches for the study of single-cell proteomics.

ANPELA: Significantly Enhanced Quantification Tool for Cytometry-Based Single-Cell Proteomics

ANPELA is widely used for quantifying traditional bulk proteomic data. Recently, there is a clear shift from bulk proteomics to the single-cell ones (SCP), for which powerful cytometry techniques demonstrate the fantastic capacity of capturing cellular heterogeneity that is completely overlooked by traditional bulk profiling. However, the in-depth and high-quality quantification of SCP data is still challenging and severely affected by the large numbers of quantification workflows and extreme performance dependence on the studied datasets. In other words, the proper selection of well-performing workflow(s) for any studied dataset is elusory, and it is urgently needed to have a significantly enhanced and accelerated tool to address this issue. However, no such tool is developed yet. Herein, ANPELA is therefore updated to its 2.0 version (https://idrblab.org/anpela/), which is unique in providing the most comprehensive set of quantification alternatives (>1000 workflows) among all existing tools, enabling systematic performance evaluation from multiple perspectives based on machine learning, and identifying the optimal workflow(s) using overall performance ranking together with the parallel computation. Extensive validation on different benchmark datasets and representative application scenarios suggest the great application potential of ANPELA in current SCP research for gaining more accurate and reliable biological insights.

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.

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.

INFLECT: an R-package for cytometry cluster evaluation using marker modality

Background

Current methods of high-dimensional unsupervised clustering of mass cytometry data lack means to monitor and evaluate clustering results. Whether unsupervised clustering is correct is typically evaluated by agreement with dimensionality reduction techniques or based on benchmarking with manually classified cells. The ambiguity and lack of reproducibility of sequential gating has been replaced with ambiguity in interpretation of clustering results. On the other hand, spurious overclustering of data leads to loss of statistical power. We have developed INFLECT, an R-package designed to give insight in clustering results and provide an optimal number of clusters. In our approach, a mass cytometry dataset is overclustered intentionally to ensure the smallest phenotypically different subsets are captured using FlowSOM. A range of metacluster number endpoints are generated and evaluated using marker interquartile range and distribution unimodality checks. The fraction of marker distributions that pass these checks is taken as a measure of clustering success. The fraction of unimodal distributions within metaclusters is plotted against the number of generated metaclusters and reaches a plateau of diminishing returns. The inflection point at which this occurs gives an optimal point of capturing cellular heterogeneity versus statistical power.

Results

We applied INFLECT to four publically available mass cytometry datasets of different size and number of markers. The unimodality score consistently reached a plateau, with an inflection point dependent on dataset size and number of dimensions. We tested both ConsenusClusterPlus metaclustering and hierarchical clustering. While hierarchical clustering is less computationally expensive and thus faster, it achieved similar results to ConsensusClusterPlus. The four datasets consisted of labeled data and we compared INFLECT metaclustering to published results. INFLECT identified a higher optimal number of metaclusters for all datasets. We illustrated the underlying heterogeneity within labels, showing that these labels encompass distinct types of cells.

Conclusion

INFLECT addresses a knowledge gap in high-dimensional cytometry analysis, namely assessing clustering results. This is done through monitoring marker distributions for interquartile range and unimodality across a range of metacluster numbers. The inflection point is the optimal trade-off between cellular heterogeneity and statistical power, applied in this work for FlowSOM clustering on mass cytometry datasets.

Single-cell profiling reveals a memory B cell-like subtype of follicular lymphoma with increased transformation risk

Follicular lymphoma (FL) is an indolent cancer of mature B-cells but with ongoing risk of transformation to more aggressive histology over time. Recurrent mutations associated with transformation have been identified; however, prognostic features that can be discerned at diagnosis could be clinically useful. We present here comprehensive profiling of both tumor and immune compartments in 155 diagnostic FL biopsies at single-cell resolution by mass cytometry. This revealed a diversity of phenotypes but included two recurrent patterns, one which closely resembles germinal center B-cells (GCB) and another which appears more related to memory B-cells (MB). GCB-type tumors are enriched for EZH2, TNFRSF14, and MEF2B mutations, while MB-type tumors contain increased follicular helper T-cells. MB-type and intratumoral phenotypic diversity are independently associated with increased risk of transformation, supporting biological relevance of these features. Notably, a reduced 26-marker panel retains sufficient information to allow phenotypic profiling of future cohorts by conventional flow cytometry.

Innate type 2 immunity controls hair follicle commensalism by Demodex mites

Demodex mites are commensal parasites of hair follicles (HFs). Normally asymptomatic, inflammatory outgrowth of mites can accompany malnutrition, immune dysfunction, and aging, but mechanisms restricting Demodex outgrowth are not defined. Here, we show that control of mite HF colonization in mice required group 2 innate lymphoid cells (ILC2s), interleukin-13 (IL-13), and its receptor, IL-4Ra-IL-13Ra1. HF-associated ILC2s elaborated IL-13 that attenuated HFs and epithelial proliferation at anagen onset; in their absence, Demodex colonization led to increased epithelial proliferation and replacement of gene programs for repair by aberrant inflammation, leading to the loss of barrier function and HF exhaustion. Humans with rhinophymatous acne rosacea, an inflammatory condition associated with Demodex, had increased HF inflammation with decreased type 2 cytokines, consistent with the inverse relationship seen in mice. Our studies uncover a key role for skin ILC2s and IL-13, which comprise an immune checkpoint that sustains cutaneous integrity and restricts pathologic infestation by colonizing HF mites.

Therapeutic Targeting of Macrophage Plasticity Remodels the Tumor-Immune Microenvironment

SIGNIFICANCE

Comprehensive single-cell proteomics analyses of lung adenocarcinoma progression reveal the role of tumor-associated macrophages in resistance to PD-1 blockade therapy. See related commentary by Lee et al., p. 2515.

Sotigalimab and/or nivolumab with chemotherapy in first-line metastatic pancreatic cancer: clinical and immunologic analyses from the randomized phase 2 PRINCE trial

Chemotherapy combined with immunotherapy has improved the treatment of certain solid tumors, but effective regimens remain elusive for pancreatic ductal adenocarcinoma (PDAC). We conducted a randomized phase 2 trial evaluating the efficacy of nivolumab (nivo; anti-PD-1) and/or sotigalimab (sotiga; CD40 agonistic antibody) with gemcitabine/nab-paclitaxel (chemotherapy) in patients with first-line metastatic PDAC ( NCT03214250 ). In 105 patients analyzed for efficacy, the primary endpoint of 1-year overall survival (OS) was met for nivo/chemo (57.7%, P = 0.006 compared to historical 1-year OS of 35%, n = 34) but was not met for sotiga/chemo (48.1%, P = 0.062, n = 36) or sotiga/nivo/chemo (41.3%, P = 0.223, n = 35). Secondary endpoints were progression-free survival, objective response rate, disease control rate, duration of response and safety. Treatment-related adverse event rates were similar across arms. Multi-omic circulating and tumor biomarker analyses identified distinct immune signatures associated with survival for nivo/chemo and sotiga/chemo. Survival after nivo/chemo correlated with a less suppressive tumor microenvironment and higher numbers of activated, antigen-experienced circulating T cells at baseline. Survival after sotiga/chemo correlated with greater intratumoral CD4 T cell infiltration and circulating differentiated CD4 T cells and antigen-presenting cells. A patient subset benefitting from sotiga/nivo/chemo was not identified. Collectively, these analyses suggest potential treatment-specific correlates of efficacy and may enable biomarker-selected patient populations in subsequent PDAC chemoimmunotherapy trials.

A cytotoxic-skewed immune set point predicts low neutralizing antibody levels after Zika virus infection

Although generating high neutralizing antibody levels is a key component of protective immunity after acute viral infection or vaccination, little is known about why some individuals generate high versus low neutralizing antibody titers. Here, we leverage the high-dimensional single-cell profiling capacity of mass cytometry to characterize the longitudinal cellular immune response to Zika virus (ZIKV) infection in viremic blood donors in Puerto Rico. During acute ZIKV infection, we identify widely coordinated responses across innate and adaptive immune cell lineages. High frequencies of multiple activated cell types during acute infection are associated with high titers of ZIKV neutralizing antibodies 6 months post-infection, while stable immune features suggesting a cytotoxic-skewed immune set point are associated with low titers. Our study offers insight into the coordination of immune responses and identifies candidate cellular biomarkers that may offer predictive value in vaccine efficacy trials aimed at inducing high levels of antiviral neutralizing antibodies.

Differential immune landscapes in appendicular versus axial skeleton

Roughly 400,000 people in the U.S. are living with bone metastases, the vast majority occurring in the spine. Metastases to the spine result in fractures, pain, paralysis, and significant health care costs. This predilection for cancer to metastasize to the bone is seen across most cancer histologies, with the greatest incidence seen in prostate, breast, and lung cancer. The molecular process involved in this predilection for axial versus appendicular skeleton is not fully understood, although it is likely that a combination of tumor and local micro-environmental factors plays a role. Immune cells are an important constituent of the bone marrow microenvironment and many of these cells have been shown to play a significant role in tumor growth and progression in soft tissue and bone disease. With this in mind, we sought to examine the differences in immune landscape between axial and appendicular bones in the normal noncancerous setting in order to obtain an understanding of these landscapes. To accomplish this, we utilized mass cytometry by time-of-flight (CyTOF) to examine differences in the immune cell landscapes between the long bone and vertebral body bone marrow from patient clinical samples and C57BL/6J mice. We demonstrate significant differences between immune populations in both murine and human marrow with a predominance of myeloid progenitor cells in the spine. Additionally, cytokine analysis revealed differences in concentrations favoring a more myeloid enriched population of cells in the vertebral body bone marrow. These differences could have clinical implications with respect to the distribution and permissive growth of bone metastases.

Variation of Immune Cell Responses in Humans Reveals Sex-Specific Coordinated Signaling Across Cell Types

Assessing the health and competence of the immune system is central to evaluating vaccination responses, autoimmune conditions, cancer prognosis, and treatment. With an increasing number of studies examining immune dysregulation, there is a growing need for a curated reference of variation in immune parameters in healthy individuals. We used mass cytometry (CyTOF) to profile blood from 86 humans in response to 15 ex vivo immune stimuli. We present reference ranges for cell-specific immune markers and highlight differences that appear across sex and age. We identified modules of immune features that suggest there exists an underlying structure to the immune system based on signaling pathway responses across cell types. We observed increased MAPK signaling in inflammatory pathways in innate immune cells and greater overall coordination of immune cell responses in females. In contrast, males exhibited stronger pSTAT1 and pTBK1 responses. These reference data are publicly available as a resource for immune profiling studies.

CyTOF® for the Masses

Mass cytometry has revolutionized immunophenotyping, particularly in exploratory settings where simultaneous breadth and depth of characterization of immune populations is needed with limited samples such as in preclinical and clinical tumor immunotherapy. Mass cytometry is also a powerful tool for single-cell immunological assays, especially for complex and simultaneous characterization of diverse intratumoral immune subsets or immunotherapeutic cell populations. Through the elimination of spectral overlap seen in optical flow cytometry by replacement of fluorescent labels with metal isotopes, mass cytometry allows, on average, robust analysis of 60 individual parameters simultaneously. This is, however, associated with significantly increased complexity in the design, execution, and interpretation of mass cytometry experiments. To address the key pitfalls associated with the fragmentation, complexity, and analysis of data in mass cytometry for immunologists who are novices to these techniques, we have developed a comprehensive resource guide. Included in this review are experiment and panel design, antibody conjugations, sample staining, sample acquisition, and data pre-processing and analysis. Where feasible multiple resources for the same process are compared, allowing researchers experienced in flow cytometry but with minimal mass cytometry expertise to develop a data-driven and streamlined project workflow. It is our hope that this manuscript will prove a useful resource for both beginning and advanced users of mass cytometry.

The Route of Vaccine Administration Determines Whether Blood Neutrophils Undergo Long-Term Phenotypic Modifications

Innate immunity modulates adaptive immunity and defines the magnitude, quality, and longevity of antigen-specific T- and B- cell immune memory. Various vaccine and administration factors influence the immune response to vaccination, including the route of vaccine delivery. We studied the dynamics of innate cell responses in blood using a preclinical model of non-human primates immunized with a live attenuated vaccinia virus, a recombinant Modified vaccinia virus Ankara (MVA) expressing a gag-pol-nef fusion of HIV-1, and mass cytometry. We previously showed that it induces a strong, early, and transient innate response, but also late phenotypic modifications of blood myeloid cells after two months when injected subcutaneously. Here, we show that the early innate effector cell responses and plasma inflammatory cytokine profiles differ between subcutaneous and intradermal vaccine injection. Additionally, we show that the intradermal administration fails to induce more highly activated/mature neutrophils long after immunization, in contrast to subcutaneous administration. Different batches of antibodies, staining protocols and generations of mass cytometers were used to generate the two datasets. Mass cytometry data were analyzed in parallel using the same analytical pipeline based on three successive clustering steps, including SPADE, and categorical heatmaps were compared using the Manhattan distance to measure the similarity between cell cluster phenotypes. Overall, we show that the vaccine per se is not sufficient for the late phenotypic modifications of innate myeloid cells, which are evocative of innate immune training. Its route of administration is also crucial, likely by influencing the early innate response, and systemic inflammation, and vaccine biodistribution.

Twin study reveals non-heritable immune perturbations in multiple sclerosis

Multiple sclerosis (MS) is a chronic inflammatory disorder of the central nervous system underpinned by partially understood genetic risk factors and environmental triggers and their undefined interactions1,2. Here we investigated the peripheral immune signatures of 61 monozygotic twin pairs discordant for MS to dissect the influence of genetic predisposition and environmental factors. Using complementary multimodal high-throughput and high-dimensional single-cell technologies in conjunction with data-driven computational tools, we identified an inflammatory shift in a monocyte cluster of twins with MS, coupled with the emergence of a population of IL-2 hyper-responsive transitional naive helper T cells as MS-related immune alterations. By integrating data on the immune profiles of healthy monozygotic and dizygotic twin pairs, we estimated the variance in CD25 expression by helper T cells displaying a naive phenotype to be largely driven by genetic and shared early environmental influences. Nonetheless, the expanding helper T cells of twins with MS, which were also elevated in non-twin patients with MS, emerged independent of the individual genetic makeup. These cells expressed central nervous system-homing receptors, exhibited a dysregulated CD25-IL-2 axis, and their proliferative capacity positively correlated with MS severity. Together, our matched-pair analysis of the extended twin approach allowed us to discern genetically and environmentally determined features of an MS-associated immune signature.

Revealing new biology from multiplexed, metal-isotope-tagged, single-cell readouts

Mass cytometry (MC) is a recent technology that pairs plasma-based ionization of cells in suspension with time-of-flight (TOF) mass spectrometry to sensitively quantify the single-cell abundance of metal-isotope-tagged affinity reagents to key proteins, RNA, and peptides. Given the ability to multiplex readouts (~50 per cell) and capture millions of cells per experiment, MC offers a robust way to assay rare, transitional cell states that are pertinent to human development and disease. Here, we review MC approaches that let us probe the dynamics of cellular regulation across multiple conditions and sample types in a single experiment. Additionally, we discuss current limitations and future extensions of MC as well as computational tools commonly used to extract biological insight from single-cell proteomic datasets.

Single-Cell Profiling of the Intrinsic Apoptotic Pathway by Mass Cytometry (CyTOF)

Mass cytometry time-of-flight (CyTOF) is a technology for the study of complex biological processes at the single-cell level. The technology enables measurement of >50 protein moieties on the surface and inside the cell. The power of CyTOF lies in the application of purpose-built panels of antibody probes that resolve features of key biological processes in a cell. Here, we describe this technology's use to profile changes in the intrinsic apoptotic (cell death) protein machinery at a single-cell level. We provide a comprehensive overview of a tailor-made set of cell survival/death antibodies, ideal staining conditions, and high-dimensional data analysis.

Mass Synaptometry: Applying Mass Cytometry to Single Synapse Analysis

Synaptic degeneration is one of the earliest and phenotypically most significant features associated with numerous neurodegenerative conditions, including Alzheimer's and Parkinson's diseases. Synaptic changes are also known to be important in neurocognitive disorders such as schizophrenia and autism spectrum disorders. Several labs, including ours, have demonstrated that conventional (fluorescence-based) flow cytometry of individual synaptosomes is a robust and reproducible method. However, the repertoire of probes needed to assess comprehensively the type of synapse, pathologic proteins (including protein products of risk genes discovered in GWAS), and markers of stress and injury far exceeds what is achievable with conventional flow cytometry. We recently developed a method that applies CyTOF (Cytometry by Time-Of-Flight mass spectrometry) to high-dimensional analysis of individual human synaptosomes, overcoming many of the multiplexing limitations of conventional flow cytometry. We call this new method Mass Synaptometry. Here we describe the preparation of synaptosomes from human and mouse brain, the generation and quality control of the "SynTOF" (Synapse by Time-Of-Flight mass spectrometry) antibody panel, the staining protocol, and CyTOF parameter setup for acquisition, post-acquisition processing, and analysis.

Single-cell profiling of immune system alterations in lymphoid, barrier and solid tissues in aged mice

Aging exerts profound and paradoxical effects on the immune system, at once impairing proliferation, cytotoxicity and phagocytosis, and inducing chronic inflammation. Previous studies have focused on individual tissues or cell types, while a comprehensive multisystem study of tissue-resident and circulating immune populations during aging is lacking. Here we reveal an atlas of age-related changes in the abundance and phenotype of immune cell populations across 12 mouse tissues. Using cytometry-based high parametric analysis of 37 mass-cytometry and 55 spectral flow-cytometry parameters, mapping samples from young and aged animals revealed conserved and tissue-type-specific patterns of both immune atrophy and expansion. We uncovered clear phenotypic changes in both lymphoid and myeloid lineages in aged mice, and in particular a contraction in natural killer cells and plasmacytoid dendritic cells. These changes correlated with a skewing towards myelopoiesis at the expense of early lymphocyte genesis in aged mice. Taken together, this atlas represents a comprehensive, systematic and thorough resource of the age-dependent alterations of the mammalian immune system in lymphoid, barrier and solid tissues.

Conditional Deletion of Pdcd1 Identifies the Cell-Intrinsic Action of PD-1 on Functional CD8 T Cell Subsets for Antitumor Efficacy

Programmed cell death-1 (PD-1) blockade has a profound effect on the ability of the immune system to eliminate tumors, but many questions remain about the cell types involved and the underlying mechanisms of immune activation. To shed some light on this, the cellular and molecular events following inhibition of PD-1 signaling was investigated in the MC-38 colon carcinoma model using constitutive (PD-1 KO) and conditional (PD1cKO) mice and in wild-type mice treated with PD-1 antibody. The impact on both tumor growth and the development of tumor immunity was assessed. In the PD-1cKO mice, a complete deletion of Pdcd1 in tumor-infiltrating T cells (TILs) after tamoxifen treatment led to the inhibition of tumor growth of both small and large tumors. Extensive immune phenotypic analysis of the TILs by flow and mass cytometry identified 20-different T cell subsets of which specifically 5-CD8 positive ones expanded in all three models after PD-1 blockade. All five subsets expressed granzyme B and interferon gamma (IFNγ). Gene expression analysis of the tumor further supported the phenotypic analysis in both PD-1cKO- and PD-1 Ab-treated mice and showed an upregulation of pathways related to CD4 and CD8 T-cell activation, enhanced signaling through costimulatory molecules and IFNγ, and non-T-cell processes. Altogether, using PD-1cKO mice, we define the intrinsic nature of PD-1 suppression of CD8 T-cell responses in tumor immunity.

Recent Advances in Developmental Hematopoiesis: Diving Deeper With New Technologies

The journey of a hematopoietic stem cell (HSC) involves the passage through successive anatomical sites where HSCs are in direct contact with their surrounding microenvironment, also known as niche. These spatial and temporal cellular interactions throughout development are required for the acquisition of stem cell properties, and for maintaining the HSC pool through balancing self-renewal, quiescence and lineage commitment. Understanding the context and consequences of these interactions will be imperative for our understanding of HSC biology and will lead to the improvement of in vitro production of HSCs for clinical purposes. The aorta-gonad-mesonephros (AGM) region is in this light of particular interest since this is the cradle of HSC emergence during the embryonic development of all vertebrate species. In this review, we will focus on the developmental origin of HSCs and will discuss the novel technological approaches and recent progress made to identify the cellular composition of the HSC supportive niche and the underlying molecular events occurring in the AGM region.

Layilin Anchors Regulatory T Cells in Skin

Regulatory T cells (Tregs) reside in nonlymphoid tissues where they carry out unique functions. The molecular mechanisms responsible for Treg accumulation and maintenance in these tissues are relatively unknown. Using an unbiased discovery approach, we identified LAYN (layilin), a C-type lectin-like receptor, to be preferentially and highly expressed on a subset of activated Tregs in healthy and diseased human skin. Expression of layilin on Tregs was induced by TCR-mediated activation in the presence of IL-2 or TGF-β. Mice with a conditional deletion of layilin in Tregs had reduced accumulation of these cells in tumors. However, these animals somewhat paradoxically had enhanced immune regulation in the tumor microenvironment, resulting in increased tumor growth. Mechanistically, layilin expression on Tregs had a minimal effect on their activation and suppressive capacity in vitro. However, expression of this molecule resulted in a cumulative anchoring effect on Treg dynamic motility in vivo. Taken together, our results suggest a model whereby layilin facilitates Treg adhesion in skin and, in doing so, limits their suppressive capacity. These findings uncover a unique mechanism whereby reduced Treg motility acts to limit immune regulation in nonlymphoid organs and may help guide strategies to exploit this phenomenon for therapeutic benefit.

Towards high throughput and high information coverage: advanced single-cell mass spectrometric techniques

Mass spectrometry (MS) is attractive for single-cell analysis because of its high sensitivity, rich information, and large dynamic ranges, especially for the single-cell metabolome and proteome analysis. Efforts have been made to deal with the throughput and information coverage problems in typical manual single-cell MS techniques. In this review, advanced techniques to improve the automation and throughput for single-cell sampling and single-cell metabolome and proteome MS detection have been discussed. Furthermore, representative MS-based strategies that can increase the in-depth cellular information coverage and achieve the more comprehensive single-cell multiomics information during high throughput detection have been highlighted, providing an ongoing perspective of the MS performance for the single-cell research.

Making the most of high-dimensional cytometry data

High-dimensional cytometry represents an exciting new era of immunology research, enabling the discovery of new cells and prediction of patient responses to therapy. A plethora of analysis and visualization tools and programs are now available for both new and experienced users; however, the transition from low- to high-dimensional cytometry requires a change in the way users think about experimental design and data analysis. Data from high-dimensional cytometry experiments are often underutilized, because of both the size of the data and the number of possible combinations of markers, as well as to a lack of understanding of the processes required to generate meaningful data. In this article, we explain the concepts behind designing high-dimensional cytometry experiments and provide considerations for new and experienced users to design and carry out high-dimensional experiments to maximize quality data collection.

Human Systemic Immune Response to Ingestion of the Oral Probiotic Streptococcus salivarius BLIS K12

Streptococcus salivarius K12 is an oral probiotic known to contribute to protection against oral pathogenic bacteria in humans. Studies of immune responses to S. salivarius K12 have focused on the oral cavity, and systemic immune responses have not yet been reported. The aim of this study was to identify acute systemic immune responses to the commercial product, S. salivarius BLIS K12, in a double-blinded, placebo-controlled human clinical trial. It was hypothesised that consumption of S. salivarius BLIS K12 would induce an anti-inflammatory response and a decrease in pro-inflammatory cytokines. Blood samples were obtained from participants prior to a single dose of S. salivarius BLIS K12 or a placebo and then secondary blood samples were obtained 24 h and 7 days post-consumption. Samples were analysed using multi-parametric flow cytometry, to quantify immune cell frequency changes, and by a LEGENDplex assay of human inflammatory cytokines. Consumption of S. salivarius BLIS K12 was associated with increased levels of IL-8 at 24 h. The frequency of Tregs increased in samples taken 7 days after probiotic consumption, and IL-10 concentrations were higher at 7 days than 24 h after consumption. There was no difference in the frequency and/or activation of CD4⁺ T cells, CD8⁺ T cells, B cells and NK cells. Interestingly, there was an increase in IL-12, 7 days after the consumption of S. salivarius BLIS K12. Collectively, this research demonstrates that ingestion of the probiotic S. salivarius K12 can induce changes in the systemic immune response. The implications of the generation and type of immune response warrant further study to determine potential health benefits.

Oropharyngeal Squamous Cell Carcinoma Treatment in the Era of Immune Checkpoint Inhibitors

While head and neck squamous cell carcinomas (HNSCC) are marginally decreasing due to the reduction in exposure to the major risk factors, tobacco and alcohol, the incidence of high-risk human papillomavirus (HPV)-positive oropharynx squamous cell carcinomas (OPSCC), especially those in the tonsil and base of tongue subsites, are increasing. Patients with the latter are younger, display a longer overall survival, and show a lower recurrence rate after standard-of-care treatment than those with HPV-negative OPSCC. This may reflect an important role for immune surveillance and control during the natural history of the virally driven tumour development. Immune deviation through acquisition of immune-suppressive factors in the tumour microenvironment (TME) is discussed in relation to treatment response. Understanding how the different immune factors are integrated in the TME battleground offers opportunities for identifying prognostic biomarkers as well as novel therapeutic strategies. OPSCC generally receive surgery or radiotherapy for early-stage tumour treatment, but many patients present with locoregionally advanced disease requiring multimodality therapies which can involve considerable complications. This review focuses on the utilization of newly emerged immune checkpoint inhibitors (PD-1/PD-L1 pathway) for treatment of HNSCC, in particular HPV-positive OPSCC, since they could be less toxic and more efficacious. PD-1/PD-L1 expression in the TME has been extensively investigated as a biomarker of patient response but is yet to provide a really effective means for stratification of treatment. Extensive testing of combinations of therapeutic approaches by types and sequencing will fuel the next evolution of treatment for OPSCC.

Single-cell profiling of myasthenia gravis identifies a pathogenic T cell signature

Myasthenia gravis (MG) is an autoimmune disease characterized by impaired neuromuscular signaling due to autoantibodies targeting the acetylcholine receptor. Although its auto-antigens and effector mechanisms are well defined, the cellular and molecular drivers underpinning MG remain elusive. Here, we employed high-dimensional single-cell mass and spectral cytometry of blood and thymus samples from MG patients in combination with supervised and unsupervised machine-learning tools to gain insight into the immune dysregulation underlying MG. By creating a comprehensive immune map, we identified two dysregulated subsets of inflammatory circulating memory T helper (Th) cells. These signature Th_CD103 and Th_GM cells populated the diseased thymus, were reduced in the blood of MG patients, and were inversely correlated with disease severity. Both signature Th subsets rebounded in the blood of MG patients after surgical thymus removal, indicative of their role as cellular markers of disease activity. Together, this in-depth analysis of the immune landscape of MG provides valuable insight into disease pathogenesis, suggests novel biomarkers and identifies new potential therapeutic targets for treatment.

Spatially visualized single-cell pathology of highly multiplexed protein profiles in health and disease

Deep molecular profiling of biological tissues is an indicator of health and disease. We used imaging mass cytometry (IMC) to acquire spatially resolved 20-plex protein data in tissue sections from normal and chronic tonsillitis cases. We present SpatialViz, a suite of algorithms to explore spatial relationships in multiplexed tissue images by visualizing and quantifying single-cell granularity and anatomical complexity in diverse multiplexed tissue imaging data. Single-cell and spatial maps confirmed that CD68+ cells were correlated with the enhanced Granzyme B expression and CD3+ cells exhibited enrichment of CD4+ phenotype in chronic tonsillitis. SpatialViz revealed morphological distributions of cellular organizations in distinct anatomical areas, spatially resolved single-cell associations across anatomical categories, and distance maps between the markers. Spatial topographic maps showed the unique organization of different tissue layers. The spatial reference framework generated network-based comparisons of multiplex data from healthy and diseased tonsils. SpatialViz is broadly applicable to multiplexed tissue biology.

SMNN: batch effect correction for single-cell RNA-seq data via supervised mutual nearest neighbor detection

Batch effect correction has been recognized to be indispensable when integrating single-cell RNA sequencing (scRNA-seq) data from multiple batches. State-of-the-art methods ignore single-cell cluster label information, but such information can improve the effectiveness of batch effect correction, particularly under realistic scenarios where biological differences are not orthogonal to batch effects. To address this issue, we propose SMNN for batch effect correction of scRNA-seq data via supervised mutual nearest neighbor detection. Our extensive evaluations in simulated and real datasets show that SMNN provides improved merging within the corresponding cell types across batches, leading to reduced differentiation across batches over MNN, Seurat v3 and LIGER. Furthermore, SMNN retains more cell-type-specific features, partially manifested by differentially expressed genes identified between cell types after SMNN correction being biologically more relevant, with precision improving by up to 841.0%.

Arid5a Promotes Immune Evasion by Augmenting Tryptophan Metabolism and Chemokine Expression

The acquisition of mesenchymal traits leads to immune evasion in various cancers, but the underlying molecular mechanisms remain unclear. In this study, we found that the expression levels of AT-rich interaction domain-containing protein 5a (Arid5a), an RNA-binding protein, were substantially increased in mesenchymal tumor subtypes. The deletion of Arid5a in tumor cell lines enhanced antitumor immunity in immunocompetent mice, but not in immunodeficient mice, suggesting a role for Arid5a in immune evasion. Furthermore, an Arid5a-deficient tumor microenvironment was shown to have robust antitumor immunity, as manifested by suppressed infiltration of granulocytic myeloid-derived suppressor cells and regulatory T cells. In addition, infiltrated T cells were more cytotoxic and less exhausted. Mechanistically, Arid5a stabilized Ido1 and Ccl2 mRNAs and augmented their expression, resulting in enhanced tryptophan catabolism and an immunosuppressive tumor microenvironment. Thus, our findings demonstrate the role of Arid5a beyond inflammatory diseases and suggest Arid5a as a promising target for the treatment of immunotolerant malignant tumors.See related Spotlight by Van den Eynde, p. 854.

Chromatin Landscape Underpinning Human Dendritic Cell Heterogeneity

Human dendritic cells (DCs) comprise subsets with distinct phenotypic and functional characteristics, but the transcriptional programs that dictate their identity remain elusive. Here, we analyze global chromatin accessibility profiles across resting and stimulated human DC subsets by means of the assay for transposase-accessible chromatin using sequencing (ATAC-seq). We uncover specific regions of chromatin accessibility for each subset and transcriptional regulators of DC function. By comparing plasmacytoid DC responses to IFN-I-producing and non-IFN-I-producing conditions, we identify genetic programs related to their function. Finally, by intersecting chromatin accessibility with genome-wide association studies, we recognize DC subset-specific enrichment of heritability in autoimmune diseases. Our results unravel the basis of human DC subset heterogeneity and provide a framework for their analysis in disease pathogenesis.

Human dendritic cells (DCs) orchestrate immune responses by a division of labor between functionally specialized subsets; however, the transcriptional basis of this heterogeneity is poorly understood. Using ATAC-seq, Leylek et al. profile the chromatin landscape of human DC subsets, providing insight into the underlying regulatory mechanisms that modulate their function.