The availability of single-cell transcriptomics has allowed the construction of reference cell atlases, but their usefulness depends on the quality of dataset integration and the ability to map new samples. Previous benchmarks have compared integration methods and suggest that feature selection improves performance but have not explored how best to select features. Here, we benchmark feature selection methods for single-cell RNA sequencing integration using metrics beyond batch correction and preservation of biological variation to assess query mapping, label transfer and the detection of unseen populations. We reinforce common practice by showing that highly variable feature selection is effective for producing high-quality integrations and provide further guidance on the effect of the number of features selected, batch-aware feature selection, lineage-specific feature selection and integration and the interaction between feature selection and integration models. These results are informative for analysts working on large-scale tissue atlases, using atlases or integrating their own data to tackle specific biological questions.

Abstract Central nervous system tumors are the most common solid tumors and the leading cause of cancer-related deaths in children. Outcomes remain dismal for a significant proportion of patients despite aggressive multimodal therapy and survivors suffer from treatment-related morbidity and neurological deficits. Despite their remarkable success in adult cancers, immunotherapies are not approved for pediatric brain tumors. Several factors like the relatively smaller patient population, a tumor location that makes them difficult to access through biopsies and poorly-studied tumor immune microenvironment (TME) have posed challenges in advancing immunotherapeutic options for these patients. Immune checkpoint blockers (ICB) tested in unstratified cohorts of children with brain tumors have not shown promising results; therefore, the role of immunotherapy remains unclear. Because the magnitude and nature of pre-existing T cell responses in the tumor has been shown to correlate with response to ICB, we asked if pediatric brain tumors show evidence of such favorable T cell responses. In the current study, we address this question using an unbiased approach to analyze the clonal relationship and transcriptional status of T cells present in tumor tissue of children with brain tumors using paired single-cell RNA-sequencing and TCR-sequencing. Our results demonstrate that a large fraction of T cells in the tumor tissue are clonally expanded with potential to recognize tumor antigens. Such clonally-expanded T cells display states linked to anti-tumor immunity, express higher levels of transcripts encoding for molecules linked to T cell activation, effector functions, immune cell recruitment, tissue-residency, immune checkpoints, and importantly, show significant enrichment of signatures linked to neoantigen-specific T cells and immunotherapy response. Notably, we identify several neoantigens in pediatric brain tumors, and show that neoantigen-specific T cell gene signatures are linked to better survival outcomes in high-grade glioma. We further show that PD1-expressing CD8+ T cells in pediatric brain tumors are indeed functional as evidenced by their capacity for cytotoxicity, cytokine production and proliferation. Among the patients in our cohort, we observe substantial heterogeneity in the degree of clonal expansion and expression of transcripts encoding immune checkpoints in tumor-infiltrating T cells. Overall, our findings suggest that accurate characterization of intra-tumoral T cell responses for features linked to immune check point therapy (ICB) response may enable selection of patients where immunotherapy is likely to be beneficial, an approach that requires prospective validation in future clinical trials. Citation Format: Anusha Preethi Ganesan, Aditi Upadhye, Kevin Meza Landeros, Ciro Ramírez-Suástegui1, Benjamin Schmiedel, Edwin Woo, Serena J. Chee, Denise Malicki, Nicole Coufal, David Gonda, Michael Levy, Jason Greenbaum, Gregory Seumois, John Crawford, William Roberts, Stephen Schoenberger, Hilde Cheroutre, Christian Ottensmeier, Pandurangan Vijayanand. Single-cell transcriptomics and TCR clonality reveal a role for immunotherapy in pediatric brain tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB310.

