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.

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.

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.

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.

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.

Disease-Associated Microglia: A Universal Immune Sensor of Neurodegeneration.

A major challenge in the field of neurodegenerative diseases and brain aging is to identify the body's intrinsic mechanism that could sense the central nervous system (CNS) damage early and protect the brain from neurodegeneration. Accumulating evidence suggests that disease-associated microglia (DAM), a recently identified subset of CNS resident macrophages found at sites of neurodegeneration, might play such a protective role. Here, we propose that microglia are endowed with a dedicated sensory mechanism, which includes the Trem2 signaling pathway, to detect damage within the CNS in the form of neurodegeneration-associated molecular patterns (NAMPs). Combining data from transcriptional analysis of DAM at single-cell level and from human genome-wide association studies (GWASs), we discuss potential function of different DAM pathways in the diseased brain and outline how manipulating DAM may create new therapeutic opportunities.

A Unique Microglia Type Associated with Restricting Development of Alzheimer's Disease.

Alzheimer's disease (AD) is a detrimental neurodegenerative disease with no effective treatments. Due to cellular heterogeneity, defining the roles of immune cell subsets in AD onset and progression has been challenging. Using transcriptional single-cell sorting, we comprehensively map all immune populations in wild-type and AD-transgenic (Tg-AD) mouse brains. We describe a novel microglia type associated with neurodegenerative diseases (DAM) and identify markers, spatial localization, and pathways associated with these cells. Immunohistochemical staining of mice and human brain slices shows DAM with intracellular/phagocytic Aß particles. Single-cell analysis of DAM in Tg-AD and triggering receptor expressed on myeloid cells 2 (Trem2)-/- Tg-AD reveals that the DAM program is activated in a two-step process. Activation is initiated in a Trem2-independent manner that involves downregulation of microglia checkpoints, followed by activation of a Trem2-dependent program. This unique microglia-type has the potential to restrict neurodegeneration, which may have important implications for future treatment of AD and other neurodegenerative diseases. VIDEO ABSTRACT.

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.

Dissecting Immune Circuits by Linking CRISPR-Pooled Screens with Single-Cell RNA-Seq.

In multicellular organisms, dedicated regulatory circuits control cell type diversity and responses. The crosstalk and redundancies within these circuits and substantial cellular heterogeneity pose a major research challenge. Here, we present CRISP-seq, an integrated method for massively parallel single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-pooled screens. We show that profiling the genomic perturbation and transcriptome in the same cell enables us to simultaneously elucidate the function of multiple factors and their interactions. We applied CRISP-seq to probe regulatory circuits of innate immunity. By sampling tens of thousands of perturbed cells in vitro and in mice, we identified interactions and redundancies between developmental and signaling-dependent factors. These include opposing effects of Cebpb and Irf8 in regulating the monocyte/macrophage versus dendritic cell lineages and differential functions for Rela and Stat1/2 in monocyte versus dendritic cell responses to pathogens. This study establishes CRISP-seq as a broadly applicable, comprehensive, and unbiased approach for elucidating mammalian regulatory circuits.

The Spectrum and Regulatory Landscape of Intestinal Innate Lymphoid Cells Are Shaped by the Microbiome.

Innate lymphoid cells (ILCs) are critical modulators of mucosal immunity, inflammation, and tissue homeostasis, but their full spectrum of cellular states and regulatory landscapes remains elusive. Here, we combine genome-wide RNA-seq, ChIP-seq, and ATAC-seq to compare the transcriptional and epigenetic identity of small intestinal ILCs, identifying thousands of distinct gene profiles and regulatory elements. Single-cell RNA-seq and flow and mass cytometry analyses reveal compartmentalization of cytokine expression and metabolic activity within the three classical ILC subtypes and highlight transcriptional states beyond the current canonical classification. In addition, using antibiotic intervention and germ-free mice, we characterize the effect of the microbiome on the ILC regulatory landscape and determine the response of ILCs to microbial colonization at the single-cell level. Together, our work characterizes the spectrum of transcriptional identities of small intestinal ILCs and describes how ILCs differentially integrate signals from the microbial microenvironment to generate phenotypic and functional plasticity.

