2.1.

Surgical Procedures

All animal experiments were performed in accordance with the guidelines of the White Oak Institutional Animal Care and Use Committee. Optical and electrophysiological measurements were performed on 10 male C57BL/6J mice (12 to 24 weeks). Anesthesia was induced by an initial exposure to 4% isoflurane (vaporized in medical grade compressed oxygen) for $<30\mathrm{s}$. Animals were additionally administered an injection of xylazine ($18\mathrm{mg}/\mathrm{kg}$, IP) to provide a stable plane of anesthesia at low isoflurane concentrations (0.1% to 0.25%) for the remainder of the experimental session, which typically lasted 3 to 6 h. Maintenance doses of xylazine ($6\mathrm{mg}/\mathrm{kg}$) were administered once every $\sim 2.5\mathrm{h}$. Following initial anesthesia induction, animals were positioned in a stereotaxic apparatus (David Kopf Instruments, California). Their body temperature was measured and maintained at 37°C with a closed-loop temperature-controlled heating pad (Model TC-1000, CWE). Respiratory rate was also monitored and maintained at $\sim 100\text{breaths}/\min$ during the surgical and experimental procedures. Skin incisions were infused with lidocaine, and the eyes were covered with ointment (Lacri-Lube) to prevent drying. A midline sagittal incision was made in the skin, which exposed the coronal and lambdoid sutures on the skull. The intersection of these sutures with the midline (i.e., bregma and lambda) served as landmarks for recording locations and was also used as a guide when drilling burr holes. The area of the skull under the probe was cleaned with 70% ethanol, and the optical probe was secured to the skull by cyanoacrylate glue (Loctite 454, Hankel, Australia) [Fig. 1(a)]. Burr holes (diameter: $\sim 0.5\mathrm{mm}$) were made for placement of two silver wire electrodes: a recording electrode, placed 2.5-mm lateral and 1-mm posterior to bregma, and a reference electrode, placed 1-mm lateral and 1-mm posterior to lambda. The recording electrode was embedded within the optical probe tip, and the reference electrode was secured into position with Kwik-Sil adhesive (World Precision Instruments, Florida). A ground needle electrode was placed subcutaneously on the back of the animal.

Fig. 1

Experimental design and measurement apparatus. (a) Schematic of the recording and stimulation configuration. A DCS probe tip containing optical fibers is placed on the region of the skull above the primary somatosensory cortex, at a location coinciding with the forepaw’s representation. A recording EEG electrode is inserted through the center of the probe. Pulses of electrical current are delivered to the contralateral limb with respect to the optical and electrophysiological measurements. Closed-skull impact is delivered by a 5-mm impactor tip $\sim 4\text{-}\mathrm{mm}$ anterior and $\sim 1\text{-}\mathrm{mm}$ lateral of bregma contralateral to the optical measurement (see Sec. 2). (b) Overview of synchronization between stimulus, measurements of CBF, and SSEPs. The red traces depict simulated hemodynamic signals, which are broad compared with electrophysiological potentials (note the time scale difference in the insets). The blue pulses indicate electrical stimuli delivered to the median nerve. (c) Block diagram of the DCS device used in this study. (d) Photograph and schematic of the 3-D-printed DCS probe tip. Light is delivered and collected from the tissue through fiber optics that are coupled to the skin with microprisms.

Although epidural measurements of SSEPs are invasive, we have previously compared the effects of injury on SSEPs measured both epidurally and epidermally.^23,24 In mouse experiments, SSEP waveforms measured using both approaches differ only in signal-to-noise ratio, which is roughly an order of magnitude higher when measured epidurally. In this study, the use of epidural recordings permitted us to rapidly acquire clear SSEPs at a high sampling rate following injury owing to the fact that less time was needed to be spent averaging, compared with epidermal recordings. This enhanced the temporal resolution with which injury-induced changes could be tracked.

2.2.

Sensory Stimuli and Electrophysiological Recordings

