Time-resolved single-cell transcriptomics defines immune trajectories in glioblastoma.

Deciphering the cell-state transitions underlying immune adaptation across time is fundamental for advancing biology. Empirical in vivo genomic technologies that capture cellular dynamics are currently lacking. We present Zman-seq, a single-cell technology recording transcriptomic dynamics across time by introducing time stamps into circulating immune cells, tracking them in tissues for days. Applying Zman-seq resolved cell-state and molecular trajectories of the dysfunctional immune microenvironment in glioblastoma. Within 24 hours of tumor infiltration, cytotoxic natural killer cells transitioned to a dysfunctional program regulated by TGFB1 signaling. Infiltrating monocytes differentiated into immunosuppressive macrophages, characterized by the upregulation of suppressive myeloid checkpoints Trem2, Il18bp, and Arg1, over 36 to 48 hours. Treatment with an antagonistic anti-TREM2 antibody reshaped the tumor microenvironment by redirecting the monocyte trajectory toward pro-inflammatory macrophages. Zman-seq is a broadly applicable technology, enabling empirical measurements of differentiation trajectories, which can enhance the development of more efficacious immunotherapies.

Targeted delivery of tumor necrosis factor in combination with CCNU induces a T cell-dependent regression of glioblastoma.

Glioblastoma is the most aggressive primary brain tumor with an unmet need for more effective therapies. Here, we investigated combination therapies based on L19TNF, an antibody-cytokine fusion protein based on tumor necrosis factor that selectively localizes to cancer neovasculature. Using immunocompetent orthotopic glioma mouse models, we identified strong anti-glioma activity of L19TNF in combination with the alkylating agent CCNU, which cured the majority of tumor-bearing mice, whereas monotherapies only had limited efficacy. In situ and ex vivo immunophenotypic and molecular profiling in the mouse models revealed that L19TNF and CCNU induced tumor DNA damage and treatment-associated tumor necrosis. In addition, this combination also up-regulated tumor endothelial cell adhesion molecules, promoted the infiltration of immune cells into the tumor, induced immunostimulatory pathways, and decreased immunosuppression pathways. MHC immunopeptidomics demonstrated that L19TNF and CCNU increased antigen presentation on MHC class I molecules. The antitumor activity was T cell dependent and completely abrogated in immunodeficient mouse models. On the basis of these encouraging results, we translated this treatment combination to patients with glioblastoma. The clinical translation is ongoing but already shows objective responses in three of five patients in the first recurrent glioblastoma patient cohort treated with L19TNF in combination with CCNU (NCT04573192).

Dual ontogeny of disease-associated microglia and disease inflammatory macrophages in aging and neurodegeneration.

Brain macrophage populations include parenchymal microglia, border-associated macrophages, and recruited monocyte-derived cells; together, they control brain development and homeostasis but are also implicated in aging pathogenesis and neurodegeneration. The phenotypes, localization, and functions of each population in different contexts have yet to be resolved. We generated a murine brain myeloid scRNA-seq integration to systematically delineate brain macrophage populations. We show that the previously identified disease-associated microglia (DAM) population detected in murine Alzheimer's disease models actually comprises two ontogenetically and functionally distinct cell lineages: embryonically derived triggering receptor expressed on myeloid cells 2 (TREM2)-dependent DAM expressing a neuroprotective signature and monocyte-derived TREM2-expressing disease inflammatory macrophages (DIMs) accumulating in the brain during aging. These two distinct populations appear to also be conserved in the human brain. Herein, we generate an ontogeny-resolved model of brain myeloid cell heterogeneity in development, homeostasis, and disease and identify cellular targets for the treatment of neurodegeneration.

LGR5 expressing skin fibroblasts define a major cellular hub perturbed in scleroderma.

Systemic sclerosis (scleroderma, SSc) is an incurable autoimmune disease with high morbidity and mortality rates. Here, we conducted a population-scale single-cell genomic analysis of skin and blood samples of 56 healthy controls and 97 SSc patients at different stages of the disease. We found immune compartment dysfunction only in a specific subtype of diffuse SSc patients but global dysregulation of the stromal compartment, particularly in a previously undefined subset of LGR5+-scleroderma-associated fibroblasts (ScAFs). ScAFs are perturbed morphologically and molecularly in SSc patients. Single-cell multiome profiling of stromal cells revealed ScAF-specific markers, pathways, regulatory elements, and transcription factors underlining disease development. Systematic analysis of these molecular features with clinical metadata associates specific ScAF targets with disease pathogenesis and SSc clinical traits. Our high-resolution atlas of the sclerodermatous skin spectrum will enable a paradigm shift in the understanding of SSc disease and facilitate the development of biomarkers and therapeutic strategies.

