- Open Access
An overlooked subset of Cx3cr1wt/wt microglia in the Cx3cr1CreER-Eyfp/wt mouse has a repopulation advantage over Cx3cr1CreER-Eyfp/wt microglia following microglial depletion
Journal of Neuroinflammation volume 19, Article number: 20 (2022)
Fluorescent reporter labeling and promoter-driven Cre-recombinant technologies have facilitated cellular investigations of physiological and pathological processes, including the widespread use of the Cx3cr1CreER-Eyfp/wt mouse strain for studies of microglia.
Immunohistochemistry, Flow Cytometry, RNA sequencing and whole-genome sequencing were used to identify the subpopulation of microglia in Cx3cr1CreER-Eyfp/wt mouse brains. Genetically mediated microglia depletion using Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice and CSF1 receptor inhibitor PLX3397 were used to deplete microglia. Primary microglia proliferation and migration assay were used for in vitro studies.
We unexpectedly identified a subpopulation of microglia devoid of genetic modification, exhibiting higher Cx3cr1 and CX3CR1 expression than Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia in Cx3cr1CreER-Eyfp/wt mouse brains, thus termed Cx3cr1highCre−Eyfp− microglia. This subpopulation constituted less than 1% of all microglia under homeostatic conditions, but after Cre-driven DTA-mediated microglial depletion, Cx3cr1highCre−Eyfp− microglia escaped depletion and proliferated extensively, eventually occupying one-third of the total microglial pool. We further demonstrated that the Cx3cr1highCre−Eyfp− microglia had lost their genetic heterozygosity and become homozygous for wild-type Cx3cr1. Therefore, Cx3cr1highCre−Eyfp− microglia are Cx3cr1wt/wtCre−Eyfp−. Finally, we demonstrated that CX3CL1–CX3CR1 signaling regulates microglial repopulation both in vivo and in vitro.
Our results raise a cautionary note regarding the use of Cx3cr1CreER-Eyfp/wt mouse strains, particularly when interpreting the results of fate mapping, and microglial depletion and repopulation studies.
Microglia are derived from the yolk sac during early embryonic development and represent approximately 10% of the healthy adult brain's total cell population [1, 2]. They play critical roles in maintaining brain development and function [3, 4]. In their homeostatic state, they display a highly ramified morphology, efficiently surveying the central nervous system (CNS) microenvironment, recognizing and clearing cell debris . Microglia play pivotal roles in diverse CNS diseases including neurodevelopmental disorders, neurodegenerative disorders, and high-grade glioma [6,7,8,9,10,11,12,13]. Targeting microglia has thus emerged as an attractive strategy to modulate neuroinflammation in the context of various CNS diseases [14, 15].
The Cx3cr1 gene is mainly expressed in microglia in the CNS parenchyma, and the CX3CL1/CX3CR1 axis plays a crucial role in microglia–neuron communication [16, 17]. Cx3cr1 is widely used to genetically label microglia by inserting Gfp [18, 19] since this enables microglia tracing, visualization and sorting [1, 20,21,22,23]. Cx3cr1Cre/wt and Cx3cr1CreER-Eyfp/wt mouse strains are potent, commonly used tools for studies of microglial fate mapping [1, 24], microglial depletion [25,26,27,28] and modification of the microglial genome by leveraging floxed target genes [29,30,31]. These methods have greatly increased our understanding of microglia in CNS homeostasis and disease conditions.
In this study, we report the existence of an unexpected, small population of Cx3cr1highCre−Eyfp− microglia in Cx3cr1CreER-Eyfp/wt mice. After genetically mediated microglial depletion using Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice, Cx3cr1highCre−Eyfp− microglia escape depletion, display a repopulation advantage and eventually constitute one-third of the total repopulated microglial pool. We further determined that the Cx3cr1highCre−Eyfp− microglia are Cx3cr1wt/wtCre−Eyfp−. Finally, we demonstrate the vital role of CX3CL1–CX3CR1 signaling in regulating microglial repopulation post-depletion. Not being aware of this population may result in misinterpretation of the results generated since these cells escape detection (not carrying the Eyfp or Gfp) and cannot be modified (lacking Cre expression) as expected.
Breeding pairs of Rosa26DTA/DTA and Cx3cr1CreER-Eyfp/CreER-Eyfp mice were bought from the Jackson Laboratory (Bar Harbor, ME, USA) with stock numbers 009669 and 021160, respectively. Second generation pups with genotypes Cx3cr1CreER-Eyfp/wtRosa26wt/wt, Cx3cr1CreER-Eyfp/wtRosa26DTA/wt, Cx3cr1CreER-Eyfp/CreER-EyfpRosa26wt/wt and Cx3cr1CreER-Eyfp/CreER-EyfpRosa26DTA/wt were used for all experiments. Cx3cr1Gfp/Gfp mice breeding pairs were bought from the Jackson Laboratory (005582, Bar Harbor, ME, USA), C57BL/6 mice were bought from Charles River (Sulzfeld, Germany), and the pups of Cx3cr1Gfp/Gfp and C57BL/6 mice were used for the experiments. All mice were housed with a 12:12 h light–dark cycle and had free access to food and water.
Tamoxifen (T5648-1G, Sigma-Aldrich, Merck, Germany) was dissolved in corn oil (C8267, Sigma-Aldrich, Merck, Germany), and 125 mg/kg or 62.5 mg/kg was injected i.p. every 24 h for 3 (postnatal days 18, 19 and 20) or 10 consecutive days (postnatal days 18–27).
PLX3397 chow treatment
PLX3397 (HY-16749, MedChemExpress) was formulated with standard chow at 290 mg/kg by SAFE Nutrition Service, France. Four-week-old mice were treated with PLX3397 chow for 3 weeks and sacrificed at different timepoints afterwards.
DNA was extracted from ear punches and microglia. The following primer sets were used for the genotyping, Cx3cr1CreER-Eyfp: common 5′-AAG ACT CAC GTG GAC CTG CT-3′, Wild-type reverse 5′-AGG ATG TTG ACT TCC GAG TTG-3′, Mutant reverse 5′-CGG TTA TTC AAC TTG CAC CA-3′, Long targeting primer 5′-TGTCCGTATAGGTTGGGAAA-3′; DTA: common 5′-AAA GTC GCT CTG AGT TGT TAT-3′, Wild-type reverse 5′-GGA GCG GGA GAA ATG GAT ATG-3′, Mutant reverse 5′-GCG AAG AGT TTG TCC TCA ACC-3′.
Immunohistochemistry, cell counting, and image analysis
After perfusion with cold PBS, the brain was dissected, immersed in 4% paraformaldehyde (PFA) for 2 days, followed by immersion in 30% sucrose in 0.1 M phosphate buffer until fully saturated. Brain tissues were sagittally sliced into 25 µm sections and stored in tissue cryoprotectant solution (25% ethylene glycol and 25% glycerin in 0.1 M phosphate buffer). The sections were rinsed in Tris-buffered saline (TBS, 50 mM Tris–HCl in 150 mM NaCl, pH 7.5) and then blocked with 3% donkey serum in TBS with 0.1% Triton X-100. After blocking, the sections were incubated overnight in primary antibody solution (Goat anti-Iba-1, 1:500, ab5076, Abcam; Chicken anti-GFP, 1:800, ab13970 Abcam; Rabbit anti-Tmem119, 1:500, ab209064, Abcam; Rabbit anti-P2ry12, 1:500, generated at Harvard Medical School; Rat anti-Ki67, 1:500, 14-5698-82, eBioscience). After washing, the sections were incubated for 2 h at room temperature with a secondary antibody (Donkey anti-goat, 1:1,000, Alexa Fluor 633, 20127, BIOTIUM; Donkey anti-chicken, 1:1000, Alexa Fluor 488, 703-545-155, Jakson Immuno Research; Donkey anti-rabbit, 1:1000, Alexa Fluor 555, 1945911, Invitrogen; Donkey anti-Rat, 1:1,000, Alexa Fluor 555, ab150154, Abcam) together with Hoechst (1:1000, 33342 Life Technologies). The sections were mounted onto glass slides with ProLong® Gold antifade reagent (1925239, Invitrogen). For Iba-1+ counting, the sections were blocked with 3% donkey serum in TBS with 0.1% Triton X-100, and then incubated in Rabbit anti-Iba-1 1:1,000 (Wako) overnight. After that, sections were incubated with the secondary antibody (donkey anti-rabbit 1:500, 20215, BIOTIUM) for 2 h, and then endogenous peroxidase activity was blocked with 3% H2O2 for 10 min. Visualization was performed using Vectastain ABC Elite with 3,3-diaminobenzidine.
Iba-1+ microglia were counted using an Axio Imager M2 microscope with Apotome attachment (Carl Zeiss). Cells were counted exhaustively on day 1, day 3, day 7 and day 10 after the final Tam injection. While a fractionator was used to estimate the positive cells on day 42 after the final Tam injection, the counting frame width and height are 150 μm, sample grid X and Y are 300 μm. Thus, all the data has a 2nd estimated coefficient of errors (Schmitz-Hof) less than 0.1 . Then, relative microglia number ratios were calculated with reference to the control group.
All the fluorescence images were captured using an LSM700 laser scanning confocal microscope (Axio-observer Z1; CarlZeiss microscopy, Germany) and analyzed using ZEN software (the black edition; Zeiss).
Confocal microscopy was used to acquire images at 3 µm intervals, with a 40× objective lens used for Cx3cr1 depletion analysis (Fig. 2g and Additional file 1: Fig. S4G) and a 20× lens used for Cx3cr1CreER-Eyfp/wt physiological analysis (Additional file 1: Fig. S5A) (Plan-Apochromat lens, 20×/0.8 and 40×/1.3 Oil objective, Carl Zeiss). Microglial morphology was analyzed with the skeleton analysis method using FIJI open-source image analysis software (Image J 1.51s, NIH) [33, 34]. Briefly, the maximum intensity projections of z-series stacks were created. The Iba-1+ channel and the CX3CR1-EYFP channel were enhanced in brightness to visualize all processes, followed by de-speckling to eliminate single-pixel backgrounds. Then, each channel image was processed into a binary image, and a topological skeleton image was created (Additional file 1: Fig. S2A). The Analyze Skeleton plugin feature was applied to the whole image frame, and the process length and the number of endpoints were normalized by the cell number per frame. Three morphological parameters were used: (1) total length of the processes, (2) number of endpoints, and (3) microglial process area. The process length and the number of endpoints of CX3CR1+/+EYFP−Iba-1+ cells were calculated by subtracting the CX3CR1-EYFP channel data from the Iba-1 channel data. Each cell's process area was represented as the convex hull area by connecting the process ends using the polygon tool (Additional file 1: Fig. S2B, C) . The selection criteria for the process area analysis were relatively isolated from the processes of the surrounding cells, and whole processes were not truncated and are within the image.
Single-cell suspensions were prepared using a Neural Tissue Dissociation Kit (T, Miltenyi Biotec), and myelin was removed using 38% Percoll. After passing through a 40 μm cell strainer, cells were incubated with the CD11b magnetic bead for 20 min. CD11b+ cells were isolated and collected using a magnetic field and MS column (Miltenyi Biotec, Germany).
After microglia isolation, cells were transferred into 96-well V-bottom plates and incubated at 4 °C for 20 min with the following antibodies: yellow dead cell stain kit (L34959, Thermo Fisher Scientific), Near-IR dead cell stain kit (L34975, Thermo Fisher Scientific), CD11b (PerCP-Cyanine5.5, Clone M1/70, BioLegend), CD45 (PE-Cyanine7, Clone 30F11, BioLegend), Ly6C (PE, Clone HK1.4, BioLegend), Ly-6G (V450, Clone 1A8, BD Biosciences), F4/80 (APC, Clone BM8, BioLegend), MHCII (Alexa Fluor 700, Clone M5/114.15.2, BioLegend), Cx3cr1 (Alexa Fluor 700, Clone SA011F11, Biolegend), Cre Recombinase (PE, Clone D7L7L, Cell Signaling Technology) and Ki67 (V450, Clone B56, BD). Cells were run on a Gallios flow cytometer (Beckman Coulter), and the data were analyzed with Kaluza software (Beckman Coulter).
Cell sorting and bulk RNA sequencing
Cells were sorted into CD11b+CD45+Ly6C−CX3CR1+EYFP+ and CD11b+CD45+Ly6C−CX3CR1+EYFP− populations using a BD Influx Cell Sorter. Briefly, 100 cells were collected directly into lysis buffer, and the library was built by the smart-seq method. Gene expression data were analyzed using Qlucore Omics Explorer 3.4. For two-group and multi-group comparisons of candidate genes, expression was considered as significantly different if p < 0.05 using a heteroscedastic two-tailed Student’s t test or an F test, respectively. For transcriptome-wide analyses, the Benjamini–Hochberg method was used to correct multiple tests.
Whole-genome sequencing was performed essentially as in DNTR-seq  with minor modifications. Cell nuclei (1,500–5,000 nuclei per sample) were treated with 15 μl proteinase K (0.4 ml/ml) for 2 h at 50 °C followed by heat inactivation at 80 °C for 30 min. 2 μl of the solution was taken for tagmentation with 1 μl Tn5 stock solution and 1.6 μl reaction buffer (50 mM TAPS, 25 mM MgCl2, 40% PEG 8k; final concentration 8% PEG, 5 mM MgCl2, 10 mM TAPS) and 3.4 ml H2O (final volume 8 μl), incubated at 55 °C for 10 min, and then inactivated by adding 2 μl 0.2% SDS and incubating at 55 °C for 10 min. For barcoding PCR, each sample was split into 4 individual reaction, each barcoded with a unique molecular barcode (fastq files for the replicate reactions were merged after sequencing). 1 μl of the tagmentated DNA solution was used as the template for each barcoding PCR reaction. The PCR program was 72 °C/3 min, 95 °C/30 s, [95 °C/15 s, 67 °C/30 s, 72 °C/1 min] × 15 cycles, 72 °C/5 min, and then 10 °C hold. After PCR amplification, the reactions were cleaned-up by SPRI-beads (at 0.9× volume) for size selection and were pooled together according in equal parts according to their DNA concentration (qubit dsDNA quantification assays, ThermoFisher). The pooled library was then cleaned-up one more time by SPRI-beads (at 0.9× volume) and sequenced on an Illumina NextSeq 550 using the high output 75 bp kit.
Genomic reads were trimmed for adaptor sequences using TrimGalore and mapped to human genome reference build hg38 with the non-genomic parts of the Cx3cr1-CreERT2 targeting vector added as a separate template. Duplicates were removed using Picard MarkDuplicate, along with read pairs with MAPQ < 20. Reads mapping to the part of the vector which is inserted into the genome (position 1173–6936), and to the region surrounding the insertion site (chr9:120,040,000–120,080,000) were extracted from the bam files. Physical coverage (eg. the number of times a base is spanned by paired read mapping positions) was visualized using a custom R script.
Primary microglia proliferation and migration assay
Primary microglia were cultured as follows: cerebrum tissues were collected from postnatal days 3 mice and used for generating mixed glia cultures in T75 flasks. After 2 weeks of primary culture (DMEM F-12 with 20 ng/ml M-CSF), the microglia were isolated with CD11b magnetic bead sorting (MACS, Miltenyi Biotec) according to the manufacturer’s protocol. For proliferation assay, isolated CD11b positive cells were seeded at a density of 0.2 × 105 on the coverslips (83.1840.002, Sarstedt) in 24-well plate for 2 days with DMEM F-12 medium containing 10% FCS and 20 ng/ml M-CSF, followed by replacement with serum-free medium (DMEM F-12) and serum-free medium (DMEM F-12) with 100 ng/ml CX3CL1 (472-FF-025/CF, R&D systems). 24 h later the cells were fixed by 4% paraformaldehyde, then incubated overnight in primary antibody solution (Goat anti-Iba-1, 1:500, ab5076, Abcam; Rat anti-Ki67, 1:500, 14-5698-82, eBioscience). After washing, the coverslips were incubated for 2 h at room temperature with a secondary antibody (Donkey anti-goat, 1:1000, Alexa Fluor 488, 1942238, invitrogen; Donkey anti-Rat, 1:1000, Alexa Fluor 555, ab150154, Abcam). The rate of ki67 positive cells in Iba-1 positive cells were calculated as the proliferating rate of microglia. For migration assay, isolated CD11b positive cells were seeded at a density of 0.5 × 105 in an Incucyte ImageLock 96-well plate (4379, Sartorius) for overnight with the addition of 10% FCS and 20 ng/ml M-CSF. On the following day, a standardized scratched wound was induced in each well simultaneously using an Incucyte WoundMaker kit (Cat. No. 4493, Sartorius), followed by replacement with serum-free, M-CSF free DMEM/F-12 medium containing indicated concentration of CX3CL1 (472-FF-025/CF, R&D systems). The plate was then placed in an Incucyte Zoom System, and live images were taken every 2 h. The images were analyzed according to the manual of Incucyte Zoom.
Identification of a subset of Cx3cr1 high Cre − Eyfp − parenchymal microglia in Cx3cr1 CreER-Eyfp/wt mice
In Cx3cr1CreER-Eyfp/wt mice, the protein-coding exon of one of the Cx3cr1 alleles is replaced by the Cre-ERT2 fusion gene and the enhanced yellow fluorescent protein gene (Eyfp), implying that all CX3CR1+ cells (mainly microglia in the CNS parenchyma) should express both Cre-ERT2 and EYFP. Unexpectedly, we observed a subpopulation of Iba-1+Tmem119+EYFP− cells with microglia-like morphology in the CNS parenchyma of Cx3cr1CreER-Eyfp/wt mice under homeostatic conditions (Fig. 1a and Additional file 1: Fig. S1a). These Iba-1+Tmem119+EYFP− cells and Iba-1+Tmem119+EYFP+ microglia exhibited similar morphologies as measured by process complexity (numbers of endpoints, length) and area occupied (Fig. 1b and Additional file 1: Fig. S2). Flow cytometric analysis (gating strategy depicted in Additional file 1: Fig. S3a) confirmed the existence of EYFP− cells in Cx3cr1CreER-Eyfp/wt mice, accounting for less than 1% of total CD11b+CD45+ Ly6C–Ly6G– cells at 3–15 weeks, with no significant differences between ages. One outlier sample in 9-week age mice was excluded from the analysis, having 6% of EYFP− cells of the total CD11b+CD45+ Ly6C−Ly6G− pool (Fig. 1c and Additional file 1: Fig. S3b). Furthermore, we observed that both EYFP+ microglia and EYFP− microglia-like cells were CX3CR1+ (Fig. 1d). These data indicate that CX3CR1+EYFP− cells are resident microglia.
To further confirm this finding, we performed fluorescence-activated cell (FAC) sorting of CX3CR1+EYFP+ and CX3CR1+EYFP− cells (gated as CD11b+CD45+Ly6C−Ly6G−CX3CR1+) followed by RNA sequencing (Fig. 1d). To investigate the Cre and Eyfp expression, the RNA sequencing reads were mapped on the mouse genome integrated with the CreERT2-Eyfp knockin-allele. We determined that both CX3CR1+EYFP+ and CX3CR1+EYFP− cells expressed similar levels of specific markers for myeloid cells (Hexb, CD68, Fcgr1, Adgre1, Itgam, Sparc, C1qa, Aif1, and Fcer1g) and microglia (P2ry12, Tmem119, and Sall1) (data not included), further confirming that both CX3CR1+EYFP− and CX3CR1+EYFP+ cells were CNS parenchymal microglia. Next, the transcriptomes of CX3CR1+EYFP− and CX3CR1+EYFP+ microglia were compared, indicating 32 differently expressed genes (DEGs) with a false discovery rate (FDR) < 0.05 (Additional file 2: Table S1). The top 5 downregulated genes in CX3CR1+EYFP+ microglia were Wfdc21, Camp, Chil3, Ltf, and Cd177, while all the upregulated genes in CX3CR1+EYFP+ microglia were Eyfp, Cre, H2-Aa, Ifi47, Zfp940, Ms4a4a, and Itsn1 (Fig. 1e). Furthermore, we determined that CX3CR1+EYFP− microglia expressed twice as much Cx3cr1 mRNA compared with CX3CR1+EYFP+ microglia (Fig. 1f). Moreover, CX3CR1+EYFP− microglia did not express Eyfp or Cre mRNA (Fig. 1g, h). We thus defined that CX3CR1+EYFP− microglia are Cx3cr1highCre−Eyfp−.
To assess the variation of Cx3cr1 expression in C57BL/6 mice, single-cell RNA sequencing data from our previous publication  were reanalyzed in CX3CR1+ microglia from the hippocampus. No clusters with high or low Cx3cr1 expression were detected (Additional file 1: Fig. S3c). The Cx3cr1 expression was normally distributed in CX3CR1+ microglia (Additional file 1: Fig. S3d), supporting that the difference in Cx3cr1 expression levels between Cx3cr1highCre−Eyfp− and Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia in Cx3cr1CreER-Eyfp/wt mice was attributed to the genetic modification.
To confirm that Cx3cr1highCre−Eyfp−microglia were not an anomaly of one specific mouse strain, we investigated the microglial populations in Cx3cr1GFP/wt mice (JAX, 005582) by flow cytometry (gating strategy depicted in Additional file 1: Fig. S3a) and immunohistochemistry. The results revealed that Iba-1+Tmem119+GFP− cells with microglial morphology could also be detected in the Cx3cr1GFP/wt adult mouse brain, constituting 0.59% ± 0.27% of total microglia (Fig. 1i, j).
Cx3cr1 high Cre − Eyfp − microglia repopulated and became a dominant subgroup of the repopulated microglia pool following microglial depletion
To investigate whether Cx3cr1highCre−Eyfp− microglia are capable of filling the vacant microglial niche following microglial depletion, we first crossed Cx3cr1CreER-Eyfp/CreER-Eyfp and Rosa26DTA/DTA mouse strains (Fig. 2a, JAX, 009669) to breed Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice. Following tamoxifen (Tam) injections, Cx3cr1CreER-Eyfp/wt cells in the brain parenchyma (mainly Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia) can be depleted through the intracellular release of diphtheria toxin . To track dynamic microglial alterations in the brain, Cx3cr1CreER-Eyfp/wt Rosa26DTA/wt and Cx3cr1CreER-Eyfp/wt mice were sacrificed at different timepoints (days 1, 3, 7, 8, 10, 21, and 42) after three consecutive Tam injections on postnatal days 18, 19, and 20, respectively (Fig. 2b). Iba-1+ microglia were quantified by immunohistochemistry, revealing that approximately 89.7 ± 2.4%, 90.9 ± 0.3%, and 94.5 ± 1.4% of the microglia were depleted by the first day after the final Tam injection in the hippocampus, cortex, and cerebellum, respectively, (Fig. 2c, d, Additional file 1: Fig. S4a, b). Moreover, newly repopulated microglia were noted 1 week later, amounting to 21.9 ± 8.5% in the cortex and 34.9 ± 5.9% in the cerebellum compared to control Cx3cr1CreER-Eyfp/wt mice (Additional file 1: Fig. S4a, b). We also recorded an overabundance of microglia 10 days after the final Tam injection, 123.3 ± 2.7% in the hippocampus and 139.8 ± 18.6% in the cerebellum compared to the control group (Fig. 2d and Additional file 1: Fig. S4a). No significant differences were observed in the numbers of newly repopulated microglia 42 days after the final Tam injection, compared with the baseline level (Fig. 2c, d, Additional file 1: Fig. S4a, b).
The numbers of Cx3cr1highCre−Eyfp− (EYFP−Iba1+) microglia increased from day 1 to day 3 after the final Tam injection (Additional file 1: Fig. S4c, d). Distinct clusters of Cx3cr1highCre−Eyfp− microglia were formed by day 7, further expanded on days 8 and 10, and then maintained 42 days after the final Tam injection (Fig. 2e, f). Interestingly, newly repopulated Cx3cr1highCre−Eyfp− microglial clusters were distributed in a stochastic manner with a considerable regional variation between different mice, not co-localizing with Cx3cr1CreER-Eyfp/wt Cre+Eyfp+ microglia clusters (Fig. 2e, f). We have previously demonstrated that high and low expression of F4/80 can be used to distinguish resident repopulated microglia and peripheral infiltrating microglia-like cells following experimental microglial depletion . Flow cytometric analysis revealed 3 main subsets of microglia-like cells after repopulation in Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice (Fig. 2g). These included EYFP+F4/80low, EYFP−F4/80low, and EYFP+F4/80hi, likely representing resident Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia, resident Cx3cr1highCre−Eyfp− microglia, and peripherally derived microglia-like cells, respectively (Fig. 2g). The finding was further confirmed using Immunohistochemistry (IHC), again demonstrating mainly 3 repopulated microglial subgroups expressing EYFP+Tmem119+, EYFP−Tmem119+, and EYFP+Tmem119−, respectively (Fig. 2h).
We next addressed if these newly repopulated Cx3cr1highCre−Eyfp− microglia were derived from the periphery. We and others have previously demonstrated that the empty microglial niche can be repopulated within weeks through resident microglia proliferation and concomitant infiltration of monocytes [38, 39]. Two novel specific surface markers (Tmem119 and P2ry12) were used to identify CNS-resident microglia [40, 41]. Our results demonstrated that all newly repopulated Cx3cr1highCre−Eyfp− microglia expressed both Tmem119 (Fig. 2h) and P2ry12 (Fig. 3a). Furthermore, the vast majority of Cx3cr1highCre−Eyfp− microglia had low expression of F4/80 (Fig. 2g).
To confirm the CNS origin of Cx3cr1highCre−Eyfp− microglia, we reanalyzed data from our previous publication . Specifically, Cx3cr1CreER-Eyfp/wtRosa26DTA/wt and Cx3cr1CreER-Eyfp/wt mice (CD45.2) were exposed to head-protected irradiation and transplanted with congenic CD45.1 bone marrow, Tam was injected 8 weeks after reconstitution, and the chimeras were analyzed another 8 weeks later (Fig. 3b). We determined that the vast majority of the F4/80high cells were derived from CD45.1+ donor cells, while all Cx3cr1CreER-Eyfp/wtCre+Eyfp+ (EYFP+F4/80low) microglia originated from CD45.2+ residential microglia in both Cx3cr1CreER-Eyfp/wt and Cx3cr1CreER-Eyfp/wtRosa26DTA/wt chimeras, and that the vast majority of the Cx3cr1highCre−Eyfp− (EYFP−F4/80low) microglia were derived from CD45.2+ residential microglia, indicating their CNS origin.
This notion was further confirmed by RNA sequencing of the 4 groups of repopulated microglia expressing EYFP+F4/80low, EYFP−F4/80low, EYFP+F4/80high, and EYFP−F4/80high, respectively (Fig. 2g). Specifically, our results demonstrated that the F4/80low groups had lower expression of signature genes of engrafted macrophages, including Ms4a7 and Lyz2, compared to the F4/80high groups (Fig. 3c, d). Principle component analysis (PCA) indicated that the greatest difference was evident between F4/80low and F4/80high groups (Fig. 3e). We further investigated the expression levels of 14 macrophage- and microglia-related genes , and EYFP+F4/80low and EYFP−F4/80low microglia expressed similar patterns, while EYFP+F4/80high and EYFP−F4/80high microglia expressed similar patterns (Fig. 3f). For example, the F4/80high groups expressed higher levels of Adgre1 (F4/80) and Ptprc (CD45) than did the F4/80low groups (Fig. 3f), supporting the efficiency of our sorting strategy and the peripheral source of the F4/80high groups. F4/80high and F4/80low groups expressed similar levels of Fcer1g, Aif1 (Iba-1), and Cd68. Moreover, F4/80high groups expressed lower levels of Fcgr1, C1qa, Cx3cr1, Sall1, Tmem119, sparc1, P2ry12, Hexb and Itgam (CD11b) (Fig. 3f). Taken together, newly repopulated Cx3cr1highCre−Eyfp−microglia originate from the resident Cx3cr1highCre−Eyfp− microglia pool.
Collectively, although only less than 1% Cx3cr1highCre−Eyfp− microglia exist in the Cx3cr1CreER-Eyfp/wt mouse brain, they can account for one-third of the total newly repopulated microglia pool following conditional genetic depletion in Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice.
Cx3cr1 high Cre − Eyfp − microglia have lost the CreERT2-Eyfp fusion gene and reveal a wild-type Cx3cr1 genotype
We next investigated the newly repopulated CNS-derived (F4/80low) microglia. Consistent with our findings in naive mice, EYFP− microglia sorted after repopulation lacked Cre and Eyfp expression but expressed twice the levels of Cx3cr1 compared to EYFP+ microglia (Additional file 1: Fig. S5a–c). Moreover, flow cytometric analysis indicated a higher expression of CX3CR1 protein in Cx3cr1highCre−Eyfp− (EYFP−) microglia compared to Cx3cr1CreER-Eyfp/wtCre+Eyfp+ (EYFP+) microglia (Additional file 1: Fig. S5d, e).
Given the lack of Cre and Eyfp expression on the mRNA and protein levels, we reasoned that the CreER-Eyfp knockin-allele was either epigenetically/genetically silenced or deleted. We thus investigated Cx3cr1highCre−Eyfp− microglia at the DNA level. Genotyping of DNA extracted from Cx3cr1highCre−Eyfp− and Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia revealed that Cx3cr1highCre−Eyfp− microglia apparently lacked the CreERT2-Eyfp insert (Fig. 4a). This analysis was performed using the primers and protocols Jax lab recommended, which only target part of the CreERT2-Eyfp knockin-allele (Additional file 1: Fig. S6a). We, therefore, further designed primers and optimized protocols to target longer sequences, including the whole Cx3cr1 exon 2 containing the inserted CreERT2-Eyfp fusion gene (Additional file 1: Fig. S6a). The mutant band was not detected in Cx3cr1highCre−Eyfp− microglia (Fig. 4b), again indicating that they lack the CreERT2-Eyfp insert. Whole-genome sequencing of the DNA extracted from Cx3cr1highCre−Eyfp− and Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia further confirmed the complete absence of the CreERT2-Eyfp knockin-allele in Cx3cr1highCre−Eyfp− microglia. Interestingly, the coverage of the integrated Cx3cr1 coding exon in Cx3cr1highCre−Eyfp− microglia is twice that of the coverage in the Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia (Fig. 4c–e, Additional file 1: Fig. S6b, e), indicating that Cx3cr1highCre−Eyfp− microglia are Cx3cr1wt/wtCre−Eyfp−.
We next addressed the origin of Cx3cr1wt/wtCre−Eyfp− microglia. The experimental Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice were bred from male Cx3cr1CreER-Eyfp/CreER-EyfpRosa26DTA/DTA and female Cx3cr1wt/wtpRosa26wt/wt mice. We first excluded the possibility of maternally derived microglia since the Cx3cr1wt/wtCre−Eyfp− microglia inherited the paternal DTA knockin-allele (Additional file 1: Fig. S6d). We then hypothesized that homologous recombination during microglial mitosis leads to loss of the CreERT2-Eyfp insert. To address this hypothesis, we tested whether breeding Cx3cr1CreER-Eyfp/CreER-EyfpRosa26DTA/wt mice would result in the loss of the CreERT2-Eyfp insert, given the fact that there is no wild type allele to recombine with. Our results revealed that the EYFP−F4/80low microglia, here denoted as Cx3cr1wt/wtCre−Eyfp− microglia, were not detected in Cx3cr1CreER-Eyfp/CreER-EyfpRosa26DTA/wt mice after genetically microglial depletion and repopulation (Additional file 1: Fig. S6e). These results indicate that homologous recombination is a possible mechanism leading to loss of heterozygosity (LOH) of Cx3cr1wt/wtCre−Eyfp− microglia in Cx3cr1CreER-Eyfp/wt mice.
CX3CL1–CX3CR1 signaling regulates resident microglia repopulation
We investigated whether the Cx3cr1wt/wtCre−Eyfp− microglia have a competitive advantage during repopulation following depletion due to their higher expression of Cx3cr1, analogous to what was reported in a study of the retina . Our results revealed that after depletion, the Cx3cr1wt/wtCre−Eyfp− (EYFP−F4/80low) microglia indeed had an advantage over Cx3cr1CreER-Eyfp/wtCre+Eyfp+ (EYFP+F4/80low) microglia during repopulation. Their numbers increased more quickly, reaching similar numbers as the Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia despite originating from a smaller number of surviving cells (Fig. 5a). To investigate whether this repopulation advantage is due to higher Cx3cr1 expression, their proliferative capacity was assayed in Cx3cr1wt/wt, Cx3cr1CreER-Eyfp/wt, and Cx3cr1CreER-Eyfp/CreER-Eyfp primary culture microglia. We observed that microglial proliferation rates decreased dose-dependently with decreasing CX3CR1 expression, but independently of the presence of the ligand, CX3CL1 (Fig. 5b, c). Furthermore, a scratch assay revealed that microglial migration was decreased in Cx3cr1CreER-Eyfp/CreER-Eyfp microglia compared to Cx3cr1wt/wt and Cx3cr1CreER-Eyfp/wt microglia. After adding exogenous CX3CL1, microglial migration rates increased in Cx3cr1wt/wt and Cx3cr1CreER-Eyfp/wt microglia, more so in Cx3cr1wt/wt than Cx3cr1CreER-Eyfp/wt, but not in Cx3cr1CreER-Eyfp/CreER-Eyfp microglia (Fig. 5d). These data indicate that the CX3XL1/CX3CR1 axis is important in microglial repopulation following depletion.
To further confirm the role of CX3CR1 in regulating microglial repopulation. We depleted microglia using PLX3397 chow in Cx3cr1wt/wt, Cx3cr1CreER-Eyfp/wt, and Cx3cr1CreER-Eyfp/CreER-Eyfp mice. Previous studies have demonstrated that blood-borne infiltrating cells do not contribute to microglia repopulation following microglial depletion using PLX3397 administered in chow . Approximately 70% of the microglia were depleted after 21 days PLX3397 treatment and there was no difference in depletion efficiency between the 3 genotypes. Microglia repopulated quickly 3 days after ceasing treatment, without any differences between the 3 genotypes. However, 7 days after repopulation the Cx3cr1wt/wt group exhibited more microglia than did the Cx3cr1CreER-Eyfp/CreER-Eyfp group, reaching 20% higher numbers than the control group, while the Cx3cr1CreER-Eyfp/CreER-Eyfp group stopped at the control, pre-depletion level and the Cx3cr1CreER-Eyfp/wt group displayed intermediate numbers (Fig. 5e), further indicating that CX3CR1 plays a role in repopulation when the niche is filling up.
Collectively, our results show that Cx3cr1highCre−Eyfp− microglia are Cx3cr1wt/wtCre−Eyfp−. They have an intrinsic, competitive repopulation advantage, which is at least partly due to higher Cx3cr1 expression than in Cx3cr1CreER-Eyfp/wt microglia.
The Cx3cr1CreER-Eyfp mouse line constitutes a central tool for studying microglia biology . In the present study, we unexpectedly identified the presence of Cx3cr1wt/wtCre−Eyfp− microglia lacking CreERT2-Eyfp locus in the Cx3cr1CreER-Eyfp/wt mouse brain, which can be repopulated after microglial depletion following Tam injections in Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice.
One previous study reported that all repopulated microglia were derived from CX3CR1+ cells, as measured by fate mapping in Cx3cr1CreER::Ai14 mice. All CX3CR1+ microglia were tdTomato+, and all Iba-1+ cells were tdTomato+, leading to the interpretation that all the repopulated Iba-1+ cells arose from the surviving CX3CR1+ cells following microglial depletion . However, in both our Cx3cr1CreER-Eyfp/wt and Cx3cr1GFP/wt transgenic mice, the presence of Cre−EYFP− and GFP− microglia was recorded. Proliferating GFP− microglia in Cx3cr1GFP/+ mice have also been reported following microglia depletion using the CSF1 receptor inhibitor PLX3397 . In that study, the authors showed that proliferating GFP− cells could become Iba-1+ microglia; however, whether these Iba-1+GFP− microglia were also CX3CR1+ was not addressed . Our results indicated that these cells could be de facto microglia expressing all thus far established canonical markers and that they existed before microglial depletion and repopulation. Another study reported that most newly repopulated microglia were tdTomato+ using the CSF1 receptor inhibitor PLX5622 in Cx3cr1CreER/+:tdTomato mice, and the Iba1+tdTomato− cells were interpreted as being infiltrating peripheral monocytes . However, based on our study, we consider that the Iba1+tdTomato− cells could be derived from the local microglial cell pool, which does not carry the transgene Cx3cr1CreER. Furthermore, another study demonstrated that EYFP− microglia escaped Cre-mediated recombination and could repopulate the CNS following microglial depletion in Rosa26−STOP−EyfpCx3cr1CreER:iDTR mice . Phagocytosis and activation of microglia is, at least partly, dependent on the CX3CR1/CX3CRL axis [44,45,46,47,48,49]. Hence, the EYFP–microglia, expressing twice the amount of CX3CR1, might have a different phagocytosis capacity. In light of these findings, interpretation of results from Cx3cr1CreER/wt and Cx3cr1GFP/wt transgenic mice should be made with caution.
We consider two possible explanations for the presence of Cx3cr1wt/wtCre−Eyfp− microglia in the Cx3cr1CreER-Eyfp/wt mouse brain: (i) maternal-derived macrophages or microglia, or (ii) LOH of microglia during mitosis. The possibility of maternal-derived microglia or macrophages was excluded by our results, as the Cx3cr1wt/wtCre−Eyfp− microglia inherited the paternal DTA gene. Alternatively, LOH due to homologous recombination could explain the existence of Cx3cr1wt/wtCre−Eyfp− microglia in the Cx3cr1CreER-Eyfp/wt mouse brain. LOH is a phenomenon whereby the cells only possess the genetic information from one of the parental chromosomes, as previously described in the cancer cells  and mammalian cells in vivo and in vitro [51,52,53]. Our results showing a similar percentage of Cx3cr1wt/wtCre−Eyfp− microglia in each of the heterozygous Cx3cr1CreER-Eyfp/wt mouse brains and the absence of Cx3cr1wt/wtCre−Eyfp− microglia in the homozygous Cx3cr1CreER-Eyfp/CreER-Eyfp mouse brain support the theory of LOH through homologous recombination during mitosis. The mechanisms of microglial LOH in Cx3cr1CreER-Eyfp/wt and Cx3cr1GFP/wt mouse strains and whether the phenomena of LOH also occurs in other analogous genetic modified heterozygous strains need to be further investigated.
Microglial depletion and repopulation studies have expanded our knowledge of microglia in physiological and pathological states, exerting favorable effects in different preclinical disease models [54, 55]. Moreover, new microglia rapidly repopulate the brain parenchyma following microglial depletion, although the origin of these newly repopulated microglia has been debated. Elmore and colleagues reported that CX3CR1+GFP− cells were potential microglial precursor cells during the microglial repopulation period. However, this viewpoint has been challenged, and newly repopulated microglia are proposed to only arise from the surviving microglia . This is further supported using fate mapping approaches showing that The new forming microglia only temporary expression nestin, and no microglia were derived from None microglia nestin+ cells . Our RNA sequencing data revealed no increased expression of precursor of stem cell genes in Cx3cr1wt/wtCre−Eyfp− microglia. Moreover, microglial single-cell data from wild type mice revealed no high Cx3cr1-expressing microglia that clustered separately from other microglia. Rather, the cells expressing higher levels of Cx3cr1 were evenly distributed (Additional file 1: Fig. S3C, D). Taken together, these data support that Cx3cr1wt/wtCre−Eyfp− microglia are not precursors or stem cells. How these surviving microglia could survive depletion and whether they were unaffected by the depletion period remain unclear. Microglia cannot actually be fully depleted using currently available depletion models  and the existence of a microglial subset that may be resistant to depletion has been previously proposed . Cx3cr1wt/wtCre−Eyfp− microglia could be one such subset, at least in the Cx3cr1-Cre derived depletion setting. Our results indicate that both unaffected Cx3cr1wt/wtCre−Eyfp− microglia and surviving Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia repopulate the brain competitively, with contribution from peripherally derived macrophages. Moreover, it is widely known that individual microglia occupy non-overlapping spatial territories [58, 59]. Herein, we report a novel finding that non-overlapping niches exist between repopulated Cx3cr1wt/wtCre−Eyfp− and Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia; however, related signal pathways need to be further investigated.
Tissue macrophages with a self-renewing capacity are seeded during embryonic development, and some macrophages can be replaced or renewed postnatally by peripheral monocytes [60,61,62,63,64,65]. In the brain, the resident microglial pool is self-renewing without contribution from peripheral monocytes during the whole life span . However, in disease states, such as CNS injury and neurodegenerative diseases, monocytes can infiltrate the brain and become microglia-like cells, albeit with different functionalities . Here we used a mouse model in which peripheral monocytes entered the brain after microglial depletion . Our results demonstrate that the main transcriptomic difference is noted between the F4/80high and F4/80low groups, but not between the EYFP− and EYFP+ groups after competitive microglial repopulation, which indicates that the infiltrating monocytes, imprinted by the CNS microenvironment, are different from the resident repopulated microglia, and confirming our previous publication .
Microglia replacement therapy has been proposed for CNS diseases linked to microglial dysfunctions or gene mutations. Replacing microglia by genetically modified or engineered cells may hold promise for distinct CNS diseases, yet the factors regulating competitive engraftment of different populations including microglia, monocytes and engineered cells are poorly understood [67, 68]. CX3CR1 regulates microglia colonization and distribution in the brain , and Cx3cr1 gene-deleted mice exhibited lower microglial density in the developing brain . The CX3CL1–CX3CR1 axis also regulates microglial repopulation following microglial depletion in the mouse retina , but whether this axis regulates competitive microglial repopulation in the brain has not previously been addressed. Our results indicate that repopulating microglia originate from three different predominant sources, resident microglia, including both Cx3cr1wt/wtCre−Eyfp− and Cx3cr1CreER-Eyfp/wt microglia, competing with infiltrating peripheral-derived microglia-like cells. Our data demonstrate that Cx3cr1wt/wtCre−Eyfp− microglia have a competitive advantage over Cx3cr1CreER-Eyfp/wt microglia. Furthermore, resident microglia lacking Cx3cr1 (from Cx3cr1CreER-Eyfp/CreER-Eyfp mice) were unable to compete with the peripheral-derived microglia-like cells following microglial depletion. This indicates that the resident microglia repopulation, but not peripheral derived microglia-like cell repopulation relies on CX3CR1. The proliferation rate of microglia in vitro was decreased in Cx3cr1 gene-depleted primary microglia, and this was not affected when the cells were challenged with CX3CL1 protein. This likely explains the repopulation advantage of Cx3cr1wt/wtCre−Eyfp− over Cx3cr1CreER-Eyfp/wt microglia post-depletion, at least partly. Cx3cr1 deficiency decreased the microglial migration rate, consistent with previous studies reporting that Cx3cr1 deficiency impaired microglia migration in vivo [71, 72] and in vitro . However, microglial migration was not impaired in Cx3cr1CreER-Eyfp/wt microglia. CX3CL1 increased migration rates of Cx3cr1wt/wtCre−Eyfp− microglia but not of Cx3cr1CreER-Eyfp/wt and Cx3cr1CreER-Eyfp/CreER-Eyfp microglia, which indicates that CX3CL1–CX3CR1 regulates microglial migration. The migration rate is higher in Cx3cr1wt/wtCre−Eyfp− than Cx3cr1wt/wtCre−Eyfp− microglia after adding CX3CL1, which further explains the repopulation advantage of Cx3cr1wt/wtCre−Eyfp− over Cx3cr1CreER-Eyfp/wt microglia post-depletion. Taken together, we conclude that the CX3CL1–CX3CR1 axis is important for the resident microglial repopulation and for competition with peripheral monocyte-derived microglia-like cells. Thus, limiting residential microglia repopulation by inhibiting CX3CL1–CX3CR1 signaling improves the microglial replacement efficiency by peripheral derived monocytes.
A small portion (less than 1%) of Cx3cr1wt/wtCre−Eyfp− microglia do not carry genetic labels in the widely used Cx3cr1GFP/wt or Cx3cr1CreER-Eyfp/wt mouse models. Not being aware of this population may lead to significant data misinterpretation since these cells may escape detection (not carrying the Eyfp or Gfp) and cannot be modified (lacking Cre expression) as expected. We further demonstrate the important role of the CX3CL1–CX3CR1 axis in regulation of microglial repopulation post-depletion. These findings raise an important cautionary note, not only when using the strains mentioned above but also for other strains that might display similar phenomena.
Availability of data and materials
The data supporting the findings are available upon request from the corresponding author. The RNA and DNA sequencing data generated during the current study are available in the public repository GEO with accession numbers GSE 186700 and PRJNA 778933, respectively.
Diphtheria toxin subunit A
Loss of heterozygosity
Central nervous system
- Eyfp :
Enhanced yellow fluorescent protein gene
Ginhoux F, Greter M, Leboeuf M, Nandi S, See P, Gokhan S, Mehler MF, Conway SJ, Ng LG, Stanley ER, et al. Fate mapping analysis reveals that adult microglia derive from primitive macrophages. Science. 2010;330:841–5.
Alliot F, Godin I, Pessac B. Microglia derive from progenitors, originating from the yolk sac, and which proliferate in the brain. Brain Res Dev Brain Res. 1999;117:145–52.
Tremblay ME, Stevens B, Sierra A, Wake H, Bessis A, Nimmerjahn A. The role of microglia in the healthy brain. J Neurosci. 2011;31:16064–9.
Nayak D, Roth TL, McGavern DB. Microglia development and function. Annu Rev Immunol. 2014;32:367–402.
Kierdorf K, Prinz M. Microglia in steady state. J Clin Invest. 2017;127:3201–9.
Salter MW, Stevens B. Microglia emerge as central players in brain disease. Nat Med. 2017;23:1018–27.
Spangenberg EE, Green KN. Inflammation in Alzheimer’s disease: lessons learned from microglia-depletion models. Brain Behav Immun. 2017;61:1–11.
Mondelli V, Vernon AC, Turkheimer F, Dazzan P, Pariante CM. Brain microglia in psychiatric disorders. Lancet Psychiatry. 2017;4:563–72.
Bisht K, Sharma KP, Lecours C, Sanchez MG, El Hajj H, Milior G, Olmos-Alonso A, Gomez-Nicola D, Luheshi G, Vallieres L, et al. Dark microglia: a new phenotype predominantly associated with pathological states. Glia. 2016;64:826–39.
Wolf SA, Boddeke HW, Kettenmann H. Microglia in physiology and disease. Annu Rev Physiol. 2017;79:619–43.
Colonna M, Butovsky O. Microglia function in the central nervous system during health and neurodegeneration. Annu Rev Immunol. 2017;35:441–68.
Pyonteck SM, Akkari L, Schuhmacher AJ, Bowman RL, Sevenich L, Quail DF, Olson OC, Quick ML, Huse JT, Teijeiro V, et al. CSF-1R inhibition alters macrophage polarization and blocks glioma progression. Nat Med. 2013;19:1264–72.
Quail DF, Bowman RL, Akkari L, Quick ML, Schuhmacher AJ, Huse JT, Holland EC, Sutton JC, Joyce JA. The tumor microenvironment underlies acquired resistance to CSF-1R inhibition in gliomas. Science. 2016;352: aad3018.
Wes PD, Sayed FA, Bard F, Gan L. Targeting microglia for the treatment of Alzheimer’s disease. Glia. 2016;64:1710–32.
Michell-Robinson MA, Touil H, Healy LM, Owen DR, Durafourt BA, Bar-Or A, Antel JP, Moore CS. Roles of microglia in brain development, tissue maintenance and repair. Brain. 2015;138:1138–59.
Mecca C, Giambanco I, Donato R, Arcuri C. Microglia and aging: the role of the TREM2-DAP12 and CX3CL1-CX3CR1 axes. Int J Mol Sci. 2018;19:318.
Harrison JK, Jiang Y, Chen S, Xia Y, Maciejewski D, McNamara RK, Streit WJ, Salafranca MN, Adhikari S, Thompson DA, et al. Role for neuronally derived fractalkine in mediating interactions between neurons and CX3CR1-expressing microglia. Proc Natl Acad Sci USA. 1998;95:10896–901.
Elmore MR, Najafi AR, Koike MA, Dagher NN, Spangenberg EE, Rice RA, Kitazawa M, Matusow B, Nguyen H, West BL, Green KN. Colony-stimulating factor 1 receptor signaling is necessary for microglia viability, unmasking a microglia progenitor cell in the adult brain. Neuron. 2014;82:380–97.
Jung S, Aliberti J, Graemmel P, Sunshine MJ, Kreutzberg GW, Sher A, Littman DR. Analysis of fractalkine receptor CX(3)CR1 function by targeted deletion and green fluorescent protein reporter gene insertion. Mol Cell Biol. 2000;20:4106–14.
Wolf Y, Yona S, Kim KW, Jung S. Microglia, seen from the CX3CR1 angle. Front Cell Neurosci. 2013;7:26.
Tay TL, Sagar, Dautzenberg J, Grun D, Prinz M. Unique microglia recovery population revealed by single-cell RNAseq following neurodegeneration. Acta Neuropathol Commun. 2018;6:87.
Matcovitch-Natan O, Winter DR, Giladi A, Vargas Aguilar S, Spinrad A, Sarrazin S, Ben-Yehuda H, David E, Zelada Gonzalez F, Perrin P, et al. Microglia development follows a stepwise program to regulate brain homeostasis. Science. 2016;353: aad8670.
Wieghofer P, Prinz M. Genetic manipulation of microglia during brain development and disease. Biochim Biophys Acta. 2016;1862:299–309.
Tay TL, Mai D, Dautzenberg J, Fernandez-Klett F, Lin G, Sagar, Datta M, Drougard A, Stempfl T, Ardura-Fabregat A, et al. A new fate mapping system reveals context-dependent random or clonal expansion of microglia. Nat Neurosci. 2017;20:793–803.
Parkhurst CN, Yang G, Ninan I, Savas JN, Yates JR, Lafaille JJ, Hempstead BL, Littman DR, Gan WB. Microglia promote learning-dependent synapse formation through brain-derived neurotrophic factor. Cell. 2013;155:1596–609.
Zhang Y, Zhao L. Repopulating retinal microglia restore endogenous organization and function under CX3CL1–CX3CR1 regulation. Sci Adv. 2018;4:8492.
Lund H, Pieber M. Fatal demyelinating disease is induced by monocyte-derived macrophages in the absence of TGF-beta signaling. Nat Immunol. 2018;19:1–7.
Xu Z, Rao Y, Huang Y, Zhou T, Feng R, Xiong S, Yuan TF, Qin S, Lu Y, Zhou X, et al. Efficient strategies for microglia replacement in the central nervous system. Cell Rep. 2020;32: 108041.
Fonseca MI, Chu SH, Hernandez MX, Fang MJ, Modarresi L, Selvan P, MacGregor GR, Tenner AJ. Cell-specific deletion of C1qa identifies microglia as the dominant source of C1q in mouse brain. J Neuroinflamm. 2017;14:48.
Wieghofer P, Knobeloch KP, Prinz M. Genetic targeting of microglia. Glia. 2015;63:1–22.
Merlini M, Rafalski VA, Ma K, Kim KY, Bushong EA. Microglial G(i)-dependent dynamics regulate brain network hyperexcitability. Nat Neurosci. 2020;24(1):19–23.
Schmitz C, Hof PR. Recommendations for straightforward and rigorous methods of counting neurons based on a computer simulation approach. J Chem Neuroanat. 2000;20:93–114.
Morrison HW, Filosa JA. A quantitative spatiotemporal analysis of microglia morphology during ischemic stroke and reperfusion. J Neuroinflamm. 2013;10:4.
Fontainhas AM, Wang M, Liang KJ, Chen S, Mettu P, Damani M, Fariss RN, Li W, Wong WT. Microglial morphology and dynamic behavior is regulated by ionotropic glutamatergic and GABAergic neurotransmission. PLoS ONE. 2011;6: e15973.
Fernández-Arjona MDM, Grondona JM, Granados-Durán P, Fernández-Llebrez P, López-Ávalos MD. Microglia morphological categorization in a rat model of neuroinflammation by hierarchical cluster and principal components analysis. Front Cell Neurosci. 2017;11:235.
Zachariadis V, Cheng H, Andrews N, Enge M. A highly scalable method for joint whole-genome sequencing and gene-expression profiling of single cells. Mol Cell. 2020;80:541-553.e545.
Osman AM, Sun Y, Burns TC, He L, Kee N, Oliva-Vilarnau N, Alevyzaki A, Zhou K, Louhivuori L, Uhlén P, Hedlund E, Betsholtz C, Lauschke VM, Kele J, Blomgren K. Radiation triggers a dynamic sequence of transient microglial alterations in juvenile brain. Cell Rep. 2020. https://doi.org/10.1016/j.celrep.2020.107699.
Lund H, Pieber M, Parsa R, Han J, Grommisch D, Ewing E, Kular L, Needhamsen M, Espinosa A, Nilsson E, et al. Competitive repopulation of an empty microglial niche yields functionally distinct subsets of microglia-like cells. Nat Commun. 2018;9:4845.
Bruttger J, Karram K, Wortge S, Regen T, Marini F, Hoppmann N, Klein M, Blank T, Yona S, Wolf Y, et al. Genetic cell ablation reveals clusters of local self-renewing microglia in the mammalian central nervous system. Immunity. 2015;43:92–106.
Butovsky O, Jedrychowski MP, Moore CS, Cialic R, Lanser AJ, Gabriely G, Koeglsperger T, Dake B, Wu PM, Doykan CE, et al. Identification of a unique TGF-beta-dependent molecular and functional signature in microglia. Nat Neurosci. 2014;17:131–43.
Bennett ML, Bennett FC, Liddelow SA, Ajami B, Zamanian JL, Fernhoff NB, Mulinyawe SB, Bohlen CJ, Adil A, Tucker A, et al. New tools for studying microglia in the mouse and human CNS. Proc Natl Acad Sci USA. 2016;113:E1738–46.
Van Hove H, Martens L, Scheyltjens I, De Vlaminck K, Pombo Antunes AR, De Prijck S, Vandamme N, De Schepper S, Van Isterdael G, Scott CL, et al. A single-cell atlas of mouse brain macrophages reveals unique transcriptional identities shaped by ontogeny and tissue environment. Nat Neurosci. 2019;22:1021–35.
Huang Y, Xu Z, Xiong S, Sun F, Qin G, Hu G, Wang J, Zhao L, Liang YX, Wu T, et al. Repopulated microglia are solely derived from the proliferation of residual microglia after acute depletion. Nat Neurosci. 2018;21:530–40.
Zabel MK, Zhao L, Zhang Y, Gonzalez SR, Ma W, Wang X, Fariss RN, Wong WT. Microglial phagocytosis and activation underlying photoreceptor degeneration is regulated by CX3CL1-CX3CR1 signaling in a mouse model of retinitis pigmentosa. Glia. 2016;64:1479–91.
Lee S, Varvel NH, Konerth ME, Xu G, Cardona AE, Ransohoff RM, Lamb BT. CX3CR1 deficiency alters microglial activation and reduces beta-amyloid deposition in two Alzheimer’s disease mouse models. Am J Pathol. 2010;177:2549–62.
Thome AD, Standaert DG, Harms AS. Fractalkine signaling regulates the inflammatory response in an α-synuclein model of Parkinson disease. PLoS ONE. 2015;10: e0140566.
Castro-Sánchez S, García-Yagüe ÁJ, Kügler S, Lastres-Becker I. CX3CR1-deficient microglia shows impaired signalling of the transcription factor NRF2: implications in tauopathies. Redox Biol. 2019;22: 101118.
Gyoneva S, Hosur R, Gosselin D, Zhang B, Ouyang Z, Cotleur AC. Cx3cr1-deficient microglia exhibit a premature aging transcriptome. Life Sci. 2019. https://doi.org/10.26508/lsa.201900453.
Corona AW, Huang Y, O’Connor JC, Dantzer R, Kelley KW, Popovich PG, Godbout JP. Fractalkine receptor (CX3CR1) deficiency sensitizes mice to the behavioral changes induced by lipopolysaccharide. J Neuroinflamm. 2010;7:93.
Nichols CA, Gibson WJ, Brown MS, Kosmicki JA. Loss of heterozygosity of essential genes represents a widespread class of potential cancer vulnerabilities. Nat Commun. 2020;11:2517.
Shao C, Deng L, Henegariu O, Liang L, Raikwar N, Sahota A, Stambrook PJ, Tischfield JA. Mitotic recombination produces the majority of recessive fibroblast variants in heterozygous mice. Proc Natl Acad Sci USA. 1999;96:9230–5.
Fisher G, Stephenson DA, West JD. Investigation of the potential for mitotic recombination in the mouse. Mutat Res. 1986;164:381–8.
Mortensen RM, Conner DA, Chao S, Geisterfer-Lowrance AA, Seidman JG. Production of homozygous mutant ES cells with a single targeting construct. Mol Cell Biol. 1992;12:2391–5.
Han J, Zhu K, Zhang XM, Harris RA. Enforced microglial depletion and repopulation as a promising strategy for the treatment of neurological disorders. Glia. 2019;67:217–31.
Willis EF, MacDonald KPA, Nguyen QH, Garrido AL, Gillespie ER, Harley SBR, Bartlett PF, Schroder WA, Yates AG, Anthony DC, et al. Repopulating microglia promote brain repair in an IL-6-dependent manner. Cell. 2020;180:833-846.e816.
Bruttger J, Karram K, Wörtge S, Regen T, Marini F, Hoppmann N, Klein M, Blank T, Yona S, Wolf Y, et al. Genetic cell ablation reveals clusters of local self-renewing microglia in the mammalian central nervous system. Immunity. 2015;43:92–106.
Miron VE, Priller J. Investigating microglia in health and disease: challenges and opportunities. Trends Immunol. 2020;41:785–93.
Salter MW, Beggs S. Sublime microglia: expanding roles for the guardians of the CNS. Cell. 2014;158:15–24.
Parkhurst CN, Gan WB. Microglia dynamics and function in the CNS. Curr Opin Neurobiol. 2010;20:595–600.
Hoeffel G, Chen J, Lavin Y, Low D, Almeida FF, See P, Beaudin AE, Lum J, Low I, Forsberg EC, et al. C-Myb(+) erythro-myeloid progenitor-derived fetal monocytes give rise to adult tissue-resident macrophages. Immunity. 2015;42:665–78.
Gomez Perdiguero E, Klapproth K, Schulz C, Busch K, Azzoni E, Crozet L, Garner H, Trouillet C, de Bruijn MF, Geissmann F, Rodewald HR. Tissue-resident macrophages originate from yolk-sac-derived erythro-myeloid progenitors. Nature. 2015;518:547–51.
Bain CC, Bravo-Blas A, Scott CL, Perdiguero EG, Geissmann F, Henri S, Malissen B, Osborne LC, Artis D, Mowat AM. Constant replenishment from circulating monocytes maintains the macrophage pool in the intestine of adult mice. Nat Immunol. 2014;15:929–37.
Molawi K, Wolf Y, Kandalla PK, Favret J, Hagemeyer N, Frenzel K, Pinto AR, Klapproth K, Henri S, Malissen B, et al. Progressive replacement of embryo-derived cardiac macrophages with age. J Exp Med. 2014;211:2151–8.
Mossadegh-Keller N, Gentek R, Gimenez G, Bigot S, Mailfert S, Sieweke MH. Developmental origin and maintenance of distinct testicular macrophage populations. J Exp Med. 2017;214:2829–41.
Ginhoux F, Guilliams M. Tissue-resident macrophage ontogeny and homeostasis. Immunity. 2016;44:439–49.
Herz J, Filiano AJ, Smith A, Yogev N, Kipnis J. Myeloid cells in the central nervous system. Immunity. 2017;46:943–56.
Prinz M, Jung S, Priller J. Microglia biology: one century of evolving concepts. Cell. 2019;179:292–311.
Priller J, Prinz M. Targeting microglia in brain disorders. Science. 2019;365:32–3.
Combadiere C, Feumi C, Raoul W, Keller N, Rodero M, Pezard A, Lavalette S, Houssier M, Jonet L, Picard E, et al. CX3CR1-dependent subretinal microglia cell accumulation is associated with cardinal features of age-related macular degeneration. J Clin Invest. 2007;117:2920–8.
Paolicelli RC, Bolasco G, Pagani F, Maggi L, Scianni M, Panzanelli P, Giustetto M, Ferreira TA, Guiducci E, Dumas L, et al. Synaptic pruning by microglia is necessary for normal brain development. Science. 2011;333:1456–8.
Arnoux I, Audinat E. Fractalkine signaling and microglia functions in the developing brain. Neural Plast. 2015;2015: 689404.
Liang KJ, Lee JE, Wang YD, Ma W, Fontainhas AM, Fariss RN, Wong WT. Regulation of dynamic behavior of retinal microglia by CX3CR1 signaling. Invest Ophthalmol Vis Sci. 2009;50:4444–51.
Castro-Sanchez S, Garcia-Yague AJ, Kugler S, Lastres-Becker I. CX3CR1-deficient microglia shows impaired signalling of the transcription factor NRF2: implications in tauopathies. Redox Biol. 2019;22: 101118.
We would like to thank: The biomedical flow cytometry core facility, Karolinska Institutet and Center for molecular medicine flow cytometry facility, Karolinska Institutet for cell sorting; the core facility Bioinformatics and Expression Analysis (BEA), Karolinska Institutet, for the bulk RNA sequencing and whole-genome sequencing; the Eukaryotic Single Cell Genomics Facility (ESCG, Karolinska Institutet) for the single-cell RNA sequencing. The data handling was enabled by resources provided by the Swedish National Infrastructure for Computing (SNIC) at UPPMAX, partially funded by the Swedish Research Council through Grant Agreement No. 2018-05973.
Open access funding provided by Karolinska Institutet. The Swedish Research Council (2019-01684 to K.B.), the Swedish Childhood Cancer Fund (PR2018-0158, PROF20/0001 to K.B.), Swedish Cancer Foundation (20 1297 PjF to K.B.), Swedish Government grants to scientists working in healthcare (LS 2017-1341 to K.B.), Hjärnfonden (FO2018-0257 to K.B.), Cancer Research KI Translational Seed Funding (to K.B.), the Radiumhemmet Research Funds (201072 to K.B.), the Gunnar and Märta Philipson Foundation (to K.B.), and the Frimurare Barnhuset Foundation of Stockholm (to K.B.).
All experiments in this study were approved and performed following the Swedish National Board for Laboratory Animals’ guidelines and the European Community Council Directive (86/609/EEC) under the ethical permits N163/15, N127/16, and 13676-2020.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Additional data related to Fig. 1. A Representative 8 sections of Hoechst staining of Cx3cr1CreER-Eyfp/wt mouse brain sequential sagittal slices, the numbers representing the location and number of EYFP− microglia in that location, each color represents one mouse, n = 6. Figure S2. Morphology analysis. A The process to prepare topological skeleton from the original photomicrographs. The full-size maximum intensity projection images (top) were processed into binary images (middle) and then skeletonized (bottom) following the ImageJ plugin protocol. Cropped images (right) from the original full-size images (left) were shown to improve the visualization. Scale bar, 50 μm. B The Analyze Skeleton plugin was applied, and skeletonized endpoints were tagged purple, the slab is orange, and the junction is pink. The tagged data are summarized as total length (sum of endpoints, slab, and junction) and the total number of endpoints. C The process area is represented as the convex hull area by connecting the process ends using the ImageJ polygon tool. Figure S3. Additional data related to Fig. 1. A Gating strategies of flow cytometry analyses. B EYFP− microglia ratio of total microglia at 3 weeks, 4 weeks, 6 weeks, 9 weeks and 15 weeks old mice, n = 4, 4, 6, 4, respectively, mean ± s.d. No significant difference by one-way ANOVA. C UMAP plots from single-cell sequencing, each dot represents one cell, the color represents the expression levels of Cx3cr1. D Histogram plot of Cx3cr1. The x-axis represents Cx3cr1 expression read counts; the y-axis represents cell numbers. Figure S4. Additional data related to Fig. 2. A Graph showing relative Iba-1+ microglia at 1 day, 3 days and 7 days after the final Tam injection, respectively; n = 3–4, mean ± s.d. *p < 0.05, **p < 0.01, ***p < 0.001 by Student’s two-tailed unpaired t test. B Representative images of Iba-1 microglia staining of Cx3cr1CreER-Eyfp/wt mice and Cx3cr1CreER-Eyfp/wt Rosa26DTA/wt mice at 1 day after the final Tam injection in cortex and cerebellum. Scale bar, 200 μm. C Representative 3 sequential sagittal sections (25 µm/section, with 36 sections interval) after Hoechst staining of Cx3cr1CreER-Eyfp/wt Rosa26DTA/wt mice brains. The numbers represent the location and number of Cx3cr1highCre−Eyfp− microglia at 1 day and 3 days after the final Tam injection, each color represents one mouse, n = 3. D Quantitative data showing the total number of microglia in the 3 sections in Ctrl group, 1 day and 3 days after the final Tam injection, n = 4, 3, 3, mean ± s.d. Figure S5. Additional data related to Fig. 4. A–C Bar graphs showing RPKM of Cre, EYFP, and Cx3cr1, respectively, n = 5, mean ± s.d. ***p < 0.001 by Student’s two-tailed unpaired t test. D The expression of CX3CR1 in both Cx3cr1CreER-Eyfp/wtCre+Eyfp+ and Cx3cr1highCre−Eyfp− microglia in Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice at day 42 after the final Tam injection. E Bar graph showing the mean fluorescence intensity (MFI) of CX3CR1, n = 4, mean ± s.d. ***p < 0.001 by Student’s two-tailed unpaired t test. Figure S6. Additional data related to Fig. 4. A Schematic diagram showed how the primers were designed. The red color indicates the common primers used in all the PCR reactions; the green and blue color indicate the WT primer and Mutant primer used in PCR reaction for Fig. 4f; the orange color indicates the long PCR product targeting primer used in PCR reaction for Fig. 4g. B Whole-genome sequencing of Cx3cr1highCre−Eyfp− microglia in mouse 2. The top panel is the sequencing mapping on the CreERT2-Eyfp vector; the middle panel is the sequencing mapping on chromosome 9; the bottom panel is the schematic graph showing the homozygosity of the exon 2 of Cx3cr1. C Whole-genome sequencing of Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia. The top panel is the sequencing mapping CreERT2-Eyfp vector; the middle panel is the sequencing mapping the chromosome 9; the bottom panel is the schematic graph showing the heterogeneity of the exon 2 of Cx3cr1. D DNA electrophoresis of PCR product from EYFP− and EYFP+ microglia. E Representative dot plots showing the percentage of EYFP+F4/80low, EYFP−F4/80low, EYFP+F4/80high, and EYFP−F4/80high repopulated microglia at day 42 after Tam in Cx3cr1CreER-Eyfp/ CreER-Eyfp Rosa26DTA/wt mice, mean ± s.d.
Additional file 2: Table S1.
Gene list including Gene symbols, Log2 fold changes, P-values and FDR values. The table shows FDR values and P-values in ascending order.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhou, K., Han, J., Lund, H. et al. An overlooked subset of Cx3cr1wt/wt microglia in the Cx3cr1CreER-Eyfp/wt mouse has a repopulation advantage over Cx3cr1CreER-Eyfp/wt microglia following microglial depletion. J Neuroinflammation 19, 20 (2022). https://doi.org/10.1186/s12974-022-02381-6
- Microglial depletion
- Microglial repopulation
- Diphtheria toxin subunit A (DTA)
- Loss of heterozygosity (LOH)
- Homologous recombination