Behavioral, blood, and magnetic resonance imaging biomarkers of experimental mild traumatic brain injury

Abstract

Repeated mild traumatic brain injuries (mTBI) may lead to serious neurological consequences, especially if re-injury occurs within the period of increased cerebral vulnerability (ICV) triggered by the initial insult. MRI and blood proteomics might provide objective measures of pathophysiological changes in mTBI, and indicate when the brain is no longer in a state of ICV. This study assessed behavioral, MRI, and blood-based markers in a rat model of mTBI. Rats were given a sham or mild fluid percussion injury (mFPI), and behavioral testing, MRI, and blood collections were conducted up to 30 days post-injury. There were cognitive impairments for three days post-mFPI, before normalizing by day 5 post-injury. In contrast, advanced MRI (i.e., tractography) and blood proteomics (i.e., vascular endothelial growth factor) detected a number of abnormalities, some of which were still present 30 days post-mFPI. These findings suggest that MRI and blood proteomics are sensitive measures of the molecular and subtle structural changes following mTBI. Of particular significance, this study identified novel tractography measures that are able to detect mTBI and may be more sensitive than traditional diffusion-tensor measures. Furthermore, the blood and MRI findings may have important implications in understanding ICV and are translatable to the clinical setting.

Behavioral, blood, and magnetic resonance imaging biomarkers of experimental mild traumatic brain injury

OPEN

R A

P

David K.Wright,, JackTrezise, Alaa Kamnaksh, Ramsey Bekdash, Leigh A. Johnston,, Roger Ordidge, Bridgette D. Semple, Andrew J. Gardner, Peter Stanwell,Terence J. OBrien, DenesV. Agoston & Sandy R. Shultz

Repeated mild traumatic brain injuries (mTBI) may lead to serious neurological consequences, especially if re-injury occurs within the period of increased cerebral vulnerability (ICV) triggered bythe initial insult. MRI and blood proteomics might provide objective measures of pathophysiological changes in mTBI, and indicate when the brain is no longer in a state of ICV. This study assessed days post-injury. There were cognitive impairments for three days post-mFPI, before normalizing by novel tractography measures that are able to detect mTBI and may be more sensitive than traditional in understanding ICV and are translatable to the clinical setting.

A mild traumatic brain injury (mTBI), oen referred to as a concussion, rarely has lasting eects and is oen presumed to cause only transient disturbances to brain function15. However, repeated mTBIs, particularly those occurring in the sports and military settings, have been associated with cumulative and chronic neurological impairments68, and the development of neurodegenerative diseases such as chronic traumatic encephalopathy (CTE)69. There is evidence that these long-term adverse eects of repeated mTBIs are in part due to the recurring insults occurring before the brain has recovered from the initial mTBI and is in a period of increased cerebral vulnerability (ICV)5,10,11. There is increasing evidence that mTBI triggers complex biological changes including inammatory, metabolic, neuronal, vascular and axonal abnormalities14. It is believed that such changes are responsible for ICV and therefore, the identication of reliable markers that indicate when the brain is no longer in a state of ICV might allow them to be used to guide medical decisions.

The current clinical management of mTBI is largely guided by the presence or absence of neuropsychological symptoms, and typically evaluated by subjective and/or self-reported methods24. Symptoms may include physical, cognitive, co-ordination, emotional, and sleep abnormalities24. The onset of symptoms, although typically rapid, can take minutes or hours to occur, and symptoms are usually mild, or may even go unrecognized24,12,13.

Recovery is determined to have occurred aer all post-injury symptoms have resolved, at which point patients

The Florey Institute of Department of Medicine, The Royal Melbourne Department of Anatomy, Physiology, and Department of Centre for

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

are commonly cleared to return to pre-injury activity24,12. However, there is now evidence that the resolution of symptoms might not accurately indicate that the brain has recovered from the neuropathophysiological changes induced by mTBI15,10.

Considering the possible long-term consequences of repetitive mTBIs, and the limitations of current mTBI management approaches, research is required to guide and facilitate more informed medical decisions pertaining to return to pre-injury activity. In particular, it is critical that objective markers sensitive to the brains changes and recovery aer an mTBI are identied. Magnetic resonance imaging (MRI) is a non-invasive and common clinical tool that may be capable of providing objective and quantitative indicators of mTBI pathophysiology to help guide clinical management15,11. Although conventional structural MRI is oen unable to identify pathological changes due to the absence of macroscopic changes1, initial studies applying advanced MRI techniques are gaining prominence as potentially highly sensitive indicators of mTBI. In particular, there is growing evidence that diusion weighted imaging (DWI) and proton magnetic resonance spectroscopy (1H-MRS) may be sensitive to the subtle pathophysiological changes that occur in the mildly injured brain1418. Blood-based biomarkers also hold great promise in the mTBI eld, as there are several candidate protein biomarkers that may be indicative of neuronal and glial cell loss, metabolic abnormalities, vascular changes, neuroinammation, axonal injury, and other pathophysiological mechanisms associated with mTBI11,1928.

Although these initial MRI and blood-based protein ndings demonstrate the potential of these methods to provide insight into brain abnormalities post-mTBI, more detailed studies are required to characterize and validate these methods as well as correlate the microstructural and molecular changes to the commonly used neurobehavioral outcomes. Animal models allow for the control of confounding factors, as well as the rigorous investigation of biomarkers. Previous studies from our laboratory, and others, report that a mild uid percussion injury (mFPI) in the rat induces transient behavioral and pathophysiological changes that occur in the absence of signicant neuronal loss or structural brain damage2931, which is consistent with what may occur aer a single mTBI in humans24. Therefore, here we employed the mFPI model in conjunction with serial multi-modal MRI, blood-based proteomics, and behavioral analyses in order to assess the ability of these methods to detect changes and estimate recovery aer experimental mTBI.

Methods

Subjects. 46 male Long-Evans rats were purchased from Monash animal research services (Melbourne, Australia). All rats were 812 weeks of age, weighed 250300g, and were experimentally nave prior to surgical procedures. Aer surgery, rats were housed individually under a 12h:12h light/dark cycle with ad libitum access to water and food for the duration of the study. All experimental procedures were approved by The University of Melbourne and The Florey Institute of Neuroscience and Mental Health animal ethics committees, and complied with the guidelines of the Australian Code of Practice for the Care and Use of Animals for Scientic Purposes.

All of the surgical procedures used were based on standard protocols as previously described and used by our group29,30. Briey, rats were placed in a sealed Plexiglas box into which 4% Isourane and 2L/min oxygen ow was introduced for anesthesia. Rats were then placed in a stereo-taxic frame via ear bars, with anesthesia maintained at 2% isourane and 1L/min oxygen, and given a subcutaneous injection of analgesic (carprofen 5mg/kg). A craniotomy (5mm diameter) was performed (anterior/posterior 3.0mm, medial/lateral 4.0mm relative to Bregma) to expose the intact dura and a hollow injury cap was affixed to the skull over the craniotomy window. Rats were then attached to the uid percussion device via the headcap. At the rst sign of hind paw withdrawal to a pinch test, rats received a uid percussion pulse of 11.5 atm before being disconnected from the device. Acute injury measures were monitored and included: duration of apnea, dened as the time from injury to the return of spontaneous breathing (sham/mFPI= 0/0s); loss of consciousness, dened as the time to hind limb withdrawal to a pinch test (sham/mFPI=0/1.9 0.9s); and time to self-righting, defined as the time from injury to return of an upright position (sham/mFPI = 82.8 4.0/149.3 9.2 s, p< 0.05). The protocol was identical for sham injuries but no uid percussion pulse was given.