Brain tumors in children are a devastating disease in a high proportion of patients. Owing to inconsistent results in clinical trials in unstratified patients, the role of immunotherapy remains unclear. We performed an in-depth survey of the single-cell transcriptomes and clonal relationship of intra-tumoral T cells from children with brain tumors. Our results demonstrate that a large fraction of T cells in the tumor tissue are clonally expanded with the potential to recognize tumor antigens. Such clonally expanded T cells display enrichment of transcripts linked to effector function, tissue residency, immune checkpoints and signatures of neoantigen-specific T cells and immunotherapy response. We identify neoantigens in pediatric brain tumors and show that neoantigen-specific T cell gene signatures are linked to better survival outcomes. Notably, among the patients in our cohort, we observe substantial heterogeneity in the degree of clonal expansion and magnitude of T cell response. Our findings suggest that characterization of intra-tumoral T cell responses may enable selection of patients for immunotherapy, an approach that requires prospective validation in clinical trials. Ganesan and colleagues characterize T cell clonality and transcriptomes at the single-cell level in pediatric brain tumor samples, providing insights into existing tumor neoantigens and T cell responses and the potential for effective immunotherapy.

The rapid adoption of single-cell technologies has created an opportunity to build single-cell ‘atlases’ integrating diverse datasets across many laboratories. Such atlases can serve as a reference for analyzing and interpreting current and future data. However, it has become apparent that atlasing approaches differ, and the impact of these differences are often unclear. Here we review the current atlasing literature and present considerations for building and using atlases. Importantly, we find that no one-size-fits-all protocol for atlas building exists, but rather we discuss context-specific considerations and workflows, including atlas conceptualization, data collection, curation and integration, atlas evaluation and atlas sharing. We further highlight the benefits of integrated atlases for analyses of new datasets and deriving biological insights beyond what is possible from individual datasets. Our overview of current practices and associated recommendations will improve the quality of atlases to come, facilitating the shift to a unified, reference-based understanding of single-cell biology. This Review provides a comprehensive and detailed discussion about how to build and use single-cell atlases.

To understand the role of T cells in the pathogenesis of ulcerative colitis (UC), we analyzed colonic T cells isolated from patients with UC and controls. Here we identified colonic CD4+ and CD8+ T lymphocyte subsets with gene expression profiles resembling stem-like progenitors, previously reported in several mouse models of autoimmune disease. Stem-like T cells were increased in inflamed areas compared to non-inflamed regions from the same patients. Furthermore, TCR sequence analysis indicated stem-like T cells were clonally related to proinflammatory T cells, suggesting their involvement in sustaining effectors that drive inflammation. Using an adoptive transfer colitis model in mice, we demonstrated that CD4+ T cells deficient in either BCL-6 or TCF1, transcription factors that promote T cell stemness, had decreased colon T cells and diminished pathogenicity. Our results establish a strong association between stem-like T cell populations and UC pathogenesis, highlighting the potential of targeting this population to improve clinical outcomes. Li et al. identify a human TCF1-expressing stem-like T cell population that is associated with immune-mediated pathogenesis in ulcerative colitis.

ABSTRACTOrgan- and body-scale cell atlases have the potential to transform our understanding of human biology. To capture the variability present in the population, these atlases must include diverse demographics such as age and ethnicity from both healthy and diseased individuals. The growth in both size and number of single-cell datasets, combined with recent advances in computational techniques, for the first time makes it possible to generate such comprehensive large-scale atlases through integration of multiple datasets. Here, we present the integrated Human Lung Cell Atlas (HLCA) combining 46 datasets of the human respiratory system into a single atlas spanning over 2.2 million cells from 444 individuals across health and disease. The HLCA contains a consensus re-annotation of published and newly generated datasets, resolving under- or misannotation of 59% of cells in the original datasets. The HLCA enables recovery of rare cell types, provides consensus marker genes for each cell type, and uncovers gene modules associated with demographic covariates and anatomical location within the respiratory system. To facilitate the use of the HLCA as a reference for single-cell lung research and allow rapid analysis of new data, we provide an interactive web portal to project datasets onto the HLCA. Finally, we demonstrate the value of the HLCA reference for interpreting disease-associated changes. Thus, the HLCA outlines a roadmap for the development and use of organ-scale cell atlases within the Human Cell Atlas.

