- Department of Physiology and Pharmacology, College of Veterinary Medicine, University of Georgia, Athens, GA, United States
Gulf War Illness (GWI) affects nearly 30% of veterans from the 1990–1991 Gulf War (GW) and is a multi-symptom illness with many neurological effects attributed to in-theater wartime chemical overexposures. Brain-focused studies have revealed persistent structural and functional alterations in veterans with GWI, including reduced volumes, connectivity, and signaling that correlate with poor cognitive and motor performance. GWI symptomology components have been recapitulated in rodent models as behavioral, neurochemical, and neuroinflammatory aberrations. However, preclinical structural imaging studies remain limited. This study aimed to characterize the progression of brain structural alterations over the course of 12 months in two established preclinical models of GWI. In the PB/PM model, male C57BL/6 J mice (8–9 weeks) received daily exposure to the nerve agent prophylactic pyridostigmine bromide (PB) and the pyrethroid insecticide permethrin (PM) for 10 days. In the PB/DEET/CORT/DFP model, mice received daily exposure to PB and the insect repellent DEET (days 1–14) and corticosterone (CORT; days 7–14). On day 15, mice received a single injection of the sarin surrogate diisopropylfluorophosphate (DFP). Using a Varian 7 T Bore MRI System, structural (sagittal T2-weighted) scans were performed at 6-, 9-, and 12-months post GWI exposures. Regions of interest, including total brain, ventricles, cortex, hippocampus, cerebellum, and brainstem were delineated in the open source Aedes Toolbox in MATLAB, followed by brain volumetric and cortical thickness analyses in ImageJ. Limited behavioral testing 1 month after the last MRI was also performed. The results of this study compare similarities and distinctions between these exposure paradigms and aid in the understanding of GWI pathogenesis. Major similarities among the models include relative ventricular enlargement and reductions in hippocampal volumes with age. Key differences in the PB/DEET/CORT/DFP model included reduced brainstem volumes and an early and persistent loss of total brain volume, while the PB/PM model produced reductions in cortical thickness with age. Behaviorally, at 13 months, motor function was largely preserved in both models. However, the GWI mice in the PB/DEET/CORT/DFP model exhibited an elevation in anxiety-like behavior.
1 Introduction
Gulf War Illness (GWI), a persisting multi-symptom malady predominantly affecting the nervous and immune systems, presents in about a third of the veterans from the 1990–1991 Gulf War (White et al., 2016). Major symptoms include fatigue, aberrations in cognition/memory, reduced motor function, and mood disturbances. Precise etiology of GWI is still unknown but it is largely attributed to war-time chemical overexposures including to insecticides/repellents [permethrin (PM), chlorpyrifos, DEET], nerve agents (sarin/cyclosarin), and a nerve agent prophylactic (pyridostigmine bromide; PB). While the understanding of GWI has improved, symptoms remain heterogenous among affected veterans and curative treatments do not exist; however, therapeutic interventions are being explored (White et al., 2016).
Several clinical imaging studies have evaluated the neuroanatomical correlates to the progressive, long-term aberrations in GWI symptomology. These structural magnetic resonance imaging (MRI) studies found specific reductions in white and gray matter volumes in multiple cortical and subcortical areas (e.g., hippocampus, hypothalamus, brainstem, cerebellum) in GWI veterans (Chao et al., 2011; Rayhan et al., 2013; Chao et al., 2015; Christova et al., 2017; James et al., 2017; Van Riper et al., 2017; Chao, 2020; Zhang et al., 2020; Zhang et al., 2021). Functional MRI (fMRI) studies corroborated these structural abnormalities by detecting significant signal changes (e.g., activity) in regions governing cognitive, memory, motor and mood function such as the prefrontal, somatosensory, and motor cortices (Gopinath et al., 2012; Wylie et al., 2019), as well as the thalamus, caudate, and hippocampus (Calley et al., 2010; Li et al., 2011; Cooper et al., 2016). Across these studies, the patient populations included veterans diagnosed with GWI based on the Haley criteria which groups GWI symptoms into three syndromes: (1) impaired cognition (i.e., problems attention, memory, reasoning, insomnia, depression, daytime sleepiness, and headaches); (2) confusion-ataxia (i.e., problems thinking, disorientation, balance, vertigo, and impotence); and (3) arthromyoneuropathy (i.e., joint and muscle pain, muscle fatigue, difficulty lifting, and extremity paresthesia) (Haley et al., 1997; IOM, 2014). These studies found marked behavior aberrations that correlated with functional abnormalities (Calley et al., 2010; Li et al., 2011; Gopinath et al., 2012; Cooper et al., 2016; Wylie et al., 2019). Additionally, an earlier clinical study detected neuronal loss in the basal ganglia with corresponding dopamine dysfunction (Haley et al., 2000). The structural and functional abnormalities observed in these areas are consistent with reported symptomology by veterans and correlate strongly with performance deficits in working memory, attention, and motor tasks as well as with mood impairments, pain, and sleep quality (Calley et al., 2010; Hubbard et al., 2014; Clarke et al., 2019; Chao, 2020; Washington et al., 2020; Zhang et al., 2021). However, neurological abnormalities have not been correlated with types and/or severity of GWI exposures.
Many of the symptoms exhibited by GWI veterans have been recapitulated in GWI rodent models. Studies utilizing these models have demonstrated multiple behavioral deficits (i.e., cognition/memory, motor, and mood) and biological parameters (i.e., neuroinflammation) consistent with symptomology in GWI veterans (Parihar et al., 2013; Megahed et al., 2014; O’Callaghan et al., 2015; Zakirova et al., 2015; Zakirova et al., 2016; Miller et al., 2018; Macht et al., 2019; Dickey et al., 2021; Ribeiro and Deshpande, 2021). To date, our group has explored the acute and long-term neurobiological and behavioral alterations (Carpenter et al., 2020; Brown et al., 2021a; Brown et al., 2021b; Carpenter et al., 2021; Carpenter et al., 2022) in two established GWI models, one utilizing subacute PB and PM exposure (Zakirova et al., 2015; Zakirova et al., 2016) and another employing PB, DEET, corticosterone (CORT), and diisopropylfluorophosphate (DFP) exposure (O’Callaghan et al., 2015) In these studies, we observed: (1) acute monoaminergic disbalance in multiple brain regions associated with motor, memory, and mood function (Carpenter et al., 2020), (2) hippocampal neuroinflammation (Carpenter et al., 2020; Carpenter et al., 2021; Carpenter et al., 2022) and (3) changes in hippocampal synaptic plasticity and transmission (Brown et al., 2021a; Brown et al., 2021b). Further, our long-term studies (6–10 months post GWI exposures) revealed multiple GWI-related behavioral deficits in motor function, cognition, and mood (Carpenter et al., 2021; Carpenter et al., 2022) and alterations in hippocampal electrophysiology (Brown et al., 2021a; Carpenter et al., 2021).
While rodent GWI models produce GWI-like behavioral and pathological phenotypes, functional and structural imaging studies investigating the neuroanatomical changes within GWI models remain limited (Koo et al., 2018; Wu et al., 2021). This may be due to the complexity of conducting rodent imaging studies, especially studies that involve repeated imaging over time. Nevertheless, results from the few available GWI preclinical studies provide highly translational and valuable insights into this chronic illness (Koo et al., 2018; Wu et al., 2021). In particular, high order diffusion MRI revealed distinct cortical and subcortical (hippocampus and hypothalamus) alterations that corresponded to neuroinflammation following exposure to the GWI-relevant chemicals, DFP and CORT (Koo et al., 2018). Additionally, structurally significant increases in lateral ventricle volume and decreases in hippocampal and thalamic volumes were observed in another GWI model (PB, PM, DEET, and restraint stress) 10 months post exposure that aligned with neurobehavioral impairments in cognition and mood function (Wu et al., 2021).
Due to the paucity of longitudinal information surrounding GWI preclinical neuroanatomical changes and possible GWI exposure specific differences, the present study aimed to characterize the progression of neuroanatomical alterations over the course of 12 months following two, distinct GWI chemical exposure paradigms. Utilizing the PB/PM (Zakirova et al., 2015) and PB/DEET/CORT/DFP (O’Callaghan et al., 2015) models, we evaluated GWI neuroanatomical volume changes at 6-, 9-, and 12-months post exposure in various brain regions including the total brain, ventricles, cortex, hippocampus, cerebellum, and brainstem. Cortical thickness measurements were also conducted over the course of the study, alongside limited behavioral and stress responsivity evaluations at 13 months.
2 Materials and methods
2.1 Materials
The following chemicals were used for animal treatments: pyridostigmine bromide (PB; Sigma Aldrich, St. Louis, MO), permethrin (PM; 29.5% cis/69.5% trans isomer; Chem Service Inc., West Chester, PA), diisopropylfluorophosphate (DFP; Sigma Aldrich), N-Diethyl-3-methylbenzamide (DEET; Sigma Aldrich), and corticosterone (CORT; Steraloids, Newport, RI). All additional chemicals and reagents used in this study, unless otherwise noted, were of analytical or higher grade and were obtained from Sigma Aldrich or Fisher Scientific (Hampton, NH).
2.2 Animals and GWI models
2.2.1 Animals
Male C57BL/6 J mice (8–9 weeks old; Jackson Laboratories, Bar Harbor, ME) were housed 4 per cage in an environmentally controlled room (22–24°C) and maintained on a 12 h light/dark cycle (0700–1900 lights on) throughout the study. Mice were handled daily for 1 week prior to the start of the study to minimize experimenter-induced stress. Food and water were available ad libitum. All procedures were approved in advance by the University of Georgia Institutional Animal Care and Use Committee (IACUC), as well as by DoD ACURO, and were in accordance with the latest National Institutes of Health guidelines.
Body weights were measured daily during the GWI chemicals exposure and then biweekly following the last GWI exposure until study completion. All mice were euthanized 13.5 months post GWI chemicals exposure (CO2, followed by decapitation). When the mice were euthanized, body weight and length were measured to determine morphometric-related parameters including body mass index [BMI; Body weight (g)/length2 (cm2)] and the Lee index [body weight (g)/length(cm)] as in Novelli et al. (2007). Brains were extracted, weighed, and split sagittally before one half was quickly frozen on dry ice and the other half was immersion fixed in 4% paraformaldehyde as in Carpenter et al. (2021) for future analyses.
2.2.2 GWI models
Two established, chemically different GWI treatment paradigms were utilized for this study (Figure 1). For both models, mice were selected for treatment randomly. Following the Zakirova et al. model (Zakirova et al., 2015), mice were co-administered the nerve agent prophylactic pyridostigmine bromide [PB; 0.7 mg/kg body weight (BW); IP] and the pyrethroid insecticide permethrin (PM; 200 mg/kg BW; IP) or DMSO vehicle over 10 days. Following the O’Callaghan et al. model (O’Callaghan et al., 2015), mice received daily administration of PB (2 mg/kg BW; SC) and the insect repellent DEET (30 mg/kg BW; SC) or saline control for 14 days with concurrent stress [corticosterone (CORT): 200 mg/kg in 1.2% EtOH drinking water] exposure on days 8–14. Control mice during this period (day 8–14) received 1.2% EtOH drinking water. On day 15, mice received a single dose of the nerve agent surrogate diisopropylfluorophosphate (DFP; 3.75 mg/kg BW; IP) or saline control. This resulted in 4 treatment groups (N = 24, n = 6/group/per model).
Figure 1. Experimental timeline and MRI collection. Experimental timeline for the study investigating two GWI exposure paradigms (PB/PM: 10 days; PB/DEET/CORT/DFP: 15 days). MRI evaluations were conducted at 6-, 9-, and 12-months post-exposure. Open field and temperature response stress tests were performed at 13 months post-exposure. Figure created with BioRender.com. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; GWI, Gulf War Illness; MRI, magnetic resonance imaging; PB, pyridostigmine bromide; PM, permethrin.
2.3 Structural magnetic resonance imaging
Structural brain images were acquired using a 7 Tesla magnet (Agilent, Santa Clara, CA), and MRI sequences were conducted on mice at 6, 9, and 12 months post GWI chemicals exposures (Figure 1). Within a model, order of imaging was random. For all imaging processing post collection, individual animals were coded, and MRI data were processed and analyzed by a treatment-blinded experimenter. Prior to imaging, mice were anesthetized with isoflurane (3% for induction, 1.0–1.5% for maintenance) in a 30%:70% O2:N2 gas mixture with a flow rate 0.8–1.0 L/min. Respiratory rate was monitored using a small animal monitoring system (Small Animal Instruments, Inc., Stony Brook, NY) throughout the imaging period. Axial, two dimensional (2D) T1 and T2 weighted images were obtained using a spin echo sequence with the following acquisition parameters: TR 500 ms, TE 17 ms, 8 averages, FOV 25, data matrix 256 × 256, 17 slices, thickness 1.00 mm with no gap (Krishna et al., 2014). Sagittal, 3D T2 weighted images were obtained using a fast spin echo sequence with the following parameters: TR 4,000 ms, TE 33.38 ms, 2 averages, FOV 35, data matrix 256 × 256, 15 slices, thickness 1.00 mm with no gap.
2.3.1 Volume measurements
Images were displayed and masked using the open-source MATLAB toolbox AEDES using methods outlined in (Minkeviciene et al., 2019; Grant et al., 2020). Regions of interest (ROIs) including total brain, cortex, hippocampus, cerebellum, medulla, pons, and ventricles (i.e., lateral, 3rd, and 4th) were manually drawn according to the Allen Brain Atlas.1 ROIs were saved individually as masks, and volumetric analysis was conducted on the masks by using the pixel count generated with the AEDES toolbox. Volume (mm3) was calculated by first calculating the voxel volume and then multiplying the measured output with the calculated voxel volume. Refer to Figures 2A–E for volume analysis process and representative masks.
Figure 2. MRI image processing for volume and cortical thickness analyses. (A) For volume analysis of brain regions, sagittal T2 images were processed and brain regions of interest (ROIs) were delineated manually using the MATLAB Aedes toolbox. Representative ROIs depicted include the total brain (B,D), cortex (C, yellow; E, green), hippocampus (C, purple; E, orange), cerebellum (C, teal; E, purple), medulla (C, orange; E, dark blue), pons (C, green; E, yellow), lateral ventricle (C, light green), 3rd ventricle (E, pink), and 4th ventricle (C, red; E, light blue). (F) For the cortical thickness analysis, cortical ROIs were delineated using the straight tool in ImageJ on T2-weighted axial images, as represented in (G) showing in the somatosensory cortex.
2.3.2 Thickness measurements
Axial images were exported using the open-source MATLAB toolbox AEDES. Thickness measurements were done using the straight tool in ImageJ by manually drawing the thickness of ROIs including the frontal, motor, auditory, somatosensory, and dorsolateral entorhinal cortices as in Minkeviciene et al. (2019). Six measurement lines were drawn per slice, and line thickness in pixels (e.g., the number of pixels for cortical width) was obtained using the measure function in ImageJ. Cortical thickness was calculated by converting the number of pixels into mm. Refer to Figures 2F,G for the cortical thickness analysis process and a representative image.
2.4 Behavioral analyses
Limited behavioral assays were conducted at 13 months to evaluate sensitivity to a mild stressor, locomotor activity, and anxiety-like activity. The experimenter was treatment-blinded during the behavioral tests and behavioral sequence was randomized.
2.4.1 Open field
Locomotor activity and anxiety-like behavior were assessed as previously described (Carpenter et al., 2021). Briefly, mice were individually placed into an open field arena (25 cm × 25 cm × 40 cm; Coulbourn Instruments, Whitehall, PA) and allowed to freely explore for 30 min. Locomotor parameters (i.e., distance traveled) and anxiety-like behaviors (i.e., entries and time spent in the center and corners) were scored using AnyMaze software (Stoelting) for the total 30 min and per 5 min intervals.
2.4.2 Temperature response stress test
The temperature response stress test was used to assess stress-induced hyperthermia following methodology outlined by Borsini et al. (1989) and Bouwknecht et al. (2000). Rectal temperature measurements were obtained utilizing a thermocouple meter (Digi-Sense Single-Input thermometer) with a mouse-specific probe. Two temperature readings were recorded with a 15-min interval between each reading. After the first reading, the mouse was individually housed until the second reading. The change in temperature was calculated and analyzed.
2.5 Statistical analysis
All data were analyzed first by using preplanned comparisons within a model and a time point using Student’s t-test. To determine the effects of treatment over time or treatment x time interactions, a two-way repeated analysis of variance (RM-ANOVA) within a GWI model was used. Here, if a RM-ANOVA was significant (p ≤ 0.05), treatment means were separated by Student–Newman–Keuls (SNK) post hoc test. All data were analyzed using SigmaPlot 12.5 (San Jose, CA), and graphs were generated using GraphPad Prism (San Diego, CA).
3 Results
3.1 Longitudinal weights
Body weights were not different in either model at the beginning of the study (PB/PM model: control, 25.7 ± 0.76 vs. PB/PM, 25.9 ± 0.62, p ≥ 0.89; PB/DEET/CORT/DFP model: control 26.1 ± 0.77 vs. PB/DEET/CORT/DFP, 25.1 ± 0.70 p ≥ 0.38). In both models, weights significantly increased over time [PB/PM model: F(14, 140) = 162.24, p ≤ 0.001, Figure 3A; PB/DEET/CORT/DFP model: F(14, 137) = 76.20, p ≤ 0.001, Figure 3C]. In the PB/PM model, body weights were not significantly different between treatments throughout the study except for a numerical increase by the control group at 12 months (p = 0.10) (Figure 3A). An opposite numerical increase was observed in the PB/DEET/CORT/DFP treated mice (p = 0.13) (Figure 3C). Further, in the PB/PM model, there were no treatment differences for the calculated BMI or Lee indices (p’s ≥ 0.33) (Figure 3B). The Lee index, albeit numerically greater, was unaffected by treatment in the PB/DEET/CORT/DFP model (p ≥ 0.29). However, there was a strong trend for increased BMI in the PB/DEET/CORT/DFP treated mice [t(9) = 2.06, p = 0.07] (Figure 3D). At euthanasia, brain weights, both absolute and normalized to body weight, did not differ in either model (p’s ≥ 0.24, data not shown).
Figure 3. Monthly weights and calculated morphometrical parameters. (A,C) Weights were monitored and are presented for the start of the study (pre-treatment), directly after treatment administration cessation (PB/PM: 10 days; PB/DEET/CORT/DFP: 15 days), and monthly until the end of the study at 13.5 months. (B,D) Morphometrical parameters were evaluated by calculating the body mass index [BMI; Body weight (g)/length2 (cm2)] and the Lee index [body weight (g)/length (cm)]. Data are presented at ± SEM. Sample sizes were n = 6/group/timepoint for the PB/PM model and n = 5–6/group/timepoint for the PB/DEET/CORT/DFP model. ***Indicates a significant effect of time, p ≤ 0.001. ^Indicates a trending effect of treatment, p ≤ 0.10. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; g, grams; PB, pyridostigmine bromide; PM, permethrin.
3.2 Prior PB/DEET/CORT/DFP exposure, but not PB/PM exposure, leads to overall brain volume reduction
In the PB/PM model, there were no significant effects of time or treatment on overall brain volume (p’s ≥ 0.42) (Figure 4A). Over the course of the study, a slight 0.42% increase and 0.21% decrease in total brain volume were observed in the GWI and control groups, respectively (Supplementary Table S1). In the PB/DEET/CORT/DFP model, significant changes in global brain volume were observed for treatment [F(1, 19) = 8.55, p ≤ 0.05] (Figure 4B) and time [F(2, 19) = 7.22, p ≤ 0.01] (Figure 4B). Here, there were significant reductions in total brain volume for the PB/DEET/CORT/DFP group compared to control at the 6, 9, and 12-month time points (p’s ≤ 0.05), and this translated to an overall reduction of 0.08 and 0.96% in brain volume over the course of the study in the control and GWI groups, respectively (Supplementary Table S1). Within treatment, total brain volume was reduced significantly from 9 to 12 months in the PB/DEET/CORT/DFP (p ≤ 0.05) (Figure 4B) group; in the control mice, there was a trend for the age-related decrease of brain volume at 12 months (vs. 9 months, p = 0.07) (Figure 4B). This suggests GWI-related early loss of total brain volume, especially in the PB/DEET/CORT/DFP model, that is accelerated by age.
Figure 4. Longitudinal total brain volumes in two models of Gulf War Illness. Longitudinal measurement of total brain volume in (A) PB/PM and (B) PB/DEET/CORT/DFP models of GWI. Data are presented as mean ± SEM. Sample sizes were: PB/PM model, n = 6/group/timepoint; PB/DEET/CORT/DFP model, n = 5–6/group/timepoint. **Indicates significant main effect of time where p ≤ 0.01. †Indicates significant main effect of treatment, p ≤ 0.05, and ‡ denotes a significant pairwise comparison of treatment within a time point. a and a^ denote a significant (p ≤ 0.05) and trending (p ≤ 0.10) pairwise comparison from the 9-month time point, respectively. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin.
3.3 Ventricular enlargement is present in both models over time
Total ventricular size (3rd, 4th, and lateral ventricles) was increased in the PB/PM [F(2, 30) = 10.87, p ≤ 0.001] (Figure 5A) model over time. An early increase was observed in the GWI group starting at 9 months (vs. 6 months, p ≤ 0.05), and ventricles became larger at 12 months (vs. 6 months, p’s ≤ 0.01). Additionally, an increase in ventricular size was observed in the control group at 12 months (vs. 6 months, p ≤ 0.05; vs. 9 months p = 0.08). By study end, this resulted in a 17 and 25% increase in volume from the 6-month timepoint in the control and GWI groups, respectively (Supplementary Table S1). Ventricular enlargement was also observed over time in the PB/DEET/CORT/DFP model [F(2, 19) = 7.41, p ≤ 0.01] (Figure 5B) in addition to a trending effect of treatment [F(1, 19) = 3.73, p = 0.08] (Figure 5B). Here, this treatment trend was driven by initially smaller ventricular volumes in the GWI group within the 6-month time point [t(10) = −3.25, p ≤ 0.01] (Figure 5B), an effect likely resulting from the decreased total brain volume in this model (Figure 4B). However, with age, significant increases in total ventricular size were seen in GWI mice at the 9-and 12-month time points (vs. 6, p’s ≤ 0.05). At the end of the study, the ventricular volume in the GWI group increased 20% from the 6-month time point, whereas the control increased 11% (Supplementary Table S1).
Figure 5. Longitudinal MRI analysis of ventricular volume in two models of Gulf War Illness. Longitudinal measurement of total ventricular volume (mm3) in (A) PB/PM and (B) PB/DEET/CORT/DFP models of GWI. Data are presented as mean ± SEM. Sample sizes were: PB/PM model, n = 6/group/timepoint; PB/DEET/CORT/DFP model, n = 5–6/group/timepoint. ** and *** indicate a significant main effect of time, p ≤ 0.01 and 0.001, respectively. ^Indicates a trending main effect of treatment, p = 0.08. a and b denote a significant pairwise comparison from the 6- and 9-month time points (p ≤ 0.05), respectively, whereas b^ denotes a trending pairwise comparison from the 9-month time point (p = 0.08). ‡Denotes a significant pairwise comparison of treatment within a time point. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin.
Individual ventricles were delineated to determine where these observed increases originated. In the PB/PM model, there was a significant volume increase in both the 4th ventricle [F(2, 30) = 4.25, p ≤ 0.05] and lateral ventricle [F(2, 30) = 12.93, p ≤ 0.001] over time, but not the 3rd ventricle (Table 1). Within the 4th ventricle, there was a significant volume increase in both groups at 9 months that was driven by the GWI group (vs. 6 months, p ≤ 0.05) (Table 1). Interestingly, this increase at 9 months was transient, as it did not persist to month 12 (p = 0.16) (Table 1). Lateral ventricle size gradually increased over time in both groups (p’s ≤ 0.05), and this was driven by early, significant increases in the GWI group at 9 months that persisted to 12 months (vs. 6 months, p’s ≤ 0.01) (Table 1) and a significant increase in the control group at 12 months (vs. 6 and 9 months, p’s ≤ 0.05) (Table 1).
Table 1. Longitudinal MRI assessment of ventricular and hindbrain area volumes in the PB/PM model of Gulf War Illness.
Within the PB/DEET/CORT/DFP model, the total ventricular enlargement was largely driven by significant treatment-related increases in the lateral ventricle over time {Treatment: [F(1, 19) = 7.62, p ≤ 0.05] and Time: [F(2, 19) = 15.33, p ≤ 0.001]; Table 2}. Lateral ventricle size gradually increased over time in both groups (control, 6 vs. 9 and 12 months, p’s ≤ 0.05), however, these increases were more pronounced within the GWI group (6 vs. 9 and 12 months, p’s ≤ 0.01 and 0.001, respectively). Interestingly, likely due to the decreased brain volume overall, lateral ventricle size within the GWI group was significantly smaller than the control group at all-time points (6 and 9 months, p ≤ 0.05, 12 months, p = 0.06). While there were numerical increases in the 3rd and 4th ventricles over time, none were significant (p’s ≥ 0.16).
Table 2. Longitudinal MRI assessment of the ventricular and hindbrain area volumes in the PB/DEET/CORT/DFP model of Gulf War Illness.
3.4 Cortical volume and thickness are decreased in both models over time
There were no significant time or GWI treatment effects for cortical volume observed in the PB/PM model (p’s ≥ 0.36) (Figure 6A); over time, cortical volume decreased by 0.09 and 0.18% in the control and GWI groups, respectively (Supplementary Table S1). Similarly, in the PB/DEET/CORT/DFP model, the age effect on cortical volume was not significant (p’s ≥ 0.64). The percent change from 6 to 12 months was 0.56 and 2% for the control and GWI groups, respectively (Supplementary Table S1). Notably, a significant treatment effect was observed at 9 months where there was a reduction in cortical volume in the GWI mice compared to control [t(10) = −2.20, p ≤ 0.05] (Figure 6B); a similar, non-significant trend was observed at 12 months (Figure 6B).
Figure 6. Longitudinal MRI analysis of cortical volume and thickness in two models of Gulf War Illness. Longitudinal measurement of cortical volumes (mm3) in (A) PB/PM and (B) PB/DEET/CORT/DFP models of GWI. Longitudinal measurement of thickness (in mm) of the somatosensory and dorsolateral entorhinal cortices were conducted in (C,E) PB/PM and (D,F) PB/DEET/CORT/DFP models of GWI, respectively, over the course of 12 months. Data are presented as mean ± SEM. Sample sizes were: PB/PM model, n = 6/group/timepoint; PB/DEET/CORT/DFP model, n = 5–6/group/timepoint. **Indicates a significant main effect of time where p ≤ 0.01. †Indicates significant main effect of treatment, p ≤ 0.05. ‡Denotes a significant pairwise comparison of treatment within a time point. a and b denote a significant pairwise comparison from the 6- and 9-month time points, respectively, p ≤ 0.05. a^Denotes a trending effect from the 6-month time point, p ≤ 0.10. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin.
To complement the volume analyses, cortical thickness measurements were conducted across the frontal, motor, auditory, somatosensory, and dorsolateral entorhinal cortices in both models. In the PB/PM model, there were no significant effects of treatment observed in the somatosensory cortex (Figure 6C). In the frontal cortex, there was a transient decrease in thickness within the GWI group at 9 months (vs. 6 months, [t(10) = 2.68, p ≤ 0.05]; Table 3) that did not persist to 12 months. A trending decrease in motor cortex thickness within the control group at 12 months {vs. 6 months, [t(10) = 2.08, p = 0.06; Table 3]} also resulted in a treatment trend when controls were compared with the GWI group [t(10) = −2.00, p = 0.07] (Table 3). Similar reductions were apparent within the auditory cortex over time [F(2, 30) = 4.62, p ≤ 0.05]; these effects were driven largely by numerical decreases in both groups at month 12 (vs 6 months, p’s ≤ 0.10) (Table 3). Significant thinning of the dorsolateral entorhinal cortex was apparent with age [F(2, 30) = 7.78, p ≤ 0.01] (Figure 6E) and treatment [F(2, 30) = 6.11, p ≤ 0.05] (Figure 6E). Interestingly, there was an increase in thickness at 9 months in the GWI group (vs. 6 months, p ≤ 0.05) that was significantly higher than the control group (p ≤ 0.01). This increase was transient, as thickness decreased by 12 months (vs. 9 months, p ≤ 0.01 and vs. 6 months, p = 0.07; Figure 6E) and was similar to the control group.
In the PB/DEET/CORT/DFP model, there were no significant changes in cortical thickness observed with treatment or over time in the frontal, motor, auditory, or dorsolateral entorhinal cortices (Table 3; Figure 6F). However, there was a trending decrease in thickness over the course of 12 months within the somatosensory cortex [F(2, 19) = 2.95, p = 0.08] (Figure 6D); this numerical decrease was driven by significant decreases in thickness within the GWI group for this region at 12 months {vs. 6 months: [t(10) = −2.96, p ≤ 0.05] and vs. 9 months: [t(10) = −2.48, p ≤ 0.05]; Figure 6D).
3.5 Reductions in hippocampal volumes are present in both models over time
Reduced hippocampal volume was present over time in the PB/PM [F(2, 30) = 7.39, p ≤ 0.01] (Figure 7A) and PB/DEET/CORT/DFP [F(2, 19) = 5.24, p ≤ 0.05] (Figure 7B) models. Pairwise comparisons revealed that hippocampal volume in the PB/PM model began to decline in the PB/PM group from month 6 to 9 (p ≤ 0.05) (Figure 7A) and continued to month 12 (p ≤ 0.01) (Figure 7A). In the control group, this effect was only observed at 12 months (9 vs. 12 months, p ≤ 0.05) (Figure 7A), highlighting an earlier onset of the decrease driven by GWI treatment. In the PB/DEET/CORT/DFP model, both groups had reduced hippocampal volumes at 12 months (9 vs. 12 months, p’s ≤ 0.05) (Figure 7B). Within the PB/PM model, the overall percent change in volume from 6 to 12 months was 7.75 and 11.40% in the control and GWI groups, respectively (Supplementary Table S1). Within the PB/DEET/CORT/DFP model, the percentage change in volume from 6 to 12 months was 5.03 and 2.03% in the control and GWI groups, respectively (Supplementary Table S1).
Figure 7. Longitudinal MRI analysis of hippocampal volume in two models of Gulf War Illness. Longitudinal measurement of hippocampal volumes in (A) PB/PM and (B) PB/DEET/CORT/DFP models of GWI. Data are presented as mean ± SEM. Sample sizes were: PB/PM model, n = 6/group/timepoint; PB/DEET/CORT/DFP model, n = 5–6/group/timepoint. * and ** indicate a significant main effect of time where p ≤ 0.05 and 0.01, respectively. a and b denote a significant pairwise comparison from the 6- and 9-month time point, respectively, p ≤ 0.05. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin.
3.6 Alterations in the hindbrain over time in both models
Overall, no significant effects for time or treatment were observed in the hindbrain (cerebellum, medulla, and pons) in the PB/PM model (p’s ≥ 0.22; Figure 8A). However, the percentage change in hindbrain from 6 to 12 months indicated an overall reduction of 3.95% in the GWI group and an increase of 0.21% in the control group. A numerical (p = 0.12) reduction in hindbrain volume was observed over time within the PB/DEET/CORT/DFP model; this translated to an overall reduction of 2.15 and 1.27% in the control and GWI groups, respectively from 6 to 12 months. Notably, the hindbrain volume of the GWI group was significantly smaller than the control group [F(1, 19) = 8.46, p ≤ 0.05] (Figure 8B), particularly within 6- and 9-month time points (p’s ≤ 0.05) (Figure 8B), which might partly explain the total brain volume difference between groups in this model.
Figure 8. Longitudinal MRI analysis of the hindbrain volume in two models of Gulf War Illness. Longitudinal measurement of the hindbrain (cerebellum, pons, and medulla) volumes in the (A) PB/PM and (B) PB/DEET/CORT/DFP models of GWI. Data are presented as mean ± SEM. Sample sizes were n = 6/group/timepoint for the PB/PM model and n = 5–6/group/timepoint for the PB/DEET/CORT/DFP model. †Denotes a significant main effect of treatment (p ≤ 0.05), and ‡ indicates a significant pairwise comparison of treatment within a time point (p ≤ 0.05). CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin.
Similar to the ventricular volumetric analysis, the individual areas of the hindbrain (cerebellum, medulla, and pons) were analyzed to gather more specific information on these brain regions over time. Interestingly, cerebellar volumes fluctuated slightly in the PB/PM model over the course of 12 months (Table 1) and were significantly decreased over time in the PB/DEET/CORT/DFP model [F(2, 19) = 5.37, p ≤ 0.05] (Table 2). A significant reduction in cerebellar volume was observed at 12 months within the control group of the PB/DEET/CORT/DFP model (vs. 6 months, p ≤ 0.05); no significant reductions were observed in the GWI group with age (Table 2). Notably, a significant main effect of treatment was present within the cerebellum [F(1, 19) = 6.30, p ≤ 0.05] where cerebellar volume was smaller within the GWI group compared to control, particularly at months 6 and to a lesser extent, 9 (p’s ≤ 0.05 and 0.08, respectively).
Within the brainstem, both medullary and pontine volumes numerically decreased over the course of 12 months in the PB/PM model (Table 1). In the pons, there was a trending effect for reduction in volume over time [F(2, 30) = 2.70, p = 0.08]. However, this reduction was not significant in the medulla (p = 0.66). Within the PB/DEET/DFP/CORT model, there were significant and trending main effects of treatment for the medulla [F(1, 19) = 8.02, p ≤ 0.05] and pons [F(1, 19) = 3.40, p = 0.09], respectively (Table 2). Here, medullary volume was significantly smaller in the GWI group than controls at the 6-and 9-month time points (p ≤ 0.05) and trended at the 12-month time point (p = 0.07). Similar strong trends were observed in the pons where volume was smaller in the GWI group compared to control at 6 and 9 months (p’s = 0.08 and 0.06, respectively).
3.7 Behavior analyses
3.7.1 Open field test
The open field test was used to gather insights into motor function and anxiety-like behavior in these mice at 13 months post study initiation. As expected, all mice habituated to the OF arena as evident by decreases in distance traveled over time [PB/PM model, F(5, 50) = 9.09, p ≤ 0.001; data not shown] and PB/DEET/CORT/DFP model, [F(5,50) = 9.90, p ≤ 0.001, data not shown]. In the PB/PM model, both groups traveled similar distances during the first 5 and total 30 minutes; at 30 min, the slightly greater distance covered by the GWI mice was not significant (p = 0.44) (Figure 9A). In the PB/DEET/CORT/DFP model, the numerically shorter distance that the GWI mice traveled during the first 5 min and the total 30 min was also not significant (p’s ≥ 0.56) (Figure 9C).
Figure 9. Open field test. Motor and mood effects were evaluated using the open field test 13 months post GWI chemicals exposure. Motor function was examined by the distance traveled (A,C) for the first 5 min and total 30 min. Anxiety-like behavior was evaluated by evaluating the (B,D) time spent in the corners of the arena. Data are presented as mean ± SEM. Sample sizes were n = 6/group/timepoint for the PB/PM model and n = 5–6/group/timepoint for the PB/DEET/CORT/DFP model. †Indicates a significant treatment difference, p ≤ 0.05.
Anxiety-like behavior was measured by examining the time spent in the corners of the OF arena. There were no significant impacts on anxiety-like behavior within the first 5 min of the test for either model (p’s ≥ 0.87; Figures 9B,D). Interestingly at 30 min, prior GWI exposure led to more anxiety-like behavior compared to controls in the PB/DEET/CORT/DFP model as these mice spent more time in the corners [t(9) = 2.72, p ≤ 0.05] (Figure 9D); this was not the case in the PB/PM model (p = 0.25) (Figure 9B).
3.7.2 Temperature response stress test
All mice in both models exhibited increases in body temperature indicative of stress-induced (probe insertion) hyperthermia. Within the PB/PM model, there were no significant treatment differences between groups, just significant increases in temperature from the baseline reading for both groups {control [t(10) = −2.73, p ≤ 0.05]; PB/PM [t(10) = −7.28, p ≤ 0.001]; Figure 10A}. Body temperature increased by 1.42 ± 0.28°C and 1.67 ± 0.16°C, respectively, for the control and GWI mice; no significant treatment difference was observed for the change in temperature as well (p = 0.46) (Figure 10A).
Figure 10. Temperature response stress test conducted 13 months post GWI chemicals exposure. The temperature response stress test was used to assess stress-induced hyperthermia following two temperature measurements (baseline and post 15 min). Data are presented as mean ± SEM. Sample sizes were: PB/PM model, n = 6/group/timepoint; PB/DEET/CORT/DFP model, n = 5–6/group/timepoint. *,**, and *** indicate a significant effect of time where p ≤ 0.05, 0.01, and 0.001, respectively. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin; T, temperature.
Similarly, within the PB/DEET/CORT/DFP model, there were no significant treatment differences among groups. Although not significant, at baseline, the PB/DEET/CORT/DFP group had a numerically lower temperature compared to the control group. The control mice exhibited a significant increase in temperature [t(8) = −3.58, p ≤ 0.01], but the GWI group, despite having a lower baseline temperature, had only a numerical elevation in body temperature (p = 0.19) (Figure 10B). Body temperature increased by 1.30 ± 0.21°C and 1.03 ± 0.26°C, respectively, for the control and GWI mice; no significant difference was observed for the change in temperature between treatment groups (p = 0.46) (Figure 10B).
4 Discussion
This study sought to evaluate the progression of structural brain alterations in two GWI preclinical models over the course of 12 months. The results, summarized in Figure 11, highlight similarities and distinctions between these exposure paradigms and aid in the understanding of GWI pathogenesis and endophenotypes. Major similarities among the models include ventricular enlargement and reductions in global brain and hippocampal volumes with age. Key differences were model specific, in which the PB/DEET/CORT/DFP model led to reduced brainstem volumes and an early and persistent loss of total brain volume, while the PB/PM model produced reductions in cortical volume and thickness with age. Behaviorally, at 13 months, motor function was largely preserved in both models. However, the GWI mice in the PB/DEET/CORT/DFP model exhibited an elevation in anxiety-like behavior.
Figure 11. Summary of study results. The effects of two Gulf War Illness models on brain volume over the course of 12 months are shown. Model-specific and shared effects are presented as either increases or decreases in volume at the 6-, 9-, and 12-month time points. *Indicates a treatment difference within the GWI model and is plotted in the direction of the GWI effect. CORT, corticosterone; DEET, N,N-diethyl-meta-toluamide; DFP, diisopropylfluorophosphate; PB, pyridostigmine bromide; PM, permethrin; T, temperature.
Weights were measured for overall health monitoring throughout this study. There were no significant treatment differences over the course of the study (13.5 months). However, PB/PM treated mice weighed numerically less than controls, while PB/DEET/CORT/DFP mice gained more weight over time and trended higher in BMI than controls at the end of the study. These results align with the higher rates of overweight and obesity observed in GWI veterans (Coughlin et al., 2011; Breland et al., 2017), and suggest that prior exposure to this model’s specific chemicals (PB/DEET/DFP) in addition to stress (CORT in this study) may impact weight (e.g., metabolic function) later in life. In fact, it was observed with the temperature response stress test that the GWI mice in the PB/DEET/CORT/DFP model had lower baseline body temperatures and had the smallest elevation in temperature compared to the other groups. This suggests that these mice may have an altered metabolism due to weight gain or potential structural changes within the hypothalamus, including blunted HPA axis (Carmo-Silva and Cavadas, 2017). While this study was MRI-focused and did not dive into metabolic mechanisms or evaluate the hypothalamus, studies show that higher BMIs and obesity are associated with brain structural changes (e.g., reduced volumes; Ward et al., 2005), neuroinflammation, changes in cognitive function, and increased risk of neurological diseases, such as dementias or Alzheimer’s Disease (Whitmer et al., 2008; O’Brien et al., 2017). Further investigations into fat distribution and metabolic dysregulation in a model-specific context are warranted with increased sample sizes appropriate for such investigations.
Many of the effects observed in this study occurred over time, indicating that normal aging is an important factor to consider. In fact, numerous clinical studies have shown that with age, particularly beginning around middle age, global and regional brain volumes decrease in size (Gur et al., 1991; Murphy et al., 1992; Scahill et al., 2003). At the time of study completion (13.5 months post exposures), mice were approximately 15 months old, an age that corresponds to late middle-age in humans and the current status of many GWI veterans. Of note, at the end of the study, there were no differences in either model for gross wet brain weight. However, MRI analysis did reveal longitudinal alterations in global brain volume over time that were most pronounced at the latest scanning (12 months), particularly within the PB/DEET/CORT/DFP model. It has been shown that with age and age-related diseases, such as Alzheimer’s and Parkinson’s, increases in brain water are associated with gray matter atrophy and the breakdown of white matter (Maier-Hein et al., 2015; Ofori et al., 2015; Chad et al., 2018; Gullett et al., 2020). Thus, the lack of changes in brain weight at the end of this study might be an indicator of increased extracellular water compensating for parenchymal decreases. Indeed, significant treatment differences were observed for overall brain volume in the PB/DEET/CORT/DFP model, in which GWI mice exhibited smaller brain volumes compared to control at each time point examined; however, reductions in brain volume within controls did not appear until the 12-month time point. In the PB/PM model, the enlargement of ventricles over time may explain the lack of any discernible change in total brain volume; in fact, within the GWI group, there was a 0.42% increase in total brain volume that may stem from the 25% increase in ventricular size at 12 months. Further, compared to the PB/DEET/CORT/DFP model, the hippocampus was the only brain area examined that significantly decreased over time in the PB/PM model. This suggests the combination of chemicals utilized in this model do not globally impact brain structures, at least at the time points examined. In contrast, in the PB/DEET/CORT/DFP model, the effects observed over the course of the study were largely driven by volume reductions in the GWI group, indicating a treatment dependent effect that may also be compounded by normal aging. Investigations at later stages of life in these models may shed more light on the impact of aging on brain structure in the context of GWI.
It has been shown that ventricular size increases with age (Scahill et al., 2003; Preul et al., 2006) and is associated with declines in brain structure (e.g., cortical thinning) and consequently, neurobehavioral function like cognitive performance (Preul et al., 2006; Carmichael et al., 2007a, 2007b). Here, total ventricular volume was increased at 12 months in both groups of the PB/PM model, albeit to a less degree in the control group (17% control vs. 25% GWI). Similarly, ventricular volume enlarged with age in the PB/DEET/CORT/DFP model. While GWI mice in this model initially had significantly smaller ventricular volume compared to control mice, a pronounced enlargement with age was observed as early as month 9 that persisted to month 12. In fact, ventricular size increased by 20% in the GWI group, almost twice more than the controls. To further characterize where these increases in ventricular volume occurred, ventricles were delineated by the lateral, 3rd, and 4th ventricles and examined. Notably, these measures revealed similarities and distinctions across the evaluated models. Neither model produced significant effects of the 3rd ventricle. However, lateral ventricle size significantly increased with age in both models starting at 9 months post exposures, and was more pronounced in both GWI groups. The changes observed PB/PM model at 12 months align with another GWI study investigating volumetric changes within this region at 10 months (Wu et al., 2021). Additionally, there was an increase in 4th ventricle size within the PB/PM model at the 9-month time point for GWI mice, which was not present at 12 months. This, coupled with the enlargement of the lateral ventricle, suggests a potential effect of neuroinflammation within this region and model. Studies have linked inflammation with transient overproduction of cerebral spinal fluid (CSF) and ventricular enlargement (Lepore et al., 2013; Millward et al., 2020). An earlier study utilizing the PB/PM model found GWI-specific increases in neuroinflammation at this time point as well as behavioral alterations (Carpenter et al., 2021). It is notable that areas containing CSF may contribute to structural changes in gray/white matter (Resnick et al., 2000; Scahill et al., 2003; Preul et al., 2006).
Regional alterations in gray matter were also apparent within these models over time. The PB/DEET/CORT/DFP model produced slight decreases in cortical volume and marginal cortical thinning was observed in both models with age. In the PB/PM model, the GWI increase in ventricular volume and age-associated cortical thinning, particularly in the frontal and dorsolateral entorhinal cortices, could explain the cognitive aberrations seen previously in this model (Zakirova et al., 2015; Zakirova et al., 2016; Carpenter et al., 2021) and GWI veterans. There was a transient decrease in thickness for the frontal cortex at 9 months in the GWI group, which aligns to behavioral aberrations in cognitive performance seen previously in this model at this time point (Carpenter et al., 2021). Interestingly, this effect did not persist beyond 9 months, and warrants further evaluation. In the PB/DEET/CORT/DFP model, thinning of the somatosensory cortex was observed with age and was more prominent in the GWI treated group. This aligns with behavioral outcomes from our earlier studies that identified significant impairments in the sticker removal test (Carpenter et al., 2021; Carpenter et al., 2022). The sticker removal test assesses sensorimotor function, which is negatively affected in veterans with GWI (Axelrod and Milner, 1997; Proctor et al., 2006; Toomey et al., 2009). The present results strengthen the connection between GWI chemicals exposure, cortical impact, and behavior.
Both models saw reductions in hippocampal volumes over time and in all groups, suggestive of an age-related volume reduction. Among models, there was a greater reduction in hippocampal volume for the GWI group in the PB/PM model (11%) than the PB/DEET/CORT/DFP model (2%) at 12 months. In the PB/DEET/CORT/DFP model, both GWI and control groups decreased at similar rates at 12 months. However, in the PB/PM model this reduction was observed at 9 months only in the GWI treated group, indicating that this prior PB/PM exposure impacts this region earlier. These data are in line with findings from GWI studies depicting volume reductions in rats (Wu et al., 2021) and in the GWI veteran population (Chao et al., 2010). Additionally, the significant thinning of the dorsolateral entorhinal cortex over time in this model, particularly within the GWI group, may further compound hippocampal effects, such as cognitive, memory, and synaptic transmission deficits that have been observed in a previous studies utilizing this paradigm (Brown et al., 2021b; Carpenter et al., 2021). Hindbrain atrophy has also been observed in GWI veterans (Christova et al., 2017) and may contribute to the various symptoms experienced by these individuals. Atrophy in the areas of the hindbrain could potentially impact vital functions such as motor control, sleep patterns, and respiration (Moini and Piran, 2020). Here, there were numerical reductions in the hindbrain with age in both models, but effects were more pronounced in the PB/DEET/CORT/DFP model. In the PB/DEET/CORT/DFP model there was a significant treatment difference in volume compared to control at all-time points examined. The hindbrain was delineated into its key regions (e.g., cerebellum, pons, and medulla) to identify any structural changes that may be in line with previously observed behavioral changes in these models and clinical symptomology of veterans with GWI. Interestingly, cerebellar volumes fluctuated within the PB/PM model, and ultimately were numerically reduced at 12 months; a similar reduction was observed in the PB/DEET/CORT/DFP model and these effects are attributed to aging. Further within the PB/DEET/CORT/DFP model, the GWI group started out with significantly smaller cerebellar volumes compared to control. This reduction is in line with clinical observations of cerebellar atrophy in GWI veterans and indicate that this model’s chemicals may contribute to this pathology (Christova et al., 2017). Age-associated declines in volume were also apparent for the pons and medulla within both models, albeit not significant for the PB/PM model. These results are consistent with findings in veterans with GWI (James et al., 2017; Zhang et al., 2020; Zhang et al., 2021). Like the cerebellum, both areas of the brainstem were smaller across the time points examined in a model-specific manner. Here, the GWI group of the PB/DEET/CORT/DFP model had smaller pons and medulla volumes compared to control, as early as 6–9 months. As the brainstem is an important regulatory area for vital functions such as respiration and sleep (Moini and Piran, 2020), further characterization of these functions within these models and in veterans with GWI is warranted.
Overall, this study provides insight into the progressive impact of GWI chemical exposures on brain structures that may be associated with previously reported behavioral and neurological effects observed in these models (Zakirova et al., 2016; Brown et al., 2021a; Carpenter et al., 2021). The data obtained here provide a comprehensive summary of the impact of these models’ chemical exposures, especially when considered alongside our earlier studies, which identified significant behavioral impairments, neuroinflammation, and aberrations in synaptic transmission, particularly at the 9-month time point (Brown et al., 2021a; Brown et al., 2021b; Carpenter et al., 2021; Carpenter et al., 2022). Structural alterations we observed here, especially the markedly earlier total brain volume reduction in PB/DEET/CORT/DFP model are indicative of accelerated brain aging. This is in line with multiple recent reports where (1) the rate of mild cognitive impairment (MCI) was greater in veterans with GWI in two studies (Chao, 2020; Chao et al., 2023), (2) Gulf War veterans not only have greater rates of MCI, but also experience olfactory impairments (Chao, 2024), and (3) veterans with severe GWI symptoms might be impacted by accelerated aging the most (Petry et al., 2024; Thompson et al., 2024). GWI is heterogeneous among veterans; this heterogeneity is attributed not only to varying symptom presentation and severity, but also to diverse chemical exposures during the Gulf War, including exposure to sarin/cyclosarin among a subset of individuals (White et al., 2016). This study provides valuable insights into the structural differences between distinct chemical exposures, such as the earlier PB/PM hippocampal and PB/DEET/CORT/DFP brainstem deficits that may benefit future investigations into targeted neuroprotective interventions based on GWI exposure and symptom endophenotypes.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by the University of Georgia Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
JC: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. SH: Writing – original draft, Writing – review & editing. NF: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Department of Defense (grant nos. W81XWH1610586 and W81XWH2110661) to NMF and the Franklin Foundation Imaging Award to JMC.
Acknowledgments
We would like to extend a special thanks to Kim Mason and William Oliver at the University of Georgia’s Bio-Imaging Research Center for their assistance in this study. Additionally, we thank members of the Filipov Lab, Danielle Ludwig and Alexandria Goodyear, for their assistance in transporting mice to the imaging facility.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2024.1465701/full#supplementary-material
Footnotes
References
Axelrod, B. N., and Milner, I. B. (1997). Neuropsychological findings in a sample of operation desert storm veterans. J. Neuropsychiatry Clin. Neurosci. 9, 23–28. doi: 10.1176/jnp.9.1.23
Borsini, F., Lecci, A., Volterra, G., and Meli, A. (1989). A model to measure anticipatory anxiety in mice? Psychopharmacology 98, 207–211. doi: 10.1007/BF00444693
Bouwknecht, J. A., Hijzen, T. H., Van Der Gugten, J., Maes, R. A. A., and Olivier, B. (2000). Stress-induced hyperthermia in mice: effects of flesinoxan on heart rate and body temperature. Eur. J. Pharmacol. 400, 59–66. doi: 10.1016/S0014-2999(00)00387-3
Breland, J. Y., Phibbs, C. S., Hoggatt, K. J., Washington, D. L., Lee, J., Haskell, S., et al. (2017). The obesity epidemic in the veterans health administration: prevalence among key populations of women and men veterans. J. Gen. Intern. Med. 32, 11–17. doi: 10.1007/s11606-016-3962-1
Brown, K. A., Carpenter, J. M., Preston, C. J., Ludwig, H. D., Clay, K. B., Harn, D. A., et al. (2021a). Lacto-N-fucopentaose-III ameliorates acute and persisting hippocampal synaptic plasticity and transmission deficits in a gulf war illness mouse model. Life Sci. 279:119707. doi: 10.1016/j.lfs.2021.119707
Brown, K. A., Preston, C. J., Carpenter, J. M., Ludwig, H. D., Norberg, T., Harn, D. A., et al. (2021b). Lacto-N-fucopentaose-III (LNFPIII) ameliorates acute aberrations in hippocampal synaptic transmission in a gulf war illness animal model. Brain Res. 1766:147513. doi: 10.1016/j.brainres.2021.147513
Calley, C. S., Kraut, M. A., Spence, J. S., Briggs, R. W., Haley, R. W., and Hart, J. Jr. (2010). The neuroanatomic correlates of semantic memory deficits in patients with gulf war illnesses: a pilot study. Brain Imaging Behav. 4, 248–255. doi: 10.1007/s11682-010-9103-2
Carmichael, O. T., Kuller, L. H., Lopez, O. L., Thompson, P. M., Dutton, R. A., Lu, A., et al. (2007a). Cerebral ventricular changes associated with transitions between normal cognitive function, mild cognitive impairment, and dementia. Alzheimer Dis. Assoc. Disord. 21, 14–24. doi: 10.1097/WAD.0b013e318032d2b1
Carmichael, O. T., Kuller, L. H., Lopez, O. L., Thompson, P. M., Dutton, R. A., Lu, A., et al. (2007b). Ventricular volume and dementia progression in the cardiovascular health study. Neurobiol. Aging 28, 389–397. doi: 10.1016/j.neurobiolaging.2006.01.006
Carmo-Silva, S., and Cavadas, C. (2017). “Hypothalamic dysfunction in obesity and metabolic disorders” in Obesity and brain function. eds. L. Letra and R. Seiça (Cham: Springer International Publishing), 73–116.
Carpenter, J. M., Brown, K. A., Diaz, A. N., Dockman, R. L., Benbow, R. A., Harn, D. A., et al. (2021). Delayed treatment with the immunotherapeutic LNFPIII ameliorates multiple neurological deficits in a pesticide-nerve agent prophylactic mouse model of gulf war illness. Neurotoxicol. Teratol. 87:107012. doi: 10.1016/j.ntt.2021.107012
Carpenter, J. M., Brown, K. A., Veltmaat, L., Ludwig, H. D., Clay, K. B., Norberg, T., et al. (2022). Evaluation of delayed LNFPIII treatment initiation protocol on improving long-term behavioral and neuroinflammatory pathology in a mouse model of gulf war illness. Brain Behav. Immun. Health 26:100553. doi: 10.1016/j.bbih.2022.100553
Carpenter, J. M., Gordon, H. E., Ludwig, H. D., Wagner, J. J., Harn, D. A., Norberg, T., et al. (2020). Neurochemical and neuroinflammatory perturbations in two gulf war illness models: modulation by the immunotherapeutic LNFPIII. Neurotoxicology 77, 40–50. doi: 10.1016/j.neuro.2019.12.012
Chad, J. A., Pasternak, O., Salat, D. H., and Chen, J. J. (2018). Re-examining age-related differences in white matter microstructure with free-water corrected diffusion tensor imaging. Neurobiol. Aging 71, 161–170. doi: 10.1016/j.neurobiolaging.2018.07.018
Chao, L. L. (2020). The prevalence of mild cognitive impairment in a convenience sample of 202 gulf war veterans. Int. J. Environ. Res. Public Health 17:7158. doi: 10.3390/ijerph17197158
Chao, L. L. (2024). Olfactory and cognitive decrements in 1991 gulf war veterans with gulf war illness/chronic multisymptom illness. Environ. Health 23:14. doi: 10.1186/s12940-024-01058-2
Chao, L. L., Abadjian, L., Hlavin, J., Meyerhoff, D. J., and Weiner, M. W. (2011). Effects of low-level sarin and cyclosarin exposure and gulf war illness on brain structure and function: a study at 4T. Neurotoxicology 32, 814–822. doi: 10.1016/j.neuro.2011.06.006
Chao, L. L., Rothlind, J. C., Cardenas, V. A., Meyerhoff, D. J., and Weiner, M. W. (2010). Effects of low-level exposure to sarin and cyclosarin during the 1991 gulf war on brain function and brain structure in US veterans. Neurotoxicology 31, 493–501. doi: 10.1016/j.neuro.2010.05.006
Chao, L. L., Sullivan, K., Krengel, M. H., Killiany, R. J., Steele, L., Klimas, N. G., et al. (2023). The prevalence of mild cognitive impairment in gulf war veterans: a follow-up study. Front. Neurosci. 17:1301066. doi: 10.3389/fnins.2023.1301066
Chao, L. L., Zhang, Y., and Buckley, S. (2015). Effects of low-level sarin and cyclosarin exposure on white matter integrity in gulf war veterans. Neurotoxicology 48, 239–248. doi: 10.1016/j.neuro.2015.04.005
Christova, P., James, L. M., Engdahl, B. E., Lewis, S. M., Carpenter, A. F., and Georgopoulos, A. P. (2017). Subcortical brain atrophy in gulf war illness. Exp. Brain Res. 235, 2777–2786. doi: 10.1007/s00221-017-5010-8
Clarke, T., Jamieson, J. D., Malone, P., Rayhan, R. U., Washington, S., Vanmeter, J. W., et al. (2019). Connectivity differences between gulf war illness (GWI) phenotypes during a test of attention. PLoS One 14:e0226481. doi: 10.1371/journal.pone.0226481
Cooper, C. M., Briggs, R. W., Farris, E. A., Bartlett, J., Haley, R. W., and Odegard, T. N. (2016). Memory and functional brain differences in a national sample of U.S. veterans with gulf war illness. Psychiatry Res. Neuroimaging 250, 33–41. doi: 10.1016/j.pscychresns.2016.03.004
Coughlin, S. S., Kang, H. K., and Mahan, C. M. (2011). Selected health conditions among overweight, obese, and non-obese veterans of the 1991 gulf war: results from a survey conducted in 2003-2005. Open Epidemiol J 4, 140–146. doi: 10.2174/1874297101104010140
Dickey, B., Madhu, L. N., and Shetty, A. K. (2021). Gulf war illness: mechanisms underlying brain dysfunction and promising therapeutic strategies. Pharmacol. Ther. 220:107716. doi: 10.1016/j.pharmthera.2020.107716
Gopinath, K., Gandhi, P., Goyal, A., Jiang, L., Fang, Y., Ouyang, L., et al. (2012). FMRI reveals abnormal central processing of sensory and pain stimuli in ill gulf war veterans. Neurotoxicology 33, 261–271. doi: 10.1016/j.neuro.2012.01.014
Grant, M. K., Bobilev, A. M., Rasys, A. M., Branson Byers, J., Schriever, H. C., Hekmatyar, K., et al. (2020). Global and age-related neuroanatomical abnormalities in a Pax6-deficient mouse model of aniridia suggests a role for Pax6 in adult structural neuroplasticity. Brain Res. 1732:146698. doi: 10.1016/j.brainres.2020.146698
Gullett, J. M., O'shea, A., Lamb, D. G., Porges, E. C., O'shea, D. M., Pasternak, O., et al. (2020). The association of white matter free water with cognition in older adults. Neuroimage 219:117040. doi: 10.1016/j.neuroimage.2020.117040
Gur, R. C., Mozley, P. D., Resnick, S. M., Gottlieb, G. L., Kohn, M., Zimmerman, R., et al. (1991). Gender differences in age effect on brain atrophy measured by magnetic resonance imaging. Proc. Natl. Acad. Sci. USA 88, 2845–2849. doi: 10.1073/pnas.88.7.2845
Haley, R. W., Fleckenstein, J. L., Marshall, W. W., Mcdonald, G. G., Kramer, G. L., and Petty, F. (2000). Effect of basal ganglia injury on central dopamine activity in gulf war syndrome: correlation of proton magnetic resonance spectroscopy and plasma homovanillic acid levels. Arch. Neurol. 57, 1280–1285. doi: 10.1001/archneur.57.9.1280
Haley, R. W., Kurt, T. L., and Hom, J. (1997). Is there a gulf war syndrome? Searching for syndromes by factor analysis of symptoms. JAMA 277, 215–222. doi: 10.1001/jama.1997.03540270041025
Hubbard, N. A., Hutchison, J. L., Motes, M. A., Shokri-Kojori, E., Bennett, I. J., Brigante, R. M., et al. (2014). Central executive dysfunction and deferred prefrontal processing in veterans with gulf war illness. Clin. Psychol. Sci. 2, 319–327. doi: 10.1177/2167702613506580
IOM (2014). Chronic multisymptom illness in gulf war veterans: case definitions reexamined. Washington, DC: The National Academies Press.
James, L. M., Christova, P., Engdahl, B. E., Lewis, S. M., Carpenter, A. F., and Georgopoulos, A. P. (2017). Human leukocyte antigen (HLA) and gulf war illness (GWI): HLA-DRB1*13:02 spares subcortical atrophy in gulf war veterans. EBioMedicine 26, 126–131. doi: 10.1016/j.ebiom.2017.11.005
Koo, B. B., Michalovicz, L. T., Calderazzo, S., Kelly, K. A., Sullivan, K., Killiany, R. J., et al. (2018). Corticosterone potentiates DFP-induced neuroinflammation and affects high-order diffusion imaging in a rat model of gulf war illness. Brain Behav. Immun. 67, 42–46. doi: 10.1016/j.bbi.2017.08.003
Krishna, S., Dodd, C. A., Hekmatyar, S. K., and Filipov, N. M. (2014). Brain deposition and neurotoxicity of manganese in adult mice exposed via the drinking water. Arch. Toxicol. 88, 47–64. doi: 10.1007/s00204-013-1088-3
Lepore, S., Waiczies, H., Hentschel, J., Ji, Y., Skodowski, J., Pohlmann, A., et al. (2013). Enlargement of cerebral ventricles as an early indicator of encephalomyelitis. PLoS One 8:e72841. doi: 10.1371/journal.pone.0072841
Li, X., Spence, J. S., Buhner, D. M., Hart, J. Jr., Cullum, C. M., Biggs, M. M., et al. (2011). Hippocampal dysfunction in gulf war veterans: investigation with ASL perfusion MR imaging and physostigmine challenge. Radiology 261, 218–225. doi: 10.1148/radiol.11101715
Macht, V. A., Woodruff, J. L., Maissy, E. S., Grillo, C. A., Wilson, M. A., Fadel, J. R., et al. (2019). Pyridostigmine bromide and stress interact to impact immune function, cholinergic neurochemistry and behavior in a rat model of gulf war illness. Brain Behav. Immun. 80, 384–393. doi: 10.1016/j.bbi.2019.04.015
Maier-Hein, K. H., Westin, C. F., Shenton, M. E., Weiner, M. W., Raj, A., Thomann, P., et al. (2015). Widespread white matter degeneration preceding the onset of dementia. Alzheimers Dement. 11, 485–493.e482. doi: 10.1016/j.jalz.2014.04.518
Megahed, T., Hattiangady, B., Shuai, B., and Shetty, A. K. (2014). Parvalbumin and neuropeptide Y expressing hippocampal GABA-ergic inhibitory interneuron numbers decline in a model of gulf war illness. Front. Cell. Neurosci. 8:447. doi: 10.3389/fncel.2014.00447
Miller, J. V., Lebouf, R. F., Kelly, K. A., Michalovicz, L. T., Ranpara, A., Locker, A. R., et al. (2018). The Neuroinflammatory phenotype in a mouse model of gulf war illness is unrelated to brain regional levels of acetylcholine as measured by quantitative HILIC-UPLC-MS/MS. Toxicol. Sci. 165, 302–313. doi: 10.1093/toxsci/kfy130
Millward, J. M., Ramos Delgado, P., Smorodchenko, A., Boehmert, L., Periquito, J., Reimann, H. M., et al. (2020). Transient enlargement of brain ventricles during relapsing-remitting multiple sclerosis and experimental autoimmune encephalomyelitis. JCI Insight 5:e140040. doi: 10.1172/jci.insight.140040
Minkeviciene, R., Hlushchenko, I., Virenque, A., Lahti, L., Khanal, P., Rauramaa, T., et al. (2019). MIM-deficient mice exhibit anatomical changes in dendritic spines, cortex volume and brain ventricles, and functional changes in motor coordination and learning. Front. Mol. Neurosci. 12:276. doi: 10.3389/fnmol.2019.00276
Moini, J., and Piran, P. (2020). “Chapter 9 – brainstem” in Functional and clinical neuroanatomy. eds. J. Moini and P. Piran (London: Academic Press), 293–317.
Murphy, D. G., Decarli, C., Schapiro, M. B., Rapoport, S. I., and Horwitz, B. (1992). Age-related differences in volumes of subcortical nuclei, brain matter, and cerebrospinal fluid in healthy men as measured with magnetic resonance imaging. Arch. Neurol. 49, 839–845. doi: 10.1001/archneur.1992.00530320063013
Novelli, E. L., Diniz, Y. S., Galhardi, C. M., Ebaid, G. M., Rodrigues, H. G., Mani, F., et al. (2007). Anthropometrical parameters and markers of obesity in rats. Lab. Anim. 41, 111–119. doi: 10.1258/002367707779399518
O’Brien, P. D., Hinder, L. M., Callaghan, B. C., and Feldman, E. L. (2017). Neurological consequences of obesity. Lancet Neurol. 16, 465–477. doi: 10.1016/S1474-4422(17)30084-4
O’Callaghan, J. P., Kelly, K. A., Locker, A. R., Miller, D. B., and Lasley, S. M. (2015). Corticosterone primes the neuroinflammatory response to DFP in mice: potential animal model of gulf war illness. J. Neurochem. 133, 708–721. doi: 10.1111/jnc.13088
Ofori, E., Pasternak, O., Planetta, P. J., Li, H., Burciu, R. G., Snyder, A. F., et al. (2015). Longitudinal changes in free-water within the substantia nigra of Parkinson's disease. Brain 138, 2322–2331. doi: 10.1093/brain/awv136
Parihar, V. K., Hattiangady, B., Shuai, B., and Shetty, A. K. (2013). Mood and memory deficits in a model of gulf war illness are linked with reduced neurogenesis, partial neuron loss, and mild inflammation in the hippocampus. Neuropsychopharmacology 38, 2348–2362. doi: 10.1038/npp.2013.158
Petry, S. E., Thompson, A. D. Jr., Hauser, E. R., Lynch, S. M., Boyle, S. H., Upchurch, J., et al. (2024). Characterizing deficit accumulation among gulf war era veterans. J. Frailty Aging 13, 300–306. doi: 10.14283/jfa.2024.44
Preul, C., Hund-Georgiadis, M., Forstmann, B. U., and Lohmann, G. (2006). Characterization of cortical thickness and ventricular width in normal aging: a morphometric study at 3 tesla. J. Magn. Reson. Imaging 24, 513–519. doi: 10.1002/jmri.20665
Proctor, S. P., Heaton, K. J., Heeren, T., and White, R. F. (2006). Effects of sarin and cyclosarin exposure during the 1991 gulf war on neurobehavioral functioning in US army veterans. Neurotoxicology 27, 931–939. doi: 10.1016/j.neuro.2006.08.001
Rayhan, R. U., Stevens, B. W., Timbol, C. R., Adewuyi, O., Walitt, B., Vanmeter, J. W., et al. (2013). Increased brain White matter axial diffusivity associated with fatigue, pain and hyperalgesia in gulf war illness. PLoS One 8:e58493. doi: 10.1371/journal.pone.0058493
Resnick, S. M., Goldszal, A. F., Davatzikos, C., Golski, S., Kraut, M. A., Metter, E. J., et al. (2000). One-year age changes in MRI brain volumes in older adults. Cereb. Cortex 10, 464–472. doi: 10.1093/cercor/10.5.464
Ribeiro, A. C. R., and Deshpande, L. S. (2021). A review of pre-clinical models for gulf war illness. Pharmacol. Ther. 228:107936. doi: 10.1016/j.pharmthera.2021.107936
Scahill, R. I., Frost, C., Jenkins, R., Whitwell, J. L., Rossor, M. N., and Fox, N. C. (2003). A longitudinal study of brain volume changes in normal aging using serial registered magnetic resonance imaging. Arch. Neurol. 60, 989–994. doi: 10.1001/archneur.60.7.989
Thompson, A. D. Jr., Petry, S. E., Hauser, E. R., Boyle, S. H., Pathak, G. A., Upchurch, J., et al. (2024). Longitudinal patterns of multimorbidity in gulf war era veterans with and without gulf war illness. J. Aging Health. doi: 10.1177/08982643241245163 [Epub ahead of print].
Toomey, R., Alpern, R., Vasterling, J. J., Baker, D. G., Reda, D. J., Lyons, M. J., et al. (2009). Neuropsychological functioning of U.S. gulf war veterans 10 years after the war. J. Int. Neuropsychol. Soc. 15, 717–729. doi: 10.1017/S1355617709990294
Van Riper, S. M., Alexander, A. L., Koltyn, K. F., Stegner, A. J., Ellingson, L. D., Destiche, D. J., et al. (2017). Cerebral white matter structure is disrupted in gulf war veterans with chronic musculoskeletal pain. Pain 158, 2364–2375. doi: 10.1097/j.pain.0000000000001038
Ward, M. A., Carlsson, C. M., Trivedi, M. A., Sager, M. A., and Johnson, S. C. (2005). The effect of body mass index on global brain volume in middle-aged adults: a cross sectional study. BMC Neurol. 5:23. doi: 10.1186/1471-2377-5-23
Washington, S. D., Rayhan, R. U., Garner, R., Provenzano, D., Zajur, K., Addiego, F. M., et al. (2020). Exercise alters cerebellar and cortical activity related to working memory in phenotypes of gulf war illness. Brain Commun 2:fcz039. doi: 10.1093/braincomms/fcz039
White, R. F., Steele, L., O'callaghan, J. P., Sullivan, K., Binns, J. H., Golomb, B. A., et al. (2016). Recent research on gulf war illness and other health problems in veterans of the 1991 gulf war: effects of toxicant exposures during deployment. Cortex 74, 449–475. doi: 10.1016/j.cortex.2015.08.022
Whitmer, R. A., Gustafson, D. R., Barrett-Connor, E., Haan, M. N., Gunderson, E. P., and Yaffe, K. (2008). Central obesity and increased risk of dementia more than three decades later. Neurology 71, 1057–1064. doi: 10.1212/01.wnl.0000306313.89165.ef
Wu, X., Shetty, A. K., and Reddy, D. S. (2021). Long-term changes in neuroimaging markers, cognitive function and psychiatric symptoms in an experimental model of gulf war illness. Life Sci. 285:119971. doi: 10.1016/j.lfs.2021.119971
Wylie, G. R., Genova, H., Dobryakova, E., Deluca, J., Chiaravalloti, N., Falvo, M., et al. (2019). Fatigue in gulf war illness is associated with tonically high activation in the executive control network. Neuroimage Clin 21:101641. doi: 10.1016/j.nicl.2018.101641
Zakirova, Z., Crynen, G., Hassan, S., Abdullah, L., Horne, L., Mathura, V., et al. (2016). A chronic longitudinal characterization of neurobehavioral and neuropathological cognitive impairment in a mouse model of gulf war agent exposure. Front. Integr. Neurosci. 9:71. doi: 10.3389/fnint.2015.00071
Zakirova, Z., Tweed, M., Crynen, G., Reed, J., Abdullah, L., Nissanka, N., et al. (2015). Gulf war agent exposure causes impairment of long-term memory formation and neuropathological changes in a mouse model of gulf war illness. PLoS One 10:e0119579. doi: 10.1371/journal.pone.0119579
Zhang, Y., Avery, T., Vakhtin, A. A., Mathersul, D. C., Tranvinh, E., Wintermark, M., et al. (2020). Brainstem atrophy in gulf war illness. Neurotoxicology 78, 71–79. doi: 10.1016/j.neuro.2020.02.006
Keywords: Gulf War Illness, magnetic resonance imaging (MRI), pyridostigmine bromide, permethrin, DEET, diisopropylfluorophosphate (DFP), corticosterone, longitudinal
Citation: Carpenter JM, Hughes SN and Filipov NM (2024) Longitudinal evaluation of structural brain alterations in two established mouse models of Gulf War Illness. Front. Neurosci. 18:1465701. doi: 10.3389/fnins.2024.1465701
Edited by:
Ashok K. Shetty, Texas A&M University School of Medicine, United StatesReviewed by:
Kimberly A. Kelly, National Institute for Occupational Safety and Health, United StatesLindsay T. Michalovicz, National Institute for Occupational Safety and Health, United States
Copyright © 2024 Carpenter, Hughes and Filipov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nikolay M. Filipov, ZmlsaXBvdkB1Z2EuZWR1