5 December 2017 Longitudinal, transcranial measurement of functional activation in the rat brain by diffuse correlation spectroscopy
Author Affiliations +
Neurophotonics, 4(4), 045006 (2017). doi:10.1117/1.NPh.4.4.045006
Neural activity is an important biomarker for the presence of neurodegenerative diseases, cerebrovascular alterations, and brain trauma; furthermore, it is a surrogate marker for treatment effects. These pathologies may occur and evolve in a long time-period, thus, noninvasive, transcutaneous techniques are necessary to allow a longitudinal follow-up. In the present work, we have customized noninvasive, transcutaneous, diffuse correlation spectroscopy (DCS) to localize changes in cerebral blood flow (CBF) induced by neural activity. We were able to detect changes in CBF in the somatosensory cortex by using a model of electrical forepaw stimulation in rats. The suitability of DCS measurements for longitudinal monitoring was demonstrated by performing multiple sessions with the same animals at different ages (from 6 to 18 months). In addition, functional DCS has been cross-validated by comparison with functional magnetic resonance imaging (fMRI) in the same animals in a subset of the time-points. The overall results obtained with transcutaneous DCS demonstrates that it can be utilized in longitudinal studies safely and reproducibly to locate changes in CBF induced by neural activity in the small animal brain.
Blanco, Zirak, Dragojević, Castellvi, Durduran, and Justicia: Longitudinal, transcranial measurement of functional activation in the rat brain by diffuse correlation spectroscopy



Alterations in neuronal activity drive local changes in blood flow in order to increase blood supply to the metabolically active regions.1 The neurovascular coupling is one of the primary biomarkers utilized in the study of functional brain activity.

Functional neuroimaging techniques are especially interesting for the study of the progression of different pathologies23.4.5.6 affecting the central nervous system, as well as different mechanisms of recovery, for instance, therapy-potentiated functional recovery.7 This requires systems that do not interfere with the progression of either the pathology or the possible therapies, and allow the same individuals to be studied over extended periods.

The most common noninvasive neuroimaging method is functional magnetic resonance imaging (fMRI). fMRI is based on the blood oxygen level dependent (BOLD) signal,8 which is associated with deoxyhemoglobin concentration. Unfortunately, changes in the BOLD signal can be influenced by vascular architecture, intrinsic hemodynamic responses, and other factors, such as spin density, volume fraction, the levels of different neurotransmitters, and perfusion.9 Even though fMRI is considered as one of the most applicable techniques for the study of neuronal activity, it has some drawbacks, such as availability, cost, and robustness.

On the other hand, there are optical imaging modalities that are relatively inexpensive, which do not need a special infrastructure or use of an exogenous contrast. Optical techniques can provide advantages in the study of the functional activation in the brain by using intrinsic changes that can detect variations in fluorescence light, absorption, or light scattering. Techniques, such as optical intrinsic signal imaging (OISI), functional near-infrared spectroscopy, or photoacoustics study changes in the blood volume and blood oxygen saturation due to changes in the concentration of oxy- and deoxyhemoglobin induced by the neuronal activity.1011.12.13 Indirect information about cerebral blood flow (CBF) can be obtained by looking at the changes of the blood volume over time by using contrast agents. There are optical techniques that directly measure blood flow, such as diffuse correlation spectroscopy (DCS), laser speckle imaging (LSI), and laser Doppler flowmetry (LDF).1415.16.17.18 Other MRI modalities, such as cerebral blood volume (CBV)-weighted fMRI, where an intravascular contrast agent is necessary to monitor CBV changes induced by local neural activity, or arterial spin labeling (ASL), that localizes local changes in CBF, have a relatively low temporal resolution and low signal-to-noise ratio (SNR).19

Optical techniques are often used for imaging the hemodynamic response to different functional stimuli in the small animal brain. LDF is used to quantify the blood flow changes during the stimulation in small animal models, however, it is limited to point and superficial tissue measurements (<1  mm).2021.22 Another technique that is often used is LSI. LSI provides two dimensional images and can measure absolute blood flow,23,24 but it is also limited to superficial imaging (<1  mm) of hemodynamics.16,25,26 The downside of both LDF and LSI is that it is necessary to clear, remove, or thin the skull depending on the species.23,25 On the other hand, diffuse optical techniques, such as near-infrared spectroscopy (NIRS), DCS,27 and speckle contrast optical tomography (SCOT)28 allow noninvasive measurements of the cerebral hemodynamics in the deep tissues.12,25,29,30 Both NIRS and DCS have been used for measurement of rodent brain transcranially31 with scalp retraction27,28 and with implanted probes in mice heads for longitudinal studies.32

In the present study, we combine a mild sedative with transcranial DCS to measure CBF in longitudinal (over months) brain functional activation due to forepaw stimulation and to cross-validate the activation pattern with results obtained by fMRI in the same rats by employing the routinely used BOLD signal.



All experimental procedures involving the use of laboratory animals have been carried out following the ARRIVE guidelines, the local and European legislation, and have been approved by the Ethical Committee of Animal Experimentation (approval number 7883) of the Generalitat de Catalunya.

Eleven male Wistar rats with a body weight of (320±50)  g at the beginning of the experiment (3 to 4 months old) were initially used for fMRI, and the same rats were measured with DCS at the age of 6 and 18 months. Six DCS experiments were carried out at the age of 6 months, with a resting period of 15 days between sessions, and an additional four experiments were done with the same rats at the age of 18 months [Fig. 1(a)].

Fig. 1