SSEPs were recorded in single-ended configuration (RZ5D processor, PZ2 preamplifier, ZC16 headstage, Tucker-Davis Technologies, Florida) with a shared common reference electrode, at a sampling frequency of 3 kHz. Optical measurements of $\mathrm{\Delta }\mathrm{CBF}$ were performed concurrently, driven by a separate computer. The median nerve contralateral to the recording locations was stimulated via a pair of 27-gauge stainless-steel needles inserted subcutaneously into the forelimb of the animal. The electrical stimulus consisted of a train of 12 current pulses generated by a constant current stimulator (DS7A, Digitimer) (amplitude 4 mA, frequency 3 Hz, pulse duration $200\mu \mathrm{s}$, and total pulse-train duration 4 s). Individual measurements of sensory-evoked $\mathrm{\Delta }\mathrm{CBF}$ were separated by 45 s, which we empirically found to be the minimum intertrial duration that did not elicit alterations in the steady-state CBF. An overview of the acquisition and stimulus timing is shown in Fig. 1(b). The major peaks associated with the SSEP waveform were identified with MATLAB^® as maxima and minima within boundaries set using published values as a reference ^23,25,26 as follows: P1 corresponded to the time of maximum voltage within the window of 15 to 30 ms, N1 corresponded to the time of minimum voltage within the window of 21 to 50 ms, and P2 corresponded to the time of maximum voltage at times later than the identified N1 time.

2.3.

Controlled Cortical Impact

CCI is a well-established and highly reproducible brain injury model.^{20,21,27–30} Closed-skull impact was delivered with an Impact One™ stereotaxic impactor (Leica Microsystems). Using a 5-mm-diameter metal impact tip, we used the following settings: velocity $4\mathrm{m}/\mathrm{s}$, dwell time 100 ms, and impact depth 0.6 mm. Strike velocities $>5\mathrm{m}/\mathrm{s}$ and/or $\text{strike depths}>0.8\mathrm{mm}$ resulted in skull fracture. Additionally, impactor tips with smaller diameter tended to cause fractures. The metal impactor tip was positioned over the exposed skull $\sim 4\text{-}\mathrm{mm}$ anterior and $\sim 1\text{-}\mathrm{mm}$ lateral of bregma, contralateral to the optical measurement.

2.4.

Optical Measurement of Cerebral Blood Flow

The DCS signal can be used to measure blood flow by using speckle correlation techniques. The basic theory underlying DCS has been extensively described in previous publications.^16,31 Briefly, when laser light migrates through tissue, the emerging intensity pattern, called a speckle pattern, is composed of bright and dark spots, which are caused by constructive and destructive interference of photons that traverse different path lengths. The intensity fluctuations of a single region, which are caused by interactions of scattered light with moving particles (i.e., red blood cells), can be utilized to extract information about blood flow.

A block diagram of our DCS recording system is shown in Fig. 1(c). It consists of a long coherence length continuous-wave near-infrared laser (785 nm, CrystaLaser, Nevada), a photon-counting avalanche photodiode (APD) (SPCM-AQRH-12-FC, Excelitas, Quebec, Canada), and an autocorrelator signal processing board (Ref. 32, New Jersey). NIR excitation light was delivered to the brain with a multimode optical fiber ($200\text{-}\mu \mathrm{m}$ core diameter, Thorlabs, New Jersey), and the scattered light was detected with a single-mode fiber ($5\text{-}\mu \mathrm{m}$ core diameter) connected to the APD. The source and detection fibers were separated by 5 mm. The APD signal was sent to the correlator board, which computed the intensity of the autocorrelation function. The correlation board streams output signals continuously via USB to a laptop PC for further data analysis.