Design. After their assigned sham or mFPI injury, rats were randomly allocated to separate behavior (sham=8; mFPI= 8) or MRI (sham=10; mFPI= 10) cohorts. The behavior cohort underwent behavioral testing on days 1 (D1), 3 (D3), 5 (D5), and 7 (D7) post-injury before being euthanized on D7 for blood collection. The MRI cohort underwent MRI at D1, D3, D5, D7 and 30 days post-injury (D30) before being euthanized on D30 for blood collection. Note that the MRI cohort also underwent behavioral testing on the day prior to their D30 MRI. This study design allowed us to maximize the amount of data collected on corresponding post-injury days, and also avoided the confounding factor of anesthesia (required for MRI) on behavioral outcomes. An additional 10 rats (sham=5; mFPI= 5) were euthanized at D1 for blood analysis. Thus, blood samples were available for D1, D7 and D30 post-injury (5/group/day).

Behavioral Testing. All behavioral testing was conducted by a researcher blinded to experimental conditions. Cognitive function was assessed using the water maze29,30. The water maze apparatus consisted of a circular tank (163 cm in diameter) lled with water (29 1 C), and a hidden acrylic escape platform (10 cm diameter) submerged 2cm below the water surface in one of the quadrants of the pool. Unique visual cues were positioned at the north, south, east and west quadrants. Both the position of the escape platform, and the visual clues, were randomly assigned each day to maintain novelty throughout the course of the experiment. Rats completed four trials per testing day, with the entry position randomized between north, south, east and west for each trial. If the rat failed to nd the platform in the allotted minute, then they were placed on the platform for 15s. Ethovision behavioural tracking soware analyzed video from an overhead camera positioned at the center of the tank. Search time and direct and circle swims (cognition), as well as swim speed (locomotion), were recorded for each trial.

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

Sensorimotor function was assessed using a 1m long, 2 cm wide elevated wooden beam as previously described29,30. Baseline training was performed prior to surgery and consisted of ve successful trials on a 4cm wide beam, followed by ve successful trials on a 2 cm wide beam. On each day of testing, rats completed ten trials on the 2 cm wide beam with the time taken to traverse the beam, as well as the numbers of slips and falls, recorded. The maximum time allowed per trial was 60s, which was also assigned to rats that fell.

As previously described29,30, anxiety-like behavior was assessed using an elevated plus maze. The maze consists of both open and closed arms, 110 cm long and 12 cm wide, perpendicular to each other, and dividing each in half to form a cross. 50cm high walls enclosed the closed arms while the open arms contained no walls. A single 5minute trial was performed on each day of testing with the rat initially positioned at the center of the cross facing an open arm. Ethovision behavioural tracking soware analyzed video from an overhead camera positioned at the center of the cross. The percentage of time spent within the open arm (anxiety), as well as the total number of entries into the closed and open arms (locomotion), were calculated. An entry was only completed once all four paws had entered one of the arms.

Locomotor activity and anxiety-like behavior were assessed using a circular open eld (diameter 100 cm, 40 cm high wall) as previously described29,30. Ethovision behavioural tracking soware (Noldus, Netherlands) analyzed video from an overhead camera positioned at the center of the eld. A single 5 minute trial was performed on each day of testing with the rat positioned in the center of the eld. The total distance travelled (locomotion), as well as the number of entries into and time spent within the center 66 cm diameter area of the eld (anxiety), were recorded.

MRI Acquisition. MRI scanning was performed on D1, D3, D5, D7 and D30 using a 4.7 T Bruker Avance III scanner with 30cm horizontal bore. The magnet was tted with a BGA12S2 actively shielded gradient set capable of 440mT/m and actively decoupled volume transmit and 4-channel surface receive coils were used for imaging. Rats were anesthetized in a clear plastic box with 5% Isourane in a 1:1 mixture of medical grade air and oxygen. Aer loss of consciousness, rats were positioned supinely on a purpose-built rat cradle with stereotactic xation and a nose cone positioned over the rats snout to maintain anesthesia with approximately 2% Isourane in a 1:1 mixture of medical grade air and oxygen. Body temperature was maintained throughout the experiment with a hot water circulation system built into the cradle. Arterial oxygen saturation, heart and breath rates and pulse distention were monitored with a MouseOxTM pulse oximeter (STARR Life Sciences Corp. Pittsburgh PA, USA) and anesthesia adjusted accordingly.

The scanning protocol consisted of a 3-plane localizer sequence followed by multi-slice axial, coronal and sagittal scout images to accurately determine the position of the rat brain. A multi-echo T2*-weighted image was acquired using a 2D gradient echo sequence with the following imaging parameters: repetition time (TR)=4.4s; 8 echoes with the rst echo at 7.5ms and an echo spacing of 7.5 ms; eld of view (FOV)=28.828.8mm2; matrix size=160 160; number of slices= 64; and slice thickness=180m.

Diusion-weighted imaging (DWI) was performed using a 2D echo planar, spin-echo sequence with the following imaging parameters: TR= 6 s; echo time (TE)= 35 ms; FOV=25.625.6mm2; matrix size=128128; number of slices= 24; and slice thickness=600 m. Diusion weighting was performed with diusion duration ()= 3.5ms, diusion gradient separation ()= 14 ms and b-value=1200 s/mm2 in 81 non-collinear directions with 8 non-diusion-weighted (b0) images.

Point resolved spectroscopy (PRESS) was acquired in the ipsilateral cortex with VAPOR water suppression and outer volume saturation. Acquisition parameters were: TR = 2.5 s; TE = 20 ms; number of excitations (NEX)=700; spectral width= 6ppm; number of points=2048; oset=2 ppm; and voxel size=245mm3.

Non-water suppressed spectra were also acquired for eddy current correction and water scaling for metabolite quantication.

