In the United States, each year there are over 1.5 million traumatic brain injuries (TBIs), resulting in 50,000 deaths.1 TBI induces temporary or permanent impairment of cognition, physical function, and psychosocial behavior.2 TBI severity ranges widely depending on the nature of the injury. The majority of these injuries are classified as “mild” TBI (mTBI), which are challenging to diagnose and track because quantitative biomarkers for mTBI are lacking. In fact, one of the defining characteristics of mTBI is that it cannot be validated through standard methods of clinical imaging.3 Given the difficulty of rapid detection, mTBI poses a particular challenge to public health because repeated injuries such as concussions have a cumulative effect on brain health.4–6
Among a wide spectrum of TBI sequelae, sensory and cognitive deficits are some of the most common.7 A difficulty in assessing these deficits is that they are frequently mutually confounding. For example, auditory dysfunction following injury includes both peripheral deficits, e.g., increased hearing thresholds, as well as difficulties with more sophisticated auditory tasks, such as discriminating sounds in noisy environments or comprehending speech.8 Visual deficits following TBI are similarly complex; the most common complaints, for example, are accommodative deficiencies,9 which include blurred vision, headache, motion sickness, or loss of concentration during visual task performance.10 While auditory and visual performance may be difficult to dissociate from cognitive deficits, olfaction, which does not involve significant feedback with subcortical structures, is also impacted by TBI, and the degree of anosmia is correlated with the severity of injury.11
The integrity of sensory systems can be probed electrophysiologically via evoked potentials (EPs). EPs can be measured noninvasively and comprise a series of positive and negative voltage deflections, which reflect the afferent relay and processing of sensory information. Auditory-, visual-, and somatosensory-evoked potentials (SSEPs) are routinely used in multiple clinical contexts to aid neurological assessment of TBI and to monitor functional recovery over time.12 Sensory-evoked neural activity can also be inferred based on cerebral hemodynamic responses, including changes in cerebral blood flow (), volume, and oxygenation that are evoked by brief sensory stimuli. Damage to any portion of the cerebral vasculature fundamentally alters the ability of the network to supply neurons with energy.13 In fact, blockage of even single capillaries can cause larger-scale changes in blood flow14 and may result in microvascular ischemia.15 Repercussions of these injuries, even if mild, may continue to progress following the primary trauma, ultimately leading to more global sequelae.
We explored the possibility of deriving innovative indicators for mTBI based on intrinsic correlations between hemodynamic and neuronal activity. Specifically, we investigated signals that can be measured rapidly and using portable technology, which is a prerequisite for expedited assessment. To noninvasively monitor sensory-evoked hemodynamics, we used diffuse correlation spectroscopy (DCS) that takes advantage of the dynamic scattering properties of red blood cells to directly measure cerebral blood flow (CBF). DCS is particularly sensitive to flow in the cortical microvasculature due to the high absorption (and thus low probability of photon escape) in larger blood vessels.16 DCS has been used to measure functional hemodynamics associated with sensory stimuli and motor tasks17,18 and has been used to track baseline CBF following brain trauma.19 We supplemented our ongoing optical recordings of the sensory-evoked with concomitant measurements of SSEPs and applied this multimodal approach to a mouse model of TBI that employed controlled cortical impact (CCI) as the source of primary injury.20–22 The ability to noninvasively monitor both aspects of neural response enabled us to obtain the first detailed, in vivo portrait of the effects of acute injury on sensory processing in the brain.
Materials and Methods
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 . Animals were additionally administered an injection of xylazine (, 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 () were administered once every . 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 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: ) 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.
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.
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 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 , and total pulse-train duration 4 s). Individual measurements of sensory-evoked 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.
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 , dwell time 100 ms, and impact depth 0.6 mm. Strike velocities and/or resulted in skull fracture. Additionally, impactor tips with smaller diameter tended to cause fractures. The metal impactor tip was positioned over the exposed skull anterior and lateral of bregma, contralateral to the optical measurement.
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 ( core diameter, Thorlabs, New Jersey), and the scattered light was detected with a single-mode fiber ( 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.33
Modeling the Hemodynamic Response Function
To extract quantitative features from the sensory-evoked , 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)],34 i.e.,
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., ).
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 EEGLAB35 to produce spectrograms that quantified power as well as intertrial coherence (ITC), which indicates the degree of phase locking in the observed . 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 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.
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: ; 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 (, i.p.) and perfused transcardially with saline and then formalin. Mice were then decapitated and the brain removed and postfixed in formalin overnight. coronal brain sections were acquired with a vibrating microtome at anteroposterior (AP) levels spanning bregma through . 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 () 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 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.
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 displayed an initial increase in blood flow of () that peaked at , followed by a subsequent undershoot of magnitude , which reached a negative peak after the stimulus ( 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 poststimulus, a negative peak (N1) at , and a broader positive peak at (P2) [Fig. 2(b)].
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 response peak decreased by , and the undershoot amplitude decreased as well, though to a lesser extent ( to 40%). The peak latencies additionally displayed changes on the order of , though with opposite trends; the initial increase in CBF peaked earlier than preinjury, and the undershoot peaked 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 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 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.
SSEPs displayed similarly profound alterations within the first 5 min after CCI. As shown in Fig. 5, P1 fell to 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 , N1 and P2 both increased by up to and . 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 .
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, , which is defined as the ratio of the SD to the mean value. In the first 5 min following CCI, 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.
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.
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 ), this is consistent with baseline expression patterns of GFAP,36 and there was no apparent difference relative to the sham condition.
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 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 flow37–39 as well as following functional activation.40 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 autoregulation41 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.35 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 ( 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 undershoot (region three) demonstrated significant alterations from baseline.
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 (, 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 peak was highly correlated with all aspects of the SSEP peaks.
A major repercussion of traumatic injury in the brain is an ensuing mismatch between oxygen supply and neural metabolic demand.42,43 Changes in CBF and local tissue oxygenation dynamics ultimately reflect an interplay between CBF and changes in the cerebral metabolic rate of oxygen (), both of which are susceptible to alteration following acute injury. The primary mechanical injury, for example, can alter the vascular network that delivers oxygen to the parenchyma and, simultaneously, can alter neural activity and metabolic rate through a wide spectrum of mechanisms, including mechanical shearing44 or excitotoxicity due to damaged glia.45 Our concomitant measurements of CBF and electrical activity provide a picture of the repercussions of injury on both “supply and demand” aspects of the networks underlying sensation. While it is possible for injury to globally affect brain oxygenation (and, indeed, systemic physiology) as a result of compression of cardiorespiratory brainstem nuclei, our histological assessment does not find evidence of this scenario, and the changes in physiological signals likely reflect more subtle, transient perturbations.
While the general trends in SSEP alterations following injury are consistent with previous work in animal models,23,25,26,46 in the present study, we were able to identify time-domain alterations in the SSEP in finer detail. The most prominent SSEP peaks—N1, P1, and N2—exhibited differing alterations after injury, some of which are consistent with the known anatomical correlates of time-domain waveform. For instance, considering that the site of primary injury in CCI is the cortical gray matter,47 it is not entirely surprising that P1 ( onset delay preinjury) exhibited the mildest decrease following injury given the fact that this component of the SSEP reflects subcortical activity related to thalamocortical relay.48 Likewise, the high sensitivity of P2 likely reflects the fact that its underlying generators are derived from cortico-cortical processing, which is sensitive to lateral cortical perturbations. For example, modulating cortical activity with topical application of pharmacological agents generally affects N1, but largely spares P1.49 Consistent with this, in our previous work using a high-intensity focused ultrasound model for explosive blast, the mechanical insult penetrated deeper beneath the surface of the cortex and P1 was more significantly altered.23 The high correlation between hemodynamics and P2 may also reflect an intrinsic, disproportionately large influence of superficial cortical injury on electrical activity and hemodynamics within the cortical microvasculature. The superficial layers of the cortex feature extensive lateral networking with other cortical areas;50 thus, local disruption due to impact injury would likely have diffuse and far-reaching repercussions.
Within the temporal span of individual hemodynamic responses, the major features of the preinjury that we observed—including both an initial peak and subsequent undershoot—have been observed in other DCS studies51,52 as well as studies using laser speckle imaging.17 While evidence of a significant undershoot is lacking in laser Doppler measurements, a clear inflection point between the CBF peak and return to baseline has been observed.53 A significant undershoot has also been reported in some fMRI studies, measured through arterial spin labeling,54–56 but not all.57 A major contributing factor underlying discrepancies in CBF results is likely the anatomical specificity of the measurement modality. For example, DCS measurements and speckle flow imaging are highly sensitive to flow in smaller vessels, whereas laser Doppler, which yields a minimal undershoot, is biased toward signals in large vessels, which generally exhibit smaller fractional changes in flow during functional activation.
The initial CBF peak reflects the influx of blood flow recruited by local changes in neural activity,58 and, on the timescale of days, its amplitude has been found to be reduced following traumatic injury59 and stroke.60 To our knowledge, however, the impact of injury on sensory-evoked CBF, both in terms of amplitude and temporal waveform, has not been explored within the first moments following TBI. It remains to be seen whether our observed acute dynamics, which depict significant recovery of some aspects of the sensory-evoked , continuously evolve into the previously reported long-term deficits. In terms of injury-induced alterations in peak CBF timing, the same prior investigations reported a delayed flow onset that is observable on the timescale of days following injury. Our observation of a reduced onset latency is innovative; however, it may be unique to the immediate aftermath of acute injury. In general, the fact that most aspects of the response and SSEPs substantially recover within a relatively short period of time is remarkable given the known chronic effects that TBI exerts on baseline hemodynamics, notably, changes in cerebrovascular autoregulation.61,62 The findings concur with previous reports that functional CBF is independent of baseline values;53,63,64 however, our results indicate that this independence extends to injury-induced global hemodynamic changes.
It has long been assumed that the sensory-evoked hemodynamic response undershoot is primarily a vascular phenomenon, representing a poststimulus period of increased cerebral blood volume even after CBF and have returned to baseline.65 However, more recent work strongly indicates that the poststimulus undershoot is highly dependent on neural activity.56 Our results depict the postinjury undershoot as chronically attenuated and more variable. Similarly, SSEP components also displayed a marked increase in variability in both time and frequency domain. These observations are, therefore, consistent with the hypothesis that neural activity plays a significant role in shaping the undershoot. However, the fact that the variability and amplitude of the undershoot does not recover as fully as other CBF waveform components even after SSEP attributes have returned to baseline status suggests that the prolonged undershoot alterations may not be primarily shaped by neural activity. Future studies involving systematic variation of stimulus parameters may be able to more directly probe this hypothesis.
Using a variety of techniques, low-frequency () CBF oscillations similar to the ones we detected following injury (Figs. 9 and 10) have been previously observed, both spontaneously and following stimulus-evoked CBF, in the form of vasomotion.37,40,66,67 While the physiological basis of these oscillations is still currently unclear, multiple mechanisms have been suggested, including emergent properties of molecular signaling pathways68 or more fundamental fluid mechanical principles.69,70 Such fluid mechanical models derive an oscillatory component based on inclusion of an inertial term, owing to blood volume in larger, more distal venules and arterioles.71 A simplistic, general interpretation of the results is that the emergence of postundershoot oscillations in CBF reflects perturbation of an intrinsically oscillatory system, regardless of whether the affected components are mechanical or phenomenological (e.g., related to molecular signaling).
In emergency medicine, there is a well-known principle of the “golden hour,” the period of time following a traumatic injury within which medical intervention can have a maximal impact on patient outcome.72 Although the precise duration is debated,73,74 timely transport to a hospital is essential. Our experimental investigation focused on disruptions in neurovascular coupling that occurs within this relatively short time window, and we identified several metrics that were profoundly altered by TBI. It is, however, both possible and likely that at later time points, the trends we observed will change or reverse. The rapid energy depletion associated with disrupted CBF and mechanical shearing of neurons dramatically increases extracellular and glutamate concentrations.75,76 In turn, resulting excitotoxicity triggers a variety of rapidly evolving adverse events through the “ischemic cascade.” For example, initial cytotoxic edema is followed by the onset of vasogenic edema within the first few hours.77,78 These various sequelae will likely be reflected in the optical and electrophysiological signals. Future experimental work focusing on the longer-term evolution of neurovascular coupling following CCI could provide a fuller physiological context for the initial effects that we characterized in this present study. On a technical level, focusing on physiological signals during the early timeframe of injury provided a way to mitigate potential confounds owing to anesthesia pharmacodynamics, which could be affected by injury.
A general challenge for parsing the bulk signals obtained through diffuse optical neuromonitoring techniques is that the spatial resolution is limited because the detected signals include light that has traversed a large span of optical paths in tissue.79 In the case of DCS, although the signal origin is putatively weighted toward microvasculature, the fidelity with which the bulk signal aligns with a microscale portrait of blood flow in the layers of the cerebral cortex remains to be quantified. Future work in animal models of TBI that supplements DCS measurements with high-resolution imaging of activity in microvasculature, such as optical coherence angiography80,81 or two-photon microscopy,82 may better inform interpretation.
While significant progress has been made toward understanding the repercussions and mechanisms of mTBI at the cellular level using invasive, high-resolution techniques, biomarkers for rapidly assessing human subjects in active settings outside the clinic are constrained to signals that can be obtained with noninvasive technology. Such measurement modalities, however, are limited in their ability to spatially resolve signals and biomarkers observed in preclinical work. As we have demonstrated, functional neurovascular coupling is a property that can be measured with low-resolution techniques and contains information about microscale physiological function. This information can be extracted from variable, high-noise physiological recordings because they are temporally synchronized with external stimuli that can be delivered in a controlled manner. The event-related optical and electrophysiological changes that we have discovered, albeit correlative, are significantly indicative of TBI in an animal model. Notably, the dynamic alterations we observed in stimulus-evoked blood flow and electrophysiological responses were not associated with significant tissue damage; rather, our histological examination revealed only a small increase in astrogliosis associated with CCI. We expect that these functional indicators will ultimately add informative dimensions to current multimodal approaches for diagnosing mTBI in human subjects, given the heterogeneity of baseline values across populations.
The authors report no financial interests or potential conflicts of interest related to this work.
The mention of commercial products, their sources, or their use in connection with material reported herein is not to be construed as either an actual or implied endorsement of such products by the Department of Health and Human Services. The authors thank Dr. A. Parthasarathy, Dr. W. Baker, Dr. T. Durduran, Dr. A. G. Yodh, Dr. G. Yu, Dr. D. Busch, Dr. A. Koller, Dr. D. Ress, and Dr. J. Pfefer for helpful discussions, and Dr. S. Vasudevan for technical assistance. The work was supported by the National Science Foundation (Award Nos. 1541612 and 1641133) as well as intramural recruitment funds from New York Medical College.
Hyounguk Jang received his BS degree and PhD in physics from the University of Sungkyunkwan, Korea, in 2002 and Kansas State University in 2010, respectively. He performed postdoctoral research in the Department of Physics, University of Maryland, from 2013 and the Department of Bioengineering, University of Maryland, from 2015. Currently, he is a postdoctoral researcher in the Department of Physiology, New York Medical College, and a scholar-in-residence at the U.S. Food and Drug Administration (FDA).
Stanley Huang received his BS and MS degrees in biomedical engineering from Johns Hopkins University, Baltimore, Maryland, in 2003 and 2005, respectively. He received his PhD in biomedical engineering from Columbia University, New York City in 2011. Currently, he is a research fellow at the U.S. FDA. His research interests include biosignal processing, bioelectronics, medical imaging, and neurological devices.
Daniel X. Hammer has over 20 years of experience in the invention, development, and application of biomedical optics to a broad range of interdisciplinary research programs. His current research interests include optical coherence tomography, adaptive optics, and neuroimaging. He is currently the deputy director of the Division of Biomedical Physics (DBP) within the Center for Devices and Radiological Health’s (CDRH) Office of Science and Engineering Laboratories (OSEL) at the U.S. Food and Drug Administration (FDA) in Silver Spring, Maryland. He is a member of SPIE (senior member), the Optical Society of America, the Association for Research in Vision and Ophthalmology, and the American Institute for Medical and Biological Engineering (fellow).
Lin Wang was a graduate student at New York Medical College and is currently a medical student at the Medical College of Wisconsin. Her current research interests include improving medical diagnostic tools and modalities, particularly for mild traumatic brain injury (TBI). She has also been involved in the design of a portable device for diagnosing and quantifying neurovascular alternations after mild TBI at New York Medical College.
Harmain Rafi received her BS degree in neuroscience from George Mason University, Fairfax, Virginia, in 2016. Her previous work includes retrospective research on spinal cord stimulators at INOVA Fairfax hospital and locomotor behavioral studies in rats at the Krasnow Institute. Currently, she is an ORISE fellow at the U.S. FDA in Silver Spring, Maryland, working in the Neural Implants Lab.
Meijun Ye received her PhD in neuroscience from the University of Arkansas for Medical Sciences. After that she did her postdoctoral research in visual sciences at Yale University. She joined the Neurological Devices Laboratory at the U.S. FDA in 2013. Her research interests include neurophysiological signatures of TBI, safety and effectiveness of neural diagnostic and intervention devices, cellular and circuitry mechanisms of brain injury, retinal neural circuits, and neural interactions in the reticular activating system.
Cristin G. Welle received her BS degree from the College of William and Mary in 2004 and her PhD from the University of Pennsylvania in 2010. She is an assistant professor at the University of Colorado, Denver. Previously, she was a staff scientist at the U.S. FDA, where she led the Neural Interfaces Laboratory. Her research interests center on using small animal models to explore long-term safety and performance of medical devices that interface with the nervous system.
Jonathan A. N. Fisher received his BA degree and PhD in physics from the University of Pennsylvania in 2000 and 2007, respectively. He performed his postdoctoral research in neuroscience at The Rockefeller University. From 2013 to 2014, he was an ORISE faculty fellow at the U.S. FDA/CDRH/OSEL/Division of Biomedical Physics. Currently, he is an assistant professor of physiology at New York Medical College and his awards include a Blavatnik Award for Young Scientists.