To identify the cells at the origin of melanoma, we combined single-cell lineage-tracing and transcriptomics approaches with time-lapse imaging. A mouse model that recapitulates key histopathological features of human melanomagenesis was created by inducing a BRafV600E-driven melanomagenic program in tail interfollicular melanocytes. Most targeted mature, melanin-producing melanocytes expanded clonally within the epidermis before losing their differentiated features through transcriptional reprogramming and eventually invading the dermis. Tumors did not form within interscales, which contain both mature and dormant amelanotic melanocytes. The hair follicle bulge, which contains melanocyte stem cells, was also refractory to melanomagenesis. These studies identify varying tumor susceptibilities within the melanocytic lineage, highlighting pigment-producing cells as the melanoma cell of origin, and indicate that regional variation in tumor predisposition is dictated by microenvironmental cues rather than intrinsic differences in cellular origin. Critically, this work provides in vivo evidence that differentiated somatic cells can be reprogrammed into cancer initiating cells.

The phenotypic and functional dichotomy between IRF8+ type 1 and IRF4+ type 2 conventional dendritic cells (cDC1s and cDC2s, respectively) is well accepted; it is unknown how robust this dichotomy is under inflammatory conditions, when additionally monocyte-derived cells (MCs) become competent antigen-presenting cells (APCs). Using single-cell technologies in models of respiratory viral infection, we found that lung cDC2s acquired expression of the Fc receptor CD64 shared with MCs and of IRF8 shared with cDC1s. These inflammatory cDC2s (inf-cDC2s) were superior in inducing CD4+ T helper (Th) cell polarization while simultaneously presenting antigen to CD8+ T cells. When carefully separated from inf-cDC2s, MCs lacked APC function. Inf-cDC2s matured in response to cell-intrinsic Toll-like receptor and type 1 interferon receptor signaling, upregulated an IRF8-dependent maturation module, and acquired antigens via convalescent serum and Fc receptors. Because hybrid inf-cDC2s are easily confused with monocyte-derived cells, their existence could explain why APC functions have been attributed to MCs.

Modulation and depletion strategies of regulatory T cells (Tregs) constitute valid approaches in antitumor immunotherapy but suffer from severe adverse effects due to their lack of selectivity for the tumor-infiltrating (ti-)Treg population, indicating the need for a ti-Treg specific biomarker.

Cancers acquire several capabilities to survive the multistep process in carcinogenesis. Resisting cell death is one of them. Silencing of the necroptosis initiator Ripk3 occurs in a wide variety of cancer types including melanoma. Little is known about the role of the necroptosis executioner MLKL in tumor development. Studies often indicate opposing roles for MLKL as a tumor-suppressing or a tumor-promoting protein. This study investigates the role of MLKL during melanoma initiation and progression using a tamoxifen-inducible melanoma mouse model driven by melanocyte-specific overexpression of mutated Braf and simultaneous deletion of Pten (BrafV600EPten-/-). In this model we observed a clear sex difference: melanoma initiation and progression were faster in females mice. Mlkl deficiency in male mice resulted in a modest but significant reduction of nevi growth rate compared to the littermate control. In these mice, infiltration and expansion of melanoma cells in the inguinal lymph node were also modestly decreased. This is likely to be a consequence of the delay in nevi development. No significant difference was observed in the Mlkl-deficient condition in female mice in which melanoma development was faster. Overall, our results indicate that in this genetic model MLKL has a minor role during melanoma initiation and progression.