Patients with severe uncontrolled asthma represent a distinct endotype with persistent airway inflammation and remodeling that is refractory to corticosteroid treatment. To determine T cell subsets and effector molecules that drive pathogenesis of severe asthma, we performed single-cell transcriptome analysis of >50,000 airway CD4+ T cells isolated from bronchoalveolar lavage (BAL) samples from 30 patients with mild and severe asthma. We observed striking heterogeneity in the nature of CD4+ T cells present in asthmatics' airways with tissue-resident memory (TRM) cells making a dominant contribution. Notably, in severe asthmatics a subset of CD4+ TRM cells (CD103-expressing) was significantly increased, comprising nearly 65% of all CD4+ T cells in the airways of male patients with severe asthma when compared to mild asthma (13%). This subset was enriched for transcripts linked to T cell receptor (TCR) activation (HLA-DRB1, HLA-DPA1, CD40LG) and cytotoxicity (GZMB, GZMH), and following stimulation expressed high levels of transcripts encoding for pro-inflammatory non-TH2 cytokines (CCL3, CCL4, CCL5, TNF, LIGHT) that could fuel persistent airway inflammation and remodeling. Our findings indicate the need to look beyond the traditional T2 model of severe asthma to better understand the heterogeneity of this disease. This work was supported by (i) NIH research grant R01HL114093 (PV) and equipment grants (S10RR027366 - BD FACSAria Fusion, and S10OD025052 - Illumina Novaseq 6000); the William K. Bowes Jr. Foundation (P.V.). The WATCH study is supported by the Southampton NIHR Biomedical Research Centre and the Southampton NIHR Clinical Research Facility which are funded by the NIHR and are a partnership between the University of Southampton and University Hospital Southampton NHS Foundation Trust.

Mucosal-associated invariant T (MAIT) cells are a subset of T lymphocytes that respond to microbial metabolites. We defined MAIT cell populations in different organs and characterized the developmental pathway of mouse and human MAIT cells in the thymus using single-cell RNA sequencing and phenotypic and metabolic analyses. We showed that the predominant mouse subset, which produced IL-17 (MAIT17), and the subset that produced IFN-γ (MAIT1) had not only greatly different transcriptomes but also different metabolic states. MAIT17 cells in different organs exhibited increased lipid uptake, lipid storage, and mitochondrial potential compared with MAIT1 cells. All these properties were similar in the thymus and likely acquired there. Human MAIT cells in lung and blood were more homogeneous but still differed between tissues. Human MAIT cells had increased fatty acid uptake and lipid storage in blood and lung, similar to human CD8 T resident memory cells, but unlike mouse MAIT17 cells, they lacked increased mitochondrial potential. Although mouse and human MAIT cell transcriptomes showed similarities for immature cells in the thymus, they diverged more strikingly in the periphery. Analysis of pet store mice demonstrated decreased lung MAIT17 cells in these so-called "dirty" mice, indicative of an environmental influence on MAIT cell subsets and function.

Junctional adhesion molecule-like protein (JAML) serves as a co-stimulatory molecule in γδ T cells. While it has recently been described as a cancer immunotherapy target in mice, its potential to cause toxicity, specific mode of action with regard to its cellular targets, and whether it can be targeted in humans remain unknown. Here, we show that JAML is induced by T cell receptor engagement, reveal that this induction is linked to cis-regulatory interactions between the CD3D and JAML gene loci. When compared with other immunotherapy targets plagued by low target specificity and end-organ toxicity, we find JAML to be mostly restricted to and highly expressed by tissue-resident memory CD8+ T cells in multiple cancer types. By delineating the key cellular targets and functional consequences of agonistic anti-JAML therapy in a murine melanoma model, we show its specific mode of action and the reason for its synergistic effects with anti-PD-1.