The probe that coupled the optical fibers and recording electrode to the head was fabricated from semiflexible acrylate polymer using a three-dimensional (3-D) printer (Objet 260 Connex 3 printer, Stratasys, Minnesota). Embedded microprisms directed light from the fibers down to the head. The recording electrode was situated at the midpoint between the DCS optical source–detector separation. The midpoint of the probe [indicated in green in Fig. 1(d)] was thinner and permitted the probe to bend and conform to the curvature of the mouse’s skull, providing tight contact to the surface. The probe was positioned on a relatively flat region of intact skull above the forelimb’s representation in primary somatosensory cortex, 2.5-mm lateral of the midline at the same coordinate as bregma on the rostrocaudal axis.³³

2.5.

Modeling the Hemodynamic Response Function

To extract quantitative features from the sensory-evoked $\mathrm{\Delta }\mathrm{CBF}$, we fit the observed waveform to a canonical hemodynamic response function (HRF) that uses two-gamma density functions to approximate the hemodynamic response [Fig. 2(a)],³⁴ i.e.,

Fig. 2

Somatosensory-evoked hemodynamic and electrophysiological responses. (a) Optically measured changes in relative CBF following electrical stimulation of the median nerve. The gray trace represents the average of 10 trials in one experiment; blood flow values are normalized to the average of CBF values recorded in 4 s prior to the stimulus onset, which is indicated by the arrow. The dotted red line that is superimposed on the trace represents a fit to the observed hemodynamic response using a conventional HRF model (two-gamma functions), as described in Sec. 2. (b) SSEP recorded concurrently during acquisition of the optical data depicted in (a). The prominent peaks P1, N1, and P2 occur at temporal latencies that are well-established in the literature. The trace represents the average obtained from the same 10 trials in (a); however, for each trial, there are 12 SSEPs recorded in response to individual stimulus pulses (4 s of pulses delivered at 3 Hz); thus, the trace represents the average of 120 SSEPs. The green lines represent the widths of positive and negative peaks, measured at half of the maximum/minimum values. Note that the timescale in (a) is roughly $1000\times$ longer than in (b).

2.6.

Data Quality Criteria for Inclusion

Animals varied in the magnitude of noise associated with CBF and electrophysiological measurements. Sources of variability included optical coupling efficiency, electrode impedances, and physiological noise, among others. To best quantify injury-related changes, we only included data from experiments in which preinjury optical and electrophysiological noise levels were low enough to permit the observation of clear stimulus-evoked CBF changes and SSEPs. Quantitatively, our data inclusion criterion for optical measurements was that the initial CBF peak had an amplitude that exceeded twice the standard deviation (SD) of the baseline signal (i.e., $z\text{-}\text{score}>2$).

2.7.

Time-Frequency Analysis of the Hemodynamic Response Function

In the interest of parsing apparent oscillatory behavior associated with postinjury modifications in the hemodynamic response following the main positive and negative peaks, we performed spectrotemporal analysis on the CBF data throughout the entire duration of six experiments. Briefly, the continuous DCS data were divided into six 5-min epochs, and analysis of the evoked responses within those epochs was used to inform single “frames.” We utilized time-frequency transform functions from the MATLAB^® toolbox EEGLAB³⁵ to produce spectrograms that quantified power as well as intertrial coherence (ITC), which indicates the degree of phase locking in the observed $\mathrm{\Delta }\mathrm{CBF}$. Note that for each data point in Figs. 9 and 10, the spectrograms were based on only the trials within the 5-min epochs; although a larger sample size generally improves signal-to-noise, in this case because the $\mathrm{\Delta }\mathrm{CBF}$ was dynamically changing in time following injury, there was a trade-off between signal-to-noise and the temporal resolution with which the time evolution of power and ITC could be quantified. We empirically converged on 5-min epochs as a compromise based on observations from the continuous CBF that indicated that recovery generally occurred on a timescale slower than 5 min following a rapid decrement due to CCI.

2.8.

Histological Procedures