(a) Timing of the experimental procedures. (b) Experimental protocol for induction of anesthesia and stimulation session. Fifteen minutes under isoflurane (ISO) were dedicated to positioning the rats’ heads on the stereotaxic frame (or on the cradle in the fMRI experiments). Fifteen minutes after medetomidine bolus injection (MED), the infusion started. Thirty minutes were allowed for sedative stabilization before the first electrical stimulation. Each session consisted in six consecutive experiments, alternating right and left forepaw stimulation. Rectangular pulses (0.3 ms; 6 Hz; 2 mA) were used in five blocks (15-s activation and 45-s resting period each), and 5 minutes were left for recovery between experiments.




Rats were anesthetized with 4% isoflurane in O2:N2O (30:70). Animals were placed in a stereotaxic frame (or the MRI cradle in the fMRI experiments) with a face mask delivering the same gas mixture with 1.5% isoflurane, and a subcutaneous bolus (0.05  mg/Kg) of medetomidine (Domtor, Pfizer) was injected. Isoflurane was slowly discontinued over the course of the next 15 min at a rate of 0.25% every 1 min starting 3 min after the bolus. At this point, a continuous subcutaneous infusion of medetomidine (1  ml/h and 0.1  mg/Kg) was initiated for the whole duration of the experiment.

The scalp of the rat was shaved using an electric razor and depilatory cream to avoid side effects in the optical signal produced by fur. Two subdermal electrode needles were inserted in each forepaw for stimulation. Respiration rate (number of breaths per minute) was constantly registered and the body temperature was monitored with a rectal probe and maintained at (37±0.5)°C with a feedback controlled electrical blanket. Stimulation experiments started 30 min after isoflurane was discontinued and rats presented a stable condition with a respiration rate of 40% to 50% compared to the initial values.

Once the experiment was concluded, animals received an intraperitoneal injection (0.1  mg/Kg) of atipemazole (Antisedan, Pfizer) to reverse the effect of medetomidine.


Electrical Forepaw Stimulation

Forepaw stimulation consisted in rectangular pulses (2.0 mA, 6 Hz, 0.3 ms) in a paradigm of five consecutive blocks of 15-s stimulation followed by 45 s of recovery [Fig. 1(b)] with a total time of 5 minutes. Functional activation monitoring was conducted alternately three times for each forepaw. A resting period of 5 min was allowed between forepaw stimulation experiments [Fig. 1(b)].


Functional Magnetic Resonance Imaging

The fMRI experiments were conducted on a 7.0 T BioSpec 70/30 horizontal scanner (Bruker BioSpin, Ettlingen, Germany) equipped with an actively shielded gradient system (400  mT/m, 12-cm inner diameter). The receiver coil was a four-channel phased-array surface coil for the rat brain. Functional activation imaging was achieved with BOLD contrast MRI. Coronal multislice spin-echo (SE) EPI images were acquired using the following parameters: TE/TR=30=3000  ms; BW=150  kHz; five consecutive slices of 2 mm thickness; field-of-view=2.56×2.56  cm2; matrix of 64×64  pixels.


Data analysis

Statistical parametric activation maps from BOLD fMRI were constructed with the software STIMULATE.33 The time course was examined for each pixel during forepaw stimulation using a paired Student’s t-test with p<0.01 as the significance level, where only statically significant pixels were used to calculate median BOLD signal intensity [ΔSI BOLD (%)].

fMRI maps were calculated without any previous filtering or signal baseline thresholding. Clustering was done considering at least four adjacent activated pixels.

DCS signal was normalized to the 5-min baseline prior to the corresponding stimulation experiment.

Statistical analyses were conducted using GraphPad Prism v4.0. All data are presented as (mean±SEM)%, where SEM is the standard error of the mean. To check if the ΔSI BOLD signal and ΔrCBF are correlated, we have carried out a linear regression analysis and calculated the squared Pearson’s correlation coefficient (r2).


Diffuse Correlation Spectroscopy


Optical method

DCS is a noninvasive diffuse optical method that quantifies deep tissues blood flow. DCS measures temporal speckle fluctuations of the scattered light that are sensitive to the movement of the scattering particles, such as red blood cells.14,34 The dynamics of the scatterers is determined by the measurement of the intensity autocorrelation function, from which the normalized electric autocorrelation function is derived. This is then fitted with the solution of the correlation diffusion equation for the semi-infinite geometry to obtain the effective, scattering number weighted Brownian motion coefficient (DB), which was shown to be proportional to the blood flow index (BFI).14,35,36 The relative cerebral blood flow (rCBF) is then obtained by normalizing the measured BFI with a given baseline.37


Optical instrumentation

Light was injected into the tissue by using a long coherence laser source at 785 nm (120 mW, Crystalaser, Reno, Nevada) through a multimode fiber. The light was collected using eight single-mode fibers and sent into two arrays of four single-photon counting avalanche photodiodes (SPCM-AQ4C, Dumberry, Vaudreuil, Canada) conforming eight (2×4) detectors. The output of each detector was used to build the normalized intensity autocorrelation function by using an eight-channel correlator (Correlator.com, New Jersey).


Diffuse correlation spectroscopy probe

The DCS probe consisted of eight detectors and two source fibers (Fig. 2), with three source–detector separations 2.6, 5.5, and 9.1 mm, respectively. The probe was placed in a holder that was designed to be merged with the manipulator arm of the stereotaxic frame, ensuring a sufficiently precise positioning and good contact between the head of the rat and the transcutaneous probe. External landmarks (eyes, ears, and nose) were used for the orientation and the placement of the probe since the skull and scalp were intact, the Bregma was not visible and could not be used as a reference point.

Fig. 2

(a) A schematic of the transcutaneous probe and its design with the source–detector separations and (b) its placement on the rat head.