Phosphoinositide 3-kinase δ (PI3Kδ) has a key role in lymphocytes, and inhibitors that target this PI3K have been approved for treatment of B cell malignancies 1–3 . Although studies in mouse models of solid tumours have demonstrated that PI3Kδ inhibitors (PI3Kδi) can induce anti-tumour immunity 4,5 , its effect on solid tumours in humans remains unclear. Here we assessed the effects of the PI3Kδi AMG319 in human patients with head and neck cancer in a neoadjuvant, double-blind, placebo-controlled randomized phase II trial (EudraCT no. 2014-004388-20). PI3Kδ inhibition decreased the number of tumour-infiltrating regulatory T (T reg ) cells and enhanced the cytotoxic potential of tumour-infiltrating T cells. At the tested doses of AMG319, immune-related adverse events (irAEs) required treatment to be discontinued in 12 out of 21 of patients treated with AMG319, suggestive of systemic effects on T reg cells. Accordingly, in mouse models, PI3Kδi decreased the number of T reg cells systemically and caused colitis. Single-cell RNA-sequencing analysis revealed a PI3Kδi-driven loss of tissue-resident colonic ST2 T reg cells, accompanied by expansion of pathogenic T helper 17 (T H 17) and type 17 CD8 + T (T C 17) cells, which probably contributed to toxicity; this points towards a specific mode of action for the emergence of irAEs. A modified treatment regimen with intermittent dosing of PI3Kδi in mouse models led to a significant decrease in tumour growth without inducing pathogenic T cells in colonic tissue, indicating that alternative dosing regimens might limit toxicity.