We assessed the effects of CCI on the brain by performing vital staining and immunohistological analysis. Two mice were subjected to CCI using the same parameters as all of the functional monitoring experiments (velocity: $4\mathrm{m}/\mathrm{s}$; dwell time: 100 ms; and impact depth: 0.6 mm). Another two animals were subjected to CCI of higher severity (impact depth was increased to 1.1 mm), and an additional two sham animals received no impact yet were positioned in the stereotaxic apparatus and the cortical impactor was fixed above the animal’s head. Animals were permitted to recover following the procedure. 24 h after CCI or sham treatment, animals were anesthetized with sodium pentobarbital ($100\mathrm{mg}/\mathrm{kg}$, i.p.) and perfused transcardially with saline and then formalin. Mice were then decapitated and the brain removed and postfixed in formalin overnight. $50\text{-}\mu \mathrm{m}$ coronal brain sections were acquired with a vibrating microtome at anteroposterior (AP) levels spanning bregma $+2\mathrm{mm}$ through $-2\mathrm{mm}$. For glial fibrillary acidic protein (GFAP) immunohistological analysis, floating slices were incubated overnight with monoclonal antimouse GFAP antibody (1:500 dilution, Invitrogen). The slices were then incubated with Peroxidase-conjugated AffiniPure Goat Antirat IgG ($\mathrm{H}+\mathrm{L}$) antibody (1:100 dilution, Jackson ImmunoResearch Laboratories) and stained with a 3,3-diaminobenzidine horseradish peroxidase substrate kit (Vector Laboratories), then mounted onto gelatin subbed slides and coverslipped. For Nissl staining, sections from the same brain, yet AP $+50\mu \mathrm{m}$ to the locations at which slices for GFAP analysis, were mounted and dried on gelatin subbed slides. They were then immersed in a heated solution of 0.1% cresyl violet (CV) acetate (Electron Microscopy Sciences) for 4 min and dehydrated in alcohol. Slides were then cover slipped with Permount (Fisher Scientific) and dried before imaging.

3.1.

Concurrent Measurements of Sensory-Evoked Cerebral Hemodynamics and SSEPs

Using the DCS component of our multimodal apparatus, we recorded sensory-evoked changes in CBF concurrently with electrophysiological measurements of SSEPs. Figure 2(a) shows an example of the preinjury CBF response. Averaged over 10 animals, the $\mathrm{\Delta }\mathrm{CBF}$ displayed an initial increase in blood flow of $14\%\pm 4.8\%$ ($\mathrm{mean}\pm \mathrm{SD}$) that peaked at $9\pm 0.3\mathrm{s}$, followed by a subsequent undershoot of magnitude $-12\%\pm 2.5\%$, which reached a negative peak $17\pm 1.1\mathrm{s}$ after the stimulus ($12\pm 1.1\mathrm{s}$ after the first peak). CBF returned to baseline values roughly 10 s following maximum undershoot. In terms of electrophysiological recordings, SSEPs displayed well-defined peaks, most prominently a positive deflection (P1) at $\sim 12\text{-}\mathrm{ms}$ poststimulus, a negative peak (N1) at $\sim 21\mathrm{ms}$, and a broader positive peak at $\sim 44\mathrm{ms}$ (P2) [Fig. 2(b)].

3.2.

TBI Acutely Alters Hemodynamic and Electrophysiological Functional Responses

Cortical impact acutely reduced the amplitude and altered the waveforms of sensory-evoked hemodynamic and electrophysiological responses (Fig. 3). Within the first 5 min following impact, the $\mathrm{\Delta }\mathrm{CBF}$ response peak decreased by $\sim 50\%$, and the undershoot amplitude decreased as well, though to a lesser extent ($\sim 20\%$ to 40%). The peak latencies additionally displayed changes on the order of $\sim 10\%$, though with opposite trends; the initial increase in CBF peaked $\sim 1\mathrm{s}$ earlier than preinjury, and the undershoot peaked $\sim 1\mathrm{s}$ later than in the preinjury case. The alterations in these parameters, based on the results of 10 mice, are summarized in Fig. 4. The peak-to-peak amplitude recovered to $\sim 80\%$ of its baseline value within 30 min, and during this period, the shifts in peak latencies resolved completely. It should be noted that despite being more sensitive to injury than other features of the $\mathrm{\Delta }\mathrm{CBF}$ waveform, the initial peak also recovered at a relatively fast rate, and within 5 to 10 min following injury, the peak amplitude had essentially recovered completely. In contrast, the undershoot alterations following injury, though more difficult to quantify because of the variability in that portion of the waveform post-CCI, displayed no recovery trend within the same time period.

