Heterogeneity of human bone marrow and blood natural killer cells defined by single-cell transcriptome

HIV/HBV coinfection remodels the immune landscape and natural killer cell ADCC functional responses

BACKGROUND AND AIMS

HBV and HIV coinfection is a common occurrence globally, with significant morbidity and mortality. Both viruses lead to immune dysregulation including changes in natural killer (NK) cells, a key component of antiviral defense and a promising target for HBV cure strategies. Here we used high-throughput single-cell analysis to explore the immune cell landscape in people with HBV mono-infection and HIV/HBV coinfection, on antiviral therapy, with emphasis on identifying the distinctive characteristics of NK cell subsets that can be therapeutically harnessed.

APPROACH AND RESULTS

Our data show striking differences in the transcriptional programs of NK cells. HIV/HBV coinfection was characterized by an over-representation of adaptive, KLRC2 -expressing NK cells, including a higher abundance of a chemokine-enriched ( CCL3/CCL4 ) adaptive cluster. The NK cell remodeling in HIV/HBV coinfection was reflected in enriched activation pathways, including CD3ζ phosphorylation and ZAP-70 translocation that can mediate stronger antibody-dependent cellular cytotoxicity responses and a bias toward chemokine/cytokine signaling. By contrast, HBV mono-infection imposed a stronger cytotoxic profile on NK cells and a more prominent signature of "exhaustion" with higher circulating levels of HBsAg. Phenotypic alterations in the NK cell pool in coinfection were consistent with increased "adaptiveness" and better capacity for antibody-dependent cellular cytotoxicity compared to HBV mono-infection. Overall, an adaptive NK cell signature correlated inversely with circulating levels of HBsAg and HBV-RNA in our cohort.

CONCLUSIONS

This study provides new insights into the differential signature and functional profile of NK cells in HBV and HIV/HBV coinfection, highlighting pathways that can be manipulated to tailor NK cell-focused approaches to advance HBV cure strategies in different patient groups.

Pan-cancer profiling of tumor-infiltrating natural killer cells through transcriptional reference mapping

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions.

High-dimensional single-cell analysis of human natural killer cell heterogeneity

