Abstract
Repeated mild head injuries due to sports, or domestic violence and military service are increasingly linked to debilitating symptoms in the long term. Although symptoms may take decades to manifest, potentially treatable neurobiological alterations must begin shortly after injury. Better means to diagnose and treat traumatic brain injuries requires an improved understanding of the mechanisms underlying progression and means through which they can be measured. Here, we employ a repetitive mild traumatic brain injury (rmTBI) and chronic variable stress mouse model to investigate emergent structural and functional brain abnormalities. Brain imaging is achieved with [18F]SynVesT-1 positron emission tomography, with the synaptic vesicle glycoprotein 2A ligand marking synapse density and BOLD (blood-oxygen-level-dependent) functional magnetic resonance imaging (fMRI). Animals were scanned six weeks after concluding rmTBI/Stress procedures. Injured mice showed widespread decreases in synaptic density coupled with an increase in local BOLD-fMRI synchrony detected as regional homogeneity. Injury-affected regions with higher synapse density showed a greater increase in fMRI regional homogeneity. Taken together, these observations may reflect compensatory mechanisms following injury. Multimodal studies are needed to provide deeper insights into these observations.
Introduction
Repetitive mild traumatic brain injury (rmTBI), sustained through sport, military service or domestic violence, is associated with a myriad of long-term symptoms including impaired cognition, problems with impulse control, aggression, dementia, and suicide (see reviews, Ramalho & Castillo et al. 2015; Ozga et al. 2018; Maresca et al. 2023; Chen et al. 2022). 1 –4 Symptoms can take a decade or more to manifest. Yet, the neuropathological processes which give rise to deleterious long-term outcomes must begin shortly after injury. This latency makes the management of patients with rmTBIs especially challenging. Research efforts are needed to uncover methods for assessing rmTBI outcomes and the efficacy of novel intervention strategies. To this end, comprehensive multimodal approaches using animal models are essential for identifying the neuropathological changes which contribute to neurodegeneration after head injury that are accessible by clinically applicable modalities.
Modeling the biophysical and mechanical forces experienced during rmTBI by humans poses challenges in animal studies. Consequently, a number of models have been developed 5 –10 which, to varying degrees, recapitulate various aspects of rmTBI pathology. 5,9 Given the significance of pathological tau species (herein, p-tau) in some human TBIs, 11 it is imperative for animal models to explore recapitulating this phenotype. Here, we use a mouse model which combines rmTBI, using a closed-head injury (CHI), and stress using a chronic variable stress paradigm (CVS), referred to as rmTBI/Stress henceforth. 12,13 Critically, this combination results in persistent neurocognitive and behavioral deficits, neuroinflammation, and p-tau (which we and others have shown previously, see Tang et al. 2020; Wen et al. 2019; Ojo et al. 2014; Fesharaki-Zadeh et al. 2020). 12,14 –16 While this may not be of direct relevance to injuries sustained through sport, it is likely of high relevance to veterans, and victims of assault and domestic violence. Elevated p-tau in particular has been linked to more severe long-term outcomes. 11,17 –21 Furthermore, CHI does not involve a craniotomy, which enables the inclusion of multiple impacts, and is, overall, a better model of the injuries sustained by humans. 22,23
Typically, studies of rmTBI in animal models use a battery of behavioral assessments and/or immunohistological (IHC) measurements 9 to characterize injury and/or recovery. Here, we use the same model as our group’s previous work, 24 where it was shown that the microlesion at the impact site is restricted to superficial cortical layers, with the perilesional and distal sites remaining unaltered by routine histological stains. The neuronal marker Neuronal Nuclei revealed limited or no decline in neuronal density relative to sham controls. In the present study, we elect to focus on in vivo assessments of injury as these are more readily translatable to patient populations. The use of in vivo imaging modalities, such as magnetic resonance imaging (MRI) and positron emission tomography (PET), in animal models is not as widespread as is behavior and IHC studies, despite their use within the patient population.
Functional MRI (fMRI) analyses from human rmTBI patient studies, such as measures of regional homogeneity (ReHo), at various times after injury, show markedly heterogeneous results in terms of the brain regions and circuits affected. 25 –29 Similarly, heterogeneous alterations in fMRI functional connectivity (FC) have been reported. 30 –34 Although significantly less studied, animal models of rmTBI, at various times after injury (either single or multiple concussions) also report diverse alterations in FC and ReHo. 35 –38 Combining fMRI with another clinically accessible modality, such as PET, may provide further insights into heterogeneous fMRI changes observed across rmTBI animal models.
Combining advanced neuroimaging modalities, such as PET and fMRI, may yield biomarkers that can depict tissue status, structural and functional brain alterations, ultimately aiding in diagnosis of the complex heterogeneous disease process of TBI. Improved characterization of post-injury mechanisms stands to improve rmTBI diagnosis, classification, and prognosis. In an effort to establish increased translational relevance of preclinical studies, we employ two clinically accessible modalities-PET with the synaptic vesicle glycoprotein 2A ligand marking synapse density and fMRI-to assess abnormalities in a mouse model where rmTBI and CVS are combined. To the best of our knowledge, no studies report using the synaptic vesicle glycoprotein 2A (SV2A) marker to assess changes in synaptic density in either humans, or animal models of TBI. Moreover, there are no studies that combine this measure with fMRI.
Our results show widespread decreases in synaptic density (PET) and detect an increase in ReHo (fMRI) in brain regions with decreased synaptic density. SV2A-PET standardized uptake value ratios (SUVRs) and BOLD (blood-oxygen-level-dependent) fMRI ReHo show a positive correlation in healthy mice (across the whole brain). With injury, this relationship (SV2A-PET SUVRs vs. BOLD-fMRI ReHo) strengthens. Moreover, separate from the impact foci and in areas without any physical damage (based on structural MRI), the affected brain regions in injured mice with higher synaptic density show a greater proportional increase in ReHo than do regions with lower synapse density. FC changes between regions with decreased synapse density are also reported. Taken together, these observations may reflect partial spontaneous neural circuit compensation for injuries sustained. Further multimodal studies are needed to provide clarification and deeper insight into these observations.
Methods
Our study exclusively examined male mice. Thus, due to the known differences between males and females in TBI (see review, Biegon 2021), 39 it is unknown whether the findings presented here are relevant for female mice. All experiments and procedures are conducted following the Institutional animal care and use committee (IACUC) approved protocols at Yale University. Wildtype male C57BL/6J mice (n = 20) from Jackson Labs (stock number: 000664) at 16 weeks of age are used. All mice are kept in standard housing under a 12 h light/dark cycle with food and water provided ad libitum throughout the experiment (with the exception of overnight food deprivation which was one of the possible stress exposures included in our CVS paradigm, see below).
Repetitive mild closed-head injury combined with chronic variable stress
We combine rmTBI, using a CHI, and CVS (stress), referred to as rmTBI/Stress hereafter. 12 CVS is induced as described previously (n = 10). 12 Briefly, CVS consists of 2-weeks of daily exposure to 3/5 aversive stimuli: (1) 3-min cold water swim (16–18°C), (2) overnight food deprivation with access to water ad libitum, (3) 3-h in a cage with 300 ml of water added (4) 3-hour exposure to a cage tilted at a 45° angle, and (5) 15-min immobilization in a flat-bottomed restraint chamber (Braintree Scientific Instruments). Aversive stimuli are randomized in each session for each mouse to simulate the unpredictable nature of psychological trauma, and limit habituation. Sham-stress mice are transferred to a single occupancy cage for the duration of the stressor. Daily rmTBI occurs within 1–3 h of stress exposure.
Prior to each injury, anesthesia was induced with 3.5% isoflurane (in O2, 1.0L/min) and maintained at 3% until immediately after impact when anesthesia is discontinued (<10 min). The head is shaved, the animal placed in a stereotaxic frame (Fig. 1), and rmTBI induced using a 5 mm diameter tip operated by an electromagnetic impactor (Leica Microsystems, Buffalo Grove, IL). The impactor tip is placed 5 mm lateral of the sagittal line, 5 mm caudal of the eye, and at 20° from vertical (Fig. 1). The impact is at a velocity of 5 m/sec, depth 1 mm and dwell time 100 ms. Impacts alternate between the right and left hemisphere on consecutive days to minimize the probability of hemorrhage.