Fig. 3

Sensory-evoked blood flow responses from forepaw stimulation for a representative experiment pre- and postinjury. The DCBF and SSEP traces shown here represent the average of multiple stimulations (average of 10 trials preinjury and 5 trials for each time point postinjury). The black line is the average, and the red error bars represent a SD.

Fig. 4

(a)–(d) The grand average summary of the effects of CCI on the amplitude and latency of prominent time-domain features of the $\mathrm{\Delta }\mathrm{CBF}$ waveform. The red dots represent the average values over all animals ($n=10$). Values are normalized to the baseline value (blue dot). The time of impact is designated here as $t=0\mathrm{s}$. Error bars represent SD. Significant changes relative to baseline are indicated by asterisks, which refer to $P$-values obtained from a two-tailed, paired student’s $t$-test ($*=P<0.05$, $**=P<0.01$). (e) A montage of sensory-evoked $\mathrm{\Delta }\mathrm{CBF}$ waveforms taken repeatedly over the course of an hour in a sham experiment.

SSEPs displayed similarly profound alterations within the first 5 min after CCI. As shown in Fig. 5, P1 fell to $\sim 40\%$ of its preinjury amplitude; however, N1 and P2 exhibited greater reductions of 81% and 78%, respectively. CCI additionally altered the peak latencies. Although the P1 latency increased by less than $\sim 14\%$, N1 and P2 both increased by up to $\sim 13\%$ and $\sim 50\%$. In general, P1 appeared to be the least sensitive to injury, in terms of both amplitude and latency. In the ensuing 30 min, the peak latencies and amplitudes largely—although not entirely—returned to their baseline values, yet at differing rates. The amplitude of P2, for example, recovered more rapidly than N1 or P1 in the first 15 min. Overall, however, following CCI, the initial decrease in SSEP amplitude (peak-to-peak, N1 to P2) was significantly more profound than the drop in CBF. For example, SSEP peak-to-peak amplitude was reduced by nearly 75%, while CBF peak-to-peak dropped by only $\sim 45\%$.

Fig. 5

(a)–(f) The grand average summary of the effects of CCI on major SSEP peaks. As in Fig. 4, red data points indicate postinjury averages over all animals $(n=10)\pm \mathrm{SD}$ and are normalized to baseline values (blue data points). The time of impact is designated here as $t=0\mathrm{s}$, and significance relative to baseline is indicated by asterisks, as in Fig. 4. (g) A montage of SSEP waveforms taken repeatedly over the course of an hour in a sham experiment.

In addition to alterations in average amplitude and latency of SSEP and CBF features, injury acutely increased their trial-to-trial variability of these parameters (Fig. 6). In our time-domain measurements, we quantified this variability in terms of the coefficient of variation, $C_v$, which is defined as the ratio of the SD to the mean value. In the first 5 min following CCI, $C_v$ increased by over 200% for all electrophysiological and blood flow parameters. Additionally, while the variability of other parameters returned to near baseline values within 10 min after injury, the CBF undershoot amplitude maintained an elevated variability over 30 min after injury, on average.

Fig. 6

Alterations in the variability of evoked responses following injury. (a) Dynamics of the $\mathrm{\Delta }\mathrm{CBF}$ variability are plotted over time, normalized to a baseline value (black dot). Injury was induced at 0 min in these plots. Variability is quantified here in terms of the coefficient of variation, $C_v$, defined as the ratio of the SD to the mean. $C_v$ is plotted for both the initial $\mathrm{\Delta }\mathrm{CBF}$ peak (in red) and the subsequent undershoot (purple). (b) Evolution of $C_v$ for SSEP peaks N1 (red), P1 (purple), and P2 (blue), normalized to baseline values (black dot). Error bars represent standard error of the mean, and $P$-value indications are as described in the figure insets.