Natural killer (NK) cells are innate lymphoid cells (ILCs) contributing to immune responses to microbes and tumors. Historically, their classification hinged on a limited array of surface protein markers. Here, we used single-cell RNA sequencing (scRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to dissect the heterogeneity of NK cells. We identified three prominent NK cell subsets in healthy human blood: NK1, NK2 and NK3, further differentiated into six distinct subgroups. Our findings delineate the molecular characteristics, key transcription factors, biological functions, metabolic traits and cytokine responses of each subgroup. These data also suggest two separate ontogenetic origins for NK cells, leading to divergent transcriptional trajectories. Furthermore, we analyzed the distribution of NK cell subsets in the lung, tonsils and intraepithelial lymphocytes isolated from healthy individuals and in 22 tumor types. This standardized terminology aims at fostering clarity and consistency in future research, thereby improving cross-study comparisons.

Single-cell profiling identifies a CD8bright CD244bright Natural Killer cell subset that reflects disease activity in HLA-A29-positive birdshot chorioretinopathy

Birdshot chorioretinopathy is an inflammatory eye condition strongly associated with MHC-I allele HLA-A29. The striking association with MHC-I suggests involvement of T cells, whereas natural killer (NK) cell involvement remains largely unstudied. Here we show that HLA-A29-positive birdshot chorioretinopathy patients have a skewed NK cell pool containing expanded CD16 positive NK cells which produce more proinflammatory cytokines. These NK cells contain populations that express CD8A which is involved in MHC-I recognition on target cells, display gene signatures indicative of high cytotoxic activity (GZMB, PRF1 and ISG15), and signaling through NK cell receptor CD244 (SH2D1B). Long-term monitoring of a cohort of birdshot chorioretinopathy patients with active disease identifies a population of CD8bright CD244bright NK cells, which rapidly declines to normal levels upon clinical remission following successful treatment. Collectively, these studies implicate CD8bright CD244bright NK cells in birdshot chorioretinopathy.

Myeloid-T cell interplay and cell state transitions associated with checkpoint inhibitor response in melanoma

BACKGROUND

The treatment of melanoma, the deadliest form of skin cancer, has greatly benefited from immunotherapy. However, many patients do not show a durable response, which is only partially explained by known resistance mechanisms.

METHODS

We performed single-cell RNA sequencing of tumor immune infiltrates and matched peripheral blood mononuclear cells of 22 checkpoint inhibitor (CPI)-naive stage III-IV metastatic melanoma patients. After sample collection, the same patients received CPI treatment, and their response was assessed.

FINDINGS

CPI responders showed high levels of classical monocytes in peripheral blood, which preferentially transitioned toward CXCL9-expressing macrophages in tumors. Trajectories of tumor-infiltrating CD8+ T cells diverged at the level of effector memory/stem-like T cells, with non-responder cells progressing into a state characterized by cellular stress and apoptosis-related gene expression. Consistently, predicted non-responder-enriched myeloid-T/natural killer cell interactions were primarily immunosuppressive, while responder-enriched interactions were supportive of T cell priming and effector function.

CONCLUSIONS

Our study illustrates that the tumor immune microenvironment prior to CPI treatment can be indicative of response. In perspective, modulating the myeloid and/or effector cell compartment by altering the described cell interactions and transitions could improve immunotherapy response.

FUNDING

This research was funded by Roche Pharma Research and Early Development.

Comparative single-cell transcriptional and proteomic atlas of clinical-grade injectable mesenchymal source tissues

Bone marrow aspirate concentrate (BMAC) and adipose-derived stromal vascular fraction (ADSVF) are the most marketed stem cell therapies to treat a variety of conditions in the general population and elite athletes. Both tissues have been used interchangeably clinically even though their detailed composition, heterogeneity, and mechanisms of action have neither been rigorously inventoried nor compared. This lack of information has prevented investigations into ideal dosages and has facilitated anecdata and misinformation. Here, we analyzed single-cell transcriptomes, proteomes, and flow cytometry profiles from paired clinical-grade BMAC and ADSVF. This comparative transcriptional atlas challenges the prevalent notion that there is one therapeutic cell type present in both tissues. We also provide data of surface markers that may enable isolation and investigation of cell (sub)populations. Furthermore, the proteome atlas highlights intertissue and interpatient heterogeneity of injected proteins with potentially regenerative or immunomodulatory capacities. An interactive webtool is available online.

HLA and KIR genetic association and NK cells in anti-NMDAR encephalitis

Introduction

Genetic predisposition to autoimmune encephalitis with antibodies against N-methyl-D-aspartate receptor (NMDAR) is poorly understood. Given the diversity of associated environmental factors (tumors, infections), we hypothesized that human leukocyte antigen (HLA) and killer-cell immunoglobulin-like receptors (KIR), two extremely polymorphic gene complexes key to the immune system, might be relevant for the genetic predisposition to anti-NMDAR encephalitis. Notably, KIR are chiefly expressed by Natural Killer (NK) cells, recognize distinct HLA class I allotypes and play a major role in anti-tumor and anti-infection responses.

Methods

We conducted a Genome Wide Association Study (GWAS) with subsequent control-matching using Principal Component Analysis (PCA) and HLA imputation, in a multi-ethnic cohort of anti-NMDAR encephalitis (n=479); KIR and HLA were further sequenced in a large subsample (n=323). PCA-controlled logistic regression was then conducted for carrier frequencies (HLA and KIR) and copy number variation (KIR). HLA-KIR interaction associations were also modeled. Additionally, single cell sequencing was conducted in peripheral blood mononuclear cells from 16 cases and 16 controls, NK cells were sorted and phenotyped.

Results

Anti-NMDAR encephalitis showed a weak HLA association with DRB1*01:01~DQA1*01:01~DQB1*05:01 (OR=1.57, 1.51, 1.45; respectively), and DRB1*11:01 (OR=1.60); these effects were stronger in European descendants and in patients without an underlying ovarian teratoma. More interestingly, we found increased copy number variation of KIR2DL5B (OR=1.72), principally due to an overrepresentation of KIR2DL5B*00201. Further, we identified two allele associations in framework genes, KIR2DL4*00103 (25.4% vs. 12.5% in controls, OR=1.98) and KIR3DL3*00302 (5.3% vs. 1.3%, OR=4.44). Notably, the ligands of these KIR2DL4 and KIR3DL3, respectively, HLA-G and HHLA2, are known to act as immune checkpoint with immunosuppressive functions. However, we did not find differences in specific KIR-HLA ligand interactions or HLA-G polymorphisms between cases and controls. Similarly, gene expression of CD56dim or CD56bright NK cells did not differ between cases and controls.

Discussion

Our observations for the first time suggest that the HLA-KIR axis might be involved in anti-NMDAR encephalitis. While the genetic risk conferred by the identified polymorphisms appears small, a role of this axis in the pathophysiology of this disease appears highly plausible and should be analyzed in future studies.

Inflammatory and Immune Mechanisms for Atherosclerotic Cardiovascular Disease in HIV

Atherosclerotic vascular disease disproportionately affects persons living with HIV (PLWH) compared to those without. The reasons for the excess risk include dysregulated immune response and inflammation related to HIV infection itself, comorbid conditions, and co-infections. Here, we review an updated understanding of immune and inflammatory pathways underlying atherosclerosis in PLWH, including effects of viral products, soluble mediators and chemokines, innate and adaptive immune cells, and important co-infections. We also present potential therapeutic targets which may reduce cardiovascular risk in PLWH.

Distinct developmental pathways generate functionally distinct populations of natural killer cells

Natural killer (NK) cells function by eliminating virus-infected or tumor cells. Here we identified an NK-lineage-biased progenitor population, referred to as early NK progenitors (ENKPs), which developed into NK cells independently of common precursors for innate lymphoid cells (ILCPs). ENKP-derived NK cells (ENKP_NK cells) and ILCP-derived NK cells (ILCP_NK cells) were transcriptionally different. We devised combinations of surface markers that identified highly enriched ENKP_NK and ILCP_NK cell populations in wild-type mice. Furthermore, Ly49H+ NK cells that responded to mouse cytomegalovirus infection primarily developed from ENKPs, whereas ILCP_NK cells were better IFNγ producers after infection with Salmonella and herpes simplex virus. Human CD56dim and CD56bright NK cells were transcriptionally similar to ENKP_NK cells and ILCP_NK cells, respectively. Our findings establish the existence of two pathways of NK cell development that generate functionally distinct NK cell subsets in mice and further suggest these pathways may be conserved in humans.

Transcriptomic diversity of innate lymphoid cells in human lymph nodes compared to BM and spleen

Innate lymphoid cells (ILCs) are largely tissue-resident, mostly described within the mucosal tissues. However, their presence and functions in the human draining lymph nodes (LNs) are unknown. Our study unravels the tissue-specific transcriptional profiles of 47,287 CD127+ ILCs within the human abdominal and thoracic LNs. LNs contain a higher frequency of CD127+ ILCs than in BM or spleen. We define independent stages of ILC development, including EILP and pILC in the BM. These progenitors exist in LNs in addition to naïve ILCs (nILCs) that can differentiate into mature ILCs. We define three ILC1 and four ILC3 sub-clusters in the LNs. ILC1 and ILC3 subsets have clusters with high heat shock protein-encoding genes. We identify previously unrecognized regulons, including the BACH2 family for ILC1 and the ATF family for ILC3. Our study is the comprehensive characterization of ILCs in LNs, providing an in-depth understanding of ILC-mediated immunity in humans.

Venetoclax acts as an immunometabolic modulator to potentiate adoptive NK cell immunotherapy against leukemia

Natural killer (NK) cell-based immunotherapy holds promise for cancer treatment; however, its efficacy remains limited, necessitating the development of alternative strategies. Here, we report that venetoclax, an FDA-approved BCL-2 inhibitor, directly activates NK cells, enhancing their cytotoxicity against acute myeloid leukemia (AML) both in vitro and in vivo, likely independent of BCL-2 inhibition. Through comprehensive approaches, including bulk and single-cell RNA sequencing, avidity measurement, and functional assays, we demonstrate that venetoclax increases the avidity of NK cells to AML cells and promotes lytic granule polarization during immunological synapse (IS) formation. Notably, we identify a distinct CD161lowCD218b+ NK cell subpopulation that exhibits remarkable sensitivity to venetoclax treatment. Furthermore, venetoclax promotes mitochondrial respiration and ATP synthesis via the NF-κB pathway, thereby facilitating IS formation in NK cells. Collectively, our findings establish venetoclax as a multifaceted immunometabolic modulator of NK cell function and provide a promising strategy for augmenting NK cell-based cancer immunotherapy.

Dynamic Transcriptional Programs During Single NK Cell Killing: Connecting Form to Function in Cellular Immunotherapy

Introduction

Natural killer (NK) cell-based therapies are a promising new method for treating indolent cancer, however engineering new therapies is complex and progress towards therapy for solid tumors is slow. New methods for determining the underlying intracellular signaling driving the killing phenotype would significantly improve this progress.

Methods

We combined single-cell RNA sequencing with live cell imaging of a model system of NK cell killing to correlate transcriptomic data with functional output. A model of NK cell activity, the NK-92 cell line killing of HeLa cervical cancer cells, was used for these studies. NK cell killing activity was observed by microscopy during co-culture with target HeLa cells and killing activity subsequently manually mapped based on NK cell location and Annexin V expression. NK cells from this culture system were profiled by single-cell RNA sequencing using the 10× Genomics platform, and transcription factor activity inferred using the Viper and DoRothEA R packages. Luminescent microscopy of reporter constructs in the NK cells was then used to correlate activity of inferred transcriptional activity with killing activity.

Results

NK cells had heterogeneous killing activity during 10 h of culture with target HeLa cells. Analysis of the single cell sequencing data identified Nuclear Factor Kappa B (NF-κB), Signal Transducer and Activator of Transcription 1 (STAT1) and MYC activity as potential drivers of NK cell functional phenotype in our model system. Live cell imaging of the transcription factor activity found NF-κB activity was significantly correlated with past killing activity. No correlation was observed between STAT1 or MYC activity and NK cell killing.

Conclusions

Combining luminescent microscopy of transcription factor activity with single-cell RNA sequencing is an effective means of assigning functional phenotypes to inferred transcriptomics data.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-024-00812-3.

Spatial transcriptomic brain imaging reveals the effects of immunomodulation therapy on specific regional brain cells in a mouse dementia model

Increasing evidence of brain-immune crosstalk raises expectations for the efficacy of novel immunotherapies in Alzheimer's disease (AD), but the lack of methods to examine brain tissues makes it difficult to evaluate therapeutics. Here, we investigated the changes in spatial transcriptomic signatures and brain cell types using the 10x Genomics Visium platform in immune-modulated AD models after various treatments. To proceed with an analysis suitable for barcode-based spatial transcriptomics, we first organized a workflow for segmentation of neuroanatomical regions, establishment of appropriate gene combinations, and comprehensive review of altered brain cell signatures. Ultimately, we investigated spatial transcriptomic changes following administration of immunomodulators, NK cell supplements and an anti-CD4 antibody, which ameliorated behavior impairment, and designated brain cells and regions showing probable associations with behavior changes. We provided the customized analytic pipeline into an application named STquantool. Thus, we anticipate that our approach can help researchers interpret the real action of drug candidates by simultaneously investigating the dynamics of all transcripts for the development of novel AD therapeutics.

A comprehensive single-cell breast tumor atlas defines epithelial and immune heterogeneity and interactions predicting anti-PD-1 therapy response

We present an integrated single-cell RNA sequencing atlas of the primary breast tumor microenvironment (TME) containing 236,363 cells from 119 biopsy samples across eight datasets. In this study, we leverage this resource for multiple analyses of immune and cancer epithelial cell heterogeneity. We define natural killer (NK) cell heterogeneity through six subsets in the breast TME. Because NK cell heterogeneity correlates with epithelial cell heterogeneity, we characterize epithelial cells at the level of single-gene expression, molecular subtype, and 10 categories reflecting intratumoral transcriptional heterogeneity. We develop InteractPrint, which considers how cancer epithelial cell heterogeneity influences cancer-immune interactions. We use T cell InteractPrint to predict response to immune checkpoint inhibition (ICI) in two breast cancer clinical trials testing neoadjuvant anti-PD-1 therapy. T cell InteractPrint was predictive of response in both trials versus PD-L1 (AUC = 0.82, 0.83 vs. 0.50, 0.72). This resource enables additional high-resolution investigations of the breast TME.

Persistence of a Skewed Repertoire of NK Cells in People with HIV-1 on Long-Term Antiretroviral Therapy

HIV-1 infection greatly alters the NK cell phenotypic and functional repertoire. This is highlighted by the expansion of a rare population of FcRγ- NK cells exhibiting characteristics of traditional immunologic memory in people with HIV (PWH). Although current antiretroviral therapy (ART) effectively controls HIV-1 viremia and disease progression, its impact on HIV-1-associated NK cell abnormalities remains unclear. To address this, we performed a longitudinal analysis detailing conventional and memory-like NK cell characteristics in n = 60 PWH during the first 4 y of ART. Throughout this regimen, a skewed repertoire of cytokine unresponsive FcRγ- memory-like NK cells persisted and accompanied an overall increase in NK surface expression of CD57 and KLRG1, suggestive of progression toward immune senescence. These traits were linked to elevated serum inflammatory biomarkers and increasing Ab titers to human CMV, with human CMV viremia detected in approximately one-third of PWH at years 1-4 of ART. Interestingly, 40% of PWH displayed atypical NK cell subsets, representing intermediate stages of NK-poiesis based on single-cell multiomic trajectory analysis. Our findings indicate that NK cell irregularities persist in PWH despite long-term ART, underscoring the need to better understand the causative mechanisms that prevent full restoration of immune health in PWH.

Differences in nature killer cell response and interference with mitochondrial DNA induced apoptosis in moxifloxacin environment

OBJECTIVES

As antibiotics become more prevalent, accuracy and safety are critical. Moxifloxacin (MXF) have been reported to have immunomodulatory effects on a variety of immune cells and even anti-proliferative and pro-apoptotic effects, but the mechanism of action is not fully clear.

METHODS

Peripheral blood mononuclear cells (PBMC) from experimental groups of healthy adults (n = 3) were treated with MXF (10ug/ml) in vitro for 24 h. Single-cell sequencing was performed to investigate differences in the response of each immune cell to MXF. Flow cytometry determined differential gene expression in subsets of most damaged NK cells. Pseudo-time analysis identified drivers that influence MXF-stimulated cell differentiation. Detection of mitochondrial DNA and its involvement in the mitochondrial respiratory chain pathway clarifies the origin of MXF-induced stress injury.

RESULTS

Moxifloxacin-environmental NK cells are markedly reduced: a new subset of NK cells emerges, and immediate-early-response genes in this subset indicate the presence of an early activation response. The inhibitory receptor-dominant subset shows enhanced activation, leading to increased expression of cytokines and chemokines. The near-mature subset showed greater cytotoxicity and the most pronounced cellular damage. CD56bright cells responded by antagonizing the regulation of activation and inhibitory signals, demonstrating a strong cleavage capacity. The severe depletion of mitochondrial genes was focused on apoptosis induced by the mitochondrial respiratory chain complex.

CONCLUSION

NK cells exhibit heightened sensitivity to the MXF environment. Different NK subsets upregulate the expression of cytokines and chemokines through different activation pathways. Concurrently, MXF induces impairment of the mitochondrial oxidative phosphorylation system, culminating in apoptosis.

Dynamic immune recovery process after liver transplantation revealed by single-cell multi-omics analysis

Elucidating the temporal process of immune remodeling under immunosuppressive treatment after liver transplantation (LT) is critical for precise clinical management strategies. Here, we performed a single-cell multi-omics analysis of peripheral blood mononuclear cells (PBMCs) collected from LT patients (with and without acute cellular rejection [ACR]) at 13 time points. Validation was performed in two independent cohorts with additional LT patients and healthy controls. Our study revealed a four-phase recovery process after LT and delineated changes in immune cell composition, expression programs, and interactions along this process. The intensity of the immune response differs between the ACR and non-ACR patients. Notably, the newly identified inflamed NK cells, CD14+RNASE2+ monocytes, and FOS-expressing monocytes emerged as predictive indicators of ACR. This study illuminates the longitudinal evolution of the immune cell landscape under tacrolimus-based immunosuppressive treatment during LT recovery, providing a four-phase framework that aids the clinical management of LT patients.

Single-cell, single-nucleus, and spatial transcriptomics characterization of the immunological landscape in the healthy and PSC human liver

BACKGROUND & AIMS

Primary sclerosing cholangitis (PSC) is an immune-mediated cholestatic liver disease for which there is an unmet need to understand the cellular composition of the affected liver and how it underlies disease pathogenesis. We aimed to generate a comprehensive atlas of the PSC liver using multi-omic modalities and protein-based functional validation.

METHODS

We employed single-cell and single-nucleus RNA sequencing (47,156 cells and 23,000 nuclei) and spatial transcriptomics (one sample by 10x Visium and five samples with Nanostring GeoMx DSP) to profile the cellular ecosystem in 10 PSC livers. Transcriptomic profiles were compared to 24 neurologically deceased donor livers (107,542 cells) and spatial transcriptomics controls, as well as 18,240 cells and 20,202 nuclei from three PBC livers. Flow cytometry was performed to validate PSC-specific differences in immune cell phenotype and function.

RESULTS

PSC explants with parenchymal cirrhosis and prominent periductal fibrosis contained a population of cholangiocyte-like hepatocytes that were surrounded by diverse immune cell populations. PSC-associated biliary, mesenchymal, and endothelial populations expressed chemokine and cytokine transcripts involved in immune cell recruitment. Additionally, expanded CD4+ T cells and recruited myeloid populations in the PSC liver expressed the corresponding receptors to these chemokines and cytokines, suggesting potential recruitment. Tissue-resident macrophages, by contrast, were reduced in number and exhibited a dysfunctional and downregulated inflammatory response to lipopolysaccharide and interferon-γ stimulation.

CONCLUSIONS

We present a comprehensive atlas of the PSC liver and demonstrate an exhaustion-like phenotype of myeloid cells and markers of chronic cytokine expression in late-stage PSC lesions. This atlas expands our understanding of the cellular complexity of PSC and has potential to guide the development of novel treatments.

IMPACT AND IMPLICATIONS

Primary sclerosing cholangitis (PSC) is a rare liver disease characterized by chronic inflammation and irreparable damage to the bile ducts, which eventually results in liver failure. Due to a limited understanding of the underlying pathogenesis of disease, treatment options are limited. To address this, we sequenced healthy and diseased livers to compare the activity, interactions, and localization of immune and non-immune cells. This revealed that hepatocytes lining PSC scar regions co-express cholangiocyte markers, whereas immune cells infiltrate the scar lesions. Of these cells, macrophages, which typically contribute to tissue repair, were enriched in immunoregulatory genes and demonstrated a lack of responsiveness to stimulation. These cells may be involved in maintaining hepatic inflammation and could be a target for novel therapies.

HLA-DR Expression in Natural Killer Cells Marks Distinct Functional States, Depending on Cell Differentiation Stage

HLA-DR-positive NK cells, found in both healthy individuals and patients with different inflammatory diseases, are characterized as activated cells. However, data on their capacity for IFNγ production or cytotoxic response vary between studies. Thus, more precise investigation is needed of the mechanisms related to the induction of HLA-DR expression in NK cells, their associations with NK cell differentiation stage, and functional or metabolic state. In this work, HLA-DR-expressing NK cell subsets were investigated using transcriptomic analysis, metabolic activity assays, and analysis of intercellular signaling cascades. We demonstrated that HLA-DR+CD56bright NK cells were characterized by a proliferative phenotype, while HLA-DR+CD56dim NK cells exhibited features of adaptive cells and loss of inhibitory receptors with increased expression of MHC class II trans-activator CIITA. The activated state of HLA-DR-expressing NK cells was confirmed by higher levels of ATP and mitochondrial mass observed in this subset compared to HLA-DR- cells, both ex vivo and after stimulation in culture. We showed that HLA-DR expression in NK cells in vitro can be induced both through stimulation by exogenous IL-2 and IL-21, as well as through auto-stimulation by NK-cell-produced IFNγ. At the intracellular level, HLA-DR expression depended on the activation of STAT3- and ERK1/2-mediated pathways, with subsequent activation of isoform 3 of the transcription factor CIITA. The obtained results broaden the knowledge about HLA-DR-positive NK cell appearance, diversity, and functions, which might be useful in terms of understanding the role of this subset in innate immunity and assessing their possible implications in NK cell-based therapy.

Single-Cell Informatics for Tumor Microenvironment and Immunotherapy

Cancer comprises malignant cells surrounded by the tumor microenvironment (TME), a dynamic ecosystem composed of heterogeneous cell populations that exert unique influences on tumor development. The immune community within the TME plays a substantial role in tumorigenesis and tumor evolution. The innate and adaptive immune cells "talk" to the tumor through ligand-receptor interactions and signaling molecules, forming a complex communication network to influence the cellular and molecular basis of cancer. Such intricate intratumoral immune composition and interactions foster the application of immunotherapies, which empower the immune system against cancer to elicit durable long-term responses in cancer patients. Single-cell technologies have allowed for the dissection and characterization of the TME to an unprecedented level, while recent advancements in bioinformatics tools have expanded the horizon and depth of high-dimensional single-cell data analysis. This review will unravel the intertwined networks between malignancy and immunity, explore the utilization of computational tools for a deeper understanding of tumor-immune communications, and discuss the application of these approaches to aid in diagnosis or treatment decision making in the clinical setting, as well as the current challenges faced by the researchers with their potential future improvements.

Single-cell data revealed exhaustion of characteristic NK cell subpopulations and T cell subpopulations in hepatocellular carcinoma

BACKGROUND

The treatment and prognosis of patients with advanced hepatocellular carcinoma (HCC) have been a major medical challenge. Unraveling the landscape of tumor immune infiltrating cells (TIICs) in the immune microenvironment of HCC is of great significance to probe the molecular mechanisms.

METHODS

Based on single-cell data of HCC, the cell landscape was revealed from the perspective of TIICs. Special cell subpopulations were determined by the expression levels of marker genes. Differential expression analysis was conducted. The activity of each subpopulation was determined based on the highly expressed genes. CTLA4+ T-cell subpopulations affecting the prognosis of HCC were determined based on survival analysis. A single-cell regulatory network inference and clustering analysis was also performed to determine the transcription factor regulatory networks in the CTLA4+ T cell subpopulations.

RESULTS

10 cell types were identified and NK cells and T cells showed high abundance in tumor tissues. Two NK cells subpopulations were present, FGFBP2+ NK cells, B3GNT7+ NK cells. Four T cells subpopulations were present, LAG3+ T cells, CTLA4+ T cells, RCAN3+ T cells, and HPGDS+ Th2 cells. FGFBP2+ NK cells, and CTLA4+ T cells were the exhaustive subpopulation. High CTLA4+ T cells contributed to poor prognostic outcomes and promoted tumor progression. Finally, a network of transcription factors regulated by NR3C1, STAT1, and STAT3, which were activated, was present in CTLA4+ T cells.

CONCLUSION

CTLA4+ T cell subsets in HCC exhibited functional exhaustion characteristics that probably inhibited T cell function through a transcription factor network dominated by NR3C1, STAT1, and STAT3.

Therapeutic potential of natural killer cells in neuroimmunological diseases

Natural killer (NK) cells, a major component of the innate immune system, have prominent immunoregulatory, antitumor proliferation, and antiviral activities. NK cells act as a double-edged sword with therapeutic potential in neurological autoimmunity. Emerging evidence has identified NK cells are involved in the development and progression of neuroimmunological diseases such as multiple sclerosis, neuromyelitis optica spectrum disorders, autoimmune encephalitis, Guillain-Barré Syndrome, chronic inflammatory demyelinating polyneuropathy, myasthenia gravis, and idiopathic inflammatory myopathy. However, the regulatory mechanisms and functional roles of NK cells are highly variable in different clinical states of neuroimmunological diseases and need to be further determined. In this review, we summarize the evidence for the heterogenic involvement of NK cells in the above conditions. Further, we describe cutting-edge NK-cell-based immunotherapy for neuroimmunological diseases in preclinical and clinical development and highlight challenges that must be overcome to fully realize the therapeutic potential of NK cells.

Extracellular cancer‑associated fibroblasts: A novel subgroup in the cervical cancer microenvironment that exhibits tumor‑promoting roles and prognosis biomarker functions

Tumor invasion and metastasis are the processes that primarily cause adverse outcomes in patients with cervical cancer. Cancer-associated fibroblasts (CAFs), which participate in cancer progression and metastasis, are novel targets for the treatment of tumors. The present study aimed to assess the heterogeneity of CAFs in the cervical cancer microenvironment through single-cell RNA sequencing. After collecting five cervical cancer samples and obtaining the CAF-associated gene sets, the CAFs in the cervical cancer microenvironment were divided into myofibroblastic CAFs and extracellular (ec)CAFs. The ecCAFs appeared with more robust pro-tumorigenic effects than myCAFs according to enrichment analysis. Subsequently, through combining the ecCAF hub genes and bulk gene expression data for cervical cancer obtained from The Cancer Genome Atlas and Gene Ontology databases, univariate Cox regression and least absolute shrinkage and selection operator analyses were performed to establish a CAF-associated risk signature for patients with cancer. The established risk signature demonstrated a stable and strong prognostic capability in both the training and validation cohorts. Subsequently, the association between the risk signature and clinical data was evaluated, and a nomogram to facilitate clinical application was established. The risk score was demonstrated to be associated with both the tumor immune microenvironment and the therapeutic responses. Moreover, the signature also has predictive value for the prognosis of head and neck squamous cell carcinoma, and bladder urothelial carcinoma, which were also associated with human papillomavirus infection. In conclusion, the present study assessed the heterogeneity of CAFs in the cervical cancer microenvironment, and a subgroup of CAFs that may be closely associated with tumor progression was defined. Moreover, a signature based on the hub genes of ecCAFs was shown to have biomarker functionality in terms of predicting survival rates, and therefore this CAF subgroup may become a therapeutic target for cervical cancer in the future.

Exploring NK cell receptor dynamics in paediatric leukaemias: implications for immunotherapy and prognosis

Objectives

Immunotherapies targeting natural killer (NK) cell receptors have shown promise against leukaemia. Unfortunately, cancer immunosuppressive mechanisms that alter NK cell phenotype prevent such approaches from being successful. The study utilises advanced cytometry to examine how cancer immunosuppressive pathways affect NK cell phenotypic changes in clinical samples.

Methods

In this study, we conducted a high-dimensional examination of the cell surface expression of 16 NK cell receptors in paediatric patients with acute myeloid leukaemia and acute lymphoblastic leukaemia, as well as in samples of non-age matched adult peripheral blood (APB) and umbilical cord blood (UCB). An unsupervised analysis was carried out in order to identify NK cell populations present in paediatric leukaemias.

Results

We observed that leukaemia NK cells clustered together with UCB NK cells and expressed relatively higher levels of the NKG2A receptor compared to APB NK cells. In addition, CD56dimCD16+CD57- NK cells lacking NKG2A expression were mainly absent in paediatric leukaemia patients. However, CD56br NK cell populations expressing high levels of NKG2A were highly represented in paediatric leukaemia patients. NKG2A expression on leukaemia NK cells was found to be positively correlated with the expression of its ligand, suggesting that the NKG2A-HLA-E interaction may play a role in modifying NK cell responses to leukaemia cells.

Conclusion

We provide an in-depth analysis of NK cell populations in paediatric leukaemia patients. These results support the development of immunotherapies targeting immunosuppressive receptors, such as NKG2A, to enhance innate immunity against paediatric leukaemia.

Single-cell RNA sequencing of cervical exfoliated cells reveals potential biomarkers and cellular pathogenesis in cervical carcinogenesis

Cervical cancer (CC) is a common gynecological malignancy. Despite the current screening methods have been proved effectively and significantly decreased CC morbidity and mortality, deficiencies still exist. Single-cell RNA sequencing (scRNA-seq) approach can identify the complex and rare cell populations at single-cell resolution. By scRNA-seq, the heterogeneity of tumor microenvironment across cervical carcinogenesis has been mapped and described. Whether these alterations could be detected and applied to CC screening is unclear. Herein, we performed scRNA-seq of 56,173 cervical exfoliated cells from 15 samples, including normal cervix, low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), and malignancy. The present study delineated the alteration of immune and epithelial cells derived during the cervical lesion progression. A subset of lipid-associated macrophage was identified as a tumor-promoting element and could serve as a biomarker for predicting the progression of LSIL into HSIL, which was then verified by immunofluorescence. Furthermore, cell-cell communication analysis indicated the SPP1-CD44 axis might exhibit a protumor interaction between epithelial cell and macrophage. In this study, we investigated the cervical multicellular ecosystem in cervical carcinogenesis and identified potential biomarkers for early detection.