Timeline. Mice (n = 10) are subjected to daily random stressors for 2 weeks (chronic variable stress, CVS), followed immediately by a mild brain injury (repeated daily) over a 2-week period. Six-weeks after concluding rmTBI/Stress procedures, PET data and MRI data are collected (1-week apart). For all methodological details (e.g., sequences used and anesthesia), see Methods. MRI, magnetic resonance imaging; PET, positron emission tomography; rmTBI, Repetitive mild traumatic brain injury.
AZD0530 treatment
A subset of 5 injured mice receive AZD0530 treatment (saracatinib), while the rest are given the vehicle. Amongst tau-interacting proteins is the neuronally-enriched cytoplasmic tyrosine kinase, Fyn, a member of the Src family. In preclinical models of Amyloid-ß-triggered deficits, Fyn inhibition with the orally available kinase inhibitor AZD0530 (saracatinib) reduces synaptic loss, memory deficits and p-tau. 40 Mice began treatment (or vehicle) 24-h following the last day of rmTBI/Stress. 12 Treatment doses of 5 mg/kg/day are given in two equal aliquots by oral gavage until the study end-point. The vehicle is 0.5% weight/volume (wt/vol) hydroxypropyl-methylcellulose per 0.1% wt/vol polysorbate 80. Each dose is 250 µl. Treatment effects are not detectable (possibly due to a subtherapeutic low dosage). All 10 injured mice (regardless of treatment-status) are grouped together in our analyses (but Figures include shaded dots to indicate treated vs. vehicle). Independent two-tailed t-tests are used to compare the synaptic density changes between treated and untreated animals for a number of different brain regions (Supplementary Fig. S2).
Control group
n = 3 mice undergo sham-rmTBI and sham-stress. Sham-stress mice are transferred to a single occupancy cage for the duration of the stressor. Sham-CHI mice are anesthetized and have their heads shaved but do not undergo impaction. To address this small sample size, n = 7 mice are added which underwent identical imaging procedures but did not undergo sham-CHI or sham-CVS. To assess the similarity of sham and naïve mouse datasets, we compare whole-brain FC between each mouse and the group average using Spearman’s rho (Supplementary Fig. S1). No notable deviations are found. Thus, all 10 uninjured mice (regardless of sham-status) are grouped together in our analyses (but Figures. include shaded dots to indicate sham vs. naïve).
Head-post surgery
To reduce motion and susceptibility artifacts during MRI with low-dose (0.5% isoflurane) anesthesia, all mice undergo a minimally invasive surgery where a custom head-plate is affixed to the skull as we have described previously. 41 –43 Mice undergo head-plate implantation following the last day of rmTBI/Stress. Briefly, mice are initially anesthetized with 5% isoflurane (70/30 medical air/O2) and secured in a stereotaxic frame (KOPF) using ear bars and an incisor bar (∼3 min). Isoflurane is then reduced to 2%. Ointment is applied to the eyes to prevent dryness; meloxicam (2 mg/kg body weight) is administered subcutaneously, and bupivacaine (0.1%) injected locally at the incision site. After hair removal (Nair, chemical hair remover), the scalp is washed 3 times with betadine followed by ethanol 70%. The skin and soft tissue overlying the skull is surgically removed. The skull is cleaned and dried. Antibiotic powder (Neo-Predef) is applied to the incision site, and isoflurane reduced further to 1.5% (∼10 min). Superglue (Locite) is applied to the exposed skull, followed by transparent dental cement C&B Metabond (Parkell). The pre-cut head-plate is attached to the dental cement prior to its solidification. The entire procedure is concluded in ∼20 min, recovery takes <10 min.
PET data acquisition, processing, and analyses
Data acquisition and processing
All mice undergo [18F]SynVesT-1 PET scans after completing treatment, or 6-weeks after the last TBI (Fig. 1), on a Focus 220 scanner (Siemens Medical Solutions). A 57Co transmission scan is collected for attenuation correction. During acquisition, mice are maintained on isoflurane anesthesia (1.5%–2.5%). Three mice are placed in the scanner during each PET acquisition. [18F]SynVesT-1 is administered via intramuscular injection (<0.2 mL). The injected radioactive doses were 25.0 ± 10.1 MBq. After acquisition, anesthesia is discontinued, and recovery monitored closely (mice recover in <10 min). Thus, the entire imaging procedure lasts <90 min.
Images are reconstructed with a 3D ordered subset expectation maximization method (3D-OSEM; 2 iterations, 9 subsets) with maximum a posteriori probability algorithm (MAP; 18 iterations) with corrections for decay, random effects, attenuation, and scatter. Images of SUVs (i.e., activity normalized to injected dose and body weight, or voxel value in the reconstructed image (Bq/mL) · body weight (g)/injected dose in Bq) are generated using data from 30–60 min post-injection. SUVRs are normalized to striatum. This region was selected based on the effect size for group differences (controls vs. injured) was the lowest in this region. One of the major challenges in rodent PET studies is quantification. In this study, we utilized the averaged SUV from 30–60 min normalized by the striatum. The optimal reference region in this study would be an area unaffected by CVS or rmTBI. Although previous literature indicates that TBI might induce presynaptic terminal degeneration in the striatum, 44 our findings show that the striatum SUV (normalized by injected activity and body weight only) exhibited the smallest group differences on average and lowest effect size. Based on those findings, we hypothesize that normalization by striatum values has the smallest impact on the group differences while minimizing the data variability.
Analyses
FSL General Linear Model (GLM) is used to perform a statistical comparison of SUVRs across the entire brain between injured and control mice. We perform nonparametric permutation testing with 5,000 permutations, using family-wise error correction with threshold-free cluster enhancement (TFCE). Statistical significance is defined as p < 0.05. Based on the significant results obtained from this group comparison, we create a group of Region of interest (ROIs) (Fig. 3B). Specifically, we define regions as the intersect between voxels we determine to show synapse losses (SV2A-PET) and regions within the Allen Brain atlas. 45 The cerebellum, olfactory bulb and medulla are excluded from fMRI analyses based on limited coverage. This approach results in 10 regions, which we use in all other analyses (unless specified).
MRI data acquisition, processing, and analyses
Data acquisition
Data are acquired on an 11.7T system (Bruker, Billerica, MA), using ParaVision 6.0.1 software. Body temperature is continuously monitored (Neoptix fiber) and maintained with a circulating water bath at 36.6–37°C. BOLD-fMRI data are acquired using a gradient-echo, echo-planar-imaging (EPI) sequence with a repetition time (TR) of 1.8 sec, and echo time (TE) of 9.1 ms. Data are isotropic 0.275 × 0.275 × 0.275 mm3, with 35 slices yielding near whole-brain coverage (no gaps). Each functional imaging run is comprised of 334 repetitions (∼10 min). In total, 30-min of resting-fMRI data are acquired during each imaging session (i.e., from each mouse).
Structural MRI
These data are used for assessment of physical damage (Fig. 2) and registration to a common space as we have described previously. 41,43 (1) A high in-plane resolution image of the fMRI field-of-view. Using a MSME (multi-spin, multi-echo) sequence, in 7 min 20 s, using a TR/TE of 2500/20 ms, we obtain 35 slices (0.275 mm thick) with an in-plane resolution of 0.1 × 0.1 mm2. The slice prescription matches those of the fMRI data (i.e., they are of the same anatomy). (2) Isotropic 3D anatomy of the whole brain. Using a MSME imaging sequence, in 5 min 20 s, using a TR/TE of 5500/15 ms, we obtained a 0.2 × 0.2 × 0.2 mm3 (single average) image of the whole brain.

