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.

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.

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.

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.

Microbiota Diurnal Rhythmicity Programs Host Transcriptome Oscillations.

The intestinal microbiota undergoes diurnal compositional and functional oscillations that affect metabolic homeostasis, but the mechanisms by which the rhythmic microbiota influences host circadian activity remain elusive. Using integrated multi-omics and imaging approaches, we demonstrate that the gut microbiota features oscillating biogeographical localization and metabolome patterns that determine the rhythmic exposure of the intestinal epithelium to different bacterial species and their metabolites over the course of a day. This diurnal microbial behavior drives, in turn, the global programming of the host circadian transcriptional, epigenetic, and metabolite oscillations. Surprisingly, disruption of homeostatic microbiome rhythmicity not only abrogates normal chromatin and transcriptional oscillations of the host, but also incites genome-wide de novo oscillations in both intestine and liver, thereby impacting diurnal fluctuations of host physiology and disease susceptibility. As such, the rhythmic biogeography and metabolome of the intestinal microbiota regulates the temporal organization and functional outcome of host transcriptional and epigenetic programs.

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.

The 'TranSeq' 3'-end sequencing method for high-throughput transcriptomics and gene space refinement in plant genomes.

High-throughput RNA sequencing has proven invaluable not only to explore gene expression but also for both gene prediction and genome annotation. However, RNA sequencing, carried out on tens or even hundreds of samples, requires easy and cost-effective sample preparation methods using minute RNA amounts. Here, we present TranSeq, a high-throughput 3'-end sequencing procedure that requires 10- to 20-fold fewer sequence reads than the current transcriptomics procedures. TranSeq significantly reduces costs and allows a greater increase in size of sample sets analyzed in a single experiment. Moreover, in comparison with other 3'-end sequencing methods reported to date, we demonstrate here the reliability and immediate applicability of TranSeq and show that it not only provides accurate transcriptome profiles but also produces precise expression measurements of specific gene family members possessing high sequence similarity. This is difficult to achieve in standard RNA-seq methods, in which sequence reads cover the entire transcript. Furthermore, mapping TranSeq reads to the reference tomato genome facilitated the annotation of new transcripts improving >45% of the existing gene models. Hence, we anticipate that using TranSeq will boost large-scale transcriptome assays and increase the spatial and temporal resolution of gene expression data, in both model and non-model plant species. Moreover, as already performed for tomato (ITAG3.0; www.solgenomics.net), we strongly advocate its integration into current and future genome annotations.

Deterministic direct reprogramming of somatic cells to pluripotency.

Somatic cells can be inefficiently and stochastically reprogrammed into induced pluripotent stem (iPS) cells by exogenous expression of Oct4 (also called Pou5f1), Sox2, Klf4 and Myc (hereafter referred to as OSKM). The nature of the predominant rate-limiting barrier(s) preventing the majority of cells to successfully and synchronously reprogram remains to be defined. Here we show that depleting Mbd3, a core member of the Mbd3/NuRD (nucleosome remodelling and deacetylation) repressor complex, together with OSKM transduction and reprogramming in naive pluripotency promoting conditions, result in deterministic and synchronized iPS cell reprogramming (near 100% efficiency within seven days from mouse and human cells). Our findings uncover a dichotomous molecular function for the reprogramming factors, serving to reactivate endogenous pluripotency networks while simultaneously directly recruiting the Mbd3/NuRD repressor complex that potently restrains the reactivation of OSKM downstream target genes. Subsequently, the latter interactions, which are largely depleted during early pre-implantation development in vivo, lead to a stochastic and protracted reprogramming trajectory towards pluripotency in vitro. The deterministic reprogramming approach devised here offers a novel platform for the dissection of molecular dynamics leading to establishing pluripotency at unprecedented flexibility and resolution.

Each cell counts: Hematopoiesis and immunity research in the era of single cell genomics.

Hematopoiesis and immunity are mediated through complex interactions between multiple cell types and states. This complexity is currently addressed following a reductionist approach of characterizing cell types by a small number of cell surface molecular features and gross functions. While the introduction of global transcriptional profiling technologies enabled a more comprehensive view, heterogeneity within sampled populations remained unaddressed, obscuring the true picture of hematopoiesis and immune system function. A critical mass of technological advances in molecular biology and genomics has enabled genome-wide measurements of single cells - the fundamental unit of immunity. These new advances are expected to boost detection of less frequent cell types and fuzzy intermediate cell states, greatly expanding the resolution of current available classifications. This new era of single-cell genomics in immunology research holds great promise for further understanding of the mechanisms and circuits regulating hematopoiesis and immunity in both health and disease. In the near future, the accuracy of single-cell genomics will ultimately enable precise diagnostics and treatment of multiple hematopoietic and immune related diseases.