Integrated Analysis of Single-Cell and Bulk RNA Sequencing Data Reveals Memory-like NK Cell Subset Associated with Mycobacterium tuberculosis Latency

Natural killer (NK) cells are innate-like lymphocytes that belong to the family of type-1 innate lymphoid cells and rapidly respond to virus-infected and tumor cells. In this study, we have combined scRNA-seq data and bulk RNA-seq data to define the phenotypic and molecular characteristics of peripheral blood NK cells. While the role of NK cells in immune surveillance against virus infections and tumors has been well established, their contribution to protective responses to other intracellular microorganisms, such as Mycobacterium tuberculosis (Mtb), is still poorly understood. In this study, we have combined scRNA-seq data and bulk RNA-seq data to illuminate the molecular characteristics of circulating NK cells in patients with active tuberculosis (TB) disease and subjects with latent Mtb infection (LTBI) and compared these characteristics with those of healthy donors (HDs) and patients with non-TB other pulmonary infectious diseases (ODs). We show here that the NK cell cluster was significantly increased in LTBI subjects, as compared to patients with active TB or other non-TB pulmonary diseases and HD, and this was mostly attributable to the expansion of an NK cell population expressing KLRC2, CD52, CCL5 and HLA-DRB1, which most likely corresponds to memory-like NK2.1 cells. These data were validated by flow cytometry analysis in a small cohort of samples, showing that LTBI subjects have a significant expansion of NK cells characterized by the prevalence of memory-like CD52+ NKG2C+ NK cells. Altogether, our results provide some new information on the role of NK cells in protective immune responses to Mtb.