The pancreas is comprised of exocrine and endocrine compartments releasing digestive enzymes into the duodenum and regulating blood glucose levels by insulin and glucagon release. Tissue homeostasis is depending on transcription factor networks, involving Ptf1α, Ngn3, Nkx6.1, and Sox9, which are already activated during organogenesis. However, proper organ function is challenged by diets of high sugar and fat content, increasing the risk of type 2 diabetes and other disorders. A detailed understanding of processes that are important for homeostasis and are impaired during type 2 diabetes is lacking. Here, we show that Zeb1-a transcription factor known for its pivotal role in epithelial-mesenchymal transition, cell plasticity, and metastasis in cancer-is expressed at low levels in epithelial cells of the pancreas and is crucial for organogenesis and pancreas function. Loss of Zeb1 in these cells result in an increase of islet mass, impaired glucose tolerance, and sensitizes to develop liver and pancreas steatosis during diabetes and obesity. Interestingly, moderate overexpression of Zeb1 results in severe pancreas agenesis and lethality after birth, due to islet insufficiency and lack of acinar structures. We show that Zeb1 induction interferes with proper differentiation, cell survival, and proliferation during pancreas formation, due to deregulated expression of endocrine-specific transcription factors. In summary, our analysis suggests a novel role of Zeb1 for homeostasis in epithelial cells that is indispensable for pancreas morphogenesis and proper organ function involving a tight regulation of Zeb1 expression.

Stromal cells support epithelial cell and immune cell homeostasis and play an important role in inflammatory bowel disease (IBD) pathogenesis. Here, we quantify the stromal response to inflammation in pediatric IBD and reveal subset-specific inflammatory responses across colon segments and intestinal layers. Using data from a murine dynamic gut injury model and human ex vivo transcriptomic, protein and spatial analyses, we report that PDGFRA+CD142-/low fibroblasts and monocytes/macrophages co-localize in the intestine. In primary human fibroblast-monocyte co-cultures, intestinal PDGFRA+CD142-/low fibroblasts foster monocyte transition to CCR2+CD206+ macrophages through granulocyte-macrophage colony-stimulating factor (GM-CSF). Monocyte-derived CCR2+CD206+ cells from co-cultures have a phenotype similar to intestinal CCR2+CD206+ macrophages from newly diagnosed pediatric IBD patients, with high levels of PD-L1 and low levels of GM-CSF receptor. The study describes subset-specific changes in stromal responses to inflammation and suggests that the intestinal stroma guides intestinal macrophage differentiation.

Natural killer (NK) cell maturation is a tightly controlled process that endows NK cells with functional competence and the capacity to recognize target cells. Here, we found that the transcription factor (TF) Zeb2 was the most highly induced TF during NK cell maturation. Zeb2 is known to control epithelial to mesenchymal transition, but its role in immune cells is mostly undefined. Targeted deletion of Zeb2 resulted in impaired NK cell maturation, survival, and exit from the bone marrow. NK cell function was preserved, but mice lacking Zeb2 in NK cells were more susceptible to B16 melanoma lung metastases. Reciprocally, ectopic expression of Zeb2 resulted in a higher frequency of mature NK cells in all organs. Moreover, the immature phenotype of Zeb2(-/-) NK cells closely resembled that of Tbx21(-/-) NK cells. This was caused by both a dependence of Zeb2 expression on T-bet and a probable cooperation of these factors in gene regulation. Transgenic expression of Zeb2 in Tbx21(-/-) NK cells partially restored a normal maturation, establishing that timely induction of Zeb2 by T-bet is an essential event during NK cell differentiation. Finally, this novel transcriptional cascade could also operate in human as T-bet and Zeb2 are similarly regulated in mouse and human NK cells.

Heterogeneity between different macrophage populations has become a defining feature of this lineage. However, the conserved factors defining macrophages remain largely unknown. The transcription factor ZEB2 is best described for its role in epithelial to mesenchymal transition; however, its role within the immune system is only now being elucidated. We show here that Zeb2 expression is a conserved feature of macrophages. Using Clec4f-cre, Itgax-cre, and Fcgr1-cre mice to target five different macrophage populations, we found that loss of ZEB2 resulted in macrophage disappearance from the tissues, coupled with their subsequent replenishment from bone-marrow precursors in open niches. Mechanistically, we found that ZEB2 functioned to maintain the tissue-specific identities of macrophages. In Kupffer cells, ZEB2 achieved this by regulating expression of the transcription factor LXRα, removal of which recapitulated the loss of Kupffer cell identity and disappearance. Thus, ZEB2 expression is required in macrophages to preserve their tissue-specific identities.