Distinct biological events generated by ECM proteolysis by two homologous collagenases.

It is well established that the expression profiles of multiple and possibly redundant matrix-remodeling proteases (e.g., collagenases) differ strongly in health, disease, and development. Although enzymatic redundancy might be inferred from their close similarity in structure, their in vivo activity can lead to extremely diverse tissue-remodeling outcomes. We observed that proteolysis of collagen-rich natural extracellular matrix (ECM), performed uniquely by individual homologous proteases, leads to distinct events that eventually affect overall ECM morphology, viscoelastic properties, and molecular composition. We revealed striking differences in the motility and signaling patterns, morphology, and gene-expression profiles of cells interacting with natural collagen-rich ECM degraded by different collagenases. Thus, in contrast to previous notions, matrix-remodeling systems are not redundant and give rise to precise ECM-cell crosstalk. Because ECM proteolysis is an abundant biochemical process that is critical for tissue homoeostasis, these results improve our fundamental understanding its complexity and its impact on cell behavior.

Chronic exposure to TGFß1 regulates myeloid cell inflammatory response in an IRF7-dependent manner.

Tissue microenvironment influences the function of resident and infiltrating myeloid-derived cells. In the central nervous system (CNS), resident microglia and freshly recruited infiltrating monocyte-derived macrophages (mo-MΦ) display distinct activities under pathological conditions, yet little is known about the microenvironment-derived molecular mechanism that regulates these differences. Here, we demonstrate that long exposure to transforming growth factor-ß1 (TGFß1) impaired the ability of myeloid cells to acquire a resolving anti-inflammatory phenotype. Using genome-wide expression analysis and chromatin immunoprecipitation followed by next-generation sequencing, we show that the capacity to undergo pro- to anti-inflammatory (M1-to-M2) phenotype switch is controlled by the transcription factor interferon regulatory factor 7 (IRF7) that is down-regulated by the TGFß1 pathway. RNAi-mediated perturbation of Irf7 inhibited the M1-to-M2 switch, while IFNß1 (an IRF7 pathway activator) restored it. In vivo induction of Irf7 expression in microglia, following spinal cord injury, reduced their pro-inflammatory activity. These results highlight the key role of tissue-specific environmental factors in determining the fate of resident myeloid-derived cells under both physiological and pathological conditions.

Modeling T cell temporal response to cancer immunotherapy rationalizes development of combinatorial treatment protocols.

Successful immunotherapy relies on triggering complex responses involving T cell dynamics in tumors and the periphery. Characterizing these responses remains challenging using static human single-cell atlases or mouse models. To address this, we developed a framework for in vivo tracking of tumor-specific CD8+ T cells over time and at single-cell resolution. Our tools facilitate the modeling of gene program dynamics in the tumor microenvironment (TME) and the tumor-draining lymph node (tdLN). Using this approach, we characterize two modes of anti-programmed cell death protein 1 (PD-1) activity, decoupling induced differentiation of tumor-specific activated precursor cells from conventional type 1 dendritic cell (cDC1)-dependent proliferation and recruitment to the TME. We demonstrate that combining anti-PD-1 therapy with anti-4-1BB agonist enhances the recruitment and proliferation of activated precursors, resulting in tumor control. These data suggest that effective response to anti-PD-1 therapy is dependent on sufficient influx of activated precursor CD8+ cells to the TME and highlight the importance of understanding system-level dynamics in optimizing immunotherapies.

Physically interacting beta-delta pairs in the regenerating pancreas revealed by single-cell sequencing.

OBJECTIVES: Until recently, communication between neighboring cells in islets of Langerhans was overlooked by genomic technologies, which require rigorous tissue dissociation into single cells. METHODS: We utilize sorting of physically interacting cells (PICs) with single-cell RNA-sequencing to systematically map cellular interactions in the endocrine pancreas after pancreatectomy. RESULTS: The pancreas cellular landscape features pancreatectomy associated heterogeneity of beta-cells, including an interaction-specific program between paired beta and delta-cells. CONCLUSIONS: Our analysis suggests that the particular cluster of beta-cells that pairs with delta-cells benefits from stress protection, implying that the interaction between beta- and delta-cells might safeguard against pancreatectomy associated challenges. The work encourages testing the potential relevance of physically-interacting beta-delta-cells also in diabetes mellitus.

Dissection of floral transition by single-meristem transcriptomes at high temporal resolution.