Single-cell analysis of immune recognition in chronic myeloid leukemia patients following tyrosine kinase inhibitor discontinuation

Immunological control of residual leukemia cells is thought to occur in patients with chronic myeloid leukemia (CML) that maintain treatment-free remission (TFR) following tyrosine kinase inhibitor (TKI) discontinuation. To study this, we analyzed 55 single-cell RNA and T cell receptor (TCR) sequenced samples (scRNA+TCRαβ-seq) from patients with CML (n = 13, N = 25), other cancers (n = 28), and healthy (n = 7). The high number and active phenotype of natural killer (NK) cells in CML separated them from healthy and other cancers. Most NK cells in CML belonged to the active CD56dim cluster with high expression of GZMA/B, PRF1, CCL3/4, and IFNG, with interactions with leukemic cells via inhibitory LGALS9-TIM3 and PVR-TIGIT interactions. Accordingly, upregulation of LGALS9 was observed in CML target cells and TIM3 in NK cells when co-cultured together. Additionally, we created a classifier to identify TCRs targeting leukemia-associated antigen PR1 and quantified anti-PR1 T cells in 90 CML and 786 healthy TCRβ-sequenced samples. Anti-PR1 T cells were more prevalent in CML, enriched in bone marrow samples, and enriched in the mature, cytotoxic CD8 + TEMRA cluster, especially in a patient maintaining TFR. Our results highlight the role of NK cells and anti-PR1 T cells in anti-leukemic immune responses in CML.