Bispecific antibodies increase the therapeutic window of CD40 agonists through selective dendritic cell targeting.

Therapeutic use of agonistic anti-CD40 antibodies is a potentially powerful approach for activation of the immune response to eradicate tumors. However, the translation of this approach to clinical practice has been substantially restricted due to the severe dose-limiting toxicities observed in multiple clinical trials. Here, we demonstrate that conventional type 1 dendritic cells are essential for triggering antitumor immunity but not the toxicity of CD40 agonists, while macrophages, platelets and monocytes lead to toxic events. Therefore, we designed bispecific antibodies that target CD40 activation preferentially to dendritic cells, by coupling the CD40 agonist arm with CD11c-, DEC-205- or CLEC9A-targeting arms. These bispecific reagents demonstrate a superior safety profile compared to their parental CD40 monospecific antibody while triggering potent antitumor activity. We suggest such cell-selective bispecific agonistic antibodies as a drug platform to bypass the dose-limiting toxicities of anti-CD40, and of additional types of agonistic antibodies used for cancer immunotherapy.

Alzheimer's disease modification mediated by bone marrow-derived macrophages via a TREM2-independent pathway in mouse model of amyloidosis.

Microglia and monocyte-derived macrophages (MDM) are key players in dealing with Alzheimer's disease. In amyloidosis mouse models, activation of microglia was found to be TREM2 dependent. Here, using Trem2-/-5xFAD mice, we assessed whether MDM act via a TREM2-dependent pathway. We adopted a treatment protocol targeting the programmed cell death ligand-1 (PD-L1) immune checkpoint, previously shown to modify Alzheimer's disease via MDM involvement. Blockade of PD-L1 in Trem2-/-5xFAD mice resulted in cognitive improvement and reduced levels of water-soluble amyloid beta1-42 with no effect on amyloid plaque burden. Single-cell RNA sequencing revealed that MDM, derived from both Trem2-/- and Trem2+/+5xFAD mouse brains, express a unique set of genes encoding scavenger receptors (for example, Mrc1, Msr1). Blockade of monocyte trafficking using anti-CCR2 antibody completely abrogated the cognitive improvement induced by anti-PD-L1 treatment in Trem2-/-5xFAD mice and similarly, but to a lesser extent, in Trem2+/+5xFAD mice. These results highlight a TREM2-independent, disease-modifying activity of MDM in an amyloidosis mouse model.

Tumor cells in light-chain amyloidosis and myeloma show distinct transcriptional rewiring of normal plasma cell development.

Although light-chain amyloidosis (AL) and multiple myeloma (MM) are characterized by tumor plasma cell (PC) expansion in bone marrow (BM), their clinical presentation differs. Previous attempts to identify unique pathogenic mechanisms behind such differences were unsuccessful, and no studies have investigated the differentiation stage of tumor PCs in patients with AL and MM. We sought to define a transcriptional atlas of normal PC development in secondary lymphoid organs (SLOs), peripheral blood (PB), and BM for comparison with the transcriptional programs (TPs) of tumor PCs in AL, MM, and monoclonal gammopathy of undetermined significance (MGUS). Based on bulk and single-cell RNA sequencing, we observed 13 TPs during transition of normal PCs throughout SLOs, PB, and BM. We further noted the following: CD39 outperforms CD19 to discriminate newborn from long-lived BM-PCs; tumor PCs expressed the most advantageous TPs of normal PC differentiation; AL shares greater similarity to SLO-PCs whereas MM is transcriptionally closer to PB-PCs and newborn BM-PCs; patients with AL and MM enriched in immature TPs had inferior survival; and protein N-linked glycosylation-related TPs are upregulated in AL. Collectively, we provide a novel resource to understand normal PC development and the transcriptional reorganization of AL and other monoclonal gammopathies.

XCR1+ type 1 conventional dendritic cells drive liver pathology in non-alcoholic steatohepatitis.