The vegetative-to-floral transition is a dramatic developmental change of the shoot apical meristem, promoted by the systemic florigen signal. However, poor molecular temporal resolution of this dynamic process has precluded characterization of how meristems respond to florigen induction. Here, we develop a technology that allows sensitive transcriptional profiling of individual shoot apical meristems. Computational ordering of hundreds of tomato samples reconstructed the floral transition process at fine temporal resolution and uncovered novel short-lived gene expression programs that are activated before flowering. These programs are annulled only when both florigen and a parallel signalling pathway are eliminated. Functional screening identified genes acting at the onset of pre-flowering programs that are involved in the regulation of meristem morphogenetic changes but dispensable for the timing of floral transition. Induced expression of these short-lived transition-state genes allowed us to determine their genetic hierarchies and to bypass the need for the main flowering pathways. Our findings illuminate how systemic and autonomous pathways are integrated to control a critical developmental switch.

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.

Inquiring into the differential action of interferons (IFNs): an IFN-alpha2 mutant with enhanced affinity to IFNAR1 is functionally similar to IFN-beta.

Alpha and beta interferons (IFN-alpha and IFN-beta) are multifunctional cytokines that exhibit differential activities through a common receptor composed of the subunits IFNAR1 and IFNAR2. Here we combined biophysical and functional studies to explore the mechanism that allows the alpha and beta IFNs to act differentially. For this purpose, we have engineered an IFN-alpha2 triple mutant termed the HEQ mutant that mimics the biological properties of IFN-beta. Compared to wild-type (wt) IFN-alpha2, the HEQ mutant confers a 30-fold higher binding affinity towards IFNAR1, comparable to that measured for IFN-beta, resulting in a much higher stability of the ternary complex as measured on model membranes. The HEQ mutant, like IFN-beta, promotes a differentially higher antiproliferative effect than antiviral activity. Both bring on a down-regulation of the IFNAR2 receptor upon induction, confirming an increased ternary complex stability of the plasma membrane. Oligonucleotide microarray experiments showed similar gene transcription profiles induced by the HEQ mutant and IFN-beta and higher levels of gene induction or repression than those for wt IFN-alpha2. Thus, we show that the differential activities of IFN-beta are directly related to the binding affinity for IFNAR1. Conservation of the residues mutated in the HEQ mutant within IFN-alpha subtypes suggests that IFN-alpha has evolved to bind IFNAR1 weakly, apparently to sustain differential levels of biological activities compared to those induced by IFN-beta.

MARS-seq2.0: an experimental and analytical pipeline for indexed sorting combined with single-cell RNA sequencing.

Human tissues comprise trillions of cells that populate a complex space of molecular phenotypes and functions and that vary in abundance by 4-9 orders of magnitude. Relying solely on unbiased sampling to characterize cellular niches becomes infeasible, as the marginal utility of collecting more cells diminishes quickly. Furthermore, in many clinical samples, the relevant cell types are scarce and efficient processing is critical. We developed an integrated pipeline for index sorting and massively parallel single-cell RNA sequencing (MARS-seq2.0) that builds on our previously published MARS-seq approach. MARS-seq2.0 is based on >1 million cells sequenced with this pipeline and allows identification of unique cell types across different tissues and diseases, as well as unique model systems and organisms. Here, we present a detailed step-by-step procedure for applying the method. In the improved procedure, we combine sub-microliter reaction volumes, optimization of enzymatic mixtures and an enhanced analytical pipeline to substantially lower the cost, improve reproducibility and reduce well-to-well contamination. Data analysis combines multiple layers of quality assessment and error detection and correction, graphically presenting key statistics for library complexity, noise distribution and sequencing saturation. Importantly, our combined FACS and single-cell RNA sequencing (scRNA-seq) workflow enables intuitive approaches for depletion or enrichment of cell populations in a data-driven manner that is essential to efficient sampling of complex tissues. The experimental protocol, from cell sorting to a ready-to-sequence library, takes 2-3 d. Sequencing and processing the data through the analytical pipeline take another 1-2 d.

DC Respond to Cognate T Cell Interaction in the Antigen-Challenged Lymph Node.

Dendritic cells (DC) are unrivaled in their potential to prime naive T cells by presenting antigen and providing costimulation. DC are furthermore believed to decode antigen context by virtue of pattern recognition receptors and to polarize T cells through cytokine secretion toward distinct effector functions. Diverse polarized T helper (TH) cells have been explored in great detail. In contrast, studies of instructing DC have to date largely been restricted to in vitro settings or adoptively transferred DC. Here we report efforts to unravel the DC response to cognate T cell encounter in antigen-challenged lymph nodes (LN). Mice engrafted with antigen-specific T cells were immunized with nanoparticles (NP) entrapping adjuvants and absorbed with antigen to study the immediate DC response to T cell encounter using bulk and single cell RNA-seq profiling. NP induced robust antigen-specific TH1 cell responses with minimal bystander activation. Fluorescent-labeled NP allowed identification of antigen-carrying DC and focus on transcriptional changes in DC that encounter T cells. Our results support the existence of a bi-directional crosstalk between DC and T cells that promotes TH1 responses, including involvement of the ubiquitin-like molecule Isg15 that merits further study.