Single-cell functional genomics reveals determinants of sensitivity and resistance to natural killer cells in blood cancers

Cancer cells can evade natural killer (NK) cell activity, thereby limiting anti-tumor immunity. To reveal genetic determinants of susceptibility to NK cell activity, we examined interacting NK cells and blood cancer cells using single-cell and genome-scale functional genomics screens. Interaction of NK and cancer cells induced distinct activation and type I interferon (IFN) states in both cell types depending on the cancer cell lineage and molecular phenotype, ranging from more sensitive myeloid to less sensitive B-lymphoid cancers. CRISPR screens in cancer cells uncovered genes regulating sensitivity and resistance to NK cell-mediated killing, including adhesion-related glycoproteins, protein fucosylation genes, and transcriptional regulators, in addition to confirming the importance of antigen presentation and death receptor signaling pathways. CRISPR screens with a single-cell transcriptomic readout provided insight into underlying mechanisms, including regulation of IFN-γ signaling in cancer cells and NK cell activation states. Our findings highlight the diversity of mechanisms influencing NK cell susceptibility across different cancers and provide a resource for NK cell-based therapies.

Abundant binary promoter switches in lineage-determining transcription factors indicate a digital component of cell fate determination

Previous studies of the murine Ly49 and human KIR gene clusters implicated competing sense and antisense promoters in the control of variegated gene expression. In the current study, an examination of transcription factor genes defines an abundance of convergent and divergent sense/antisense promoter pairs, suggesting that competing promoters may control cell fate determination. Differentiation of CD34+ hematopoietic progenitors in vitro shows that cells with GATA1 antisense transcription have enhanced GATA2 transcription and a mast cell phenotype, whereas cells with GATA2 antisense transcription have increased GATA1 transcripts and an erythroblast phenotype. Detailed analyses of the AHR and RORC genes demonstrate the ability of competing promoters to act as binary switches and the association of antisense transcription with an immature/progenitor cell phenotype. These data indicate that alternative cell fates generated by promoter competition in lineage-determining transcription factors contribute to the programming of cell differentiation.

Single cell transcriptomics shows that malaria promotes unique regulatory responses across multiple immune cell subsets

Plasmodium falciparum malaria drives immunoregulatory responses across multiple cell subsets, which protects from immunopathogenesis, but also hampers the development of effective anti-parasitic immunity. Understanding malaria induced tolerogenic responses in specific cell subsets may inform development of strategies to boost protective immunity during drug treatment and vaccination. Here, we analyse the immune landscape with single cell RNA sequencing during P. falciparum malaria. We identify cell type specific responses in sub-clustered major immune cell types. Malaria is associated with an increase in immunosuppressive monocytes, alongside NK and γδ T cells which up-regulate tolerogenic markers. IL-10-producing Tr1 CD4 T cells and IL-10-producing regulatory B cells are also induced. Type I interferon responses are identified across all cell types, suggesting Type I interferon signalling may be linked to induction of immunoregulatory networks during malaria. These findings provide insights into cell-specific and shared immunoregulatory changes during malaria and provide a data resource for further analysis.

Understanding cervical cancer at single-cell resolution

Cervical cancer is now the fourth most prevalent malignancy in women worldwide, representing a tremendous burden of cancer. The heterogeneity of complex tumor ecosystem impacts tumorigenesis, malignant progression, and response to treatment; thus, a thorough understanding of the tumor ecosystem is vital for enhancing the prognosis of patients with cervical cancer. The rapid development and widespread use of single-cell sequencing have generated a new paradigm of cancer research, providing a comprehensive and in-depth understanding of cancers. In this review, we give an overview of the recent advances made by leveraging single-cell sequencing studies in the dissection of cervical cancer ecosystem heterogeneity. We highlight the evolution of the cervical cancer ecosystem during tumor initiation, progression, and treatment. High-resolution dissection of cervical cancer at the single-cell level has the potential to drive the development of targeted therapies and enable the realization of personalized medicine.

MAGNETO: Cell type marker panel generator from single-cell transcriptomic data

Single-cell RNA sequencing experiments produce data useful to identify different cell types, including uncharacterized and rare ones. This enables us to study the specific functional roles of these cells in different microenvironments and contexts. After identifying a (novel) cell type of interest, it is essential to build succinct marker panels, composed of a few genes referring to cell surface proteins and clusters of differentiation molecules, able to discriminate the desired cells from the other cell populations. In this work, we propose a fully-automatic framework called MAGNETO, which can help construct optimal marker panels starting from a single-cell gene expression matrix and a cell type identity for each cell. MAGNETO builds effective marker panels solving a tailored bi-objective optimization problem, where the first objective regards the identification of the genes able to isolate a specific cell type, while the second conflicting objective concerns the minimization of the total number of genes included in the panel. Our results on three public datasets show that MAGNETO can identify marker panels that identify the cell populations of interest better than state-of-the-art approaches. Finally, by fine-tuning MAGNETO, our results demonstrate that it is possible to obtain marker panels with different specificity levels.