The body temperature was kept in the physiological range of (37.4±0.2)°C under the isoflurane and at (37.1±0.6)°C with medetomidine while the respiratory rate decreased during the experimental procedure when using medetomidine as a sedative agent. The initial respiratory rate, which, under isoflurane, was (87±21)  bpm dropped to (55±12)  bpm after shifting to medetomidine sedation and remained stable for the next 1.5 to 2 h increasing gradually thereafter. This temporal window was considered as the stable period to carry out the stimulation experiments.


Functional Magnetic Resonance Imaging

All rats (n=11) were initially subjected to fMRI and showed neuronal activation in the primary somatosensory cortex during unilateral electrical forepaw stimulation.

In Fig. 3(a), we show an example of the activation map of one animal stimulated in the right forepaw. The color map shows the activation on the cortical somatosensory region (SIfl). Averaged BOLD signal increase (%), of the 11 animals over all six experiments (three left and three right forepaw stimulation), is shown in Fig. 3(b), resulting in a mean BOLD signal increase of (4.3±1.1)%, where the profile of five consecutive time-series of activation periods composed in each experiment can be observed. In Table 1, we show data from all the animals.

Table 1

Relative rCBF from all animals at 6 months of age (week 0, week 2, and week 4) and at 18 months of age (week 52 and week 54) presented as ΔrCBF (%) (mean±SEM). fMRI results are presented as [ΔSIBOLD (%)].

AnimalΔSI  BOLD  (%)ΔrCBF  (%)
Week 0Week 0Week 2Week 4Grouped (6 months)Week 52Week 54Grouped (18 months)
Rat 14.3±1.322.8±3.221.9±4.820.2±3.621.6±2.918.4±3.616.3±3.717.3±3.1
Rat 23.9±0.817.0±4.316.1±4.915.9±3.416.2±3.516.1±3.314.2±3.715.1±2.9
Rat 34.5±2.027.8±4.424.7±2.625.9±5.526.1±4.123.6±4.221.0±3.422.2±1.7
Rat 44.0±0.820.4±3.920.9±2.823.4±4.721.6±1.521.2±4.716.9±3.019.0±2.2
Rat 54.5±1.222.4±3.827.9±5.826.7±6.125.6±5.225.8±5.924.6±5.425.2±5.6
Rat 64.2±1.222.4±4.616.3±3.824.7±5.621.1±3.719.5±2.817.1±3.718.2±2.1
Rat 74.1±0.420.3±4.016.7±3.417.2±3.718.0±3.219.3±3.917.4±2.618.3±3.3
Rat 84.6±0.725.6±4.220.4±3.823.7±5.023.2±1.923.2±4.419.1±3.421.1±1.8

Fig. 3

(a) BOLD fMRI activation map during right forepaw stimulation overlapped on the anatomical T2-w images. The activated region corresponds to the forelimb area of the somatosensory cortex (SIfl). (b) BOLD signal increase averaged from all rats (six experiments each). Change in signal increase was (4.2±1.6)% [(mean±SEM)%].



Diffuse Correlation Spectroscopy

After the administration of the medetomidine bolus, and the discontinuation of isoflurane, the corneal and tail-pinch reflexes were absent, except in two animals, in which medetomidine did not cause a sedative effect, probably due to a high level of stress, and, therefore, were excluded. Furthermore, one animal died due to respiratory difficulties. After excluding these animals, a total of eight rats completed all the five sessions of DCS measurements.

In each stimulation experiment, the strongest signal was obtained from the second pair (L-2 or R-2) except for one animal, where the signal was better in L-3 and R-3 source–detector pairs. In Fig. 4, simultaneously measured rCBF change from both hemispheres is shown for a representative animal. The cortex contralateral to the stimulus showed changes in rCBF, whereas the ipsilateral cortex results remained unaltered.

Fig. 4

The stimulation pattern from the cortex: (a) contralateral to the stimulus and (b) ipsilateral to the stimulus. The plots correspond to the source/detector pair L-2 and R-2, respectively. The profile of five consecutive trains of stimuli is evident in the contralateral hemisphere, whereas the ipsilateral hemisphere shows no significant changes in CBF.


Figure 5 shows five different sessions from a representative animal, three sessions at the age of 6 months (week 0, week 2, and week 4), and two at the age of 18 months (week 52 and week 54). The study was performed with a resting period of 2 weeks between each experiment cluster. Each graph shows changes in rCBF averaged from six consecutive experiments (three right and three left forepaw stimulation), from L-2 and R-2 source–detector pairs and are represented as mean±SEM. Relative CBF was (22.4±3.8)%; (27.9±5.8)%; and (26.7±6.1)% for week 0, week 2, and week 4, respectively, at the age of 6 months, and (25.8±5.9)%; (24.6±5.4)% for the last sessions at the age of 18 months. Statistical analysis [one-way analysis of variance (ANOVA), p<0.05] showed no significant difference between the sessions from the same animal. In Table 1, we show data for individual weeks and animals, averaged over all six consecutive experiments. The variability in the fMRI signal is 2% to 4%38 and in the optics, the variation in the signal is around 10% in accordance with the literature.26,39

Fig. 5

The complete set of experimental sessions from a representative animal. Each plot is composed of six experiments of forepaw stimulation (three right forepaw and three left forepaw). (a) Sessions at the age of 6 month (adulthood, week 0), with a rest period of 2 weeks between sessions (weeks 2 and 4), and (b) sessions at 18 month (elderly, weeks 52 and 54). Changes in CBF measured by DCS did not show statistical significant differences over time (one-way ANOVA; p<0.05).