The impact of genetic variants on cells challenged in biologically relevant contexts has not been fully explored. Here, we activated CD4+ T cells from 89 healthy donors and performed a single-cell RNA sequencing assay with >1 million cells to examine cell type-specific and activation-dependent effects of genetic variants. Single-cell expression quantitative trait loci (sc-eQTL) analysis of 19 distinct CD4+ T cell subsets showed that the expression of over 4000 genes is significantly associated with common genetic polymorphisms and that most of these genes show their most prominent effects in specific cell types. These genes included many that encode for molecules important for activation, differentiation, and effector functions of T cells. We also found new gene associations for disease-risk variants identified from genome-wide association studies and highlighted the cell types in which their effects are most prominent. We found that biological sex has a major influence on activation-dependent gene expression in CD4+ T cell subsets. Sex-biased transcripts were significantly enriched in several pathways that are essential for the initiation and execution of effector functions by CD4+ T cells like TCR signaling, cytokines, cytokine receptors, costimulatory, apoptosis, and cell-cell adhesion pathways. Overall, this DICE (Database of Immune Cell Expression, eQTLs, and Epigenomics) subproject highlights the power of sc-eQTL studies for simultaneously exploring the activation and cell type-dependent effects of common genetic variants on gene expression (https://dice-database.org).

Parkinson's disease (PD) as a progressive neurodegenerative disorder arises from multiple genetic and environmental factors. However, underlying pathological mechanisms remain poorly understood. Using multiplexed single-cell transcriptomics, we analyze human neural precursor cells (hNPCs) from sporadic PD (sPD) patients. Alterations in gene expression appear in pathways related to primary cilia (PC). Accordingly, in these hiPSC-derived hNPCs and neurons, we observe a shortening of PC. Additionally, we detect a shortening of PC in PINK1-deficient human cellular and mouse models of familial PD. Furthermore, in sPD models, the shortening of PC is accompanied by increased Sonic Hedgehog (SHH) signal transduction. Inhibition of this pathway rescues the alterations in PC morphology and mitochondrial dysfunction. Thus, increased SHH activity due to ciliary dysfunction may be required for the development of pathoetiological phenotypes observed in sPD like mitochondrial dysfunction. Inhibiting overactive SHH signaling may be a potential neuroprotective therapy for sPD.

NK cells are the predominant innate lymphocyte subset that contribute to anti-tumor immunity. The molecular profile of human tumor-infiltrating NK cells, and particularly, their characteristics in relation to CD8 + cytotoxic T lymphocyte (CTL) density within tumors is poorly understood. Here, we performed single-cell and bulk transcriptomic analysis of NK cells within tumor and adjacent normal lung tissue from patients with treatment-naïve lung cancer. We found divergent NK cell subsets with distinct functional features that were associated with the magnitude of CD8 + CTL response in tumors. CD16 pos NK cells present in tumors with high density of CD8 + CTLs were enriched for transcripts linked to cytotoxicity, pro-inflammatory chemokines and effector function. CD16 neg NK cells predominantly infiltrated tumors relative to lung tissue and harbored sub-populations with unique gene expression programs. Notably, tumors with high density of CD8 + CTLs were enriched for a CD16 neg NK cell subset with features of tissue residency and NK cell memory as well as a CD16 neg ILC3-like subset that displayed transcriptional features linked to activation, NFκB signaling, dendritic cell recruitment and lymphoid-tissue inducing properties. These data reveal that a potentially synergistic NK cell innate immune program may co-operate with CD8 + CTLs to mediate tumor clearance in tumors with a high tumor-infiltrating lymphocyte (TIL) density.

Immune-checkpoint blockade (ICB) has shown remarkable clinical success in boosting antitumor immunity. However, the breadth of its cellular targets and specific mode of action remain elusive. We find that tumor-infiltrating follicular regulatory T (TFR) cells are prevalent in tumor tissues of several cancer types. They are primarily located within tertiary lymphoid structures and exhibit superior suppressive capacity and in vivo persistence as compared with regulatory T cells, with which they share a clonal and developmental relationship. In syngeneic tumor models, anti-PD-1 treatment increases the number of tumor-infiltrating TFR cells. Both TFR cell deficiency and the depletion of TFR cells with anti-CTLA-4 before anti-PD-1 treatment improve tumor control in mice. Notably, in a cohort of 271 patients with melanoma, treatment with anti-CTLA-4 followed by anti-PD-1 at progression was associated with better a survival outcome than monotherapy with anti-PD-1 or anti-CTLA-4, anti-PD-1 followed by anti-CTLA-4 at progression or concomitant combination therapy. Vijayanand and colleagues show highly suppressive CD4+CTLA-4+PD-1+ follicular regulatory T (TFR) cells reside within tumor microenvironments. Depleting TFR cells or blocking their activity with CTLA-4-depleting antibodies before anti-PD-1 checkpoint blockade therapy improved the efficacy of anti-PD-1 treatment in mouse tumor models and was also associated with better survival outcomes in a large cohort of patients with melanoma.

Phosphoinositide 3-kinase δ (PI3Kδ) plays a key role in lymphocytes and inhibitors targeting this PI3K have been approved for hematological malignancies. While studies in hematological and solid tumor models in mice have demonstrated that PI3Kδ inhibitors (PI3Kδi) can induce anti-tumor immunity, the impact of PI3Kδi on solid tumors in humans remains unclear. Here, we assessed the effects of the PI3Kδi AMG319 in patients with resectable head and neck cancer in a neoadjuvant, double-blind, placebo-controlled randomised phase-II trial. We find that PI3Kδ inhibition decreases tumor-infiltrating immunosuppressive TREG cells and causes heightened cytotoxic potential of tumor-infiltrating CD8+ and CD4+ T cells. Loss of intratumoral TREG cells and an increase in the frequency of activated TREG cells in the blood post-treatment are indicative of systemic effects on TREG tissue retention and maintenance. At the tested AMG319 doses, immune-related adverse events caused treatment discontinuation in 12/21 of AMG319-treated patients, further suggestive of systemic effects on TREG cells. Consistent with this notion, in a murine syngeneic tumor model, PI3Kδi decreased TREG cells in both tumor and non-malignant tissues and affected TREG subtype composition, maintenance and functionality. Our data demonstrate the cancer-immunotherapy potential of PI3Kδ inhibition in humans, but its modulation will need to be carefully balanced to harness its anti-tumor capacity while minimizing immune related toxicity.

The molecular properties of CD8+ T cells that respond to SARS-CoV-2 infection are not fully known. Here, we report on the single-cell transcriptomes of >80,000 virus-reactive CD8+ T cells, obtained using a modified Antigen-Reactive T cell Enrichment (ARTE) assay, from 39 COVID-19 patients and 10 healthy subjects. COVID-19 patients segregated into two groups based on whether the dominant CD8+ T cell response to SARS-CoV-2 was 'exhausted' or not. SARS-CoV-2-reactive cells in the exhausted subset were increased in frequency and displayed lesser cytotoxicity and inflammatory features in COVID-19 patients with mild compared to severe illness. In contrast, SARS-CoV-2-reactive cells in the dominant non-exhausted subset from patients with severe disease showed enrichment of transcripts linked to co-stimulation, pro-survival NF-κB signaling, and anti-apoptotic pathways, suggesting the generation of robust CD8+ T cell memory responses in patients with severe COVID-19 illness. CD8+ T cells reactive to influenza and respiratory syncytial virus from healthy subjects displayed polyfunctional features and enhanced glycolysis. Cells with such features were largely absent in SARS-CoV-2-reactive cells from both COVID-19 patients and healthy controls non-exposed to SARS-CoV-2. Overall, our single-cell analysis revealed substantial diversity in the nature of CD8+ T cells responding to SARS-CoV-2.

Immune checkpoint blockade (ICB) has shown remarkable clinical success in boosting anti-tumor immunity. However, the breadth of its cellular targets and specific mode of action remain elusive. We find that T follicular regulatory (TFR) cells are present in high numbers in multiple tumors, inhibit anti-tumor immunity and are responsive to ICB. TCR-seq data, trajectory analyses and adoptive transfer studies indicate intratumoral TREG to TFR cell conversion. When compared to TREG cells, TFR cells exhibited enhanced suppressive capacity. In syngeneic tumor models, anti-PD-1 treatment increased the number of tumor-infiltrating TFR cells. Conditional knockout of TFR cells or depletion of TFR cells with anti-CTLA-4 antibody prior to anti-PD-1 treatment, improved tumor control in mice. Notably, in a cohort of melanoma patients, treatment with anti-CTLA-4 followed by anti-PD-1 at progression was associated with better long-term survival outcomes than anti-PD-1 or anti-CTLA-4 monotherapy, anti-PD-1 followed by anti CTLA-4 at progression or concomitant combination therapy.

The contribution of CD4(+) T cells to protective or pathogenic immune responses to SARS-CoV-2 infection remains unknown. Here, we present single-cell transcriptomic analysis of >100,000 viral antigen-reactive CD4(+) T cells from 40 COVID-19 patients. In hospitalized patients compared to non-hospitalized patients, we found increased proportions of cytotoxic follicular helper cells and cytotoxic T helper (T-H) cells (CD4-CTLs) responding to SARS-CoV-2 and reduced proportion of SARS-CoV-2-reactive regulatory T cells (T-REG). Importantly, in hospitalized COVID-19 patients, a strong cytotoxic T-FH response was observed early in the illness, which correlated negatively with antibody levels to SARS-CoV-2 spike protein. Polyfunctional T(H)1 and T(H)17 cell subsets were underrepresented in the repertoire of SARS-CoV-2-reactive CD4(+) T cells compared to influenza-reactive CD4(+) T cells. Together, our analyses provide insights into the gene expression patterns of SARS-CoV-2-reactive CD4(+) T cells in distinct disease severities.