Deciphering the localization and trajectory of human natural killer cell development

Innate immune cells represent the first line of cellular immunity, comprised of both circulating and tissue-resident natural killer cells and innate lymphoid cells. These innate lymphocytes arise from a common CD34+ progenitor that differentiates into mature natural killer cells and innate lymphoid cells. The successive stages in natural killer cell maturation are characterized by increased lineage restriction and changes to phenotype and function. Mechanisms of human natural killer cell development have not been fully elucidated, especially the role of signals that drive the spatial localization and maturation of natural killer cells. Cytokines, extracellular matrix components, and chemokines provide maturation signals and influence the trafficking of natural killer cell progenitors to peripheral sites of differentiation. Here we present the latest advances in our understanding of natural killer and innate lymphoid cell development in peripheral sites, including secondary lymphoid tissues (i.e. tonsil). Recent work in the field has provided a model for the spatial distribution of natural killer cell and innate lymphoid cell developmental intermediates in tissue and generated further insights into the developmental niche. In support of this model, future studies using multifaceted approaches seek to fully map the developmental trajectory of human natural killer cells and innate lymphoid cells in secondary lymphoid tissues.

Emergence of human CMV-induced NKG2C+ NK cells is associated with CD8+ T-cell recovery after allogeneic HCT

Cytomegalovirus (CMV) infection is associated with the expansion of a mature NKG2C+FcεR1γ- natural killer (NK) cell population. The exact mechanism underlying the emergence of NKG2C+ NK cells, however, remains unknown. Allogeneic hematopoietic cell transplantation (HCT) provides an opportunity to longitudinally study lymphocyte recovery in the setting of CMV reactivation, particularly in patients receiving T-cell-depleted (TCD) allografts. We analyzed peripheral blood lymphocytes from 119 patients at serial time points after infusion of their TCD allograft and compared immune recovery with that in samples obtained from recipients of T-cell-replete (T-replete) (n = 96) or double umbilical cord blood (DUCB) (n = 52) allografts. NKG2C+ NK cells were detected in 92% (45 of 49) of recipients of TCD HCT who experienced CMV reactivation. Although NKG2A+ cells were routinely identifiable early after HCT, NKG2C+ NK cells were identified only after T cells could be detected. T-cell reconstitution occurred at variable times after HCT among patients and predominantly comprised CD8+ T cells. In patients with CMV reactivation, recipients of TCD HCT expressed significantly higher frequencies of NKG2C+ and CD56neg NK cells compared with patients who received T-replete HCT or DUCB transplantation. NKG2C+ NK cells after TCD HCT were CD57+FcεR1γ+ and degranulated significantly more in response to target cells compared with the adaptive the NKG2C+CD57+FcεR1γ- NK cell population. We conclude that the presence of circulating T cells is associated with the expansion of a CMV-induced NKG2C+ NK cell population, a potentially novel example of developmental cooperation between lymphocyte populations in response to viral infection.

Dual Immune Checkpoint Blockade Induces Analogous Alterations in the Dysfunctional CD8+ T-cell and Activated Treg Compartment

To dissect the effect of neoadjuvant PD-1 and CTLA4 blockade on intratumoral T cells in treatment-naive head and neck squamous cell carcinoma, we analyzed primary tumor immune infiltrates from responding and nonresponding patients. At baseline, a higher ratio between active (4-1BB/OX40+) and inactive regulatory CD4+ T cells was associated with immunotherapy response. Furthermore, upon therapy, this active regulatory T-cell (Treg) population showed a profound decrease in responding patients. In an analogous process, intratumoral dysfunctional CD8+ T cells displayed decreased expression of activity and dysfunction-related genes in responding patients, whereas in clinical nonresponders, natural killer cells showed an increased cytotoxic profile early upon treatment. These data reveal immunologic changes in response to dual PD-1/CTLA4 blockade, including a parallel remodeling of presumed tumor-reactive Treg and CD8+ T-cell compartments in responding patients, and indicate that the presence of activated Tregs at baseline may be associated with response.

SIGNIFICANCE

In head and neck squamous cell carcinoma, neoadjuvant PD-1/CTLA4 blockade has shown substantial response rates (20%-35%). As recognition of tumor antigens by T cells appears to be a critical driver of therapy response, a better understanding of alterations in T-cell state that are associated with response and resistance is of importance. This article is featured in Selected Articles from This Issue, p. 2109.

Integrative Magneto-Microfluidic Separation of Immune Cells Facilitates Clinical Functional Assays

Accurately analyzing the functional activities of natural killer (NK) cells in clinical diagnosis remains challenging due to their coupling with other immune effectors. To address this, an integrated immune cell separator is required, which necessitates a streamlined sample preparation workflow including immunological cell isolation, removal of excess red blood cells (RBCs), and buffer exchange for downstream analysis. Here, a self-powered integrated magneto-microfluidic cell separation (SMS) chip is presented, which outputs high-purity target immune cells by simply inputting whole blood. The SMS chip intensifies the magnetic field gradient using an iron sphere-filled inlet reservoir for high-performance immuno-magnetic cell selection and separates target cells size-selectively using a microfluidic lattice for RBC removal and buffer exchange. In addition, the chip incorporates self-powered microfluidic pumping through a degassed polydimethylsiloxane chip, enabling the rapid isolation of NK cells at the place of blood collection within 40 min. This chip is used to isolate NK cells from whole blood samples of hepatocellular cancer patients and healthy volunteers and examined their functional activities to identify potential abnormalities in NK cell function. The SMS chip is simple to use, rapid to sort, and requires small blood volumes, thus facilitating the use of immune cell subtypes for cell-based diagnosis.

Mapping the chemotactic landscape in NK cells reveals subset-specific synergistic migratory responses to dual chemokine receptor ligation

BACKGROUND

Natural killer (NK) cells have a unique capability of spontaneous cytotoxicity against malignant cells and hold promise for off-the-shelf cell therapy against cancer. One of the key challenges in the field is to improve NK cell homing to solid tumors.

METHODS

To gain a deeper understanding of the cellular mechanisms regulating trafficking of NK cells into the tumor, we used high-dimensional flow cytometry, mass cytometry, and single-cell RNA-sequencing combined with functional assays, creating a comprehensive map of human NK cell migration phenotypes.

FINDINGS

We found that the chemokine receptor repertoire of peripheral blood NK cells changes in a coordinated manner becoming progressively more diversified during NK cell differentiation and correlating tightly with the migratory response of the distinct NK cell subsets. Simultaneous ligation of CXCR1/2 and CX3CR1, synergistically potentiated the migratory response of NK cells. Analysis of 9471 solid cancers from publicly available TCGA/TARGET repositories revealed dominant chemokine patterns that varied across tumor types but with no tumor group expressing ligands for more than one chemokine receptor present on mature NK cells.

INTERPRETATION

The finding that chemokine stimulation can elicit a synergistic migratory response in NK cells combined with the identified lack of naturally occurring pairs of chemokines-chemokine receptors in human cancers may explain the systematic exclusion of NK cells from the tumor microenvironment and provides a basis for engineering next-generation NK cell therapies against malignancies.

FUNDING

The Polish Ministry of Science and Higher Education, the National Science Centre, Poland, The Norwegian Cancer Society, the Norwegian Research Council, the South-Eastern Norway Regional Health Authority, The Swedish Cancer Society, the Swedish Children's Cancer Foundation, The Swedish Research Council, The Center of Excellence: Precision Immunotherapy Alliance, Knut and Alice Wallenberg Foundation and National Cancer Institute.

The immunome of mobilized peripheral blood stem cells is predictive of long-term outcomes and therapy-related myeloid neoplasms in patients with multiple myeloma undergoing autologous stem cell transplant