In order to look at the reproducibility of results in the long-term study, the signal was averaged over L-2 and R-2 source–detector pairs (except for one rat, where the signal was stronger for L-3 and R-3), followed by the intrasubject average for all rats and grouped by age. These results give a similar profile of activation at adulthood and after aging (Fig. 6). The averaged mean of ΔrCBF in the 6 months group was (21.7±3.2)% and 1 year later was (19.5±2.8)% in the same animals. Statistical analysis did not show significant differences between groups (p<0.05).

Fig. 6

Changes in rCBF (averaged from eight rats), induced by functional brain activity, measured in adulthood (6 months) and in old age (18 months). The percentage change in rCBF did not show significant difference (unpaired t-test; p<0.05).


Linear regression (Fig. 7) shows a positive correlation with the squared Pearson’s correlation coefficient (r2) of 0.7 between the BOLD signal and ΔrCBF at the age of 6 months [Fig. 7(a)] and 0.6 at the age of 18 months [Fig. 7(b)] with the p<0.05 for both cases. The correlation between ΔrCBF at 6 and 18 months is shown in Fig. 7(c) with r2=0.8 and p<0.05.

Fig. 7

Linear correlation analysis between BOLD and DCS measures as well as between DCS measures at two time-points. (a) ΔrCBF at 6 months is positively correlated to changes in BOLD signal intensity (Pearson’s correlation r2=0.7, p<0.05). (b) ΔrCBF at 18 months is positively correlated to BOLD (%SI) with Pearson’s correlation of r2=0.6, and p<0.05). (c) Correlation between ΔrCBF (%) at 6 and 18 months (r2=0.8, p<0.05).




The present work demonstrates that DCS is a useful tool to study functional brain activity longitudinally, from a few weeks up to a year, in a noninvasive manner in rats. The data obtained by the transcutaneous DCS probe were reliable and reproducible. The results were comparable to fMRI.

Although DCS has been used in the past to measure tissue blood flow and its alterations in animal models of brain-related pathologies, such as stroke32,37 or postnatal hypoxia,40 it has not been used to monitor changes in blood flow longitudinally.

To cross-validate our system, we have compared DCS measurements with fMRI, a well-characterized technique to measure brain activity. DCS measures changes of CBF in the microvasculature,14 while fMRI (BOLD) is sensitive to local changes in deoxyhemoglobin concentration produced by the hemodynamic response to a stimulus.41,42 Although DCS and fMRI measure different parameters, both provide similar information about hemodynamic activity reflecting the neurovascular coupling, and are a good indicator of the hemodynamic response promoted by neuronal activity in noninvasive manner.

Most of the methods in optics measuring CBF, such as LDF21,22,43, LSI,23,24 or OISI,26 require the removal of the scalp or/and thinning of the skull to avoid the possibility of masking effects from external tissues on the measurements, or to increase the SNR. All these surgical interventions carry a risk of infection and/or inflammation, which can cause changes in the vasculature of the extracranial tissues as well as the formation of the scar tissue. It can also lead to localized inflammation of the area of the brain closest to the thinned skull because of the heating generated by the drill.44 Any of these reactions can cause changes locally in the CBF. In addition, one of the targets of the small animal experiments is to obtain data comparable to those obtained in patients with whom it is generally not possible to use invasive probes.

DCS has been validated in different animal models to measure rCBF27,45 under different anesthesia protocols,46 but to our knowledge, no studies have been done in rats with the intact scalp and skull to perform longitudinal studies over months.

Hypoxia, stroke, some neurodegenerative diseases, among other pathologies, cause chronic alterations in brain function. Thus, experimental animal models that mimic such pathologies need to be studied in a longitudinal manner, from days, weeks up to years. In addition, many therapeutic strategies can present their effects in a gradual way, promoting small functional improvements that require a follow-up over an extended time-period. Furthermore, longitudinal studies allow tracking of individual subjects over time, rather than relying on cross-sectional studies with multiple individuals at different time points. In this work, we have shown results on eight rats measured repeatedly at the adult age (6 months old), and the same animals were measured 1 year later (18 months old).

Our experimental paradigm consisted of electrical forepaw stimulation performed under sedation. The optimal neuronal response to nonpainful forepaw electrical stimulation was detected at a stimulating frequency of 6 Hz and is in agreement with previous publications,47 thus, DCS experiments were performed at the same frequency, yielding an easily detectable signal. Both BOLD and rCBF results agree with published results in the literature.44,48,49

In experimental studies with animals, it is essential to use anesthetics or sedatives that immobilize the animal to avoid motion artifacts, reduce stress, and to increase reproducibility.50,51 The most frequently used inhalatory anesthetics agents (halothane, isoflurane) can interfere with normal neural activity, and can affect the neurovascular coupling.44,50,52,53 Other anesthetics that preserve such activity often require the use of paralyzing agents, with the need for intubation and mechanical ventilation of the rat, and often with certain degree of toxicity, thus making longitudinal studies impossible. In the last decade, different methods for anesthesia/sedation have been implemented that allow performing functional experiments in a noninvasive and longitudinal manner.44,48,54

In the present work, a widely validated sedation system for fMRI has been used.44,48,49 Medetomidine, an α2-adrenoreceptor agonist, produces analgesia, sedation, muscle relaxation, and anxiolysis. However, this sedative has some drawbacks to be considered: it promotes a decrease in the respiratory rate, which may alter the basal levels of CBF and CBV,44 and an increase in variability among different individuals due to the depth of sedation. Stress level at the onset of the experiment could make the achievement of an optimal sedation more difficult. Therefore, some animals were excluded from the study. However, all those individuals who achieved satisfactory sedation showed significant changes in CBF evoked by forepaw stimulation. New anesthetics, such as those used recently in fMRI55 that minimally alter physiological parameters or neurovascular coupling, could be used in future studies.