3.3.

Histological Assessment of the Effects of CCI

To assess the repercussions of our implementation of cortical impact on the brain, we explored both macroscopic, tissue-level effects as well as GFAP expression in brains of animals that had been exposed to CCI (Figs. 7 and 8, respectively). We assessed effects at locations (1) close to the site of impact, which was anterior to the site of DCS and SSEP measurements, (2) at the same AP location as the optical and electrophysiological measurements, and (3) at a location distal (posterior) to the site of impact.

Fig. 7

CV staining reveals negligible macroscopic anatomical damage 24 h after CCI. Columns (a) and (b) depict CV-stained coronal brain sections from mice that received “mild” and “moderate” CCI treatment, respectively, and (c) depicts CV-stained sections from sham experiments. Here, “mild” CCI corresponds to the parameter settings employed for all of the functional data in this study (velocity: $4\mathrm{m}/\mathrm{s}$; dwell time: 100 ms; impact depth: 0.6 mm), while “moderate” corresponds to an increased impact depth (1.1 mm), with other parameters remaining the same. Slices from each condition are shown at locations on the AP axis corresponding to $+2\mathrm{mm}$, 0 mm, and $-2\mathrm{mm}$ relative to bregma. Square subregions that are indicated in the slices at bregma $+0\mathrm{mm}$ (middle row) are shown on an expanded scale (imaged at $40\times$, $\text{scale bar}=100\mu \mathrm{m}$) in (d)–(f).

Fig. 8

GFAP immunoreactivity 24 h after CCI indicates mild astrogliosis, but no macroscopic damage. Following the layout of Fig. 7, columns (a) and (b) depict GFAP immunoreactivity in coronal slices from mice that received mild and moderate CCI treatment, respectively, and (c) depicts results from sham experiments. Definitions of mild and moderate are as described in the caption for Fig. 7. Mild CCI corresponds to the parameter settings employed for all of the functional data in this study. The small notches on the ventrolateral corners of the brain slices indicate the hemisphere contralateral to the central axis of the impactor tip. Subregions that are indicated in the slices at bregma $+0\mathrm{mm}$ (middle row) are shown on an expanded scale in (d)–(f).

Figure 7 shows CV Nissl staining of brain sections from animals that had been exposed to CCI. Figure 7(a) shows CV-stained slices from an animal subjected to the same CCI parameters used during acquisition of all functional data (i.e., Figs. 3Fig. 4Fig. 5–6 and 9Fig. 10–11); we termed these settings “mild” for the sake of comparison with the histological results obtained following a slightly more intense impact (i.e., impact depth of 1.1 mm versus 0.8), which we termed “moderate.” Compared with sham animals, CV staining revealed no evidence of significant alterations in tissue morphology, such as persisting edema, major hematoma, or direct lesions resulting from the CCI settings used throughout this study (two animals per condition). GFAP immunoreactivity (Fig. 8), however, revealed a proliferation of reactive astrocytes biased toward the hemisphere that was exposed to a greater area of the impactor tip (which was 5 mm in diameter). This astrogliosis effect was evident at the site on the AP axis where functional measurements were centered and was more pronounced at moderate CCI (impact depth 1.1 mm) compared with mild (0.8 mm). Although GFAP immunoreactivity appeared to be more globally intense at the most distal region assessed (AP $-2\mathrm{mm}$), this is consistent with baseline expression patterns of GFAP,³⁶ and there was no apparent difference relative to the sham condition.

Fig. 9

