- Research
- Open access
- Published:
Interaction of high-fat diet and brain trauma alters adipose tissue macrophages and brain microglia associated with exacerbated cognitive dysfunction
Journal of Neuroinflammation volume 21, Article number: 113 (2024)
Abstract
Obesity increases the morbidity and mortality of traumatic brain injury (TBI). Detailed analyses of transcriptomic changes in the brain and adipose tissue were performed to elucidate the interactive effects between high-fat diet-induced obesity (DIO) and TBI. Adult male mice were fed a high-fat diet (HFD) for 12 weeks prior to experimental TBI and continuing after injury. High-throughput transcriptomic analysis using Nanostring panels of the total visceral adipose tissue (VAT) and cellular components in the brain, followed by unsupervised clustering, principal component analysis, and IPA pathway analysis were used to determine shifts in gene expression patterns and molecular pathway activity. Cellular populations in the cortex and hippocampus, as well as in VAT, during the chronic phase after combined TBI-HFD showed amplification of central and peripheral microglia/macrophage responses, including superadditive changes in selected gene expression signatures and pathways. Furthermore, combined TBI and HFD caused additive dysfunction in Y-Maze, Novel Object Recognition (NOR), and Morris water maze (MWM) cognitive function tests. These novel data suggest that HFD-induced obesity and TBI can independently prime and support the development of altered states in brain microglia and VAT, including the disease-associated microglia/macrophage (DAM) phenotype observed in neurodegenerative disorders. The interaction between HFD and TBI promotes a shift toward chronic reactive microglia/macrophage transcriptomic signatures and associated pro-inflammatory disease-altered states that may, in part, underlie the exacerbation of cognitive deficits. Thus, targeting of HFD-induced reactive cellular phenotypes, including in peripheral adipose tissue immune cell populations, may serve to reduce microglial maladaptive states after TBI, attenuating post-traumatic neurodegeneration and neurological dysfunction.
Background
Traumatic brain injury (TBI) is a debilitating disorder that affects more than 1.5 million individuals annually in the United States [1]. Obesity is among the most prevalent pre-existing conditions that can negatively impact TBI outcomes. Currently, approximately 13% of the world population is considered clinically obese [body mass index BMI greater than 30 kg/m2] [2] placing a large burden on healthcare systems [3]. Emerging evidence indicates that obese patients have increased TBI complications and higher mortality rates (36%) [4,5,6,7].
Activation of brain and peripheral inflammation pathways is a key pathophysiological feature in both TBI [8,9,10,11,12] and obesity [12,13,14,15,16,17,18,19,20] and has been implicated in the development of associated neurological dysfunction [8, 21, 22]. Notably, evidence suggests that pre-existing diet-induced obesity is associated with an amplification of post-TBI pro-inflammatory responses and increased microglia-altered states in several brain regions [23,24,25]. Additional studies report diet-induced exacerbations in TBI-induced cognitive decline [24, 25] and central (brain) insulin resistance [26].
Microglia play a key role in mediating the inflammatory responses after TBI. In both clinical and experimental TBI, microglia undergo chronic reactive changes that may contribute to long-term neurodegenerative processes and cognitive decline [10, 11]. In animal models, pharmacological depletion of microglia or other interventions that reduce chronic microglial activation decrease TBI-induced neuroinflammation and associated neurological deficits [8, 27, 28]. Microglia also mediate hypothalamic inflammation and neuronal stress responses in diet-induced obesity [19]. Selective pharmacological inhibition of microglial phagocytosis limits diet-induced obesity hyperphagia and weight gain, decreases dendritic spine loss, and reduces cognitive dysfunction [21].
Obesity triggers peripheral inflammation that includes reactive changes in adipose tissue macrophages (ATM), which account for > 50% of the cells in the increased visceral adipose tissue (VAT) [29]. VAT is an immunogenic tissue and the resulting state of obesity-induced chronic low-grade inflammation and increased secretion of pro-inflammatory factors into the circulation contribute to brain neuroinflammatory responses [29, 30]. Obesity is also associated with increased cognitive decline and dementia [31,32,33]. Inflammation may play a role in driving the negative pathological consequences of diet-induced obesity, including behavioral deficits [34,35,36,37]. NLR family, pyrin domain-containing 3 (NLRP3) contributes to obesity-induced inflammation [36, 38] through NLRP3-induced activation of microglial IL-1 receptor 1 (IL-1R1) [36]. Targeted knockout of NLRP3 at the level of the VAT attenuates diet-induced obesity cognitive deficits, thus implicating a potential role for NLRP3/IL-1β signaling in brain-visceral adipose tissue interactions related to obesity [36]. Furthermore, the depletion of myeloid NOX-2, decreases VAT inflammation, VAT macrophage infiltration, and cognitive dysfunction [37].
HFD-induced obesity may promote the development of VAT macrophage reactive states that prime the development of deleterious brain microglia phenotypes, increasing vulnerability to secondary insults including TBI. Thus, the amplification of chronic microglial maladaptive states in combined TBI-HFD may contribute to the observed exacerbation of cognitive dysfunction [13, 23]. Although the individual effects of TBI on the brain and of diet on the adipose tissue and brain have been well studied, the present study examines several important issues that have not been addressed. These include: the effect of isolated TBI on the adipose tissue inflammatory environment; and how interactive effects of diet-induced obesity and TBI affect transcriptomic changes associated with microglia/macrophage states centrally and systemically to impact neurological function.
Materials and methods
Animals
Studies were performed using adult male C57Bl/6 J mice (8–10 weeks old; Jackson Laboratories). Mice were housed in shoebox cages (5 mice in each cage) at least 1 week prior to any procedures in a room (22–23 °C) with a 12-h/12-h light–dark cycle. Food and water were provided ad libitum. The mice were handled briefly before use. Procedures were conducted from 10:00 to 17:00 in a quiet room. All protocols involving the use of animals complied with the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health (NIH) (DHEW publication NIH 85-23-2985) and were approved by the University of Maryland Animal Use Committee.
Experimental design
To induce our model of DIO, C57Bl/6 mice (n = 8–12/group) were placed on either a high-fat diet (HFD) (60 kcal% Fat; D12492, Research Diets) or a standard diet (SD) (10 kcal% Fat; D12450B, Research Diets) for a period of 12 weeks prior to exposure to controlled cortical impact (CCI) or Sham surgery. Group sizes were determined based on 25 ± 10% (TBI) and 50 ± 10% (diet), anticipated differences between groups on functional recovery and metabolic parameters, respectively, with a α = 0.05, statistical power of 80%. Standard deviation used in sample size calculation was obtained from previous/pilot experiments within our laboratory.
Cohort 1: Following exposure to CCI/Sham surgery, mice were anesthetized (100 mg/kg sodium pentobarbital, I.P.) and transcardially perfused with ice-cold 0.9% saline (100 ml) at 28 days post-injury (dpi). VAT, brain tissue, and blood samples were collected and stored at − 80 °C until processed for RNA/protein analysis.
Cohort 2: To assess more chronic effects of diet and/or TBI, on neurobehavioral, neuroinflammatory and neurodegenerative outcomes, a second cohort of mice was fed a SD or HFD diet for 12 weeks prior to induction of CCI/Sham surgery. Mice underwent a battery of neurobehavioral tasks to assess cognitive behaviors throughout 90 dpi. Mice were anesthetized (100 mg/kg sodium pentobarbital, I.P.) and transcardially perfused with ice-cold 0.9% saline (100 ml) at 90 dpi. Cluster of differentiation 11b (CD11b) positively selected cells were isolated from hippocampal and cortical tissue using Miltenyi MACS Cell separation. RNA was extracted from both the CD11b positively selected cells (primarily microglia but also resident/infiltrating macrophages) as well as flowthrough (the remaining brain cells, including neurons, astrocytes, and oligodendrocytes) and analyzed using a Nanostring Glial Panel. RNA was also extracted from the VAT and analyzed using a Nanostring Neuroinflammation panel.
Controlled cortical impact
Our custom-designed controlled cortical impact (CCI) injury device consists of a microprocessor-controlled pneumatic impactor with a 3.5 mm diameter tip. Mice were anesthetized with isoflurane evaporated in a gas mixture containing 70% N2O and 30% O2 and administered through a nose mask (induction at 4% and maintenance at 2%). The depth of anesthesia was assessed by monitoring respiration rate and pedal withdrawal reflexes. Mice were placed on a heated pad, and core body temperature was maintained at 37 °C. The head was mounted on a stereotaxic frame, and the surgical site was clipped and cleaned with Nolvasan and ethanol scrubs. A 10-mm midline incision was made over the skull, the skin and fascia were reflected, and a 5-mm craniotomy was made on the central aspect of the left parietal bone. The impounder tip of the injury device was then extended to its full stroke distance (44 mm), positioned to the surface of the exposed dura, and reset to impact the cortical surface. Moderate-level CCI (n = 12) was induced using an impactor velocity of 6 m/s and deformation depth of 2 mm as previously described [8, 39,40,41,42]. After injury, the incision was closed with interrupted 6–0 silk sutures, anesthesia was terminated, and the animal was placed into a heated cage to maintain normal core temperature for 45 min post-injury. Sham animals underwent the same procedure as TBI mice except for the craniotomy and impact. All animals were monitored daily post-injury.
Neurobehavioral testing
Y-maze spontaneous alternation
The Y-maze was carried out at 70 days dpi to access spatial working memory and was essentially performed as previously described [8, 41]. Briefly, the Y-maze (Stoelting Co., Wood Dale, IL) consisted of three identical arms, each arm 35 cm long, 5 cm wide, and 10 cm high, at an angle of 120° with respect to the other arms. One arm was randomly selected as the “start” arm, and the mouse was placed within and allowed to explore the maze freely for 5 min. Arm entries (arms A–C) were recorded by analyzing mouse activity using ANY-maze software (Stoelting Co., Wood Dale, IL). An arm entry was attributed when all four paws of the mouse entered the arm, and an alternation was designated when the mouse entered three different arms consecutively. The percentage of alternation was calculated as follows: total alternations × 100/ (total arm entries − 2). If a mouse scored > 50% alternations (the chance level for choosing the unfamiliar arm), this was indicative of proper spatial working memory.
Novel object recognition task
The Novel object recognition (NOR) was carried out to assess non-spatial hippocampal-mediated memory as previously described with slight modifications, on 77–78 dpi [8]. Mice were placed in the NOR chamber where two identical objects were placed near the left and right corners of the open field for training (familiar phase) and allowed to freely explore until they spent a total of 20 s exploring the identical objects (exploration was recorded when the front paws or nose contacted the object). After 24 h, object recognition was tested by substituting a novel object for a familiar training object (the novel object location was counterbalanced across mice). The time spent with each object was recorded using the computed Any-Maze automated software. Because mice inherently prefer to explore novel objects, a preference for the novel object (more time than chance [15 s] spent with the novel object) indicates intact memory for the familiar object. Mice that failed to explore and remained immobile throughout the NOR test were excluded from the analysis.
Morris water maze (MWM)
Spatial learning and memory were assessed using the Morris Water Maze (MWM) task carried out at 80–84 dpi, as previously described [8]. The MWM protocol included two phases: (1) standard hidden platform training (acquisition) and (2) the twenty-four-hour probe test. Briefly, a circular tank (100 cm in diameter) was filled with water (23 ± 2 °C) and was surrounded by various extra-maze cues on the wall of the testing area. A transparent platform (10 cm in diameter) was submerged 0.5 cm below the surface of the water. Starting at 80 dpi, the mice were trained to find the hidden submerged platform located in the northeast (NE) quadrant of the tank for 4 consecutive days (80-83dpi). The mice underwent four trials per day, starting from a randomly selected release point (east, south, west, and north). Each mouse was allowed a maximum of 90 s to find the hidden submerged platform. The latency to the platform was recorded by the Any-Maze automated video tracking system. Reference memory was assessed by a probe test carried out at 24 h following the last acquisition day, on 84 dpi. The platform was removed, and the mice were released from the southwest (SW) position, and the time in the target quadrant was recorded. Additionally, search strategy analysis was performed as previously described [8]. Mice that failed to swim and remained immobile throughout the 90-s trial in the MWM test were excluded from the analysis.
Isolation of CD11b-positive cells
A magnetic bead-conjugated anti-cluster of differentiation 11b (CD11b) was used to isolate microglia/macrophages from ipsilateral (injured hemisphere) cortical and hippocampal tissue using Miltenyi MACS Separation Technology (Miltenyi Biotec, Auburn, CA) as per manufacturer’s instructions. Briefly, ipsilateral perilesional cortex and hippocampus from Sham and CCI mice were rapidly microdissected and a single cell suspension was prepared from the combined tissues (pooled tissue from two mice) using enzymatic digestion (Neural Tissue Dissociation Kit; Miltenyi Biotec) in combination with a gentleMACS Dissociator. Myelin was removed using Debris Remval Solution step (Miltenyi Biotec) and cells were incubated with anti-CD11b MicroBeads (Miltenyi Biotec) and loaded onto an LS column (Miltenyi Biotec) placed in the magnetic field of a MACS separator. The negative fraction (flow-through) was collected, and the column was washed 3 times with MACS buffer (Miltenyi Biotech). CD11b-positive cells were eluted by removing the magnetic field, resulting in the isolation of viable CD11b-positve cells (microglia/macrophages) from Sham and CCI mice. Cells were snap-frozen on liquid N2 for RNA extraction performed using Direct-zol RNA MicroPrep kit (Zymo Research).
Real-time PCR
Quantitative gene expression analysis in the adipose tissue as well as ipsilateral hippocampus and peri-lesional cortex of Sham and CCI mice was performed using Taqman technology as previously described [8, 41]. Real-time PCR for target mRNAs was performed using TaqMan gene expression assays (NADPH oxidase 2 (NOX2), Mm01287743_m1; human neutrophil cytochrome blight chain (p22phox), Mm00514478_m1; Interleukin-1beta (IL-1β), Mm01336189_m1; Tumor necrosis factor-alpha (TNF-α), Mm00443258_m1; Interleukin-10 (IL-10), Mm01288386_m1; NLR family pyrin domain containing 3 (NLRP3), Mm00840904_m1; glial fibrillary acidic protein (GFAP), Mm01253033_m1; Integrin Subunit Alpha M (ITGAM; CD11b). Mm00434455_m1; and GAPDH, Mm99999915_g1; Applied Biosystems, Carlsbad, CA) on a QUANTSTUDIO 5 Real Time PCR machine (Applied Biosystems). Samples were assayed in duplicate in one run (40 cycles), which was composed of 3 stages, 50 °C for 2 min, 95 °C for 10 s for each cycle (denaturation) and finally the transcription step at 60 °C for 1 min. Gene expression was normalized by GAPDH and compared to the control sample to determine relative expression values by 2−ΔΔCt method.
Flow cytometry analysis
Immediately following euthanasia, a 30 mg piece of adipose tissue was dissected from each mouse, minced into small pieces, and incubated with 150U/ml Collagenase IV (Worthington Biochemical Corporation Lakewood, NJ) and 10 mg/ml DNase II (Sigma) for 1 h at 37 °C in a rotational shaker. The suspension was passed through a 70 µm filter to mechanically dissociate adipose tissue. For immune cell surface markers, leukocytes were washed with FACS buffer (5% fetal bovine serum in 1 × HBSS) with sodium azide (NaN3) and blocked with 1:50 mouse Fc Block (anti-CD16/32) prior to staining with primary antibody-conjugated fluorophores at 1:50 concentration, including CD45-eF450, and CD11b-APCeF780. All antibodies were commercially purchased from Biolegend. For live/dead discrimination, a fixable viability dye, Zombie Aqua (Biolegend), was diluted at 1:100 in Hank's balanced salt solution (HBSS; Gibco). Cells were briefly fixed in 2% paraformaldehyde (PFA). Data were acquired on a LSRII using FACSDiva 6.0 (BD Biosciences) and analyzed using FlowJo (Treestar Inc.).
To measure reactive oxygen species (ROS) levels, leukocytes were incubated with dihydrorhodamine (DHR)123 (5 mM; 1:500 in RPMI; Ex/Em: 500/536), a cell-permeable fluorogenic probe (Life Technologies/Invitrogen, Waltham, MA). DHR123 passively diffuses into cells and is oxidized by peroxide and peroxynitrite, causing a reaction that produces a green fluorescence that can be measured by flow cytometry. Cells were loaded for 20 min in a 37 °C water bath, washed twice with FACS buffer (without NaAz), and then stained for surface markers including viability dye, and subsequently fixed in PFA.
Leptin, insulin, and monocyte chemoattractant protein (MCP)-1 determination
At either 28 or 90 dpi, serum Monocyte Chemoattractant Protein (MCP)-1, Leptin and Insulin were measured using the Mouse CCL2/JE/MCP-1 DuoSet ELISA (Cat No. DY479), Mouse/Rat Leptin Quantikine ELISA Kit (Cat No. MOB00B, R&D Bio systems, Minneapolis, MN, USA) and the Insulin Mouse ELISA Kit (Cat No. EMINS; Thermo-Fisher, Waltham, MA, USA). The assays were performed according to the manufacturer’s instructions.
Determination of resting blood glucose levels
Prior to euthanasia at 28 dpi, mice were fasted for 12 h. Using the AlphaTRAK Blood Glucose Monitoring System (Zoetis Inc. Parsippany, NJ, USA) resting blood glucose levels were measured via tail snip samples according to the manufacturer’s instructions.
Nanostring gene expression data analysis
The raw Nanostring files were analyzed using Rosalind software to generate the normalized gene expression data (used for all later analyses) and log2 fold changes. The average gene expression across groups as well as the expression Z-scores were calculated in Excel; where indicated, we used Z-score addition to create a composite score (Z-score sum) [43]. Two-way ANOVA in R was used to study the gene expression variability among groups and identify interaction effects between injury and diet, as well as main effects—genes which are affected only by the injury or only by diet, independent of each other—as previously described [44]. The resulting p values and log2 fold changes were imported into IPA package to analyze the activation of various molecular pathways including Canonical, Upstream regulators (genes, RNAs, and proteins); abbreviated throughout as Upstream regulators – Genes, Upstream regulators (drugs and chemicals); abbreviated throughout as Upstream regulators – Drugs, and Diseases and Bio Functions; abbreviated throughout as Diseases.
Two-way ANOVA in GraphPad Prism was used to determine the significance of differences between cumulative gene expression Z-scores and the pathway activation Z-scores generated by IPA software. PCA analysis of the Rosalind differentially expressed genes was performed in R using the pcaExplorer software package [45]. The principal components were derived from the Z-score values of the differentially expressed subset of genes.
Statistical analysis
Blinding was achieved by ensuring that the individual who carried out behavioral and stereological analyses were blinded to injury or diet groups. Quantitative data were expressed as mean ± SEM or mean ± STDEV with individual data points as indicated. Normality testing was performed; as datasets met normality requirements (D'Agostino and Pearson omnibus normality test), parametric analyses were used. Statistical analysis was performed using a two-way ANOVA with Tukey post-hoc tests with post-hoc statistics shown in the representative figures. When comparisons were made between two conditions, an unpaired Student's t test was performed. Statistical analyses utilized Prism v8 for Windows (GraphPad Software) or in R. Significance level was set at p < 0.05.
The consistent application of two-way ANOVA model permits rigorous determination of the effects of each of the two factors (diet and injury), as indicated by the separate factor significance as well as the direction, and magnitude of their combined effects as evidenced by the interaction and post-hoc test significance [46]. A non-significant interaction denotes a combined effect that is simply additive whereas, in a significant interaction, the effects of one of the factors of change depend on the condition of the other factor, and the combined effect reflects either a super additive amplification (synergy) or the opposite. Only significant comparisons evident in post-hoc tests are indicated in respective figures.
We used two-way ANOVA significant post-hoc p values to select genes for IPA analysis rather that t-test p-value (or multiple comparisons-adjusted p values). Although this approach is more stringent and may exclude some biologically significant changes, it is also more appropriate considering the multiple groups evaluated. Furthermore, the application of adjusted p values is not appropriate here in which we focus on Nanostring panels that include genes selected based on their contribution to similar pathways (e.g., inflammation) and may therefore be modulated in a similar fashion. The analysis of IPA data employed repeated two-way ANOVA to account for changes in specific pathways across groups. To perform the analyses, we introduced the following conventions: a baseline Sham-SD vs Sham-SD comparison with a z-score “0” across all pathways and replaced all IPA-generated N/A classifications (not able to detect activation level) by a score of “0”.
The absence of additive effects can also result in a significant interaction factor and may reflect a “ceiling effect” in which the relatively severe TBI model used in our studies prevents detection of the additive impairments after combined TBI and HFD insults, as was observed for cognitive deficits and lesion volume (an indicator of post-traumatic neurodegeneration, data not shown).
Results
High fat diet (HFD) increases body weight and elevates markers of metabolic dysfunction and systemic inflammation, effects of which are partially altered in the presence of TBI
Mice fed a chronic HFD diet starting at 12 weeks prior to TBI/Sham surgery exhibited a time-dependent increase in body weight when compared to mice fed on a SD (Fig. 1A). In diet-induced obesity models, mice are considered obese when their body weight is 25 percent greater than that of their counterparts; severe obese models have greater than 50 percent increases in body weight (g) [21]. Beginning at 11 weeks prior to injury, we observed a significant effect of time (F(3.268, 117.6) = 47.77, p < 0.0001; two-way repeated ANOVA; Fig. 1A), treatment (F(30,36) = 281.5, p < 0.0001) and an interaction between time * treatment (F(30,360) = 8.538, p < 0.001) on body weight (g). Post-hoc analysis revealed a significant effect of HFD on body weight (p < 0.05 vs Sham-SD mice; Fig. 1A); at 11- (p < 0.01), 10- (p < 0.001), 9- (p < 0.0001), 8- (p < 0.0001), 7- (p < 0.001), 6- (p < 0.001), 5- (p < 0.01), 4- (p < 0.01), 3- (p < 0.01), 2- (p < 0.001), and 1- (p < 0.001) weeks prior to induction of TBI. After Sham or TBI, mice fed on a HFD showed significant increases in body weight at 0- (p < 0.001), 1- (p < 0.001), 3- (p = 0.000), 7- (p = 0.000), 14- (p = 0.000), 21- (p = 0.000) and 28- (p = 0.000) dpi, compared to their Sham-SD and TBI-SD counterparts (Fig. 1A). TBI alone did not induce significant changes in body weight (g), nor did it influence HFD-induced increases in body weight, with no significant factor effect, interaction, or difference between Sham-HFD vs. TBI-HFD in post hoc analysis (Fig. 1A). In addition, HFD-induced increases in body weight were independent of total food intake in grams (data not shown). Both diet (F(1,14) = 284.2, p < 0.0001) and TBI (F(1,14) = 5.175, p < 0.05) had significant factor effects inducing significant increases in percent serum leptin at 28 dpi (Fig. 1B), but no significant interaction was detected. Post-hoc analysis revealed that Sham-HFD mice had significantly increased percent serum leptin when compared to Sham-SD counterparts (p < 0.0001), at 28 dpi. TBI-SD mice had significantly increased percent leptin when compared to Sham-SD counterparts (p < 0.05) although the magnitude of the increase was much smaller compared to HFD groups. While TBI did not further exacerbate the HFD-induced increases in percent serum leptin, TBI-HFD mice displayed significant increases in percent serum leptin when compared to TBI-SD counterparts (p < 0.0001), at 28 dpi (Fig. 1B).
Statistical analysis revealed that chronic HFD, but not TBI, resulted in increases in resting glucose (F(1,19) = 30.02, p < 0.0001) and insulin (F(1,17) = 19.47, p < 0.001) compared to their SD counterparts at 28 dpi. No significant interaction between diet and TBI factors was detected. Post-hoc analysis determined that Sham-HFD mice had significantly increased resting blood glucose and insulin when compared to their Sham-SD counterparts (p < 0.01; Fig. 1B). In addition, TBI-HFD mice had significantly increased levels of resting blood glucose when compared to TBI-SD mice (p < 0.05).
Statistical analysis revealed that diet (F(1,18) = 11.76, p < 0.01), but not TBI, had a significant effect on serum protein levels of the inflammatory mediator, MCP-1, at 90 dpi. Post-hoc analysis revealed that while HFD alone did not significantly alter serum protein levels of MCP-1, TBI-HFD mice had significantly increased protein levels of MCP-1 when compared to Sham-SD fed counterparts at 90 dpi (p < 0.05) (Fig. 1C).
HFD increases visceral adipose tissue weight and elevates cellular and molecular markers of inflammation; TBI exacerbates selective HFD effects
Chronic HFD induced significant increases in VAT weight (g/g body weight) (F(1,18) = 15.57, p < 0.001), when compared to Sham-SD mice at 28 dpi; no significant TBI factor or interaction effects were detected (Fig. 2A).
Flow cytometry was used to characterize the VAT immune cells environment. Diet but not TBI had significant factor effects on CD45+ (F(1,16) = 8.596, p < 0.01) and CD11b+ (F(1,16) = 14.55, p < 0.01) cell populations at 28 dpi (Fig. 2B–D). Diet had a significant effect on Bead+ CD11b+ ATM/myeloid cell population cell number (F(1,16) = 14.18, p < 0.01; Fig. 2E). In addition, a 2-way ANOVA revealed an effect of diet on % Red Bead+ of CD11b+ ATM/myeloid cells (F(1,18) = 6.664, p < 0.05), while post-hoc analysis revealed that TBI-HFD mice displayed a significant increase in MFI of Red Bead+ of CD11b+ ATM/myeloid cells when compared to TBI-SD counterparts (p < 0.05; Additional file 1: Fig. S1). CD11b− infiltrating lymphocytes were not impacted by diet or TBI (Fig. 2F). Post hoc analysis revealed that while HFD alone did not significantly alter ATM/myeloid cell infiltration, TBI-HFD mice had increased infiltration of CD45+ (p < 0.05), CD11b+ ATM/myeloid (p < 0.01), and Bead+ CD11b+, ATM/myeloid cells (p < 0.05) when compared to Sham-SD counterparts (p < 0.05) at 28 dpi (Fig. 2B, C, E). Reported changes in cell infiltration were not associated with altered intracellular ROS in either CD11b+ or CD11b− cell populations (data not shown).
To further assess the effect of chronic HFD and TBI on the VAT environment, we also used sensitive qRT-PCR measurement of mRNA expression of several molecular markers of inflammatory responses at 28 dpi. Diet had a significant factor effect on expression of IL-1β (F(1,23) = 29.89, p < 0.0001), TNF-α (F(1,23) = 22.57, p < 0.0001), NLRP3 (F(1,23) = 72.88, p < 0.0001), NOX-2 (F(1,21) 25.08 = , p < 0.0001), p22phox (F(1,21) = 64.21, p < 0.0001), and IL-10 (F(1,23) = 31.95, p < 0.0001) at 28 dpi (Fig. 2G). TBI had a significant effect on the expression of IL-1β (F(1,23) = 4.679, p < 0.05) and NLRP3 (F(1,23) = 4.925, p < 0.05). Furthermore, a significant interaction effect (diet * TBI) was observed on expression of IL-1β (F(1,23) = 5.083, p < 0.05), and NLRP3 (F(1,23) = 7.799, p < 0.05). Post hoc analysis revealed that chronic HFD feeding was associated with significant increases in TNF-α (p < 0.05), NLRP3 (p < 0.01), p22phox (p < 0.01), and IL-10 (p < 0.001) when compared to Sham-SD counterparts, at 28 dpi (Fig. 2G). Furthermore, TBI-HFD mice displayed significantly increased expression of IL-1β (p < 0.0001), and NLRP3 (p < 0.0001) when compared to TBI-SD mice, at 28 dpi, while NLRP3 was significantly increased in TBI-HFD mice when compared to Sham-HFD mice (p < 0.01) (Fig. 2G).
HFD and TBI interact to amplify VAT transcriptome changes that are characteristic of macrophage reactive states and inflammation pathways
We extracted the RNA from VAT in the experimental animals in the second cohort at 90 dpi and examined the gene expression changes using the Nanostring Neuroinflammation panel. The Nanostring data were analyzed as described above detecting 483 genes out of a total of 777 with average expression > 50. The general heat map generated by the Morpheus unsupervised hierarchical clustering displays a consistent order in the gradient of differential expression (Sham-SD, TBI-SD, Sham-HFD, TBI-HFD) which in contrast to the brain changes, has diet as the primary factor (Fig. 3A). The cell-phenotype-directed heatmaps highlight the abundance of changes impacting the VAT immune/inflammation responses, representing the majority component among activated genes, that maintain the same clustering order (Fig. 3B). The general heat map based only on the 369 genes with significant ANOVA outcomes display a similar order while bringing into sharper distinction the upregulated inflammation markers (Additional file 1: Fig. S2). Statistical analysis reveals a significant effect of diet on the Z-score sum on adipose macrophage homeostatic (F(1,8) = 56.16, p < 0.0001), adipose macrophage DAM (F(1,8) = 39.55, p < 0.001), adipose tissue function (F(1,8) = 35.34, p < 0.001), adipose inflammatory signaling (F(1,8) = 20.38, p < 0.01), and macrophage pro-inflammatory (F(1,8) = 85.04, p < 0.0001) but not adipose innate immunity (Fig. 4; Additional file 1: Fig. S16A). There was a significant injury effect reported for adipose macrophage function (F(1,8) = 11.30, p < 0.01), adipose innate immunity (F(1,8) = 14.57, p < 0.01) and macrophage pro-inflammatory (F(1,8) = 7.55, p < 0.05). Furthermore, an interaction effect of diet*injury was reported for adipose macrophage homeostatic (F(1,8) = 7.12, p < 0.05) and adipose innate immunity (F(1,8) = 5.81, p < 0.05; Additional file 1: Fig. S16A). Post-hoc analysis revealed a significant increase in the Z-score sum of adipose macrophage homeostatic and adipose macrophage DAM (p < 0.05) as well as pro-inflammatory associated macrophage phenotype (p < 0.01) of Sham-HFD vs Sham-SD diet (Fig. 4). Notably, TBI-HFD mice displayed a significant increase in adipose macrophage homeostatic (p < 0.05), adipose macrophage function (p < 0.05), and adipose innate immunity (p < 0.01) vs Sham-HFD mice. The pathways-centered analysis in IPA demonstrates that while diet (F(1,52) = 1293, p < 0.001) is the dominant factor, the injury (F(1,52) = 119.2, p < 0.0001) factor also reaches significance and that a significant diet*injury interaction (F(1,52) = 90.09, p < 0.0001) is associated with the robust superadditive activation of multiple Canonical pathways(majority are inflammation-related), by combined TBI-HFD (p < 0.0001) versus Sham-HFD (Fig. 5; Additional file 1: Fig. S16A). The Upstream regulators—Genes (Additional file 1: Fig. S3), Upstream regulators—Drugs (Additional file 1: Fig. S4) and Disease pathways (Additional file 1: Fig. S5) also display consistent HFD-driven activation as a principal factor accompanied by injury factor significance and post-hoc tests as well as significant diet*injury interaction associated with amplification of both activated (the majority are inflammation-related) and inhibited pathways by combined TBI-HFD, including upstream genes upregulated (p < 0.0001), upregulated (p < 0.0001) and downregulated (p < 0.0001) upstream drugs, and diseases upregulated and downregulated (p < 0.0001), when compared to Sham-HFD.
TBI increases neuroinflammation in the injured hippocampus and cortex; HFD exacerbates selective TBI effects
We assessed inflammatory gene changes in key brain regions including the ipsilateral (injured hemisphere) hippocampus, and cortex at 28 dpi. In the hippocampus, we detected a significant factor effect of TBI on IL-1β (F(1,20) = 19.02, p < 0.001), NOX-2 (F(1,20) = 5.781, p < 0.05), p22phox (F(1,20) = 16.71, p < 0.001), and GFAP (F(1,20) = 8.408, p < 0.01) (Additional file 1: Fig. S1a). We also observed a significant effect of diet on hippocampal expression of p22phox (F(1,20) = 16.71, p < 0.001), and CD11b (F(1,20) = 6.993, p < 0.05). A significant interaction effect (injury * diet) was noted on hippocampal expression of p22phox (F(1,20) = 12.55, p < 0.01) at 28 dpi. Post-hoc analysis revealed that TBI resulted in a significant increase in hippocampal expression of IL-1β (p < 0.05), when compared to Sham-SD counterparts (Additional file 1: Fig. S6A). Although chronic HFD alone (Sham-HFD), did not significantly alter hippocampal expression of any of the inflammatory genes assessed, post-hoc analysis detected a significant elevation in p22phox in TBI-HFD mice when compared to Sham-HFD (p < 0.001) and TBI-SD (p < 0.001) counterparts (Additional file 1: Fig. S6A). In the cortex, there was a significant effect of TBI, but not HFD, on expression of IL-1β (F(1,18) = 11.05, p < 0.01), TNF-a (F(1,17) = 30.04, p < 0.0001), NOX-2 (F(1,17) = 22.74, p < 0.001), and p22phox (F(1,17) = 32.13, p < 0.0001). No significant interaction (TBI * diet) was detected on cortical gene expression at 28 dpi (Additional file 1: Fig. S6B). Post-hoc analysis revealed that TBI resulted in a significant increase in cortical expression of TNF-α (p < 0.05), NOX-2 (p < 0.05) and p22phox (p < 0.01), but not IL-1β, when compared to Sham-SD counterparts (Additional file 1: Fig. S6B). These injury effects were not altered by chronic HFD.
TBI results in chronic memory/learning deficits; HFD exacerbates elements of cognitive dysfunction
To examine the effects of chronic HFD and TBI on long-term neurological outcomes, a second cohort of mice underwent a battery of neurobehavioral tasks probing cognitive functions up to 90 dpi. In the Y-maze, a task assessing hippocampal-dependent working memory, both TBI (F(1,32) = 9.799, p < 0.01) and diet (F(1,32) = 4.276, p < 0.05) factors had significant effects on % spontaneous alterations, but no significant interaction effect (TBI * diet) was observed, at 70 dpi (Fig. 6A). Post-hoc analysis determined that TBI-SD mice had significantly less % spontaneous alterations when compared with their Sham-SD mice (p < 0.05; Fig. 6A). All groups had equivalent numbers of arm entries in the Y-maze task (Fig. 6A).
Non-spatial hippocampal-mediated memory was assessed using the NOR task at 77–78 dpi (Fig. 6B). During the familiar phase of the task, no significant differences were detected between any of the four groups with regards to time spent with either the right or left familiar object. Twenty-four hours later, mice were retested with a novel object and analysis revealed a significant effect of injury (F(1,31) = 62.41, p < 0.001) and diet (F(1,31) = 14.81, p < 0.001) factors on percent time spent with the novel object, but no significant interaction was detected (Fig. 6B). Post-hoc analysis showed that, compared with Sham-SD mice, TBI-SD mice spent significantly less time with the novel object (42.00 ± 3.27% versus 69.29 ± 3.27%; p < 0.0001). Sham-HFD mice did not spend significantly less time with the novel object compared to the Sham-SD counterparts, but TBI-HFD mice spent significantly less time with the novel object when compared to TBI-SD fed mice (29.55 ± 3.19% versus 42.00 ± 3.27%, p < 0.05), indicating that combined pre-existing chronic HFD and TBI further increases TBI-induced deficits in nonspatial hippocampal-mediated memory.
The MWM was designed based on the concept that the hippocampus creates a spatial cognitive map and is responsible for spatial navigation [47]. Spatial learning and memory were examined using the MWM task at 80–84 dpi. During the acquisition trials, two factor repeated measures ANOVA revealed an effect of time (F(3,140) = 7.309, p < 0.001), treatment (F(9,140) = 17.94, p < 0.001), but no interaction effect between time (acquisition day) * treatment group on latency (s) to the escape platform (Fig. 6C). Twenty-four hours later, retention memory was assessed using the probe trial. There was a significant effect of injury (F(1,33) = 29.69, p < 0.0001) and diet (F(1,33) = 10.17, p < 0.01) factors as well as a significant interaction injury * diet (F(1,33) = 7.59, p < 0.01) on the time to the hidden platform (Fig. 6D). Post-hoc analysis revealed that TBI-SD (p < 0.0001), Sham-HFD (p < 0.01) and TBI-HFD (p < 0.0001) had significantly increased latency to the escape platform compared to Sham-SD counterparts. No significant differences were observed between TBI-HFD and TBI-SD (Fig. 6D). There was no significant effect of diet or TBI factors on swim speed (Fig. 6E). The swim pattern used to find the hidden platform was analyzed and assigned a search strategy, classified according to the type of learning taking place, with the least precise/efficient strategies being repetitive looping, followed by non-spatial systematic and the more precise strategies being spatial as described previously [8]. Search strategy analysis permits differentiation between hippocampus-dependent spatial (allocentric) and hippocampus-independent non-spatial (egocentric) strategies including systematic and repetitive looping search strategies [48]. A spatial pattern indicates an intact hippocampus and adequate spatial memory, whereas non-spatial systematic and repetitive looping patterns indicate increasing hippocampal dysfunction and growing spatial memory impairments, respectively [48]. Sham-SD mice demonstrated proper memory while increasing memory deficits were noted in TBI-SD, Sham-HFD and TBI-HFD (p < 0.05, x = 14.60) (Fig. 6F, G).
TBI causes brain microglia transcriptome changes driving increased reactive states and decreased homeostatic phenotype; HFD amplifies posttraumatic activation of inflammation pathways
To identify microglia transcriptional changes, we performed CD11b-positive cells selection (using Miltenyi MACS technology; microglia are the principal component of this cellular fraction) from the cortex and hippocampus in the second cohort at 90 dpi followed by high-throughput gene expression examination using Nanostring panels. RNA was extracted and analyzed using the nCounter Mouse Glial panel (Nanostring Technologies) that profiles more than 750 genes related to glia functions and metabolism. Gene expression data was normalized using the Rosalind Nanostring analysis package and all further analyses were limited to genes whose average expression across all experimental groups was at least 50, to reduce the impact of large variations common in low expressor genes and increase biological significance. In the microglia fraction, 526 genes out of a total of 757 genes met the expression requirement and underwent a further Z-score normalization that permits examination of expression changes independently of absolute expression values. The per-gene Z-scores list was analyzed using the Morpheus heatmap program, an unsupervised hierarchical clustering package (https://software.broadinstitute.org/morpheus) [49]. The generated heat map shows that while TBI is the primary determining factor for gene expression changes with HFD having a lesser impact; a consistent gradient of changes is observed with the following group order: Sham-SD, Sham-HFD, TBI-SD, TBI-HFD (Fig. 7A). The gene function description is based on Nanostring annotations that were updated based on recent data. The general heat map (Fig. 7A) and the microglia phenotype-specific heat maps (Fig. 7B) identifies homeostatic microglia markers along with other microglia markers as among the most consistently downregulated phenotypes while the disease-associated microglia (DAM) and ‘M1’-like micorglia which are associated with the classical pro-inflammatory markers (pro-inflammatory-associated microglia) [50], were among the most consistently upregulated phenotypes. Markers of autophagy-phagocytosis and interferon signaling were also elevated (Fig. 7A, B).
To quantify the gene expression changes reflected by the heatmaps we summated in each sample the Z-scores that corresponded to genes defining specific cellular phenotypes. A baseline Sham-SD vs Sham-SD comparison with a z-score “0” was generated across pathways. The generated Z-score sum is an overall indicator of the relative change in expression of the phenotype-specific genes in each sample and the values across groups were analyzed using two-way ANOVA (Fig. 8). Statistical analysis revealed a significant effect of injury on homeostatic-associated microglia (F(1,8) = 39.63, p < 0.0001), autophagic-associated microglia (F(1,8) = 78.60, p < 0.0001), interferon-associated microglia (F(1,8) = 11.86, p < 0.01), pro-inflammatory-associated microglia (F(1,8) = 27.90, p < 0.001), DAM (F(1,8) = 225.4, p < 0.0001) and microglia markers (F(1,8) = 34.43, p < 0.001) (Fig. 8; Additional file 1: Fig. S6B). No significant effect of diet or diet*injury interaction was observed, except a modest diet effect in microglia markers. Post-hoc analysis revealed that TBI-SD mice had a significant decrease in homeostatic-associated microglia (p < 0.01) and microglial markers (p < 0.01) while displaying a significant increase in autophagic-associated microglia (p < 0.01), pro-inflammatory-associated microglia (p < 0.01), and DAM (p < 0.0001), when compared to Sham-SD counterparts. TBI-HFD showed similar changes compared to Sham-HFD. No significant changes were observed in Sham-HFD vs Sham-SD or TBI-HFD vs TBI-SD confirming that the injury (TBI) factor but not the diet (HFD) factor was associated with significant F ratios and p value in post-hoc comparisons (Fig. 8).
We also analyzed the normalized expression levels (all genes with average expression > 50) using a two-way ANOVA. Genes that demonstrated a significant p value (p ≤ 0.05) in at least one factor, interaction or post-hoc test were selected for further analysis (n = 312 genes). The Morpheus-generated unsupervised hierarchical clustering heat map (Additional file 1: Fig. S7; annotations include all ANOVA outcomes) demonstrates the same patterns as in Fig. 7A, although gene expression patterns are easier to detect due to the elimination of genes that do not significantly change with injury/diet.
Ingenuity Pathway Analysis (IPA) software (QIAGEN) is an advanced bioinformatics tool, based on built-in scientific literature databases (www.ingenuity.com) that identifies the most impacted molecular pathways and predicts cellular responses. For the 312 genes with significant ANOVA results we imported into IPA the Rosalind-generated log2FC of Sham-HFD vs Sham-SD, TBI-SD- vs Sham-SD and TBI-HFD vs Sham-SD and the corresponding post-hoc p values generated by the repeated two-way ANOVA analysis. We determined the activation of IPA canonical pathways based on the inclusion of genes with p ≤ 0.05 and pathway z-score > 2 or < − 2 (the minimal requirements for IPA to conclude that a pathway is activated or inhibited, respectively) (Fig. 9). The data suggest that TBI-HFD tends to cause more pathway changes including stronger activation of multiple inflammatory pathways. We also quantified these changes after separating the activated and inhibited pathways (based on the outcome in the TBI-HFD group) and analyzing the pathway activity z-scores with repeated two-way ANOVA. To assist analysis, we introduced a baseline Sham-SD vs Sham-SD comparison with a z-score “0” across all pathways and replaced any IPA-generated N/A classification (not able to detect activation level) with a score of “0”. Statistical analysis shows injury (F(1,16) = 199.8, p < 0.0001), diet (F(1,16) = 10.29, p < 0.01), and an injury*diet interaction (F(1,16) = 10.29, p < 0.01) associated with activation/upregulation of a group of Canonical pathways mostly composed by pro-inflammatory pathways (Fig. 9; Additional file 1: Fig. S16B). For the less numerous Canonical pathways that were downregulated/inhibited, significance was only observed for an injury effect (F(1,4) = 458.9, p < 0.0001). Post-hoc analysis revealed that TBI-SD mice had significantly increased Z-score (p < 0.0001) vs Sham-SD mice. TBI-HFD showed similar changes when compared to Sham-HFD. Sham-HFD was not associated with pathway activation or inhibition. Notably, TBI-HFD mice displayed significantly increased Z-scores (p = 0.0017) vs TBI-SD, for canonical upregulated pathways (Fig. 9). For the downregulated pathways, post-hoc analysis revealed that TBI-SD mice had a significantly decreased Z-score (p < 0.0001) vs Sham-SD and that TBI-HFD showed similar changes when compared to Sham-HFD.
We also examined changes in IPA Upstream regulators—Genes (Additional file 1: Fig. S8), Upstream regulators—Drugs (Additional file 1: Fig. S9) and Diseases (Additional file 1: Fig. S10). The injury factor and related post-hoc tests were significant in all cases, but significant diet factor was only observed in Upstream Genes and Drugs and significant interaction effects with significant TBI-HFD vs TBI-SD post-hoc test was only observed in Upstream Drugs upregulated pathways which included several inflammatory pathways (e.g., lipopolysaccharide (LPS), polyinosinic:polycytidylic (polyI-C), etc.)).
TBI is the primary factor driving changes in non-microglia cellular transcriptomes
To determine the specificity of the CD11b-positive selection and the resulting selective separation of microglia and CD11b-negative (non-microglia) populations (neurons, astrocytes, oligodendrocytes, etc.), we averaged the gene expression levels across microglia and non-microglia fraction in each experimental group followed by Z-score transformation. The ratio between averaged microglia and flowthrough normalized expression data (560 genes out of a total of 769 genes had average expression > 50 in the flowthrough) was used to determine the top 50 overexpressed and under expressed genes in microglia vs flowthrough. The Z-scores for these genes were uploaded into Morpheus for unsupervised hierarchical clustering (Additional file 1: Fig. S11). The data confirm the quality of the separation as the CD11b-positive fraction demonstrates over 100-fold enrichment of key microglia markers. The reverse is true for the flowthrough in which neuronal, astrocytes and oligodendrocytes markers are highly enriched. Notably, even at this reduced scale the gradient of injury/diet changes is maintained and a decline in the microglia homeostatic phenotype and neurogenesis as well as increase in microglia and astrocyte altered phenotypes after injury/HFD are observed in the appropriate compartments.
Performing the Morpheus heat map on all 560 genes in the flow-through did not result in an ordered clustering across the groups highlighting the smaller number of genes with consistent and significant changes (Additional file 1: Fig. S12). Conversely, when a further selection restricted the list to genes demonstrating a p ≤ 0.05 in at least one of the two-way ANOVA outcomes (as described above; n = 121 genes; Fig. 10) the resulting heat map displayed a consistent gradient of changes across groups that matched the clustering observed in the microglia (Sham-SD, Sham-HFD, TBI-SD, TBI-HFD). Genes downregulated included neurogenesis and neuronal signaling markers while genes upregulated included astrocyte activation phenotype markers (Fig. 10; Additional file 1: Fig. S13 includes the heatmap with ANOVA data). The generated Z-score sum demonstrates that injury is the only factor driving significant downregulation of neuronal markers (F(1,8) = 30.95, p < 0.001), an upregulation of astrocytes (F(1,8) = 18.79, p < 0.01) and oligodendrocytes markers (F(1,8) = 23.60, p < 0.01); injury*diet interaction was significant for neuronal (F(1,8) = 8.78, p < 0.05) and oligodendrocyte (F(1,8) = 5.41, p < 0.05) markers but it was indicative of absence of additive effects as post-hoc statistical analysis revealed that TBI-SD (p < 0.01) and TBI-HFD (p < 0.01) had similar effects vs. Sham-SD while diet alone (Sham-HFD) had effects in the same direction as injury although these reached significance only in the case of neuronal markers (Fig. 10; Additional file 1: Fig. S16C). These changes may in part reflect declines in neuronal populations and increases in astrocytic and oligodendrocytes populations. The IPA examination of molecular responses failed to detect any Canonical pathways with significant changes (z > 2 or z < − 2). The few activated Upstream genes pathways showed significance only for the TBI factor and no significant interaction (Additional file 1: Fig. S14). The Diseases pathways displayed multiple inhibited or activated genes; the majority of the activated pathways were inflammatory with only TBI as a significant factor and a modest but significant increase in this response was observed in the TBI-HFD group vs. TBI-SD (Additional file 1: Fig. S15).
HFD and TBI contribute to global separation of cellular phenotypes
We used principal component analysis, a linear dimensionality reduction technique to highlight the transcriptomic differences between specific groups. The gene expression changes for every cell fraction, based the Z-scores of genes that have been identified as differentially expressed in Rosalind, were analyzed by PCA. Consistent with previous data, the microglia and flowthrough show TBI as a primary factor for changes with HFD playing a secondary role, whereas the reverse was true for the VAT. Notably, when PCA was applied to an integrated data set (brain microglia and flowthrough as well as VAT) for every experimental animal, we generated a composite score that resulted in robust separation of all groups with combined TBI-HFD located at the union of individual factors (TBI and HFD) directions of change and most distinct from Sham-SD (Fig. 11).
Discussion
The present study examined the individual and combined effects of HFD and TBI on VAT, systemic circulation, and brain, with a specific focus on cellular transcriptomic programs and cognitive changes. At 28 dpi (cohort 1), HFD is the major driver of time-dependent changes in body weight, metabolic markers, VAT weight, and VAT immune cell activation. Although HFD alone was not associated with a significant increase in ATM/myeloid cell numbers, HFD in the presence of TBI modified the VAT immune environment with significantly increased numbers of CD45+ and CD11b+ ATM/myeloid cells at 28 dpi; effects of which were associated with increased ATM/myeloid cell phagocytic activity. Furthermore, HFD increased VAT expression of the pro-inflammatory mediators TNF-a, NLRP3, NOX-2, and p22phox as well as the anti-inflammatory mediator IL-10. Notably, HFD in the presence of TBI showed significant diet*injury interactions with exacerbation of HFD-dependent increases in VAT pro-inflammatory molecules IL-1β and NLRP3, indicating a brain trauma-dependent amplification of diet-induced adipose tissue inflammation. Systemically, HFD resulted in evident metabolic dysfunction including increased leptin, insulin, and resting glucose levels; effects of which were independent of TBI, except for evident TBI-induced increases in circulating leptin. Notably, circulating levels of the inflammatory mediator, MCP-1, was significantly increased in the presence of -TBI-HFD. At the level of the brain, although TBI was shown to be the primary driver of microglial inflammatory changes; HFD exacerbated, at least in part, TBI-induced deficits in cognitive function. Overall, this study demonstrates a potential bi-directional neuroimmune relationship between the adipose tissue and brain in the presence of co-morbid diet-induced obesity and TBI; that may be targetable.
ATM are a diverse population of cells with distinct developmental lineages including yolk-sac derived and self-maintaining tissue resident macrophages as well as macrophages resulting from recruitment of bone-marrow derived monocytes, a process that is accelerated in the adipose tissue of obese animals [51]. Monocytes infiltrate the adipose tissue [52, 53] and differentiate in macrophages beginning at 8 weeks following initiation of diet-induced obesity [54, 55] and these bone-marrow (myeloid) derived monocytes/macrophages may be primarily involved in adipose tissue inflammation [51]. Infiltrating ATM secrete pro-inflammatory mediators including TNF-α, IL-1β, MCP-1, and IL-6 which are suggested to play an important role in driving a pro-inflammatory state within the VAT microenvironment and subsequent low grade chronic inflammation in the systemic circulation [56,57,58,59]. Other myeloid immune cells, namely neutrophils, have been shown to transiently infiltrate the adipose tissue as early as 3 days following initiation of HFD [60] and are identified as key drivers of obesity-induced increases in IL-1β [61]. In the present study, we demonstrate that HFD has a significant factor effect to elevate CD45+ and CD11b+ cell numbers in the VAT microenvironment and that HFD in the presence of TBI, results in a significant increase in the numbers of these cell populations at 28 dpi; these effects are associated with increased phagocytic activity. While we cannot distinguish the specific functions of the various lineages of the CD45+/CD11b+ cells including the contribution of yolk-sac vs. myeloid infiltrating ATM/monocytes or that of other infiltrating myeloid cells such as neutrophils, our data suggest that overall, the ATM/myeloid populations contribute to driving the reported inflammatory changes in the adipose tissue microenvironment.
To further explore the VAT immune cell microenvironment, we performed nanostring analysis on the VAT at 90 dpi. The data show an ordered clustering of differential gene expression, which establishes diet as the primary factor driving changes. Disease associated macrophages (DAM) molecules, as well as markers of various pro-inflammatory states, likely expressed in ATM, are among the most numerous HFD/TBI upregulated genes, whereas other immune and epigenetic regulators are downregulated. Importantly, the VAT shows a significant effect of both TBI and HFD, and/or the presence of significant interactions driving super-additive pro-inflammatory effects. The TBI/HFD-dependent elevation was evident not only for the inflammatory states but also for the homeostatic state, suggesting a complex response in which HFD increases some macrophage sub-populations, whereas the combined TBI/HFD increases all macrophage sub-populations, typically at higher levels than HFD alone. IPA pathway analysis extends these findings, identifying the inflammatory pathways as the largest component of upregulated pathways, with both diet and injury as significant factors (although diet remains the dominant factor), and a significant diet-injury interaction driving a synergistic increase in the pathway activation levels after combined TBI-HFD compared to Sham-HFD.
A recent report provides evidence that ATM accumulate lipids in obese mice and become polarized toward the lipid-associated macrophages (LAM) state [62]. LAM development is driven by Trem2 and characterized by upregulation of components of an enzymatic machinery that recognizes, scavenges, and catabolizes lipids such as fatty acids transporter Cd36, fatty acid binding proteins 5 (Fabp5) and lipoprotein lipase (Lpl). Notably, these and other molecules defining LAM establish a profile that is virtually identical with the conserved signature for DAM [63]. The LAM/DAM signature in VAT samples is elevated by HFD and appears to be a Trem2-dependent uniform response and responsible, at least initially, for the clearance of injured cells and damaged tissue, resolution of inflammation, and improvement of metabolic changes in obesity [64, 65]. However, ATM also transform to pro-inflammatory phenotypes, including pro-inflammatory associated genes, changes that may underlie the disease-driving role played by ATM in metabolic conditions and obesity-associated inflammation [66].
Circulating levels of MCP-1 are increased in animal models of diet-induced obesity [67, 68] where it appears to play a role in HFD-induced insulin resistance, glucose intolerance [52, 69] and ATM infiltration [52, 69, 70]. In the present study, we show that whereas HFD alone does not significantly increase serum levels of MCP-1, HFD in the presence of TBI causes a significant increase in serum levels of MCP-1 at 90 dpi. This observation may have significant implications, given that modulating circulating MCP-1 signaling reduces microglial activation and improves neurological recovery in hepatic encephalopathy [71]. Furthermore, limiting MCP-1 signaling suppresses LPS-induced seizures [72]. The adipokine leptin is another factor reported to play a key role in the communication between peripheral immune responses and the and brain [73, 74]. Notably, both HFD and TBI, albeit to a lesser extent, increase circulating levels of leptin at 28 dpi, suggesting a potential role for HFD-induced increases in circulating leptin in driving brain changes.
Diet-induced obesity is associated with brain neuroinflammation [16], including in the hypothalamus [75, 76]. There are conflicting reports regarding the effect of HFD and obesity on extra-hypothalamic inflammation [21, 37, 77,78,79], with some studies reporting no change [78, 79] and others identifying neuroinflammation in the hippocampus and cortex [37, 77]. Isolated cortical microglia from mice fed on a HFD regime release more TNF-α compared to SD-fed counterparts [80]. Diet-induced obesity, unlike short-term HFD feeding, leads to altered hippocampal microglial states, reductions in dendritic spines at sites of excitatory synapses and memory impairments; the cognitive deficits were reversed by partial knockdown of microglia and associated reduction in microglial phagocytosis [21]. In TBI patients and experimental models, microglia undergo a chronic transformation to reactive phenotypes that are associated with neurodegenerative processes and cognitive decline [10, 11]. Delayed transient removal of microglia after experimental TBI reduced the inflammatory lesion microenvironment and improved long-term neurological recovery [8]. Although previous studies have reported increased microglial altered states in the hypothalamus and the prefrontal cortex in the context of co-morbid diet-induced obesity and TBI [23, 26], the present study is the first to investigate hippocampal and cortical inflammatory responses at 28dpi (cohort 1). Unlike TBI, HFD significantly altered only two pro-inflammatory genes (p22phox and CD11b) in the hippocampus. However, combined TBI-HFD displayed a significant diet*injury interaction, including an elevation in hippocampal p22phox compared to TBI-SD mice. These results suggest that HFD can prime development of TBI-induced, pro-inflammatory responses.
To better characterize the interaction between HFD and TBI on secondary posttraumatic injury, we examined cognitive (memory/learning) functions as well as large-scale transcriptomic changes in specific cellular compartments after chronic brain trauma as late as 90 dpi (cohort 2). Y-Maze and NOR tests showed significant effects for both TBI and HFD, although no significant diet*injury interaction was observed, suggesting a simple additive association. In contrast, HFD and TBI had significant factor effects in MWM but not an additive one. However, combined TBI-HFD did show additive deficits in the more sensitive search strategy measurement in the MWM test.
Transcriptome signatures have been proposed to define the complex and dynamic microglial states [81, 82]. To examine the impact of diet and injury on the molecular pathways potentially related to neurodegeneration and neurological dysfunction, we isolated and analyzed transcriptomic profiles of microglial populations from the perilesional cortex and hippocampus. The heatmap of microglia gene differential expression at 90 dpi displayed a consistent and ordered clustering of the experimental groups: Sham-SD, Sham-HFD, TBI-SD and TBI-HFD. The fact that Sham-HFD diet is consistently distinct from Sham-SD suggests that, at least at a global level, HFD influences microglia responses. The secondary grouping of Sham-SD with Sham-HFD and TBI-SD with TBI-HFD demonstrate that TBI is the dominant factor responsible for change, with HFD having a secondary, albeit persistent, impact. Major components of differential expression in response to TBI and/or HFD include downregulation of homeostatic microglia that occur in response to pathological conditions (injury and/or disease) [81, 83] and may serve an adaptive role when acting to promote repair, but may also represent a maladaptive response when associated with dysfunctional proinflammatory states that perpetuate tissue damage [83]. The transcriptomic signature of the DAM and related microglia neurodegenerative phenotype (MGnD) state represents one of the most consistent differentially expressed profiles. However, as with other microglia pro-inflammatory signatures it reflected TBI, but not HFD, as a significant factor. The collective magnitude for overexpression of virtually all DAM/MGnD markers- including ApoE, Spp1, Trem2, Axl, Clec7a and Lpl [81, 83]- is paralleled by opposite directional changes for homeostatic markers, indicating a strong shift of microglial populations toward the DAM reactive phenotype after TBI/HFD. Downregulation of the homeostatic signature is one of the most common microglial changes shared across various neurodegenerative conditions [81]. The DAM/MGnD state, driven by the Trem2-APOE signaling pathway [84], is found in many neurodegenerative disorders, but its function in this regard remains to be elucidated [83]. Trem2, which along with APOE is elevated in microglia after TBI in our study, plays a key role in microglial responses to pathological signals such as amyloid accumulation, promoting their exit from the homeostatic state and shift to the DAM reactive state [85], especially DAM stage 2 [83]. Trem2 may be commonly required for microglia to switch from the homeostatic state to all reactive states [83]. Conversely, Trem2 absence impairs phagocytosis of key substrates, including APOE, resulting in an impaired response to beta-amyloid; thus, in certain circumstances the shift of Trem2-dependent microglia to the DAM state may reflect an adaptive response [85]. In contrast to the differential expression approach, which identified only TBI as a significant factor, IPA analysis using a pathway-centric approach detected a significant effect for TBI and HFD, as well as the significant interaction underlying the super-additive activation of multiple microglial reactive inflammatory pathways in TBI-HFD compared to TBI-SD.
Overall, our data indicate that in the more chronic phase (90 dpi), microglia remain strongly shifted toward heterogenous reactive states/responses primarily driven by TBI, but with HFD having an effect. Thus microglia, due to inability to promote effective repair, especially in the presence of an additional insult (HFD), enter persistent altered states that are maladaptive and may propagate secondary injury.
The gene expression changes in the non-microglia populations also show TBI-dependent down-regulation of neuronal responses and neurogenesis, as well as upregulation of astrocyte and oligodendrocyte signatures, that are consistent with neurodegeneration [86, 87]. Notably, HFD is also able to decrease neuronal signature. Pro-inflammatory changes were also found in non-microglia populations, suggesting that microglia are not the only cell population involved in post-traumatic inflammatory responses. Overall, principal component analysis (PCA) illustrates the predominant importance of injury for brain changes and diet for VAT changes, but also shows clear separation of responses across the various experimental groups when a composite score is generated that includes both brain and VAT changes.
Evidence demonstrates that pre-existing obesity is associated with worsened outcomes following TBI [23,24,25,26]. A prior study implicated a pivotal role for VAT NLRP3-induced signaling in driving obesity-induced cognitive deficits; these effects reflected IL-1-mediated microglial activation [36]. Furthermore, an important role for the IL-1/IL-1R1 signaling axis in driving TBI-induced cognitive deficits has been reported [88, 89]. Thus, supporting the notion that targeting the NLRP3/IL-1 signaling axis may be a feasible target in the context of co-morbid TBI-HFD. Notably, the present study demonstrates an interaction effect of TBI-HFD on VAT expression of both NLRP3 and IL-1β mRNA expression and a combined HFD-TBI effect on levels of systemic inflammatory mediators. These findings support the conclusion that HFD-TBI induces a reinforcing loop in which pro-inflammatory molecules released by VAT into the circulation promote a ‘primed’ state within the brain, wherein a subsequent insult (i.e. TBI) serves to exacerbate neuroinflammatory responses and related cognitive decline.
However, there are important limitations in the present studies that should be recognized:
-
1.
HFD was initiated prior to TBI and continued after injury. Thus, it is not possible to distinguish between the pre- and post-injury effects of diet. Furthermore, HFD and resultant obesity are not equivalent because HFD is maintained throughout the experiment, even after obesity was established. Therefore, we cannot distinguish the effects of continuous HFD on the VAT (and indirectly on the brain) versus more direct actions on the brain.
-
2.
The Nanostring transcriptomic analysis, although including more than 750 genes, is less comprehensive than RNAseq. More importantly, the panels used are highly enriched in inflammatory molecules, and as such it is not possible to apply gene set enrichment analysis, which restricted our IPA work to pathway activation score analysis.
-
3.
Sampling limitations include the probable presence of brain-resident macrophages and infiltrating monocytes as a smaller component in the CD11b positively selected fraction, the inclusion of combined injured hippocampus and cortex regions for cell isolation, and the absence of cell-specific isolation (macrophages) in the VAT. In addition, the presence of heterogeneous cellular populations in samples such as brain flowthrough (neuronal, astrocytic, oligodendrocytes, etc.) and adipose tissue (adipocytes, macrophages, etc.) complicate the analysis and may represent an important confounding factor when interpreting reported molecular/cellular changes in this data set. Future studies should also examine specific cell populations in the adipose tissue, including adipocytes, ATM, and neutrophils.
-
4.
This study focused on high throughput transcriptomic analysis at level of the adipose tissue and brain. Detailed proteomic analysis was not performed; therefore, the biological relevance of identified gene expression changes remain speculative. Future studies should include protein changes related to key identified pathways in both the brain and adipose tissue, effects of modulating such pathways, and elucidating more detailed region-specific brain changes.
In summary, our findings further support the concept of bi-directional communications between the adipose tissue and brain, which may contribute to brain injury and be potentially targetable. In-depth transcriptomic analysis at both the tissue and cellular level revealed that HFD enhances TBI-induced brain pro-inflammatory responses and that TBI augments HFD-induced VAT inflammatory pathway changes. We hypothesize that this interaction creates a self-propagating secondary inflammatory injury loop with an exacerbated microglial response, which contributes to chronic posttraumatic disability including cognitive decline.
Availability of data and materials
Data will be made available upon request.
Abbreviations
- TBI:
-
Traumatic brain injury
- CCI:
-
Controlled cortical impact
- DIO:
-
Diet-induced obesity
- HFD:
-
High fat diet
- SD:
-
Standard diet
- VAT:
-
Visceral adipose tissue
- DAM:
-
Disease associated microglia
- MWM:
-
Morris water maze
- NOR:
-
Novel object recognition
References
Corrigan JD, Selassie AW, Orman JA. The epidemiology of traumatic brain injury. J Head Trauma Rehabil. 2010;25(2):72–80.
Chooi YC, Ding C, Magkos F. The epidemiology of obesity. Metabolism. 2019;92:6–10.
Winfield RD, Bochicchio GV. The critically injured obese patient: a review and a look ahead. J Am Coll Surg. 2013;216(6):1193–206.
Chabok SY, et al. The impact of body mass index on treatment outcomes among traumatic brain injury patients in intensive care units. Eur J Trauma Emerg Surg. 2014;40(1):51–5.
Christmas AB, et al. Morbid obesity impacts mortality in blunt trauma. Am Surg. 2007;73(11):1122–5.
Ditillo M, et al. Morbid obesity predisposes trauma patients to worse outcomes: a National Trauma Data Bank analysis. J Trauma Acute Care Surg. 2014;76(1):176–9.
Mishra R, et al. Obesity as a predictor of outcome following traumatic brain injury: a systematic review and meta-analysis. Clin Neurol Neurosurg. 2022;217: 107260.
Henry RJ, et al. Microglial depletion with CSF1R inhibitor during chronic phase of experimental traumatic brain injury reduces neurodegeneration and neurological deficits. J Neurosci. 2020;40(14):2960–74.
Johnson VE, et al. Inflammation and white matter degeneration persist for years after a single traumatic brain injury. Brain. 2013;136(Pt 1):28–42.
Loane DJ, et al. Progressive neurodegeneration after experimental brain trauma: association with chronic microglial activation. J Neuropathol Exp Neurol. 2014;73(1):14–29.
Simon DW, et al. The far-reaching scope of neuroinflammation after traumatic brain injury. Nat Rev Neurol. 2017;13(3):171–91.
Chong AJ, et al. Effects of a high-fat diet on neuroinflammation and apoptosis in acute stage after moderate traumatic brain injury in rats. Neurocrit Care. 2020;33(1):230–40.
Aguilar-Valles A, et al. Obesity, adipokines and neuroinflammation. Neuropharmacology. 2015;96(Pt A):124–34.
Bousquet M, et al. High-fat diet exacerbates MPTP-induced dopaminergic degeneration in mice. Neurobiol Dis. 2012;45(1):529–38.
Bruce-Keller AJ, Keller JN, Morrison CD. Obesity and vulnerability of the CNS. Biochim Biophys Acta. 2009;1792(5):395–400.
Miller AA, Spencer SJ. Obesity and neuroinflammation: a pathway to cognitive impairment. Brain Behav Immun. 2014;42:10–21.
Purkayastha S, Cai D. Neuroinflammatory basis of metabolic syndrome. Mol Metab. 2013;2(4):356–63.
Campillo BW, et al. Short-term high-fat diet alters the mouse brain magnetic resonance imaging parameters consistently with neuroinflammation on males and metabolic rearrangements on females. A pre-clinical study with an optimized selection of linear mixed-effects models. Front Neurosci. 2022;16:1025108.
Valdearcos M, et al. Microglial inflammatory signaling orchestrates the hypothalamic immune response to dietary excess and mediates obesity susceptibility. Cell Metab. 2017;26(1):185-197.e3.
Valdearcos M, et al. Microglia dictate the impact of saturated fat consumption on hypothalamic inflammation and neuronal function. Cell Rep. 2014;9(6):2124–38.
Cope EC, et al. Microglia play an active role in obesity-associated cognitive decline. J Neurosci. 2018;38(41):8889–904.
Henn RE, et al. Obesity-induced neuroinflammation and cognitive impairment in young adult versus middle-aged mice. Immun Ageing. 2022;19(1):67.
Sherman M, et al. Adult obese mice suffer from chronic secondary brain injury after mild TBI. J Neuroinflammation. 2016;13(1):171.
Thomson S, et al. Impact of high fat consumption on neurological functions after traumatic brain injury in rats. J Neurotrauma. 2022;39(21–22):1547–60.
Ibeh S, et al. High fat diet exacerbates long-term metabolic, neuropathological, and behavioral derangements in an experimental mouse model of traumatic brain injury. Life Sci. 2023;314: 121316.
Karelina K, et al. Traumatic brain injury and obesity induce persistent central insulin resistance. Eur J Neurosci. 2016;43(8):1034–43.
Byrnes KR, et al. Delayed mGluR5 activation limits neuroinflammation and neurodegeneration after traumatic brain injury. J Neuroinflammation. 2012;9:43.
Piao CS, et al. Late exercise reduces neuroinflammation and cognitive dysfunction after traumatic brain injury. Neurobiol Dis. 2013;54:252–63.
Catrysse L, van Loo G. Adipose tissue macrophages and their polarization in health and obesity. Cell Immunol. 2018;330:114–9.
Kane H, Lynch L. Innate immune control of adipose tissue homeostasis. Trends Immunol. 2019;40(9):857–72.
Debette S, et al. Midlife vascular risk factor exposure accelerates structural brain aging and cognitive decline. Neurology. 2011;77(5):461–8.
Fitzpatrick AL, et al. Midlife and late-life obesity and the risk of dementia: cardiovascular health study. Arch Neurol. 2009;66(3):336–42.
Whitmer RA, et al. Body mass index in midlife and risk of Alzheimer disease and vascular dementia. Curr Alzheimer Res. 2007;4(2):103–9.
Fourrier C, et al. Brain tumor necrosis factor-alpha mediates anxiety-like behavior in a mouse model of severe obesity. Brain Behav Immun. 2019;77:25–36.
Gaudet AD, et al. miR-155 deletion in female mice prevents diet-induced obesity. Sci Rep. 2016;6:22862.
Guo DH, et al. Visceral adipose NLRP3 impairs cognition in obesity via IL-1R1 on CX3CR1+ cells. J Clin Invest. 2020;130(4):1961–76.
Pepping JK, et al. Myeloid-specific deletion of NOX2 prevents the metabolic and neurologic consequences of high fat diet. PLoS ONE. 2017;12(8): e0181500.
Vandanmagsar B, et al. The NLRP3 inflammasome instigates obesity-induced inflammation and insulin resistance. Nat Med. 2011;17(2):179–88.
Barrett JP, et al. Interferon-beta plays a detrimental role in experimental traumatic brain injury by enhancing neuroinflammation that drives chronic neurodegeneration. J Neurosci. 2020;40(11):2357–70.
Barrett JP, et al. NOX2 deficiency alters macrophage phenotype through an IL-10/STAT3 dependent mechanism: implications for traumatic brain injury. J Neuroinflammation. 2017;14(1):65.
Henry RJ, et al. Inhibition of miR-155 limits neuroinflammation and improves functional recovery after experimental traumatic brain injury in mice. Neurotherapeutics. 2019;16(1):216–30.
Henry RJ, et al. Longitudinal assessment of sensorimotor function after controlled cortical impact in mice: comparison of beamwalk, rotarod, and automated gait analysis tests. J Neurotrauma. 2020;37(24):2709–17.
Andrade C. Z scores, standard scores, and composite test scores explained. Indian J Psychol Med. 2021;43(6):555–7.
Pavlidis P. Using ANOVA for gene selection from microarray studies of the nervous system. Methods. 2003;31(4):282–9.
Marini F, Binder H. pcaExplorer: an R/Bioconductor package for interacting with RNA-seq principal components. BMC Bioinformatics. 2019;20(1):331.
Toussaint PA, et al. Explainable artificial intelligence for omics data: a systematic mapping study. Brief Bioinform. 2023. https://doi.org/10.1093/bib/bbad453.
Rogers J, et al. Search strategy selection in the Morris water maze indicates allocentric map formation during learning that underpins spatial memory formation. Neurobiol Learn Mem. 2017;139:37–49.
Curdt N, et al. Search strategy analysis of Tg4-42 Alzheimer Mice in the Morris Water Maze reveals early spatial navigation deficits. Sci Rep. 2022;12(1):5451.
Chun YS, et al. MEK1/2 inhibition rescues neurodegeneration by TFEB-mediated activation of autophagic lysosomal function in a model of Alzheimer’s Disease. Mol Psychiatry. 2022;27(11):4770–80.
Colonna M, Butovsky O. Microglia function in the central nervous system during health and neurodegeneration. Annu Rev Immunol. 2017;35:441–68.
Cox N, Geissmann F. Macrophage ontogeny in the control of adipose tissue biology. Curr Opin Immunol. 2020;62:1–8.
Cranford TL, et al. Role of MCP-1 on inflammatory processes and metabolic dysfunction following high-fat feedings in the FVB/N strain. Int J Obes (Lond). 2016;40(5):844–51.
Pepping JK, et al. NOX2 deficiency attenuates markers of adiposopathy and brain injury induced by high-fat diet. Am J Physiol Endocrinol Metab. 2013;304(4):E392-404.
Strissel KJ, et al. Adipocyte death, adipose tissue remodeling, and obesity complications. Diabetes. 2007;56(12):2910–8.
Weisberg SP, et al. Obesity is associated with macrophage accumulation in adipose tissue. J Clin Invest. 2003;112(12):1796–808.
Makki K, Froguel P, Wolowczuk I. Adipose tissue in obesity-related inflammation and insulin resistance: cells, cytokines, and chemokines. ISRN Inflamm. 2013;2013: 139239.
Russo L, Lumeng CN. Properties and functions of adipose tissue macrophages in obesity. Immunology. 2018;155(4):407–17.
Salas-Venegas V, et al. The obese brain: mechanisms of systemic and local inflammation, and interventions to reverse the cognitive deficit. Front Integr Neurosci. 2022;16: 798995.
Sam S, Mazzone T. Adipose tissue changes in obesity and the impact on metabolic function. Transl Res. 2014;164(4):284–92.
Elgazar-Carmon V, et al. Neutrophils transiently infiltrate intra-abdominal fat early in the course of high-fat feeding. J Lipid Res. 2008;49(9):1894–903.
Watanabe Y, et al. Bidirectional crosstalk between neutrophils and adipocytes promotes adipose tissue inflammation. FASEB J. 2019;33(11):11821–35.
Reich T, et al. TREM2 has a significant, gender-specific, effect on human obesity. Sci Rep. 2023;13(1):482.
Jaitin DA, et al. Lipid-associated macrophages control metabolic homeostasis in a trem2-dependent manner. Cell. 2019;178(3):686-698.e14.
Kim K, et al. Characteristics of plaque lipid-associated macrophages and their possible roles in the pathogenesis of atherosclerosis. Curr Opin Lipidol. 2022;33(5):283–8.
Worthmann AH. TREM2-positive lipid-associated macrophages (LAMs) control white adipose tissue remodeling and metabolic adaptation in obesity. Immunometabolism. 2020;2(2):e200014.
Florance I, Ramasubbu S. current understanding on the role of lipids in macrophages and associated diseases. Int J Mol Sci. 2022;24(1):589.
Tsai VW, et al. Treatment with the TGF-b superfamily cytokine MIC-1/GDF15 reduces the adiposity and corrects the metabolic dysfunction of mice with diet-induced obesity. Int J Obes (Lond). 2018;42(3):561–71.
van der Heijden RA, et al. Obesity-induced chronic inflammation in high fat diet challenged C57BL/6J mice is associated with acceleration of age-dependent renal amyloidosis. Sci Rep. 2015;5:16474.
Kanda H, et al. MCP-1 contributes to macrophage infiltration into adipose tissue, insulin resistance, and hepatic steatosis in obesity. J Clin Invest. 2006;116(6):1494–505.
Dommel S, Bluher M. Does C-C motif chemokine ligand 2 (CCL2) link obesity to a pro-inflammatory state? Int J Mol Sci. 2021;22(3):1500.
McMillin M, et al. Neuronal CCL2 is upregulated during hepatic encephalopathy and contributes to microglia activation and neurological decline. J Neuroinflammation. 2014;11:121.
Cerri C, et al. The chemokine CCL2 mediates the seizure-enhancing effects of systemic inflammation. J Neurosci. 2016;36(13):3777–88.
Da Re C, et al. Neuroinflammation induced by lipopolysaccharide leads to memory impairment and alterations in hippocampal leptin signaling. Behav Brain Res. 2020;379: 112360.
Rummel C, et al. Leptin regulates leukocyte recruitment into the brain following systemic LPS-induced inflammation. Mol Psychiatry. 2010;15(5):523–34.
De Souza CT, et al. Consumption of a fat-rich diet activates a proinflammatory response and induces insulin resistance in the hypothalamus. Endocrinology. 2005;146(10):4192–9.
Posey KA, et al. Hypothalamic proinflammatory lipid accumulation, inflammation, and insulin resistance in rats fed a high-fat diet. Am J Physiol Endocrinol Metab. 2009;296(5):E1003–12.
Jeon BT, et al. Resveratrol attenuates obesity-associated peripheral and central inflammation and improves memory deficit in mice fed a high-fat diet. Diabetes. 2012;61(6):1444–54.
Milanski M, et al. Saturated fatty acids produce an inflammatory response predominantly through the activation of TLR4 signaling in hypothalamus: implications for the pathogenesis of obesity. J Neurosci. 2009;29(2):359–70.
Thaler JP, et al. Obesity is associated with hypothalamic injury in rodents and humans. J Clin Invest. 2012;122(1):153–62.
Puig KL, et al. Amyloid precursor protein and proinflammatory changes are regulated in brain and adipose tissue in a murine model of high fat diet-induced obesity. PLoS ONE. 2012;7(1): e30378.
Paolicelli RC, et al. Microglia states and nomenclature: a field at its crossroads. Neuron. 2022;110(21):3458–83.
Chintamen S, et al. Distinct microglial transcriptomic signatures within the hippocampus. PLoS ONE. 2024;19(1): e0296280.
Chen Y, Colonna M. Microglia in Alzheimer’s disease at single-cell level. Are there common patterns in humans and mice? J Exp Med. 2021;218(9):e20202717.
Krasemann S, et al. The TREM2-APOE Pathway Drives the Transcriptional Phenotype of Dysfunctional Microglia in Neurodegenerative Diseases. Immunity. 2017;47(3):566-581.e9.
McQuade A, et al. Gene expression and functional deficits underlie TREM2-knockout microglia responses in human models of Alzheimer’s disease. Nat Commun. 2020;11(1):5370.
Brandebura AN, et al. Astrocyte contribution to dysfunction, risk and progression in neurodegenerative disorders. Nat Rev Neurosci. 2023;24(1):23–39.
Pandey S, et al. Disease-associated oligodendrocyte responses across neurodegenerative diseases. Cell Rep. 2022;40(8): 111189.
Bachstetter AD, et al. Attenuation of traumatic brain injury-induced cognitive impairment in mice by targeting increased cytokine levels with a small molecule experimental therapeutic. J Neuroinflammation. 2015;12:69.
Vincent JC, et al. Correction : IL-1R1 signaling in TBI: assessing chronic impacts and neuroinflammatory dynamics in a mouse model of mild closed-head injury. J Neuroinflammation. 2023;20(1):287.
Acknowledgements
We are grateful for the technical support provided by Xiaoyi Lin and Cera Huffman.
Funding
This work was supported by National Institutes of Health (NIH) grants R01 NS110756 to A.I.F. and B.A.S. and R01 NS129094 to B.A.S. and R.J.H.; VA grants I01BX004866 and I01RX003379 to B.A.S; MPOWER seed grants to R.J.H. and B.A.S; Royal Society Science Foundation Ireland (RS-SFI) University Research Fellowship (URF) grant URF\R1\221653 to R.J.H
Author information
Authors and Affiliations
Contributions
R.J.H. and B.A.S designed research. R.J.H., and J.P.B carried out research. R.J.H., J.P.B., M.V., N.Z.K., O.M., R.M.R., and B.A.S., analyzed data. R.J.H. wrote the first draft of the manuscript, and B.A.S., and A.I.F., edited the paper to produce final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethics approval and consent to participate
All protocols involving the use of animals complied with the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health (NIH) (DHEW publication NIH 85- 23–2985) and were approved by the University of Maryland Animal Use Committee.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Gating strategy for CD45+/CD11b+ and CD11b+ Red Bead+ ATM/myeloid cell populations. A Gating strategy for CD45+/CD11b+ ATM/myeloid cell populations in the adipose tissue at 28 dpi. B Gating strategy for CD11b+ Red Bead+ ATM/myeloid cell population in the adipose tissue at 28 dpi. Although there was a factor effect of diet on % Red Bead+ of CD11b+ ATM/myeloid cells, there was no reported differences between groups in post-hoc analysis (C). Post-hoc analysis revealed that TBI-HFD had a significant increase in MFI of Red Bead+ of CD11b+ ATM/myeloid cells when compared to TBI-SD counterparts (D). Figure S2. Morpheus unsupervised hierarchical clustering heat map on VAT cellular transcriptome. Only genes that showed evidence of significant changes (p < 0.05 in at least one of the two-way ANOVA outcomes) were included (369 genes). The heat map demonstrated the same ordered clustering across the groups observed Fig. 9A. Figure S3. IPA detects Upstream Regulator-Genes upregulated and downregulated genes in VAT. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05; all pathways with z-score > 2.8 or < − 2.8 in at least one group are presented. Analysis of the Z-score of adipose upstream upregulated genes pathways including multiple inflammatory-related, identify significant TBI, HFD, and TBI-HFD interaction effects. Notably, TBI-HFD significantly increases upregulated pathways compared to TBI-SD. Only diet was a significant factor for downregulated genes. Figure S4. IPA detects activated and inhibited Upstream Regulator-Drugs pathways in VAT. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05; all pathways with z-score > 2.8 or < − 2.8 in at least one group are presented. Analysis of the Z-score of adipose upstream upregulated and downregulated drugs pathways, including multiple inflammatory-related, identify significant TBI, HFD, and TBI-HFD interaction effects. Notably, TBI-HFD significantly increases upregulated pathways and decreases downregulated pathways compared to TBI-SD. Figure S5. IPA detects upregulated and downregulated Diseases and Bio Functions pathways in VAT. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05; all pathways with z-score > 3.5 or < − 3.5 in at least one group are presented. Analysis of the Z-score of adipose disease pathways, including multiple inflammatory-related, identify significant TBI, HFD, and TBI-HFD interaction effects (except lack of significant injury effect in downregulated pathways). Notably, TBI-HFD significantly increases upregulated pathways and decreases downregulated pathways compared to TBI-SD. Figure S6. HFD selectively exacerbates TBI-induced increases in the expression of hippocampal, but not cortical neuroinflammatory genes. A. TBI-SD mice display a significant increase in hippocampal expression of IL-1β, NOX2, and p22phox when compared to Sham-SD counterparts, at 28 dpi. TBI-HFD mice exhibit a significant increase in hippocampal expression of p22phox when compared to Sham-HFD and TBI-SD. In addition, TBI-HFD displays a significant increase in hippocampal expression of CD11b, when compared to Sham-HFD counterparts. B. TBI-SD mice display an increase in cortical expression of TNF-α, NOX-2, and p22phox when compared to Sham-SD counterparts. *p < 0.05, **p < 0.01 vs Sham-SD; + + + p < 0.001 vs TBI-SD; ###p < 0.001 vs Sham-HFD. N = 6–7 per group. Data expressed as mean ± SEM. Figure S7. Morpheus unsupervised hierarchical clustering heat map on isolated microglia (CD11b+) cells. Although only the genes that showed evidence of significant changes (p < 0.05 in at least one of the two-way ANOVA outcomes) were included (as for Fig. 8) the pattern remained to Fig. 5A. Figure S8. IPA detects activated and inhibited Upstream Regulator-Genes pathways in microglia. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05; all pathways with z-score > 2.5 or < -2.5 in at least one group are presented. TBI, but not HFD (except for a modest significant factor effect in upregulated) had strong significant factor effects in both upregulated (including several inflammatory-related pathways) and downregulated pathways; no significant TBI-HFD interaction was observed. Figure S9. IPA detects activated and inhibited Upstream Regulator-Drugs pathways in microglia. IPA determined the activation of pathways based on the inclusion of genes showing changes in each group vs Sham-SD with p < 0.05; all pathways with z-score > 2.5 or <− 2.5 in at least one group are presented. There was a significant effect of TBI, HFD, and TBI-HFD interaction in upregulated upstream drug pathways which included several inflammatory pathways. Notably, TBI-HFD showed increased effects compared to TBI-SD. However, only TBI effects were reported on downregulated upstream drugs in microglia. Figure S10. IPA detects upregulated and downregulated Diseases and Bio Functions pathways in microglia. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05; all pathways with z-score > 2.5 or < − 2.5 in at least one group are presented. TBI, but not HFD, or TBI-HFD effects are reported on upregulated and downregulated diseases in microglia. Figure S11. Morpheus unsupervised hierarchical clustering heat map on microglia vs. flowthrough cellular transcriptome. The ratio between global microglia and flowthrough normalized expression data (560 genes out of a total of 769 genes had average expression > 50) was used to determine the top 50 overexpressed and under-expressed genes in microglia vs flowthrough. Figure S12. Morpheus unsupervised hierarchical clustering heat map on flowthrough cellular transcriptome. Morpheus heat map on all 560 genes shows no ordered clustering across the groups indicative of the smaller number of genes with consistent and significant changes. Figure S13. Morpheus unsupervised hierarchical clustering heat map on flowthrough cellular transcriptome. When only genes showing significant changes (p < 0.05 in at least one of the two-way ANOVA outcomes) were selected the heat map demonstrated the same ordered clustering across the groups observed in microglia. Genes downregulated included neurogenesis and neuronal signaling markers while genes upregulated included oligodendrocyte phenotype markers. Figure S14. IPA detects Upstream Regulator-Genes upregulated in flowthrough cellular transcriptome. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05. TBI, but no HFD or TBI-HFD showed significant effects on upstream genes upregulated in flowthrough cellular transcriptome. Figure S15. IPA detects upregulated and downregulated Diseases and Bio Functions pathways in flowthrough cellular transcriptome. IPA determined the activation of pathways based on the inclusion of genes with p < 0.05; all pathways with z-score > 2.0 or < − 2.0 in at least one group are presented. Significant TBI, but no HFD or TBI-HFD interaction effects are observed in diseases upregulated and downregulated pathways in flowthrough cellular transcriptome. Notably, TBI-HFD showed increased effects compared to TBI-SD in upregulated pathways. Figure S16. Statistical analysis tables for the main figures. We present the statistical tables for all the transcriptomic (Nanostring) and pathways analysis (IPA) data included in the main figures. The corresponding figure numbers are indicated.
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
Henry, R.J., Barrett, J.P., Vaida, M. et al. Interaction of high-fat diet and brain trauma alters adipose tissue macrophages and brain microglia associated with exacerbated cognitive dysfunction. J Neuroinflammation 21, 113 (2024). https://doi.org/10.1186/s12974-024-03107-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12974-024-03107-6