The depth penetration in DCS depends on the source–detector separation on the tissue surface. Our probe was designed to detect CBF in the rat cortex at roughly three different depths for both hemispheres [Fig. 2(b)]. The best signal was achieved with the second source–detector pair combination (5.5 mm), which probes a depth of around 2 mm from the surface, except in one session, where the best response was with the third source–detector pair (9.1 mm) for one animal. If we consider the thickness of the scalp and skull, the actual depth of scanning corresponded to the cortex, in agreement with the activation zone in the somatosensory cortex detected by fMRI. The number of source–detector pairs in our probe was limited to improve the temporal resolution. In the future, probes could be implemented by increasing the number of sources and detectors, encompassing a greater area and depth in the cerebral cortex, allowing the composition of accurate activation maps.14,28,31,56,57

We have shown the reproducibility of the results, with little variability when comparing the same individual at different time points, or when we compare between different individuals. Despite this, a certain degree of variability in CBF measurements should be mentioned. The contribution of extracerebral tissues to the measured CBF over time should be considered. Aging increases the skin thickness and the fat content in soft tissues throughout the duration of the experiment (1 year). Alterations in the mechanism of neurovascular coupling can take place with aging, due to effects, such as the thickening of the vessel wall and a decrease in elasticity.58 Other factors related to age (hypertension, diabetes) have to be taken into consideration when performing longitudinal studies regardless of the pathology of interest. In the future, we believe that with the transcutaneous DCS, it will be possible to study aging and neurovascular coupling in a noninvasive manner in elderly animals.

Furthermore, the case of a slight offset in placement with the transcutaneous probe from the correct place may cause differences in the measured CBF values, since electrical forepaw stimulation elicits a discrete hemodynamic response in a small region of the cerebral cortex. Therefore, the accurate placement of the probe is an issue to consider. External landmarks, such as the position of the eyes, ears, and nose, have been used to localize Bregma as a reference point for the probe placement. However, these landmarks may vary slightly among individuals, as well as with head growing with aging, which may add a component of variability to the measurements.

There is still space for more improvements. We note here that we have followed the literature and used for absorption 0.1  cm1 and for reduced scattering coefficients 15  cm1.30 On the other hand, it has been shown previously in the literature59 that errors in baseline optical parameters minimally affect high relative changes in CBF especially since the reduced scattering coefficient does not change significantly due to functional stimuli. Furthermore, the current study was a proof-of-principle to show that we can see the changes in the CBF over a long time (1 year). It is possible to use more complex systems (tomography) or models, but those would require different protocols and probes. In the future, measuring changes with a new method SCOT28 would give advantages by measuring the CBF through the intact scalp with a large number of detectors.

The used algorithm assumes the homogeneous background of the medium, which may be problematic in the case of high absorptions and for small source–detector separation. Higher order approximations or Monte Carlo models can provide more accurate estimation of the light transport in the tissue,60 introducing the heterogeneities for the upper layers, such as the skull.



The present work demonstrates the ability of DCS to detect neuronal activity through the intact scalp of the rat in a longitudinal manner (over months). This methodology could be useful for the study, in animal models, of the alterations in neuronal activity, which can be compromised in several cerebral pathologies.


ICFO has equity ownership in the spin-off company HemoPhotonics S.L. that commercializes DCS technology. Potential financial conflicts of interest and objectivity of research have been monitored by ICFO’s Knowledge & Technology Transfer Department. No financial conflicts of interest were identified or declared by the authors.


The project was funded by Fundació CELLEX Barcelona, Ministerio de Economía y Competitividad (PHOTODEMENTIA, DPI2015-64358-C2-1-R and DPI2015-64358-C2-2-R, PHOTOSTROKE), the “Severo Ochoa” Programme for Centres of Excellence in R&D (SEV-2015-0522), the Obra social “la Caixa” Foundation (LlumMedBcn), Instituto de Salud Carlos III / FEDER (Medphotage, DTS16/00087), AGAUR-Generalitat (2014-SGR725, 2014SGR-1555), LaserLab-Europe IV. We are indebted to the Experimental MRI 7T Unit of the IDIBAPS for the technical help.


1. C. Roy and C. Sherrington, “On the regulation of the blood-supply of the brain,” J. Physiol. 11(1–2), 85–158 (1890).JPHYA70022-3751 http://dx.doi.org/10.1113/jphysiol.1890.sp000321 Google Scholar

2. C. Tang et al., “Decreased functional connectivity of homotopic brain regions in chronic stroke patients: a resting state fMRI study,” PLoS One 11(4), e0152875 (2016).POLNCL1932-6203 http://dx.doi.org/10.1371/journal.pone.0152875 Google Scholar

3. C.-M. Lin et al., “Arterial spin labeling perfusion study in the patients with subacute mild traumatic brain injury,” PLoS One 11(2), e0149109 (2016).POLNCL1932-6203 http://dx.doi.org/10.1371/journal.pone.0149109 Google Scholar

4. R. Watson and S. J. Colloby, “Imaging in dementia with Lewy bodies: an overview,” J. Geriatr. Psychiatry Neurol. 29(5), 254–260 (2016). http://dx.doi.org/10.1177/0891988716654984 Google Scholar

5. A. R. Switzer et al., “Longitudinal decrease in blood oxygenation level dependent response in cerebral amyloid angiopathy,” Neuroimage Clin. 11, 461–467 (2016). http://dx.doi.org/10.1016/j.nicl.2016.02.020 Google Scholar