Mapping Nucleosome Resolution Chromosome Folding in Yeast by Micro-C.

We describe a Hi-C-based method, Micro-C, in which micrococcal nuclease is used instead of restriction enzymes to fragment chromatin, enabling nucleosome resolution chromosome folding maps. Analysis of Micro-C maps for budding yeast reveals abundant self-associating domains similar to those reported in other species, but not previously observed in yeast. These structures, far shorter than topologically associating domains in mammals, typically encompass one to five genes in yeast. Strong boundaries between self-associating domains occur at promoters of highly transcribed genes and regions of rapid histone turnover that are typically bound by the RSC chromatin-remodeling complex. Investigation of chromosome folding in mutants confirms roles for RSC, "gene looping" factor Ssu72, Mediator, H3K56 acetyltransferase Rtt109, and the N-terminal tail of H4 in folding of the yeast genome. This approach provides detailed structural maps of a eukaryotic genome, and our findings provide insights into the machinery underlying chromosome compaction.

Transcriptional Heterogeneity and Lineage Commitment in Myeloid Progenitors.

Within the bone marrow, stem cells differentiate and give rise to diverse blood cell types and functions. Currently, hematopoietic progenitors are defined using surface markers combined with functional assays that are not directly linked with in vivo differentiation potential or gene regulatory mechanisms. Here, we comprehensively map myeloid progenitor subpopulations by transcriptional sorting of single cells from the bone marrow. We describe multiple progenitor subgroups, showing unexpected transcriptional priming toward seven differentiation fates but no progenitors with a mixed state. Transcriptional differentiation is correlated with combinations of known and previously undefined transcription factors, suggesting that the process is tightly regulated. Histone maps and knockout assays are consistent with early transcriptional priming, while traditional transplantation experiments suggest that in vivo priming may still allow for plasticity given strong perturbations. These data establish a reference model and general framework for studying hematopoiesis at single-cell resolution.

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.

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.

High-resolution chromatin dynamics during a yeast stress response.

Covalent histone modifications are highly conserved and play multiple roles in eukaryotic transcription regulation. Here, we mapped 26 histone modifications genome-wide in exponentially growing yeast and during a dramatic transcriptional reprogramming-the response to diamide stress. We extend prior studies showing that steady-state histone modification patterns reflect genomic processes, especially transcription, and display limited combinatorial complexity. Interestingly, during the stress response we document a modest increase in the combinatorial complexity of histone modification space, resulting from roughly 3% of all nucleosomes transiently populating rare histone modification states. Most of these rare histone states result from differences in the kinetics of histone modification that transiently uncouple highly correlated marks, with slow histone methylation changes often lagging behind the more rapid acetylation changes. Explicit analysis of modification dynamics uncovers ordered sequences of events in gene activation and repression. Together, our results provide a comprehensive view of chromatin dynamics during a massive transcriptional upheaval.

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.

A high-throughput chromatin immunoprecipitation approach reveals principles of dynamic gene regulation in mammals.

Understanding the principles governing mammalian gene regulation has been hampered by the difficulty in measuring in vivo binding dynamics of large numbers of transcription factors (TF) to DNA. Here, we develop a high-throughput Chromatin ImmunoPrecipitation (HT-ChIP) method to systematically map protein-DNA interactions. HT-ChIP was applied to define the dynamics of DNA binding by 25 TFs and 4 chromatin marks at 4 time-points following pathogen stimulus of dendritic cells. Analyzing over 180,000 TF-DNA interactions we find that TFs vary substantially in their temporal binding landscapes. This data suggests a model for transcription regulation whereby TF networks are hierarchically organized into cell differentiation factors, factors that bind targets prior to stimulus to prime them for induction, and factors that regulate specific gene programs. Overlaying HT-ChIP data on gene-expression dynamics shows that many TF-DNA interactions are established prior to the stimuli, predominantly at immediate-early genes, and identified specific TF ensembles that coordinately regulate gene-induction.