MRI Analysis. T2*-weighted Imaging. A study-specic T2*-weighted template image was created by rst generating templates for each cohort at each time point using Advanced Normalization Tools (ANTs, http:// stnava.github.io/ANTs/). ANTs derives an optimal template based on symmetric normalization unbiased by an individual subjects topographical idiosyncrasy32. The resulting templates were then combined, again with ANTs, to create the study-specic template. T2*-weighted images for each rat were then registered to the study-specic template (T2*-to-template) using standard symmetric normalization (SyN), a dieomorphic algorithm included in ANTs32.

Eight a-priori regions of interest (ROIs) were traced using FSLView, a component of FMRIBs Soware Library (FSL, www.fmrib.ox.ac.uk/fsl), on the study-specic T2*-weighted template. Tracing started at the anterior tip of eld CA3 of the hippocampus, approximately 1.72 mm posterior to bregma and continued for 16 slices. ROIs included the ipsilateral and contralateral cortex, hippocampus, corpus callosum and ventricles. The ROIs were registered to subject space using the inverse T2*-to-template dieomorphisms and the total volumes for each structure were calculated using MATLAB (The MathWorks, Natick, MA, USA).

Diusion-Weighted Imaging. Diusion tensor images (DTI) were calculated using DTIFit, part of the FMRIB Diusion Toolbox (FDT, part of FSL) and registered to a study-specic DTI template constructed using DTI-TK (http://dti-tk.sourceforge.net). DTI-TK uses the local bre orientations to guide the alignment of white-matter tracts, which has been shown to result in improved registration over standard fractional anisotropy (FA)-based approaches33. Templates for each cohort at each time point were then used to generate a study template. Each rats DTI was warped directly to the study template (DTI-to-template) and resampled to have a resolution of 0.10.10.3mm3.

The ipsilateral and contralateral cortex, hippocampus and corpus callosum were traced on the study template using FSLview for a-priori ROI-based analyses. Tracing started at approximately 1.8mm bregma and continued

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

for 10 slices. Template ROIs were warped back to individual space using the inverse DTI-to-template dieo-morphisms. For each ROI, the median of FA, radial diusivity (RD), axial diusivity (AD) and trace (TR) were calculated using MATLAB.

Tractography. DWI preprocessing for tractography was performed as described previously34. Briey, spatial intensity inhomogeneity was corrected by estimating the bias eld across the mean b0 image and images were normalized by the contralateral white matter. The bre orientation distribution (FOD) within each voxel was estimated using constrained spherical deconvolution, with a group average single-bre response function, to a maximum harmonic degree of 6.

Whole brain tractograms of 2 million streamlines were generated using the iFOD2 tractography algorithm, part of the MRtrix package (www.brain.org.au/soware)35. Tracking was performed with the following parameters: FOD amplitude for initiating tracks= 0.2; FOD termination threshold= 0.1; FOD samples per step=4; step size=100m; minimum track length=200m; maximum angle between steps=22.5.

Tractograms were registered to the nal DTI template using the DTI-to-template dieomorphisms ensuring that both the length and spatial location of the streamlines were normalized36. The tractograms were then used to generate two tract-weighted images: average pathlength map, representing the average length of tracks passing through each voxel36; and track-weighted curvature map, representing the average curvature of all tracks passing through each voxel (TWI-MC). Tractogram derived average pathlength and TWI-MC images show novel contrast not observed in diusion tensor images such as FA (see Fig.1). As with the tensor measures, the median of each was calculated for the a-priori ROIs.

Magnetic Resonance Spectroscopy. MRS data was processed using LCModel with a metabolite basis set matched to eld strength and TE. The combined absolute concentrations of glutamate + glutamine (Glu+Gln), creatine + phosphocreatine (Cr+PCr), N-acetylaspartate + N-acetylaspartylglutamate (NAA+NAAG), glycerophosphocholine + phosphocholine (GPC+PCh) as well as the single peaks myo-inositol (Ins) and glutathione (GSH), were quantied. Maximum uncertainty in the concentrations (Cramr-Rao lower bounds) was 5% for the combined concentrations and 10% for Ins and GSH.

Plasma Analysis. Blood collection and plasma separation. Following the completion of experiments on D1, D7, or D30 (see Design above), rats were deeply anaesthetized by an intraperitoneal injection of 0.5 ml sodium pentobarbital (Lethabarb, Virbac, Australia). The chest cavity was then opened and blood was collected into a BD Vacutainer K2 EDTA (K2E) Plus Blood collection tube. The tubes were then centrifuged at 6,000g for 15min at room temperature, and plasma was pipetted into 0.5 ml aliquots, ash-frozen in liquid nitrogen, and stored at 80C until use.

Proteomics. Plasma levels of ceruloplasmin, neurolament heavy chain (NF-H), tau protein, vascular endothelial growth factor (VEGF), 4-hydroxynonenal Michael adducts (4-HNE), neuron-specic enolase (NSE), glial brillary acidic protein (GFAP), and S100 calcium binding protein subunit (S100) were assayed using reverse phase protein microarray (RPPM). Sample preparation, printing, scanning, and data analysis for RPPM were performed as described previously1923,3740. Frozen plasma samples were thawed on ice and diluted 1:10 with Dilution Buer (3 parts Lysis Buer [TPER, 10% Glycerol, 1 HALT] and 1 part 4 SDS Sample Buer [35% Glycerol, 0.8% SDS, 10 TBS, 10 TCEP, 1 HALT, 0.0035% NaN3]). Samples were then transferred to a JANUS

Varispan Integrator and Expanded Platform Workstation (PerkinElmer, Waltham, MA) to perform 1:1 serial dilutions using Dilution Buer in 384-well microarray plates (product no. X7022; Molecular Devices, Sunnyvale, CA). An equal amount of 2 PBS Buer (80% Glycerol, 2 PBS) was added to the volume of liquid in each well. The microarray plates were subsequently transferred into an Aushon 2470 Arrayer (Aushon Biosystems, Billerica, MA) where plasma samples were printed onto ONCYTE AVID nitrocellulose lm slides (product no. 305177;

Grace Bio-Labs, Bend, OR). The Aushon 2470 Arrayer was set up with 16 pins and programmed for 2 depositions per spot. The spot diameter was set to 250nm with spacing between dots at 500nm on the x-axis and 375nm on the y-axis. Wash time was set to 2s without delays.

Immunochemical detection. Aer overnight desiccation at 4C, the slides were blocked with a solution of 1 TBS, 5% non-fat dry milk, and 0.1% Tween-20. Slides were then incubated with the primary antibody solutions and a cover slip (product no. 25x 60I-M-5439-001-LS; mSeries LierSlip; Thermo Fisher Scientic, Waltham, MA) overnight at 4C. All primary antibodies were pretested for specicity by Western blot analysis using positive controls and neutralizing peptides when available37. The primary antibodies were diluted in antibody incubation buer (0.1% BSA, EDTA-free Halt Protease and Phosphatase Inhibitor Cocktail [Thermo Fisher Scientic], 1 TBS, and 0.5% Tween-20) and used in the following dilutions: ceruloplasmin (1:20; Santa Cruz Biotechnology, sc-21240), NF-H (1:100; SigmaAldrich, N4142), tau (1:100; Cell Signaling Technology, 4019), VEGF (1:50; Abcam, ab53465), 4-HNE (1:1000; EMD Millipore, 393207), NSE (1:50; Abcam, ab53025), GFAP (1:20; Abcam, ab48050), and S100 (1:20; Abcam, ab41548).

The following day, slides were washed three times with TBST and then incubated with the appropriate secondary antibody solutions for 1h at room temperature. The secondary antibodies Alexa Fluor 790 goat anti-rabbit

(catalog no. A-11369), 680 rabbit anti-goat (catalog no. A-21088), and 680 goat anti-mouse (catalog no. A-21058) (Invitrogen, Eugene, OR) were used at a 1:20,000 dilution. Aer three thorough washes with TBST followed by a single wash with 1 TBS, the slides were air dried and subsequently scanned in an Innopsys InnoScan 710 IR microarray scanner (Innopsys, Carbonne, France).

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

Figure 1. Tract-weighted imaging provides additional contrast to traditional diusion tensor imaging.

Example FA, tractogram, average pathlength and TWI-MC template images for sham rats at D30. Tractogram images were generated with one million streamlines, colour-encoded according to orientation: red, lateral/ medial; green, superior/inferior; blue, anterior/posterior. Average pathlength and TWI-MC images calculated from the tractograms show novel contrast compared to the FA map. TWI-MC=track-weighted mean curvature.

Data acquisition and analysis. Scanner uorescence data were imported into a Microso Excel-based bioinformatics program. After correcting for local background noise, points indiscernible from background were excluded (SNR < 2, Net Fluorescence < 10). Net intensity vs. dilution was plotted on a log-log scale. A ve-parameter logistic (5PL) master curve was tted to the slide. Using this master curve, outliers were excluded using the False Discovery Rate method (Q = 0.01) and the master curve is further rened. Each local block of samples is t individually, using master curve values as the initial conditions for the local curve t. The asymptotic maximum and minimum, Hill-Slope, and asymmetry constant were shared between samples. The slope of the linear portion of the logistic curve was calculated and the line extrapolated back to zero (i.e., the y-intercept), assessing the amount of protein expressed. Data (y-intercept values) are presented as the meanS.E.M.

Statistical analysis. Between group dierences in behavior, MRI and plasma protein levels were analyzed with independent t-tests using SPSS Statistics 22 soware (IBM, New York, NY). Statistical signicance was set at p< 0.05 in all cases.

Results

Rats were tested in the water maze to assess cognitive function. Rats given an mFPI showed signicantly longer search times in the water maze on D3 compared to sham-injured rats (t= 2.215, p= 0.044, Fig.2a). Consistent with these ndings, mFPI rats also performed fewer direct and circle swims compared to sham-injured rats on D3 (t = 2.160, p = 0.047, Fig.2b). Notably, by D5 there were no signicant dierences on either of these measures of cognitive ability between mFPI and sham groups. There were no signicant dierences on the measure of swim speed in the water maze on any of the days, suggesting that motor abnormalities were not a confounding factor in cognitive outcomes. There

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

Figure 2. mFPI results in transient cognitive decits. Rats given an mFPI had longer total search times(a) and fewer direct and circle swims (b) during water maze testing on D3. There were no group dierences on Beam Task (c,d), elevated-plus maze (e,f) or open eld (g,h). Means.e.m., *mFPI group signicantly dierent than sham group, p<0.05.

were no signicant dierences (all p > 0.05) between mFPI and sham rats on any of the beam task (sensorimotor function; Fig.2c,d), elevated-plus maze measures (anxiety-like behavior; Fig.2e,f) or open eld (locomotion, anxiety-like behavior; Fig.2g,h).

mFPI does not result in macroscopic brain damage. To assess for any overt structural brain damage we conducted a volumetric analysis on T2*-weighted images. As shown in Fig.3, mFPI did not result in signicant volumetric changes to any of the eight ROIs, including the ipsilateral and contralateral cortex, corpus callosum, hippocampus, and lateral ventricles (all p>0.05).

mFPI induces transient H-MRS abnormalities. 1H-MRS analyses revealed that mFPI rats had signicantly lower Ins in the ipsilateral cortex on D3 compared to sham-injured rats (t=2.659, p= 0.016; Fig.4a). There were no signicant dierences between mFPI and shams on levels of NAA+NAAG, Glu+Gln, GPC+PCh, GSH, or Cr+PCr (all p>0.05; Fig.4bf).

mFPI induces transient DTI abnormalities. ROI-based DTI analyses revealed that mFPI rats had signicantly decreased FA in the ipsilateral corpus callosum on D3 (t = 2.240, p = 0.038) and D5 (t = 2.620, p = 0.017) when compared to shams (Fig.5b). Additionally, mFPI rats had signicantly increased RD in the ipsilateral corpus callosum on D5 (t = 2.462, p = 0.024, Fig.5c); increased AD in the ipsilateral cortex on D3 (t=2.898, p= 0.009, Fig.5d) and D5 (t=2.354, p= 0.030, Fig.5d); and increased TR in the ipsilateral cortex on D3 (t=2.547 p= 0.020, Fig.5e) compared to sham-injured rats. There were no signicant dierences in any DTI measure in the contralateral hemisphere (Fig.5fi).

mFPI induces persistent tractogram abnormalities. As shown in Fig.6, a-priori ROI-based analyses revealed that mFPI rats had signicantly shorter average pathlength in the ipsilateral corpus callosum on D1 (t = 3.744, p = 0.002), D3 (t = 3.204, p = 0.005) and D5 (t = 2.548, p = 0.020) than their sham-injured counterparts (Fig.6b). Additionally, rats given an mFPI also had signicantly shorter streamlines in the ipsilateral cortex on D1 (t= 2.868, p= 0.010, Fig.6c) and D3 (t = 2.902, p= 0.010, Fig.6c) and in the ipsilateral hippocampus on D3 (t=2.897, p= 0.010, Fig.6d) and D5 (t=2.765, p= 0.013, Fig.6d) compared to shams. No signicant dierences were observed in the contralateral hemisphere at any time point (Fig.6eg).

Rats given an mFPI also had signicantly reduced mean curvature in the ipsilateral corpus callosum on D1 (t = 3.246, p = 0.005, Fig.7b), D3 (t = 2.994, p = 0.008, Fig.7b), D5 (t = 2.491, p = 0.023, Fig.7b) and D30 (t= 2.397, p= 0.028, Fig.7b), and the contralateral corpus callosum on D1 (t= 2.439, p= 0.025), D3 (t= 2.144, p = 0.031) and D30 (t = 2.346, p = 0.046) compared to sham-injured rats (Fig.7e). Additionally, rats given an mFPI also had signicantly reduced mean curvature in the ipsilateral cortex on D1 (t=2.156, p=0.045, Fig.7c) and in the ipsilateral hippocampus on D3 (t=3.525, p= 0.002, Fig.7d) and D5 (t=2.590, p=0.019, Fig.7d).

mFPI alters the levels of select protein biomarkers in the plasma. Rats given an mFPI had significantly elevated levels of ceruloplasmin (t=2.617, p= 0.047, Fig.8a) and NF-H (t=3.435, p= 0.006, Fig.8b) on D1 compared to sham rats. In addition, mFPI rats had decreased levels of tau on D7 (t=2.509, p=0.029, Fig.8c), and decreased VEGF on D7 (t=3.186, p= 0.006, Fig.8d) and D30 (t=2.407, p= 0.037, Fig.8d). There were no statistically signicant dierences in 4-HNE, NSE, GFAP, or S100 plasma levels (all p>0.05; Fig.8eh).

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

Figure 3. mFPI did not result in signicant volumetric changes. (a) T2*-weighted template images for sham (top row) and mFPI (bottom row) groups. There were no signicant dierences observed in any of the eight a-priori ROIs (bi). Means.e.m.

Discussion

This study employed the rat mFPI model in conjunction with advanced MRI, blood proteomics, and behavioral analyses to assess the ability of these dierent measures to detect changes aer experimental mTBI. We report that an experimental mTBI resulted in transient cognitive abnormalities that normalized by D5. Numerous MRI and blood-based outcomes detected changes beyond the resolution of behavioral symptoms, and in some cases these changes were present at D30. Although all three of the biomarker platforms assessed were able to detect abnormalities aer experimental mTBI, our ndings suggest that advanced MRI and blood-based proteomics are more sensitive in identifying the duration of mTBI-induced molecular and ne structural changes than the behavioural measures examined in this study.

Transient cognitive decits were observed in the mFPI group, with no evidence of decits by D5. This is consistent with previously published ndings with the mFPI model2931,41,

as well as clinical concussion where post-concussive symptoms typically resolve within seven days of injury4.

Although no group dierences were observed on D1, this was likely due to the sham animals requiring some exposure to the water maze before they learned how to complete the task. Because there was no group dierence on the measure of swim speed, or on other motor tasks, it is likely that the increased search times and fewer direct and circle swims performed by the mFPI rats were indicative of cognitive impairment and not confounded by motor decits. Although we found no evidence for behavioral impairments beyond D3, it is important to consider whether other behavioral tasks that were not used in this study may have revealed decits at other post-injury time points.

Interestingly, there were a number of abnormalities indicated by MRI that were present at the time of the cognitive decits in mFPI rats. 1H-MRS revealed a signicant decrease in Ins that corresponded with the onset and resolution of cognitive decits. As discussed in greater detail below, changes in Ins may reect an astroglial response, and astrogliosis has previously been associated with cognitive impairments in the days following mFPI29. DTI identied abnormalities in the corpus callosum present at D3 indicative of axonal injury, and previous ndings have also associated DTI-detected corpus callosum injury with cognitive decits aer mTBI14,15.

SCIENTIFIC REPORTS

7

www.nature.com/scientificreports/

Figure 4. mFPI induces transient 1H-MRS abnormalities. mFPI rats had signicantly decreased myo-inositol at D3 compared to sham animals (a). There were no signicant dierences observed in any other metabolites (bf). NAA+NAAG=N-acetylaspartate + N-acetylaspartylglutamate; GPC+PCh=glycerophosphocholine + phosphocholine. Means.e.m., *mFPI group signicantly dierent than sham group, p<0.05.

The biological underpinnings of MRI and blood-based biomarkers. Consistent with neuroimaging results in clinical concussion where conventional structural MRI sequences typically report negative ndings4,

here we found no evidence of overt structural brain damage aer mFPI on conventional T2*-weighted images. However, 1H-MRS and DWI were both able to identify mFPI-induced abnormalities. 1H-MRS revealed a significant decrease in Ins in mFPI rats at D3, however there were no statistically signicant dierences between the groups on 1H-MRS measures from D5 onwards. Ins is a molecule mostly located within astrocytes and is thus widely considered an astroglia marker42. Consistent with our nding here, Ins has previously been found to be decreased in other rat models of TBI43, in the motor cortex and hippocampus aer sports-related concussion18, and in other forms of TBI44. It has been speculated that decreased Ins early aer TBI may be indicative of cell damage and / or death, or be a consequence of an efflux of Ins from astrocytes in order to regulate volume under conditions of edema42. However, it should also be noted that other studies have found increased Ins levels aer TBI45. The exact reasons for these contradictory ndings regarding Ins levels aer TBI are not yet known and require further investigation, though temporal and spatial aspects as well as the severity of the injury are all critical factors to consider.

Signicant abnormalities in the corpus callosum and cortex of mFPI rats were detected by DTI measures on D3, with some of these changes persisting aer cognitive symptoms had normalized. In the ipsilateral corpus callosum there was a signicant decrease in FA on D3 and D5, and an increase in RD on D5, in rats given an mFPI. Whereas in the ipsilatereal cortex there was increased AD on D3 and D5, as well as increased TR on D3.

SCIENTIFIC REPORTS

8

www.nature.com/scientificreports/

Figure 5. mFPI induces transient DTI abnormalities. (a) FA template images for sham (top row) and mFPI (bottom row) rats at each time point. A section of the ipsilateral hemisphere that includes the cortex, corpus callosum and hippocampus is shown enlarged (red box). Qualitatively there appears to be reduced FA in the ipsilateral corpus callosum that persists from D1 to D30. However, ROI-based quantication revealed that mFPI only induced statistically signicant reductions in FA at D3 and D5 in the ipsilateral corpus callosum (b) with a corresponding increase in RD at D5 (c). Additionally, mFPI rats had increased AD at D3 and D5 (d) and increased TR on D3 (e) in the ipsilateral cortex. There were no dierences in any DTI measure in the contralateral hemisphere (fi). Means.e.m., *mFPI group signicantly dierent than sham group, p<0.05.

Traditional DTI measures, such as FA, AD, RD, and TR may reect various pathological events that occur aer brain injury, including axonal injury, edema, gliosis, and blood barrier disruption, as each of these processes are capable of aecting the diusivity of water molecules in the brain46,47. Decreased FA in the corpus callosum has previously been reported aer mTBI15,18,48,49, and this is widely considered a marker for axonal injury within this major white matter tract. The nding of increased AD in the injured cortex aer TBI is consistent with our previous ndings50 and may be related to gliosis46.

Track-weighted imaging methods are relatively new. Derived from properties of the streamlines, such as length, curvature and density, these maps give additional contrasts that may be sensitive to the subtle changes

SCIENTIFIC REPORTS

9

www.nature.com/scientificreports/

Figure 6. mFPI induces transient reductions in average streamline pathlength. (a) Average pathlength template images for sham (top row) and mFPI (bottom row) rats at each time point. Rats given an mFPI had shorter streamlines in the ipsilateral corpus callosum at D1, D3 and D5 (b), in the ipsilateral cortex at D1 and D3 (c), and in the ipsilateral hippocampus at D3 and D5 (d). No dierences were observed contralaterally (eg). Means.e.m., *mFPI group signicantly dierent than sham group, p<0.05.

occurring in the brain aer mTBI. Here, we found that the average pathlength was reduced in the ipsilateral corpus callosum to D5. This is consistent with the original average pathlength article of Pannek and colleagues who reported a signicant decrease in pathlength in the corpus callosum of a TBI patient compared to a cohort of controls36. Additionally, we also observed reduced average pathlength in the ipsilateral hippocampus to D5 and in the ipsilateral cortex to D3. These changes likely reect underlying DTI changes within these regions, which would be expected to cause the premature termination of streamlines. While the average pathlength did not appear to be any more sensitive than traditional DTI measures, we found changes in streamline curvature at D30 in both the

SCIENTIFIC REPORTS

10

www.nature.com/scientificreports/

Figure 7. mFPI induces persistent reductions in average streamline curvature. (a) TWI-MC template images for sham (top row) and mFPI (bottom row) rats at each time point. Rats given an mFPI had signicantly reduced curvature in the ipsilateral (b) and contralateral (e) corpus callosum that persisted to D30 post-injury. Additionally, mFPI rats had signicantly reduced mean curvature in the ipsilateral cortex at D1 (c) and the ipsilateral hippocampus on D3 and D5 (d). Means.e.m., *mFPI group signicantly dierent than sham group, p<0.05.

ipsilateral and contralateral corpus callosum, suggesting that streamline curvature may be more sensitive to the subtle changes occurring in the brain aer mFPI. Reduced streamline curvature may also reect the early termination of streamlines; however, further research is required to determine how changes in streamline properties, such as length and curvature, relate to the underlying pathology. Additionally, as we ceased this study at D30 it remains to be seen whether these changes are permanent, or whether there is still the potential for full recovery.

SCIENTIFIC REPORTS

11

www.nature.com/scientificreports/

Figure 8. mFPI alters protein biomarker levels in plasma. Rats given an mFPI had signicantly increased levels of ceruloplasmin (a) and NF-H (b) in plasma on D1 compared to sham-injured rats. Additionally, mFPI rats had signicantly decreased levels of tau on D7 (c), and VEGF on D7 and D30 (d) compared to their sham counterparts. There were no statistically signicant dierences between the groups on any of the other plasma biomarkers (eh). Means.e.m., *mFPI group signicantly dierent than sham group, p<0.05.

RPPM analysis of blood-based protein biomarkers revealed that mFPI resulted in increased levels of ceruloplasmin and NF-H at D1, decreased levels of tau at D7, and decreased levels of VEGF at D7 and D30. Ceruloplasmin is considered a marker for metabolic abnormalities post-TBI19, and our nding of elevated circulating ceruloplasmin at D1 is consistent with previous studies that also report increased ceruloplasmin in the early stages post-TBI26. Ceruloplasmin, a positive acute-phase protein, is involved in exporting free iron and in iron metabolism in the brain26,51. TBI causes elevated levels of iron, which can trigger oxidative stress and neuroinammation26,51. Thus, the increased ceruloplasmin levels observed in this study may be in response to an mFPI-induced acute increase in free iron.

Increased levels of circulating NF-H have been reported aer TBI20,52,53, including mTBI26, and is considered a marker of axonal injury. Consistent with previous ndings, the current study also found elevated levels of NF-H in the plasma of mFPI rats at D1. Combined, the DWI and NF-H ndings support the notion that mFPI induced axonal injury. Though it is interesting to note that the plasma changes in NF-H manifested earlier (i.e., D1) but were not present at the later time points, whereas DWI changes emerged later (i.e., D3) and were present at D30.

Tau was decreased in the plasma of mFPI rats at D7, which is consistent with previous studies from our laboratory indicating decreased levels of total tau in brain tissue following TBI54. However, our plasma ndings contradict those from other studies that have reported increased levels of circulating tau, which is postulated to be a marker of axonal injury19, aer mTBI25,26. The reasons for these discrepancies should be further explored, and the mechanism/pattern of injury (i.e., blast versus focal versus diuse) may be an important factor to consider as this could aect the consequent pathophysiological cascade. Unfortunately, the absence of a specic phospho-tau antibody suitable for RPPM at the time of the study prevented us from determining the eect of the injury on the ratio between total tau and phosphorylated tau. In our earlier studies, we have found increased levels of phosphorylated tau, an important cofactor in neurodegenerative diseases such as Alzheimers and CTE, in brain tissue following TBI54,55. It is possible that the reduced circulating tau in this study may be related to an increase in phosphorylated-tau. Because tau abnormalities are implicated in the cumulative and long-term eects of mTBIs, it is important that further studies are conducted to gain a more comprehensive understanding of the nature of these changes in the brain, how they relate to other pathology (e.g., axonal injury, neuroinammation), and how blood-protein levels reect these changes.

VEGF is known to be a potent stimulator of angiogenesis (i.e., the formation of new blood vessels), but VEGF also has direct eects on other types of neural cells56. In fact, inhibition of VEGF signaling results in increased cell death aer TBI57. Conversely, treatment with VEGF has been found to stimulate neurogenesis and is neuro-protective aer experimental TBI in rats58,59. In this study, VEGF was decreased at D7 and D30 in the plasma of rats given an mFPI. These ndings may have important implications, particularly with regards to susceptibility to repeated mTBI and neurodegeneration. For example, angiogenesis post-TBI may be an important recovery process that restores adequate blood and oxygen supply to damaged tissue60. If this process has been compromised in rats given an mFPI, the injured tissue may be at increased vulnerability if exposed to a consequent insult. Notably, reduced VEGF levels have been shown to cause neurodegeneration by impairing neural tissue perfusion56, and mice with reduced VEGF levels develop motor neuron degeneration reminiscent of amyotrophic lateral sclerosis61, a neurodegenerative disease that has been epidemiologically linked to brain trauma8. Taken together, this novel nding of decreased VEGF in the sub-acute to long-term stages aer a single mTBI may have important implications in understanding the underlying mechanisms involved in ICV, the degenerative eects of repeated mTBI, as well as potential treatment strategies, and should be further investigated in future studies.

Although other studies have reported that mTBI induces changes in S100, GFAP, 4-HNE, and NSE2426,62,63,

we did not nd any statistically signicant dierences pertaining to these markers in our study. Variable ndings between laboratories/studies with regards to specic blood-based protein biomarkers is not uncommon and may be related to methodological differences24. Our findings also indicate temporal complexities with the various markers, which is consistent with previous ndings suggesting that a single point assay of a single

SCIENTIFIC REPORTS

12

www.nature.com/scientificreports/

protein biomarker is not able to capture the entire pathophysiological cascade of mTBI23,26. However, when taken together, RPPM analysis of plasma proteins did reveal mFPI-induced changes that emerged on D1 (i.e., before signicant behavioral or MRI-related changes), and were present at D7 and D30 (i.e., aer behavioral symptoms normalized).

Future direction and hurdles. While the use of symptom-based tools can aid in mTBI diagnosis, their subjective nature, limited sensitivity, need for accurate baseline measures, practice eects, and lack of any patho-physiology quantication are problematic4,6466. It is for these reasons that objective biomarker platforms capable of assessing biological changes are being pursued to either supplement or replace symptom-based management approaches. This study provides evidence that DTI and tractography methods detect changes for up to 30 days post-injury, long aer behavioral changes had normalized, and should continue to be developed and employed in future investigations of mTBI. In particular, clinical studies with serial neuroimaging aer mTBI are required to better characterize the temporal progression of injury and recovery, and how these changes relate to symptoms and neuropsychological alterations. Unfortunately, due to the high cost, need for specialized equipment/personnel, and the additional research required validating their use as reliable and eective biomarkers, it is unlikely that neuroimaging will be widely accessible as a clinical mTBI management tool for routine clinical use in the near future.

Alternatively, blood-based biomarkers may represent a more readily available strategy to provide information about ongoing pathophysiological changes aer mTBI that aids in diagnosis and monitoring injury progression/ recovery24. In this study tau and VEGF detected changes that were present aer behavioral recovery, and there are other protein biomarkers that may be sensitive to the neuronal and glial cell loss, metabolic abnormalities, neuroinammation, axonal injury, and other pathophysiological changes associated with mTBI19,24. In addition, microRNAs represent another source of promising blood-based biomarkers for mTBI that should continue to be explored67. However, to our knowledge there are no detailed longitudinal studies that have assessed blood bio-markers in mTBI patients in order to validate their use as diagnostic and prognostic tools24. Furthermore, the lack of standardized methods to quantify blood biomarkers, and the need for baseline measures for each individual and/or validated cut-o points are additional limitations to be addressed before the use of blood-based biomarkers in mTBI management24.

While MRI and blood-based biomarker platforms may have the capacity to help guide medical decisions, to date, there is no direct evidence that the implementation of any of these methods will actually prevent or reduce the long-term eects of repeated mTBIs. Because the studies required to assess this would take decades to complete in humans, and would be confounded by a number of factors, animal model studies that administer repeated mTBIs with the inter-injury times guided by a specic biomarker would be useful in providing evidence that management based on particular biomarkers is eective.

There are some limitations with this study that should also be considered. The mFPI requires a craniotomy and the administration of anesthetic, both of which are not typically present in the clinical mTBI setting and have the potential to confound results6870. Nonetheless, previous studies have found that a single mFPI can result in transient motor and cognitive decits29,31,41,71, as well as abnormalities in other behaviors29, sleep72, and electrophysiology73,74. An mFPI also induces transient neuroinammation29,30, axonal injury29, the hyperphosphorylation of tau55, and reduced cerebral blood ow73, but does not result in signicant neuronal loss, visible brain contusion, or focal lesion2931,71,73. These changes are consistent with those that might occur in clinical mTBI. Taken together, the use of the mFPI to model mTBI in this and future studies is warranted, particularly while newer models of mTBI/concussion that avoid limitations such as craniotomy and anesthetic are developed, characterized, validated, and become widely available to other laboratories. Another limitation is the lack of histological/immunohistochemical analyses of brain tissue that may have complemented the other outcomes. For example, additional information pertaining to the nature of the brain damage and other neuropathological changes observed in rats given an mFPI (e.g., conrming DTI ndings of axonal injury in the corpus callosum; relationship between Ins and neuroinammation; localization of blood biomarker abnormalities such as VEGF) would have been informative and should be examined in future studies. Though it is worth noting that previous studies of neuropathology aer an mFPI in the rat have found that it causes increased silver staining and amyloid precursor protein immunoreactivity (i.e., markers of axonal injury) in the corpus callosum and other white matter tracts, as well as neuroinammation in the frontal and parietal cortex, corpus callosum, and thalamus29,31, all of which may be of relevance to the current ndings.

Conclusion

This study investigated the use of MRI, blood proteomics, and behavioral methods as markers to detect changes and estimate recovery aer experimental mTBI. All three of the platforms were able to detect abnormalities aer mFPI; however, advanced MRI and blood-based proteomics proved to be more sensitive for the detection of persisting and/or long-term changes than the behavioral measures examined. In particular, our ndings related to blood biomarkers, such as VEGF and tau, might have important implications in the context of ICV and the cumulative/chronic consequences of repeated mTBI. Furthermore, to our knowledge this is the rst study to report that advanced track-weighted imaging measures, such as mean curvature and average pathlength, are able to detect mTBI and may be more sensitive to persisting white matter abnormalities than traditional DTI measures (e.g., FA). It is important to note that further studies are required to determine the biological signicance/ underpinnings of these markers and their relationship to neurobehavioral changes, and there are a number of challenging hurdles before these markers can be reliably applied in the clinical mTBI setting. Nonetheless, the characterization and utilization of these and other objective biomarkers should continue to be pursued in both the laboratory and the clinic due to their potential value in guiding management decisions, including return to play, following concussion.

SCIENTIFIC REPORTS

13

www.nature.com/scientificreports/

References

1. Barkhoudarian, G., Hovda, D. A. & Giza, C. C. The molecular pathophysiology of concussive brain injury. Clin. Sports Med. 30, 3348 viiiii (2011).

2. Khurana, V. G. & Kaye, A. H. An overview of concussion in sport. J. Clin. Neurosci. 19, 111 (2012).3. Jordan, B. D. The clinical spectrum of sport-related traumatic brain injury. Nat. Rev. Neurol. 9, 222230 (2013).4. McCrory, P. et al. Consensus statement on concussion in sport-The 4th international conference on concussion in sport held in Zurich, November 2012. J. Sci. Med. Sport 16, 178189 (2013).

5. Prins, M. L., Alexander, D., Giza, C. C. & Hovda, D. A. Repeated mild traumatic brain injury: mechanisms of cerebral vulnerability. J. Neurotrauma 30, 3038 (2013).6. McKee, A. C. et al. Chronic traumatic encephalopathy in athletes: progressive tauopathy aer repetitive head injury. J. Neuropathol. Exp. Neurol. 68, 709735 (2009).

7. McKee, A. C. et al. The spectrum of disease in chronic traumatic encephalopathy. Brain 136, 4364 (2013).8. McKee, A. C. et al. TDP-43 proteinopathy and motor neuron disease in chronic traumatic encephalopathy. J. Neuropathol. Exp. Neurol. 69, 918929 (2010).

9. Blennow, K., Hardy, J. & Zetterberg, H. The neuropathology and neurobiology of traumatic brain injury. Neuron 76, 886899 (2012).

10. Povlishock, J. T. The window of risk in repeated head injury. J. Neurotrauma 30, 1 (2013).11. Jeter, C. B. et al. Biomarkers for the diagnosis and prognosis of mild traumatic brain injury/concussion. J. Neurotrauma 30, 657670 (2013).

12. Moser, R. S. et al. Neuropsychological evaluation in the diagnosis and management of sports-related concussion. Arch. Clin. Neuropsychol. 22, 909916 (2007).

13. Vagnozzi, R. et al. Assessment of metabolic brain damage and recovery following mild traumatic brain injury: a multicentre, proton magnetic resonance spectroscopic study in concussed patients. Brain 133, 32323242 (2010).

14. Cubon, V. A., Putukian, M., Boyer, C. & Dettwiler, A. A diusion tensor imaging study on the white matter skeleton in individuals with sports-related concussion. J. Neurotrauma 28, 189201 (2011).

15. Niogi, S. N. et al. Extent of microstructural white matter injury in postconcussive syndrome correlates with impaired cognitive reaction time: a 3T diusion tensor imaging study of mild traumatic brain injury. AJNR Am. J. Neuroradiol. 29, 967973 (2008).16. Henry, L. C. et al. Acute and chronic changes in diusivity measures aer sports concussion. J. Neurotrauma 28, 20492059 (2011).

17. Gasparovic, C. et al. Neurometabolite concentrations in gray and white matter in mild traumatic brain injury: an 1H-magnetic resonance spectroscopy study. J. Neurotrauma 26, 16351643 (2009).

18. Chamard, E. et al. Neurometabolic and microstructural alterations following a sports-related concussion in female athletes. Brain Inj. 27, 10381046 (2013).

19. Agoston, D. V., Gyorgy, A., Eidelman, O. & Pollard, H. B. Proteomic biomarkers for blast neurotrauma: targeting cerebral edema, inammation, and neuronal death cascades. J. Neurotrauma 26, 901911 (2009).

20. Gyorgy, A. et al. Time-dependent changes in serum biomarker levels after blast traumatic brain injury. J. Neurotrauma 28, 11211126 (2011).

21. Agoston, D. V. & Elsayed, M. Serum-based protein biomarkers in blast-induced traumatic brain injury spectrum disorder. Front. Neurol. 3, 107 (2012).

22. Ahmed, F. et al. Time-dependent changes of protein biomarker levels in the cerebrospinal uid aer blast traumatic brain injury. Electrophoresis 33, 37053711 (2012).

23. Ahmed, F. A., Kamnaksh, A., Kovesdi, E., Long, J. B. & Agoston, D. V. Long-term consequences of single and multiple mild blast exposure on select physiological parameters and blood-based biomarkers. Electrophoresis 34, 22292233 (2013).

24. Zetterberg, H., Smith, D. H. & Blennow, K. Biomarkers of mild traumatic brain injury in cerebrospinal uid and blood. Nat. Rev. Neurol. 9, 201210 (2013).

25. Shahim, P. et al. Blood biomarkers for brain injury in concussed professional ice hockey players. JAMA Neurol. 71, 684692 (2014).

26. Ahmed, F., Plantman, S. & Cernak, I. The temporal pattern of changes in serum biomarker levels reveals complex and dynamically changing pathologies aer exposure to a single low-intensity blast in mice. Front. Neurol. 6, 114 (2015).

27. Siman, R. et al. Serum SNTF increases in concussed professional ice hockey players and relates to the severity of postconcussion symptoms. J. Neurotrauma 32, 12941300 (2015).

28. Shahim, P. et al. Serum tau fragments predict return to play in concussed professional ice hockey players. J. Neurotrauma, doi: 10.1089/neu.2014.3741 (2015).

29. Shultz, S. R., MacFabe, D. F., Foley, K. A., Taylor, R. & Cain, D. P. A single mild uid percussion injury induces short-term behavioral and neuropathological changes in the Long-Evans rat: support for an animal model of concussion. Behav. Brain Res. 224, 326335 (2011).

30. Shultz, S. R. et al. Repeated mild lateral uid percussion brain injury in the rat causes cumulative long-term behavioral impairments, neuroinammation, and cortical loss in an animal model of repeated concussion. J. Neurotrauma 29, 281294 (2012).

31. Hylin, M. J. et al. Behavioral and histopathological alterations resulting from mild uid percussion injury. J. Neurotrauma 30, 702715 (2013).

32. Avants, B. B. et al. A reproducible evaluation of ANTs similarity metric performance in brain image registration. Neuroimage 54, 20332044 (2011).

33. Keihaninejad, S. et al. An unbiased longitudinal analysis framework for tracking white matter changes using diusion tensor imaging with application to Alzheimers disease. Neuroimage 72, 153163 (2013).

34. Willats, L. et al. Quantication of track-weighted imaging (TWI): characterisation of within-subject reproducibility and between-subject variability. Neuroimage 87, 1831 (2014).

35. Tournier, J.-D., Calamante, F. & Connelly, A. Robust determination of the bre orientation distribution in diusion MRI: nonnegativity constrained super-resolved spherical deconvolution. Neuroimage 35, 14591472 (2007).

36. Pannek, K. et al. The average pathlength map: a diusion MRI tractography-derived index for studying brain pathology. Neuroimage 55, 133141 (2011).

37. Gyorgy, A. et al. Reverse phase protein microarray technology in traumatic brain injury. J. Neurosci. Methods 192, 96101 (2010).38. Kovesdi, E. et al. Acute minocycline treatment mitigates the symptoms of mild blast-induced traumatic brain injury. Front. Neurol. 3, 111 (2012).

39. Rostami, E. et al. Time-dependent changes in serum level of protein biomarkers after focal traumatic brain injury. Int. J. Neurorehabil. 2, 26 (2015).

40. Kwon, S. K. et al. Stress and traumatic brain injury: a behavioral, proteomics, and histological study. Front. Neurol. 2, 12 (2011).41. DeRoss, A. L. et al. Multiple head injuries in rats: eects on behavior. J. Trauma 52, 708714 (2002).42. Harris, J. L., Choi, I.-Y. & Brooks, W. M. Probing astrocyte metabolism in vivo: proton magnetic resonance spectroscopy in the injured and aging brain. Front. Aging Neurosci. 7, 202 (2015).

43. Xu, S. et al. Early microstructural and metabolic changes following controlled cortical impact injury in rat: a magnetic resonance imaging and spectroscopy study. J. Neurotrauma 28, 20912102 (2011).

SCIENTIFIC REPORTS

14

www.nature.com/scientificreports/

44. Harris, J. L. et al. Altered neurochemical prole aer traumatic brain injury: (1)H-MRS biomarkers of pathological mechanisms. J. Cereb. Blood Flow Metab. 32, 21222134 (2012).45. Ashwal, S. et al. Proton spectroscopy detected myoinositol in children with traumatic brain injury. Pediatr. Res. 56, 630638 (2004).46. Budde, M. D., Janes, L., Gold, E., Turtzo, L. C. & Frank, J. A. The contribution of gliosis to diffusion tensor anisotropy and tractography following traumatic brain injury: validation in the rat using Fourier analysis of stained tissue sections. Brain 134, 22482260 (2011).

47. Shultz, S. R., OBrien, T. J., Stefanidou, M. & Kuzniecky, R. I. Neuroimaging the epileptogenic process. Neurotherapeutics 11, 347357 (2014).

48. Arfanakis, K. et al. Diusion tensor MR imaging in diuse axonal injury. AJNR Am J Neuroradiol 23, 794802 (2002).49. Inglese, M. et al. Diuse axonal injury in mild traumatic brain injury: a diusion tensor imaging study. J. Neurosurg. 103, 298303 (2005).

50. Johnstone, V. P. A. et al. Experimental traumatic brain injury results in long-term recovery of functional responsiveness in sensory cortex but persisting structural changes and sensorimotor, cognitive, and emotional decits. J. Neurotrauma 32, 13331346 (2015).

51. Ayton, S. et al. Ceruloplasmin and -amyloid precursor protein confer neuroprotection in traumatic brain injury and lower neuronal iron. Free Radic. Biol. Med. 69, 331337 (2014).

52. Anderson, K. J. et al. The phosphorylated axonal form of the neurolament subunit NF-H (pNF-H) as a blood biomarker of traumatic brain injury. J. Neurotrauma 25, 10791085 (2008).

53. urek, J., Bartlov, L. & Fedora, M. Hyperphosphorylated neurolament NF-H as a predictor of mortality aer brain injury in children. Brain Inj. 25, 221226 (2011).

54. Shultz, S. R. et al. Sodium selenate reduces hyperphosphorylated tau and improves outcomes aer traumatic brain injury. Brain 138, 12971313 (2015).

55. Tan, X.L., et al. Sodium selenate, a protein phosphatase 2A activator, mitigates hyperphosphorylated tau and improves repeated mild traumatic brain injury outcomes. Neuropharmacology, doi: 10.1016/ j.neuropharm.2016.05.001. (2016)

56. Storkebaum, E. & Carmeliet, P. VEGF: a critical player in neurodegeneration. J. Clin. Invest. 113, 1418 (2004).57. Lee, C. & Agoston, D. V. Inhibition of VEGF receptor 2 increased cell death of dentate hilar neurons aer traumatic brain injury. Exp. Neurol. 220, 400403 (2009).

58. Thau-Zuchman, O., Shohami, E., Alexandrovich, A. G. & Leker, R. R. Subacute treatment with vascular endothelial growth factor aer traumatic brain injury increases angiogenesis and gliogenesis. Neuroscience 202, 334341 (2012).

59. Thau-Zuchman, O., Shohami, E., Alexandrovich, A. G. & Leker, R. R. Vascular endothelial growth factor increases neurogenesis aer traumatic brain injury. J. Cereb. Blood Flow Metab. 30, 10081016 (2010).

60. Xiong, Y., Mahmood, A. & Chopp, M. Angiogenesis, neurogenesis and brain recovery of function following injury. Curr. Opin. Investig. Drugs 11, 298308 (2010).

61. Oosthuyse, B. et al. Deletion of the hypoxia-response element in the vascular endothelial growth factor promoter causes motor neuron degeneration. Nat. Genet. 28, 131138 (2001).

62. Mannix, R., Eisenberg, M., Berry, M., Meehan, W. P. & Hayes, R. L. Serum biomarkers predict acute symptom burden in children aer concussion: a preliminary study. J. Neurotrauma 31, 10721075 (2014).

63. Papa, L. et al. GFAP out-performs S100 in detecting traumatic intracranial lesions on computed tomography in trauma patients with mild traumatic brain injury and those with extracranial lesions. J. Neurotrauma 31, 18151822 (2014).

64. Comper, P., Hutchison, M., Magrys, S., Mainwaring, L. & Richards, D. Evaluating the methodological quality of sports neuropsychology concussion research: a systematic review. Brain Inj. 24, 12571271 (2010).

65. Echemendia, R. J. et al. Advances in neuropsychological assessment of sport-related concussion. Br. J. Sports Med. 47, 294298 (2013).

66. Randolph, C., McCrea, M. & Barr, W. B. Is neuropsychological testing useful in the management of sport-related concussion? J. Athl. Train. 40, 139154 (2005).

67. Redell, J. B., Moore, A. N., Ward, N. H., Hergenroeder, G. W. & Dash, P. K. Human traumatic brain injury alters plasma microRNA levels. J. Neurotrauma 27, 21472156 (2010).

68. Statler, K. D. et al. Comparison of seven anesthetic agents on outcome aer experimental traumatic brain injury in adult, male rats. J. Neurotrauma 23, 97108 (2006).69. Cole, J. T. et al. Craniotomy: true sham for traumatic brain injury, or a sham of a sham? J. Neurotrauma 28, 359369 (2011).70. Plog, B. A. et al. Biomarkers of traumatic injury are transported from brain to blood via the glymphatic system. J. Neurosci. 35, 518526 (2015).

71. Gurko, G. G., Giza, C. C. & Hovda, D. A. Lateral uid percussion injury in the developing rat causes an acute, mild behavioral dysfunction in the absence of signicant cell death. Brain Res. 1077, 2436 (2006).

72. Aungst, S. L., Kabadi, S. V., Thompson, S. M., Stoica, B. A. & Faden, A. I. Repeated mild traumatic brain injury causes chronic neuroinammation, changes in hippocampal synaptic plasticity, and associated cognitive decits. J. Cereb. Blood Flow Metab. 34, 12231232 (2014).

73. Lim, M. M. et al. Dietary therapy mitigates persistent wake decits caused by mild traumatic brain injury. Sci. Transl. Med. 5, 215ra173 (2013).

74. Johnstone, V. P. A., Shultz, S. R., Yan, E. B., OBrien, T. J. & Rajan, R. The acute phase of mild traumatic brain injury is characterized by a distance-dependent neuronal hypoactivity. J. Neurotrauma 31, 18811895 (2014).

Acknowledgements

The authors thank Mr. John Olson for his assistance with the blood biomarker assays and acknowledge the facilities, and the scientific and technical assistance of the National Imaging Facility at the Florey Institute of Neuroscience and Mental Health, Parkville. Research reported in this publication was supported by the Australian National Health and Medical Research Council (grant numbers 1062653, 1087172) and the Center for Neuroscience and Regenerative Medicine (grant number G170AZ).

Author Contributions

All authors contributed to the design of the study, the interpretation of the results, and writing the manuscript. D.K.W. carried out the acquisition and analysis of MRI data. J.T. performed behavioural tests and tissue collection. A.K., R.B. and D.V.A. designed and performed the blood biomarker assays. S.R.S. conceived the study and performed all surgeries.

Additional Information

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Wright, D. K. et al. Behavioral, blood, and magnetic resonance imaging biomarkers of experimental mild traumatic brain injury. Sci. Rep. 6, 28713; doi: 10.1038/srep28713 (2016).

SCIENTIFIC REPORTS

15

www.nature.com/scientificreports/

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

16