6. H. Lockau et al., “Structural and functional magnetic resonance imaging,” PET Clin. 8(4), 407–430 (2013). http://dx.doi.org/10.1016/j.cpet.2013.08.004 Google Scholar

7. P. Ramos-Cabrer et al., “Stem cell mediation of functional recovery after stroke in the rat,” PLoS One 5(9), e12779 (2010).POLNCL1932-6203 http://dx.doi.org/10.1371/journal.pone.0012779 Google Scholar

8. S. Ogawa et al., “Brain magnetic resonance imaging with contrast dependent on blood oxygenation,” Proc. Natl. Acad. Sci. U. S. A. 87(24), 9868–9872 (1990). http://dx.doi.org/10.1073/pnas.87.24.9868 Google Scholar

9. S.-G. Kim and S. Ogawa, “Biophysical and physiological origins of blood oxygenation level-dependent fMRI signals,” J. Cereb. Blood Flow Metab. 32(7), 1188–1206 (2012). http://dx.doi.org/10.1038/jcbfm.2012.23 Google Scholar

10. M. J. P. Barrett and V. Suresh, “Improving estimates of the cerebral metabolic rate of oxygen from optical imaging data,” Neuroimage 106, 101–110 (2015).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2014.11.041 Google Scholar

11. M. S. Andrew et al., “Temporal comparison of functional brain imaging with diffuse optical tomography and fMRI during rat forepaw stimulation,” Phys. Med. Biol. 48(10), 1391–1403 (2003).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/48/10/311 Google Scholar

12. J. P. Culver et al., “Volumetric diffuse optical tomography of brain activity,” Opt. Lett. 28(21), 2061–2063 (2003).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.28.002061 Google Scholar

13. P. J. van den Berg, K. Daoudi and W. Steenbergen, “Review of photoacoustic flow imaging: Its current state and its promises,” Photoacoustics 3(3), 89–99 (2015). http://dx.doi.org/10.1016/j.pacs.2015.08.001 Google Scholar

14. T. Durduran et al., “Diffuse optics for tissue monitoring and tomography,” Rep. Prog. Phys. 73(7), 76701 (2010). http://dx.doi.org/10.1088/0034-4885/73/7/076701 Google Scholar

15. R. Bonner and R. Nossal, “Model for laser Doppler measurements of blood flow in tissue,” Appl. Opt. 20(12), 2097–2107 (1981).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.20.002097 Google Scholar

16. D. A. Boas and A. K. Dunn, “Laser speckle contrast imaging in biomedical optics,” J. Biomed. Opt. 15(1), 011109 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3285504 Google Scholar

17. D. A. Boas and A. G. Yodh, “Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation,” J. Opt. Soc. Am. A 14(1), 192 (1997).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.14.000192 Google Scholar

18. A. K. Dunn, “Laser speckle contrast imaging of cerebral blood flow,” Ann. Biomed. Eng. 40(2), 367–377 (2012).ABMECF0090-6964 http://dx.doi.org/10.1007/s10439-011-0469-0 Google Scholar

19. A. Borogovac and I. Asllani, “Arterial spin labeling (ASL) fMRI: advantages, theoretical constrains and experimental challenges in neurosciences,” Int. J. Biomed. Imaging 2012, 818456 (2012). http://dx.doi.org/10.1155/2012/818456 Google Scholar

20. J. Mayhew et al., “Increased oxygen consumption following activation of brain: theoretical footnotes using spectroscopic data from barrel cortex,” Neuroimage 13(6), 975–987 (2001).NEIMEF1053-8119 http://dx.doi.org/10.1006/nimg.2001.0807 Google Scholar

21. H. Piilgaard and M. Lauritzen, “Persistent increase in oxygen consumption and impaired neurovascular coupling after spreading depression in rat neocortex,” J. Cereb. Blood Flow Metab. 29(9), 1517–1527 (2009). http://dx.doi.org/10.1038/jcbfm.2009.73 Google Scholar

22. U. Dirnagl et al., “Continuous measurement of cerebral cortical blood flow by laser-Doppler flowmetry in a rat stroke model,” J. Cereb. Blood Flow Metab. 9(5), 589–596 (1989). http://dx.doi.org/10.1038/jcbfm.1989.84 Google Scholar

23. A. B. Parthasarathy, S. M. S. Kazmi and A. K. Dunn, “Quantitative imaging of ischemic stroke through thinned skull in mice with multi exposure speckle imaging,” Biomed. Opt. Express 1(1), 246–259 (2010).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.1.000246 Google Scholar

24. S. M. S. Kazmi et al., “Chronic imaging of cortical blood flow using multi-exposure speckle imaging,” J. Cereb. Blood Flow Metab. 33(6), 798–808 (2013). http://dx.doi.org/10.1038/jcbfm.2013.57 Google Scholar

25. T. Durduran et al., “Spatiotemporal quantification of cerebral blood flow during functional activation in rat somatosensory cortex using laser-speckle flowmetry,” J. Cereb. Blood Flow Metab. 24, 518–525 (2004). http://dx.doi.org/10.1097/00004647-200405000-00005 Google Scholar

26. W. B. Baker et al., “Neurovascular coupling varies with level of global cerebral ischemia in a rat model,” J. Cereb. Blood Flow Metab. 33(1), 97–105 (2012). http://dx.doi.org/10.1038/jcbfm.2012.137 Google Scholar

27. C. Cheung et al., “In vivo cerebrovascular measurement combining diffuse near-infrared absorption and correlation spectroscopies,” Phys. Med. Biol. 46(8), 2053–2065 (2001).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/46/8/302 Google Scholar