Non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH) are prevalent liver conditions that underlie the development of life-threatening cirrhosis, liver failure and liver cancer. Chronic necro-inflammation is a critical factor in development of NASH, yet the cellular and molecular mechanisms of immune dysregulation in this disease are poorly understood. Here, using single-cell transcriptomic analysis, we comprehensively profiled the immune composition of the mouse liver during NASH. We identified a significant pathology-associated increase in hepatic conventional dendritic cells (cDCs) and further defined their source as NASH-induced boost in cycling of cDC progenitors in the bone marrow. Analysis of blood and liver from patients on the NAFLD/NASH spectrum showed that type 1 cDCs (cDC1) were more abundant and activated in disease. Sequencing of physically interacting cDC-T cell pairs from liver-draining lymph nodes revealed that cDCs in NASH promote inflammatory T cell reprogramming, previously associated with NASH worsening. Finally, depletion of cDC1 in XCR1DTA mice or using anti-XCL1-blocking antibody attenuated liver pathology in NASH mouse models. Overall, our study provides a comprehensive characterization of cDC biology in NASH and identifies XCR1+ cDC1 as an important driver of liver pathology.

NASH limits anti-tumour surveillance in immunotherapy-treated HCC.

Hepatocellular carcinoma (HCC) can have viral or non-viral causes1-5. Non-alcoholic steatohepatitis (NASH) is an important driver of HCC. Immunotherapy has been approved for treating HCC, but biomarker-based stratification of patients for optimal response to therapy is an unmet need6,7. Here we report the progressive accumulation of exhausted, unconventionally activated CD8+PD1+ T cells in NASH-affected livers. In preclinical models of NASH-induced HCC, therapeutic immunotherapy targeted at programmed death-1 (PD1) expanded activated CD8+PD1+ T cells within tumours but did not lead to tumour regression, which indicates that tumour immune surveillance was impaired. When given prophylactically, anti-PD1 treatment led to an increase in the incidence of NASH-HCC and in the number and size of tumour nodules, which correlated with increased hepatic CD8+PD1+CXCR6+, TOX+, and TNF+ T cells. The increase in HCC triggered by anti-PD1 treatment was prevented by depletion of CD8+ T cells or TNF neutralization, suggesting that CD8+ T cells help to induce NASH-HCC, rather than invigorating or executing immune surveillance. We found similar phenotypic and functional profiles in hepatic CD8+PD1+ T cells from humans with NAFLD or NASH. A meta-analysis of three randomized phase III clinical trials that tested inhibitors of PDL1 (programmed death-ligand 1) or PD1 in more than 1,600 patients with advanced HCC revealed that immune therapy did not improve survival in patients with non-viral HCC. In two additional cohorts, patients with NASH-driven HCC who received anti-PD1 or anti-PDL1 treatment showed reduced overall survival compared to patients with other aetiologies. Collectively, these data show that non-viral HCC, and particularly NASH-HCC, might be less responsive to immunotherapy, probably owing to NASH-related aberrant T cell activation causing tissue damage that leads to impaired immune surveillance. Our data provide a rationale for stratification of patients with HCC according to underlying aetiology in studies of immunotherapy as a primary or adjuvant treatment.

Identification of resistance pathways and therapeutic targets in relapsed multiple myeloma patients through single-cell sequencing.

Multiple myeloma (MM) is a neoplastic plasma-cell disorder characterized by clonal proliferation of malignant plasma cells. Despite extensive research, disease heterogeneity within and between treatment-resistant patients is poorly characterized. In the present study, we conduct a prospective, multicenter, single-arm clinical trial (NCT04065789), combined with longitudinal single-cell RNA-sequencing (scRNA-seq) to study the molecular dynamics of MM resistance mechanisms. Newly diagnosed MM patients (41), who either failed to respond or experienced early relapse after a bortezomib-containing induction regimen, were enrolled to evaluate the safety and efficacy of a daratumumab, carfilzomib, lenalidomide and dexamethasone combination. The primary clinical endpoint was safety and tolerability. Secondary endpoints included overall response rate, progression-free survival and overall survival. Treatment was safe and well tolerated; deep and durable responses were achieved. In prespecified exploratory analyses, comparison of 41 primary refractory and early relapsed patients, with 11 healthy subjects and 15 newly diagnosed MM patients, revealed new MM molecular pathways of resistance, including hypoxia tolerance, protein folding and mitochondria respiration, which generalized to larger clinical cohorts (CoMMpass). We found peptidylprolyl isomerase A (PPIA), a central enzyme in the protein-folding response pathway, as a potential new target for resistant MM. CRISPR-Cas9 deletion of PPIA or inhibition of PPIA with a small molecule inhibitor (ciclosporin) significantly sensitizes MM tumor cells to proteasome inhibitors. Together, our study defines a roadmap for integrating scRNA-seq in clinical trials, identifies a signature of highly resistant MM patients and discovers PPIA as a potent therapeutic target for these tumors.