Time-frequency analysis reveals significant changes in spectrogram after injury. (a) Representative time-frequency spectrograms from a single animal showing the alterations in spectrotemporal content following CCI. The white arrows indicate the time point at which forelimb stimulation was applied. (b) Spectrotemporal changes over time, averaged over all animals. Each image depicts, in grayscale, the difference in spectrotemporal content via subtraction of the preinjury spectrogram (darker shades indicate decreases, brighter shades indicate increases). Superimposed in color on the grayscale difference spectrogram are points that, when analyzed over all experiments, demonstrated a difference from the preinjury values with $P<0.05$ using a two-tailed student’s $t$-test. The colors are defined as in the colorbar of (a), and note that nearly all of the spectrotemporal power changes are focal increases. (c) Comparison of a representative postinjury $\mathrm{\Delta }\mathrm{CBF}$ waveform (here, at a latency of 25 min following injury) with a plot showing the sum, over all experiments, of pixels showing alterations from baseline with $P<0.05$. Note that the significantly changed spectrotemporal features lined up with the cycles of apparent “oscillations” following injury. The gray shading projected onto the postinjury $\mathrm{\Delta }\mathrm{CBF}$ waveform (a representative plot) highlights temporal periods that contained statistically significant alterations in spectral power, assessed over all experiments ($n=10$).

Fig. 10

ITC alterations following controlled cortical injury. (a) Representative spectrograms depicting the absolute magnitude of ITC before (top) and after (bottom) injury. These particular spectrograms are derived from the same animal as in Fig. 9(a). The color definitions are as indicated in the ITC colorbar; a value of 1 represents perfect, in-phase spectral responses. While the depicted values are real values, they represent the absolute magnitude of complex values that originally contain phase information. The black arrows indicate the time point at which forelimb stimulation was applied. (b) A series of ITC spectrograms, averaged over all animals, showing features that differ significantly from preinjury with $P<0.05$ obtained from a two-tailed student’s t-test. As in Fig. 9(b), the statistically significant pixels, which are color-coded according to the colorbar of (a), are superimposed on a grayscale spectrogram representing the postinjury ITC spectrogram from which preinjury spectrogram values have been subtracted. While most of the changes represent increases in ITC at time points following the initial $\mathrm{\Delta }\mathrm{CBF}$ waveform peak and undershoot, there was a significant decrease in ITC in the region of the undershoot itself, roughly the area of frequencies below 0.4 Hz between 12- and 16-s poststimulus (indicated by the white arrow in the top left superimposed image). (c) Regional assessment of ITC changes following injury. A sample $\mathrm{\Delta }\mathrm{CBF}$ waveform fit divided into three regions for assessing changes in ITC at particular phases of the functional blood flow response is shown on top. The specific regions were defined as (1) an initial “rise” extending from the time directly after sensory stimulation cessation up through the point where the second derivative of the $\mathrm{\Delta }\mathrm{CBF}$ fit passes zero; (2) the undershoot period, extending from the end of (1) through the latest point of intersection of the $\mathrm{\Delta }\mathrm{CBF}$ fit with the undershoot half-maximum (or minimum, here) line; (3) the period of time following the main two phases of the $\mathrm{\Delta }\mathrm{CBF}$ waveform and extending until the subsequent stimulus period. Asterisks in the right-most bar chart on the bottom indicate values that differ from baseline, preinjury values with $P<0.05$ according to a two-tailed Student’s t-test.

Fig. 11

Correlation matrix comparing the post-CCI recovery trends for the prominent features of the CBF and SSEP responses. Each matrix element represents the Pearson’s correlation coefficient obtained by comparing the post-CCI time course (as in Figs. 4 and 5) of the major peaks and latencies. A correlation coefficient value of 1 represents identical waveforms, $-1$ indicates that the two traces are perfectly anticorrelated. The matrix is symmetrical about the diagonal. Correlation coefficients with an absolute value higher than 0.8 are indicated with white dots. The dotted rectangular perimeter delimits the matrix region representing correlation between the two modalities, optical and electrophysiological. The yellow perimeter highlights the high correlation between the P2 amplitude and postinjury hemodynamics, whereas the green perimeter highlights the fact that the $\mathrm{\Delta }\mathrm{CBF}$ peak latency is highly correlated with all aspects of SSEPs.