Physical damage (edema) from repetitive mTBI/stress.
After structural and functional MRI data have been acquired, anesthesia is discontinued, and recovery monitored closely (mice recover in <10 min). Thus, the entire imaging procedure lasts <60 min.
Data processing
Functional MRI data are processed using the open-source RABIES (Rodent Automated BOLD Improvement of EPI sequences) software package (https://github.com/CoBrALab/RABIES). 46 A volumetric EPI is derived from trimmed average EPI frames after motion realignment. Using this volumetric EPI as a target, the head motion parameters are estimated by realigning each EPI frame to the target using a rigid registration. To conduct common space alignment, structural images, which are acquired alongside the EPI scans, are initially corrected for inhomogeneities and then registered together to allow the alignment of different MRI acquisitions. This registration is conducted by generating an unbiased data-driven template through the iterative nonlinear registration of each image to the dataset consensus average. The average is updated at each iteration to provide an increasingly representative dataset template. 47 The finalized template, after the last iteration, provides a representative alignment of each MRI session to a template that shares the acquisition properties of the dataset (e.g., brain shape, Field of view, and anatomical contrast), which makes it a stable registration target for cross-subject alignment. After aligning the MRI sessions, this newly generated unbiased template is then itself registered, using a nonlinear registration, to an external reference atlas to provide both an anatomical segmentation and a common space comparable across studies defined from the provided reference atlas. To correct for EPI susceptibility distortions, the volumetric EPI is also subjected to inhomogeneity correction, and then registered using a nonlinear registration to the anatomical scan from the same MRI session, which allows the calculation of the required geometrical transformations for recovering brain anatomy. 48 Finally, after calculating the transformations required to correct head motion and susceptibility distortions, slice timing correction is applied to the timeseries (3dTshift). 49 Preprocessed timeseries in common space are generated by concatenating the transforms allowing resampling to the reference atlas (the Allen Brain Atlas), 45 at a voxel resolution of 0.2 × 0.2 × 0.2 mm3.
Confound correction of EPI timeseries data is conducted in common space. First, frames with prominent corruption are censored (i.e., scrubbing, see Power et al. 2012). 50 Framewise displacement 50 is measured across time and each frame surpassing 0.075 mm of motion, together with 1 backward and 2 forward frames, are removed. Next, voxelwise detrending is applied to remove first-order drifts and the average image. Selected nuisance regressors are then used for confound regression. More specifically, using ordinary least square regression, the 6 rigid motion parameters, the mean signal from the cerebral spinal fluid mask and the global signal are modeled at each voxel and regressed from the data. Finally, a spatial Gaussian smoothing filter(nilearn.image.smooth_img) 50,51 is applied at 0.3 mm full-width at half maximum (FWHM).
Analyses
Physical damage
Lesion volume is assessed by manual segmentation of MSME MRI data. Voxels with signal intensities above or below average grey/white matter are delineated in native space using tools developed for thispurpose within BioImage Suite (BIS, https://bioimagesuiteweb.github.io/alphaapp/index.html).
ReHo
Voxelwise ReHo maps are computed using AFNI. 49,52 ReHo is calculated for a given voxel and its 27 nearest neighbors and smoothed by a Gaussian kernel with FWHM equal to 10.
fALFF
Voxelwise fractional Amplitude of Low Frequency Fluctuations (fALFF) is assessed as a ratio of BOLD-fMRI activity in the low frequency range (0.01–0.075 Hz) to the entire frequency range (0.01–0.03 Hz). 53,54
FC strength
Correlation matrices are obtained using the whole-brain Allen atlas (Supplementary Fig. S2) or ROIs generated from SV2A-PET data (Fig. 3B). Pearson’s correlation is computed between all region pairs.
Statistics
FSL GLM (Functional Magnetic Resonance Imaging of the Brain Software Library) 55 is used to perform statistical comparison analyses between injured and uninjured mice. We perform nonparametric permutation testing with 5,000 permutations, using familywise error correction with TFCE. Statistical significance is defined as p < 0.05. Additionally, independent t-tests are used to ascertain differences between groups. Specifically, average SUVRs (Fig. 3C), ReHo (Fig. 4C) and fALFF (Supplementary Fig. S5), which are adjusted using FDR corrections for multiple comparisons. 56

Widespread decreases in synapse density in “spared” brain areas.

Increased BOLD-fMRI ReHo in regions with depleted synapse densities.
Data and code are available upon request.
Results
We used daily repeated mTBI in combination with stress exposure over 14 days weeks as a mouse model of combined brain injury and chronic stress (Fig. 1). 16 Six-weeks after injury, mice underwent two scan sessions, one using the PET tracer, [18F]SynVesT-1, a ligand of SV2A as a pan-neural synapse density marker, 57 –63 and a second scan session using MRI (Fig. 1). Our goal was to investigate the emergent functional abnormalities that follow injury and stress using clinically accessible in vivo imaging modalities—and the intermodal relationship of their complementary contrasts.
Mild and heterogeneous tissue damage or edema in cortical areas at or near the impact sites
Our goal was to induce a mild injury with limited destruction of brain tissue. To assess our success, we quantified changes in T2-contrast using MSME-MRI (Methods). Voxels with high intensities (likely indicative of edema 64 ) were segmented manually by a blinded observer (no hypo-intensities were noted). Figure 2A displays 3D reconstructions (BIS, https://bioimagesuiteweb.github.io/alphaapp/index.html) of each injured brain. An example coronal image, from each mouse highlights the area where the most evidence of damage was observed. Overall, these data indicated minimal destruction of brain tissue, and an appreciable amount of heterogeneity within the group — akin to the human TBI patient population. Figure 2B summarizes these data. A maximum of four out of 10 mice had overlapping physical damage (i.e., the same voxels were identified across mice), confirming the heterogeneity within the dataset and the mild nature of the induced injuries (on average). Mostly cortical regions at or proximal to the impact sites were affected (as expected).
A widespread decrease in synapse density followed injury in ‘spared’ brain areas
Six-weeks after injury (Fig. 1), we assessed changes in synapse density via SV2A-PET. Specifically, we used a PET tracer, [18F]SynVesT-1, which targets SV2A (Methods). This marker correlated well with in vitro histopathological synapse markers and thus allowed for an in vivo evaluation of synapse density. 57 –63 PET SUVRs, using striatum as the reference region, from injured and control mice were compared at the voxel-level (FSL, GLM). Figure 3A shows widespread decreases in synapse density in regions that showed no physical damage (on MRI) in the injured relative to the control group (p < 0.005; TFCE; corrected; uncorrected results are shown in Supplementary Fig. S3). Based on the intersection of affected voxels and the Allen Brain Atlas ontology 45 (Fig. 3B), we compared average SUVRs from each mouse (Fig. 3C, independent t-test two-tailed, FDR corrected, p < 0.005). Reduced synapse density was observed bilaterally in olfactory areas, thalamus, and midbrain. Unilateral decrements were more prominent on the right side and observed in visual and visceral areas, hippocampus, agranular insular dorsal/ventral areas, and the cortical amygdala area.
Increased BOLD-fMRI ReHo in areas with depleted synapse densities
We assessed the effects of rmTBI/Stress on fMRI measures of ReHo, a metric of local BOLD synchrony (correlation strength among neighboring voxels). 52 Specifically, within regions with synapse loss (as defined in Fig. 3B), ReHo was computed for each voxel and its 27 neighbors (Fig. 4A), and compared between injured and control groups (FSL, GLM). The olfactory bulb was excluded due to a lack of coverage in fMRI data. Significant but uncorrected changes in ReHo are shown in Supplementary Figure S4. Across injured mice, a significant (p < 0.05, corrected) increase in ReHo was observed relative to the control group within the thalamus, midbrain, left hippocampus and right visual area (Fig. 4B, C).
To assess changes in regional excitability, we computed voxelwise the fALFF 53,54 as the ratio between BOLD activity measured in the low frequency (0.01–0.075 Hz) and entire frequency (0.01–0.03 Hz) range (Supplementary Fig. S5). This metric is thought to capture neural excitability 53,54,65,66 and was thus one way to assess regional (hyper)excitability. Our results did not show a significant change in fALFF, although a statistically non-significant trend towards an increase in fALFF was observed in the hippocampus and thalamus (p < 0.05, uncorrected) of injured mice (Supplementary Fig. S5).
In injured animals, regions with higher synapse densities showed a greater increase in local BOLD-fMRI synchrony
Our multimodal data offered us the unique opportunity to examine inter-contrast relationships. We first assessed how synapse density (SUVR) and ReHo were related in control mice, i.e., without any exposure to rmTBI/Stress. Average SUVR and ReHo were extracted from regions following the Allen Brain Atlas. 45 Results indicated a moderate but significant positive correlation (Pearson’s rho = 0.582; p < 0.001) between ReHo and SUVR in the healthy condition (Fig. 5A). This means that the regions with higher SUVR showed greater ReHo.

Relationship between synapse density (SV2A-PET SUVRs) and ReHo (BOLD-fMRI).
Next, we assessed the relationship in injured mice (n = 10). Average SV2A-PET SUVRs and BOLD-fMRI ReHo values were extracted from n = 10 regions (Methods) as defined in Figure 3B. Between measures, a positive correlation (Pearson’s rho =0.582; p < 0.001) was uncovered (Fig. 5A). This was compared, using the same regions, to the control group. A strong and statistically significant correlation (Pearson’s rho = 0.812; p = 0.004) was found in the injured group (Fig. 5B). A statistically insignificant trend towards a moderate positive correlation (Pearson’s rho = 0.547; p = 0.101) was uncovered in controls (likely due to our small ‘N’ and effect size). Further, a difference between slopes across groups was observed (t-value = 2.15; Dfd = 16; p value, two-tailed = 0.047). In summary, SV2A-PET SUVRs, and BOLD-fMRI ReHo showed a positive correlation in healthy mice (across the whole brain). With injury, this relationship (between SV2A-PET SUVRs and BOLD-fMRI ReHo) strengthened in “spared” brain regions. Moreover, separate from the impact foci and in areas without any physical damage on structural MRI, variation in synapse density was highly correlated with an increase in ReHo, such that regions with higher synapse densities showed a greater increase in ReHo after injury.
Sparse indications of increases in BOLD-fMRI FC were uncovered following injury
In addition to assessing local BOLD-fMRI synchrony within regions, we assessed injury-induced changes in inter-regional FC strength. 67 We assessed brain-wide FC calculated based on the Allen Brain atlas 45 (Supplementary Fig. S6A). Permutation testing (5,000 permutations) indicated no significant (uncorrected) differences between groups. However, among ‘spared’ regions showing decreases in synapse densities (Fig. 3), FC strength (z-scored) in injured relative to control mice (Fig. 6A) showed a decrease between the right agranular and visual areas, as well as more widespread increases including between gustatory and visual areas, thalamus and agranular areas and left and right thalamus (FSL, GLM) (Fig. 6B). These analyses were exploratory (p < 0.005, uncorrected), as our observations did not survive FDR correction for multiple comparisons.

Patterns of injury-induced changes in BOLD-fMRI FC strength.
Discussion
rmTBI can lead to debilitating protracted neurological and behavioral symptoms that can have a delayed onset. Clinical standardized treatment methods to properly manage patients, and effectively monitor and treat rmTBI symptoms and potential progression are lacking. To address these gaps, preclinical studies can play a crucial translational role. Here, we use a mouse model which combines rmTBI and CVS to recapitulate key phenotypes of rmTBI/Stress. To this end, we and others have shown that stress exacerbates mild mechanical trauma 15,16 and subsequent tauopathy. 12
Six-weeks post-injury, mice undergo PET with an SV2A ligand and MRI. With a focus on translational relevance, we use clinically assessable imaging modalities to characterize injury-induced brain abnormalities. Injured mice show limited and heterogeneous physical damage (MRI), but widespread decrease in synapse density—including in regions “spared” from physical damage (on structural MRI)—coupled with an increase in BOLD-fMRI ReHo. In uninjured mice, a positive correlation is found between SV2A-PET SUVRs and BOLD-fMRI ReHo. This relationship is enhanced following injury, suggesting that brain areas with higher synapse densities show a greater increase in ReHo post-injury. Additionally, sparse increases in FC strength are observed. Taken together, these results may reflect endogenous functional compensation for injury (more below).
A widespread decrease in synapse density accompanies an increase in BOLD-fMRI ReHo
SV2A is ubiquitously expressed in presynaptic vesicles, making it a suitable in vivo biomarker for synapse density measurements. 68,69 Recent validation studies in humans 61,70 –73 and nonhuman primates 74,75 include pioneering work measuring synapse losses in patients with temporal lobe epilepsy, 71 and severe depression. 73 Studies in rodents, validated with in vitro immunoreactivity 40,69 show synapse losses in a mouse model of Alzheimer’s disease (AD), 57,58 and in rat cortical and hippocampal regions following a single head injury. 76 Our findings align well with previous SV2A-PET studies. After rmTBI/Stress, mice show a decrease in synapse density, not only in cortical and hippocampal areas, but also in the cerebellum, midbrain, olfactory areas, and thalamus.
To assess how decreases in synapse density affect local connectivity we measure BOLD-fMRI synchrony using ReHo. 52 In patients with TBI, ReHo, at various times after injury, shows markedly heterogeneous results in terms of the brain regions and circuits affected, as well as whether ReHo goes up or down (Table 1). 25 –30,77,78 Given the heterogeneity in time after injury, mechanism of injury (e.g., military service or sports-related), and injury location, there is a clear need for further research.
Summary of BOLD-fMRI ReHo Measurements in TBI Patients
BOLD, blood-oxygen-level-dependent; FC, functional connectivity; fMRI, functional magnetic resonance imaging; L, left; R, right; mid, middle; ReHo, Regional homogeneity; Sup, supplementary; OFC, orbital frontal cortex; TBI, traumatic brain injury.
Experimental animal models of TBI are likely to improve our understanding of brain injury progression, and its imaging correlates, given our ability to control injury mechanism, location, frequency, and severity. Studies on rodent TBI models measuring BOLD-fMRI ReHo are just beginning to emerge. We identified two in the literature. The first is a longitudinal study of mice following a single CHI. 79 Despite the difference in model (rmTBI/Stress vs. a single CHI), and timepoint (6 weeks vs. 2 days and 2 weeks) our findings align. To et al. observe an increase in ReHo in similar brain regions including the hippocampus, thalamus, and midbrain. 79 On the other hand, a more recent study on drug seeking mice at a later timepoint (7 weeks post-injury), using a different injury mechanism (3 repeated blast injuries), uncovers very sparse (or negligible) increases in ReHo, but more widespread increases when the interaction between injury and drug seeking behavior is considered. 80 Given the potential for translation to the clinic, future work is needed to determine the effects of injury mechanism, and time, on BOLD-fMRI measures of ReHo in animal models of (rm)TBI.
Increases in ReHo could be related to p-tau, neuroinflammation, and/or neural hyperexcitability. Our previous work in a CHI mouse model 16,18,24 uncovered p-tau and neuroinflammation similar to observations made in animal models of AD. 12,58,81 In AD, we and others have established a link between p-tau, neuroinflammation, and synapse losses measured via SV2A-PET. 58,72,82 P-tau and neural hyperexcitability have also been reported in mouse models of AD 83,84 and brain injury, 37 alongside increases in astrocytic activity. 85 In AD, an increase in hippocampal ReHo, a positive correlation between ReHo and tau-immunoreactivity, and evidence of neural hyperexcitability, have been reported. 86 Taken together, accumulating evidence seems to suggest that brain injury or AD can elicit p-tau coupled with an immune response and neural hyperexcitability which in turn manifests as an increase in BOLD-fMRI ReHo. To & Nasrallah proposed that hyperexcitability is responsible for their observations of increased ReHo and BOLD-fMRI FC, in a single concussion TBI model. 38 One interpretation of increased ReHo and neural hyperexcitability is that it reflects evidence of endogenous compensation for injury (or disease)—neurons attempting compensate which could result in greater local synchrony. Future studies are needed to dig deeper into these speculations including in vivo and ex vivo markers of brain activity and neurodegeneration (ideally in the same subjects).
To provide preliminary exploration of the hypothesis that increases in ReHo reflect neural hyperexcitability, we compute voxel wise fALFF 53,54 (Supplementary Fig. S5). This metric is thought to capture neural excitability 53,54,65,66 and is thus one way to probe our hypothesis. Encouragingly, TBI patients at three months post-injury, show increases in fALFF within regions that that show increases in ReHo, and there is a significant correlation between the two metrics. 30 However, to the best of our knowledge, fALFF has not yet been reported in an animal study of brain injury. While our results do not show a significant change in fALFF, a statistically non-significant trend towards an increase in fALFF is observed in the hippocampus and thalamus (p < 0.05, uncorrected) of injured mice (Supplementary Fig. S5). Further studies are warranted to reliably assess changes in fALFF caused by brain injury and to explore relationships between BOLD-fMRI measures (ReHo, FC, and fALFF).
To the best of our knowledge, this is the first study to compare SV2A-PET and BOLD-fMRI measurements in the same brain. Thus, we begin by quantifying the relationship between these metrics in healthy (uninjured) mice. We find a robust positive correlation between SV2A-PET SUVs and BOLD-fMRI ReHo across the whole-brain. This suggests that BOLD-fMRI ReHo is (to some extent) reflective of underlying synapse densities. In brain regions susceptible to injury-induced synapse loss, this correlation strengthens with injury and ReHo changes are greatest in regions with greater synapse density. The brain regions with high excitatory synapse densities include hippocampus, midbrain, cerebellum, and thalamus. 87 The dependence of injury-induced ReHo on baseline synapse density is consistent with the hypothesis of increased post-injury excitability. This hypothesis requires further exploration which takes into account both pre- and post-injury measures.
Although we observe sparse changes in BOLD-fMRI FC strength that do not survive corrections for multiple comparisons, the circuit showing the greatest numerical change aligns well with underlying anatomical connectivity. The gustatory cortex projects to the hippocampus, amygdala, visceral and agranular (insular) cortex, thalamus, and visual cortex. Insular areas project to the amygdala, while agranular areas project to parts of the thalamus and midbrain (Supplementary Fig. S7). Collectively this circuit is responsible for touch, taste, and olfaction. 88 –90 Encouragingly, a similar pattern of increased FC strength at eight weeks after a single controlled cortical impact has recently been reported. 91 Given the collective role of the injury-affected circuit in sensory processes, it would be intriguing to study further via coupled behavioral and multimodal imaging experiments.
For human TBI, BOLD-fMRI measurements of FC strength show pronounced heterogeneity, matching the variability of ReHo (Table 1). 25 –30,77,78 This likely reflects the diverse nature of clinical (rm)TBI studies. Yet, relative to BOLD-fMRI ReHo, the FC literature is considerably more developed (see reviews, Lunkova et al. 2021; Mayer et al. 2015; Jantzen et al. 2010; Eierud et al. 2014). 31 –34 Of particular relevance to the present study, brain-wide increases in FC strength, alongside increases in ReHo and fALFF have been recently reported. 30
Study limitations
Despite employing rigorous statistical approaches to mitigate the impact of a small sample size, our ability to draw definitive conclusions, especially in fALFF and FC analyses, is limited. Our control group is comprised of both sham and naïve mice (Methods). Although we took careful steps to ensure that sham and naïve mice show no differences across measures (see Methods and Supplementary Fig. S1), a group of sham-only mice would have been preferable. A subset of our injured mice received AZD0530 treatment (saracatinib) or a vehicle as part of a pilot study. Despite the potential neuroprotective effects of AZD0530 in a mouse model of AD, 40 we do not observe any differences between treated mice and those given a vehicle (Supplementary Fig. S2). Consequently, we consolidate all injured mice into one group. This study assesses mice with clear yet mild physical brain damage (on MRI), while some studies only include animals with no evidence of structural alterations. 35,68 Our view is that the heterogeneity and presence of structural damage here provide a closer approximation to the human condition. 92
While this variability may mirror the inherent diversity in patients, it is critical to acknowledge that such variability may also indicate potential limitations within our model of injury—especially given our small sample size. Further, to address potential irregularities, we exclude all areas with evidence of physical damage from our analyses. Finally, this study is conducted on male mice only. Given the known differences between males and females in TBI, 39 future work will include both male and female animals.
Conclusions
This study utilizes a mouse model of rmTBI/Stress 93 over a course of two-weeks. Six-weeks post-injury, mice undergo SV2A-PET and MRI imaging sessions. Our focus is on clinically accessible modalities (PET and fMRI) for assessing emergent abnormalities indicative of risk of neurodegeneration. Our findings reveal widespread decreases in synapse density and increases in BOLD-fMRI ReHo in areas without physical damage. A moderate positive correlation is found between SV2A-PET SUVRs and BOLD-fMRI ReHo measures in uninjured mice. In mice with brain injuries, greater increases in BOLD-fMRI ReHo are observed in regions with higher synapse density. Taken together, these observations may reflect compensatory mechanisms following injury. Further multimodal studies are needed to provide deeper insights into these observations.
Transparency, rigor and reproducibility statement
All experiments and procedures adhere to the protocols approved by the Institutional Animal Care and Use Committee (IACUC-2023–20284) at Yale University. Twenty-one wildtype male mice are randomly allocated into control and injured groups, guided by prior research from the laboratory, <10 animals per group is sufficient to detect behavioral and brain changes. Repetitive mild traumatic brain injury (rmTBI) with stress is administered daily in a randomized manner. Specifically, rmTBI comprises a two-week regimen involving daily exposure to three out of five aversive stimuli: (1) a 3-min cold water swim (16–18°C), (2) overnight food deprivation with access to water ad libitum, (3) 3 h in a cage with 300 ml of water added, (4) 3 h exposed to a cage tilted at a 45° angle, and (5) 15 min of immobilization in a flat-bottomed restraint chamber (Braintree Scientific Instruments). The aversive stimuli are randomized to mimic the unpredictable nature of psychological trauma and to limit habituation. Sham-stressed mice are transferred to a single-occupancy cage during the stressor. Daily rmTBI occurs within 1–3 h of stress exposure (further details in the Methods under “Repetitive Mild Closed-Head Injury with Chronic Variable Stress”). A subset of five randomly selected mice receive a sub-therapeutic dose of saracatinib treatment daily for six weeks at doses of 5 mg/kg/day, administered in two equal aliquots via oral gavage until the study end-point. Independent t-tests, followed by false discovery rate (FDR) correction, show no significant effect of treatment (Supplementary Fig. S2) when comparing treated and untreated injured mice. Consequently, treated and untreated mice are combined and treated as the injured group. Structural damage assessments using MRI scans are conducted by a blinded experimenter. Cortical regions showing structural damage are excluded from further analysis, and PET and fMRI data are preprocessed according to standardized protocols and analyzed as outlined in the Methods sections (“PET/MRI Data Acquisition, Processing, and Analyses”). Statistical comparisons between control and injured groups are conducted using nonparametric testing with 5000 permutations, corrected using family-wise error correction with TFCE, setting statistical significance at p < 0.05, unless otherwise specified. Replication studies involving both sexes are currently being planned, and efforts to secure funding are ongoing.
Footnotes
Acknowledgment
The authors thank the Yale PET Center staff for their professional technical assistance on the PET tracer production and animal PET scans.
Authors’ Contributions
M.M. —Conceptualization, Data curation, Formal analysis, Visualization, Writing—original draft, review and editing. F.M. Conceptualization, Methodology, Data curation, Data acquisition, Writing—review and editing. T.T.—Conceptualization, Methodology, Data acquisition, Formal analysis, Writing—review and editing. J.C.—Supervision, Funding acquisition, Writing—review and editing. A. F.-Z. —Conceptualization, Methodology, Data acquisition, Writing—review and editing, X.S.—Formal analysis, Writing—review and editing. S.M.S.—Conceptualization, Methodology, Formal analysis, Supervision, Funding acquisition, Writing—review and editing. E.L.—Conceptualization, Methodology, Data curation, Formal analysis, Supervision, Funding acquisition, Writing—review and editing.
Author Disclosure Statement
The authors have declared that no conflict of interest exists.
Funding Information
M.M. is funded by SNSF Postdoc.Mobility grant (214392). Z.C. is funded by NIH grants
Supplementary Material
Supplementary Figure S1
Supplementary Figure S2
Supplementary Figure S3
Supplementary Figure S4
Supplementary Figure S5
Supplementary Figure S6
Supplementary Figure S7
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