The metabolic principles underlying the differences between follicular and marginal zone B cells (FoB and MZB, respectively) are not well understood. Here we show, by studying mice with B cell-specific ablation of the catalytic subunit of glutamate cysteine ligase (Gclc), that glutathione synthesis affects homeostasis and differentiation of MZB to a larger extent than FoB, while glutathione-dependent redox control contributes to the metabolic dependencies of FoB. Specifically, Gclc ablation in FoB induces metabolic features of wild-type MZB such as increased ATP levels, glucose metabolism, mTOR activation, and protein synthesis. Furthermore, Gclc-deficient FoB have a block in the mitochondrial electron transport chain (ETC) due to diminished complex I and II activity and thereby accumulate the tricarboxylic acid cycle metabolite succinate. Finally, Gclc deficiency hampers FoB activation and antibody responses in vitro and in vivo, and induces susceptibility to viral infections. Our results thus suggest that Gclc is required to ensure the development of MZB, the mitochondrial ETC integrity in FoB, and the efficacy of antiviral humoral immunity.

ZEB1 and ZEB2 play pivotal roles in solid cancer metastasis by allowing cancer cells to invade and disseminate through the transcriptional regulation of epithelial-to-mesenchymal transition. ZEB expression is also associated with the acquisition of cancer stem cell properties and therapy resistance. Consequently, expression levels of ZEB1/2 and of their direct target genes are widely seen as reliable prognostic markers for solid tumor aggressiveness and cancer patient outcome. Recent loss-of-function mouse models demonstrated that both ZEBs are also essential hematopoietic transcription factors governing blood lineage commitment and fidelity. Interestingly, both gain- and loss-of-function mutations have been reported in multiple hematological malignancies. Combined with emerging functional studies, these data suggest that ZEB1 and ZEB2 can act as tumor suppressors and/or oncogenes in blood borne malignancies, depending on the cellular context. Here, we review these novel insights and discuss how balanced expression of ZEB proteins may be essential to safeguard the functionality of the immune system and prevent leukemia.

GM-CSF promotes myelopoiesis and inflammation, and GM-CSF blockade is being evaluated as a treatment for COVID-19-associated hyperinflammation. Alveolar GM-CSF is, however, required for monocytes to differentiate into alveolar macrophages (AMs) that control alveolar homeostasis. By mapping cross-species AM development to clinical lung samples, we discovered that COVID-19 is marked by defective GM-CSF-dependent AM instruction and accumulation of pro-inflammatory macrophages. In a multi-center, open-label RCT in 81 non-ventilated COVID-19 patients with respiratory failure, we found that inhalation of rhu-GM-CSF did not improve mean oxygenation parameters compared with standard treatment. However, more patients on GM-CSF had a clinical response, and GM-CSF inhalation induced higher numbers of virus-specific CD8 effector lymphocytes and class-switched B cells, without exacerbating systemic hyperinflammation. This translational proof-of-concept study provides a rationale for further testing of inhaled GM-CSF as a non-invasive treatment to improve alveolar gas exchange and simultaneously boost antiviral immunity in COVID-19. This study is registered at ClinicalTrials.gov (NCT04326920) and EudraCT (2020-001254-22).