3.4.

Spectrotemporal Analysis the Hemodynamic Response Function Following TBI

Modulation of prominent aspects of the hemodynamic response such as peak amplitude and latency are straightforward to quantify because the waveforms can fit well to a conventional HRF model (as described in Sec. 2). However, injury also induced new temporal features into the CBF response. For example, in the $\sim 30\mathrm{s}$ following each evoked peak in CBF, there emerged slow oscillations in flow at a frequency of roughly 0.1 to 0.3 Hz. Such low-frequency oscillations have been reported both as spontaneous fluctuations in ongoing flow^37–39 as well as following functional activation.⁴⁰ These emerging oscillations introduce additional waveform peaks that cannot be effectively fit using the standard, two-gamma function model. As an alternative approach for quantifying these new components of the CBF response, we performed time-frequency spectral analysis to explore and more directly visualize these effects. While not typically applied to evoked responses, such spectral analysis has been used to explore dynamic, low-frequency aspects of ongoing cerebral autoregulation⁴¹ and provides an approach for characterizing the putative oscillatory behavior. As shown in Fig. 9, injury elicited changes in spectral power throughout the duration of the major peaks (increase and undershoot) of the CBF response, as well as in the interim between stimulus sets. Although the apparent spectrotemporal changes following injury resemble the development of temporally separated “notches” in the spectral power for frequencies below 1 Hz, we explored whether these changes were significant across animals by assessing the patterns of spectrotemporal differences from preinjury values [Fig. 9(b)]. The regions of statistically significant spectrotemporal changes in power are relatively sparse; however, when summed over frequency and for all experiments, significant changes coincided with individual peaks in the “oscillatory” region, indicating a degree of phase locking [Fig. 9(c)].

To quantify this apparent coherence effect, we explored the ITC (Fig. 10), which can reveal spectrotemporal phase locking and variability among evoked responses.³⁵ Similar to the results of spectral power analysis shown in Fig. 9, injury qualitatively changed ITC most prominently for regions following the initial CBF peak. A large fraction of the alterations in ITC were increases, most prominently in the period after the CBF undershoot and for frequencies above 0.4 Hz. The one region that displayed a significant decrease in ITC was that corresponding to the spectrotemporal representation of the CBF undershoot, indicating greater variability of that portion of the HRF after injury. Although significant spectrotemporal increases in ITC ($P<0.05$ via two-tailed student’s t-test) were apparent throughout the time course of the CBF response, we investigated the temporal distribution of ITC alterations [Fig. 10(c)] by analyzing flow as a series of three “regions” defined by physiologically relevant landmarks. When summed over frequencies, the period following the $\mathrm{\Delta }\mathrm{CBF}$ undershoot (region three) demonstrated significant alterations from baseline.

3.5.

Correlations in the Dynamics of Cerebral Blood Flow and Electrophysiological Responses Following TBI

Given the diversity of signals that were sensitive to injury, we sought to merge the findings by exploring correlations among their recovery kinetics. The matrix in Fig. 11 presents correlation coefficients comparing the post-CCI quantitative dynamics shown in Figs. 4 and 5. Overall, SSEP peak amplitudes and latencies were highly intercorrelated. In comparison, the recovery kinetics of the CBF components was relatively weakly intercorrelated. Among the cross-correlation terms indicative of coupling between hemodynamic and neural activity, we identified two major trends that are highlighted in Fig. 11. Notably, of all time-domain electrophysiological features that we explored, the amplitude of the CBF peak was most significantly correlated with the amplitude of P2 ($z\text{-}\text{score}>3$, based on the SD of the distribution of Pearson’s correlation coefficients). Additionally, while P2 represented the strongest cross-modal correlation in terms of amplitude, the latency of the initial $\mathrm{\Delta }\mathrm{CBF}$ peak was highly correlated with all aspects of the SSEP peaks.