Coupled scRNA-Seq and Intracellular Protein Activity Reveal an Immunosuppressive Role of TREM2 in Cancer.

Cell function and activity are regulated through integration of signaling, epigenetic, transcriptional, and metabolic pathways. Here, we introduce INs-seq, an integrated technology for massively parallel recording of single-cell RNA sequencing (scRNA-seq) and intracellular protein activity. We demonstrate the broad utility of INs-seq for discovering new immune subsets by profiling different intracellular signatures of immune signaling, transcription factor combinations, and metabolic activity. Comprehensive mapping of Arginase 1-expressing cells within tumor models, a metabolic immune signature of suppressive activity, discovers novel Arg1+ Trem2+ regulatory myeloid (Mreg) cells and identifies markers, metabolic activity, and pathways associated with these cells. Genetic ablation of Trem2 in mice inhibits accumulation of intra-tumoral Mreg cells, leading to a marked decrease in dysfunctional CD8+ T cells and reduced tumor growth. This study establishes INs-seq as a broadly applicable technology for elucidating integrated transcriptional and intra-cellular maps and identifies the molecular signature of myeloid suppressive cells in tumors.

The Physiology, Pathology, and Potential Therapeutic Applications of the TREM2 Signaling Pathway.

Alzheimer's disease, obesity-related metabolic syndrome, and cancer are the leading causes of death and among the most costly medical conditions in the Western world. In all three cases, recent discoveries establish the TREM2 receptor as a major pathology-induced immune signaling hub that senses tissue damage and activates robust immune remodeling in response to it. In this review, we summarize and question what is known and remains to be discovered about TREM2 signaling pathway, track the consequences of its activation in physiological niches and pathological contexts, and highlight the promising potential of therapeutic manipulation of TREM2 signaling.

Cross-Species Single-Cell Analysis Reveals Divergence of the Primate Microglia Program.

Microglia, the brain-resident immune cells, are critically involved in many physiological and pathological brain processes, including neurodegeneration. Here we characterize microglia morphology and transcriptional programs across ten species spanning more than 450 million years of evolution. We find that microglia express a conserved core gene program of orthologous genes from rodents to humans, including ligands and receptors associated with interactions between glia and neurons. In most species, microglia show a single dominant transcriptional state, whereas human microglia display significant heterogeneity. In addition, we observed notable differences in several gene modules of rodents compared with primate microglia, including complement, phagocytic, and susceptibility genes to neurodegeneration, such as Alzheimer's and Parkinson's disease. Our study provides an essential resource of conserved and divergent microglia pathways across evolution, with important implications for future development of microglia-based therapies in humans.

Lipid-Associated Macrophages Control Metabolic Homeostasis in a Trem2-Dependent Manner.

Immune cells residing in white adipose tissue have been highlighted as important factors contributing to the pathogenesis of metabolic diseases, but the molecular regulators that drive adipose tissue immune cell remodeling during obesity remain largely unknown. Using index and transcriptional single-cell sorting, we comprehensively map all adipose tissue immune populations in both mice and humans during obesity. We describe a novel and conserved Trem2+ lipid-associated macrophage (LAM) subset and identify markers, spatial localization, origin, and functional pathways associated with these cells. Genetic ablation of Trem2 in mice globally inhibits the downstream molecular LAM program, leading to adipocyte hypertrophy as well as systemic hypercholesterolemia, body fat accumulation, and glucose intolerance. These findings identify Trem2 signaling as a major pathway by which macrophages respond to loss of tissue-level lipid homeostasis, highlighting Trem2 as a key sensor of metabolic pathologies across multiple tissues and a potential therapeutic target in metabolic diseases.

A niche-dependent myeloid transcriptome signature defines dormant myeloma cells.