28. T. Dragojevic et al., “High-density speckle contrast optical tomography (SCOT) for three dimensional tomographic imaging of the small animal brain,” Neuroimage 153, 283–292 (2017).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2017.04.003 Google Scholar

29. H. Tanaka, T. Katura and H. Sato, “Task-related oxygenation and cerebral blood volume changes estimated from NIRS signals in motor and cognitive tasks,” Neuroimage 94, 107–119 (2014).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2014.02.036 Google Scholar

30. C. Zhou et al., “Diffuse optical correlation tomography of cerebral blood flow during cortical spreading depression in rat brain,” Opt. Express 14(3), 1125–1144 (2006).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.14.001125 Google Scholar

31. G. Yu et al., “Frequency-domain multiplexing system for in vivo diffuse light measurements of rapid cerebral hemodynamics,” Appl. Opt. 42(16), 2931–2939 (2003).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.42.002931 Google Scholar

32. Y. Shang et al., “Diffuse optical monitoring of repeated cerebral ischemia in mice,” Opt. Express 19(21), 20301–20315 (2011).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.19.020301 Google Scholar

33. J. Strupp, “Stimulate: a GUI based fMRI analysis software package,” Neuroimage 3(3), S607 (1996).NEIMEF1053-8119 http://dx.doi.org/10.1016/S1053-8119(96)80609-4 Google Scholar

34. D. A. Boas and Y. Campbell, “Scattering and Imaging with diffusing temporal field correlations,” Phys. Rev. Lett. 75, 1855 (1995).PRLTAO0031-9007 http://dx.doi.org/10.1103/PhysRevLett.75.1855 Google Scholar

35. D. A. Boas et al., “Establishing the diffuse correlation spectroscopy signal relationship with blood flow,” Neurophotonics 3(3), 031412 (2016). http://dx.doi.org/10.1117/1.NPh.3.3.031412 Google Scholar

36. R. C. Mesquita et al., “Direct measurement of tissue blood flow and metabolism with diffuse optics,” Philos. Trans. A Math. Phys. Eng. Sci. 369(1955), 4390–4406 (2011). http://dx.doi.org/10.1098/rsta.2011.0232 Google Scholar

37. M. G. Burnett et al., “Electrical forepaw stimulation during reversible forebrain ischemia decreases infarct volume,” Stroke 37(5), 1327–1331 (2006).SJCCA70039-2499 http://dx.doi.org/10.1161/01.STR.0000217305.82123.d8 Google Scholar

38. K. A. Williams et al., “Comparison of α-chloralose, medetomidine and isoflurane anesthesia for functional connectivity mapping in the rat,” Magn. Reson. Imaging 28(7), 995–1003 (2010).MRIMDQ0730-725X http://dx.doi.org/10.1016/j.mri.2010.03.007 Google Scholar

39. C. Ayata et al., “Laser speckle flowmetry for the study of cerebrovascular physiology in normal and ischemic mouse cortex,” J. Cereb. Blood Flow Metab. 24, 744–755 (2004). http://dx.doi.org/10.1097/01.WCB.0000122745.72175.D5 Google Scholar

40. E. M. Buckley et al., “In vivo monitoring of cerebral hemodynamics in the immature rat: effects of hypoxia-ischemia and hypothermia,” Dev. Neurosci. 37(4–5), 407–416 (2015). Google Scholar

41. S. Ogawa et al., “Intrinsic signal changes accompanying sensory stimulation: functional brain mapping with magnetic resonance imaging,” Proc. Natl. Acad. Sci. U. S. A. 89(13), 5951–5955 (1992). http://dx.doi.org/10.1073/pnas.89.13.5951 Google Scholar

42. A. C. Silva et al., “Simultaneous blood oxygenation level-dependent and cerebral blood flow functional magnetic resonance imaging during forepaw stimulation in the rat,” J. Cereb. Blood Flow Metab. 19(8), 871–879 (1999). http://dx.doi.org/10.1097/00004647-199908000-00006 Google Scholar

43. J. A. Detre et al., “Signal averaged laser Doppler measurements of activation-flow coupling in the rat forepaw somatosensory cortex,” Brain Res. 796(1–2), 91–98 (1998).BRREAP0006-8993 http://dx.doi.org/10.1016/S0006-8993(98)00322-9 Google Scholar

44. M. Fukuda et al., “Effects of the a2-adrenergic receptor agonist dexmedetomidine on neurovascular responses in somatosensory cortex,” J. Neurosci. 37(1), 80–95 (2013). http://dx.doi.org/10.1111/ejn.12024 Google Scholar

45. K. Verdecchia et al., “Quantifying the cerebral metabolic rate of oxygen by combining diffuse correlation spectroscopy and time-resolved near-infrared spectroscopy,” J. Biomed. Opt. 18(2), 027007 (2013).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.18.2.027007 Google Scholar

46. M. A. Franceschini et al., “The effect of different anesthetisc on neurovascular coupling,” Neuroimage 51(4), 1367–1377 (2010).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2010.03.060 Google Scholar

47. N. Van Camp, M. Verhoye and A. Van der Linden, “Stimulation of the rat somatosensory cortex at different frequencies and pulse widths,” NMR Biomed. 19(1), 10–17 (2006). http://dx.doi.org/10.1002/(ISSN)1099-1492 Google Scholar

48. R. Weber et al., “A fully noninvasive and robust experimental protocol for longitudinal fMRI studies in the rat,” Neuroimage 29(4), 1303–1310 (2006).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2005.08.028 Google Scholar