In Arabidopsis thaliana, brassinosteroid (BR) signaling and stomatal development are connected through the SHAGGY/GSK3-like kinase BR INSENSITIVE2 (BIN2). BIN2 is a key negative regulator of BR signaling but it plays a dual role in stomatal development. BIN2 promotes or restricts stomatal asymmetric cell division (ACD) depending on its subcellular localization, which is regulated by the stomatal lineage-specific scaffold protein POLAR. BRs inactivate BIN2, but how they govern stomatal development remains unclear. Mapping the single-cell transcriptome of stomatal lineages after triggering BR signaling with either exogenous BRs or the specific BIN2 inhibitor, bikinin, revealed that the two modes of BR signaling activation generate spatiotemporally distinct transcriptional responses. We established that BIN2 is always sensitive to the inhibitor but, when in a complex with POLAR and its closest homolog POLAR-LIKE1, it becomes protected from BR-mediated inactivation. Subsequently, BR signaling in ACD precursors is attenuated, while it remains active in epidermal cells devoid of scaffolds and undergoing differentiation. Our study demonstrates how scaffold proteins contribute to cellular signal specificity of hormonal responses in plants.

Epithelial mesenchymal transition (EMT) is tightly regulated by a network of transcription factors (EMT-TFs). Among them is the nuclear factor ZEB2, a member of the zinc-finger E-box binding homeobox family. ZEB2 nuclear localization has been identified in several cancer types, and its overexpression is correlated with the malignant progression. ZEB2 transcriptionally represses epithelial genes, such as E-cadherin (CDH1), by directly binding to the promoter of the genes it regulates and activating mesenchymal genes by a mechanism in which there is no full agreement. Recent studies showed that EMT-TFs interact with epigenetic regulatory enzymes that alter the epigenome, thereby providing another level of control. The role of epigenetic regulation on ZEB2 function is not well understood. In this study, we aimed to characterize the epigenetic effect of ZEB2 repressive function on the regulation of a small Rab GTPase RAB25.

Metabolic-associated fatty liver disease (MAFLD) represents a spectrum of disease states ranging from simple steatosis to non-alcoholic steatohepatitis (NASH). Hepatic macrophages, specifically Kupffer cells (KCs), are suggested to play important roles in the pathogenesis of MAFLD through their activation, although the exact roles played by these cells remain unclear. Here, we demonstrated that KCs were reduced in MAFLD being replaced by macrophages originating from the bone marrow. Recruited macrophages existed in two subsets with distinct activation states, either closely resembling homeostatic KCs or lipid-associated macrophages (LAMs) from obese adipose tissue. Hepatic LAMs expressed Osteopontin, a biomarker for patients with NASH, linked with the development of fibrosis. Fitting with this, LAMs were found in regions of the liver with reduced numbers of KCs, characterized by increased Desmin expression. Together, our data highlight considerable heterogeneity within the macrophage pool and suggest a need for more specific macrophage targeting strategies in MAFLD.

Brain-immune cross-talk and neuroinflammation critically shape brain physiology in health and disease. A detailed understanding of the brain immune landscape is essential for developing new treatments for neurological disorders. Single-cell technologies offer an unbiased assessment of the heterogeneity, dynamics and functions of immune cells. Here we provide a protocol that outlines all the steps involved in performing single-cell multi-omic analysis of the brain immune compartment. This includes a step-by-step description on how to microdissect the border regions of the mouse brain, together with dissociation protocols tailored to each of these tissues. These combine a high yield with minimal dissociation-induced gene expression changes. Next, we outline the steps involved for high-dimensional flow cytometry and droplet-based single-cell RNA sequencing via the 10x Genomics platform, which can be combined with cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) and offers a higher throughput than plate-based methods. Importantly, we detail how to implement CITE-seq with large antibody panels to obtain unbiased protein-expression screening coupled to transcriptome analysis. Finally, we describe the main steps involved in the analysis and interpretation of the data. This optimized workflow allows for a detailed assessment of immune cell heterogeneity and activation in the whole brain or specific border regions, at RNA and protein level. The wet lab workflow can be completed by properly trained researchers (with basic proficiency in cell and molecular biology) and takes between 6 and 11 h, depending on the chosen procedures. The computational analysis requires a background in bioinformatics and programming in R.