The era of targeted therapies has seen significant improvements in depth of response, progression-free survival, and overall survival for patients with multiple myeloma. Despite these improvements in clinical outcome, patients inevitably relapse and require further treatment. Drug-resistant dormant myeloma cells that reside in specific niches within the skeleton are considered a basis of disease relapse but remain elusive and difficult to study. Here, we developed a method to sequence the transcriptome of individual dormant myeloma cells from the bones of tumor-bearing mice. Our analyses show that dormant myeloma cells express a distinct transcriptome signature enriched for immune genes and, unexpectedly, genes associated with myeloid cell differentiation. These genes were switched on by coculture with osteoblastic cells. Targeting AXL, a gene highly expressed by dormant cells, using small-molecule inhibitors released cells from dormancy and promoted their proliferation. Analysis of the expression of AXL and coregulated genes in human cohorts showed that healthy human controls and patients with monoclonal gammopathy of uncertain significance expressed higher levels of the dormancy signature genes than patients with multiple myeloma. Furthermore, in patients with multiple myeloma, the expression of this myeloid transcriptome signature translated into a twofold increase in overall survival, indicating that this dormancy signature may be a marker of disease progression. Thus, engagement of myeloma cells with the osteoblastic niche induces expression of a suite of myeloid genes that predicts disease progression and that comprises potential drug targets to eradicate dormant myeloma cells.

PD-1/PD-L1 checkpoint blockade harnesses monocyte-derived macrophages to combat cognitive impairment in a tauopathy mouse model.

Alzheimer's disease (AD) is a heterogeneous disorder with multiple etiologies. Harnessing the immune system by blocking the programmed cell death receptor (PD)-1 pathway in an amyloid beta mouse model was shown to evoke a sequence of immune responses that lead to disease modification. Here, blocking PD-L1, a PD-1 ligand, was found to have similar efficacy to that of PD-1 blocking in disease modification, in both animal models of AD and of tauopathy. Targeting PD-L1 in a tau-driven disease model resulted in increased immunomodulatory monocyte-derived macrophages within the brain parenchyma. Single cell RNA-seq revealed that the homing macrophages expressed unique scavenger molecules including macrophage scavenger receptor 1 (MSR1), which was shown here to be required for the effect of PD-L1 blockade in disease modification. Overall, our results demonstrate that immune checkpoint blockade targeting the PD-1/PD-L1 pathway leads to modification of common factors that go awry in AD and dementia, and thus can potentially provide an immunotherapy to help combat these diseases.

Single cell dissection of plasma cell heterogeneity in symptomatic and asymptomatic myeloma.

Multiple myeloma, a plasma cell malignancy, is the second most common blood cancer. Despite extensive research, disease heterogeneity is poorly characterized, hampering efforts for early diagnosis and improved treatments. Here, we apply single cell RNA sequencing to study the heterogeneity of 40 individuals along the multiple myeloma progression spectrum, including 11 healthy controls, demonstrating high interindividual variability that can be explained by expression of known multiple myeloma drivers and additional putative factors. We identify extensive subclonal structures for 10 of 29 individuals with multiple myeloma. In asymptomatic individuals with early disease and in those with minimal residual disease post-treatment, we detect rare tumor plasma cells with molecular characteristics similar to those of active myeloma, with possible implications for personalized therapies. Single cell analysis of rare circulating tumor cells allows for accurate liquid biopsy and detection of malignant plasma cells, which reflect bone marrow disease. Our work establishes single cell RNA sequencing for dissecting blood malignancies and devising detailed molecular characterization of tumor cells in symptomatic and asymptomatic patients.

Single-cell characterization of haematopoietic progenitors and their trajectories in homeostasis and perturbed haematopoiesis.

The dynamics of haematopoietic stem cell differentiation and the hierarchy of oligopotent stem cells in the bone marrow remain controversial. Here we dissect haematopoietic progenitor populations at single cell resolution, deriving an unbiased reference model of transcriptional states in normal and perturbed murine bone marrow. We define the signature of the naive haematopoietic stem cell and find a continuum of core progenitor states. Core cell populations mix transcription of pre-myeloid and pre-lymphoid programs, but do not mix erythroid or megakaryocyte programs with other fates. CRISP-seq perturbation analysis confirms our models and reveals that Cebpa regulates entry into all myeloid fates, while Irf8 and PU.1 deficiency block later differentiation towards monocyte or granulocyte fates. Our transcriptional map defines a reference network model for blood progenitors and their differentiation trajectories during normal and perturbed haematopoiesis.