Upfront autologous stem cell transplant (ASCT) is the standard of care for newly diagnosed multiple myeloma (MM) patients. However, relapse is ubiquitous and therapy-related myeloid neoplasms (t-MN) post-ASCT are commonly associated with poor outcomes. We hypothesized that the enrichment of abnormal myeloid progenitors and immune effector cells (IEC) in the peripheral blood stem cells (PBSCs) is associated with a higher risk of relapse and/or development of t-MN. We performed a comprehensive myeloid and lymphoid immunophenotyping on PBSCs from 54 patients with MM who underwent ASCT. Median progression-free (PFS), myeloid neoplasm-free (MNFS), and overall survival (OS) from ASCT were 49.6 months (95% CI: 39.5-Not Reached), 59.7 months (95% CI: 55-74), and 75.6 months (95% CI: 62-105), respectively. Abnormal expression of CD7 and HLA-DR on the myeloid progenitor cells was associated with an inferior PFS, MNFS, and OS. Similarly, enrichment of terminally differentiated (CD27/CD28-, CD57/KLRG1+) and exhausted (TIGIT/PD-1+) T-cells, and inhibitory NK-T like (CD159a+/CD56+) T-cells was associated with inferior PFS, MNFS, and OS post-transplant. Our observation of abnormal myeloid and IEC phenotype being present even before ASCT and maintenance therapy suggests an early predisposition to t-MN and inferior outcomes for MM, and has the potential to guide sequencing of future treatment modalities.

A pan-cancer single-cell panorama of human natural killer cells

Natural killer (NK) cells play indispensable roles in innate immune responses against tumor progression. To depict their phenotypic and functional diversities in the tumor microenvironment, we perform integrative single-cell RNA sequencing analyses on NK cells from 716 patients with cancer, covering 24 cancer types. We observed heterogeneity in NK cell composition in a tumor-type-specific manner. Notably, we have identified a group of tumor-associated NK cells that are enriched in tumors, show impaired anti-tumor functions, and are associated with unfavorable prognosis and resistance to immunotherapy. Specific myeloid cell subpopulations, in particular LAMP3+ dendritic cells, appear to mediate the regulation of NK cell anti-tumor immunity. Our study provides insights into NK-cell-based cancer immunity and highlights potential clinical utilities of NK cell subsets as therapeutic targets.

CD56bright natural killer cells preferentially kill proliferating CD4+ T cells

Human CD56br natural killer (NK) cells represent a small subset of CD56+ NK cells in circulation and are largely tissue-resident. The frequency and number of CD56br NK cells in blood has been shown to increase following administration of low-dose IL-2 (LD-IL2), a therapy aimed to specifically expand CD4+ regulatory T cells (Tregs). Given the potential clinical application of LD-IL-2 immunotherapy across several immune diseases, including the autoimmune disease type 1 diabetes, a better understanding of the functional consequences of this expansion is urgently needed. In this study, we developed an in vitro co-culture assay with activated CD4+ T cells to measure NK cell killing efficiency. We show that CD56br and CD56dim NK cells show similar efficiency at killing activated CD4+ conventional T (Tconv) and Treg cell subsets. However, in contrast to CD56dim cells, CD56br NK cells preferentially target highly proliferative cells. We hypothesize that CD56br NK cells have an immunoregulatory role through the elimination of proliferating autoreactive CD4+ Tconv cells that have escaped Treg suppression. These results have implications for the interpretation of current and future trials of LD-IL-2 by providing evidence for a new, possibly beneficial immunomodulatory mechanism of LD-IL-2-expanded CD56br NK cells.

Immune cell identifier and classifier (ImmunIC) for single cell transcriptomic readouts

Single cell RNA sequencing has a central role in immune profiling, identifying specific immune cells as disease markers and suggesting therapeutic target genes of immune cells. Immune cell-type annotation from single cell transcriptomics is in high demand for dissecting complex immune signatures from multicellular blood and organ samples. However, accurate cell type assignment from single-cell RNA sequencing data alone is complicated by a high level of gene expression heterogeneity. Many computational methods have been developed to respond to this challenge, but immune cell annotation accuracy is not highly desirable. We present ImmunIC, a simple and robust tool for immune cell identification and classification by combining marker genes with a machine learning method. With over two million immune cells and half-million non-immune cells from 66 single cell RNA sequencing studies, ImmunIC shows 98% accuracy in the identification of immune cells. ImmunIC outperforms existing immune cell classifiers, categorizing into ten immune cell types with 92% accuracy. We determine peripheral blood mononuclear cell compositions of severe COVID-19 cases and healthy controls using previously published single cell transcriptomic data, permitting the identification of immune cell-type specific differential pathways. Our publicly available tool can maximize the utility of single cell RNA profiling by functioning as a stand-alone bioinformatic cell sorter, advancing cell-type specific immune profiling for the discovery of disease-specific immune signatures and therapeutic targets.

Proteogenomic and V(D)J Analysis of Human Decidual T Cells Highlights Unique Transcriptional Programming and Clonal Distribution

Immunological tolerance toward the semiallogeneic fetus is one of many maternal adaptations required for a successful pregnancy. T cells are major players of the adaptive immune system and balance tolerance and protection at the maternal-fetal interface; however, their repertoire and subset programming are still poorly understood. Using emerging single-cell RNA sequencing technologies, we simultaneously obtained transcript, limited protein, and receptor repertoire at the single-cell level, from decidual and matched maternal peripheral human T cells. The decidua maintains a tissue-specific distribution of T cell subsets compared with the periphery. We find that decidual T cells maintain a unique transcriptome programming, characterized by restraint of inflammatory pathways by overexpression of negative regulators (DUSP, TNFAIP3, ZFP36) and expression of PD-1, CTLA-4, TIGIT, and LAG3 in some CD8 clusters. Finally, analyzing TCR clonotypes demonstrated decreased diversity in specific decidual T cell populations. Overall, our data demonstrate the power of multiomics analysis in revealing regulation of fetal-maternal immune coexistence.

Unwinding the Role of the CMG Helicase in Inborn Errors of Immunity

Inborn errors of immunity (IEI) are a collection of diseases resulting from genetic causes that impact the immune system through multiple mechanisms. Natural killer cell deficiency (NKD) is one such IEI where natural killer (NK) cells are the main immune lineage affected. Though rare, the deficiency of several genes has been described as underlying causes of NKD, including MCM4, GINS1, MCM10 , and GINS4 , all of which are involved in the eukaryotic CMG helicase. The CMG helicase is made up of C DC45 – M CM – G INS and accessory proteins including MCM10. The CMG helicase plays a critical role in DNA replication by unwinding the double helix and enabling access of polymerases to single-stranded DNA, and thus helicase proteins are active in any proliferating cell. Replication stress, DNA damage, and cell cycle arrest are among the cellular phenotypes attributed to loss of function variants in CMG helicase proteins. Despite the ubiquitous function of the CMG helicase, NK cells have an apparent susceptibility to the deficiency of helicase proteins. This review will examine the role of the CMG helicase in inborn errors of immunity through the lens of NKD and further discuss why natural killer cells can be so strongly affected by helicase deficiency.

Immunologic Characterization and T cell Receptor Repertoires of Expanded Tumor-infiltrating Lymphocytes in Patients with Renal Cell Carcinoma

The successful use of expanded tumor-infiltrating lymphocytes (TIL) in adoptive TIL therapies has been reported, but the effects of the TIL expansion, immunophenotype, function, and T cell receptor (TCR) repertoire of the infused products relative to the tumor microenvironment (TME) are not well understood. In this study, we analyzed the tumor samples (n = 58) from treatment-naïve patients with renal cell carcinoma (RCC), "pre-rapidly expanded" TILs (pre-REP TIL, n = 15) and "rapidly expanded" TILs (REP TIL, n = 25) according to a clinical-grade TIL production protocol, with single-cell RNA (scRNA)+TCRαβ-seq (TCRαβ sequencing), TCRβ-sequencing (TCRβ-seq), and flow cytometry. REP TILs encompassed a greater abundance of CD4⁺ than CD8⁺ T cells, with increased LAG-3 and low PD-1 expressions in both CD4⁺ and CD8⁺ T cell compartments compared with the pre-REP TIL and tumor T cells. The REP protocol preferentially expanded small clones of the CD4⁺ phenotype (CD4, IL7R, KLRB1) in the TME, indicating that the largest exhausted T cell clones in the tumor do not expand during the expansion protocol. In addition, by generating a catalog of RCC-associated TCR motifs from >1,000 scRNA+TCRαβ-seq and TCRβ-seq RCC, healthy and other cancer sample cohorts, we quantified the RCC-associated TCRs from the expansion protocol. Unlike the low-remaining amount of anti-viral TCRs throughout the expansion, the quantity of the RCC-associated TCRs was high in the tumors and pre-REP TILs but decreased in the REP TILs. Our results provide an in-depth understanding of the origin, phenotype, and TCR specificity of RCC TIL products, paving the way for a more rationalized production of TILs.

