- Research
- Open access
- Published:
Repeated LPS induces training and tolerance of microglial responses across brain regions
Journal of Neuroinflammation volume 21, Article number: 233 (2024)
Abstract
Background
Neuroinflammation is involved in the pathogenesis of almost every central nervous system disorder. As the brain’s innate immune cells, microglia fine tune their activity to a dynamic brain environment. Previous studies have shown that repeated bouts of peripheral inflammation can trigger long-term changes in microglial gene expression and function, a form of innate immune memory.
Methods and results
In this study, we used multiple low-dose lipopolysaccharide (LPS) injections in adult mice to study the acute cytokine, transcriptomic, and microglia morphological changes that contribute to the formation of immune memory in the frontal cortex, hippocampus, and striatum, as well as the long-term effects of these changes on behavior. Training and tolerance of gene expression was shared across regions, and we identified 3 unique clusters of DEGs (2xLPS-sensitive, 4xLPS-sensitive, LPS-decreased) enriched for different biological functions. 2xLPS-sensitive DEG promoters were enriched for binding sites for IRF and NFkB family transcription factors, two key regulators of innate immune memory. We quantified shifts in microglia morphological populations and found that while the proportion of ramified and rod-like microglia mostly remained consistent within brain regions and sexes with LPS treatment, there was a shift from ameboid towards hypertrophic morphological states across immune memory states and a dynamic emergence and resolution of events of microglia aligning end-to-end with repeated LPS.
Conclusions
Together, findings support the dynamic regulation of microglia during the formation of immune memories in the brain and support future work to exploit this model in brain disease contexts.
Background
Innate immune cells dynamically adjust their responses to immune challenges through innate immune memory, a process influenced by their past encounters with inflammatory events [1,2,3,4]. As the brain’s resident immune cells, microglia rapidly respond to inflammation and engage in classic immune functions including cytokine and chemokine production and phagocytosis of pathogens and cellular debris. In addition to their responsibilities as immune cells, microglia play critical roles in establishing and maintaining brain homeostasis including the regulation of neuron and oligodendrocyte cell numbers [5], shaping of brain circuitry [6], and fine-tuning of neuronal connections [3, 7, 8]. Many of these processes are disrupted in brain disorders, supporting a key role for microglia in health and disease. Several clinical and post-mortem studies point to a neuroinflammatory component of brain disorder etiologies, including Schizophrenia, Autism Spectrum Disorder, and Alzheimer’s Disease, driven by increased activation of microglia [9,10,11,12]. Microglia release and respond to cytokines, chemokines, and neurotransmitters in the brain [13] to communicate with other cell types and regulate their functions. Thus, disrupted microglial responses can impact communication across cell types in the brain environment and potentially lead to downstream alterations in brain function.
Innate immune cells update their responses to inflammatory events based on previous exposures to immune stimuli [4, 14, 15]. After the resolution of an initial immune activating event, innate immune cells can remain in a ‘trained’ state, characterized by exaggerated responses to subsequent immune challenges, or in a ‘tolerized’ state characterized by blunted responses [4]. Both forms of innate immune memory serve as mechanisms that enhance an organism’s ability to respond effectively to recurring infections or damage, but can become detrimental if dysregulated. While immune training allows proactive adaptation of immune responses that have a protective benefit to the organism, this process can become maladaptive in the context of disease or injury when inflammation cannot be resolved [4]. Similarly, tolerized immune responses play a crucial role in averting chronic inflammatory conditions when encountering commensal microbes [4, 16]. However, these responses can be counterproductive when facing a new pathogen. As the resident innate immune cells of the brain, microglia are unique in that they are long-lived cells that locally self-renew over the lifespan [17, 18], suggesting that formation of innate immune memories in microglia could have long-lasting impacts on brain function and disease development throughout life.
Innate immune cells, including microglia, express receptors that respond to microbe-associated molecular patterns. Microglia Toll-like receptor 4 (TLR4) binds lipopolysaccharide (LPS), a major structural component of gram-negative bacteria. Repeated peripheral injections of low-dose LPS has recently been demonstrated as a robust model to study microglial immune memory responses in the brain [4, 19]. Repeated LPS induces long-lasting immune training and tolerance in brain microglia that has been shown to persist for at least 6 months [19]. Similarly, LPS preconditioning has also been shown to reduce damage in traumatic brain injury by preventing tolerance [20]. Interestingly, Wendeln et al. [19] showed that pathological hallmarks of Alzheimer’s disease and stroke in mouse models are exacerbated by LPS-induced immune training and alleviated by immune tolerance. Long-term functional changes in microglia were paralleled by long-term reprogramming of microglial enhancers and gene expression in these models. Genes involved in hypoxia-induced factor-a (HIF-1α) regulation, metabolic glycolysis, and Rap1-mediated phagocytosis were increased several months after training but not tolerance with repeated LPS. Together, findings in the literature support a model in which long-term reprogramming of gene expression underlies long-term innate immune memory in the brain relevant to disease contexts. However, the mechanisms involved in the initial formation of innate immune memories are not as well characterized in microglia and hence form the major focus of this study.
Here, we used a mouse model of repeated LPS to identify how gene expression, cytokine expression, and microglia morphology changes drive the formation of training and tolerance in the brain. Unlike prior studies using the same repeated low-dose LPS paradigm, ours focused specifically on characterizing immediate, acute microglial changes across multiple brain regions including the frontal cortex, striatum, and hippocampus, which have known differences in microglial morphology and function [21, 22]. Furthermore, we assessed long-term impacts of innate immune memory induction in a comprehensive battery of anxiety-like, depressive-like, repetitive, and learning and memory behaviors and performed a systematic evaluation of LPS-induced shifts in microglia morphological states across brain regions. We describe shared patterns of gene expression across brain regions with shared and unique regulatory pathways. We also describe the dynamic emergence and resolution of events of microglia aligning end-to-end with the formation of immune priming and tolerance, a morphological phenotype which has not previously been described in the context of repeated LPS. Together, these findings offer new insight into microglial changes in the formation of immune memory and support future work to further explore the described mechanisms in the context of brain disease.
Methods
Animals
All experiments were conducted in accordance with the National Institutes of Health Guidelines for the Care and Use of Laboratory Animals, with approval from the Institutional Animal Care and Use Committee at the University of California, Davis and the Canadian Council on Animal Care guidelines, with approval from the Animal Care Committee at the University of British Columbia. Mice in experiments undertaken at the University of California, Davis (RNA-seq, Luminex protein array) were housed in treatment-matched groups of two to four on a regular 12-h light/12-h dark cycle and all experiments were performed in regular light during the mouse’s regular light cycle. Mice in experiments undertaken at the University of British Columbia (RT-qPCR, microglia morphology, behavioral battery) were housed in treatment-matched groups of two to four on a reversed 12-h light/12-h dark cycle. Because mice are nocturnal animals, all behavior experiments were performed in red light during the mouse’s dark cycle when they are most active. All mice had ad libitum food and water and cages were changed biweekly.
In vivo injections
Adult male and female C57BL/6J mice (10–18 weeks old) were given daily intraperitoneal injections of 0.5 mg/kg LPS (from E. coli O55:B5, Sigma-Aldrich L5418) or an equal volume of vehicle solution (1xPBS; Phosphate buffered saline, Fisher Scientific BP3991) for 4 days between 09:00 and 10:00 on each day. Animals received either four PBS vehicle injections on four consecutive days (PBS), three vehicle injections on three consecutive days followed by a single LPS injection on the fourth day (1xLPS), two vehicle injections on two consecutive days followed by two LPS injections on the following two days (2xLPS), a single vehicle injection followed by three LPS injections on the following three days (3xLPS, males only), or four LPS injections (4xLPS). PBS, 1xLPS, 2xLPS, 3xLPS, and 4xLPS were the five experimental groups analyzed in this study (Fig. 1A). Body weight was measured right before administering injections each day.
Sickness behavior scoring and analysis
Changes in sickness behavior were assessed 3 h-post injection each day. Sickness behavior was scored according to a 4-point grading scale (0 to 3 points total). Lethargy signified by diminished locomotion and curled body posture, ptosis signified by drooping eyelids, and piloerection signified by ruffled and greasy fur were each assessed individually for either 0, 0.5, or 1 point. A total cumulative score of 0 indicated no symptoms of sickness behavior and a total cumulative score of 3 points indicated that all symptoms of sickness behavior were maximally present. We assessed how cumulative scores for each mouse on the last day of assessment (3 h after final injection) changed with LPS treatment by fitting a non-parametric Aligned Ranks Transformation (ART) model to measure scores as a factor of LPS Treatment (Total Score ~ Treatment) using the ARTool R package [23]. The model was assessed using ART ANOVA and tests between Treatment groups were corrected for multiple comparisons using the Benjamini-Hochberg (BH) method. Corrected p-values < 0.05 were considered statistically significant.
Behavioral battery
Animals underwent a battery of behavioral tasks measuring anxiety-like behavior, repetitive behavior, learning and memory, and depressive-like behaviors. All behavioral tests were performed during the dark cycle under red light (40–100 lx) between the hours of 09:00–17:00. Behavioral testing began on day 9, or 5 days after the last injection, and ended on day 50 (Fig. 2B). We made the choice to start behavior on day 9 of the experiment to ensure that the animals recovered as fully as possible from the LPS paradigm (stabilization in weights, resolution of sickness behaviors) before assessing changes in the behavioral battery. This was important to be able to mitigate any effects of ongoing recovery on task performance and allowed us to specifically test how exposure to repeated LPS affects behaviors beyond the acute impacts of sickness. ANY-maze tracking software was used to record and track animal movement during the task. All statistical analyses of behavior were performed using GraphPad Prism 9 software. (Supp. File S1) Behavior was analyzed using a one-way ANOVA assessing treatment as the main effect and using the BH-correction for multiple comparisons. Corrected p-values < 0.05 were considered statistically significant. While the individual points in Fig. 2B are colored by sex for visualization purposes, the statistical analysis was run on treatment groups with sex combined (n = 4–6).
Light dark box
The light–dark box test was performed to measure anxiety-like behavior. The light chamber was an open-top 25 cm × 25 cm × 30 cm (lwh) white box. The dark chamber was a closed-top 16 cm × 25 cm × 30 cm (lwh) black box that was connected to the light box through a 10 cm × 5 cm (wh) passage and sliding door. Mice were individually placed into the dark chamber and were allowed to freely explore for 10 min once the sliding door was opened. Video footage was recorded from above and ANY-maze recorded time spent in the light chamber.
Open field test
The open field test was performed to measure anxiety-like behavior. The box was 65 cm × 65 cm × 40 cm (lwh). One at a time, mice were placed in the center of the box and allowed to explore freely for 20 min. Video footage was recorded from above and ANY-maze recorded time spent in the inner and outer zones. The zone boundary was defined at 10 cm from the box wall.
Marble burying test
The marble burying test was performed to measure repetitive behavior. To prepare for the task, a 30 cm × 19.4 cm × 39.4 cm rodent cage was filled with bedding 5 cm deep and 20 marbles were equally spaced on top of the bedding in a 4 × 5 grid using a premade stencil. Mice were placed into the centre of the box one at a time and recorded for 20 min. Video footage subsequently hand scored at a later time to count for the number of marbles buried after 20 min. A marble was scored as buried if < 50% of the marble was visible.
Object location memory test
The OLM task was performed in an open-topped square box 40 cm × 40 cm × 40 cm (lwh) lined with bedding ~ 1 cm in depth. One wall of the box was lined with duct tape as a visual cue. Objects were cleaned with 1:16 Saber solution (Rododonto, Cat # 09-12400-04) between trials and allowed to fully dry before re-use. The OLM took 8 days in total with 6 days of habituation, 1 day of training and 1 day of testing [24]. Bedding was kept in the same box each day but shuffled around between mice. Mice were habituated for 10 min/day for 6 days, where they were allowed to freely explore the same box each day with no objects present. On the training day, mice were given 10 min to investigate two identical glass beakers filled with dry cement. The 5-min testing phase took place 24 h after training. During the test, one object was moved to a novel location and the other remained in the same spot. An experienced scorer blinded to experimental conditions hand-scored the testing videos for investigation time with the unmoved and moved objects. Investigation time was scored when the mouse’s nose was within 1 cm of the object and pointing directly at the object. Looking over or past the object and climbing on top of the object were not counted as investigation. For the training and testing days, a discrimination index (DI) was calculated as [(TM − TS)/(TM + TS)] * 100, where TM is time investigating the moved object and TS is time investigating the unmoved object.
Self-grooming
Self-grooming was performed to measure repetitive-like behaviour. Mice were placed in a cage without bedding and filmed straight-on for 20 min. The last 10 min of video footage were hand-scored for the number of grooming bouts.
Forced swim test
Mice were placed individually into 1L glass beakers filled with 900 mL of 23–25 °C water. Mice were recorded for 6 min and closely monitored during the task. Mice were placed in a cage with dry paper towel on top of a heating pad after the task and were returned to the home cage once they were dry and exhibiting normal motor and grooming behaviours. The last four minutes of footage was hand-scored by a researcher blind to experimental conditions for time spent immobile, defined as the absence of escape-related movements.
Weight change analysis
Body weights were measured daily (or from another perspective, 24 h after the prior day’s injection; Fig. 2A) before undergoing the behavioral battery described above. Mice received LPS or PBS injections on Days 1–4 and were left alone to recover other than weight measurements on Days 5–8 before starting the behavioral assessments (outlined in ‘Behavioral Battery’ section) on Day 9. Each animal’s weights were normalized to their weights taken before the start of injections to calculate % weight change. We assessed how % weight differences change with LPS treatment in relation to experiment day by fitting a linear mixed effects model to measure % weight change as a factor of treatment and experiment day interactions with mouseID as a repeated measure (% weight change ~ Treatment*Timepoint + (1|MouseID)). The model was assessed using ANOVA and two separate sets of posthoc analyses using the BH-method were run: (1) to correct across all comparisons made between timepoints within Treatment groups (~ Timepoint|Treatment) and (2) to correct across comparisons made against PBS (1-4xLPS vs. PBS) within each timepoint (~ Treatment|Timepoint). BH-corrected p-values < 0.05 were considered statistically significant.
Brain tissue collection
In each of the separate cohorts outlined here, mice were euthanized 3 h after the final injection immediately after the final day’s sickness behavior was scored. For the data presented in Figure S1, mice were euthanized 3 h, 24 h, 48 h, or 7 days after injections. Mice were quickly anesthetized with isofluorane and transcardially perfused with 15 mL of 1xPBS before brains were extracted for downstream experiments. For the RNAseq and Luminex cohort: one hemisphere from each brain was microdissected for brain regions (frontal cortex, striatum, hippocampus) for RNAseq and the other hemisphere was microdissected for brain regions (frontal cortex, striatum, hippocampus, cerebellum) for Luminex protein arrays, according to the brain microdissection protocol outlined in Chiu et al. [25]. Dissected brain regions were immediately flash-frozen on dry ice and stored at − 80 °C until RNA and protein extractions for RNAseq and Luminex protein arrays, respectively. For the microglia morphology & RT-qPCR cohort: one hemisphere from each brain was collected for downstream histology analysis and the other hemisphere was microdissected for brain regions (frontal cortex, striatum, hippocampus) for RT-qPCR. Extracted hemispheres for histology were immersion-fixed in 4% paraformaldehyde for 48 h before cryoprotecting in 30% sucrose for 48 h prior to cryosectioning. Cryoprotected brains were then sectioned at 30um on the cryostat, collected in 1xPBS for long-term storage, and processed for immunohistochemistry. Dissected brain regions were immediately flash-frozen on dry ice and stored at − 80 °C until RNA extractions for RT-qPCR. For the isolated microglia experiments (Fig. S6C, E): cortex from both hemispheres was collected for microglia isolations, subsequent RNA extractions, and downstream RT-qPCR.
RNA-seq library preparation
RNA was extracted from whole tissue in the hippocampus, striatum and frontal cortex using the RNA/DNA Miniprep Plus Zymo Kit (Cat D7003) with on-column DNase digestion as described in the kit manual. RNA quality was assessed by bioanalyzer and all samples had RNA Integrity Numbers > 7 (Supp. Table S3). QuantSeq 3ʹmRNA sequencing FWD (RNAseq) libraries (Lexogen, Cat No 015) were prepared per manufacturer’s instructions using 75 ng of total input RNA per sample and 17 cycles of PCR amplification. All samples were prepared with incorporation of Unique Molecular Identifiers as part of the second strand synthesis step (Lexogen, Cat No 081) and sample unique, dual index barcodes. All samples were then pooled, exonuclease VII treated and sequenced across two lanes of a HiSeq4000 to generate Single End 100 base pair reads.
RNA-Seq alignment
Raw fastq sequencing reads were assessed for quality control using FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and multiqc [26]. The UMI index was then added to the header of each read using custom scripts from Lexogen (umi2index). Reads were then trimmed to remove the first 11 bases, adapter contamination, polyA read through and low-quality reads using BBDuk from BBtools (http://www.sourceforge.net/projects/bbmap/). All samples were then re-assessed for quality using FastQC and multiqc. Each sample was then aligned to the mouse genome (GRCm38.92) using STAR [27]. Aligned bam files were then filtered for unique UMIs using custom Lexogen script (collapse_UMI_bam). Alignments for unfiltered and UMI filtered bam files were then compared using multiqc. All subsequent analysis was performed on UMI filtered bam files. These files were indexed using samtools [28] and then reads per ensembl gene (Mus_musculus.GRCm38.92.gtf) were counted using FeatureCounts [29] with options for forward strand (-s 1). Counts were then assess using multiqc and read into R for statistical analysis using EdgeR [30] and LimmaVoom [31]. Code is available on the CierniaLab github at https://github.com/ciernialab/Kim2024_1-4xLPS.
Differential expression (DE) analysis
Low expressing genes were removed by filtering for genes with more than 1 count per million (CPM) in at least 4 samples (Figure S2C). The remaining genes (n = 16,788) were then normalized using Trimmed Mean of M-values (TMM) to correct for library composition using calcNormFactors, method = TMM (Figure S2E). To evaluate the contribution of each factor to the overall gene expression profile, Multi-Dimensional Scaling (MDS) plots were constructed for Brain Region, LPS Treatment, and Hemisphere (Figure S2B). CPM values were then fed into voom using a model design matrix for LPS treatment differences in gene expression which vary with brain region (~ 0 + Treatment:BrainRegion). Individual contrast comparisons for each brain region were then called using contrasts.fit followed by eBayes. Differentially expressed genes (DEGs) were identified for each comparison of interest using Toptreat with a BH p-value correction (significance at adjusted p-values < 0.05 and fold change > ± 2). Heatmaps of normalized CPM values of significant DEGs (n = 432) were made using the pheatmap R package [32]. CPM values for each gene were first scaled within each brain region and then scaled across brain regions as input for the DEG heatmap (Fig. 3A), before hierarchical clustering. Breaks between each color of the heatmap scale were adjusted to fit the distribution of the normalized gene expression values such that 10% of the data was contained within each color break of the inferno color palette from the viridisLite R package (Fig. S2A). This allowed for the stark visualization of the strongest patterns of gene expression across brain regions and characterization of cluster IDs as 2xLPS-sensitive, 4xLPS-sensitive, and LPS-decreased.
Gene ontology (GO) enrichment analysis
Enrichment for gene ontology (GO) terms was subsequently performed on the DEG clusters using the clusterProfiler (version 4.4.4) [33, 34] and org.Mm.eg.db (version 3.15.0) [35] R packages to further explore the functional role of the identified genes. GO enrichment was run on all 3 DEG cluster lists together using the compareCluster function in clusterProfiler with the “enrichGO” option and a p-value cutoff < 0.01 and BH-adjusted q-value cutoff < 0.05 against the Biological Processes (BP) ontology. Dotplot and cnetplot functions within the clusterProfiler package was used to create visualizations of the significant GO enrichments.
MGEnrichment analysis
Gene list enrichment analysis was performed using the 2024 database of the MGEnrichmentApp [36] with all genes interrogated in the experiment as background. Significant enrichments against mouse microglial gene lists in MGEnrichment were tested by a one-tail Fischer’s exact test and BH-corrected to reach significance at an adjusted p-value < 0.01. Target % and Background % values for significant gene list enrichments were plotted using ggplot in R.
HOMER promoter motif enrichment analysis
HOMER Motif Analysis Software, v4.11 (http://homer.ucsd.edu/homer/) [37] was used to find motif enrichments within promoters (+ 1000, − 200 to transcription start site) for significant DEGs in each cluster identified (2xLPS-sensitive, 4xLPS-senstive, LPSDownregulated) and run against background genes (all genes interrogated in the experiment, n = 16,788) using the following command: findMotifs.pl [cluster-specific genes file] mm10_promoters [output directory]-bg [background genes file]-start-1000-end 200 -p 12. Enrichment output was interpreted according to recommendations on the HOMER website (http://homer.ucsd.edu/homer/). Top-ranked de novo motifs were defined as having p-values < 1e−10 and not marked as a possible false positive in the Homer de novo Motif Results output (Supp. File S2–S4). Top-ranked known motifs were defined as having BH-corrected q-values < 0.05 in the Homer Known Motif Enrichment Results output. (Supp. File S2–S4) Target % and Background % values for significant promoter enrichments were plotted using ggplot in R. All known motifs were concatenated into one file and used as input to determine target gene promoter locations with known motifs using the following command: findMotifs.pl [cluster-specific genes file] mm10_promoters [output directory]-find [all known motifs file] > [desired name of output file]-bg [background genes file]-start-1000-end 200-p 12. Identified instances of significant known motifs in each gene promoter were binarized (yes = present, no = not present) and plotted side-by-side with each gene’s normalized gene expression using the R package pheatmap.
Published ChIP-seq promoter enrichment analysis
Publicly available ChIP-seq peaks [38,39,40,41,42] from relevant microglia and bone marrow-derived macrophage (BMDM) studies were downloaded from the NCBI GEO repository and processed using sort-bed functions and -intersect and -difference options within the bedops software (version 2.4.41) [43] to isolate out ChIP-seq mm10 promoter peaks and to characterize lost vs. gained promoter peaks in any studies with treatment conditions, respectively. Final ChIP-seq peaks were curated into collections for each study to use as input for analysis, according to the requirements outlined in the LOLA R package. ChIP-seq peak promoter enrichment analysis was performed using LOLA and significant enrichments were defined as BH-corrected adjusted p-values < 0.05. Target % and Background % values for significant enrichments were plotted using ggplot in R.
Exploration of single-cell brain atlas datasets
Publicly available gene expression datasets of the mouse brain were downloaded from BrainRNASeq.org [44] and MouseBrain.org [45]. The heatmap generated for Fig. 3A was plotted alongside expression values scaled by row from each brain cell type present in the downloaded atlas datasets using the pheatmap R package (Fig. S5).
Luminex protein arrays
As described previously [46], blood was collected by cardiac punch and brain tissue for each region (frontal cortex, striatum, hippocampus, and cerebellum) was collected and flash frozen. Blood was placed on ice and subsequently processed into serum by centrifugation. Brain tissue was lysed in cell lysis buffer (Cell Signaling Technologies, Danvers, MA, USA) containing protease and phosphatase inhibitors and incubated with agitation for 20 min on ice followed by sonication for 30 s. The cell lysate was then vortexed for 30 s and centrifuged at 20,000×g for 10 min at 4 °C. Protein concentrations were measured using a Bio-Rad Benchmark Plus Spectrophotometer system and all samples were standardized to 70 µg/mL for subsequent immunoassays.
Analysis of serum cytokines was performed using a multiplex mouse 25-plex bead immunoassay (Milliplex Mouse Cytokine/Chemokine Magnetic Bead Panel #MCYTMAG70PMX25BK). The following cytokines were quantified: G-CSF, GM-CSF, IFN-γ, IL-1α, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-7, IL-9, IL-10, IL-12 (p40), IL-12 (p70), IL-13, IL-15, IL-17, IP-10, KC, MCP-1, MIP-1α, MIP-1β, MIP-2, RANTES, and TNF-α. Standards and reagents were all prepared according to the manufacturers’ recommendations. Each serum and brain sample were diluted to a standardized concentration and run in duplicates. 25 uL of sample, standards, or blanks was loaded into a 96-well plate with appropriate amounts of assay buffer and matrix solution. The plate was then incubated overnight with antibody-coupled magnetic beads. The following day, after a series of washes, the plate was incubated with a biotinylated detection antibody on a shaker for 1 h. Streptavidin–phycoerythrin was added and incubated while shaking continued for 30 min. Plates were washed using a Bio-Plex handheld magnet (Bio-Rad Laboratories, Hercules, CA, USA). After the final wash, the plate was analyzed using a Bio-Rad Bio-Plex 200 plate reader (Bio-Rad Laboratories, Hercules, CA, USA) and analyzed using Bio-Plex Manager software (Bio-Rad Laboratories, Hercules, CA, USA). The following were the minimal amounts of detectable cytokine concentrations: G-CSF: 1.7 pg/mL; GM-CSF: 10.9 pg/mL; IFNγ: 1.1 pg/mL; IL-1α: 10.3 pg/mL; IL-1β: 5.4 pg/mL; IL-2: 1.0 pg/mL; IL-4: 0.4 pg/mL; IL-5: 1.0 pg/mL; IL-6: 1.1 pg/mL; IL-7: 1.4 pg/mL; IL-9: 17.3 pg/mL; IL-10: 2.0 pg/mL; IL-12 (p40): 3.9 pg/mL; IL-12 (p70): 4.8 pg/mL; IL-13: 7.8 pg/mL; IL-15: 7.4 pg/mL; IL-17: 0.5 pg/mL; IP-10: 0.8 pg/mL; KC: 2.3 pg/mL; MCP-1: 6.7 pg/mL; MIP-1α: 7.7 pg/mL; MIP-1β: 11.9 pg/mL; MIP-2: 30.6 pg/mL; RANTES: 2.7 pg/mL; TNF-α: 2.3 pg/mL. As per previous studies [47,48,49], sample concentrations that fell below minimal detection value were given a proxy value of half the limit of detection for statistical comparisons using non-parametric Mann–Whitney U analyses and BH-correction for multiple comparisons.
Immunohistochemistry
30 um brain sections were stained for immunofluorescent analysis with a commonly used microglia marker, purinergic receptor P2Y12 (P2ry12), to analyze microglial morphology. Brain sections used for the immunofluorescent images presented in Fig. S9 were co-stained with microglia marker ionized calcium-binding adaptor molecule 1 (Iba1) and blood vessel marker platelet endothelial cell adhesion molecule (CD-31/PECAM-1). Free-floating brain sections were washed 3 times for 5 min each in 1xPBS, permeabilized in 1xPBS and 0.5% Triton (Fisher BioReagents BP151-500) for 5 min, and incubated in blocking solution made of 1xPBS, 0.03% Triton, and 1% Bovine Serum Albumin (BSA; Bio-techne Tocris 5217) for 1 h. After the blocking steps, sections were incubated overnight at 4 °C in primary antibody solution containing 2% Normal Donkey Serum (NDS; Jackson Immunoresearch Laboratories Inc. 017-000-121), 1xPBS, 0.03% Triton, and primary antibody (rabbit anti-P2RY12: 1:500, Anaspec AS-55043A or guinea pig anti-Iba1: 1:1000, Synaptic Systems HS-234-308 and mouse anti-CD31/PECAM-1: 1:100, Invitrogen 14-0311-82). After primary antibody incubation, sections were washed 3 times for 5 min each with 1xPBS and 0.03% Triton before incubating for 2 h in secondary solution containing 2% NDS, DAPI (1:1000; Biolegend 422801), and secondary antibodies (Alexafluor 647 donkey anti-rabbit: 1:500, Jackson Immunoresearch Laboratories Inc. 711-606-152 or Alexafluor 647 donkey anti-guinea pig: 1:500, Jackson Immunoresearch Laboratories Inc. 706-605-148 and Alexafluor 488 anti-mouse: 1:500, Jackson Immunoresearch Laboratories Inc. 715-546-150). Sections were washed 3 times for 5 min each with 1xPBS and 0.03% Triton and 1 time for 5 min with 1xPBS before being transferred into 1xPBS for storage before mounting. Sections were mounted onto microscope slides (Premium Superfrost Plus Microscrope Slides, VWR CA48311-703) and air dried before being coverslipped with mounting media (EverBrite Fluorescence Antifade Mounting Media, Biotium 23001).
Imaging and microglia morphology analysis
All mounted brain sections used for morphology analysis were imaged on the ZEISS Axioscan 6 microscope slide scanner at 20 × magnification with a step-size of 1um using the z-stack acquisition parameters within the imaging software (ZEISS ZEN 3.7). Mounted brain sections used for microglia and blood vessel interaction visualization in Fig. S9 were imaged on the ZEISS Axioscan 7 microscope slide scanner at 40 × magnification with a step-size of 0.5um using the z-stack acquisition parameters within the imaging software (ZEISS ZEN 3.9). During acquisition for all images, Extended Depth of Focus (EDF) images were created using maximum projection settings within the ZEN software and saved as the outputs in the final.czi files. Maximum projection EDF images only compile the pixels of highest intensity at any given position in a z-stack to construct a new 2D image which retains the 3D information. Using ImageJ, we created.tiff images of each fluorescent channel from.czi files and selected and saved.tiffs of brain regions of interest (ROIs) to use as input for downstream morphological analysis in MicrogliaMorphology and MicrogliaMorphologyR, as detailed in Kim et al., 2023 [50]. The parameters used for MicrogliaMorphology included: auto local thresholding using the Mean method and radius 100 with an area filter of 0.01613–0.1808 in^2 (or 153.3358–1718.73 um^2 with.czi conversion applied factor of 0.325 um/pixel applied to scale). After performing principal components analysis on the 27 morphology measures output from MicrogliaMorphology, the first 3 principal components were used as input for K-means clustering in MicrogliaMorphologyR to define four distinct classes of microglia morphology (ameboid, hypertrophic, rod-like, ramified). We focused our analyses on multiple brain regions including the hippocampus, frontal cortex, and striatum. Images of coronal brain sections containing these regions were aligned to the Allen Brain Atlas (mouse.brain-map.org) [51] using the ImageJ macro FASTMAP [52]. We assessed how percentages of morphological clusters for each mouse changed with LPS treatment by fitting a generalized linear mixed model to measure percentage changes as a factor of cluster and treatment interactions for each brain region and sex with mouseID as a repeated measure (percentage ~ ClusterID*Treatment + (1|MouseID)) using MicrogliaMorphologyR, which calls to the glmmTMB R package [50, 53]. We had to analyze each sex and brain region separately because the males had an additional treatment group (3xLPS) that was not present in the females and the model converged when considering ClusterID*Treatment*BrainRegion interactions together. The model was assessed using ANOVA and tests between groups were corrected across for comparisons made within each cluster (~ Treatment|Cluster) using the BH method. Corrected p-values < 0.05 were considered statistically significant. 3D reconstructions of examples of microglia labeled as each of the different morphological classes (rod-like, ramified, ameboid, and hypertrophic) were generated in ImageJ software using the 3D viewer plugin [54].
Events of microglia aligning end-to-end were manually quantified within each brain region of interest (ROI) for every image using ImageJ software. Manual counts were normalized to the areas of the ROIs to calculate densities of event observations for each brain region analyzed. Microglial density was quantified by dividing the number of microglia analyzed in each ROI by the area of each ROI. We assessed how densities of events of microglial alignment behavior and microglial densities change with LPS treatment across brain regions by fitting separate linear mixed effects models to measure either alignment behavior event density or microglial density as a factor of treatment and brain region interactions with mouseID as a repeated measure, for each sex separately (AlignmentEventDensity ~ Treatment*BrainRegion + (1|MouseID) and MicroglialDensity ~ Treatment*BrainRegion + (1|MouseID)). Both models were assessed using ANOVA and tests between groups were corrected across for comparisons made within each brain region (~ Treatment|BrainRegion) using the BH method. Corrected p-values < 0.05 were considered statistically significant.
Microglia isolations
Isolated cortical tissue from each mouse was placed on individual petri dishes on ice with 1 mL of digestion buffer, made up of 10 ul DNase I (STEMCELL Technologies; Cat. 07900) + 990 ul activated Papain solution (Worthington Biochemical, Cat. LS003124), added to each brain sample before it was minced into small (< 1 mm) pieces using a scalpel. Minced samples in digestion solution were transferred to individual wells within a 24-well cell culture plate and incubated on ice for 30 min. Digested brain solution was transferred into a 15 mL glass dounce homogenizer on ice with 5 mL ice-cold flow cytometry (FACS) buffer comprising 2.5% bovine serum albumin (BSA; Tocris Bioscience Cat. 5217/100G) and 1 M ethylenediaminetetraacetic acid (EDTA; Invitrogen Cat. 17892) in Hanks’ Balanced Salt Solution (HBSS). Brains were slowly and gently homogenized 20 times using a loose pestle, followed by 2 times using a tight pestle. Brain homogenate was strained through a 70um nylon mesh strainer (Thermo Fisher, Cat. 08-771-1) with FACS buffer into 50 mL tubes on ice. 50 mL tubes were centrifuged at 300×g for 10 min at 4C, at the machine’s minimum brake setting. After carefully and discarding the supernatant, cell pellets were gently resuspended in 30% Percoll solution and transferred to new 15 mL tubes, where 2 mL of 37% Percoll solution was gently underlaid. To make the Percoll solutions, 3.6 mL of Isotonic Percoll (Millipore Sigma, Cat. GE17-0891-02) was added to 0.4 mL 10xHBSS per brain to make a 100% Percoll solution, which was diluted down to 30% Percoll and 37% Percoll in 1xHBSS. 15 mL tubes containing the layered Percoll gradients were centrifuged at 700×g for 10 min at 4C, at the machine’s minimum brake setting. Myelin layers and supernatant layers were discarded using a transfer pipette, leaving 1 mL of the 37% layer at the bottom of the tube. Cell pellets were gently resuspended and topped up to 10 mL with ice-cold FACS buffer then centrifuged at 300×g for 10 min at 4C, at the machine’s minimum brake setting. Supernatant in tubes was discarded down to 200 ul, 5 mL of ice-cold FACS buffer added, and tubes spun again at 300×g for 10 min at 4C at the machine’s minimum brake setting. Finally, supernatant was removed and cells resuspended in 200 uL remaining FACS buffer before being counted on the hemocytometer. After counting, samples were diluted to a concentration of 2.5 × 107 cells/mL in 0.05–0.2 mL of FACS buffer and transferred to a round-bottom, non-adherent 96-well plate (Corning, Cat. 3788). CD11b+ microglial cells were isolated by magnetic bead enrichment using the EasySep™ Mouse CD11b Positive Selection Kit (STEMCELL Technologies, Cat. 18970), according to the manufacturer’s protocol.
RT-qPCR
Isolated brain region tissues and isolated microglia cells were homogenized using QIAshredder columns (Qiagen 79656). RNA from homogenized tissue was isolated using the Monarch Total RNA Miniprep kit, according to the manufacturer’s protocol (New England Biolabs T2010S). For isolated microglia, cells were homogenized and RNA extracted using the Total RNA Purification from Cultured Mammalian Cells version of the protocol for the Monarch Total RNA Miniprep kit. Isolated RNA from brain regions and purified microglia was converted to cDNA and prepared for Real Time quantitative PCR (RT-qPCR) reactions using the LunaScript® RT SuperMix Kit (New England Biolabs E3010) and Luna® Universal qPCR Master Mix (New England Biolabs M3003), according to the manufacturer’s Protocol for Two-step RT-qPCR. The primers in Table 1 were used for RT-qPCR reactions and RT-qPCR reactions were loaded into MicroAmp™ Fast Optical 96-Well Reaction plates (Thermo Fisher Scientific 43-469-07) sealed with optical adhesive film (Thermo Fisher Scientific AB-1170). All biological samples were tested in triplicates. Reactions were run on the QuantStudio 6 qPCR machine (Life Technologies, California, USA) and samples were analyzed by the comparative CT method (ΔΔCT). The ΔCT was calculated by subtracting the average CT value of the housekeeping gene hypoxanthine phosphoribosyltransferase 1 (Hprt1) from the average CT value from the gene of interest. ΔΔCT values were calculated by subtracting the ΔCT from the gene of interest by the average ΔCT from its respective control. Fold change was calculated by using the equation 2−ΔΔCT. We assessed how expression of each target gene changes with LPS treatment by fitting a linear mixed-effects model to measure fold change values as a factor of LPS Treatment for each gene independently (gene fold change ~ Treatment + (1|MouseID)) using the nlme R package. We analyzed each sex separately, as the males had an additional treatment group, 3xLPS, that was not present in the females. The model was assessed using ANOVA and tests between LPS treatments were corrected for multiple comparisons using the BH method. Corrected p-values < 0.05 were considered statistically significant.
Results
Repeated LPS produces robust training and tolerance of acute sickness behaviors, but not long-term behavioral changes
In our experimental model, mice were given daily intraperitoneal injections of low-dose LPS (0.5 mg/kg) or vehicle (1xPBS) over a span of four days and sickness behavior was scored 3 h following the final treatment of PBS, 1xLPS, 2xLPS, 3xLPS or 4xLPS (Fig. 1A, B). Repeated LPS injections produced both training and tolerance of sickness behaviors: sickness behavior peaked and the greatest drop in body weight was observed across treatment groups after the second injection of LPS relative to PBS (Figs. 1B, 2A; BH-adjusted p < 0.05, Supp. Table S2). Sickness behavior showed recovery towards baseline levels after the third and fourth LPS injections (Figs. 1B, 2A). While body weights recovered by 2 days after peak weight loss from the second exposure to LPS across treatment groups, weights in LPS-treated groups were still significantly lower than that of the PBS group for up to 6 days after peak loss (BH-adjusted p < 0.05, Supp. Table S2) and none of the mice in the LPS-treated groups fully recovered to their pre-treatment body weights prior to the start of the behavioral experiment (Fig. 2A; BH-adjusted p < 0.05, Supp. Table S2). However, there were no long-term differences in anxiety-like, depressive-like, learning and memory, nor repetitive behaviors across treatment groups when behaviors were then assessed, indicating that all treatment conditions had fully recovered and did not have long-term detrimental impacts of repeated LPS (Fig. 2B; BH-adjusted p < 0.05, Supp. File S1).
Cytokine and chemokine protein levels increase first in serum and then in brain
To examine protein levels of known inflammatory proteins, including chemokines, cytokines, and other signaling molecules, we performed a Luminex multiplex protein quantification on serum and tissue collected from four different brain regions (frontal cortex, hippocampus, striatum, cerebellum) 3 h after final injections. We identified significant changes in protein levels of both pro- and anti-inflammatory molecules (Fig. 1C; BH-adjusted p < 0.05, Supp. Table S1). In serum, we found significant increases of G-CSF, IL-6, IP-10, KC, MCP-1, MIP-1a, MIP-1b and RANTES in response to 1xLPS that was maintained in response to 2xLPS and then returned to baseline level in most cases after 3xLPS and 4xLPS. In frontal cortex we found significant increases in G-CSF, IL-6, and IP-10 to 1xLPS, but the majority of changes were in response to the 2xLPS (10 proteins), many of which again were no longer significant after 3xLPS and 4xLPS. Similar patterns were observed in hippocampus, striatum, and cerebellum with the majority of significant increases occurring after 2xLPS (Fig. 1C; BH-adjusted p < 0.05, Supp. Table S1). We further explored these findings in males and females using RT-qPCR on RNA isolated from frontal cortex tissue for Il-1β and Tnf-α, two hallmark pro-inflammatory cytokines that are known to be produced by microglia during immune activation. Similar to the patterns observed in the protein analysis, we found the greatest induction of Il-1β and Tnf-α in response to 2xLPS, which returned towards baseline levels after 3xLPS and 4xLPS (Fig. 1D; BH-adjusted p < 0.05, Supp. Table S5). Together, these findings suggest that 2xLPS produces the largest increase in inflammatory signaling molecules and that these responses are blunted after 3xLPS and 4xLPS. These patterns also parallel the peak in sickness behaviors and drop in body weight observed after 2xLPS that resolve towards baseline levels by 4xLPS.
We also measured Il-1β and Tnf-α expression 3 h vs. 24 h after 1xLPS and 2xLPS (Fig. S1A; BH-adjusted p < 0.05, Supp. Table S5). Cytokine expression significantly decreased towards baseline levels 24 h after 1xLPS, but subsequently increased above previous levels 3 h after 2xLPS, indicating that a single hit of LPS primes a heightened response to a subsequent LPS injection (immune training). Cytokine expression significantly decreased towards baseline levels 24 h after 2xLPS and the subsequent response to 3xLPS and 4xLPS were significantly blunted, indicating that immune tolerance occurred. Furthermore, we assessed how Il-1β and Tnf-α expression changed 3 h, 24 h, 48 h, and 7 days after 1xLPS in both males and females to profile cytokine dynamics over time in response to a single dose of LPS (Fig. S1B; BH-adjusted p < 0.05, Supp. Table S5). Cytokine expression peaked acutely at 3 h after 1xLPS, significantly decreased by 24 h, and reached baseline levels of expression by 48 h or 7 days. Taken together, our findings demonstrate that the differences in gene expression observed are linked to the number of LPS injections and not a byproduct of the timing of injections.
Repeated LPS produces training and tolerance of gene expression within brain regions that is likely microglia-driven
To examine changes in gene expression in response to repeated LPS stimulation, we performed RNA-sequencing on brain tissue from the hippocampus, striatum and frontal cortex 3 h after the last LPS or PBS injection (n = 4/treatment/region; Fig. 1A, Fig. S1). Differentially expressed genes were identified across brain regions for differences in expression between the PBS controls and each LPS condition (BH-adjusted p < 0.05 and fold change > ± 2, Supp. Table S3) and clustered according to their gene expression patterns across brain regions and treatments, which were the greatest explanatory variables in the data as assessed by MDS analysis (Fig. 3A, Fig. S2B). Similar to the impacts on brain cytokine protein levels, sickness behaviors, and body weight changes, changes in gene expression were most striking following 2xLPS (Fig. S3A–C). Very few genes were decreased in response to LPS across all comparisons and brain regions. Clusters that emerged from patterns of gene expression were characterized as “2xLPS-sensitive”, “4xLPS-sensitive”, and “LPS-decreased”. Across brain regions, 2xLPS-sensitive genes (n = 226, 52% of DEGs) were uniquely increased in response to 2xLPS, 4xLPS-sensitive genes (n = 120, 28% of DEGs) were uniquely increased in response to 4xLPS, and LPS-decreased genes (n = 86, 20% of DEGs) were decreased from baseline levels across LPS treatments (Fig. 3A).
The different DEG clusters were enriched for shared and unique GO terms (Fig. 3B; BH-corrected p < 0.05, Supp. Table S3). The majority of 2xLPS-sensitive genes were not engaged at 3xLPS and 4xLPS, suggesting the formation of tolerance for the majority of these genes. 2xLPS-senstive genes were enriched for biological processes related to immune regulation including response to virus, response to interferon-beta, positive regulation of defense response, defense response to virus, defense response to symbiont, positive regulation of tumor necrosis factor production, acute inflammatory response, and myeloid leukocyte activation. These findings are consistent with previous work examining LPS-activated gene expression in the brain and periphery and indicate that alterations in immune cell signaling underlie the enhancement in gene expression to 2xLPS and repression to 3xLPS and 4xLPS in the brain [4, 19, 58, 59]. The 4xLPS-sensitive genes were enriched for many of the same processes, as well as glial cell development and astrocyte development. LPS-decreased genes were enriched for processes related to neuronal regulation including positive regulation of cell junction assembly, response to cocaine, axon guidance, neuron projection guidance, and positive regulation of monoatomic ion transport. Together, findings suggest that separate gene expression signatures are engaged in response to repeated injections of LPS.
To identify whether the differential gene expression observed in our bulk RNAseq data could be driven by microglia, we performed enrichment analysis of our DEGs against curated lists in MGEnrichment, our recently developed database of published microglial gene lists [36]. We identified shared and unique enrichment between 2 and 4xLPS-sensitive but not LPS-decreased genes with published microglia-specific gene lists important for maintaining microglia homeostasis, disease-associated microglia (DAM) phenotypes, immune primed microglia, microglia development, and microglial changes in microbiome perturbations, among many others (Fig. S4; BH-corrected p < 0.01, Supp. Table S3). We also analyzed publicly available single cell atlas datasets [44, 45] (Fig. S5) to explore which brain cell types highly express the DEGs identified in our study under homeostatic conditions with the aim of understanding which cell types could potentially be driving the observed impacts of LPS on gene expression. Across both datasets, microglia, immune cells (perivascular macrophages), and endothelial cell types expressed 2xLPS-sensitive and 4xLPS-sensitive genes most highly, while neuronal cell types expressed LPS-decreased genes most highly. These findings are consistent with the biological functions identified by gene ontology analysis on these clusters (Fig. 3B). Furthermore, two LPS-sensitive genes from our DEG list that are highly expressed in microglia (Fig. S6B, E), Cxcl16 and S100a9, also displayed the same training and tolerance gene expression patterns when assessed by RT-qPCR in repeat experiments done in frontal cortex in both males and females, as well as in microglia isolated from cortical regions in males (Fig. S6C, F; BH-corrected p < 0.05, Supp. Table S5). Together, findings suggest microglia as a likely candidate cell type driving the differential expression observed in the 2xLPS-sensitive and 4xLPS-sensitive gene clusters.
Repeated LPS regulates immune memory through engagement with IRF and NFkB family transcription factors
To examine which upstream regulatory factors could be driving gene expression in response to repeated LPS, we examined promoter regions of DEGs for enrichment of transcription factor binding motifs (Fig. 4A,C). Promoters of 2xLPS-sensitive cluster DEGs were significantly enriched for interferon regulatory factor (IRF), basic helix-loop-helix/bZIP (bHLH/bZIP), Rel homology domain (RHD), C2H2 zinc finger (Zf), and signal transducer and activator of transcription (Stat) transcription family factors across de novo and known Homer motifs. (Supp. File S2) Different sets of genes within the 2xLPS-sensitive DEG cluster had differences in the presence of binding motifs for transcription factors across the enriched transcription factor families (Fig. 4D) in their promoters. Notably, 2xLPS-sensitive genes that increased in response to 1xLPS and stayed increase after 2xLPS did not have binding sites for bZIP family transcription factors, which were present throughout most of the other 2xLPS cluster genes. Moreover, only subsets of 2xLPS-sensitive genes had binding sites for IRF family transcription factors. Binding sites for RHD family transcription factor NFkB, another critical regulator of immune genes, were found in most 2xLPS cluster DEG promoters. DEG promoters were also significantly enriched for direct binding sites of RHD family (NFkB, p65) and IRF family (IRF8) from publicly available ChIP-seq data [38,39,40,41,42], further supporting the roles of these transcription factors as major regulators of LPS-induced gene expression differences in immune training and tolerance (Fig. 4B; BH-adjusted p < 0.05, Supp. Table S3). DEG promoters were also enriched for H3K27ac ChIP-seq peaks that were lost or gained with conditional IRF8 knockouts in microglia in adults [38], suggesting that IRF8 may be interacting with other transcriptional machinery such as histone acetyltransferases and deacetylases to mediate microglial gene expression at these sites.
Microglia change their morphology across brain regions in the formation of immune training and tolerance
Using MicrogliaMorphology and MicrogliaMorphologyR, as previously described [50], we quantified shifts in morphological populations within the same brain regions assessed for RNA-seq (frontal cortex, hippocampus, striatum) from male and female mice given 1-4xLPS (Fig. 5A, S7A–E). Using unbiased clustering, we identified four clusters of microglial morphologies: ramified, rod-like, hypertrophic, and ameboid. While the percentage of ramified and rod-like microglia mostly remained consistent within the different brain regions and sexes with LPS treatment, there was an LPS-mediated shift between ameboid and hypertrophic states (Fig. 5A, B, S8A-B; BH-corrected p < 0.05, Supp. Table S4). Within males, the hippocampus and striatum showed a significant shift from ameboid to hypertrophic morphological states in the formation of immune priming (2xLPS) and back with immune tolerance (4xLPS). In males the frontal cortex maintained the 2xLPS-induced hypertrophic morphology with immune tolerance. Within females, the frontal cortex, hippocampus, and striatum showed a morphological switch to a hypertrophic morphology which was sustained with immune tolerance.
Furthermore, we also observed instances of microglia aligning end-to-end with increasing exposures to LPS. Within males, the occurrence of microglia aligning end-to-end increased with 2xLPS, began to decrease after 3xLPS, and re-emerged with 4xLPS across the frontal cortex, hippocampus, and striatum (Fig. 6A, B; BH-corrected p < 0.05, Supp. Table S4). Within females, occurrence of this alignment behavior significantly increased with 4xLPS in the frontal cortex and hippocampus, but not in the striatum. There were no significant differences in microglial density with repeated LPS across brain regions and in both sexes, indicating that the changes observed in occurrences of microglia aligning end-to-end were not due to changes in cell numbers (Fig. S7F; BH-corrected p < 0.05, Supp. Table S4). Upon further investigation by histology, we found that microglia exhibiting this alignment behavior in our LPS model were wrapped around blood vessels across all three brain regions (Fig. S9A–C).
Discussion
Disruption of microglial activity in the brain has been well-documented in the pathology of a wide range of brain disorders, but the molecular mechanisms underlying how innate immune memory directly contributes to the onset and progression of brain disorders remains understudied. In our study, we used a recently described mouse model of repeated LPS [19] to investigate the immediate transcriptomic and microglia morphological changes that contribute to the formation of immune memory in the frontal cortex, hippocampus, and striatum, as well as the long-term effects of these changes on various behaviors relevant to brain disorders. Mice showed robust changes in acute sickness behavior, body weight, and pro-inflammatory cytokine expression with LPS-induced immune training, which resolved with the formation of immune tolerance. We found transcriptomic changes in a subset of genes that showed a training response after 2xLPS and a tolerized response at 3xLPS and 4xLPS. These “2xLPS-sensitive” genes were enriched for biological processes related to immune regulation and were targeted by IRF and RHD family (NFkB) transcription factors, two key families whose coordinated activity has previously been linked to the formation of innate immune memories in microglia and other immune cells [19, 58, 60,61,62,63,64,65]. Together with the characterized LPS-induced shifts in morphology, our work identifies several key changes in microglial regulation that occur during the formation of innate immune memories that inform previous findings on innate immune memory in disease contexts in the brain.
Microglia are the only tissue-resident macrophages in the brain parenchyma [66], self-renew without contribution from the periphery [17], and exhibit different properties morphologically and functionally when placed outside of the brain environment [67], suggesting that brain-intrinsic signals that drive microglial identity make them distinct from macrophages in other tissues. However, the majority of what is known about regulation of immune training and tolerance comes from work in cultured peripheral macrophages. Work in in vitro macrophage culture models has demonstrated a clear role for epigenetic regulation of gene expression programs mediating training and tolerance [68, 69]. Recent work has demonstrated that training and tolerance occur in microglia in response to repeated peripheral injections of LPS [19, 70]. Wendeln et al. [19] demonstrated alterations in microglial gene expression and enhancer activity 6 months after LPS injection in the Alzheimer’s disease model, suggesting long-term reprogramming of microglial function after repeated LPS. Zhang et al. [58] found that enrichment of permissive epigenetic marks at accessible enhancer regions could explain training of microglia to LPS and that tolerized responses were regulated by the loss of permissive epigenetic marks at these sites. In line with both studies, our study identified LPS-sensitive DEGs that were enriched for Interferon (IRF) family and NFkB transcription factor motifs in their promoters. IRFs are induced in response to interferon (IFN), LPS, and other immune stimuli, and initiate gene expression by forming regulatory complexes with IRF family members and other transcription factors, such as PU.1, STAT1, and NFkB, which further bind with the interferon-sensitive response elements (ISRE) to regulate IFN and TLR signaling pathways [63, 71]. IRF family transcription factors have been shown to be required for Il-1β expression in microglia after peripheral nerve injury [72] and regulate disease-associated microglial gene expression of markers including Apoe and Trem2 [73]. IRF8 is particularly well-characterized in microglia, and has been shown to bind to and regulate enhancer regions of microglia during postnatal development to contribute to the establishment of microglial identity and function [38]. IRF8-deficient microglia have documented changes in morphology and function marked by reduced Iba-1 expression, less elaborated processes, decreased phagocytic capacity, reduced expression of pro-inflammatory genes, and less proliferative potential in mixed glial cultures [65, 74,75,76]. Accumulation of p50/RelA, subunits of NFkB, in the nucleus has also been tied to increased transcription of pro-inflammatory cytokines including iNOS, Tnfα, Il-1β, and Il-6 in the formation of pro-inflammatory phenotypes in microglia [64, 65, 77]. We identified enrichment for direct binding sites of NfkB/p65 and IRF8 from published microglial and bone marrow-derived macrophage studies [38, 40, 41] within our DEGs, together pointing to the coordinated activity of these transcription factors as key regulators of immune priming and tolerance in the brain.
Although LPS given peripherally reliably induces a neuroinflammatory response in the brain, LPS has been shown to be largely unable to penetrate the blood brain barrier (BBB) in healthy animals, unless at very high doses [78]. Instead, impacts of peripherally administered LPS are likely conveyed to the brain parenchyma by TLR4 receptors on meningeal and endothelial cells, which then signal to other cell types in the brain including microglia [13, 79]. Microglia are likely cooperating with other cell types such as astrocytes to perpetuate the expression of pro-inflammatory cytokines, as microglial depletion alone fails to alleviate LPS-induced sickness behaviors [13, 80]. Once immune signaling molecules reach the brain parenchyma, they can bind to cytokine receptors expressed on both neurons and glial cells to directly or indirectly modulate neuronal firing properties, which could ultimately shape circuit function and downstream behavior [13]. In our analysis of publicly available single-cell brain atlas datasets [44, 45], we found that neuronal cell types expressed LPS-decreased genes most highly (Fig. S5), which were enriched for gene ontology terms related to neuronal regulation and function (Fig. 3B). Microglia; perivascular macrophages, which exhibited high expression of border-associated macrophage (BAM) marker gene Mrc1: http://mousebrain.org/genes/Mrc1.html); and endothelial cell types expressed 2xLPS-sensitive and 4xLPS-sensitive genes most highly (Fig. S5), which were enriched for gene ontologies related to immune functions. While our analysis hints at candidate brain cell types that could be driving the differences in gene expression observed in our tissue-level RNA-seq data, the atlas datasets used for analysis were derived from mice in baseline, homeostatic conditions. As such, cell type contributions that are LPS-specific could be missed in this analysis. For example, the 4xLPS-sensitive DEG cluster identified in our study was uniquely enriched for GO terms related to astrocyte development (Fig. 3B), indicating a potential role for astrocytes in our model of immune memory. Astrocytes have been shown to exhibit LPS-induced immune training and tolerance [81,82,83,84,85,86] and respond to signals from microglia after LPS exposure [87]. Future cell type-specific experiments such as single-cell sequencing in the same repeated LPS model could be used to further discriminate the contributions of microglia vs. other immune cell types such as BAMs, as well as roles of other brain cell types including neurons and astrocytes to the gene expression differences that we observed.
A recent study using the same repeated low-dose LPS model [88] defined a temporal sequence after 2xLPS by which changes in microglia morphology, density, and expression of phagocytic markers precedes GABAergic synapse loss, followed by memory impairment on the novel object recognition task. However, this study focused on short-term changes within a range of 6 days post-2xLPS, and investigation into the prolonged effects of repeated low-dose LPS on behaviors long after recovery had yet to be thoroughly examined. We failed to identify long-term effects of repeated LPS on a battery of behaviors relevant to brain disorders. Mice showed immediate sickness behaviors including lethargy, ptosis, and piloerection 3 h after 2xLPS exposure, but they did not exhibit any long-term deficits in learning and memory tasks or anxiety-like, depressive-like, and repetitive behaviors when assessed weeks after LPS treatment. While our dataset was too underpowered to separate the statistical analysis out by sex, future experiments could focus on disentangling sex-specific differences in these behaviors. Although our data suggest that our LPS model does not contribute to any obvious observable long-term behavioral deficits in healthy mice, preconditioning with repeated peripheral injections of low-dose LPS has been demonstrated to alleviate disease pathology that develops months later in mouse models of Alzheimer’s Disease and neuroinflammatory responses in stroke and traumatic brain injury [19, 20, 89,90,91]. This suggests that a secondary disease pathology or injury is necessary to reveal the long-term reprogramming effects of LPS on behavior.
Changes in pro-inflammatory gene expression with immune priming and tolerance in the brain have been well-documented to be accompanied by changes in microglial morphology [19, 92, 93]. 2xLPS-sensitive and 4xLPS-sensitive DEGs identified in our study were enriched for DEGs previously identified as dynamically expressed in the transition from embryonic to postnatal mouse microglia development in both sexes (Fig. S4; BH-corrected p < 0.01, Supp. Table S3) [94]. During embryonic development, immature microglia display an ameboid morphology which aids in phagocytosis of apoptotic debris, synaptic pruning of developing circuits, and increased motility as microglia migrate along developing white matter tracts to populate gray matter. As development progresses postnatally, microglia morphology transforms into a fully ramified, mature form that is maintained throughout life [95,96,97]. Enrichment of our LPS innate immune memory genes for genes with differential expression between the ameboid and mature ramified development states, indicates that the morphological adaptations in the adult brain to repeated LPS exposure could involve shared developmental genes and pathways.
While the proportions of ramified and rod-like microglial forms remained mostly consistent across LPS treatments, we observed a switch between ameboid and hypertrophic morphology states in the formation of immune memory that was observed in most brain regions within both sexes. Hypertrophic morphology is described by enlarged cell soma with shorter, thicker, hyper-ramified processes [21, 98], which give microglia a “bushy” appearance. This morphological form is commonly observed in regions of brain pathology in human cases of Alzheimer’s Disease and Huntington’s Disease, as well as in mouse models of Alzheimer’s Disease, stroke, accelerated aging, chronic stress, depression, and traumatic brain injury [21]. Hypertrophic microglia have been described in close vicinity to beta-amyloid plaques in Alzheimer’s Disease in both human tissue and mouse models. Hypertrophic microglia also exhibit increased surface markers for the disease-associated microglia (DAM) phenotype including CD-33 and triggering receptor expressed on myeloid cells 2 (TREM2), a key switch in the change to the DAM phenotype [21, 96, 99, 100]. TREM2 has recently been described in the context of the same repeated LPS model used in our study to mediate decreases in microglial phagocytosis of synapses and morphology changes with the formation of immune tolerance in the hippocampus [93]. TREM2 knockout prevented acute increases in both microglial soma size and changes in process length and arborization that aligned with a hypertrophic morphology [93]. While we did not identify Trem2 as a significant DEG in our bulk tissue analysis, we did observe enrichment of 2xLPS-sensitive DEGs with published microglia-specific gene lists for the DAM phenotype (Fig. S4; BH-corrected p < 0.01, Supp. Table S3), pointing to the significance of the hypertrophic microglia morphology in regulating immune memory states in the brain.
We also observed instances of microglia aligning end-to-end along blood vessels with increasing exposures to LPS (Fig. 6, Fig. S9). Spatially distinct subsets of microglia in the brain have been found to interact closely and stably over time with the vasculature to migrate within the developing brain parenchyma [101]; and to regulate capillary diameter, blood flow, and vasodilator responses [102, 103]. Microglia associated with the vasculature have also been found to change in numbers in mouse models of ischemic cortical stroke [104], Alzheimer’s disease [103], and multiple sclerosis [105]. In line with our findings, systemic inflammation induced by a single intraperitoneal injection of LPS has been shown to cause microglia to migrate towards cerebral vessels and the proportion of microglia in contact with blood vessels increases with repeated exposures to LPS [106]. Furthermore, the alignment behavior we observed has been documented previously in the literature where rod-like microglia have been shown to align end-to-end with one another adjacent to neuronal processes in diffuse traumatic brain injury (TBI) in rats [107,108,109]. These “trains” of rod-like microglia were shown to emerge and resolve dynamically with the damage and repair cascades that evolve over the first week post-TBI, where rod trains accumulated at 2 and 6 h post-injury, with a decrease in percentage microglia population at 1 and 2 days post-injury, followed by a subsequent re-accumulation at 7 days post-injury [107]. Similar to these documented patterns, the occurrence of microglial aligning end-to-end in our study increased with 2xLPS, began to decrease towards baseline levels at 3xLPS, and re-emerged with 4xLPS (Fig. 6A), suggesting that microglia may be engaged in separate, delayed waves of repair after inflammation in the brain in the switch from an immune trained to tolerized state after the resolution of cytokine responses and obvious sickness behaviors. This could explain the short-term learning and memory impairments described by Jung et al. [88], which were observed at 6 days after LPS, and why we failed to observe them in our study in which behavioral testing was performed weeks after LPS. Future work to further characterize microglial interactions with the brain vasculature in our repeated LPS model - such as their influence on blood vessel diameter and the temporal dynamics of these interactions following LPS exposure - would lend deeper insight into the significance of the LPS-induced microglia alignment behavior that we observed and its relevance to immediate and delayed repair processes following inflammation.
Together, findings from our study offer new insight into microglial changes in the formation of immune memory. Given that microglia are long-lived cells (4–20 years in humans) [110] that self-renew without contribution from the periphery [17, 18] and epigenetically tune their function to the local brain environment through differential enhancer activity [111, 112], the acute changes induced during the formation of innate immune memories could have long lasting impacts on the brain’s response to disease and injury. Future work can link our identified acute changes to the long-lasting impacts identified in other model systems to identify how innate immune memories are maintained in microglia and impact disease risk. Understanding these underlying mechanisms will allow for future development of therapies to either enhance or blunt inappropriate microglial activity in the context of brain diseases.
Availability of data and materials
The datasets generated and/or analyzed during the current study are freely available to download and all relevant materials including input images, data, and supporting code used in this study are available on the Open Science Framework (OSF) website at https://osf.io/gb7a9/ and the Ciernia Lab Github at https://github.com/ciernialab/Kim2024_1-4xLPS. The RNA sequencing data generated in this study are available through the NCBI GEO database under accession code GSE165287.
Abbreviations
- LPS:
-
Lipopolysaccharide
- PBS:
-
Phosphate buffered saline
- TLR4:
-
Toll-like receptor 4
- MDS:
-
Multi-dimensional scaling
- CPM:
-
Counts per million
- DEG:
-
Differentially expressed gene
- TMM:
-
Trimmed mean of m-values
- GO:
-
Gene ontology
- BMDM:
-
Bone marrow-derived macrophage
- P2ry12:
-
Purinergic receptor P2Y12
- Iba1:
-
Ionized calcium-binding adaptor molecule 1
- CD-31/PECAM-1:
-
Platelet endothelial cell adhesion molecule
- BAM:
-
Border-associated macrophage
- FPKM:
-
Fragments per kilobase of transcript per million
- ROI:
-
Region of interest
- DAM:
-
Disease associated microglia
- IRF:
-
Interferon regulatory factor
- bHLH/bZIP:
-
Basic helix-loop-helix/bZip
- RHD:
-
Rel homology domain
- Zf:
-
C2H2 zinc finger
- Stat:
-
Signal transducer and activator of transcription
- BBB:
-
Blood brain barrier
- TREM2:
-
Triggering receptor expressed on myeloid cells 2
- TBI:
-
Traumatic brain injury
References
Skaper SD, Facci L, Zusso M, Giusti P. An inflammation-centric view of neurological disease: beyond the neuron. Front Cell Neurosci. 2018;12:72.
Jyonouchi H. Innate Immunity and Neuroinflammation in Neuropsychiatric Conditions Including Autism Spectrum Disorders: Role of Innate Immune Memory. Cytokines. IntechOpen; 2020. doi: 10.5772/intechopen.87167
Zhan Y, Paolicelli RC, Sforazzini F, Weinhard L, Bolasco G, Pagani F, et al. Deficient neuron-microglia signaling results in impaired functional brain connectivity and social behavior. Nat Neurosci. 2014;17(3):400–6.
Neher JJ, Cunningham C. Priming microglia for innate immune memory in the brain. Trends Immunol. 2019;40(4):358–74.
Hammond TR, Robinton D, Stevens B. Microglia and the brain: complementary partners in development and disease. Annu Rev Cell Dev Biol. 2018;34(1):523–44.
Schwarz JM, Sholar PW, Bilbo SD. Sex differences in microglial colonization of the developing rat brain. J Neurochem. 2012;120(6):948–63.
Schafer DP, Lehrman EK, Kautzman AG, Koyama R, Mardinly AR, Yamasaki R, et al. Microglia sculpt postnatal neural circuits in an activity and complement-dependent manner. Neuron. 2012;74(4):691–705.
Wang C, Yue H, Hu Z, Shen Y, Ma J, Li J, et al. Microglia mediate forgetting via complement-dependent synaptic elimination. Science. 2020;367(6478):688–94.
Petrasch-Parwez E, Schöbel A, Benali A, Moinfar Z, Förster E, Brüne M, et al. Lateralization of increased density of Iba1-immunopositive microglial cells in the anterior midcingulate cortex of schizophrenia and bipolar disorder. Eur Arch Psychiatry Clin Neurosci. 2020;270(7):819–28.
Hopperton KE, Mohammad D, Trépanier MO, Giuliano V, Bazinet RP. Markers of microglia in post-mortem brain samples from patients with Alzheimer’s disease: a systematic review. Mol Psychiatry. 2018;23(2):177–98.
Suzuki K, Sugihara G, Ouchi Y, Nakamura K, Futatsubashi M, Takebayashi K, et al. Microglial activation in young adults with autism spectrum disorder. JAMA Psychiat. 2013;70(1):49–58.
Vargas DL, Nascimbene C, Krishnan C, Zimmerman AW, Pardo CA. Neuroglial activation and neuroinflammation in the brain of patients with autism. Ann Neurol. 2005;57(1):67–81.
Salvador AF, de Lima KA, Kipnis J. Neuromodulation by the immune system: a focus on cytokines. Nat Rev Immunol. 2021;21(8):526–41.
Netea MG, Domínguez-Andrés J, Barreiro LB, Chavakis T, Divangahi M, Fuchs E, et al. Defining trained immunity and its role in health and disease. Nat Rev Immunol. 2020;20(6):375–88.
Sherwood ER, Burelbach KR, McBride MA, Stothers CL, Owen AM, Hernandez A, et al. Innate immune memory and the host response to infection. J Immunol. 2022;208(4):785–92.
Graham DB, Xavier RJ. Conditioning of the immune system by the microbiome. Trends Immunol. 2023;44(7):499–511.
Ajami B, Bennett JL, Krieger C, Tetzlaff W, Rossi FMV. Local self-renewal can sustain CNS microglia maintenance and function throughout adult life. Nat Neurosci. 2007;10(12):1538–43.
Bruttger J, Karram K, Wörtge S, Regen T, Marini F, Hoppmann N, et al. Genetic cell ablation reveals clusters of local self-renewing microglia in the mammalian central nervous system. Immunity. 2015;43(1):92–106.
Wendeln AC, Degenhardt K, Kaurani L, Gertig M, Ulas T, Jain G, et al. Innate immune memory in the brain shapes neurological disease hallmarks. Nature. 2018;556(7701):332–8.
Turner RC, Naser ZJ, Lucke-Wold BP, Logsdon AF, Vangilder RL, Matsumoto RR, et al. Single low-dose lipopolysaccharide preconditioning: neuroprotective against axonal injury and modulates glial cells. Neuroimmunol Neuroinflamm. 2017;4:6–15.
Savage JC, Carrier M, Tremblay MÈ. Morphology of microglia across contexts of health and disease. Methods Mol Biol. 2019;2034:13–26.
Tan YL, Yuan Y, Tian L. Microglial regional heterogeneity and its role in the brain. Mol Psychiatry. 2020;25(2):351–67.
Kay M, Elkin LA, Higgins JJ, Wobbrock JO. ARTool: Aligned Rank Transform for Nonparametric Factorial ANOVAs. R package version 0.11.1; 2021. https://github.com/mjskay/ARTool. https://doi.org/10.5281/zenodo.594511.
Vogel-Ciernia A, Wood MA. Examining object location and object recognition memory in mice. Curr Protoc Neurosci. 2014;69:8.31.1-8.31.17.
Chiu K, Lau WM, Lau HT, So KF, Chang RCC. Micro-dissection of rat brain for RNA or protein extraction from specific brain region. J Vis Exp. 2007;7:269.
Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47.
Kolde R. pheatmap: pretty heatmaps; 2019. https://CRAN.R-project.org/package=pheatmap
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation. 2021;2(3): 100141.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Marc Carlson. org.Mm.eg.db: genome wide annotation for Mouse.
Jao J, Ciernia AV. MGEnrichment: a web application for microglia gene list enrichment analysis 2 3. PLoS Comput Biol. 2021;17(11): e1009160.
Duttke SH, Chang MW, Heinz S, Benner C. Identification and dynamic quantification of regulatory elements using total RNA. Genome Res; 2019. https://genome.cshlp.org/content/early/2019/10/24/gr.253492.119. Accessed 12 Mar 2024.
Saeki K, Pan R, Lee E, Kurotaki D, Ozato K. IRF8 configures enhancer landscape in postnatal microglia and directs microglia specific transcriptional programs. bioRxiv. 2023. https://doi.org/10.1101/2023.06.25.546453v2.
Langlais D, Barreiro LB, Gros P. The macrophage IRF8/IRF1 regulome is required for protection against infections and is associated with chronic inflammation. J Exp Med. 2016;213(4):585–603.
Barish GD, Yu RT, Karunasiri M, Ocampo CB, Dixon J, Benner C, et al. Bcl-6 and NF-kappaB cistromes mediate opposing regulation of the innate immune response. Genes Dev. 2010;24(24):2760–5.
Link VM, Duttke SH, Chun HB, Holtman IR, Westin E, Hoeksema MA, et al. Analysis of genetically diverse macrophages reveals local and domain-wide mechanisms that control transcription factor binding and function. Cell. 2018;173(7):1796-1809.e17.
Gosselin D, Skola D, Coufal NG, Holtman IR, Schlachetzki JCM, Sajti E, et al. An environment-dependent transcriptional network specifies human microglia identity. Science. 2017;356(6344):eaal3222.
Neph S, Kuehn MS, Reynolds AP, Haugen E, Thurman RE, Johnson AK, et al. BEDOPS: high-performance genomic feature operations. Bioinformatics. 2012;28(14):1919–20.
Zhang Y, Chen K, Sloan SA, Bennett ML, Scholze AR, O’Keeffe S, et al. An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci. 2014;34(36):11929–47.
Zeisel A, Hochgerner H, Lönnerberg P, Johnsson A, Memic F, van der Zwan J, et al. Molecular architecture of the mouse nervous system. Cell. 2018;174(4):999-1014.e22.
Tamayo JM, Rose D, Church JS, Schwartzer JJ, Ashwood P. Maternal allergic asthma induces prenatal neuroinflammation. Brain Sci. 2022;12(8):1041.
Ashwood P, Krakowiak P, Hertz-Picciotto I, Hansen R, Pessah I, Van de Water J. Elevated plasma cytokines in autism spectrum disorders provide evidence of immune dysfunction and are associated with impaired behavioral outcome. Brain Behav Immun. 2011;25(1):40–5.
Osman HC, Moreno R, Rose D, Rowland ME, Ciernia AV, Ashwood P. Impact of maternal immune activation and sex on placental and fetal brain cytokine and gene expression profiles in a preclinical model of neurodevelopmental disorders. J Neuroinflamm. 2024;21(1):118.
Tamayo JM, Osman HC, Schwartzer JJ, Pinkerton KE, Ashwood P. Characterizing the neuroimmune environment of offspring in a novel model of maternal allergic asthma and particulate matter exposure. J Neuroinflamm. 2023;20(1):252.
Kim J, Pavlidis P, Ciernia AV. Development of a High-Throughput Pipeline to Characterize Microglia Morphological States at a Single-Cell Resolution. eNeuro. 2024;11(7):ENEURO.0014-24.2024. https://doi.org/10.1523/ENEURO.0014-24.2024.
Pinskiy V, Jones J, Tolpygo AS, Franciotti N, Weber K, Mitra P. High-throughput method of whole-brain sectioning, using the tape-transfer technique. PLoS ONE. 2015;10(7): e0102363.
Terstege DJ, Oboh DO, Epp JR. FASTMAP: open-source flexible atlas segmentation tool for multi-area processing of biological images. eNeuro. 2022. https://doi.org/10.1523/ENEURO.0325-21.2022.
Brooks ME, Kristensen K, van Benthem KJ, Magnusson A, Berg CW, Nielsen A, et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 2017;9(2):378–400.
Schmid B, Schindelin J, Cardona A, Longair M, Heisenberg M. A high-level 3D visualization API for Java and ImageJ. BMC Bioinform. 2010;11(1):274.
Meleady L, Towriss M, Kim J, Bacarac V, Dang V, Rowland ME, et al. Histone deacetylase 3 regulates microglial function through histone deacetylation. Epigenetics. 2023;18(1):2241008.
Towriss M, MacVicar B, Ciernia AV. Modelling microglial innate immune memory in vitro: understanding the role of aerobic glycolysis in innate immune memory. Int J Mol Sci. 2023;24(10):8967.
Vogel Ciernia A, Link VM, Careaga M, LaSalle J, Ashwood P. Genetic variants drive altered epigenetic regulation of endotoxin response in BTBR macrophages. Brain Behav Immun. 2020;89:20–31.
Zhang X, Kracht L, Lerario AM, Dubbelaar ML, Brouwer N, Wesseling EM, et al. Epigenetic regulation of innate immune memory in microglia. J Neuroinflamm. 2022;19(1):111.
Naler LB, Hsieh YP, Geng S, Zhou Z, Li L, Lu C. Epigenomic and transcriptomic analyses reveal differences between low-grade inflammation and severe exhaustion in LPS-challenged murine monocytes. Commun Biol. 2022;5(1):1–17.
Kusiak A, Brady G. Bifurcation of signalling in human innate immune pathways to NF-kB and IRF family activation. Biochem Pharmacol. 2022;205: 115246.
Wang AG, Son M, Kenna E, Thom N, Tay S. NF-κB memory coordinates transcriptional responses to dynamic inflammatory stimuli. Cell Rep. 2022;40(7): 111159.
Iwanaszko M, Kimmel M. NF-κB and IRF pathways: cross-regulation on target genes promoter level. BMC Genomics. 2015;16(1):307.
Platanitis E, Decker T. Regulatory Networks Involving STATs, IRFs, and NFκB in Inflammation. Front Immunol. 2018;9:2542.
Kaminska B, Mota M, Pizzi M. Signal transduction and epigenetic mechanisms in the control of microglia activation during neuroinflammation. Biochim et Biophys Acta BBA Mol Basis Dis. 2016;1862(3):339–51.
Maurya SK, Gupta S, Mishra R. Transcriptional and epigenetic regulation of microglia in maintenance of brain homeostasis and neurodegeneration. Front Mol Neurosci. 2023. https://doi.org/10.3389/fnmol.2022.1072046.
Mass E, Nimmerjahn F, Kierdorf K, Schlitzer A. Tissue-specific macrophages: how they develop and choreograph tissue biology. Nat Rev Immunol. 2023;23(9):563–79.
Bennett FC, Bennett ML, Yaqoob F, Mulinyawe SB, Grant GA, Gephart MH, et al. A combination of ontogeny and CNS environment establishes microglial identity. Neuron. 2018;98(6):1170-1183.e8.
Novakovic B, Habibi E, Wang SY, Arts RJW, Davar R, Megchelenbrink W, et al. β-glucan reverses the epigenetic state of LPS-induced immunological tolerance. Cell. 2016;167(5):1354-1368.e14.
Foster SL, Hargreaves DC, Medzhitov R. Gene-specific control of inflammation by TLR-induced chromatin modifications. Nature. 2007;447(7147):972–8.
Chen Z, Jalabi W, Shpargel KB, Farabaugh KT, Dutta R, Yin X, et al. Lipopolysaccharide-induced microglial activation and neuroprotection against experimental brain injury is independent of hematogenous TLR4. J Neurosci. 2012;32(34):11706–15.
Balendran T, Lim K, Hamilton JA, Achuthan AA. Targeting transcription factors for therapeutic benefit in rheumatoid arthritis. Front Immunol. 2023;14:1196931.
Masuda T, Tsuda M, Inoue K. Transcriptional regulation in microglia and neuropathic pain. Pain Management. 2016;6(2):91–4.
Gao T, Jernigan J, Raza SA, Dammer EB, Xiao H, Seyfried NT, et al. Transcriptional regulation of homeostatic and disease-associated-microglial genes by IRF1, LXRβ, and CEBPα. Glia. 2019;67(10):1958–75.
Masuda T, Tsuda M, Yoshinaga R, Tozaki-Saitoh H, Ozato K, Tamura T, et al. IRF8 is a critical transcription factor for transforming microglia into a reactive phenotype. Cell Rep. 2012;1(4):334–40.
Masuda T, Nishimoto N, Tomiyama D, Matsuda T, Tozaki-Saitoh H, Tamura T, et al. IRF8 is a transcriptional determinant for microglial motility. Purinergic Signal. 2014;10(3):515–21.
Horiuchi M, Wakayama K, Itoh A, Kawai K, Pleasure D, Ozato K, et al. Interferon regulatory factor 8/interferon consensus sequence binding protein is a critical transcription factor for the physiological phenotype of microglia. J Neuroinflamm. 2012;9:227.
Orihuela R, McPherson CA, Harry GJ. Microglial M1/M2 polarization and metabolic states. Br J Pharmacol. 2016;173(4):649–65.
Banks WA, Gray AM, Erickson MA, Salameh TS, Damodarasamy M, Sheibani N, et al. Lipopolysaccharide-induced blood-brain barrier disruption: roles of cyclooxygenase, oxidative stress, neuroinflammation, and elements of the neurovascular unit. J Neuroinflamm. 2015;12(1):223.
Vargas-Caraveo A, Sayd A, Maus SR, Caso JR, Madrigal JLM, García-Bueno B, et al. Lipopolysaccharide enters the rat brain by a lipoprotein-mediated transport mechanism in physiological conditions. Sci Rep. 2017;7(1):13113.
Vichaya EG, Malik S, Sominsky L, Ford BG, Spencer SJ, Dantzer R. Microglia depletion fails to abrogate inflammation-induced sickness in mice and rats. J Neuroinflamm. 2020;17:172.
Chistyakov DV, Astakhova AA, Azbukina NV, Goriainov SV, Chistyakov VV, Sergeeva MG. Cellular model of endotoxin tolerance in astrocytes: role of interleukin 10 and oxylipins. Cells. 2019;8(12):1553.
Schirmbeck GH, Seady M, Fróes FT, Taday J, Da Ré C, Souza JM, et al. Long-term LPS systemic administration leads to memory impairment and disturbance in astrocytic homeostasis. Neurotoxicology. 2023;99:322–31.
Bian Y, Zhao X, Li M, Zeng S, Zhao B. Various roles of astrocytes during recovery from repeated exposure to different doses of lipopolysaccharide. Behav Brain Res. 2013;253:253–61.
Beurel E. HDAC6 regulates LPS-tolerance in astrocytes. PLoS ONE. 2011;6(10): e25804.
Zamanian JL, Xu L, Foo LC, Nouri N, Zhou L, Giffard RG, et al. Genomic analysis of reactive astrogliosis. J Neurosci. 2012;32(18):6391–410.
Hasel P, Rose IVL, Sadick JS, Kim RD, Liddelow SA. Neuroinflammatory astrocyte subtypes in the mouse brain. Nat Neurosci. 2021;24(10):1475–87.
Liddelow SA, Guttenplan KA, Clarke LE, Bennett FC, Bohlen CJ, Schirmer L, et al. Neurotoxic reactive astrocytes are induced by activated microglia. Nature. 2017;541(7638):481–7.
Jung H, Lee D, You H, Lee M, Kim H, Cheong E, et al. LPS induces microglial activation and GABAergic synaptic deficits in the hippocampus accompanied by prolonged cognitive impairment. Sci Rep. 2023;13(1):6547.
Marsh B, Stevens SL, Packard AEB, Gopalan B, Hunter B, Leung PY, et al. Systemic lipopolysaccharide protects the brain from ischemic injury by reprogramming the response of the brain to stroke: a critical role for IRF3. J Neurosci. 2009;29(31):9839–49.
Rosenzweig HL, Lessov NS, Henshall DC, Minami M, Simon RP, Stenzel-Poore MP. Endotoxin preconditioning prevents cellular inflammatory response during ischemic neuroprotection in mice. Stroke. 2004;35(11):2576–81.
Hickey EJ, You X, Kaimaktchiev V, Stenzel-Poore M, Ungerleider RM. Lipopolysaccharide preconditioning induces robust protection against brain injury resulting from deep hypothermic circulatory arrest. J Thorac Cardiovasc Surg. 2007;133(6):1588–96.
Paolicelli RC, Sierra A, Stevens B, Tremblay ME, Aguzzi A, Ajami B, et al. Microglia states and nomenclature: a field at its crossroads. Neuron. 2022;110(21):3458–83.
Meng J, Han L, Xu H, Zhang L, Liu Z, Zhou Y, et al. TREM2 regulates microglial phagocytosis of synapses in innate immune tolerance. Int Immunopharmacol. 2024;127: 111445.
Hanamsagar R, Alter MD, Block CS, Sullivan H, Bolton JL, Bilbo SD. Generation of a microglial developmental index in mice and in humans reveals a sex difference in maturation and immune reactivity. Glia. 2017;65(9):1504–20.
Paolicelli RC, Bolasco G, Pagani F, Maggi L, Scianni M, Panzanelli P, et al. Synaptic pruning by microglia is necessary for normal brain development. Science. 2011;333(6048):1456–8.
Cengiz P, Zafer D, Chandrashekhar JH, Chanana V, Bogost J, Waldman A, et al. Developmental differences in microglia morphology and gene expression during normal brain development and in response to hypoxia-ischemia. Neurochem Int. 2019;127:137–47.
Parakalan R, Jiang B, Nimmi B, Janani M, Jayapal M, Lu J, et al. Transcriptome analysis of amoeboid and ramified microglia isolated from the corpus callosum of rat brain. BMC Neurosci. 2012;13(1):64.
Ayoub AE, Salm AK. Increased morphological diversity of microglia in the activated hypothalamic supraoptic nucleus. J Neurosci. 2003;23(21):7759–66.
Bouvier DS, Jones EV, Quesseveur G, Davoli MA, Ferreira AT, Quirion R, et al. High resolution dissection of reactive glial nets in Alzheimer’s disease. Sci Rep. 2016;6:24544.
Bachstetter AD, Van Eldik LJ, Schmitt FA, Neltner JH, Ighodaro ET, Webster SJ, et al. Disease-related microglia heterogeneity in the hippocampus of Alzheimer’s disease, dementia with Lewy bodies, and hippocampal sclerosis of aging. Acta Neuropathol Commun. 2015;3:32.
Mondo E, Becker SC, Kautzman AG, Schifferer M, Baer CE, Chen J, et al. A Developmental analysis of juxtavascular microglia dynamics and interactions with the vasculature. J Neurosci. 2020;40(34):6503.
Bisht K, Okojie KA, Sharma K, Lentferink DH, Sun YY, Chen HR, et al. Capillary-associated microglia regulate vascular structure and function through PANX1-P2RY12 coupling in mice. Nat Commun. 2021;12(1):5289.
Morris GP, Foster CG, Courtney JM, Collins JM, Cashion JM, Brown LS, et al. Microglia directly associate with pericytes in the central nervous system. Glia. 2023;71(8):1847–69.
Jolivel V, Bicker F, Binamé F, Ploen R, Keller S, Gollan R, et al. Perivascular microglia promote blood vessel disintegration in the ischemic penumbra. Acta Neuropathol. 2015;129(2):279–95.
Davalos D, Kyu Ryu J, Merlini M, Baeten KM, Le Moan N, Petersen MA, et al. Fibrinogen-induced perivascular microglial clustering is required for the development of axonal damage in neuroinflammation. Nat Commun. 2012;3(1):1227.
Haruwaka K, Ikegami A, Tachibana Y, Ohno N, Konishi H, Hashimoto A, et al. Dual microglia effects on blood brain barrier permeability induced by systemic inflammation. Nat Commun. 2019;10(1):5816.
Ziebell JM, Taylor SE, Cao T, Harrison JL, Lifshitz J. Rod microglia: elongation, alignment, and coupling to form trains across the somatosensory cortex after experimental diffuse brain injury. J Neuroinflamm. 2012;9(1):247.
Taylor SE, Morganti-Kossmann C, Lifshitz J, Ziebell JM. Rod microglia: a morphological definition. PLoS ONE. 2014;9(5): e97096.
Giordano KR, Denman CR, Dubisch PS, Akhter M, Lifshitz J. An update on the rod microglia variant in experimental and clinical brain injury and disease. Brain Commun. 2021;3(1):fcaa227.
Réu P, Khosravi A, Bernard S, Mold JE, Salehpour M, Alkass K, et al. The lifespan and turnover of microglia in the human brain. Cell Rep. 2017;20(4):779–84.
Lavin Y, Winter D, Blecher-Gonen R, David E, Keren-Shaul H, Merad M, et al. Tissue-resident macrophage enhancer landscapes are shaped by the local microenvironment. Cell. 2014;159(6):1312–26.
Gosselin D, Link VM, Romanoski CE, Fonseca GJ, Eichenfield DZ, Spann NJ, et al. Environment drives selection and function of enhancers controlling tissue-specific macrophage identities. Cell. 2014;159(6):1327–40.
Acknowledgements
We thank the Ciernia Lab, Ashwood Lab, and Pavlidis Lab members for their thoughtful feedback and suggestions during lab meetings throughout the progression of this project. We are grateful for the resources provided by the Neuroimaging & Neurocomputation Centre and the Dynamic Brain Circuits in Health & Disease Research Cluster at the University of British Columbia, as well as the University of California, Davis DNA Technologies & Expression Analysis Core.
Funding
This work was supported by the Canadian Open Neuroscience Platform Student Scholar Award (10901 to JK); University of British Columbia Four Year Doctoral Fellowship (6569 to JK); Canadian Institutes of Health Research CGS-M to OS; Canadian Institutes of Health Research (CRC-RS 950-232402, AWD-025853, AWD-025854 to AC); Natural Sciences and Engineering Research Council of Canada (RGPIN-2019-04450, DGECR-2019-00069 to AC); Scottish Rite Charitable Foundation (21103 to AC); Brain Canada Foundation (AWD-023132 to AC); Michael Smith Health Research Foundation (AWD-005509 to AC); National Institute of Environmental Health Sciences (R21ES035492, R21ES035969 to PA); National Institutes of Child Health (R01HD090214 to PA); National Institute of Mental Health (R21MH116383, R01MH118209 to PA; K01MH116389 to AC), and the Brain Foundation to PA. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Authors and Affiliations
Contributions
P.A. and A.V.C designed the initial experiments and A.V.C. and J.K. conceptualized the project. J.K. performed the majority of the experiments, analysis, and manuscript writing. A.V.C., J.K., and J.J. performed bioinformatic analysis. J.T. and A.M.M. assisted with injections and tissue collection for the RNAseq experiment and performed the Luminex assay. O.S. and K.L. performed the behavioral battery and analysis of the behavioral battery. B.W. assisted with RT-qPCR experiments. All authors read and edited the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
12974_2024_3198_MOESM2_ESM.xlsx
Supplementary Material 2. Table S1: Kruskal-Wallis rank sum test results, Dunn-corrected posthocs, and input data for Luminex heatmap Figure 1C.
12974_2024_3198_MOESM3_ESM.xlsx
Supplementary Material 3. Table S2: ANOVA results and BH-corrected posthocs from sickness behavior and weight change analysis, and 1-way ANOVA results for each behavioral task in behavioral battery. Related to Figures 1B and 2.
12974_2024_3198_MOESM4_ESM.xlsx
Supplementary Material 4. Table S3: RNAseq data : Sample information, differential expression analysis, genes in DEG clusters, and gene ontology terms for genes in DEG clusters; TF motif analysis data : locations of Homer known motifs in promoters of 2xLPS-sensitive cluster DEGs and : LOLA ChIPseq peak enrichment results; MGEnrichment data: Enrichments against microglial gene lists results.
12974_2024_3198_MOESM5_ESM.xlsx
Supplementary Material 5. Table S4: ANOVA results and BH-corrected posthocs from MicrogliaMorphology cluster analysis, microglia alignment event occurrence analysis, and microglial density analysis related to Figures 5, 6, and S5F.
12974_2024_3198_MOESM6_ESM.xlsx
Supplementary Material 6. Table S5: ANOVA results and BH-corrected posthocs from RT-qPCR data analysis, related to Figures 1D; S1; and S6A,C,D,F.
12974_2024_3198_MOESM7_ESM.zip
Supplementary Material 7. File S1: Homer software output for transcription factor motif analysis of 2xLPS-sensitive cluster gene promoters. Related to Figures 4A and C.
12974_2024_3198_MOESM8_ESM.zip
Supplementary Material 8. File S2: Homer software output for transcription factor motif analysis of 4xLPS-sensitive cluster gene promoters.
12974_2024_3198_MOESM9_ESM.zip
Supplementary Material 9. File S3: Homer software output for transcription factor motif analysis of LPS-decreased cluster gene promoters.
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
Kim, J., Sullivan, O., Lee, K. et al. Repeated LPS induces training and tolerance of microglial responses across brain regions. J Neuroinflammation 21, 233 (2024). https://doi.org/10.1186/s12974-024-03198-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12974-024-03198-1