Notch and TLR signaling coordinate monocyte cell fate and inflammation
Abstract
Conventional Ly6Chi monocytes have developmental plasticity for a spectrum of differentiated phagocytes. Here we show, using conditional deletion strategies in a mouse model of Toll-like receptor (TLR) 7-induced inflammation, that the spectrum of developmental cell fates of Ly6Chi monocytes, and the resultant inflammation, is coordinately regulated by TLR and Notch signaling. Cell-intrinsic Notch2 and TLR7-Myd88 pathways independently and synergistically promote Ly6Clo patrolling monocyte development from Ly6Chi monocytes under inflammatory conditions, while impairment in either signaling axis impairs Ly6Clo monocyte development. At the same time, TLR7 stimulation in the absence of functional Notch2 signaling promotes resident tissue macrophage gene expression signatures in monocytes in the blood and ectopic differentiation of Ly6Chi monocytes into macrophages and dendritic cells, which infiltrate the spleen and major blood vessels and are accompanied by aberrant systemic inflammation. Thus, Notch2 is a master regulator of Ly6Chi monocyte cell fate and inflammation in response to TLR signaling.
Introduction
Infectious agents or tissue injury trigger an inflammatory response that aims to eliminate the inciting stressor and restore internal homeostasis (Bonnardel and Guilliams, 2018). The mononuclear phagocyte system (MPS) is an integral part of the inflammatory response and consists of the lineage of monocytes and macrophages (MF) and related tissue-resident cells. A key constituent of this system are monocytes of the major (classic) monocyte subtype, in mice called Ly6Chi monocytes. They originate from progenitor cells in the bone marrow (BM), circulate in peripheral blood (PB) and respond dynamically to changing conditions by differentiation into a spectrum of at least three distinct MPS effector phagocytes: Macrophages, dendritic cells (DC), and monocytes with patrolling behavior (Arazi et al., 2019; Bonnardel and Guilliams, 2018; Chakarov et al., 2019; Gamrekelashvili et al., 2016; Hettinger et al., 2013). The diversity of monocyte differentiation responses is thought to be influenced by environmental signals as well as tissue-derived and cell-autonomous signaling mechanisms to ensure context-specific response patterns of the MPS (Okabe and Medzhitov, 2016). However, the precise mechanisms underlying monocyte cell fate decisions under inflammatory conditions are still not fully understood.
When recruited to inflamed or injured tissues, Ly6Chi monocytes differentiate into MF or DC with a variety of phenotypes and function in a context-dependent-manner and regulate the inflammatory response (Krishnasamy et al., 2017; Xue et al., 2014). However, Ly6Chi monocytes can also convert to a second, minor subpopulation of monocytes with blood vessel patrolling behavior. In mice, these are called Ly6Clo monocytes and express CD43, CD11c and the transcription factors Nr4a1, Pou2f2 (Gamrekelashvili et al., 2016; Patel et al., 2017; Varol et al., 2007; Yona et al., 2013). These monocytes have a long lifespan and remain mostly within blood vessels, where they crawl along the luminal side of blood vessels to monitor endothelial integrity and to orchestrate endothelial repair (Auffray et al., 2007; Carlin et al., 2013; Getzin et al., 2018). Steady-state monocyte conversion occurs in specialized endothelial niches and is regulated by monocyte Notch2 signaling activated by endothelial Notch ligands (Avraham-Davidi et al., 2013; Bianchini et al., 2019; Gamrekelashvili et al., 2016; Varol et al., 2007). Notch signaling is a cell-contact-dependent signaling pathway regulating cell fate decisions in the innate immune system (Radtke et al., 2013). Notch signaling regulates formation of intestinal CD11c+CX3CR1+ immune cells (Ishifune et al., 2014), Kupffer cells (Bonnardel et al., 2019; Sakai et al., 2019) and macrophage differentiation from Ly6Chi monocytes in ischemia (Krishnasamy et al., 2017), but also development of conventional DCs (Caton et al., 2007; Epelman et al., 2014; Lewis et al., 2011), which is mediated by Notch2.
Toll-like receptor 7 (TLR7) is a member of the family of pathogen sensors expressed on myeloid cells. Originally identified as recognizing imidazoquinoline derivatives such as Imiquimod (R837) and Resiquimod (R848), TLR7 senses ssRNA, and immune-complexes containing nucleic acids, in a Myd88-dependent manner during virus defense, but is also implicated in tissue-damage recognition and autoimmune disorders (Kawai and Akira, 2010). TLR7-stimulation induces cytokine-production in both mouse and human patrolling monocytes and mediates sensing and disposal of damaged endothelial cells by Ly6Clo monocytes (Carlin et al., 2013; Cros et al., 2010), while chronic TLR7-stimulation drives differentiation of Ly6Chi monocytes into specialized macrophages and anemia development (Akilesh et al., 2019). Furthermore, systemic stimulation with TLR7 agonists induces progressive phenotypic changes in Ly6Chi monocytes consistent with conversion to Ly6Clo monocytes, suggesting involvement of TLR7 in monocyte conversion (Santiago-Raber et al., 2011). Here, we show that Notch signaling alters TLR-driven inflammation and modulates Ly6Clo monocyte vs. macrophage cell fate decisions in inflammation.
Results
TLR and Notch signaling promote monocyte conversion
We first studied the effects of TLR and/or Notch stimulation on monocyte conversion in a defined in vitro system (Gamrekelashvili et al., 2016). Ly6Chi monocytes isolated from the bone marrow (Figure 1—figure supplement 1A and B) of Cx3cr1gfp/+ reporter mice (GFP+) were cultured with recombinant Notch ligand Delta-like 1 (DLL1) in the presence or absence of the TLR7/8 agonist R848 and analyzed after 24 hr for the acquisition of key features of Ly6Clo monocytes (Gamrekelashvili et al., 2016; Hettinger et al., 2013). In contrast to control conditions, cells cultured with DLL1 showed an upregulation of CD11c and CD43, remained mostly MHC-II negative, and expressed transcription factors Nr4a1 and Pou2f2, markers for Ly6Clo monocytes, leading to a significant, five-fold increase of Ly6Clo cells, consistent with enhanced monocyte conversion. Cells cultured with R848 alone showed a comparable phenotype response, both qualitatively and quantitatively (Figure 1A–C). Interestingly, on a molecular level, R848 stimulation primarily acted on Pou2f2 induction and CD43 expression, while Notch stimulation primarily induced Nr4a1 and CD11c upregulation. Furthermore, the combination of DLL1 and R848 strongly and significantly increased the number of CD11c+CD43+ Ly6Clo cells above the level of individual stimulation and significantly enhanced expression levels of both transcriptional regulators Nr4a1 and Pou2f2 (Figure 1A–C), suggesting in part synergistic and/or cumulative regulation of monocyte conversion by TLR7/8 and Notch signaling. By comparison, the TLR4 ligand LPS also increased Ly6Clo cell numbers and expression levels of Nr4a1 and Pou2f2. However, the absolute conversion rate was lower under LPS and there was no synergy/cumulative effect seen with DLL1 (Figure 1D and E).
Since monocyte conversion is regulated by Notch2 in vitro and in vivo (Gamrekelashvili et al., 2016), we next tested TLR-induced conversion in Ly6Chi monocytes with Lyz2Cre-mediated conditional deletion of Notch2 (N2ΔMy). Both, littermate control (wt) and N2ΔMy monocytes showed comparable response to R848, but conversion in the presence of DLL1, and importantly, also DLL1-R848 co-stimulation was significantly impaired in knock-out cells (Figure 1F). This suggests independent contributions of TLR and Notch signaling to monocyte conversion.
To study whether the TLR stimulation requires Myd88 we next tested purified Ly6Chi monocytes (Figure 1—figure supplement 1C and D) with Myd88 loss-of-function (Myd88-/-). Compared to wt cells, Myd88-/- monocytes showed strongly impaired conversion in response to R848 but a conserved response to DLL1. The response to DLL1-R848 co-stimulation, however, was significantly impaired (Figure 1G). Furthermore, expression of Nr4a1 and Pou2f2 by R848 was strongly reduced in Myd88-/- monocytes with or without DLL1 co-stimulation, while DLL1-dependent induction was preserved (Figure 1H). Thus, Notch and TLR signaling act independently and synergistically to promote monocyte conversion.
To address the role of TLR stimulation for monocyte conversion in vivo we adoptively transferred sorted Ly6Chi monocytes from CD45.2+GFP+ mice into CD45.1+ congenic recipients, injected a single dose of R848 and analyzed transferred CD45.2+GFP+ cells in BM and Spl after 2 days (Figure 2A and Figure 1—figure supplement 1A and B). Stimulation with R848 significantly promoted conversion into Ly6Clo monocytes displaying the proto-typical Ly6CloCD43+CD11c+MHC-IIlo/- phenotype (Figure 2B and C and Figure 2—figure supplement 1A). In contrast, transfer of Myd88-/- Ly6 Chi monocytes resulted in impaired conversion in response to R848 challenge (Figure 2D and E and Figure 1—figure supplement 1D and Figure 2—figure supplement 1B). Together, these data indicate that TLR and Notch cooperate in the regulation of monocyte conversion.
Notch2-deficient mice show altered myeloid inflammatory response
To characterize the response to TLR stimulation in vivo, we applied the synthetic TLR7 agonist Imiquimod (IMQ, R837) in a commercially available crème formulation (Aldara) daily to the skin of mice (El Malki et al., 2013; van der Fits et al., 2009) and analyzed the systemic inflammatory response in control or N2ΔMy mice (Gamrekelashvili et al., 2016; Figure 3A). While treatment with IMQ-induced comparable transient weight loss and ear swelling in both genotypes (Figure 3—figure supplement 1A), splenomegaly in response to treatment was significantly more pronounced in N2ΔMy mice (Figure 3—figure supplement 1B).
To characterize the spectrum of myeloid cells in more detail, we next performed flow cytometry of PB cells with a dedicated myeloid panel (Gamrekelashvili et al., 2016) in wt or N2ΔMy Cx3cr1gfp/+ mice and subjected live Lin-CD11b+GFP+ subsets to unsupervised t-SNE analysis (Figure 3B). This analysis strategy defined five different populations, based on single surface markers: Ly6C+, CD43+, MHC-II+, F4/80hi and CD11chi (Figure 3C). Applying these five gates to samples from separate experimental conditions identified dynamic alterations in blood myeloid subsets in response to IMQ, but also alterations in N2ΔMy mice (Figure 3D). Specifically, abundance and distribution of Ly6C+ cells, containing classical monocytes, in response to IMQ were changed to the same extend in both genotypes. In contrast, the MHC-II+ and F4/80hi subsets were more abundant in N2ΔMy mice, but also showed more robust changes in response to IMQ. On the other hand, the CD43+ subset, containing the patrolling monocyte subset, showed prominent enrichment in wt mice, but was less abundant and showed diminished distribution changes after IMQ treatment in N2ΔMy mice (Figure 3D).
To analyze the initially defined subsets more precisely, we applied a multi-parameter gating strategy to define conventional cell subsets (Figure 3—figure supplement 2A and B and Supplementary file 1; Gamrekelashvili et al., 2016).
In response to IMQ, Ly6Chi monocytes in wt mice increased transiently in blood, and this response was not altered in mice with conditional Notch2 loss-of function (Figure 3E). In contrast, while Ly6Clo monocytes robustly increased over time with IMQ treatment in wt mice, their levels in N2ΔMy mice were lower at baseline (Gamrekelashvili et al., 2016) and remained significantly reduced throughout the whole observation period (Figure 3E and F and Figure 3—figure supplement 2A and B). At the same time, while untreated N2ΔMy mice showed increased levels of MHC-II+ atypical monocytes (Figure 3E and F and Figure 3—figure supplement 2A and B; Gamrekelashvili et al., 2016), IMQ treatment induced the generation of F4/80hiCD115+ monocytes in the blood and increased MF in the spleen at d5 (Figure 3E and F and Figure 3—figure supplement 3A and B). This was followed by a peak in the DC population at d7 (Figure 3E and F and Figure 3—figure supplement 3A and B). These latter changes did not occur in bone marrow but were only observed in the periphery (Figure 3—figure supplement 3C). Together, these data suggest that wt Ly6Chi monocytes convert to Ly6Clo monocytes in response to TLR stimulation, while Notch2 deficient Ly6Chi monocytes differentiate into F4/80hiCD115+ monocytes, macrophages and DC, suggesting Notch2 as a master regulator of Ly6Chi monocyte cell fate during systemic inflammation.
Global gene expression analysis identifies macrophage gene expression signatures in monocytes of Notch2-deficient mice during acute inflammation
To characterize more broadly the gene expression changes involved in monocyte differentiation during inflammation, we next subjected monocyte subsets from PB of wt and N2ΔMy mice after Sham or IMQ treatment (Figure 4—figure supplement 1A) to RNA-sequencing and gene expression analysis. After variance filtering and hierarchical clustering, 600 genes were differentially expressed between six experimental groups (Figure 4A and Figure 4—source data 1).
Principal component analysis (PCA) of differentially expressed genes (DEG) of all experimental groups revealed a clear separation between control Ly6Chi monocytes and IMQ-treated Ly6Chi or Ly6Clo monocytes. Interestingly, the effects of Notch2 loss-of-function were most pronounced in the Ly6Clo populations, which separated quite strongly depending on genotype, while Ly6Chi monocytes from wt and N2ΔMy mice over all maintained close clustering under Sham or IMQ conditions (Figure 4B and C and Figure 4—source data 1).
Furthermore, while wt Ly6Clo cells were enriched for genes characteristic of patrolling monocytes (Hes1, Nr4a1, Ace, Cd274 and Itgb3), cells in the Ly6Clo gate from N2ΔMy mice showed upregulation of genes characteristic of mature phagocytes, such as MF (Fcgr1, Mertk, C1qa, Clec7a, Maf, Cd36, Cd14, Adgre1 (encoding F4/80)) (Figure 4D–F).
Comparative gene expression analysis of Ly6Clo cell subsets during IMQ treatment identified 373 genes significantly up- or down-regulated with Notch2 loss-of-function (p-value<0.01, Figure 4C–F and Figure 4—source data 2), which were enriched for phagosome formation, complement system components, Th1 and Th2 activation pathways and dendritic cell maturation by ingenuity canonical pathway analysis (Figure 4—figure supplement 1B). Notably, signatures for autoimmune disease processes were also enriched (Supplementary file 2). Independent gene set enrichment analysis (GSEA) (Isakoff et al., 2005; Mootha et al., 2003) confirmed consistent up-regulation of gene sets in N2ΔMy Ly6Clo cells involved in several gene ontology biological processes, such as vesicle-mediated transport (GO:0016192), defense response (GO:0006952), inflammatory response (GO:0006954), response to bacterium (GO:0009617) and endocytosis (GO:0006897) (Supplementary file 3 and Figure 4—source data 2). Overall, these data suggest regulation of Ly6Chi monocyte cell fate and inflammatory responses by Notch2.
Furthermore, changes in cell populations resulted in altered systemic inflammatory response patterns. Levels of TLR-induced cytokines and chemokines, such as TNF-α, CXCL1, IL-1β, IFN-α, were elevated to the same extend in wt and N2ΔMy mice in response to IMQ treatment, suggesting normal primary TLR-activation (Figure 4—figure supplement 1C). However, circulating levels of chemokines produced by Ly6Clo monocytes (Carlin et al., 2013), such as CCL2, CCL3, CXCL10, and IL-10 were higher in wt mice compared to N2ΔMy mice, while the levels of pro-inflammatory cytokines IL-17A, IL-6 and GMCSF were significantly enhanced in N2ΔMy mice but not in wt mice as compared to untreated controls, confirming systemic alterations in addition to cellular changes in Notch2 loss-of-function mice in response to IMQ (Figure 4G and Figure 4—figure supplement 1C).
Notch2-deficiency promotes macrophage differentiation
To match the observed gene expression pattern of inflammatory Ly6Clo cells from wt and N2ΔMy mice under IMQ treatment with previously described cells of the monocyte-macrophage lineage we performed pairwise gene set enrichment analysis with defined myeloid cell transcriptomic signatures using the GSEA software and BubbleGUM stand-alone software (Isakoff et al., 2005; Mootha et al., 2003; Spinelli et al., 2015; Figure 5A and Supplementary file 4 and 5 and Figure 4—source data 2). Out of 29 transcriptomic signatures - representing tissue MF, monocyte derived DC (MoDC), conventional DC (cDC), plasmacytoid DC (pDC), classical (Ly6Chi) monocytes (cMonocyte), non-classical (Ly6Clo) monocytes (ncMonocyte) and B cells (as a reference) – significant enrichment was registered in seven signatures (normalized enrichment score (NES >1.5, FDR < 0.1)). The cell fingerprint representing ncMonocyte (#3) was highly enriched in the gene set from wt Ly6Clo monocytes (NES = 1.91, FDR = 0.039), while all other cell fingerprints showed no significant similarity (Figure 5A and Supplementary file 4), confirming a strong developmental restriction toward Ly6Clo monocytes in wt cells. In contrast, N2ΔMy gene sets showed the highest similarity (NES >1.8 and FDR < 0.01) with four cell fingerprints (#1, 5, 6, 7) representing different MF populations, and weak similarity to MoDC (NES = 1.64, FDR < 0.1) and cMonocyte (NES = 1.54, FDR < 0.1) (Figure 5A and Supplementary files 4 and 5). Phenotyping of cell populations by flow cytometry using MF markers MerTK and CD64 (Figure 5B–D) confirmed selective expansion of an F4/80hiMerTK+ (FM+) monocyte population in IMQ-treated N2ΔMy mice (Figure 5A and Supplementary files 4 and 5). Together, these data demonstrate a cell fate switch from Ly6Clo monocytes toward macrophage signatures in the absence of Notch2.
Notch2 regulates monocyte cell fate decisions during inflammation
In the steady-state, Ly6Chi monocytes differentiate into Ly6Clo monocytes and this process is regulated by Notch2 (Gamrekelashvili et al., 2016). In order to confirm that Notch2 controls differentiation potential of Ly6Chi monocytes in response to TLR stimulation, we performed adoptive transfer of CD45.2+ wt or N2ΔMy BM Ly6Chi monocytes into IMQ-treated CD45.1+ congenic recipients and analyzed the fate of donor cells after 3 days (Figure 5—figure supplement 1A). Unsupervised t-SNE analysis of flow cytometry data showed an expanded spectrum of expression patterns in cells from N2ΔMy donors compared to wt controls (Figure 5—figure supplement 1B). More precisely, Ly6Chi monocytes from wt mice converted preferentially to Ly6Clo monocytes (Ly6CloF4/80lo/-CD11c+CD43+MHC-IIlo/- phenotype) during IMQ treatment (Figure 5E and F and Figure 5—figure supplement 1C). In contrast, conversion of Notch2-deficient Ly6Chi to Ly6Clo monocytes was strongly impaired, but the development of donor-derived F4/80hi macrophages in the spleen was strongly enhanced (Figure 5G). Furthermore, expansion of macrophages was also observed in aortas of Notch2 deficient mice in vivo after IMQ treatment (Figure 5H). Adoptive transfer studies confirmed that MF in IMQ-treated aortas originated from N2ΔMy Ly6Chi monocytes (Figure 5J and K).
To evaluate Notch targeting efficiency in our system we performed flow cytometry and transcriptome analysis of transgenic monocyte subsets. In Ly6Chi monocytes, the Lyz2Cre-mediated conditional deletion strategy induced a 40–50% reduction of Notch2 expression by flow cytometry and transcriptome analysis (Figure 5—figure supplement 2A and C), and a > 50% reduction of Notch target gene Hes1, demonstrating partial targeting and functional impairment of Notch2 in Ly6Chi monocytes (Figure 5—figure supplement 2B). At the same time, in Ly6Clo monocytes at baseline there was only minor reduction of Notch2, corroborating earlier results (Gamrekelashvili et al., 2016), but more efficient reduction of Notch2 in IMQ-treated mice to levels seen in Ly6Chi monocytes (Figure 5—figure supplement 2C). Furthermore, levels of Notch target gene Hes1 were low in Ly6Chi monocytes but significantly increased in Ly6Clo monocytes (Figure 5—figure supplement 2B; Gamrekelashvili et al., 2016), suggesting low Notch signaling in Ly6Chi monocytes but higher Notch signaling activity in Ly6Clo monocytes.
Since Ly6Clo monocyte deficiency in Nr4a1-deficient mice is caused by increased apoptosis (Hanna et al., 2011), we analyzed cells by AnnexinV staining. Notch2-deficiency did not significantly increase cell death in all analyzed cell populations, neither during steady state nor during inflammation (Figure 5L and M).
Together, these data confirm that Notch2 is a master regulator of Ly6Chi monocyte differentiation potential, regulating a switch between Ly6Clo monocyte or macrophage cell fate during inflammation. These data also demonstrate that in the context of inactive myeloid Notch2 signaling, TLR-stimulation results in systemic pro-inflammatory changes and vascular inflammation.
Discussion
Together, our data present a spectrum of developmental cell fates of Ly6Chi monocytes and their coordinated regulation by TLR and Notch signaling during inflammation. TLR and Notch signaling act independently and synergistically in promoting Ly6Clo monocyte development from Ly6Chi monocytes, while impairment in either signaling axis impairs Ly6Clo monocyte development. On the other hand, TLR stimulation in the absence of functional Notch2 signaling promotes macrophage gene expression signatures in monocytes and development of MF in aorta and spleen, suggesting Notch2 as a master regulator of Ly6Chi monocyte cell fate during systemic inflammation.
Plasticity of Ly6Chi monocytes ensures adaptation to environmental signals, which trigger distinct cell developmental programs inducing context- or tissue-specific subsets of terminally differentiated phagocytes, including Ly6Clo monocytes, MF or DC (Guilliams et al., 2018). In the steady-state, a subset of Ly6Chi monocytes converts to Ly6Clo monocytes in mice and humans, which is regulated by Notch2 and the endothelial Notch ligand Delta-like 1 (Dll1) (Gamrekelashvili et al., 2016; Patel et al., 2017; Yona et al., 2013). However, when recruited into tissues, Ly6Chi monocytes can give rise to two types of monocyte-derived resident tissue macrophages (MRTM) (Chakarov et al., 2019). These Lyve1hiMHC-IIlo and Lyve1loMHC-IIhi MRTMs differ in phenotype and function as well as spatial distribution. Our gene set enrichment analysis revealed that the signature of cells from inflamed wt mice showed highest and selective similarity to non-classical Ly6Clo monocytes, while Notch2 loss-of-function cells showed the highest similarity to the gene set of Lyve1hiMHC-IIlo MRTM, but also a more general similarity to an extended spectrum of different MF signatures, suggesting Notch2 as a gate keeper of Ly6Clo monocytes vs. macrophage differentiation during inflammation. Due to the low number Ly6Clo monocytes at the steady state, we were not able to compare baseline expression profiles, which is a limitation of this analysis. Lyve1hiMHC-IIlo interstitial MFs are closely associated with blood vessels across different tissues and mediate inflammatory reactions (Chakarov et al., 2019). In line with this, Notch2 knock-out mice showed a population of FM+ monocytes with partial gene expression signatures of MF circulating in the blood and increased MF in the aorta, and adoptive transfer studies of Ly6Chi monocytes successfully recapitulated their differentiation into aortic MF. These data also have implications for the potential developmental regulation of MRTMs by TLR and Notch2. At the same time, the Notch2-deficient population showed weaker but significant enrichment of MoDC signatures, suggesting a mixture of cell subsets representing different stages of monocyte differentiation (Menezes et al., 2016; Mildner et al., 2017) or lineage commitment (Liu et al., 2019; Yanez et al., 2017) within this cell pool, although formally we cannot exclude progenitor contamination as a confounder in our adoptive cell transfer studies.
The Lyz2Cre-mediated conditional deletion strategy induced a 40–50% reduction of Notch2 expression, and a > 50% reduction of Notch target gene Hes1 in Ly6Chi monocytes, suggesting sufficient targeting and functional impairment of Notch2. Nevertheless, partial targeting might explain the small differences seen in PCA analysis in these cells. However, in light that baseline Notch signaling activity seems to be low in Ly6Chi monocytes and only significantly increases in Ly6Clo monocytes (Gamrekelashvili et al., 2016), it suggests that Notch2 influences cell fate decision in Ly6Chi monocytes or at the early stages of conversion to Ly6Clo monocytes. Furthermore, the fact that there is only a minor reduction of Notch2 in Ly6Clo monocytes at baseline, which suggests a strong selection bias against Notch2 loss-of-function, further argues for a strict Notch2-dependence of monocyte conversion.
In the case of Nr4a1 loss-of-function, the reduced numbers of Ly6Clo monocytes are due to increased apoptosis (Hanna et al., 2011). Although we did not find evidence for increased apoptosis due to Notch2-deficiency, our data do not exclude the possibility that regulation of cell survival by Notch2 contributes to the observed phenotype. In fact, two lines of the evidence suggest that regulation of cell survival might act synergistically to cell fate choices: first, in absolute numbers, there is no compensation by alternative cell fates for the number of lacking Ly6Clo monocytes in the blood or spleen of N2ΔMy mice; second, expression of Bcl2, a strong regulator of cell survival, is downregulated in IMQ-treated N2ΔMy Ly6Clo cells as compared to controls.
While our current data clearly demonstrate that Notch2 loss-of-function promotes macrophage gene expression profiles in monocytes and macrophage development from Ly6Chi monocytes during TLR stimulation, we have previously shown that Dll1-Notch signaling promotes maturation of anti-inflammatory macrophages from Ly6Chi monocytes in ischemic muscle (Krishnasamy et al., 2017). Furthermore, Dll4-Notch signaling initiated in the liver niche was recently shown to promote Kupffer cell development after injury (Bonnardel et al., 2019; Sakai et al., 2019) or to promote pro-inflammatory macrophage development (Xu et al., 2012). This suggests that the role of Notch is ligand-, cell- and context-specific, which emphasizes the differential effects of specific ligand-receptor combinations (Benedito et al., 2009). Our data demonstrate that Notch2 is a master regulator of Ly6Chi monocyte cell fate during inflammation, which contributes to the nature of the inflammatory response.
Lastly, our data also reveal a potentially important function of myeloid Notch2 for regulation of systemic and vascular inflammation with implications for autoimmune disease. When wt mice are challenged with TLR stimulation they show predominant conversion of Ly6Chi monocytes into Ly6Clo monocytes with blood vessel patrolling and repairing function (Carlin et al., 2013) and IL-10 secretion. In contrast, Notch2 knock-out mice show predominant and ectopic differentiation of Ly6Chi monocytes into FM+ monocytes and macrophages, which appear in the bloodstream and the spleen and infiltrate major blood vessels, such as the aorta, along with aberrant cytokine profiles. In addition, absence of functional Notch2 promoted a core macrophage signature and strong upregulation of canonical pathways involved in autoimmune disease. Since TLR7 has been implicated in the development of autoimmune disease (Santiago-Raber et al., 2011; Santiago-Raber et al., 2010), our data suggest Notch2 as an important modulator of this process by regulating cell differentiation and systemic inflammation. However, the relevance and possible disease context requires further studies.
Materials and methods
Mice
B6.129P-Cx3cr1tm1Litt/J (Cx3cr1GFP/+) mice (Jung et al., 2000), B6.129P2-Lyz2tm1(cre)Ifo (Lyz2Cre) mice (Clausen et al., 1999), B6.129-Notch2tm1Frad/J (Notch2lox/lox) mice (Besseyrias et al., 2007), B6.129-Lyz2tm1(cre)IfoNotch2tm1FradCx3cr1tm1Litt (N2ΔMy) (Gamrekelashvili et al., 2016 have been previously described. B6.SJL-PtprcaPepcb/BoyJ (CD45.1+) mice were from central animal facility of Hannover Medical School (ZTL, MHH). B6.129P2-Myd88tm1Hlz/J (Myd88-/-) (Gais et al., 2012) and Myd88+/+ littermate control (wt) mice were kindly provided by Dr. Matthias Lochner. Mice were housed under specific pathogen-free conditions in the animal facility of Hannover Medical School.
Tissue and cell preparation
View detailed protocolFor single cell suspension mice were sacrificed and spleen, bone marrow, blood and aortas were collected. Erythrocytes were removed by red blood cell lysis buffer (BioLegend) or by density gradient centrifugation using Histopaque 1083 (Sigma-Aldrich). Aortas were digested in DMEM medium supplemented with 500 U/ml Collagenase II (Worthington). After extensive washing, cells were resuspended in PBS containing 10%FCS and 2 mM EDTA kept on ice, stained and used for flow cytometry or for sorting.
Flow cytometry and cell sorting
View detailed protocolNon-specific binding of antibodies to Fc-receptors was blocked with anti-mouse CD16/CD32 (TruStain fcX from BioLegend) in single-cell suspensions prepared from Spl, PB or BM. After subsequent washing step, cells were labeled with primary and secondary antibodies or streptavidin-fluorochrome conjugates and were used for flow cytometry analysis (LSR-II, BD Biosciences) or sorting (FACSAria; BD Biosciences or MoFlo XDP; Beckman Coulter).
For apoptosis assay, single-cell suspensions were stained with primary and secondary antibodies, washed, re-suspended in AnnexinV binding buffer (Biolegend) and transferred into tubes. Cells were stained with AnnexinV (AnnV) and propidium iodide (PI) at room temperature for 20 min and were immediately analyzed by flow cytometry. Antibodies and fluorochromes used for flow cytometry are described in Supplementary file 6. Flow cytometry data were analyzed using FlowJo software (FlowJo LLC). Initially cells were identified based on FSC and SSC characteristics. After exclusion of doublets (on the basis of SSC-W, SSC-A), relative frequency of each subpopulation from live cell gate, or absolute number of each subset (calculated from live cell gate and normalized per Spl, per mg Spl, mg BM, mg aorta or µl PB) were determined and are shown in the graphs as mean ± SEM, unless otherwise stated. Unsupervised t-distributed stochastic neighbor embedding (t-SNE) analysis (van der Maaten and Hinton, 2008) was performed on live CD45+Lin-GFP+CD11b+ population in concatenated samples using FlowJo.
Cytokine multiplex bead-based assay
Request a detailed protocolSera were collected from control or Aldara treated mice and kept frozen at −80°C. Concentration of IFN-γ, CXCL1, TNF-α, CCL2, IL-12(p70), CCL5, IL-1β, CXCL10, GM-CSF, IL-10, IFN-β, IFN-α, IL-6, IL-1α, IL-23, CCL3, IL-17A were measured with LEGENDplex multi-analyte flow assay kits (BioLegend) according to manufacturer’s protocol on LSR-II flow cytometer. Data were processed and analyzed with LEGENDplex data analysis software (BioLegend).
RNA isolation, library construction, sequencing and analysis
Request a detailed protocolPeripheral blood monocyte subpopulations were sorted from Aldara treated mice or untreated controls (Figure 5—figure supplement 1A) and RNA was isolated using RNeasy micro kit (Qiagen). Libraries were constructed from total RNA with the ‘SMARTer Stranded Total RNA-Seq Kit v2 – Pico Input Mammalian’ (Takara/Clontech) according to manufacturer’s recommendations, barcoded by dual indexing approach and amplified with 11 cycles of PCR. Fragment length distribution was monitored using ‘Bioanalyzer High Sensitivity DNA Assay’ (Agilent Technologies) and quantified by ‘Qubit dsDNA HS Assay Kit’ (ThermoFisher Scientific). Equal molar amounts of libraries were pooled, denatured with NaOH, diluted to 1.5pM (according to the Denature and Dilute Libraries Guide (Illumina)) and loaded on an Illumina NextSeq 550 sequencer for sequencing using a High Output Flowcell for 75 bp single reads (Illumina). Obtained BCL files were converted to FASTQ files with bcl2fastq Conversion Software version v2.20.0.422 (Illumina). Pre-processing of FASTQ inputs, alignment of the reads and quality control was conducted by nfcore/rnaseq (version 1.5dev) analysis pipeline (The National Genomics Infrastructure, SciLifeLab Stockholm, Sweden) using Nextflow tool. The genome reference and annotation data were taken from GENCODE.org (GRCm38.p6 release M17). Data were normalized with DESeq2 (Galaxy Tool Version 2.11.39) with default settings and output counts were used for further analysis with Qlucore Omics explorer (Sweden). Data were log2 transformed, 1.1 was used as a threshold and low expression genes (<50 reads in all samples) were removed from the analysis. Hierarchical clustering (HC) or principal component analysis (PCA) was performed on 600 differential expressed genes (DEGs) after variance filtering (filtering threshold 0.295) selected by ANOVA with the Benjamini-Hochberg (B-H) correction (p<0.01, FDR ≤ 0.01). For two group comparisons Student’s t-test with B-H correction was used and 373 DEGs ((Variance filtering 0.117,–0.6≥Δlog2 ≥0.6 (p<0.01, FDR ≤ 0.05)) were selected for further IPA or GSEA analysis.
Ingenuity pathway analysis (IPA) was performed on 373 DEGs using IPA software (Qiagen) with default parameters. Top 20 canonical pathways and top five immunological diseases enriched in DEGs were selected for display.
Gene set enrichment analysis (GSEA) (Mootha et al., 2003; Subramanian et al., 2005) was performed on 373 DEGs using GSEA software (Broad institute) and C5 GO biological process gene sets (Liberzon et al., 2015) from MsigDB with 1000 gene set permutations for computing p-values and FDR.
BubbleGUM software, an extension of GSEA (Spinelli et al., 2015; Vu Manh et al., 2015), GSEA software (Broad Institute) and published transcriptomic signatures (Chakarov et al., 2019; Gautier et al., 2012; Schlitzer et al., 2015; Vu Manh et al., 2015) were used to assess and visualize the enrichment of obtained gene sets for myeloid populations and define the nature of the cells from which the transcriptomes were generated.
In vitro conversion studies
Request a detailed protocol96-well flat bottom plates were coated at room temperature for 3 hr with IgG-Fc or DLL1-Fc ligands (all from R and D) reconstituted in PBS. Sorted BM Ly6Chi monocytes were cultured in coated plates and were stimulated with Resiquimod (R848, 0.2 µg ml−1, Cayman Chemicals) or LPS (0.2 µg ml−1, E. coli O55:B5 Sigma-Aldrich) in the presence of M-CSF (10 ng ml−1, Peprotech) at 37°C for 24 hr. One day after culture, cells were harvested, stained and subjected to flow cytometry. Relative frequency (from total live CD11b+GFP+ cells) or absolute numbers of Ly6Clo monocyte-like cells (CD11b+GFP+Ly6Clo/-CD11cloMHC-IIlo/-CD43+) recovered from each well served as an indicator of conversion efficiency and is shown in the graphs. Alternatively, cultured cells were harvested and isolated RNA was used for gene expression analysis.
Induction of acute systemic inflammation using IMQ
View detailed protocolMice were anesthetized and back skin was shaved and depilated using depilating crème. Two days after 50 mg/mouse/day Aldara (containing 5% Imiquimod, from Meda) or Sham crème were applied on depilated skin and right ear (where indicated) for 4–5 consecutive days (El Malki et al., 2013; van der Fits et al., 2009). Mouse weight and ear thickness were monitored daily. Mice were euthanized on the indicated time points after start of treatment (day 0, 5, 7 and 12) PB, Spl, BM and aortas were collected for further analysis.
Adoptive cell transfer experiments
View detailed protocolLin-CD11b+Ly6ChiGFP+ monocytes were sorted from BM of CD45.2+ donors and injected into CD45.1+ recipients intravenously (i.v.). In separate experiment wt or Myd88-/- LinCD45.2+CD11b+Ly6ChiCX3CR1+ monocytes were used as a donor for transfer. 30 min later PBS or R848 (37.5 µg per mouse) were injected in recipient mice. Two days after transfer Spl, PB and BM were collected and single cell suspension was prepared. After blocking of Fc receptors (anti-mouse CD16/CD32, TruStain fcX from BioLegend), cells were labeled with biotin-conjugated antibody cocktail containing anti-CD45.1 and anti-Lin (anti-CD3, CD19, B220, NK1.1, Ly6G, Ter119) antibodies, anti-biotin magnetic beads and enriched on LS columns (Miltenyi Biotec) according to manufacturer’s instructions. CD45.1negLinneg fraction was collected, stained and analyzed by flow cytometry. Ly6Clo monocytes (CD45.2+CD11b+GFP+Ly6Clo/-F4/80loCD11cloCD43hi cells) and macrophages (CD45.2+CD11b+GFP+Ly6Clo/-F4/80hiCD115+ cells) were quantified in Spl, BM and PB as relative frequency of total donor derived CD45.2+CD11b+GFP+ cells. Adoptive transfer experiments in IMQ-treated mice were performed using wt or N2ΔMy Lin-CD45.2+CD11b+Ly6ChiGFP+ donor monocytes and Spl, PB and aortas of recipients were analyzed 3 days after transfer.
Quantitative real-time PCR analysis
Request a detailed protocolTotal RNA was purified from cell lysates using Nucleospin RNA II kit (Macherey Nagel). After purity and quality check, RNA was transcribed into cDNA employing cDNA synthesis kit (Invitrogen) according to manufacturer’s instructions. Quantitative real-time PCR was performed using specific primers for Nr4a1: forward, 5’-AGCTTGGGTGTTGATGTTCC-3’, reverse, 5’-AATGCGATTCTGCAGCTCTT-3’ and Pou2f2: forward, 5’-TGCACATGGAGAAGGAAGTG-3’, reverse, 5’-AGCTTGGGACAATGGTAAGG-3’ and FastStart Essential DNA Green Master on a LightCycler 96 system from Roche according to manufacturer’s instructions. Expression of each specific gene was normalized to expression of Rps9 and calculated by the comparative CT (2-ΔΔCT) method (Schmittgen and Livak, 2008).
Statistical analysis
Request a detailed protocolResults are expressed as mean ± standard error of mean (SEM). N numbers are biological replicates of experiments performed at least three times unless otherwise indicated. Significance of differences was calculated using unpaired, two-tailed Student’s t-test with confidence interval of 95%. For comparison of multiple experimental groups one-way or two-way ANOVA and Bonferroni’s multiple-comparison test was performed.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Data from RNA sequencing have been deposited to NCBI's Gene Expression Omnibus and are available under the accession number GSE147492.
-
NCBI Gene Expression OmnibusID GSE147492. Notch and TLR signaling coordinate monocyte cell fate and inflammation.
References
-
The immune cell landscape in kidneys of patients with lupus nephritisNature Immunology 20:902–914.https://doi.org/10.1038/s41590-019-0398-x
-
On-site education of VEGF-recruited monocytes improves their performance as angiogenic and arteriogenic accessory cellsThe Journal of Experimental Medicine 210:2611–2625.https://doi.org/10.1084/jem.20120690
-
Hierarchy of Notch-Delta interactions promoting T cell lineage commitment and maturationJournal of Experimental Medicine 204:331–343.https://doi.org/10.1084/jem.20061442
-
Developmental control of macrophage functionCurrent Opinion in Immunology 50:64–74.https://doi.org/10.1016/j.coi.2017.12.001
-
Notch-RBP-J signaling controls the homeostasis of CD8- dendritic cells in the spleenJournal of Experimental Medicine 204:1653–1664.https://doi.org/10.1084/jem.20062648
-
Conditional gene targeting in macrophages and granulocytes using LysMcre miceTransgenic Research 8:265–277.https://doi.org/10.1023/a:1008942828960
-
An alternative pathway of imiquimod-induced psoriasis-like skin inflammation in the absence of interleukin-17 receptor a signalingJournal of Investigative Dermatology 133:441–451.https://doi.org/10.1038/jid.2012.318
-
Regulation of monocyte cell fate by blood vessels mediated by notch signallingNature Communications 7:12597.https://doi.org/10.1038/ncomms12597
-
Origin of monocytes and macrophages in a committed progenitorNature Immunology 14:821–830.https://doi.org/10.1038/ni.2638
-
Tissue biology perspective on macrophagesNature Immunology 17:9–17.https://doi.org/10.1038/ni.3320
-
The fate and lifespan of human monocyte subsets in steady state and systemic inflammationJournal of Experimental Medicine 214:1913–1923.https://doi.org/10.1084/jem.20170355
-
Regulation of innate and adaptive immunity by notchNature Reviews Immunology 13:427–437.https://doi.org/10.1038/nri3445
-
Analyzing real-time PCR data by the comparative C(T) methodNature Protocols 3:1101–1108.https://doi.org/10.1038/nprot.2008.73
-
Imiquimod-induced psoriasis-like skin inflammation in mice is mediated via the IL-23/IL-17 AxisThe Journal of Immunology 182:5836–5845.https://doi.org/10.4049/jimmunol.0802999
-
Monocytes give rise to mucosal, but not splenic, conventional dendritic cellsJournal of Experimental Medicine 204:171–180.https://doi.org/10.1084/jem.20061011
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (GA 2443/2-1)
- Jaba Gamrekelashvili
Deutsche Forschungsgemeinschaft (Li948-7/1)
- Florian P Limbourg
Deutsche Stiftung für Herzforschung (F/17/16)
- Jaba Gamrekelashvili
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the Central Animal Facility, Research Core Facility Cell Sorting and Research Core Unit Genomics of Hannover Medical School for support. We thank Thien-Phong Vu Manh and Lionel Spinelli for consultation on BubbleGUM. Funded by grants from DFG (GA 2443/2–1), DSHF (F/17/16) to JG, and DFG (Li948-7/1) to FPL.
Ethics
Animal experimentation: All experiments were performed with 8-12 weeks old mice and age and sex matched littermate controls with approval of the local animal welfare board LAVES (Niedersächsisches Landesamt für Verbraucherschutz und Lebensmittelsicherheit), Lower Saxony, Animal Studies Committee, animal study proposals #14-1666, #16-2251, #18-2777, #2014-63, #2018-221). Mice were housed in the central animal facility of Hannover Medical School (ZTL) and were maintained and supervised as approved by the Institutional Animal Welfare Officer (Tierschutzbeauftragter).
Copyright
© 2020, Gamrekelashvili et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,527
- views
-
- 447
- downloads
-
- 58
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Immunology and Inflammation
S100A8/A9 is an endogenous alarmin secreted by myeloid cells during many acute and chronic inflammatory disorders. Despite increasing evidence of the proinflammatory effects of extracellular S100A8/A9, little is known about its intracellular function. Here, we show that cytosolic S100A8/A9 is indispensable for neutrophil post-arrest modifications during outside-in signaling under flow conditions in vitro and neutrophil recruitment in vivo, independent of its extracellular functions. Mechanistically, genetic deletion of S100A9 in mice caused dysregulated Ca2+ signatures in activated neutrophils resulting in reduced Ca2+ availability at the formed LFA-1/F-actin clusters with defective β2 integrin outside-in signaling during post-arrest modifications. Consequently, we observed impaired cytoskeletal rearrangement, cell polarization, and spreading, as well as cell protrusion formation in S100a9-/- compared to wildtype (WT) neutrophils, making S100a9-/- cells more susceptible to detach under flow, thereby preventing efficient neutrophil recruitment and extravasation into inflamed tissue.
-
- Computational and Systems Biology
- Immunology and Inflammation
Diverse antibody repertoires spanning multiple lymphoid organs (i.e., bone marrow, spleen, lymph nodes) form the foundation of protective humoral immunity. Changes in their composition across lymphoid organs are a consequence of B-cell selection and migration events leading to a highly dynamic and unique physiological landscape of antibody repertoires upon antigenic challenge (e.g., vaccination). However, to what extent B cells encoding identical or similar antibody sequences (clones) are distributed across multiple lymphoid organs and how this is shaped by the strength of a humoral response remains largely unexplored. Here, we performed an in-depth systems analysis of antibody repertoires across multiple distinct lymphoid organs of immunized mice and discovered that organ-specific antibody repertoire features (i.e., germline V-gene usage and clonal expansion profiles) equilibrated upon a strong humoral response (multiple immunizations and high serum titers). This resulted in a surprisingly high degree of repertoire consolidation, characterized by highly connected and overlapping B-cell clones across multiple lymphoid organs. Finally, we revealed distinct physiological axes indicating clonal migrations and showed that antibody repertoire consolidation directly correlated with antigen specificity. Our study uncovered how a strong humoral response resulted in a more uniform but redundant physiological landscape of antibody repertoires, indicating that increases in antibody serum titers were a result of synergistic contributions from antigen-specific B-cell clones distributed across multiple lymphoid organs. Our findings provide valuable insights for the assessment and design of vaccine strategies.