Significance

TILs are a heterogenous group of immune cells that recognize and attack the tumor, thus are utilized in various clinical trials. In our study, we explored the TILs in patients with kidney cancer by expanding the TILs using a clinical-grade protocol, as well as observed their characteristics and ability to recognize the tumor using in-depth experimental and computational tools.

The NK cell checkpoint NKG2A maintains expansion capacity of human NK cells

Human natural killer (NK) cells are cytotoxic effector cells that are increasingly harnessed in cancer immunotherapy. NKG2A/CD94 is an inhibitory receptor on NK cells that has established regulatory functions in the direct interaction with target cells when engaged with its ligand, the non-classical HLA class I molecule HLA-E. Here, we confirmed NKG2A as a checkpoint molecule in primary human NK cells and identified a novel role for NKG2A in maintaining NK cell expansion capacity by dampening both proliferative activity and excessive activation-induced cell death. Maintenance of NK cell expansion capacity might contribute to the preferential accumulation of human NKG2A+ NK cells after hematopoietic cell transplantation and enrichment of functionally impaired NK cells in human cancers. Functional silencing of NKG2A for cancer immunotherapy is highly attractive but will need to consider that this might also lead to a reduced survival by driving activation-induced cell death in targeted NK cells.

Single-cell RNA sequencing of human lung innate lymphoid cells in the vascular and tissue niche reveals molecular features of tissue adaptation

Innate lymphoid cells (ILCs) are sentinels of healthy organ function, yet it is unknown how ILCs adapt to distinct anatomical niches within tissues. Here, we used a unique humanized mouse model, MISTRG mice transplanted with human hematopoietic stem and progenitor cells (HSPCs), to define the gene signatures of human ILCs in the vascular versus the tissue (extravascular) compartment of the lung. Single-cell RNA sequencing in combination with intravascular cell labeling demonstrated that heterogeneous populations of human ILCs and natural killer (NK) cells occupied the vascular and tissue niches in the lung of HSPC-engrafted MISTRG mice. Moreover, we discovered that niche-specific cues shape the molecular programs of human ILCs in the distinct sub-anatomical compartments of the lung. Specifically, extravasation of ILCs into the lung tissue was associated with the upregulation of genes involved in the acquisition of tissue residency, cell positioning within the lung, sensing of tissue-derived signals, cellular stress responses, nutrient uptake, and interaction with other tissue-resident immune cells. We also defined a core tissue signature shared between human ILCs and NK cells in the extravascular space of the lung, consistent with imprinting by signals from the local microenvironment. The molecular characterization of human ILCs and NK cells in the vascular and tissue niches of the lung provides new knowledge on the mechanisms of ILC tissue adaptation and represents a resource for further studies.

Single-Cell Transcriptome Analysis Identifies Subclusters with Inflammatory Fibroblast Responses in Localized Scleroderma

Localized scleroderma (LS) is an autoimmune disease with both inflammatory and fibrotic components causing an abnormal deposition of collagen in the skin and underlying tissue, often leading to disfigurement and disability. Much of its pathophysiology is extrapolated from systemic sclerosis (SSc) since the histopathology findings in the skin are nearly identical. However, LS is critically understudied. Single-cell RNA sequencing (scRNA seq) technology provides a novel way to obtain detailed information at the individual cellular level, overcoming this barrier. Here, we analyzed the affected skin of 14 patients with LS (pediatric and adult) and 14 healthy controls. Fibroblast populations were the focus, since they are the main drivers of fibrosis in SSc. We identified 12 fibroblast subclusters in LS, which overall had an inflammatory gene expression (IFN and HLA-associated genes). A myofibroblast-like cluster (SFRP4/PRSS23) was more prevalent in LS subjects and shared many upregulated genes expressed in SSc-associated myofibroblasts, though it also had strong expression of CXCL9/10/11, known CXCR3 ligands. A CXCL2/IRF1 cluster identified was unique to LS, with a robust inflammatory gene signature, including IL-6, and according to cell communication analysis are influenced by macrophages. In summary, potential disease-propagating fibroblasts and associated gene signatures were identified in LS skin via scRNA seq.

Systems biology approaches to unravel lymphocyte subsets and function

Single-cell technologies have revealed the extensive heterogeneity and complexity of the immune system. Systems biology approaches in immunology have taken advantage of the high-parameter, high-throughput data and analyzed immune cell types in a 'bottom-up' data-driven method. This approach has discovered previously unrecognized cell types and functions. Especially for human immunology, in which experimental manipulations are challenging, systems approach has become a successful means to investigate physiologically relevant contexts. This review focuses on the recent findings in lymphocyte biology, from their development, differentiation into subsets, and heterogeneity in their functions, enabled by these systems approaches. Furthermore, we review examples of the application of findings from systems approach studies and discuss how now to leave the rich dataset in the curse of high dimensionality.

Natural killer cells in aging and age-related diseases

Aging leads to escalated systemic inflammation. As the sentinel of immune system, natural killer (NK) cells are early responders that sense cues and signals from target organs and swiftly orchestrate local inflammation upon their arrival. Emerging evidence indicates a profound role of NK cells in the initiation and evolution of neuroinflammation in aging and age-related diseases. Here we discuss recent advances in NK cell biology and the organ-specific features of NK cells in normal brain aging, Alzheimer's disease, Parkinson's disease and stroke. Our increasing understanding of NK cells and their unique features in aging and age-related diseases may facilitate the future design of immune therapies targeting NK cells to benefit the old population.

TIGIT Expression on Activated NK Cells Correlates with Greater Anti-Tumor Activity but Promotes Functional Decline upon Lung Cancer Exposure: Implications for Adoptive Cell Therapy and TIGIT-Targeted Therapies

Treatments targeting TIGIT have gained a lot of attention due to strong preclinical and early clinical results, particularly with anti-PD-(L)1 therapeutics. However, this combination has failed to meet progression-free survival endpoints in phase III trials. Most of our understanding of TIGIT comes from studies of T cell function. Yet, this inhibitory receptor is often upregulated to the same, or higher, extent on NK cells in cancers. Studies in murine models have demonstrated that TIGIT inhibits NK cells and promotes exhaustion, with its effects on tumor control also being dependent on NK cells. However, there are limited studies assessing the role of TIGIT on the function of human NK cells (hNK), particularly in lung cancer. Most studies used NK cell lines or tested TIGIT blockade to reactivate exhausted cells obtained from cancer patients. For therapeutic advancement, a better understanding of TIGIT in the context of activated hNK cells is crucial, which is different than exhausted NK cells, and critical in the context of adoptive NK cell therapeutics that may be combined with TIGIT blockade. In this study, the effect of TIGIT blockade on the anti-tumor activities of human ex vivo-expanded NK cells was evaluated in vitro in the context of lung cancer. TIGIT expression was higher on activated and/or expanded NK cells compared to resting NK cells. More TIGIT⁺ NK cells expressed major activating receptors and exerted anti-tumor response as compared to TIGIT^- cells, indicating that NK cells with greater anti-tumor function express more TIGIT. However, long-term TIGIT engagement upon exposure to PVR⁺ tumors downregulated the cytotoxic function of expanded NK cells while the inclusion of TIGIT blockade increased cytotoxicity, restored the effector functions against PVR-positive targets, and upregulated immune inflammation-related gene sets. These combined results indicate that TIGIT blockade can preserve the activation state of NK cells during exposure to PVR⁺ tumors. These results support the notion that a functional NK cell compartment is critical for anti-tumor response and anti-TIGIT/adoptive NK cell combinations have the potential to improve outcomes.