Reverse transcription quantitative PCR (RT-qPCR) is the gold standard method for gene expression analysis on mRNA level. To remove experimental variation, expression levels of the gene of interest are typically normalized to the expression level of stably expressed endogenous reference genes. Identifying suitable reference genes and determining the optimal number of reference genes should precede each quantification study. Popular reference genes are not necessarily stably expressed in the examined conditions, possibly leading to inaccurate results. Stably and universally expressed repetitive elements (ERE) have previously been shown to be an excellent alternative for normalization using classic reference genes in human and zebrafish samples. Here, we confirm that in mouse tissues, EREs are broadly applicable reference targets for RT-qPCR normalization, provided that the RNA samples undergo a thorough DNase treatment. We identified Orr1a0, Rltr2aiap, and Rltr13a3 as the most stably expressed mouse EREs across six different experimental conditions. Therefore, we propose this set of ERE reference targets as good candidates for normalization of RT-qPCR data in a plethora of conditions. The identification of widely applicable stable mouse RT-qPCR reference targets for normalization has great potential to facilitate future murine gene expression studies and improve the validity of RT-qPCR data.

Agonistic αCD40 therapy has been shown to inhibit cancer progression in only a fraction of patients. Understanding the cancer cell-intrinsic and microenvironmental determinants of αCD40 therapy response is therefore crucial to identify responsive patient populations and to design efficient combinatorial treatments. Here, we show that the therapeutic efficacy of αCD40 in subcutaneous melanoma relies on preexisting, type 1 classical dendritic cell (cDC1)-primed CD8+ T cells. However, after administration of αCD40, cDC1s were dispensable for antitumor efficacy. Instead, the abundance of activated cDCs, potentially derived from cDC2 cells, increased and further activated antitumor CD8+ T cells. Hence, distinct cDC subsets contributed to the induction of αCD40 responses. In contrast, lung carcinomas, characterized by a high abundance of macrophages, were resistant to αCD40 therapy. Combining αCD40 therapy with macrophage depletion led to tumor growth inhibition only in the presence of strong neoantigens. Accordingly, treatment with immunogenic cell death-inducing chemotherapy sensitized lung tumors to αCD40 therapy in subcutaneous and orthotopic settings. These insights into the microenvironmental regulators of response to αCD40 suggest that different tumor types would benefit from different combinations of therapies to optimize the clinical application of CD40 agonists.

RNA sequencing has become the gold standard for transcriptome analysis but has an inherent limitation of challenging quantification of low-abundant transcripts. In contrast to microarray technology, RNA sequencing reads are proportionally divided in function of transcript abundance. Therefore, low-abundant RNAs compete against highly abundant - and sometimes non-informative - RNA species.

Multiplexing of samples in single-cell RNA-seq studies allows a significant reduction of the experimental costs, straightforward identification of doublets, increased cell throughput, and reduction of sample-specific batch effects. Recently published multiplexing techniques using oligo-conjugated antibodies or -lipids allow barcoding sample-specific cells, a process called "hashing."

During postnatal life, thymopoiesis depends on the continuous colonization of the thymus by bone-marrow-derived hematopoietic progenitors that migrate through the bloodstream. The current understanding of the nature of thymic immigrants is largely based on data from pre-clinical models. Here, we employed single-cell RNA sequencing (scRNA-seq) to examine the immature postnatal thymocyte population in humans. Integration of bone marrow and peripheral blood precursor datasets identified two putative thymus seeding progenitors that varied in expression of CD7; CD10; and the homing receptors CCR7, CCR9, and ITGB7. Whereas both precursors supported T cell development, only one contributed to intrathymic dendritic cell (DC) differentiation, predominantly of plasmacytoid dendritic cells. Trajectory inference delineated the transcriptional dynamics underlying early human T lineage development, enabling prediction of transcription factor (TF) modules that drive stage-specific steps of human T cell development. This comprehensive dataset defines the expression signature of immature human thymocytes and provides a resource for the further study of human thymopoiesis.