49. P. Ramos-Cabrer et al., “Continuous noninvasive monitoring of transcutaneous blood gases for a stable and persistent BOLD contrast in fMRI studies in the rat,” NMR Biomed. 18(7), 440–446 (2005). http://dx.doi.org/10.1002/(ISSN)1099-1492 Google Scholar

50. V. C. Austin et al., “Confounding effects of anesthesia on functional activation in rodent brain: a study of halothane and α-chloralose anesthesia,” Neuroimage 24(1), 92–100 (2005).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2004.08.011 Google Scholar

51. R. R. Peeters et al., “Comparing BOLD fMRI signal changes in the awake and anesthetized rat during electrical forepaw stimulation,” Magn. Reson. Imaging 19(6), 821–826 (2001).MRIMDQ0730-725X http://dx.doi.org/10.1016/S0730-725X(01)00391-5 Google Scholar

52. A. Schroeter et al., “Specificity of stimulus-evoked fMRI responses in the mouse: the influence of systemic physiological changes associated with innocuous stimulation under four different anesthetics,” Neuroimage 94, 372–384 (2014).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2014.01.046 Google Scholar

53. K. Masamoto and I. Kanno, “Anesthesia and the quantitative evaluation of neurovascular coupling,” J. Cereb. Blood Flow Metab. 32(7), 1233–1247 (2012). http://dx.doi.org/10.1038/jcbfm.2012.50 Google Scholar

54. C. P. Pawela et al., “A protocol for use of medetomidine anesthesia in rats for extended studies using task-induced BOLD contrast and resting-state functional connectivity,” Neuroimage 46(4), 1137–1147 (2009).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2009.03.004 Google Scholar

55. M. M. Petrinovic et al., “A novel anesthesia regime enables neurofunctional studies and imaging genetics across mouse strains,” Sci. Rep. 6, 24523 (2016).SRCEC32045-2322 http://dx.doi.org/10.1038/srep24523 Google Scholar

56. Z.-J. Lin et al., “Interleaved imaging of cerebral hemodynamics and blood flow index to monitor ischemic stroke and treatment in rat by volumetric diffuse optical tomography,” Neuroimage 85, 566–582 (2014).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2013.07.020 Google Scholar

57. Z. J. Lin et al., “Volumetric diffuse optical tomography for small animals using a CCD-camera-based imaging system,” Int. J. Opt. 2012, 276367 (2012). http://dx.doi.org/10.1155/2012/276367 Google Scholar

58. A. Furuta et al., “Medullary arteries in aging and dementia,” Stroke 22(4), 442–446 (1991).SJCCA70039-2499 http://dx.doi.org/10.1161/01.STR.22.4.442 Google Scholar

59. W. B. Baker et al., “Modified Beer-Lambert law for blood flow,” Biomed. Opt. Express 5(11), 4053–4075 (2014).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.5.004053 Google Scholar

60. M. A. O’Leary, “Imaging with diffuse photon density waves,” Doctoral dissertation, University of Pennsylvania (1996). Google Scholar


Igor Blanco received his PhD from ICFO—The Institute of Photonic Sciences, Spain, researching within the medical optics group, the cerebral hemodynamics in experimental and clinical neurology, under the supervision of Prof. Durduran. He is currently leading a data science team in the private consulting sector.

Peyman Zirak obtained his PhD in the laboratory of Prof. Alfons Penzkofer at University of Regensburg, Germany. As postdoctoral researcher in the laboratories of Prof. Turgut Durduran at ICFO–The Institute of Photonic Sciences, Spain, he worked on (pre-) clinical application of diffuse optical modalities. Currently, as a research fellow of Prof. Andreas Zumbusch at University of Konstanz, Germany, he works on nonlinear optical microendoscopy and waveform shaping to access deeper tissue while imaging.

Tanja Dragojević is a doctoral student at ICFO—The Institue of Photonic Sciences, Spain, with the medical optics group. She received her master’s degree in medical physics at the University of Novi Sad, Serbia. Her research is focused on developing speckle-based intrumentation and protocols, for preclinical and clinical studies.

Clara Castellvi is working in IIBB-CSCI, Barcelona, on MRI and image analysis. She obtained her bachelor’s degree in biology at the University of Autonoma, Barcelona. Her master’s degree was in bioinformatics and biostatistics from Open University of Catalonia.

Turgut Durduran was trained at the University of Pennsylvania. In 2009, he moved to ICFO—The Institute of Photonic Sciences, Spain, where he leads the medical optics group. His research interests revolve around the use of diffuse light to noninvasively probe tissue function. The group develops new technologies and algorithms and routinely translates them to preclinical, clinical, and industrial applications.

Carles Justicia received his PhD in biology from the University of Barcelona in 2001. He specialized in MR imaging at the Max Planck Institute for neurological research in Cologne under the supervision of Prof. Mathias Hoehn. He is currently a senior researcher at the Biomedical Research Institute in Barcelona (IIBB-CSIC) and is part of the group of “Brain ischemia: clinical and experimental studies” in IDIBAPS, where he investigates vascular changes associated with stroke and aging.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Igor Blanco, Peyman Zirak, Tanja Dragojević, Clara Castellvi, Turgut Durduran, Carles Justicia, "Longitudinal, transcranial measurement of functional activation in the rat brain by diffuse correlation spectroscopy," Neurophotonics 4(4), 045006 (5 December 2017). http://dx.doi.org/10.1117/1.NPh.4.4.045006 Submission: Received 10 August 2017; Accepted 9 November 2017
Submission: Received 10 August 2017; Accepted 9 November 2017

Functional magnetic resonance imaging





Blood circulation


Back to Top