Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Optical brain imaging in vivo : techniques and applications from animal to man

Optical brain imaging in vivo : techniques and applications from animal to man other imaging modalities, optical contrast is often a natural feature of many small, biocompatible molecules. Additional advantages of optical imaging include its use of nonionizing radiation, and that instruments are generally inexpensive compared to other clinical imaging modalities. The major challenge of optical imaging is to overcome the effects of light scattering, which limits penetration depth and achievable imaging resolution. “Optical imaging” therefore encompasses a very wide range of measurement techniques, each using different methods to overcome the effects of scatter. These range from laser scanning microscopy of submicron structures, to diffuse optical tomography of large volumes of tissue. Functional brain imaging has also seen dramatic growth in recent years, impacted significantly by the availability of functional magnetic resonance imaging fMRI since 1990.6–9 As a result, new applications of functional brain imaging in humans and animals are continually being established, including neurosurgical planning, the investigation of the physiological basis of neurological diseases such as epilepsy, Alzheimer’s, and stroke, the development of diagnostic methods, drugs, treatments, and interventions, and the study of cognitive and perceptual responses and developmental changes.9 Positron emission tomography PET and x-ray 1083-3668/2007/12 5 /051402/28/$25.00 © 2007 SPIE Journal of Biomedical Optics 051402-1 September/October 2007 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… computed tomography x-ray CT are also valuable and established tools for brain imaging. The obvious advantage of optical imaging over other modalities is reduced cost and infrastructure requirements such as shielded rooms and synchrotrons . However, a much more important distinction is that optical imaging offers such a broad range of contrast mechanisms. While fMRI, PET, and x-ray CT are more established clinical modalities, all suffer from limited functional contrast. For example, to visualize regions of the brain exhibiting functional changes in response to stimuli, fMRI typically exploits the blood oxygen level dependent BOLD signal, widely thought to correlate with HbR concentration. However, a BOLD increase an HbR decrease may correspond to an increase in oxygenation, or a decrease in blood volume. Newer fMRI methods such as arterial spin labeling ASL can provide measures of flow, and intravenous contrast agents such as monocrystalline iron oxide nanocompounds MION can allow cerebral blood volume CBV changes to be observed.10,11 However, these measurements cannot easily be acquired simultaneously, and only the combination of these methods can provide insight into the true hemodynamic responses and oxygenation dynamics in the brain.12 PET is performed using contrast agents such as 2-fluoro2-deoxy-D-glucose containing radioactive isotopes of carbon, nitrogen, fluorine, or oxygen. Oxygen and glucose utilization in the brain can be imaged via the localization of these isotopes. However, resolution is typically poor, data acquisition is slow and hence dynamics cannot be readily evaluated , contrast agents must be manufactured locally prior to each scan, and radioactive dose limits its use in certain patient populations, often including infants. X-ray CTs as well as structural MRI have excellent spatial resolution and can be used to detect morphological changes in the brain, for example following stroke or to delineate tumors. However, CT typically cannot elicit information about the function of the tissue, and so potentially reversible changes preceding structural changes cannot be detected. Therefore, the real advantage of optical imaging over these techniques is its ability to detect and image a wide range of functional parameters. For clinical applications, noninvasive optical imaging can provide complimentary information to other modalities such as fMRI and provide a low-cost alternative in some cases. This is in addition to serving populations often unable to receive MRI or PET scans such as young infants or the critically ill. Clinical optical brain imaging is generally noninvasive and uses near-infrared NIR light to obtain improved penetration through the scalp, skull, and brain see Fig. 1 a . In addition, optical brain imaging is finding widespread applications as a research tool for both clinical and animal studies of basic brain function and diseases. Still, so little is known about the way that the normal brain functions, in part due to the difficulties of measuring such a complex organ without disturbing or damaging its in-vivo function. Optical imaging allows the living brain to be closely observed, and many functional interactions and changes to be investigated over many length scales. The effects of pathologies and treatments can also be readily investigated using fairly basic optical imaging tools. In animals, it is also relatively simple to expose the cortex, allowing high-resolution direct optical imJournal of Biomedical Optics aging of the brain’s surface. Extrinsic dyes and visible wavelengths can also be used see Fig. 1 b . In summary, optical imaging’s key advantage is the ability to measure a range of functional contrast, can readily be exploited in functional brain imaging via a wide range of approaches, from animal studies of the intricate cellular mechanisms of normal and diseased brain, to in-vivo noninvasive clinical brain imaging. In this review, we describe a variety of optical imaging approaches to functional brain imaging. The basic principles of each technique are summarized, followed by examples of use for in-vivo brain imaging. Section 2 details the use of optical techniques for exposed cortex imaging. Section 3 describes noninvasive clinical optical imaging of the living brain. Exposed-Cortex Imaging and Microscopy The use of optical imaging in animal studies is widespread. Animal imaging has become more than just a step toward clinical diagnostic imaging, but a tool with which to learn more about the basic mechanisms of brain function as well as disease processes and the effects of treatments. These findings can be utilized in their own right to aid in the development of drugs and treatment. They can also contribute to the interpretation of results from other imaging modalities such as fMRI. Recent applications of animal imaging have included studies of stroke acute and chronic ,13 Alzheimer’s disease,14 migraine,15 epilepsy,16 the effects of cocaine,17 and the mechanisms of neurovascular coupling.18,19 Imaging animals rather than humans provides significantly more flexibility, since preparations can be much more controlled and diseases and treatments can be systematically compared. Extrinsic dyes and cross-validation techniques such as electrophysiology can be utilized and developed without the need for clinical regulatory approval. The adverse effects of light scattering worsen as the size of the tissue being interrogated increases. Imaging smaller animals therefore also offers significant technical advantages, allowing higher resolution imaging and improved sensitivity and quantitation. For very high resolution imaging, the cortex can be surgically exposed, allowing direct optical imaging of the brain’s surface with only minimal disturbance to normal brain activity. Exposed cortex optical imaging is highly versatile, and most commonly performed in animals, although it has also been performed on the human brain, for example, during epilepsy surgery.20–22 Next, we provide an overview of the most established exposed-cortex, in-vivo optical neuroimaging techniques ranging from 2-D and 3-D macroscopic imaging to microscopic in-vivo imaging using two-photon microscopy. Basic principles and methods are first described, followed by specific examples of neuroimaging applications. 2.1 Macroscopic Optical Imaging of Exposed Cortex in Two Dimensions 2.1.1 Experimental method During surgical preparation for exposed-cortex imaging, the animal commonly mouse, rat, cat, or primate is anesthetized while the scalp is retracted and the skull is gently thinned over the brain area of interest. In some cases, this piece of skull September/October 2007 051402-2 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 1 Intrinsic chromophores and structures of the brain. a For small animal imaging, it is possible to use visible light to look at oxygenationdependent hemoglobin absorption as well as exogenous absorbing and fluorescing dyes . The cortex can be exposed and imaged at high resolution, or imaged through the intact skull and scalp at lower resolution. b For human brain imaging, near-infrared light will penetrate more readily through the scalp and skull to sample the brain thanks to lower scatter and absorption. The human brain has a more complex convoluted structure compared to lower mammals. Spectra: major chromophores in brain are oxy- and deoxyhemoglobin and water lipid is omitted here . For generation of these spectra: water content is assumed to be 90% quoted value in newborn brain, 71 to 85% in adult brain34 . Hemoglobin absorption is shown assuming 2-mM concentration of hemoglobin in blood and 3% content of blood in tissue 60 mM in tissue34 calculated from spectra.257 HbR spectrum represents 3% fully deoxygenated blood in tissue. Scatter spectrum is approximate using: s = A −b, where A = 1.14 10−7 m and b = 1.3.228 may be removed to provide a direct view of an area of the brain’s surface. In mice, the brain can often be visualized without even thinning the skull. Methods have been carefully developed to ensure that normal brain function is not affected by these approaches. Most often, the animal will remain anesthetized during measurements.23 In some cases, however, a permanent window can be surgically implanted and the animals can survive for over a year while an area of their brain is accessible for repeated direct optical imaging.14,24 Imaging of pathologies requires only that the animal be anesthetized. Images during acute interventions, such as adJournal of Biomedical Optics ministration of drugs or introduction of a stroke, must be carefully synchronized with image acquisition. While the brain’s functional response to stimulus is being studied, a carefully controlled stimulus is presented during imaging. Such stimuli may include visual e.g., showing patterns or pictures , somatosensory e.g., electrical or tactile stimulation of the whiskers or forepaw , or auditory stimuli different frequencies and amplitudes . Stimuli may be repeated in regular blocks, or presented in a randomized pattern.18,25 Data are usually shown as the average response to multiple stimulus repetitions to obtain September/October 2007 051402-3 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 2 A typical setup for camera-based imaging of the exposed cortex. Inset shows typical hemodynamic and calcium-sensitive responses in rat somatosensory cortex to 1-mA forepaw stimulus delivered in 3 msec pulses at 3 Hz for 4 sec. better signal to noise where noise includes normal background physiological variations . 2.1.2 Basic principles: two-dimensional camerabased imaging 2-D camera-based imaging of the exposed cortex is the most simple of optical imaging techniques and is often referred to as “optical imaging” in the neuroscience literature. We describe three approaches to 2-D imaging of the cortex: 1. intrinsic signal or hyperspectral imaging of oxy- and deoxyhemoglobin dynamics; 2. imaging of extrinsic voltage sensitive dyes providing sensitivity to neuronal activity ; and 3. speckle flow imaging, capable of imaging the blood flow dynamics in the superficial cortex. Two-dimensional intrinsic and hyperspectral imaging. Figure 2 shows a typical setup for 2-D hyperspectral imaging of exposed rat cortex. This neuroimaging approach relies on the distinctive absorption spectra of oxy- and deoxyhemoglobin HbO2 and HbR as shown in Fig. 1 a . The hemoglobin in blood is the most significant absorber in the brain at visible and NIR wavelengths, and functional activity in the brain results in changes in blood flow, blood volume, and oxygenation.26. Blood delivers oxygen to tissues by locally converting HbO2 to HbR, so changes in the oxygenation state of blood correspond to changes in the relative concentrations of HbO2 and HbR. Consequently, variations in the amount and oxygenation state of hemoglobin modulate the absorption properties of the brain. By simply shining light into the exJournal of Biomedical Optics posed cortex and taking pictures with a camera, the hemodynamic response to stimulus can be readliy observed. Further, since HbR and HbO2 have different absorption spectra see Figs. 1 a and 1 b , measurements with different wavelengths of light can produce images that are preferentially sensitive to changes in the concentration of either HbO2 or HbR. At isobestic points, where the HbR and HbO2 absorption is the same e.g., 500, 530, 570, and 797 nm , changes in total hemoglobin concentration HbT= HbR+ HbO2 can be measured independently of changes in blood oxygenation. Early approaches to optical imaging of the exposed cortex referred to “intrinsic signal” imaging, and often used illumination at an isobestic point to examine localized changes in tissue HbT.27 By adding a second wavelength often 610 nm , researchers were able to investigate oximetric changes in addition to changes in HbT.28,29 Since HbO2 absorption at 610 nm is low around 16% of the absorption of HbR , changes in absorption at 610 nm are largely attributed to, and interpreted as, changes in HbR. However, it should be noted that even though HbO2 absorption at 610 nm is lower than that of HbR, HbO2 changes can interfere with and distort inference of HbR concentration changes from measurements at 610 nm. To calculate the actual changes in concentration of HbO2, HbR, and HbT, the absorption of the brain must be measured at two or more wavelengths hyperspectral imaging . In some cases, multiple wavelengths across the visible spectrum are used. This is typically achieved by synchronizing rapid camera acquisition with a filter wheel in front of a white light September/October 2007 051402-4 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… ln I ln I ln I DPFx = DPFx DPFx 1 HbO2, 1 2 HbO2, 2 DPFx DPFx DPFx 1 HbR2, 1 2 HbR2, 2 n HbO2, n n HbR2, n cHbO2 cHbR I = exp − I0 aDPFx +G , where I is the detected intensity, I0 is the incident intensity, DPFx is the mean pathlength traveled by light in the tissue a function of baseline scattering and absorbing properties , G is a geometry-dependent factor, and a is the absorption coefficient of the tissue as given by Beer’s law: nc n , i.e., the absorption coefficient is a linear sum of the absorption of the n chromophores in the tissue, where n is the wavelength-dependent specific absorption coefficient of a particular chromophore e.g., the HbO2 absorption spectrum and cn is its concentration. Note that “differential pathlength factor” DPF is a term used to describe the approximate multiplicative increase in pathlength experienced by a photon traveling in a scattering medium between two points a linear distance x apart. Thus, the true pathlength is DPFx. I0 may be spatially varying and difficult to accurately determine. G is also very difficult to quantify. Therefore, it is common to consider changes in detected signal relative to an initial time point t = 0 or a reference state: ln It I t0 =− t − a t0 DPFx ln I = − aDPFx. If measurements at two wavelengths 1 and 2 are acquired, it is possible to use these differential measurements to determine the changes in HbO2 and HbR concentration occurring relative to time t = t0 using: ln I ln I =− =− HbO2, 1 cHbO2 + cHbO2 + HbR, 1 cHbR DPFx HbO2, 2 HbR, 2 cHbR DPFx 2 . 4 Where more than two wavelengths are available, HbO2 and HbR can be calculated using a linear model: Journal of Biomedical Optics While imaging absorption changes in superficial tissues appears quite straightforward, the effects of scatter must be carefully considered. Since scattering is wavelength dependent in tissues an approximate scatter spectrum is shown in Fig. 1 a , each measurement wavelength will experience different levels of scattering. The mean path of light through tissue DPFx will also be altered by the background absorption properties of the tissue, which are also heavily wavelength dependent, owing to the presence of hemoglobin, so DPFx 1 DPFx 2 in almost all cases. Therefore, as given by Eqs. 4 and 5 , accurate estimates of these wavelengthdependent pathlength factors DPFx n are required even to calculate differential changes in hemoglobin concentration cHbO2 and cHbR . Note, however, that if I t / I t0 is measured at an isobestic point, and DPFx does not change over time, then I t / I t0 will be directly proportional to HbT t . Estimation of this wavelength-dependent pathlength DPFx is difficult. For exposed-cortex imaging, a Monte Carlo model of light propagation is generally employed.35 Such a model requires estimates of the background absorption and oxygenation level of the brain tissue to predict the likely scattered paths of photons.31 For camera imaging, the mean pathlength calculation should simulate a widely distributed illumination, and incorporate the numerical aperture of the camera’s lens. It has been demonstrated that incorrect estimates of DPFx can directly impact the correct calculation of HbO2, HbR, and HbT changes in tissue.35 Note also that if only values of HbO2 t and HbR t are derived, that it is not possible to infer changes in absolute oxygen saturation SO2 = HbO2 t / HbO2 t + HbR t − HbO2 t0 / HbO2 t0 + HbR t0 without first estimating the baseline absolute concentrations of HbO2 and HbR. Another important issue to be considered in 2-D imaging of the cortex is the depth sensitivity of the measurement. The detected light has scattered within the tissue, and reports absorption that it has experienced along its path. The majority of detected photons have predominantly sampled the very superficial layers of the cortex, whereas fewer have traveled deeply and managed to return to the surface. This means that the image of the cortex is a superficially weighted sum of signals from different layers. Since the cortex has a complex structure of large superficial vessels and deeper capillary beds, the relative contributions of these compartments to the overall signal can be difficult to interpret. Additionally, since the mean penetration depth of the light will also depend on its wavelength, all wavelengths have not necessarily interrogated exactly the same region of tissue, which can impact the quantitative accuracy of estimates of HbO2, HbR, and HbT. September/October 2007 051402-5 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… A typical hemodynamic response in the somatosensory cortex of a rat during a four-second forepaw stimulus 3 s pulses at 3 Hz is shown inset in Fig. 2. The hemodynamic response to stimulus is fairly slow, peaking several seconds after the first stimulus pulse. A large increase in HbO2 is generally seen, along with an increase in HbT and a decrease in HbR. The overall timing characteristics of each hemoglobin type are different from each other, and contain information about the underlying oxygen delivery and consumption dynamics.19,36,37 Some researchers routinely observe a short initial increase in HbR, attributed to rapid oxygen consumption at the site of neuronal activation.29,31,33 The origin of this initial dip is somewhat contentious, and it is discussed further in Sec. 2.1.3. The data shown in Fig. 2 are an average response of 400 stimulus repetitions. It should be noted that individual stimulus responses are often affected by physiological noise, such as variations in baseline blood flow to the region and fluctuations in systemic parameters such as blood pressure. Despite all of these important considerations, 2-D optical imaging of the exposed cortex has provided very significant and well validated observations of functional cortical hemodynamics. Specific examples of applications of these techniques are given at the end of this section. Two-dimensional voltage sensitive dye imaging. As described before, intrinsic signal and hyperspectral imaging can report the hemodynamic response to functional activation in the brain. This hemodynamic response is generally assumed to spatially correlate with the underlying brain activity corresponding to neurons firing and conducting electrical signals. The interrelation between neuronal and hemodynamic activity is currently an area of intense research.18,38,39 This neurovascular coupling is thought to be an essential part of normal brain activity and is believed to be disturbed by pathologies such as Alzheimer’s disease and stroke.40,41 Understanding whether it is possible to determine detailed information about neuronal activity from hemodynamic measures is also an important area of research, driven by the fact that clinical fMRI is currently only sensitive to hemodynamic changes. Voltage sensitive dyes VSDs provide a way to optically detect neuronal activity in the in-vivo exposed-cortex. VSDs were first used in 196842 but received renewed attention in 1986 when they were demonstrated for in-vivo mammalian brain imaging by Grinvald and colleagues43–45 studies that predate the first reported intrinsic signal measurements.27 Since then, significant advances in digital detector technology, dye formulation, and staining strategies have led to established use of VSDs to visualize neuronal activity in-vivo2,24,46,47 and also for mapping electrical activity in the heart.48,49 VSDs are molecules that generally bind across a neuron’s membrane. Changes in membrane potential cause measurable changes in the fluorescence of the dye.50 Since a large percentage of the total membrane area is dendritic, VSD signals reflect both spiking and synaptic activity.42 Calcium sensitive dyes CaSDs are compounds that increase their fluorescence in response to increases in calcium ion concentration. When cells such as neurons are loaded with CaSDs, it becomes possible to monitor intracellular calcium.51 Since a calcium influx occurs when a neuron fires an action Journal of Biomedical Optics potential, the calcium dye’s fluorescence can optically report this, providing a second tool with which to optically image neuronal activity. To image these dyes, a similar system to that shown in Fig. 2 is used, with the filter wheel replaced by an excitation filter of the appropriate wavelength. A long-pass filter is placed in front of the camera to block excitation light and isolate the fluorescence emission. In some cases, a dichroic mirror is used such that the illumination light is reflected onto the brain, and emitted fluorescence passes through the dichroic to the camera via a secondary emission filter to remove residual excitation light . This is the configuration of a commercially available VSD imaging system sold by Optical Imaging Incorporated Rehovot, Israel .2 As shown in the inset of Fig. 2, the hemodynamic response to stimulus typically peaks 2 to 4 sec after the stimulus begins. The neuronal response is much faster. Measured using VSDs and CaSDs, neuronal responses typically peak within 25 to 100 msec, and are usually followed by an inhibitory or refractory period of negative signal, which returns to baseline at around 150 to 300 msec after the stimulus starts.52 The black trace in Fig. 2 shows the fluorescence changes seen in the somatosensory cortex stained with Oregon Green 488 BAPTA-1, AM calcium sensitive dye Invitrogen during the 4-sec forepaw stimulus 3- s pulses at 3 Hz . Very fast, individual responses for each stimulation pulse can be seen, in contrast to the very slow hemodynamic response. To capture these responses, a faster camera is required than for hyperspectral or intrinsic signal imaging preferably 10 ms per frame . Voltage sensitive dyes typically have quite low contrast typical F / F 0.5% . As a result, they can be affected by changes in hemoglobin absorption due to heart rate, breathing, and of course functional activation. Two methods have been developed to overcome this. 1. New dyes have been developed that excite and emit at higher wavelengths, corresponding to regions where hemoglobin absorption is less significant see Fig. 1 a . 2. Acquisition of data is often triggered on the heart beat and sometimes breathing of the animal.2 Many repetitions are generally averaged, and “blank” acquisitions where no stimuli are presented are interleaved such that the final signal is given by the stimulated – blank data. Voltage sensitive dyes can only be used in animals at present, and are usually topically applied directly to the exposed cortex for 1 to 2 h before imaging commences.53 Some VSDs have been shown to quickly photobleach or have phototoxic effects, and so illumination times are typically kept to a minimum. It is also important to note that the change in fluorescence or absorption of VSDs is usually a manifestation of a spectral shift in either the absorption or emission spectrum of the dye.54 This means that measurement wavelength selection should not necessarily be based on the peak absorption and emission wavelengths, nor should broadband emission filters spanning the whole emission peak always be selected. The excitation and emission wavelengths should be selected based on areas of the spectrum where the largest change in intensity will result when a spectral shift occurs e.g., a band on the falling side of the emission spectrum .54 Note that selection of wavelengths may therefore determine whether an increase or deSeptember/October 2007 051402-6 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… crease in absorption or fluorescence is observed in response to an increase in membrane potential. CaSDs have much higher contrast than VSDs F/F 5 % and hence, despite generally exciting and emitting in the visible spectrum, are not strongly affected by small hemodynamic changes although the small dip after the calcium response in Fig. 2 may be due to the concurrent increase in HbT and HbO2 absorption . The challenge with CaSDs is introducing them into the neurons successfully, which to date has required the addition of dimethyl sulfoxide DMSO and pluronic acid, followed by direct pressure injection via micropipette into the cortex and incubation of 1 / 2 to 1 h.3,55 However, it is possible to achieve fairly uniform loading of neurons within an area of several millimeters of cortex, and to observe responses similar to those seen with voltage sensitive dyes. Unlike VSDs, CaSDs tend to increase the intensity of their fluorescence emission in the presence of calcium, without significant shifts in the spectral position of their emission and excitation bands. Examples of the use of VSDs and CaSDs for neuroscience studies are given at the end of this section. Imaging of CaSDs with two-photon microscopy is described in Sec. 2.2 Laser speckle-flow imaging. Another recent development in 2-D optical imaging of the exposed cortex is laser speckleflow imaging. While hyperspectral and intrinsic imaging provide measurements of the changes in HbO2 and HbR concentrations, they do not explicitly measure changes in the flow rate of the blood. Flow is a very important parameter if the rate of oxygen delivery and/or consumption is to be calculated. Laser Doppler has been widely used to provide point measurements of blood flow in the brain. However, laser speckle-flow imaging of the exposed cortex can provide a spatially resolved image of the flow throughout the vascular network during functional activation.56,57 Acquiring laser speckle-flow imaging data uses a similar setup to that shown in Fig. 2. However, in this case, the light source is replaced with a divergent laser diode typically at 780 nm . The camera then simply acquires rapid and highresolution images of the exposed cortex, which appear to be low-contrast images of the speckle pattern from the laser illumination. The speckle pattern is caused by the coherent laser light scattering within the brain. If red blood cells are moving within the image, they cause the speckle pattern to vary over time. The rate at which the speckles change relates to how fast the red blood cells are moving. Laser Doppler measures this temporal modulation. However, within the camera image, this intensity fluctuation causes spatial blurring. By calculating the spatial uniformity of an area of camera pixels e.g., 5 5 or 7 7 , the amount of blurring can be quantified from the pixel group’s standard deviation of intensity / mean intensity for each time point . This speckle contrast can be mathematically related to the correlation time, which gives a measure of the velocity of flow of the red blood cells.56 Laser speckle-flow imaging has been performed simultaneously with hyperspectral imaging of the exposed cortex by adding a second camera, a dichroic, and filters to block the laser light . This has enabled the cerebral metabolic rate of oxygen consumption to be calculated.30,36,37 Examples of the use of speckle-flow imaging are given next. Journal of Biomedical Optics 2.1.3 Imaging of exposed cortex: examples Exposed-cortex imaging has been applied to a wide range of research areas. We can broadly summarize these as: 1. functional imaging to investigate the sensory and cognitive processing functions of the brain, 2. functional imaging to improve understanding of the basic mechanisms of the hemodynamic and neuronal response to stimulus, and 3. study of the effects of diseases and treatments on the normal behavior of the brain. Next we describe a selection of recent studies applying exposed-cortex optical imaging techniques within these important brain imaging fields. Functional imaging for sensory and cognitive processing research. Most regions of the brain’s cortex appear to be dedicated to consistent tasks. Particularly in the sensory cortices, there is even spatial correspondence between, for example, the location of two whiskers on the snout of a rat and the two sites or barrels on the somatosensory cortex that respond when each whisker is flicked.58,59 Different parts of the auditory cortex respond to different frequencies of sound.60 Studies of these systems are prevalent throughout the functional imaging literature, but optical imaging of the exposed cortex is currently the only technique that can provide sufficient resolution to reveal the intricate nuances of these hemodynamic and neuronal responses. Figure 3 a shows a recent exposed-cortex optical imaging result revealing the spatially distributed cortical response to visual stimulus in anesthetized rat.61 Cortical imaging data were acquired at 590 nm and the results reveal a very clear spatially correlated retinotopic map of the cortical responses to visual stimuli across the entire visual field of one eye. The purpose of this study was to determine whether it is feasible to use rats as a model to evaluate adaptation of visual cortical function in response to retinal degeneration, and for development of associated potential therapies. Figure 3 b shows an example of a highly complex study of the visual cortex of awake monkeys using both intrinsic signal hemoglobin and VSD imaging.24 Three Macaca fascicularis monkeys had ‘windows’ surgically implanted into their skulls over their visual cortex covering areas V1, V2, and V4 . The cortical response to visual stimulus was then imaged repeatedly for up to 1 year using VSD RH-1691 Optical Imaging Incorporated, Rehovot, Israel to stain the cortex 1 to 3 times per week excited at 630 nm± 5 nm, emission: 650-nm long-pass . For studies of the visual cortex, it is very common to utilize two equivalent visual stimuli and display the difference between the two responses. This is because each response on its own consists of many different changes in the vasculature and tissue, so by subtracting two or more equivalent states, only the very specific difference between the stimuli responses should be revealed. In this case, Fig. 3 b shows intrinsic signal ocular dominance maps 605-nm reflectance generated by showing the monkey a pattern with either his left or right eye blocked the baseline image of the exposed cortex is also shown . The difference between these two states therefore reveals the response of one eye as bright areas and the response of the other as dark stripes. The image shows that this part of the visual cortex processes similar areas of the image from each eye in close proximity to each other, interleaving them, rather than processing the full imSeptember/October 2007 051402-7 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 3 Examples of in-vivo exposed cortex imaging. a Visual stimuli were presented in different parts of a rat’s field of view, and the corresponding cortical hemodynamic HbT responses were mapped color code indicates corresponding location of visual stimulus . Reproduced with permission from Gias et al.61 b A cranial window was implanted to allow direct imaging of primate visual cortex. A baseline image of the exposed brain is shown above an intrinsic signal ocular dominance map, acquired by subtracting the response from visual stimulus to the right eye from the left eye. Chronic VSD imaging was also performed, and ocular dominance maps closely resembled intrinsic signals. Time courses of the voltage-sensitive signal red showed a more prompt onset that the hemodynamic response blue . Reproduced with permission from Slovin et al.24 c An optical fiber bundle was used to image the exposed, VSD stained cortex of an awake mouse while it explored its surroundings. Image series shows the mouse’s whisker gradually brushing against an obstacle and the corresponding cortical response. Reproduced with permission from Férézou et al.65 d A fiber optic bundle was used to image the exposed cortex of a rat undergoing electrical whisker stimulation during simultaneous acquisition of fMRI at 4.7 T. The hemodynamic response shows localized cortical changes in HbR, HbO2, and HbT, which correspond well to oblique-slice fMRI BOLD response study by Hillman, Devor, DeCrespigny and D’Arceuil . e Baseline speckle contrast image showing higher flow as darker contrast and strongly accentuated vessels. Image series from top to bottom show flow prior to induction of spreading depression, 2 min after which a wave of cortical flow increase moves across the field of view, and 20 min after, when hyperperfusion in the dural vessels is prominent. Plot shows time course of flow for region of interest in cortex blue and dural vessel red . Reproduced with permission from Bolay et al.15 f left Cross sectional x-z images through the cortical HbO2 response to a 4-sec forepaw stimulus in rat imaged using laminar optical tomography, top transecting a superficial vein, and bottom transecting the deeper capillary response. right 3-D rendering of arterial, capillary, and venous compartment hemodynamic responses. Ex-vivo two-photon microscopy of vasculature inset agrees with compartment discrimination. Reproduced with permission from Hillman et al.19 Journal of Biomedical Optics 051402-8 September/October 2007 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… ages from the left and right eyes separately. Very similar maps were seen for both VSD and intrinsic signal imaging, although the timing and evolutions of the neuronal and hemodynamic signals were quite different. The data shown here are just the introduction to a very broad study where similar maps were used to investigate the visual cortex’s response to rapid saccadic eye movements. There are numerous other chronic optical imaging studies of monkey visual cortex62 and monkey somatosensory cortex.63 Very recent studies have also demonstrated that plasticity of ocular dominance can even be studied in mice using exposed-cortex optical imaging.64 Another recent novel approach to exposed-cortex optical imaging is shown in Fig. 3 c . While most rodent measurements use anesthetized animals, some studies of native behavior require that the animal be both awake and able to freely move around. Férézou, Bolea, and Petersen developed a technique where a mouse’s voltage sensitive dye-stained exposed somatosensory cortex could be imaged via a flexible fiber bundle while the awake mouse explored its surroundings.65 By simultaneously acquiring video of the mouse’s behavior, they were able to directly image the neuronal response in the brain when the mouse’s whisker contacted an object. Functional imaging to investigate the mechanisms of activation. The functional imaging studies described above exploit the fact that colocalized neuronal and hemodynamic responses can be detected in the brain in response to stimuli.26 However, very little is understood about how and why these functional responses actually occur.40,41,66 Understanding the origins and characteristics of these responses, and the interplay between neuronal and hemodynamic processes, is key to developing an improved picture of healthy and impaired brain function. It is also an essential part of developing and improving all functional brain imaging modalities that rely on the contrast provided by hemodynamic changes such as fMRI. Since exposed cortex optical imaging can provide such a close view of the working brain, as well as access for invasive measurements such as electrophysiology, it has become an established tool for the study of basic functional mechanisms. A wealth of recent studies have investigated so-called neurovascular coupling NVC relationships using optical methods.18,31,38,39,58,67–69 By quantifying the hemodynamic response via optical imaging and/or laser Doppler simultaneously with the neuronal response via electrophysiology , these studies sought to identify a parametric relationship that could be used to predict one response from the other. Devor et al. demonstrated in rat that within a single whisker barrel, the neuronal response to different amplitudes of tactile whisker stimulus saturated before the hemodynamic response.18 More recently, subtleties of the spatiotemporal evolution of both responses have shown that the NVC relationship is highly complex.58,70 A recent study also demonstrated that there can be significant differences between the responses of anesthetized and awake rodents.23 An alternative to attempting to infer neuronal responses from BOLD-type measurements is to instead locate sites where oxygen consumption is increasing by determining the cerebral metabolic rate of oxygen consumption CMRO2 .71,72 In principle, these regions should spatially correlate more closely with neuronal activity than the BOLD response, which can be significant even in distant draining Journal of Biomedical Optics veins.73 CMRO2 is a parameter that can only be calculated by considering the amount of oxygen that is delivered to and extracted by the tissue. Although 2-D hyperspectral imaging can provide changes in HbO2, HbR, and HbT concentration, the data represent discrete snapshots of the current oxygenation state of the blood, and not the rate at which the blood is being replaced. It is therefore necessary to also measure blood flow, to determine the rate of delivery of oxygen to the tissue. Other estimates are also required such as the variations originating in the venous compartment. Both laser Doppler and more recently speckle-flow imaging have allowed simultaneous measurement of oxygenation, blood volume, and flow dynamics in the exposed cortex,30,36,37 thereby allowing direct calculation of CMRO2. These calculations are based on models of the hemodynamic response to stimulus.32,74–78 Interestingly, in-vivo optical imaging and microscopy of the cortex are also allowing the vascular mechanisms of the hemodynamic response underlying these vascular models to be investigated directly.19,79,80 Another approach to comparing exposed-cortex optical imaging findings to fMRI is to perform fMRI studies and optical measurements on the same subjects in response to the same stimulus. Several recent studies have demonstrated simultaneous optical imaging of the exposed cortex during fMRI acquisition.81 An example of such a study is shown in Fig. 3 d , where an MRI-compatible fiber optic endoscope was configured to allow synchronized fMRI and functional optical imaging at two wavelengths 570 and 610 nm during electrical whisker stimulus of varying amplitudes. Similar studies have been performed in humans scheduled to undergo open cortex brain surgery21,82,83: functional responses to voluntary motion of the tongue were recorded in the language area of the brain prior to surgery using fMRI, and then compared to optical imaging results obtained during surgery, where the awake patient was asked to repeat the same stimulus paradigm.83 These results were compared to electrocortical stimulation maps acquired by directly stimulating the cortex and questioning the patient about the resulting involuntary tongue movements. A final optical imaging result of significant importance, and which has prompted much debate, was the discovery of the so-called “initial dip.”33 In some studies, the hemodynamic response to stimulus is seen to begin with an increase in HbR, prior to a subsequent decrease. This finding was interpreted as evidence that neurons initially use up local oxygen more quickly than it is being supplied by the gradually increasing blood flow, an important finding for modeling the mechanisms of the hemodynamic response. It was proposed that this early signal could also allow improved spatial localization to the neuronal response. This initial HbR rise was subsequently also seen in fMRI results as an initial dip in the BOLD signal,84,85 and in many other optical imaging studies.29,37,86 However, the dip is controversial because it is not seen in all studies.87 Reasons proposed for this discrepancy have ranged from differences between species, anesthesia, and paradigms, to measurement techniques and errors in spectroscopic analysis of optical signals especially the wavelength dependence of pathlength88 . An excellent summary of the initial dip can be found in Ref. 89. September/October 2007 051402-9 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Functional imaging to investigate pathologies and treatments. A further use of exposed cortex imaging is to directly investigate the manifestation, characteristics, and treatment responses of diseases and disorders. An example of one such study is shown in Fig. 3 e , which used speckle-flow imaging to investigate migraine in a mouse model.15 Cortical spreading depression CSD is a slow moving wave of neuronal and hemodynamic activity that spreads across the cortex following a small insult to the cortex, such as a pin prick.90 CSD is thought to be the cause of the visual disturbances that often precede a migraine headache. This study demonstrated that speckle-flow imaging could not only see the wave of CSD as a transient increase in parenchymal blood flow, but also revealed a prominent subsequent increase in blood flow in the middle meningeal artery MMA , a vessel that traverses the dura. Such a dural blood flow increase is consistent with the delayed onset of a debilitating migraine headache. Supporting measurements were then made to conclude that CSDs can trigger events leading to headache, identifying new pathways that could potentially be targeted for the development of treatments for migraine suppression. Similar studies using both speckle flow imaging and hypersepctral imaging of the exposed cortex have examined the effects of middle cerebral artery occlusion as a model of stroke.13,91 Both the acute stage and chronic manifestations of stroke can be studied. Function of the area affected by a stroke can be evaluated after days or even weeks of recovery.92,93 The effects of established and experimental treatments can be tested in these models of human stroke, and even compared to the fMRI manifestation of the same outcomes.94 It is also possible to study epilepsy in animal models using exposed-cortex intrinsic signal optical imaging95 and VSDs.96 The hemodynamic response to intrasurgical cortical stimulation during epilepsy surgeries on humans has also been investigated using optical imaging.22 Fast scattering signal. Another source of intrinsic optical signal in the brain is the fast scattering signal. Historically, this signal has been the source of much debate, particularly related to the difficulties in detecting it noninvasively,104,105 although the signal is clearly observable in isolated nerves.106 Recently, exposed cortex measurements during whisker stimulation in rats have provided compelling data revealing localized signals corresponding to regions of neuronal activity.107 Changes in the dark-field reflectance of 660-nm light occurred within 50 ms of stimulation and correspond well with electrophysiological recordings. A 4-mm-diam imaging probe was placed in direct contact with the dura over the exposed somatosensory cortex. The probe consisted of an outer ring of optical fibers delivering illumination light from 660-nm light emitting diodes, and a central fiber optic imaging conduit transmitting remitted light to a sensitive charge-coupled device CCD camera. Several possible mechanisms for the origin of the fast scattering signal have been proposed, including reorientation of membrane proteins and swelling of neuronal dendrites due to water influx in response to ionic currents during depolarization.107 Note that the possibility that fast scattering changes in the brain may distort and interfere with VSD intensity signals has not yet been explored. It is also still not fully understood whether significant scattering changes occur in the brain as part of the slower hemodynamic response. Some researchers have tried to compensate for scattering changes by including a scattering spectrum term in Eq. 2 , as if it represents a third chromophore.33 Care should be taken with this approach, since it does not completely account for scattering effects, and incorporating a third chromophore which may not in fact be changing may adversely affect the fit results of the other two chromophores HbO2 and HbR . Optical coherence tomography of functional brain activation. Optical coherence tomography OCT is a well established technique for high resolution depth-resolved imaging of structures such as the eye.108 OCT overcomes the effects of light scattering to obtain high resolution depthresolved images by using coherence to isolate light that has been directly backscattered. This light therefore has the shortest possible pathlength through the tissue. OCT is also typically performed using NIR wavelengths 1 m. As a result, OCT signals carry very little information about the hemoglobin absorption properties of tissue Eq. 1 and Fig. 1 and cannot be used to measure fluorescence contrast. The contrast in OCT is therefore generally due to refractive index mismatches, and not due to changes in absorption corresponding to hemoglobin concentration or oxygenation. Nevertheless, OCT has recently been applied to brain imaging to explore depth-resolved changes in the cortex of cats and rats during functional stimulus.109–111 To date, the functional OCT changes reported have followed the slower time scale of the hemodynamic response. Both increases and decreases in OCT contrast are seen in close proximity, at times when HbT should be increasing. It is possible that OCT could also detect fast neuronal scattering signals. However, the slower signals observed to date are hypothesized to correspond to changes in the speed and distribution of scattering red blood cells, and changes in diameter and tone of pial vessels and capillaries. September/October 2007 2.1.4 Methodological extras In addition to the three basic 2-D imaging techniques described before, a number of other distinct optical exposed cortex imaging approaches have recently been demonstrated and are described briefly next. Intrinsic fluorophore imaging. Several recent studies have shown compelling evidence that it is possible to measure concentration changes in intrinsic flavoproteins FAD in the brain from their intrinsic fluorescence in vivo.97–100 A system almost identical to that used for VSD imaging is employed, but no contrast agent is added to the brain. Excitation is between 450 and 490 nm and emission is between 500 and 550 nm. FAD provides a measure of metabolic changes in cells neurons and possibly glia , and has a distinct time course compared to shorter neuronal and longer hemodynamic responses. FAD imaging was recently used to explore visual plasticity in mice in a similar way to the studies of the visual cortex described before,101 as well as plasticity in the auditory cortex102 and epileptic activity.103 Since flavoprotein fluorescence requires visible light for its detection, it is unlikely that it will be possible to detect these signals noninvasively see Fig. 1 . Journal of Biomedical Optics 051402-10 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Laminar optical tomography. While 2-D optical imaging of the exposed cortex has found many very valuable applications, there remains concern over the fact that images are a 2-D projection of a combination of signals from superficial vessels and deeper layers of the cortex. The vascular structure of the cortex is organized with large arteries and veins on the surface, and then arteriolar and venous branches diving almost straight down to deeper capillary beds. These capillary beds supply the active neurons, which themselves are also organized in specific layers within the cortex. Therefore, not only are neuronal and vascular signals likely to be quite heterogeneous as a function of depth, but the laminar variations in neuronal and hemodynamic responses are of great interest. Laminar optical tomography LOT is a recently developed imaging technique that can add a depth dimension to hyperspectral imaging of the exposed cortex.19,112,113 LOT utilizes a combination of scanning microscopy instrumentation and optical tomographic image reconstruction techniques to acquire noncontact measurements of the cortex at more than 40 frames per second. It can image to depths of 2 mm with 100 to 200- m resolution. A recent result is shown in Fig. 3 f , where LOT was used to image the exposed-cortex hemodynamic response to rat forepaw stimulus at 473 and 532 nm. The depth-resolved images were validated by comparison to the vascular architecture, and microsvascular density determined via two-photon imaging of fluorescent vascular casts from the same animals. A spatiotemporal delineation technique was also applied that exploited the distinct differences in the temporal evolution of the functional signal in the different vascular compartments: arterioles, capillaries, and veins. Figure 3 f shows two cross sectional x − z images through the HbO2 cortical response, one transecting a superficial vein, and one transecting the deeper capillary response. A 3-D rendering of the three vascular compartments, delineated via their distinctive time courses, is also shown, revealing superficial arterioles and venules, and deeper capillary responses. A two-photon image of a vascular cast from the same animal is shown alongside, with the corresponding arterioles and veins highlighted. The temporal signatures of each different vascular compartment were crossvalidated with in-vivo two-photon imaging and provided new insights into some of the basic assumptions of the hemodynamic models for fMRI interpretation described before.19 Since LOT detects scattered light, it is highly sensitive to blood oxygenation and HbT changes and is also suitable for fluorescence imaging. As a result, LOT in now being developed to allow simultaneous depth-resolved imaging of hemodynamics and VSDs or CaSDs. Depth-resolved imaging of VSDs and CaSDs would allow the laminar properties of neuronal responses to be studied in detail. Currently, such measurements can only be made using highly invasive depthresolved electrode arrays, although demonstration of a gradient index lens-based system also recently showed promising results imaging VSDs in mouse cortex to depths of 150 m.114 1990,115 and has since become an invaluable tool for imaging intact biological specimens at very high resolution, to depths of up to 600 m. Two-photon microscopy has more recently been applied to in-vivo imaging of the brain.3,19,79,116,117 While fluorescent contrast is required, two-photon microscopy provides an unprecedented view of in-vivo brain activity on a cellular and microvascular level. 2.2 Two-Photon Microscopy Two-photon microscopy is a distinctly different functional imaging technique compared to the optical imaging described so far. Two-photon microscopy was first demonstrated in Journal of Biomedical Optics 2.2.1 Two-photon microscopy: basic principles Two-photon microscopy is similar to confocal microscopy, in that it requires a focused beam of laser light to be steered within tissue, sensing the properties of each location and using them to form a 2-D or 3-D image. Confocal microscopy relies on the rejection of scattered light by isolating signals originating from the focus of the scanning beam. Confocal xcitation wavelengths are typically in the visible spectrum, where tissue scatter and absorption are high see Fig. 1 . Therefore, when using confocal microscopy to image to depths beyond 200 to 300 m, the laser beam can no longer focus, and images become blurred and lack sensitivity. Also, while it is assumed that the detected fluorescent light originates from the focus of the scanning beam, some of the detected signal is also generated by excitation and emission light that is scattered or absorbed within the tissue above and below the focus. Not only does this light contribute to image blurring, it can also result in cumulative photodamage to areas of the tissue that are not actively being imaged. Two-photon microscopy overcomes many of the disadvantages of confocal microscopy. Instead of exciting fluorophores at visible excitation wavelengths, two-photon microscopy utilizes laser light of twice the excitation wavelength. When sufficient photon flux is achieved, it is usually possible for a fluorophore to be excited by two of these lower energy photons arriving in quick succession. Once excited, the fluorophore will emit a photon at or near its usual emission wavelength. This is illustrated schematically in Fig. 4 a . This high photon flux is achievable if a pulsed laser is used, which delivers very high energy, rapid pulses typically a Ti:sapphire laser, tunable between 700 and 1000 nm, with 80-MHz pulse repetition, 500-fs pulse width, mean power 1 W, peak power 25 kW . Even with such high peak power, the nonlinear effect of two-photon excitation will only happen at the tight focus of the laser beam. As a result, two-photon microscopy can use NIR light to image the same fluorophores as confocal microscopy. However, since scattering and absorption of NIR light in tissue is much lower than for visible light see Fig. 1 b , the laser beam’s focus can be maintained to depths of 600 m, depending on the tissue type. The fact that fluorophore excitation only occurs at the very focus of the beam means that surrounding tissue does not experience significant photodamage. The NIR excitation light is generally spectrally well separated from the emission wavelength, making rejection of excitation light much simpler than for confocal microscopy. Further, since it is expected that any light emerging from the tissue at the fluorescent emission wavelength at a particular time can only have originated from the focus of the scanning beam, it is not necessary to attempt to reject scattered light, as required in confocal microscopy. Instead, sensitive detectors can be placed as close to the tissue as possible to collect all of the emerging emission light. This September/October 2007 051402-11 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 4 Two-photon microscopy of in-vivo brain function. a Basic mechanism of two-photon fluorescence. b Schematic of surgical preparation of exposed cortex, with sealed glass window and microscope objective positioning. Green dot shows location of two-photon fluorescence. c Examples of two-photon maps of the vasculature following intravenous injection of dextran-conjugated fluorescein. Black dots and stripes show red blood cell motion. d Dual-channel imaging of neuronal green and vascular red signals: left Oregon Green 488 BAPTA-1 AM calciumsensitive dye stained neurons and right transgenic mouse expressing green fluorescent protein GFP in a subpopulation of neurons mouse supplied by Friedman, Rockefeller University, New York130 . Texas dextran red is the intravascular tracer in both cases studies by Hillman, Bouchard, Ruvinskaya and Boas . e Three-channel imaging of Tg2576 APP Alzheimer’s disease mouse model with amyloid- targeting dye blue , GFP expressing neurons and dendrites green , and vasculature red . Reproduced with permission from Spires et al.14 f 2-D camera-based imaging of hemodynamics and calcium response followed by functional two-photon imaging of the calcium response on the scale of single neurons in the same rat. Oregon Green 488 BAPTA-1 AM was loaded with serial pressure injections into the cortex. Stimulus was 1 mA, 3- s pulses at 3 Hz for 4 sec studies by Hillman, Bouchard, Ruvinskaya and Boas . Journal of Biomedical Optics 051402-12 September/October 2007 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Table 1 Summary of optimal two-photon microscope characteristics for in-vivo brain imaging. Custom-built in-vivo two-photon microscopy systems Almost always upright. Commercial two-photon microscopy systems Often inverted. Feature Form factor Imaging platform Require Z translation of stage, poor clearance between stage Fixed Z stage objective moves up and down, not stage . Motorized x − y stage. Very open clearance around objective and objective. Motorized x − y stage movement very expensive. for animal and monitoring equipment. Two-photon detector may be in place of confocal detector, therefore descanned and through pinhole—very poor light efficiency. Otherwise retrofitted detectors are rarely close to the objective. Similar objectives can be chosen for commercial systems if compatible. Detection geometry Detectors positioned very close to back of objective for maximum light collection. Multiple spectrally resolved detectors may be stacked to allow two or more emission wavelengths to be simultaneously imaged. Numerical aperture Very high numerical aperture but long working distance allows good z resolution, high efficiency, and ability to image through exposed-cortex preparation see Fig. 4 b . Scan rate Software Multimodality Typically designed to allow very high frame or line-scan rates Frame rate well optimized, but limited adaptability for to image dynamics of functional activity. patterns and line scans. Custom software allows wide variety of scan patterns and routines for a wide range of experiments. Typically possible to add additional imaging and measurement modalities to the system owing to accessible light path and animal. Can be integrated into custom software. Often proprietary and difficult to modify. Fixed scan patterns and routines. Very difficult to modify or interrupt light path. Journal of Biomedical Optics impossible. As a result, many researchers have instead built their own two-photon microscopy systems, optimized for invivo imaging.122–126 Table 1 summarizes the key attributes of in-vivo two-photon microscopy systems that are specifically designed for in-vivo imaging. 2.2.2 Two-photon microscopy: examples As for 2-D optical imaging, in-vivo brain imaging research using two-photon microscopy has also found applications in areas of functional processing, functional mechanisms, and disease research. These applications utilize a variety of methods to introduce fluorescent contrast to the brain, including intravenous injection of dextran-conjugated dyes to reveal blood vessels, topical application, or pressure injection of dyes into the cortex, systemic delivery of dyes that cross the blood brain barrier, and transgenic mutation of cells to express fluorescent proteins. Examples are given below. Functional processing. A recent example of in-vivo twophoton microscopy of brain function was the use of calciumsensitive dyes to look at visual responses in cats and rats.3,127 Data were acquired using a custom-built two-photon microscope, and used Oregon Green 488 BAPTA-1 AM dye Invitrogen excited at 800 nm to provide calcium-sensitive fluorescent contrast. The dye was pressure injected into the cortex as described before for 2-D imaging of CaSDs. The complex cortical patterns of the response to visual stimuli of different orientations similar to the ocular dominance maps in Fig. 3 b were found to occur even on the scale of discrete individual neurons. September/October 2007 051402-13 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Functional mechanisms. The advent of functional twophoton microscopy has allowed the mechanisms of neurovascular coupling to be studied in vivo in unprecedented detail. Whereas 2-D optical imaging of the cortex, as described earlier, has enabled parametric comparisons between electrophysiology and the hemodynamic response,18 two-photon imaging enables the interrelations between individual neurons, glial cells, and vessels to be examined in intricate detail. Both confocal microscopy128 and two-photon microscopy79 have been used to examine the dynamics of red blood cells flowing through capillaries in the exposed cortex. In Ref. 79 and more recently in Refs. 116 and 129, the capillary response to stimulation and hypercapnia have been examined using fluorescent plasma tracers such as dextran conjugated fluorescein or Texas dextran red Invitrogen . These high molecular weight dyes cannot pass the blood brain barrier or stain red blood cells, and so provide strong contrast to vessels. Typical in-vivo image stacks of fluorescein-perfused vessels in rat brain are shown in Fig. 4 c . The dark stripes and dots are moving red blood cells whose dynamics can be quantified from their motion in sequential images or line scans.126 The flow and dilation dynamics of pial arterioles and veins was also recently investigated using two-photon microscopy.19 Figure 4 d left shows simultaneously acquired in-vivo rat brain measurements of neurons labeled with Oregon Green 488 BAPTA-1 AM dye and vessels containing Texas dextran red. Both color channels were acquired with two-photon excitation at 800 nm and 550-long-pass dichroic and emission filters to separate the two emission bands onto different photomultiplier detectors. Figure 4 d right shows similar data acquired in a transgenic mouse, engineered such that only certain neurons in the brain express green fluorescent protein.130,131 Capillaries perfused with Texas dextran red can be seen to be weaving amidst the webs of dendrites. 2-D optical imaging of hemodynamic and CaSD changes can also be acquired prior to imaging the same animal with two-photon microscopy. Figure 4 f illustrates 2-D imaging of the calcium-sensitive Oregon Green 488 BAPTA-1 AM response in the somatosensory cortex of a rat during forepaw stimulus 1.5 mA, 3- s pulses at 3 Hz for 4 sec , followed by two-photon imaging of individual neurons at a depth of 315 nm in the cortex responding to the same forepaw stimulus. The two-photon fluorescence of the CaSDs increased by up to 17% during stimulation.126 Such direct and simultaneous observations of both neuronal and hemodynamic responses and structures in in-vivo systems enable the intricate mechanisms underlying neurovascular coupling relations to be closely examined and compared to macroscopic-scale measurements and even fMRI results. Two recent studies used similar two-photon imaging techniques to present compelling evidence supporting the role of astrocytes as active mediators of blood flow control.80,132 Astrocytes make up a large fraction of the cells in the brain, and have been shown to be both metabolically active and capable of transmitting signals between neurons and blood vessels.133,134 In Ref. 132, astrocytes in mouse brain were labeled with calcium-sensitive rhod-2 AM, caged Ca2+, and DMNP-EDTA AM, while the vasculature was perfused with dextran conjugated fluorescein. Two-photon imaging of astrocytic endfeet contacting arteries was performed in vivo at 825 nm, with spectrally separated detectors for rhod-2 AM Journal of Biomedical Optics and fluorescein. A second laser 355 nm was then focused onto an endfoot and pulsed to cause photolysis of the DMNPEDTA, triggering intracellular uncaging of calcium in the endfeet. This forced increase in calcium caused the endfoot to dilate. In a subsequent study, a CaSD fluo-4 AM and astrocyte-specific dye SR101 were used to demonstrate that similar calcium increases in astrocytic endfeet could be observed during prolonged whisker stimulation in vivo.80 Functional imaging to investigate pathologies and treatments. In-vivo two-photon imaging of intravascular dextranconjugated fluorescein in rats was recently used to investigate the effects of stroke on blood flow dynamics, as well as the recovery of flow following two treatments: hemodilution and injection of a recombinant tissue-type plasminogen activator.135 The study also created the strokes using a second laser delivering targeted high fluence ultra-short NIR pulses to cause vessel rupture, intravascular clotting, or extravasation of blood components depending on the laser power. Alzheimer’s disease research has been greatly impacted by the availability of fluorescent markers that bind to beta amyloid plaques in the brain.5,136,137 The so-called Pittsburg compound B 2- 4’-methylaminophenyl -6-hydroxybenzothiazole, PIB can be used as a clinical contrast agent for PET imaging and is showing promise as a diagnostic tool in humans.138 However, since PIB is fluorescent, it can also be used as a two-photon microscopy contrast agent for exposed cortex imaging.117 Tg2576 APP transgenic mice develop beta amyloid plaques over the course of their lives, and are assumed to provide a realistic model of the progression of Alzheimer’s disease. Since Alzheimer’s in humans can currently only be confirmed postmortem, these mouse models offer an opportunity to investigate the early effects of the disease, as well as learn about its progression and to allow development of earlier treatment and diagnostic methodologies. Mice with chronically implanted cranial windows can be repeatedly imaged with two-photon microscopy for up to a year. This has allowed two-photon microscopy to be used to evaluate the effect of developing amyloid- plaques on the structure and function of surrounding neurons and blood vessels. Figure 4 e shows three-channel two-photon image stacks of the brain of a 21 to 24 month old transgenic mouse with beta amyloid deposits shown in blue, blood vessels perfused with Texas dextran red shown in red, and neurons expressing green fluorescent protein in green. The amyloid binding dye methoxy-XO4, similar to PIB was injected intraperitoneally the day before imaging. The neurons were labeled by transfection from an adeno-associated virus injected into the brain 2 to 3 weeks prior to imaging.14 Extensions of this work to noninvasive animal and perhaps even human imaging are being explored.139 2.2.3 Methodological extras Fluorescence lifetime imaging and microscopy FLIM has also been applied to neuroimaging, allowing potential enhancement of intrinsic contrast via environment-dependent lifetime shifts. FLIM has also been utilized to improve measurements of fluorescence resonance energy transfer FRET , September/October 2007 051402-14 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… which can provide measures of the proximity of two fluorescently labeled entities and therefore direct visualization of biochemical interactions.140 There is an extensive body of neuroscience literature that utilizes living brain slices to study neuronal activity in intricate detail. Mice and rats are rapidly sacrificed before their brains are removed, sliced, and placed in a carefully prepared perfusate that supplies nutrients and oxygen. Important results have been obtained by imaging brain slices with two-photon microscopy. For example, using brain slices, Kasischke et al. demonstrated that it was possible to use two-photon microscopy to measure the intrinsic fluorescence of NADH, and distinguish between cytoplasmic and mitochondrial NADH dynamics.141 It has also recently been shown that voltagesensitive responses in neurons and dendrites in brain slices can be measured via second harmonic generation SHG microscopy when neurons are loaded with a SHG contrast enhancing agent called FM4-64.142,143 Noninvasive optical brain imaging on animals is generally performed as a precursor to clinical imaging in babies and adults. It may also be used to improve understanding of clinical data, including studies of neurovascular coupling, stroke, birth asphyxia, etc.166–169 Noninvasive animal optical brain imaging has also been used for longitudinal development studies,170 and measurements where exposing the cortex is not feasible, or would detrimentally affect the experiment. There are two main technological approaches to noninvasive optical brain imaging: fairly simple topographic measurements of cortical function as shown in Figs. 5 a and 5 b , and 3-D tomographic imaging of the brains of infants as shown in Figs. 5 c and 5 d . The basic principles of these approaches are described next, followed by specific examples of their in-vivo neuroimaging applications. 3.1 Topography 3.1.1 Experimental method Noninvasive topographic optical measurements of the human head can interrogate hemodynamic changes in the cortical surface of the brain Fig. 1 b . The majority of functional processing occurs in the cortex, and so optical topography is commonly, and extensively, used to measure functional responses to stimulus in children and adults. Early classifications of the roles of the different areas of the cortex were deduced from observations of the effects of injuries to different parts of the brain.171,172 More recently, these cortical areas have been refined and investigated using functional brain imaging modalities, including fMRI and animal imaging studies such as those described in Sec. 2. Figure 5 e shows the basic regions of the cortex and their functional roles. Optical imaging studies of brain function often focus on the occipital or visual cortex at the back of the head looking at responses to visual stimulus ,147–150 and the motor and somatosensory cortices between the crown and the ear in response to movements, tactile stimulation, or pain .151–154 However, many other studies have been reported, including those on the parietal,155 frontal,156 temporal,157 and olfactory158 lobes. In most functional imaging experiments, a well-defined stimulus is devised such as tapping fingers, or watching a pattern or picture . A subject then performs or receives the stimulus at regular or carefully controlled intervals. In almost all cases, data are then processed to isolate signals that correlate with the presentation of the stimulus. In all functional imaging modalities, and especially for higher mammals such as humans, there is always some variability in baseline brain activity during an experiment that can generally only be eliminated by performing repetitive trials and averaging or filtering. fMRI is rapidly becoming a clinical tool with which to investigate the effects of diseases and drugs, to plan surgeries, and to study brain development.9 Optical topography is also finding similar clinical applications.147,163–165 3.1.2 Basic principles: near-infrared spectroscopy for topography Direct topography. Direct topography is the most widespread approach to clinical optical brain imaging and is often referred to as near-infrared spectroscopy NIRS . In its simSeptember/October 2007 Noninvasive Optical Brain Imaging For routine optical brain imaging in humans, it is clearly necessary to develop noninvasive techniques. Approaches must overcome the effects of light scattering, while maintaining the strong benefits of functional optical contrast. To date, noninvasive optical brain imaging has been applied to babies, children, and adults for studies of functional activation, and for investigation of pathologies. Since light penetration and scattering are significant obstacles for noninvasive imaging, measurements on children and babies are easier to achieve. It is often difficult to image these infants and small children using fMRI and PET, and so optical imaging provides a unique opportunity to study many aspects of functional brain development that have previously been inaccessible.144 Imaging of pathologies in infants is also highly motivated by the prevalence of perinatal brain injuries due to prematurity, birth asphyxia, and congenital conditions.34 Thanks to the portable and nonionizing nature of optical imaging, it is possible to obtain optical measurements on babies within a neonatal intensive care unit,145 and on young children who are able to sit on a parent’s lap in a familiar environment.146 Optical imaging of the adult brain can also provide valuable functional information that complements modalities such as fMRI. A particular advantage is optical imaging’s ability to simultaneously image HbO2, HbR, and HbT. Optical imaging is being widely adopted for studies of the cortical hemodynamic response to a wide range of stimuli.147–158 Multimodality imaging of the adult head during fMRI,150,159,160 electroencephalography EEG 161 and magentoencephalography MEG 162 has provided insights into the contrast mechanisms of the fMRI BOLD signal, and allowed oxygen consumption and neurovascular coupling relations to be investigated in humans. Optical imaging in adults is also being used as a tool for monitoring and detecting hematomas and stroke, as well as the effects of other pathologies and surgery on brain perfusion.147,163–165 Noninvasive optical imaging holds the potential to provide accessible monitoring of brain function in situations where other imaging modalities such as fMRI are prohibitive, e.g., due to imaging environment or cost. Journal of Biomedical Optics 051402-15 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 6 Time-domain, frequency-domain, and continuous wave measurement sensitivities to absorption and scatter. Schematic shows simplistically how absorbing and scattering structures will change TD and FD data, but that in transmission, their effect on cw data will be indistinguishable. Fig. 5 Configurations for noninvasive optical brain imaging. a and b Functional topography uses one or more source-detector pairs to create a pseudo-2-D map of the underlying cortex. Sources and detectors may be optical fibers, or light emitting diodes and photodetectors placed directly on the skin. c and d Tomography places sources and detectors over the volume of the head, and detects light that has traveled through many projections. This approach allows subsequent reconstruction of the optical properties of the entire brain. e The main functional areas of the human cortex. f Simulation of the sensitivity of a noninvasive topographic measurement between a source and detector. The simulation used Monte-Carlo modeling and accounted for the heterogeneous 3-D structures of the adult scalp and cortex. Reproduced with permission from Boas et al.202 plest form, NIRS uses light source and detector pairs positioned on the head of infants or adults. Light that passes through the skull scatters within the cortex of the brain as illustrated in Figs. 1 b , 5 b , and 5 f . Despite significant attenuation, variations in absorption due to changes in local HbT, HbO2, and HbR concentration result in detectable modulations in the intensity of the emerging light. If an array of sources and detectors is used, the changes measured between each source-detector pair can be approximately mapped or interpolated to the underlying cortex, and a low-resolution 5 mm 2-D topographical map of cortical activation can be produced.152,154,173 When many overlapping and multidistance topographic measurements are acquired, it is possible to use mathematical models of light propagation to create a depthresolved reconstruction of functional activation.174–176 While more complex, this approach can overcome several of the limitations of direct topography, as described below. Instrumentation. Noninvasive optical brain imaging instrumentation may be continuous wave cw , frequency domain FD , or time domain TD . cw is the most common approach and simply measures the intensity of light emerging from the head at two or more wavelengths. Data acquisition is often frequency encoded up to kHz , such that multiple wavelengths177,178 or multiple sources149,179–181 can be measured simultaneously. FD and TD instruments are substantially more complex and expensive, and can measure both a cw component and the time taken for light to pass through the tissue.182 FD instruments achieve this by modulating their illumination at around 100 MHz and measuring the phase delay of the emerging light.152,183 TD instruments use a pulsed laser and measure the temporal distribution of the photons emerging from the tissue.145,184–188 The main purpose of this temporal information is to provide a measure of the scattering properties of the head, as illustrated in Fig. 6. Optical fibers are generally used to transmit light to and from the head see Fig. 5 . Direct placement of small sources and detectors on the head has also been demonstrated.166,189 Calculation of HbO2 and HbR concentrations. The change in intensity measured between a topographic source and detector the attenuation will generally reflect a change in the September/October 2007 Journal of Biomedical Optics 051402-16 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… absorption of underlying tissue Eq. 1 . The basic approach for deriving changes in HbO2 and HbR concentration from changes in attenuation was outlined in Sec. 1.1.2 Eqs. 1 and 5 . Determining I0 and G for noninvasive measurements is particularly challenging, since the amount of light coupled into and out of the scalp for every source and detector position would need to be known. Differential measurements relative to time t = t0 Eq. 3 should cancel these source- and detector-specific scaling factors, although it must be assumed that these factors stay constant for the duration of the measurement which can be challenging if the subject moves and/or hair is obstructing probe placement . As with exposed cortex imaging, if measurements of the change in attenuation at multiple NIR wavelengths are to be converted into changes in hemoglobin concentrations, the effects of scattering must be accounted for. If Eqs. 4 and 5 are to be used for “direct” topography, an estimate of DPFx for each wavelength and measurement pair is required. FD and TD systems can be used to measure DPFx directly, and as an adjunct to cw imaging can provide DPFx for a given wavelength by measuring the average time taken for light to pass from the source to detector, and converting to distance based on the speed of light in tissue.190,191 However, if a subset of DPFx measurements or generalized values of DPF are to be used, it should be noted that DPFx could vary depending on age, gender, and even position on the head, owing to the complex structures of the underlying scalp, skull, and convoluted human cortex. Simulations for calculation of DPFx can become increasingly complex if efforts are made to incorporate the curvature of the head, the heterogeneous structures of the scalp, skull, and cortex, the presence of cortical convolutions and sulci, the effects of CSF see Fig. 1 b , and even the effects of isotropic scattering in the brain.88,192–199 Estimates of the absorbing and scattering optical properties of each tissue type,200,201 as well as accurate knowledge of the source-detector separations between fibers on the head, are also required. If available, an MRI scan of the head can be used to build a customized model of light propagation for each subject, allowing constrained analysis that accommodates the exact structure of the brain, skull, and scalp underlying each detector.202 However, while this is feasible in studies where fMRI is simultaneously acquired, it is prohibitive for other applications of optical imaging. Although TD and FD systems rarely have a large number of multiwavelength channels,161,185 they have the advantage that they can be used to deduce the scattering properties of the head, which can then be incorporated into more complex calculations of hemoglobin concentrations.201,203 An FD approach is to measure between several source-detector separations and, under simplifying assumptions of homogeneity, fit for scattering and absorption based on a linearized model of light propagation versus distance.183 Such approaches do not require estimates of DPFx. In general, for direct topography, simplified models and scaling factors are applied to account for DPFx and also partial volume effects, which correct for the fact that only a fraction of the detected light actually sampled the cortex.34,153,159,183,204 Several studies have explored the importance of accurate calculation/measurement of DPFx and its Journal of Biomedical Optics impact on the determination of hemoglobin concentrations for optical topography.198,205,206 It has been shown that in many cases, useful functional information can be obtained even when simplifying assumptions are applied to analysis and interpretation. However, it is therefore important to consider the possible impact of systematic errors resulting from this analysis for all NIRS measurements, particularly when imaging the adult head, where many confounding effects can obscure cortical signals and impact correct interpretation. Extracerebral effects and physiological noise. The effects of blood flow changes in the scalp, and changes in the brain that are not directly related to the cortical stimulus response, must also be considered. The analysis approaches described above assumed that changes in signal were solely due to changes in the cortex of the brain. However, it has also been shown that breathing 0.16 Hz , heart rate 1 Hz , and 0.1 Hz 207–210 can the “vasomotor” or Mayer wave signal cause large changes in measured hemodynamic signals. Substantial repetition of a stimulus and averaging should reduce the effects of unrelated physiological variations in blood flow in the scalp and brain. However, acquisition of so many repetitions is not always possible with adult or infant subjects. In addition, some responses may not average out, or may intrinsically correlate or become entrained with the stimulus period,208 and can affect both brain blood flow and the rich blood supply to the scalp. Recent evidence has suggested that certain stimuli actually invoke a systemic response possibly through heart rate or blood pressure changes that could manifest as a change in scalp or even baseline cortical blood flow that correlates with the stimulus period.153,211 Several measurement, data analysis, and filtering approaches have been proposed to reduce the influence of hemodynamic activity unrelated to the true cortical response. One approach is to simultaneously acquire measurements between close and more distant pairs of sources and detectors, assuming that the closer detector is predominantly detecting signal from the scalp. This scalp signal can then be approximately removed from the signal from the more distant detector, assuming that a similar superficial contribution is affecting both measurements. A similar approach can be achieved with TD data, where the photons arriving at a detector earlier are likely to have traveled more superficially, and therefore contain information about the scalp signal. The benefit of using TD is that later photons that have traveled more deeply between the same source detector travel through the same region of the scalp as the early photons, so cancellation should be improved relative to the multidistance approach.187,212 An extension of this is reconstructed topography, which creates a depth-resolved map that should be able to spatially separate superficial signals from deeper cortex signals, if suitable measurements are acquired. Spatiotemporal filtering and modeling have also been proposed as a way to remove these systemic signals.213,214 While rapidly becoming more sophisticated, implementation of these data processing techniques can be complex and sometimes subjective. Extensive processing could result in distortion of true cortical signals, particularly when attempting to separate scalp and cortical signals that may both have elements correlating with the stimulus response.211 September/October 2007 051402-17 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Practical measurements. A final difficulty with acquiring functional optical imaging data is reliable placement of sources and detectors on the head. The primary problem is to be able to securely position multiple probes in a pattern such that their interfiber distances are known, and they stay stable during the course of the measurement. Hair on the head is a major challenge, as it can significantly attenuate light entering and leaving the head, and can cause measurement instability over time if the probes shift. Technical advancements have been made in positioning optical fibers and probes onto the head, and relatively stable configurations have been devised.180 However, further work will be required before it becomes possible to routinely, rapidly, and repeatably position optical topography fibers on the head. Reconstructed topography. Direct topography is the simplest way to analyze topographic data, but relies on many assumptions and simplifications. An alternative approach is to reconstruct data onto a 3-D rendering of the head. Although analysis becomes more complex, by creating a depth-resolved image of the measured absorption changes, some of the confounds described earlier can be overcome. Image reconstruction. Instead of requiring a pathlength estimate, reconstruction effectively attempts to recreate the measured data by simulating alterations in the absorption distribution within the head.187,215,216 When the simulated data and the measured data match, the absorption distribution is assumed to be correct. Reconstructing still generally requires an estimate of the baseline absorbing and scattering optical properties of the head, and an accurate model of light propagation, but can help if there are depth-dependent variations in the response, including effects from the scalp. Partial volume effects, particularly where one wavelength of light effectively probes a different region, than a second wavelength, can also be accommodated in reconstructed topography. For reconstructed topography to be effective, it is necessary to measure signals between overlapping sources and detectors, and generally at several source-detector distances see Fig. 5 b . These measurements provide a set of mutually consistent data for which there should be only one possible corresponding differential absorption distribution. This in turn should also lead to improved x-y resolution.175,176,179 Modeling of light propagation in the human head typically uses either Monte Carlo modeling or the diffusion approximation DA to the radiative transfer equation Eq. 6 and a linearizing assumption such as the Born Eq. 7 or Rytov approximations. − · = q0 , M n,m = Jn,m r r , where is the photon density within the tissue, and q0 is the source function. M n,m is the set of changes in measurement from time t = 0 between source n and detector m, and Jn,m r is the Jacobian matrix derived from the Greens function solution to Eq. 6 , representing the change in measurement M n,m corresponding to a change in absorption a at position r within the tissue.215 An example of Jn,m r for a single Journal of Biomedical Optics September/October 2007 051402-18 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… rates.152,181 However, the attenuation difference between a source and detector separated by 2.5 and 4.25 cm can be up to two orders of magnitude.179 For reconstructed topography, these multidistance measurements are vital, and detectors must have excellent dynamic range and noise characteristics to be able to simultaneously distinguish between frequencyencoded signals of such varying amplitudes. Development of faster systems, with better dynamic range or highly advanced and rapid multiplexing schemes is therefore ongoing.176,179 3.2 Tomography The most complex approach to optical brain imaging involves performing a full 3-D reconstruction of the entire head using an approach similar to x-ray computed tomography. This approach has only been demonstrated on small infants, since it requires light to be transmitted across the whole head as shown in Figs. 5 d .145,221–223 However, the benefits are that in addition to functional imaging, this approach can potentially enable imaging of pathologies such as periventricular hemorrhage and ischemia in a perinatal setting, where currently no alternative diagnostic and monitoring tools are available.217 3.2.1 Basic principles Since the use of 3-D tomography for brain imaging is not yet widespread, specific details of one study are given in the next section, and only a brief overview is provided here. For tomographic imaging, it is necessary to acquire a series of measurements or projections through the tissue being imaged. For reconstructed topography, measurements between sources and detectors with different separations allow some depth-resolved information to be deduced. In tomography, the same basic image reconstruction approach is used, but a larger number of projections through the head at different angles are required to calculate the 3-D structure deep within the brain see Figs. 5 c and 5 d . Note that optical tomography has also been used to image the adult breast for the detection or monitoring of breast cancer224 and the human forearm.225 There are three main difficulties with acquiring tomographic data on the head. The first is that light must be detected through greater than 7 cm of tissue, which in general will cause around 10−7 attenuation of the incident light. The second difficulty is that the paths of light through the head are again a function of baseline anatomical structures. Third, to look at pathologies or regions of ischemia for which a baseline e.g., prestimulus state is not available, it can be valuable to attempt to reconstruct this absolute absorbing and scattering structure of the brain, and not just measure changes as with topography . 3-D tomographic imaging of the neonatal had has been demonstrated using TD instruments.145,222,223 TD measurements that use single photon counting have the advantage that very low light levels can be detected by increasing measurement integration time. Background signal and thermal noise will not be temporally correlated, and can often be subtracted from the true signal.226 Therefore, while slow, TD systems can detect usable signals across 15 cm of adult breast, or an entire neonatal infant’s head. TD and FD measurements also have the benefit that they allow the effects of scatterers and absorbJournal of Biomedical Optics ers to be distinguished as illustrated in Fig. 6,215,227 although see Ref. 228 . In fact, using TD data, it is possible to reconstruct images of both absolute baseline absorbing and scattering structures.229 These methods are currently being developed to allow imaging of the baseline optical properties of the infant brain. Differential measurements can be used to reconstruct absorption changes, under the same assumptions as topography: that scattering does not change and the effects of heterogeneous background absorbing and scattering structure on light paths is insignificant or can be incorporated into a model of the head.217,219 Differential imaging of 3-D absorption changes can be achieved with relatively simple linearized reconstruction approaches similar to those used for reconstructed topography described before.217,222 Values of Jn.m r must be calculated based on the 3-D shape of the head, and the precise positions of the sources and detectors see Fig. 5 c . Typically, this is done using finite element models implementing the diffusion approximation, and sometimes incorporating estimates of the baseline absorbing and scattering structures within the head.230,231 For absolute imaging of the baseline absorbing and scattering structure of the brain, a nonlinear algorithm should be used that iteratively seeks a 3-D distribution of absorbing and scattering that can account for the data measured on the head. A highly complex and publicly available software package developed to do both nonlinear reconstruction of absorption and scattering images, as well as forward modeling of Jn.m r using the finite element method is called TOAST timeresolved optical absorption and scattering tomography .215,232,233 The data that are actually fed into these reconstruction algorithms will have a significant effect on image sensitivity to scatter and absorption.217 FD data consists of measurements of amplitude and phase. Measurements are typically made at modulation frequencies between 100 and 200 MHz.234 TD data represents the full impulse response of the tissue, and therefore amplitude and phase at almost any effective modulation frequency can be calculated. Datatypes such as the mean-time or integrated intensity of the signal can also be extracted.226 For images of both absorption and scatter, the image reconstruction algorithm must utilize both amplitude and phase data, or at least two sufficiently different timeresolved parameters. A very accurate model of light propagation and precise calibration of data for systematic errors is also very important, particularly if absolute images of absorption and scatter are required since the data must closely match simulated data corresponding to the correct 3-D distribution of absorption and scatter .217,226 3.3 Examples of Noninvasive Optical Brain Imaging Figure 7 shows a variety of noninvasive optical brain imaging results showing human functional responses to stimulus. Figure 7 a shows results from a study using cw direct topography to examine the hemodynamic response to different kinds of stimuli: finger tapping, median nerve electrical stimulation, and tactile stimulation.153 Stimuli were presented for 10 sec, followed by 18 sec of rest, and were repeated 20 times per subject. The cw system has 16 laser diode sources, eight at 690 nm and eight at 830 nm, frequency encoded between 4 September/October 2007 051402-19 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 7 Examples of noninvasive optical imaging of functional responses in humans. a Direct functional topography images and time courses acquired for two different motor/somatosensory stimulus paradigms on one subject. A localized HbO2 and HbR response is seen in the opposite contralateral side of the brain to the hand being stimulated. Reproduced with permission from Franceschini et al.153 b Functional topography measurements on adult head during fMRI. Both the spatial and temporal characteristics of the responses can be compared between two modalities. Reproduced with permission from Huppert et al.159 c Black and white visual stimuli were presented to the left and right fields of vision. The hemodynamic response was measured using a dense topographic array of sources and detectors. Data were reconstructed using a linear perturbation model. Reproduced with permission from Culver et al.176 d Time-resolved 3-D tomography data were acquired on a 33-week gestation infant undergoing passive motor stimulus to the left arm. Image series shows successive slices through 3-D head from left ear to right ear. Reproduced with permission from Gibson et al.222 and 7.4 kHz such that all were illuminated while 16 avalanche photodiode detectors collected emerging light. Sources and detectors were configured in two panels, arranged as shown in the photograph in Fig. 7 b where four fibers span the center each delivering both 690 and 830 nm—red arrow , straddled by two rows of four detector fibers—green arrows. The detected signal was digitized at 40 kHz and analyzed to extract signals coming only from the closest sources 3 cm to each detector. Data were then low-pass filtered and downsampled to 10 Hz, and then bandpass filtered between 0.02 and 0.5 Hz to eliminate heart beat pulsations and slow drifts. Data were then block averaged over the 20 trials to reveal the average attenuation changes corresponding to each stimulus. These attenuation values were converted to absorption changes using literature values of the differential pathlength factor DPF at 690 nm 6.51 and 830 nm 5.86 , such that the pathlength is DPFx, where x is the source-detector separation of 3 cm. The resulting absorption changes were then converted to HbO2 and HbR changes using Eq. 4 . These source-detector-specific changes in hemodynamics were then projected onto a map underneath each source-detector pair using an approximate weighted spatial distribution based on the volume sampled. The final map was then interpolated onto a finer grid. The images and time courses in Fig. 7 a show the mean maps and time courses for a single subject. Contralateral responses are those seen in the side of the brain Journal of Biomedical Optics opposite the hand receiving the stimulus where responses would be expected . Figure 7 b shows a similar study, but where optical topography data were acquired during simultaneous fMRI at 3 tesla.159 The subjects performed a finger tapping exercise for 2 sec, with a randomized interstimulus interval between 4 and 20 sec repeated between 27 and 32 times per fMRI scan, with each scan repeated four to six times. fMRI scans with a 2-Hz frame rate were used to measure either BOLD contrast or ASL data, which provides blood flow information. The cw NIRS instrument and data processing approach was similar for this and the study described previously.153 The HbO2, HbR, and HbT time courses for one subject and for each source-detector pair are shown in map form inset lower right of Fig. 7 b . Oil-filled fiducials on the optical fiber heads allowed the optical probe’s location to be coregistered with the MRI scans. The time course of the fMRI BOLD signal is compared to the HbO2, HbR, and HbT temporal responses measured using NIRS. The group average of five subjects is shown, and the work concludes that BOLD most closely temporally correlates with HbR changes. ASL changes were found to correlate most closely with HbT changes data not shown . A second paper continues this analysis to a spatial comparison between the localization of HbR and BOLD responses.235 September/October 2007 051402-20 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Figure 7 c shows a reconstructed topography result, acquired on the adult head during a 10 sec visual stimulus followed by 35 sec of darkness.176 Data were acquired with a multichannel cw NIRS system consisting of 24 radio frequency encoded light-emitting diode LED sources at 760 and 850 nm, and 28 avalanche photodiode detectors. Frequency encoding allowed simultaneous acquisition at both wavelengths, but LEDs were switched on sequentially allowing detectors with 107 dynamic range to measure signals from sources at three different separations: 13.3, 30.0, and 40.0 mm. Despite this serial source switching, full images were acquired at 10 frames per sec. Data were then reconstructed by inverting the linear problem Eq. 7 and using Jn.m r functions calculated using a finite element model simulation package called NIRFAST,236 which incorporates a simplified baseline structure consisting of a 15-mm scalp/ −1 −1 surrounding the skull layer a = 0.17 cm , s = 7.6 cm −1, −1 . brain of radius= 65 mm a = 0.19 cm s = 11.5 cm Both block averaged and even trial-by-trial data show clear hemodynamic responses to unilateral visual stimulus. Figure 7 d shows 3-D optical tomography images of the functional response to passive motor stimulus of the left arm in the brain of a premature infant.222 The image series represents vertical slices through the 3-D volume of the head from left ear to right ear . A clear volume of absorption increase is seen in the region corresponding to the right motor cortex. Data were acquired using a dual wavelength, 32-channel time-resolved optical tomography system known as MONSTIR.186 The images in Fig. 7 d show the change in absorption at 780 nm, but simultaneously acquired images at 815 nm were also created and converted into maps of the change in HbO2, HbR, and HbT.Four sleeping premature infants were measured, born at a median gestation of 28 weeks and 5 days, and imaged at a median gestational age of 34 weeks 5 days. Optical fibers were positioned on the head using a custom-made foam-lined semi-rigid cap, shaped to allow careful but stable positioning of many fibers over the full cranium.217 Since time-resolved data can be slow to acquire, motor stimulus data were collected one source at a time for up to 12 source positions. Each source was illuminated, and the time-resolved transmitted light signal then collected at up to 32 detector positions for 10 sec while the infant rested. The left arm was then gently moved back and forth by an adult for 15 sec, with transmitted light being measured for the last 10 sec of this exercise to capture the resultant hemodynamic response. The next source position was then illuminated, and the 10 sec rest/15 sec exercise cycle repeated. The resulting time-resolved data were Fourier transformed and the log amplitude and phase components extracted at 100 MHz. The difference in the log amplitude and phase measurements between the active and resting states effectively cancels out errors relating to optical coupling and the system’s impulse response function.217,225,226 The resulting differential measurements were then reconstructed into a 3-D map of absorption changes. In this case, a linearized model similar to Eq. 7 was utilized with careful normalization and scaling, and compensation for errors in optical coupling incorporated into Tikhonov inversion.237 The Jn.m r functions were generated using a finite element model implementing the diffusion equation in a homogenous 3-D volume correspondJournal of Biomedical Optics ing to the contours of the infant’s head. The MONSTIR system and finite element algorithm TOAST have been used for a number of other 3-D optical tomography imaging studies on premature infants.145,223,238 3.4 Methodological Extras 3.4.1 Cytochrome oxidase Significant efforts have been made to topographically measure changes in a third intrinsic brain chromophore called cytochrome oxidase. Cytochrome oxidase is a mitochondrial enzyme that plays an important role in oxygen metabolism in the brain. Cytochrome does not change its concentration, but changes its oxidation state and absorption spectrum during the production of adenosine triphosphate ATP . The concentration of cytochrome oxidase in the brain is lower than that of hemoglobin 2 to 3 micromolar compared to 40 to 60 micromolar , so based on published spectra, the peak absorption in brain of reduced cytochrome C would be 0.06 cm−1 at 550 nm, changing to 0.017 cm−1 when oxidized, and 0.00064 to 0.0002 cm−1 at 800 nm compare to equivalent hemoglobin spectra in Fig. 1 .239 It is therefore very desirable, but quite difficult to measure changes in the redox state of cytochrome oxidase in vivo. Nevertheless, measurement methods have been studied carefully240–242 and NIRS-based cytochrome measurements have been reported in adults.209,243 Exogenous dyes for noninvasive imaging. Exogenous contrast agents can also be used for noninvasive brain imaging. Indocyanine green ICG is an FDA-approved dye commonly used for liver and cardiac function tests.244 ICG has strong NIR absorption and fluorescence, and does not cross the blood brain barrier, and so has also found several optical brain imaging applications. A bolus of ICG can be used to quantify blood flow in the brain, by observing the rate of arrival and accumulation of the dye with NIRS measurements.245 Changes in the dynamics of a dye bolus in patients with perfusion abnormalities such as thrombolitic stroke have also been studied using NIRS.246 ICG has been used to simultaneously determine cerebral blood volume and total circulating blood volume in neonates undergoing blood transfusions.247 An interesting recent study reported that it was possible to also observe the fluorescence of ICG in the brain using noninvasive NIRS.248 This result is important as it supports the possibility that more specifically targeted longwavelength dyes could feasibly be measured noninvasively in the brain. Such molecular probes may have sufficient specificity to allow localization or even diagnosis of diseases such as brain cancer and Alzheimer’s.139,249 Noninvasive neuronal imaging via fast scattering. Since it is not possible to use voltage-sensitive dyes in humans, researchers are highly motivated to isolate a strong source of contrast that can directly report changes in neuronal activity. Observation of the fast scattering signal in exposed rodent cortex was described in Sec. 2.1.4. Noninvasive measurement of this signal is highly challenging, as scattering changes are difficult to detect accurately and separate from the effects of absorption. Recent studies of the fast scattering signal have shown compelling results, and compared findings to EEG and fMRI.161,250,251 However, it is still questioned whether the fast signal is sufficiently robust for routine measurement.104,105 Diffuse correlation spectroscopy. A relatively recent addiSeptember/October 2007 051402-21 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… tion to optical topography has been the development of diffuse correlation spectroscopy DCS , which allows noninvasive measurement of cortical blood flow. The theoretical basis of this technique is similar to that of speckle-flow imaging described in Sec. 2.1.2, but DCS measures the temporal autocorrelation of laser light with a long coherence length that has scattered through tissue containing moving red blood cells. For human imaging, DCS acquires signals through single mode fibers placed on the scalp.252–254 Promising results in humans have been reported,255 and it has been shown that it is possible to extend the technique to allow 3-D imaging.256 Summary We have reviewed a wide range of approaches to optical brain imaging, and explored their many important applications to neuroscience. We have provided insights into the fundamental basis of many optical imaging techniques, and highlighted important considerations for their practical implementation. We hope that these details will enable those interested in adopting or developing optical imaging methods to understand the key aspects that should be considered when acquiring measurements or analyzing data. We have sampled the huge body of literature that has successfully used optical brain imaging in all its forms for ground-breaking and fundamentally important research and applications. It is clear that optical brain imaging is a rapidly expanding field, continually evolving to embrace new technologies, contrast mechanisms, and in-vivo applications. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Journal of Biomedical Optics SPIE

Optical brain imaging in vivo : techniques and applications from animal to man

Hillman, Elizabeth M. C.
Journal of Biomedical Optics , Volume 12 (5) – Sep 1, 2007
Download PDF
28 pages

Loading next page...
 
/lp/spie/optical-brain-imaging-in-vivo-techniques-and-applications-from-animal-M50Kl6sITf

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

Publisher
SPIE
Copyright
Copyright © 2007 Society of Photo-Optical Instrumentation Engineers
ISSN
1083-3668
eISSN
1560-2281
DOI
10.1117/1.2789693
pmid
17994863
Publisher site
See Article on Publisher Site

Abstract

other imaging modalities, optical contrast is often a natural feature of many small, biocompatible molecules. Additional advantages of optical imaging include its use of nonionizing radiation, and that instruments are generally inexpensive compared to other clinical imaging modalities. The major challenge of optical imaging is to overcome the effects of light scattering, which limits penetration depth and achievable imaging resolution. “Optical imaging” therefore encompasses a very wide range of measurement techniques, each using different methods to overcome the effects of scatter. These range from laser scanning microscopy of submicron structures, to diffuse optical tomography of large volumes of tissue. Functional brain imaging has also seen dramatic growth in recent years, impacted significantly by the availability of functional magnetic resonance imaging fMRI since 1990.6–9 As a result, new applications of functional brain imaging in humans and animals are continually being established, including neurosurgical planning, the investigation of the physiological basis of neurological diseases such as epilepsy, Alzheimer’s, and stroke, the development of diagnostic methods, drugs, treatments, and interventions, and the study of cognitive and perceptual responses and developmental changes.9 Positron emission tomography PET and x-ray 1083-3668/2007/12 5 /051402/28/$25.00 © 2007 SPIE Journal of Biomedical Optics 051402-1 September/October 2007 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… computed tomography x-ray CT are also valuable and established tools for brain imaging. The obvious advantage of optical imaging over other modalities is reduced cost and infrastructure requirements such as shielded rooms and synchrotrons . However, a much more important distinction is that optical imaging offers such a broad range of contrast mechanisms. While fMRI, PET, and x-ray CT are more established clinical modalities, all suffer from limited functional contrast. For example, to visualize regions of the brain exhibiting functional changes in response to stimuli, fMRI typically exploits the blood oxygen level dependent BOLD signal, widely thought to correlate with HbR concentration. However, a BOLD increase an HbR decrease may correspond to an increase in oxygenation, or a decrease in blood volume. Newer fMRI methods such as arterial spin labeling ASL can provide measures of flow, and intravenous contrast agents such as monocrystalline iron oxide nanocompounds MION can allow cerebral blood volume CBV changes to be observed.10,11 However, these measurements cannot easily be acquired simultaneously, and only the combination of these methods can provide insight into the true hemodynamic responses and oxygenation dynamics in the brain.12 PET is performed using contrast agents such as 2-fluoro2-deoxy-D-glucose containing radioactive isotopes of carbon, nitrogen, fluorine, or oxygen. Oxygen and glucose utilization in the brain can be imaged via the localization of these isotopes. However, resolution is typically poor, data acquisition is slow and hence dynamics cannot be readily evaluated , contrast agents must be manufactured locally prior to each scan, and radioactive dose limits its use in certain patient populations, often including infants. X-ray CTs as well as structural MRI have excellent spatial resolution and can be used to detect morphological changes in the brain, for example following stroke or to delineate tumors. However, CT typically cannot elicit information about the function of the tissue, and so potentially reversible changes preceding structural changes cannot be detected. Therefore, the real advantage of optical imaging over these techniques is its ability to detect and image a wide range of functional parameters. For clinical applications, noninvasive optical imaging can provide complimentary information to other modalities such as fMRI and provide a low-cost alternative in some cases. This is in addition to serving populations often unable to receive MRI or PET scans such as young infants or the critically ill. Clinical optical brain imaging is generally noninvasive and uses near-infrared NIR light to obtain improved penetration through the scalp, skull, and brain see Fig. 1 a . In addition, optical brain imaging is finding widespread applications as a research tool for both clinical and animal studies of basic brain function and diseases. Still, so little is known about the way that the normal brain functions, in part due to the difficulties of measuring such a complex organ without disturbing or damaging its in-vivo function. Optical imaging allows the living brain to be closely observed, and many functional interactions and changes to be investigated over many length scales. The effects of pathologies and treatments can also be readily investigated using fairly basic optical imaging tools. In animals, it is also relatively simple to expose the cortex, allowing high-resolution direct optical imJournal of Biomedical Optics aging of the brain’s surface. Extrinsic dyes and visible wavelengths can also be used see Fig. 1 b . In summary, optical imaging’s key advantage is the ability to measure a range of functional contrast, can readily be exploited in functional brain imaging via a wide range of approaches, from animal studies of the intricate cellular mechanisms of normal and diseased brain, to in-vivo noninvasive clinical brain imaging. In this review, we describe a variety of optical imaging approaches to functional brain imaging. The basic principles of each technique are summarized, followed by examples of use for in-vivo brain imaging. Section 2 details the use of optical techniques for exposed cortex imaging. Section 3 describes noninvasive clinical optical imaging of the living brain. Exposed-Cortex Imaging and Microscopy The use of optical imaging in animal studies is widespread. Animal imaging has become more than just a step toward clinical diagnostic imaging, but a tool with which to learn more about the basic mechanisms of brain function as well as disease processes and the effects of treatments. These findings can be utilized in their own right to aid in the development of drugs and treatment. They can also contribute to the interpretation of results from other imaging modalities such as fMRI. Recent applications of animal imaging have included studies of stroke acute and chronic ,13 Alzheimer’s disease,14 migraine,15 epilepsy,16 the effects of cocaine,17 and the mechanisms of neurovascular coupling.18,19 Imaging animals rather than humans provides significantly more flexibility, since preparations can be much more controlled and diseases and treatments can be systematically compared. Extrinsic dyes and cross-validation techniques such as electrophysiology can be utilized and developed without the need for clinical regulatory approval. The adverse effects of light scattering worsen as the size of the tissue being interrogated increases. Imaging smaller animals therefore also offers significant technical advantages, allowing higher resolution imaging and improved sensitivity and quantitation. For very high resolution imaging, the cortex can be surgically exposed, allowing direct optical imaging of the brain’s surface with only minimal disturbance to normal brain activity. Exposed cortex optical imaging is highly versatile, and most commonly performed in animals, although it has also been performed on the human brain, for example, during epilepsy surgery.20–22 Next, we provide an overview of the most established exposed-cortex, in-vivo optical neuroimaging techniques ranging from 2-D and 3-D macroscopic imaging to microscopic in-vivo imaging using two-photon microscopy. Basic principles and methods are first described, followed by specific examples of neuroimaging applications. 2.1 Macroscopic Optical Imaging of Exposed Cortex in Two Dimensions 2.1.1 Experimental method During surgical preparation for exposed-cortex imaging, the animal commonly mouse, rat, cat, or primate is anesthetized while the scalp is retracted and the skull is gently thinned over the brain area of interest. In some cases, this piece of skull September/October 2007 051402-2 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 1 Intrinsic chromophores and structures of the brain. a For small animal imaging, it is possible to use visible light to look at oxygenationdependent hemoglobin absorption as well as exogenous absorbing and fluorescing dyes . The cortex can be exposed and imaged at high resolution, or imaged through the intact skull and scalp at lower resolution. b For human brain imaging, near-infrared light will penetrate more readily through the scalp and skull to sample the brain thanks to lower scatter and absorption. The human brain has a more complex convoluted structure compared to lower mammals. Spectra: major chromophores in brain are oxy- and deoxyhemoglobin and water lipid is omitted here . For generation of these spectra: water content is assumed to be 90% quoted value in newborn brain, 71 to 85% in adult brain34 . Hemoglobin absorption is shown assuming 2-mM concentration of hemoglobin in blood and 3% content of blood in tissue 60 mM in tissue34 calculated from spectra.257 HbR spectrum represents 3% fully deoxygenated blood in tissue. Scatter spectrum is approximate using: s = A −b, where A = 1.14 10−7 m and b = 1.3.228 may be removed to provide a direct view of an area of the brain’s surface. In mice, the brain can often be visualized without even thinning the skull. Methods have been carefully developed to ensure that normal brain function is not affected by these approaches. Most often, the animal will remain anesthetized during measurements.23 In some cases, however, a permanent window can be surgically implanted and the animals can survive for over a year while an area of their brain is accessible for repeated direct optical imaging.14,24 Imaging of pathologies requires only that the animal be anesthetized. Images during acute interventions, such as adJournal of Biomedical Optics ministration of drugs or introduction of a stroke, must be carefully synchronized with image acquisition. While the brain’s functional response to stimulus is being studied, a carefully controlled stimulus is presented during imaging. Such stimuli may include visual e.g., showing patterns or pictures , somatosensory e.g., electrical or tactile stimulation of the whiskers or forepaw , or auditory stimuli different frequencies and amplitudes . Stimuli may be repeated in regular blocks, or presented in a randomized pattern.18,25 Data are usually shown as the average response to multiple stimulus repetitions to obtain September/October 2007 051402-3 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 2 A typical setup for camera-based imaging of the exposed cortex. Inset shows typical hemodynamic and calcium-sensitive responses in rat somatosensory cortex to 1-mA forepaw stimulus delivered in 3 msec pulses at 3 Hz for 4 sec. better signal to noise where noise includes normal background physiological variations . 2.1.2 Basic principles: two-dimensional camerabased imaging 2-D camera-based imaging of the exposed cortex is the most simple of optical imaging techniques and is often referred to as “optical imaging” in the neuroscience literature. We describe three approaches to 2-D imaging of the cortex: 1. intrinsic signal or hyperspectral imaging of oxy- and deoxyhemoglobin dynamics; 2. imaging of extrinsic voltage sensitive dyes providing sensitivity to neuronal activity ; and 3. speckle flow imaging, capable of imaging the blood flow dynamics in the superficial cortex. Two-dimensional intrinsic and hyperspectral imaging. Figure 2 shows a typical setup for 2-D hyperspectral imaging of exposed rat cortex. This neuroimaging approach relies on the distinctive absorption spectra of oxy- and deoxyhemoglobin HbO2 and HbR as shown in Fig. 1 a . The hemoglobin in blood is the most significant absorber in the brain at visible and NIR wavelengths, and functional activity in the brain results in changes in blood flow, blood volume, and oxygenation.26. Blood delivers oxygen to tissues by locally converting HbO2 to HbR, so changes in the oxygenation state of blood correspond to changes in the relative concentrations of HbO2 and HbR. Consequently, variations in the amount and oxygenation state of hemoglobin modulate the absorption properties of the brain. By simply shining light into the exJournal of Biomedical Optics posed cortex and taking pictures with a camera, the hemodynamic response to stimulus can be readliy observed. Further, since HbR and HbO2 have different absorption spectra see Figs. 1 a and 1 b , measurements with different wavelengths of light can produce images that are preferentially sensitive to changes in the concentration of either HbO2 or HbR. At isobestic points, where the HbR and HbO2 absorption is the same e.g., 500, 530, 570, and 797 nm , changes in total hemoglobin concentration HbT= HbR+ HbO2 can be measured independently of changes in blood oxygenation. Early approaches to optical imaging of the exposed cortex referred to “intrinsic signal” imaging, and often used illumination at an isobestic point to examine localized changes in tissue HbT.27 By adding a second wavelength often 610 nm , researchers were able to investigate oximetric changes in addition to changes in HbT.28,29 Since HbO2 absorption at 610 nm is low around 16% of the absorption of HbR , changes in absorption at 610 nm are largely attributed to, and interpreted as, changes in HbR. However, it should be noted that even though HbO2 absorption at 610 nm is lower than that of HbR, HbO2 changes can interfere with and distort inference of HbR concentration changes from measurements at 610 nm. To calculate the actual changes in concentration of HbO2, HbR, and HbT, the absorption of the brain must be measured at two or more wavelengths hyperspectral imaging . In some cases, multiple wavelengths across the visible spectrum are used. This is typically achieved by synchronizing rapid camera acquisition with a filter wheel in front of a white light September/October 2007 051402-4 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… ln I ln I ln I DPFx = DPFx DPFx 1 HbO2, 1 2 HbO2, 2 DPFx DPFx DPFx 1 HbR2, 1 2 HbR2, 2 n HbO2, n n HbR2, n cHbO2 cHbR I = exp − I0 aDPFx +G , where I is the detected intensity, I0 is the incident intensity, DPFx is the mean pathlength traveled by light in the tissue a function of baseline scattering and absorbing properties , G is a geometry-dependent factor, and a is the absorption coefficient of the tissue as given by Beer’s law: nc n , i.e., the absorption coefficient is a linear sum of the absorption of the n chromophores in the tissue, where n is the wavelength-dependent specific absorption coefficient of a particular chromophore e.g., the HbO2 absorption spectrum and cn is its concentration. Note that “differential pathlength factor” DPF is a term used to describe the approximate multiplicative increase in pathlength experienced by a photon traveling in a scattering medium between two points a linear distance x apart. Thus, the true pathlength is DPFx. I0 may be spatially varying and difficult to accurately determine. G is also very difficult to quantify. Therefore, it is common to consider changes in detected signal relative to an initial time point t = 0 or a reference state: ln It I t0 =− t − a t0 DPFx ln I = − aDPFx. If measurements at two wavelengths 1 and 2 are acquired, it is possible to use these differential measurements to determine the changes in HbO2 and HbR concentration occurring relative to time t = t0 using: ln I ln I =− =− HbO2, 1 cHbO2 + cHbO2 + HbR, 1 cHbR DPFx HbO2, 2 HbR, 2 cHbR DPFx 2 . 4 Where more than two wavelengths are available, HbO2 and HbR can be calculated using a linear model: Journal of Biomedical Optics While imaging absorption changes in superficial tissues appears quite straightforward, the effects of scatter must be carefully considered. Since scattering is wavelength dependent in tissues an approximate scatter spectrum is shown in Fig. 1 a , each measurement wavelength will experience different levels of scattering. The mean path of light through tissue DPFx will also be altered by the background absorption properties of the tissue, which are also heavily wavelength dependent, owing to the presence of hemoglobin, so DPFx 1 DPFx 2 in almost all cases. Therefore, as given by Eqs. 4 and 5 , accurate estimates of these wavelengthdependent pathlength factors DPFx n are required even to calculate differential changes in hemoglobin concentration cHbO2 and cHbR . Note, however, that if I t / I t0 is measured at an isobestic point, and DPFx does not change over time, then I t / I t0 will be directly proportional to HbT t . Estimation of this wavelength-dependent pathlength DPFx is difficult. For exposed-cortex imaging, a Monte Carlo model of light propagation is generally employed.35 Such a model requires estimates of the background absorption and oxygenation level of the brain tissue to predict the likely scattered paths of photons.31 For camera imaging, the mean pathlength calculation should simulate a widely distributed illumination, and incorporate the numerical aperture of the camera’s lens. It has been demonstrated that incorrect estimates of DPFx can directly impact the correct calculation of HbO2, HbR, and HbT changes in tissue.35 Note also that if only values of HbO2 t and HbR t are derived, that it is not possible to infer changes in absolute oxygen saturation SO2 = HbO2 t / HbO2 t + HbR t − HbO2 t0 / HbO2 t0 + HbR t0 without first estimating the baseline absolute concentrations of HbO2 and HbR. Another important issue to be considered in 2-D imaging of the cortex is the depth sensitivity of the measurement. The detected light has scattered within the tissue, and reports absorption that it has experienced along its path. The majority of detected photons have predominantly sampled the very superficial layers of the cortex, whereas fewer have traveled deeply and managed to return to the surface. This means that the image of the cortex is a superficially weighted sum of signals from different layers. Since the cortex has a complex structure of large superficial vessels and deeper capillary beds, the relative contributions of these compartments to the overall signal can be difficult to interpret. Additionally, since the mean penetration depth of the light will also depend on its wavelength, all wavelengths have not necessarily interrogated exactly the same region of tissue, which can impact the quantitative accuracy of estimates of HbO2, HbR, and HbT. September/October 2007 051402-5 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… A typical hemodynamic response in the somatosensory cortex of a rat during a four-second forepaw stimulus 3 s pulses at 3 Hz is shown inset in Fig. 2. The hemodynamic response to stimulus is fairly slow, peaking several seconds after the first stimulus pulse. A large increase in HbO2 is generally seen, along with an increase in HbT and a decrease in HbR. The overall timing characteristics of each hemoglobin type are different from each other, and contain information about the underlying oxygen delivery and consumption dynamics.19,36,37 Some researchers routinely observe a short initial increase in HbR, attributed to rapid oxygen consumption at the site of neuronal activation.29,31,33 The origin of this initial dip is somewhat contentious, and it is discussed further in Sec. 2.1.3. The data shown in Fig. 2 are an average response of 400 stimulus repetitions. It should be noted that individual stimulus responses are often affected by physiological noise, such as variations in baseline blood flow to the region and fluctuations in systemic parameters such as blood pressure. Despite all of these important considerations, 2-D optical imaging of the exposed cortex has provided very significant and well validated observations of functional cortical hemodynamics. Specific examples of applications of these techniques are given at the end of this section. Two-dimensional voltage sensitive dye imaging. As described before, intrinsic signal and hyperspectral imaging can report the hemodynamic response to functional activation in the brain. This hemodynamic response is generally assumed to spatially correlate with the underlying brain activity corresponding to neurons firing and conducting electrical signals. The interrelation between neuronal and hemodynamic activity is currently an area of intense research.18,38,39 This neurovascular coupling is thought to be an essential part of normal brain activity and is believed to be disturbed by pathologies such as Alzheimer’s disease and stroke.40,41 Understanding whether it is possible to determine detailed information about neuronal activity from hemodynamic measures is also an important area of research, driven by the fact that clinical fMRI is currently only sensitive to hemodynamic changes. Voltage sensitive dyes VSDs provide a way to optically detect neuronal activity in the in-vivo exposed-cortex. VSDs were first used in 196842 but received renewed attention in 1986 when they were demonstrated for in-vivo mammalian brain imaging by Grinvald and colleagues43–45 studies that predate the first reported intrinsic signal measurements.27 Since then, significant advances in digital detector technology, dye formulation, and staining strategies have led to established use of VSDs to visualize neuronal activity in-vivo2,24,46,47 and also for mapping electrical activity in the heart.48,49 VSDs are molecules that generally bind across a neuron’s membrane. Changes in membrane potential cause measurable changes in the fluorescence of the dye.50 Since a large percentage of the total membrane area is dendritic, VSD signals reflect both spiking and synaptic activity.42 Calcium sensitive dyes CaSDs are compounds that increase their fluorescence in response to increases in calcium ion concentration. When cells such as neurons are loaded with CaSDs, it becomes possible to monitor intracellular calcium.51 Since a calcium influx occurs when a neuron fires an action Journal of Biomedical Optics potential, the calcium dye’s fluorescence can optically report this, providing a second tool with which to optically image neuronal activity. To image these dyes, a similar system to that shown in Fig. 2 is used, with the filter wheel replaced by an excitation filter of the appropriate wavelength. A long-pass filter is placed in front of the camera to block excitation light and isolate the fluorescence emission. In some cases, a dichroic mirror is used such that the illumination light is reflected onto the brain, and emitted fluorescence passes through the dichroic to the camera via a secondary emission filter to remove residual excitation light . This is the configuration of a commercially available VSD imaging system sold by Optical Imaging Incorporated Rehovot, Israel .2 As shown in the inset of Fig. 2, the hemodynamic response to stimulus typically peaks 2 to 4 sec after the stimulus begins. The neuronal response is much faster. Measured using VSDs and CaSDs, neuronal responses typically peak within 25 to 100 msec, and are usually followed by an inhibitory or refractory period of negative signal, which returns to baseline at around 150 to 300 msec after the stimulus starts.52 The black trace in Fig. 2 shows the fluorescence changes seen in the somatosensory cortex stained with Oregon Green 488 BAPTA-1, AM calcium sensitive dye Invitrogen during the 4-sec forepaw stimulus 3- s pulses at 3 Hz . Very fast, individual responses for each stimulation pulse can be seen, in contrast to the very slow hemodynamic response. To capture these responses, a faster camera is required than for hyperspectral or intrinsic signal imaging preferably 10 ms per frame . Voltage sensitive dyes typically have quite low contrast typical F / F 0.5% . As a result, they can be affected by changes in hemoglobin absorption due to heart rate, breathing, and of course functional activation. Two methods have been developed to overcome this. 1. New dyes have been developed that excite and emit at higher wavelengths, corresponding to regions where hemoglobin absorption is less significant see Fig. 1 a . 2. Acquisition of data is often triggered on the heart beat and sometimes breathing of the animal.2 Many repetitions are generally averaged, and “blank” acquisitions where no stimuli are presented are interleaved such that the final signal is given by the stimulated – blank data. Voltage sensitive dyes can only be used in animals at present, and are usually topically applied directly to the exposed cortex for 1 to 2 h before imaging commences.53 Some VSDs have been shown to quickly photobleach or have phototoxic effects, and so illumination times are typically kept to a minimum. It is also important to note that the change in fluorescence or absorption of VSDs is usually a manifestation of a spectral shift in either the absorption or emission spectrum of the dye.54 This means that measurement wavelength selection should not necessarily be based on the peak absorption and emission wavelengths, nor should broadband emission filters spanning the whole emission peak always be selected. The excitation and emission wavelengths should be selected based on areas of the spectrum where the largest change in intensity will result when a spectral shift occurs e.g., a band on the falling side of the emission spectrum .54 Note that selection of wavelengths may therefore determine whether an increase or deSeptember/October 2007 051402-6 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… crease in absorption or fluorescence is observed in response to an increase in membrane potential. CaSDs have much higher contrast than VSDs F/F 5 % and hence, despite generally exciting and emitting in the visible spectrum, are not strongly affected by small hemodynamic changes although the small dip after the calcium response in Fig. 2 may be due to the concurrent increase in HbT and HbO2 absorption . The challenge with CaSDs is introducing them into the neurons successfully, which to date has required the addition of dimethyl sulfoxide DMSO and pluronic acid, followed by direct pressure injection via micropipette into the cortex and incubation of 1 / 2 to 1 h.3,55 However, it is possible to achieve fairly uniform loading of neurons within an area of several millimeters of cortex, and to observe responses similar to those seen with voltage sensitive dyes. Unlike VSDs, CaSDs tend to increase the intensity of their fluorescence emission in the presence of calcium, without significant shifts in the spectral position of their emission and excitation bands. Examples of the use of VSDs and CaSDs for neuroscience studies are given at the end of this section. Imaging of CaSDs with two-photon microscopy is described in Sec. 2.2 Laser speckle-flow imaging. Another recent development in 2-D optical imaging of the exposed cortex is laser speckleflow imaging. While hyperspectral and intrinsic imaging provide measurements of the changes in HbO2 and HbR concentrations, they do not explicitly measure changes in the flow rate of the blood. Flow is a very important parameter if the rate of oxygen delivery and/or consumption is to be calculated. Laser Doppler has been widely used to provide point measurements of blood flow in the brain. However, laser speckle-flow imaging of the exposed cortex can provide a spatially resolved image of the flow throughout the vascular network during functional activation.56,57 Acquiring laser speckle-flow imaging data uses a similar setup to that shown in Fig. 2. However, in this case, the light source is replaced with a divergent laser diode typically at 780 nm . The camera then simply acquires rapid and highresolution images of the exposed cortex, which appear to be low-contrast images of the speckle pattern from the laser illumination. The speckle pattern is caused by the coherent laser light scattering within the brain. If red blood cells are moving within the image, they cause the speckle pattern to vary over time. The rate at which the speckles change relates to how fast the red blood cells are moving. Laser Doppler measures this temporal modulation. However, within the camera image, this intensity fluctuation causes spatial blurring. By calculating the spatial uniformity of an area of camera pixels e.g., 5 5 or 7 7 , the amount of blurring can be quantified from the pixel group’s standard deviation of intensity / mean intensity for each time point . This speckle contrast can be mathematically related to the correlation time, which gives a measure of the velocity of flow of the red blood cells.56 Laser speckle-flow imaging has been performed simultaneously with hyperspectral imaging of the exposed cortex by adding a second camera, a dichroic, and filters to block the laser light . This has enabled the cerebral metabolic rate of oxygen consumption to be calculated.30,36,37 Examples of the use of speckle-flow imaging are given next. Journal of Biomedical Optics 2.1.3 Imaging of exposed cortex: examples Exposed-cortex imaging has been applied to a wide range of research areas. We can broadly summarize these as: 1. functional imaging to investigate the sensory and cognitive processing functions of the brain, 2. functional imaging to improve understanding of the basic mechanisms of the hemodynamic and neuronal response to stimulus, and 3. study of the effects of diseases and treatments on the normal behavior of the brain. Next we describe a selection of recent studies applying exposed-cortex optical imaging techniques within these important brain imaging fields. Functional imaging for sensory and cognitive processing research. Most regions of the brain’s cortex appear to be dedicated to consistent tasks. Particularly in the sensory cortices, there is even spatial correspondence between, for example, the location of two whiskers on the snout of a rat and the two sites or barrels on the somatosensory cortex that respond when each whisker is flicked.58,59 Different parts of the auditory cortex respond to different frequencies of sound.60 Studies of these systems are prevalent throughout the functional imaging literature, but optical imaging of the exposed cortex is currently the only technique that can provide sufficient resolution to reveal the intricate nuances of these hemodynamic and neuronal responses. Figure 3 a shows a recent exposed-cortex optical imaging result revealing the spatially distributed cortical response to visual stimulus in anesthetized rat.61 Cortical imaging data were acquired at 590 nm and the results reveal a very clear spatially correlated retinotopic map of the cortical responses to visual stimuli across the entire visual field of one eye. The purpose of this study was to determine whether it is feasible to use rats as a model to evaluate adaptation of visual cortical function in response to retinal degeneration, and for development of associated potential therapies. Figure 3 b shows an example of a highly complex study of the visual cortex of awake monkeys using both intrinsic signal hemoglobin and VSD imaging.24 Three Macaca fascicularis monkeys had ‘windows’ surgically implanted into their skulls over their visual cortex covering areas V1, V2, and V4 . The cortical response to visual stimulus was then imaged repeatedly for up to 1 year using VSD RH-1691 Optical Imaging Incorporated, Rehovot, Israel to stain the cortex 1 to 3 times per week excited at 630 nm± 5 nm, emission: 650-nm long-pass . For studies of the visual cortex, it is very common to utilize two equivalent visual stimuli and display the difference between the two responses. This is because each response on its own consists of many different changes in the vasculature and tissue, so by subtracting two or more equivalent states, only the very specific difference between the stimuli responses should be revealed. In this case, Fig. 3 b shows intrinsic signal ocular dominance maps 605-nm reflectance generated by showing the monkey a pattern with either his left or right eye blocked the baseline image of the exposed cortex is also shown . The difference between these two states therefore reveals the response of one eye as bright areas and the response of the other as dark stripes. The image shows that this part of the visual cortex processes similar areas of the image from each eye in close proximity to each other, interleaving them, rather than processing the full imSeptember/October 2007 051402-7 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 3 Examples of in-vivo exposed cortex imaging. a Visual stimuli were presented in different parts of a rat’s field of view, and the corresponding cortical hemodynamic HbT responses were mapped color code indicates corresponding location of visual stimulus . Reproduced with permission from Gias et al.61 b A cranial window was implanted to allow direct imaging of primate visual cortex. A baseline image of the exposed brain is shown above an intrinsic signal ocular dominance map, acquired by subtracting the response from visual stimulus to the right eye from the left eye. Chronic VSD imaging was also performed, and ocular dominance maps closely resembled intrinsic signals. Time courses of the voltage-sensitive signal red showed a more prompt onset that the hemodynamic response blue . Reproduced with permission from Slovin et al.24 c An optical fiber bundle was used to image the exposed, VSD stained cortex of an awake mouse while it explored its surroundings. Image series shows the mouse’s whisker gradually brushing against an obstacle and the corresponding cortical response. Reproduced with permission from Férézou et al.65 d A fiber optic bundle was used to image the exposed cortex of a rat undergoing electrical whisker stimulation during simultaneous acquisition of fMRI at 4.7 T. The hemodynamic response shows localized cortical changes in HbR, HbO2, and HbT, which correspond well to oblique-slice fMRI BOLD response study by Hillman, Devor, DeCrespigny and D’Arceuil . e Baseline speckle contrast image showing higher flow as darker contrast and strongly accentuated vessels. Image series from top to bottom show flow prior to induction of spreading depression, 2 min after which a wave of cortical flow increase moves across the field of view, and 20 min after, when hyperperfusion in the dural vessels is prominent. Plot shows time course of flow for region of interest in cortex blue and dural vessel red . Reproduced with permission from Bolay et al.15 f left Cross sectional x-z images through the cortical HbO2 response to a 4-sec forepaw stimulus in rat imaged using laminar optical tomography, top transecting a superficial vein, and bottom transecting the deeper capillary response. right 3-D rendering of arterial, capillary, and venous compartment hemodynamic responses. Ex-vivo two-photon microscopy of vasculature inset agrees with compartment discrimination. Reproduced with permission from Hillman et al.19 Journal of Biomedical Optics 051402-8 September/October 2007 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… ages from the left and right eyes separately. Very similar maps were seen for both VSD and intrinsic signal imaging, although the timing and evolutions of the neuronal and hemodynamic signals were quite different. The data shown here are just the introduction to a very broad study where similar maps were used to investigate the visual cortex’s response to rapid saccadic eye movements. There are numerous other chronic optical imaging studies of monkey visual cortex62 and monkey somatosensory cortex.63 Very recent studies have also demonstrated that plasticity of ocular dominance can even be studied in mice using exposed-cortex optical imaging.64 Another recent novel approach to exposed-cortex optical imaging is shown in Fig. 3 c . While most rodent measurements use anesthetized animals, some studies of native behavior require that the animal be both awake and able to freely move around. Férézou, Bolea, and Petersen developed a technique where a mouse’s voltage sensitive dye-stained exposed somatosensory cortex could be imaged via a flexible fiber bundle while the awake mouse explored its surroundings.65 By simultaneously acquiring video of the mouse’s behavior, they were able to directly image the neuronal response in the brain when the mouse’s whisker contacted an object. Functional imaging to investigate the mechanisms of activation. The functional imaging studies described above exploit the fact that colocalized neuronal and hemodynamic responses can be detected in the brain in response to stimuli.26 However, very little is understood about how and why these functional responses actually occur.40,41,66 Understanding the origins and characteristics of these responses, and the interplay between neuronal and hemodynamic processes, is key to developing an improved picture of healthy and impaired brain function. It is also an essential part of developing and improving all functional brain imaging modalities that rely on the contrast provided by hemodynamic changes such as fMRI. Since exposed cortex optical imaging can provide such a close view of the working brain, as well as access for invasive measurements such as electrophysiology, it has become an established tool for the study of basic functional mechanisms. A wealth of recent studies have investigated so-called neurovascular coupling NVC relationships using optical methods.18,31,38,39,58,67–69 By quantifying the hemodynamic response via optical imaging and/or laser Doppler simultaneously with the neuronal response via electrophysiology , these studies sought to identify a parametric relationship that could be used to predict one response from the other. Devor et al. demonstrated in rat that within a single whisker barrel, the neuronal response to different amplitudes of tactile whisker stimulus saturated before the hemodynamic response.18 More recently, subtleties of the spatiotemporal evolution of both responses have shown that the NVC relationship is highly complex.58,70 A recent study also demonstrated that there can be significant differences between the responses of anesthetized and awake rodents.23 An alternative to attempting to infer neuronal responses from BOLD-type measurements is to instead locate sites where oxygen consumption is increasing by determining the cerebral metabolic rate of oxygen consumption CMRO2 .71,72 In principle, these regions should spatially correlate more closely with neuronal activity than the BOLD response, which can be significant even in distant draining Journal of Biomedical Optics veins.73 CMRO2 is a parameter that can only be calculated by considering the amount of oxygen that is delivered to and extracted by the tissue. Although 2-D hyperspectral imaging can provide changes in HbO2, HbR, and HbT concentration, the data represent discrete snapshots of the current oxygenation state of the blood, and not the rate at which the blood is being replaced. It is therefore necessary to also measure blood flow, to determine the rate of delivery of oxygen to the tissue. Other estimates are also required such as the variations originating in the venous compartment. Both laser Doppler and more recently speckle-flow imaging have allowed simultaneous measurement of oxygenation, blood volume, and flow dynamics in the exposed cortex,30,36,37 thereby allowing direct calculation of CMRO2. These calculations are based on models of the hemodynamic response to stimulus.32,74–78 Interestingly, in-vivo optical imaging and microscopy of the cortex are also allowing the vascular mechanisms of the hemodynamic response underlying these vascular models to be investigated directly.19,79,80 Another approach to comparing exposed-cortex optical imaging findings to fMRI is to perform fMRI studies and optical measurements on the same subjects in response to the same stimulus. Several recent studies have demonstrated simultaneous optical imaging of the exposed cortex during fMRI acquisition.81 An example of such a study is shown in Fig. 3 d , where an MRI-compatible fiber optic endoscope was configured to allow synchronized fMRI and functional optical imaging at two wavelengths 570 and 610 nm during electrical whisker stimulus of varying amplitudes. Similar studies have been performed in humans scheduled to undergo open cortex brain surgery21,82,83: functional responses to voluntary motion of the tongue were recorded in the language area of the brain prior to surgery using fMRI, and then compared to optical imaging results obtained during surgery, where the awake patient was asked to repeat the same stimulus paradigm.83 These results were compared to electrocortical stimulation maps acquired by directly stimulating the cortex and questioning the patient about the resulting involuntary tongue movements. A final optical imaging result of significant importance, and which has prompted much debate, was the discovery of the so-called “initial dip.”33 In some studies, the hemodynamic response to stimulus is seen to begin with an increase in HbR, prior to a subsequent decrease. This finding was interpreted as evidence that neurons initially use up local oxygen more quickly than it is being supplied by the gradually increasing blood flow, an important finding for modeling the mechanisms of the hemodynamic response. It was proposed that this early signal could also allow improved spatial localization to the neuronal response. This initial HbR rise was subsequently also seen in fMRI results as an initial dip in the BOLD signal,84,85 and in many other optical imaging studies.29,37,86 However, the dip is controversial because it is not seen in all studies.87 Reasons proposed for this discrepancy have ranged from differences between species, anesthesia, and paradigms, to measurement techniques and errors in spectroscopic analysis of optical signals especially the wavelength dependence of pathlength88 . An excellent summary of the initial dip can be found in Ref. 89. September/October 2007 051402-9 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Functional imaging to investigate pathologies and treatments. A further use of exposed cortex imaging is to directly investigate the manifestation, characteristics, and treatment responses of diseases and disorders. An example of one such study is shown in Fig. 3 e , which used speckle-flow imaging to investigate migraine in a mouse model.15 Cortical spreading depression CSD is a slow moving wave of neuronal and hemodynamic activity that spreads across the cortex following a small insult to the cortex, such as a pin prick.90 CSD is thought to be the cause of the visual disturbances that often precede a migraine headache. This study demonstrated that speckle-flow imaging could not only see the wave of CSD as a transient increase in parenchymal blood flow, but also revealed a prominent subsequent increase in blood flow in the middle meningeal artery MMA , a vessel that traverses the dura. Such a dural blood flow increase is consistent with the delayed onset of a debilitating migraine headache. Supporting measurements were then made to conclude that CSDs can trigger events leading to headache, identifying new pathways that could potentially be targeted for the development of treatments for migraine suppression. Similar studies using both speckle flow imaging and hypersepctral imaging of the exposed cortex have examined the effects of middle cerebral artery occlusion as a model of stroke.13,91 Both the acute stage and chronic manifestations of stroke can be studied. Function of the area affected by a stroke can be evaluated after days or even weeks of recovery.92,93 The effects of established and experimental treatments can be tested in these models of human stroke, and even compared to the fMRI manifestation of the same outcomes.94 It is also possible to study epilepsy in animal models using exposed-cortex intrinsic signal optical imaging95 and VSDs.96 The hemodynamic response to intrasurgical cortical stimulation during epilepsy surgeries on humans has also been investigated using optical imaging.22 Fast scattering signal. Another source of intrinsic optical signal in the brain is the fast scattering signal. Historically, this signal has been the source of much debate, particularly related to the difficulties in detecting it noninvasively,104,105 although the signal is clearly observable in isolated nerves.106 Recently, exposed cortex measurements during whisker stimulation in rats have provided compelling data revealing localized signals corresponding to regions of neuronal activity.107 Changes in the dark-field reflectance of 660-nm light occurred within 50 ms of stimulation and correspond well with electrophysiological recordings. A 4-mm-diam imaging probe was placed in direct contact with the dura over the exposed somatosensory cortex. The probe consisted of an outer ring of optical fibers delivering illumination light from 660-nm light emitting diodes, and a central fiber optic imaging conduit transmitting remitted light to a sensitive charge-coupled device CCD camera. Several possible mechanisms for the origin of the fast scattering signal have been proposed, including reorientation of membrane proteins and swelling of neuronal dendrites due to water influx in response to ionic currents during depolarization.107 Note that the possibility that fast scattering changes in the brain may distort and interfere with VSD intensity signals has not yet been explored. It is also still not fully understood whether significant scattering changes occur in the brain as part of the slower hemodynamic response. Some researchers have tried to compensate for scattering changes by including a scattering spectrum term in Eq. 2 , as if it represents a third chromophore.33 Care should be taken with this approach, since it does not completely account for scattering effects, and incorporating a third chromophore which may not in fact be changing may adversely affect the fit results of the other two chromophores HbO2 and HbR . Optical coherence tomography of functional brain activation. Optical coherence tomography OCT is a well established technique for high resolution depth-resolved imaging of structures such as the eye.108 OCT overcomes the effects of light scattering to obtain high resolution depthresolved images by using coherence to isolate light that has been directly backscattered. This light therefore has the shortest possible pathlength through the tissue. OCT is also typically performed using NIR wavelengths 1 m. As a result, OCT signals carry very little information about the hemoglobin absorption properties of tissue Eq. 1 and Fig. 1 and cannot be used to measure fluorescence contrast. The contrast in OCT is therefore generally due to refractive index mismatches, and not due to changes in absorption corresponding to hemoglobin concentration or oxygenation. Nevertheless, OCT has recently been applied to brain imaging to explore depth-resolved changes in the cortex of cats and rats during functional stimulus.109–111 To date, the functional OCT changes reported have followed the slower time scale of the hemodynamic response. Both increases and decreases in OCT contrast are seen in close proximity, at times when HbT should be increasing. It is possible that OCT could also detect fast neuronal scattering signals. However, the slower signals observed to date are hypothesized to correspond to changes in the speed and distribution of scattering red blood cells, and changes in diameter and tone of pial vessels and capillaries. September/October 2007 2.1.4 Methodological extras In addition to the three basic 2-D imaging techniques described before, a number of other distinct optical exposed cortex imaging approaches have recently been demonstrated and are described briefly next. Intrinsic fluorophore imaging. Several recent studies have shown compelling evidence that it is possible to measure concentration changes in intrinsic flavoproteins FAD in the brain from their intrinsic fluorescence in vivo.97–100 A system almost identical to that used for VSD imaging is employed, but no contrast agent is added to the brain. Excitation is between 450 and 490 nm and emission is between 500 and 550 nm. FAD provides a measure of metabolic changes in cells neurons and possibly glia , and has a distinct time course compared to shorter neuronal and longer hemodynamic responses. FAD imaging was recently used to explore visual plasticity in mice in a similar way to the studies of the visual cortex described before,101 as well as plasticity in the auditory cortex102 and epileptic activity.103 Since flavoprotein fluorescence requires visible light for its detection, it is unlikely that it will be possible to detect these signals noninvasively see Fig. 1 . Journal of Biomedical Optics 051402-10 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Laminar optical tomography. While 2-D optical imaging of the exposed cortex has found many very valuable applications, there remains concern over the fact that images are a 2-D projection of a combination of signals from superficial vessels and deeper layers of the cortex. The vascular structure of the cortex is organized with large arteries and veins on the surface, and then arteriolar and venous branches diving almost straight down to deeper capillary beds. These capillary beds supply the active neurons, which themselves are also organized in specific layers within the cortex. Therefore, not only are neuronal and vascular signals likely to be quite heterogeneous as a function of depth, but the laminar variations in neuronal and hemodynamic responses are of great interest. Laminar optical tomography LOT is a recently developed imaging technique that can add a depth dimension to hyperspectral imaging of the exposed cortex.19,112,113 LOT utilizes a combination of scanning microscopy instrumentation and optical tomographic image reconstruction techniques to acquire noncontact measurements of the cortex at more than 40 frames per second. It can image to depths of 2 mm with 100 to 200- m resolution. A recent result is shown in Fig. 3 f , where LOT was used to image the exposed-cortex hemodynamic response to rat forepaw stimulus at 473 and 532 nm. The depth-resolved images were validated by comparison to the vascular architecture, and microsvascular density determined via two-photon imaging of fluorescent vascular casts from the same animals. A spatiotemporal delineation technique was also applied that exploited the distinct differences in the temporal evolution of the functional signal in the different vascular compartments: arterioles, capillaries, and veins. Figure 3 f shows two cross sectional x − z images through the HbO2 cortical response, one transecting a superficial vein, and one transecting the deeper capillary response. A 3-D rendering of the three vascular compartments, delineated via their distinctive time courses, is also shown, revealing superficial arterioles and venules, and deeper capillary responses. A two-photon image of a vascular cast from the same animal is shown alongside, with the corresponding arterioles and veins highlighted. The temporal signatures of each different vascular compartment were crossvalidated with in-vivo two-photon imaging and provided new insights into some of the basic assumptions of the hemodynamic models for fMRI interpretation described before.19 Since LOT detects scattered light, it is highly sensitive to blood oxygenation and HbT changes and is also suitable for fluorescence imaging. As a result, LOT in now being developed to allow simultaneous depth-resolved imaging of hemodynamics and VSDs or CaSDs. Depth-resolved imaging of VSDs and CaSDs would allow the laminar properties of neuronal responses to be studied in detail. Currently, such measurements can only be made using highly invasive depthresolved electrode arrays, although demonstration of a gradient index lens-based system also recently showed promising results imaging VSDs in mouse cortex to depths of 150 m.114 1990,115 and has since become an invaluable tool for imaging intact biological specimens at very high resolution, to depths of up to 600 m. Two-photon microscopy has more recently been applied to in-vivo imaging of the brain.3,19,79,116,117 While fluorescent contrast is required, two-photon microscopy provides an unprecedented view of in-vivo brain activity on a cellular and microvascular level. 2.2 Two-Photon Microscopy Two-photon microscopy is a distinctly different functional imaging technique compared to the optical imaging described so far. Two-photon microscopy was first demonstrated in Journal of Biomedical Optics 2.2.1 Two-photon microscopy: basic principles Two-photon microscopy is similar to confocal microscopy, in that it requires a focused beam of laser light to be steered within tissue, sensing the properties of each location and using them to form a 2-D or 3-D image. Confocal microscopy relies on the rejection of scattered light by isolating signals originating from the focus of the scanning beam. Confocal xcitation wavelengths are typically in the visible spectrum, where tissue scatter and absorption are high see Fig. 1 . Therefore, when using confocal microscopy to image to depths beyond 200 to 300 m, the laser beam can no longer focus, and images become blurred and lack sensitivity. Also, while it is assumed that the detected fluorescent light originates from the focus of the scanning beam, some of the detected signal is also generated by excitation and emission light that is scattered or absorbed within the tissue above and below the focus. Not only does this light contribute to image blurring, it can also result in cumulative photodamage to areas of the tissue that are not actively being imaged. Two-photon microscopy overcomes many of the disadvantages of confocal microscopy. Instead of exciting fluorophores at visible excitation wavelengths, two-photon microscopy utilizes laser light of twice the excitation wavelength. When sufficient photon flux is achieved, it is usually possible for a fluorophore to be excited by two of these lower energy photons arriving in quick succession. Once excited, the fluorophore will emit a photon at or near its usual emission wavelength. This is illustrated schematically in Fig. 4 a . This high photon flux is achievable if a pulsed laser is used, which delivers very high energy, rapid pulses typically a Ti:sapphire laser, tunable between 700 and 1000 nm, with 80-MHz pulse repetition, 500-fs pulse width, mean power 1 W, peak power 25 kW . Even with such high peak power, the nonlinear effect of two-photon excitation will only happen at the tight focus of the laser beam. As a result, two-photon microscopy can use NIR light to image the same fluorophores as confocal microscopy. However, since scattering and absorption of NIR light in tissue is much lower than for visible light see Fig. 1 b , the laser beam’s focus can be maintained to depths of 600 m, depending on the tissue type. The fact that fluorophore excitation only occurs at the very focus of the beam means that surrounding tissue does not experience significant photodamage. The NIR excitation light is generally spectrally well separated from the emission wavelength, making rejection of excitation light much simpler than for confocal microscopy. Further, since it is expected that any light emerging from the tissue at the fluorescent emission wavelength at a particular time can only have originated from the focus of the scanning beam, it is not necessary to attempt to reject scattered light, as required in confocal microscopy. Instead, sensitive detectors can be placed as close to the tissue as possible to collect all of the emerging emission light. This September/October 2007 051402-11 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 4 Two-photon microscopy of in-vivo brain function. a Basic mechanism of two-photon fluorescence. b Schematic of surgical preparation of exposed cortex, with sealed glass window and microscope objective positioning. Green dot shows location of two-photon fluorescence. c Examples of two-photon maps of the vasculature following intravenous injection of dextran-conjugated fluorescein. Black dots and stripes show red blood cell motion. d Dual-channel imaging of neuronal green and vascular red signals: left Oregon Green 488 BAPTA-1 AM calciumsensitive dye stained neurons and right transgenic mouse expressing green fluorescent protein GFP in a subpopulation of neurons mouse supplied by Friedman, Rockefeller University, New York130 . Texas dextran red is the intravascular tracer in both cases studies by Hillman, Bouchard, Ruvinskaya and Boas . e Three-channel imaging of Tg2576 APP Alzheimer’s disease mouse model with amyloid- targeting dye blue , GFP expressing neurons and dendrites green , and vasculature red . Reproduced with permission from Spires et al.14 f 2-D camera-based imaging of hemodynamics and calcium response followed by functional two-photon imaging of the calcium response on the scale of single neurons in the same rat. Oregon Green 488 BAPTA-1 AM was loaded with serial pressure injections into the cortex. Stimulus was 1 mA, 3- s pulses at 3 Hz for 4 sec studies by Hillman, Bouchard, Ruvinskaya and Boas . Journal of Biomedical Optics 051402-12 September/October 2007 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Table 1 Summary of optimal two-photon microscope characteristics for in-vivo brain imaging. Custom-built in-vivo two-photon microscopy systems Almost always upright. Commercial two-photon microscopy systems Often inverted. Feature Form factor Imaging platform Require Z translation of stage, poor clearance between stage Fixed Z stage objective moves up and down, not stage . Motorized x − y stage. Very open clearance around objective and objective. Motorized x − y stage movement very expensive. for animal and monitoring equipment. Two-photon detector may be in place of confocal detector, therefore descanned and through pinhole—very poor light efficiency. Otherwise retrofitted detectors are rarely close to the objective. Similar objectives can be chosen for commercial systems if compatible. Detection geometry Detectors positioned very close to back of objective for maximum light collection. Multiple spectrally resolved detectors may be stacked to allow two or more emission wavelengths to be simultaneously imaged. Numerical aperture Very high numerical aperture but long working distance allows good z resolution, high efficiency, and ability to image through exposed-cortex preparation see Fig. 4 b . Scan rate Software Multimodality Typically designed to allow very high frame or line-scan rates Frame rate well optimized, but limited adaptability for to image dynamics of functional activity. patterns and line scans. Custom software allows wide variety of scan patterns and routines for a wide range of experiments. Typically possible to add additional imaging and measurement modalities to the system owing to accessible light path and animal. Can be integrated into custom software. Often proprietary and difficult to modify. Fixed scan patterns and routines. Very difficult to modify or interrupt light path. Journal of Biomedical Optics impossible. As a result, many researchers have instead built their own two-photon microscopy systems, optimized for invivo imaging.122–126 Table 1 summarizes the key attributes of in-vivo two-photon microscopy systems that are specifically designed for in-vivo imaging. 2.2.2 Two-photon microscopy: examples As for 2-D optical imaging, in-vivo brain imaging research using two-photon microscopy has also found applications in areas of functional processing, functional mechanisms, and disease research. These applications utilize a variety of methods to introduce fluorescent contrast to the brain, including intravenous injection of dextran-conjugated dyes to reveal blood vessels, topical application, or pressure injection of dyes into the cortex, systemic delivery of dyes that cross the blood brain barrier, and transgenic mutation of cells to express fluorescent proteins. Examples are given below. Functional processing. A recent example of in-vivo twophoton microscopy of brain function was the use of calciumsensitive dyes to look at visual responses in cats and rats.3,127 Data were acquired using a custom-built two-photon microscope, and used Oregon Green 488 BAPTA-1 AM dye Invitrogen excited at 800 nm to provide calcium-sensitive fluorescent contrast. The dye was pressure injected into the cortex as described before for 2-D imaging of CaSDs. The complex cortical patterns of the response to visual stimuli of different orientations similar to the ocular dominance maps in Fig. 3 b were found to occur even on the scale of discrete individual neurons. September/October 2007 051402-13 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Functional mechanisms. The advent of functional twophoton microscopy has allowed the mechanisms of neurovascular coupling to be studied in vivo in unprecedented detail. Whereas 2-D optical imaging of the cortex, as described earlier, has enabled parametric comparisons between electrophysiology and the hemodynamic response,18 two-photon imaging enables the interrelations between individual neurons, glial cells, and vessels to be examined in intricate detail. Both confocal microscopy128 and two-photon microscopy79 have been used to examine the dynamics of red blood cells flowing through capillaries in the exposed cortex. In Ref. 79 and more recently in Refs. 116 and 129, the capillary response to stimulation and hypercapnia have been examined using fluorescent plasma tracers such as dextran conjugated fluorescein or Texas dextran red Invitrogen . These high molecular weight dyes cannot pass the blood brain barrier or stain red blood cells, and so provide strong contrast to vessels. Typical in-vivo image stacks of fluorescein-perfused vessels in rat brain are shown in Fig. 4 c . The dark stripes and dots are moving red blood cells whose dynamics can be quantified from their motion in sequential images or line scans.126 The flow and dilation dynamics of pial arterioles and veins was also recently investigated using two-photon microscopy.19 Figure 4 d left shows simultaneously acquired in-vivo rat brain measurements of neurons labeled with Oregon Green 488 BAPTA-1 AM dye and vessels containing Texas dextran red. Both color channels were acquired with two-photon excitation at 800 nm and 550-long-pass dichroic and emission filters to separate the two emission bands onto different photomultiplier detectors. Figure 4 d right shows similar data acquired in a transgenic mouse, engineered such that only certain neurons in the brain express green fluorescent protein.130,131 Capillaries perfused with Texas dextran red can be seen to be weaving amidst the webs of dendrites. 2-D optical imaging of hemodynamic and CaSD changes can also be acquired prior to imaging the same animal with two-photon microscopy. Figure 4 f illustrates 2-D imaging of the calcium-sensitive Oregon Green 488 BAPTA-1 AM response in the somatosensory cortex of a rat during forepaw stimulus 1.5 mA, 3- s pulses at 3 Hz for 4 sec , followed by two-photon imaging of individual neurons at a depth of 315 nm in the cortex responding to the same forepaw stimulus. The two-photon fluorescence of the CaSDs increased by up to 17% during stimulation.126 Such direct and simultaneous observations of both neuronal and hemodynamic responses and structures in in-vivo systems enable the intricate mechanisms underlying neurovascular coupling relations to be closely examined and compared to macroscopic-scale measurements and even fMRI results. Two recent studies used similar two-photon imaging techniques to present compelling evidence supporting the role of astrocytes as active mediators of blood flow control.80,132 Astrocytes make up a large fraction of the cells in the brain, and have been shown to be both metabolically active and capable of transmitting signals between neurons and blood vessels.133,134 In Ref. 132, astrocytes in mouse brain were labeled with calcium-sensitive rhod-2 AM, caged Ca2+, and DMNP-EDTA AM, while the vasculature was perfused with dextran conjugated fluorescein. Two-photon imaging of astrocytic endfeet contacting arteries was performed in vivo at 825 nm, with spectrally separated detectors for rhod-2 AM Journal of Biomedical Optics and fluorescein. A second laser 355 nm was then focused onto an endfoot and pulsed to cause photolysis of the DMNPEDTA, triggering intracellular uncaging of calcium in the endfeet. This forced increase in calcium caused the endfoot to dilate. In a subsequent study, a CaSD fluo-4 AM and astrocyte-specific dye SR101 were used to demonstrate that similar calcium increases in astrocytic endfeet could be observed during prolonged whisker stimulation in vivo.80 Functional imaging to investigate pathologies and treatments. In-vivo two-photon imaging of intravascular dextranconjugated fluorescein in rats was recently used to investigate the effects of stroke on blood flow dynamics, as well as the recovery of flow following two treatments: hemodilution and injection of a recombinant tissue-type plasminogen activator.135 The study also created the strokes using a second laser delivering targeted high fluence ultra-short NIR pulses to cause vessel rupture, intravascular clotting, or extravasation of blood components depending on the laser power. Alzheimer’s disease research has been greatly impacted by the availability of fluorescent markers that bind to beta amyloid plaques in the brain.5,136,137 The so-called Pittsburg compound B 2- 4’-methylaminophenyl -6-hydroxybenzothiazole, PIB can be used as a clinical contrast agent for PET imaging and is showing promise as a diagnostic tool in humans.138 However, since PIB is fluorescent, it can also be used as a two-photon microscopy contrast agent for exposed cortex imaging.117 Tg2576 APP transgenic mice develop beta amyloid plaques over the course of their lives, and are assumed to provide a realistic model of the progression of Alzheimer’s disease. Since Alzheimer’s in humans can currently only be confirmed postmortem, these mouse models offer an opportunity to investigate the early effects of the disease, as well as learn about its progression and to allow development of earlier treatment and diagnostic methodologies. Mice with chronically implanted cranial windows can be repeatedly imaged with two-photon microscopy for up to a year. This has allowed two-photon microscopy to be used to evaluate the effect of developing amyloid- plaques on the structure and function of surrounding neurons and blood vessels. Figure 4 e shows three-channel two-photon image stacks of the brain of a 21 to 24 month old transgenic mouse with beta amyloid deposits shown in blue, blood vessels perfused with Texas dextran red shown in red, and neurons expressing green fluorescent protein in green. The amyloid binding dye methoxy-XO4, similar to PIB was injected intraperitoneally the day before imaging. The neurons were labeled by transfection from an adeno-associated virus injected into the brain 2 to 3 weeks prior to imaging.14 Extensions of this work to noninvasive animal and perhaps even human imaging are being explored.139 2.2.3 Methodological extras Fluorescence lifetime imaging and microscopy FLIM has also been applied to neuroimaging, allowing potential enhancement of intrinsic contrast via environment-dependent lifetime shifts. FLIM has also been utilized to improve measurements of fluorescence resonance energy transfer FRET , September/October 2007 051402-14 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… which can provide measures of the proximity of two fluorescently labeled entities and therefore direct visualization of biochemical interactions.140 There is an extensive body of neuroscience literature that utilizes living brain slices to study neuronal activity in intricate detail. Mice and rats are rapidly sacrificed before their brains are removed, sliced, and placed in a carefully prepared perfusate that supplies nutrients and oxygen. Important results have been obtained by imaging brain slices with two-photon microscopy. For example, using brain slices, Kasischke et al. demonstrated that it was possible to use two-photon microscopy to measure the intrinsic fluorescence of NADH, and distinguish between cytoplasmic and mitochondrial NADH dynamics.141 It has also recently been shown that voltagesensitive responses in neurons and dendrites in brain slices can be measured via second harmonic generation SHG microscopy when neurons are loaded with a SHG contrast enhancing agent called FM4-64.142,143 Noninvasive optical brain imaging on animals is generally performed as a precursor to clinical imaging in babies and adults. It may also be used to improve understanding of clinical data, including studies of neurovascular coupling, stroke, birth asphyxia, etc.166–169 Noninvasive animal optical brain imaging has also been used for longitudinal development studies,170 and measurements where exposing the cortex is not feasible, or would detrimentally affect the experiment. There are two main technological approaches to noninvasive optical brain imaging: fairly simple topographic measurements of cortical function as shown in Figs. 5 a and 5 b , and 3-D tomographic imaging of the brains of infants as shown in Figs. 5 c and 5 d . The basic principles of these approaches are described next, followed by specific examples of their in-vivo neuroimaging applications. 3.1 Topography 3.1.1 Experimental method Noninvasive topographic optical measurements of the human head can interrogate hemodynamic changes in the cortical surface of the brain Fig. 1 b . The majority of functional processing occurs in the cortex, and so optical topography is commonly, and extensively, used to measure functional responses to stimulus in children and adults. Early classifications of the roles of the different areas of the cortex were deduced from observations of the effects of injuries to different parts of the brain.171,172 More recently, these cortical areas have been refined and investigated using functional brain imaging modalities, including fMRI and animal imaging studies such as those described in Sec. 2. Figure 5 e shows the basic regions of the cortex and their functional roles. Optical imaging studies of brain function often focus on the occipital or visual cortex at the back of the head looking at responses to visual stimulus ,147–150 and the motor and somatosensory cortices between the crown and the ear in response to movements, tactile stimulation, or pain .151–154 However, many other studies have been reported, including those on the parietal,155 frontal,156 temporal,157 and olfactory158 lobes. In most functional imaging experiments, a well-defined stimulus is devised such as tapping fingers, or watching a pattern or picture . A subject then performs or receives the stimulus at regular or carefully controlled intervals. In almost all cases, data are then processed to isolate signals that correlate with the presentation of the stimulus. In all functional imaging modalities, and especially for higher mammals such as humans, there is always some variability in baseline brain activity during an experiment that can generally only be eliminated by performing repetitive trials and averaging or filtering. fMRI is rapidly becoming a clinical tool with which to investigate the effects of diseases and drugs, to plan surgeries, and to study brain development.9 Optical topography is also finding similar clinical applications.147,163–165 3.1.2 Basic principles: near-infrared spectroscopy for topography Direct topography. Direct topography is the most widespread approach to clinical optical brain imaging and is often referred to as near-infrared spectroscopy NIRS . In its simSeptember/October 2007 Noninvasive Optical Brain Imaging For routine optical brain imaging in humans, it is clearly necessary to develop noninvasive techniques. Approaches must overcome the effects of light scattering, while maintaining the strong benefits of functional optical contrast. To date, noninvasive optical brain imaging has been applied to babies, children, and adults for studies of functional activation, and for investigation of pathologies. Since light penetration and scattering are significant obstacles for noninvasive imaging, measurements on children and babies are easier to achieve. It is often difficult to image these infants and small children using fMRI and PET, and so optical imaging provides a unique opportunity to study many aspects of functional brain development that have previously been inaccessible.144 Imaging of pathologies in infants is also highly motivated by the prevalence of perinatal brain injuries due to prematurity, birth asphyxia, and congenital conditions.34 Thanks to the portable and nonionizing nature of optical imaging, it is possible to obtain optical measurements on babies within a neonatal intensive care unit,145 and on young children who are able to sit on a parent’s lap in a familiar environment.146 Optical imaging of the adult brain can also provide valuable functional information that complements modalities such as fMRI. A particular advantage is optical imaging’s ability to simultaneously image HbO2, HbR, and HbT. Optical imaging is being widely adopted for studies of the cortical hemodynamic response to a wide range of stimuli.147–158 Multimodality imaging of the adult head during fMRI,150,159,160 electroencephalography EEG 161 and magentoencephalography MEG 162 has provided insights into the contrast mechanisms of the fMRI BOLD signal, and allowed oxygen consumption and neurovascular coupling relations to be investigated in humans. Optical imaging in adults is also being used as a tool for monitoring and detecting hematomas and stroke, as well as the effects of other pathologies and surgery on brain perfusion.147,163–165 Noninvasive optical imaging holds the potential to provide accessible monitoring of brain function in situations where other imaging modalities such as fMRI are prohibitive, e.g., due to imaging environment or cost. Journal of Biomedical Optics 051402-15 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 6 Time-domain, frequency-domain, and continuous wave measurement sensitivities to absorption and scatter. Schematic shows simplistically how absorbing and scattering structures will change TD and FD data, but that in transmission, their effect on cw data will be indistinguishable. Fig. 5 Configurations for noninvasive optical brain imaging. a and b Functional topography uses one or more source-detector pairs to create a pseudo-2-D map of the underlying cortex. Sources and detectors may be optical fibers, or light emitting diodes and photodetectors placed directly on the skin. c and d Tomography places sources and detectors over the volume of the head, and detects light that has traveled through many projections. This approach allows subsequent reconstruction of the optical properties of the entire brain. e The main functional areas of the human cortex. f Simulation of the sensitivity of a noninvasive topographic measurement between a source and detector. The simulation used Monte-Carlo modeling and accounted for the heterogeneous 3-D structures of the adult scalp and cortex. Reproduced with permission from Boas et al.202 plest form, NIRS uses light source and detector pairs positioned on the head of infants or adults. Light that passes through the skull scatters within the cortex of the brain as illustrated in Figs. 1 b , 5 b , and 5 f . Despite significant attenuation, variations in absorption due to changes in local HbT, HbO2, and HbR concentration result in detectable modulations in the intensity of the emerging light. If an array of sources and detectors is used, the changes measured between each source-detector pair can be approximately mapped or interpolated to the underlying cortex, and a low-resolution 5 mm 2-D topographical map of cortical activation can be produced.152,154,173 When many overlapping and multidistance topographic measurements are acquired, it is possible to use mathematical models of light propagation to create a depthresolved reconstruction of functional activation.174–176 While more complex, this approach can overcome several of the limitations of direct topography, as described below. Instrumentation. Noninvasive optical brain imaging instrumentation may be continuous wave cw , frequency domain FD , or time domain TD . cw is the most common approach and simply measures the intensity of light emerging from the head at two or more wavelengths. Data acquisition is often frequency encoded up to kHz , such that multiple wavelengths177,178 or multiple sources149,179–181 can be measured simultaneously. FD and TD instruments are substantially more complex and expensive, and can measure both a cw component and the time taken for light to pass through the tissue.182 FD instruments achieve this by modulating their illumination at around 100 MHz and measuring the phase delay of the emerging light.152,183 TD instruments use a pulsed laser and measure the temporal distribution of the photons emerging from the tissue.145,184–188 The main purpose of this temporal information is to provide a measure of the scattering properties of the head, as illustrated in Fig. 6. Optical fibers are generally used to transmit light to and from the head see Fig. 5 . Direct placement of small sources and detectors on the head has also been demonstrated.166,189 Calculation of HbO2 and HbR concentrations. The change in intensity measured between a topographic source and detector the attenuation will generally reflect a change in the September/October 2007 Journal of Biomedical Optics 051402-16 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… absorption of underlying tissue Eq. 1 . The basic approach for deriving changes in HbO2 and HbR concentration from changes in attenuation was outlined in Sec. 1.1.2 Eqs. 1 and 5 . Determining I0 and G for noninvasive measurements is particularly challenging, since the amount of light coupled into and out of the scalp for every source and detector position would need to be known. Differential measurements relative to time t = t0 Eq. 3 should cancel these source- and detector-specific scaling factors, although it must be assumed that these factors stay constant for the duration of the measurement which can be challenging if the subject moves and/or hair is obstructing probe placement . As with exposed cortex imaging, if measurements of the change in attenuation at multiple NIR wavelengths are to be converted into changes in hemoglobin concentrations, the effects of scattering must be accounted for. If Eqs. 4 and 5 are to be used for “direct” topography, an estimate of DPFx for each wavelength and measurement pair is required. FD and TD systems can be used to measure DPFx directly, and as an adjunct to cw imaging can provide DPFx for a given wavelength by measuring the average time taken for light to pass from the source to detector, and converting to distance based on the speed of light in tissue.190,191 However, if a subset of DPFx measurements or generalized values of DPF are to be used, it should be noted that DPFx could vary depending on age, gender, and even position on the head, owing to the complex structures of the underlying scalp, skull, and convoluted human cortex. Simulations for calculation of DPFx can become increasingly complex if efforts are made to incorporate the curvature of the head, the heterogeneous structures of the scalp, skull, and cortex, the presence of cortical convolutions and sulci, the effects of CSF see Fig. 1 b , and even the effects of isotropic scattering in the brain.88,192–199 Estimates of the absorbing and scattering optical properties of each tissue type,200,201 as well as accurate knowledge of the source-detector separations between fibers on the head, are also required. If available, an MRI scan of the head can be used to build a customized model of light propagation for each subject, allowing constrained analysis that accommodates the exact structure of the brain, skull, and scalp underlying each detector.202 However, while this is feasible in studies where fMRI is simultaneously acquired, it is prohibitive for other applications of optical imaging. Although TD and FD systems rarely have a large number of multiwavelength channels,161,185 they have the advantage that they can be used to deduce the scattering properties of the head, which can then be incorporated into more complex calculations of hemoglobin concentrations.201,203 An FD approach is to measure between several source-detector separations and, under simplifying assumptions of homogeneity, fit for scattering and absorption based on a linearized model of light propagation versus distance.183 Such approaches do not require estimates of DPFx. In general, for direct topography, simplified models and scaling factors are applied to account for DPFx and also partial volume effects, which correct for the fact that only a fraction of the detected light actually sampled the cortex.34,153,159,183,204 Several studies have explored the importance of accurate calculation/measurement of DPFx and its Journal of Biomedical Optics impact on the determination of hemoglobin concentrations for optical topography.198,205,206 It has been shown that in many cases, useful functional information can be obtained even when simplifying assumptions are applied to analysis and interpretation. However, it is therefore important to consider the possible impact of systematic errors resulting from this analysis for all NIRS measurements, particularly when imaging the adult head, where many confounding effects can obscure cortical signals and impact correct interpretation. Extracerebral effects and physiological noise. The effects of blood flow changes in the scalp, and changes in the brain that are not directly related to the cortical stimulus response, must also be considered. The analysis approaches described above assumed that changes in signal were solely due to changes in the cortex of the brain. However, it has also been shown that breathing 0.16 Hz , heart rate 1 Hz , and 0.1 Hz 207–210 can the “vasomotor” or Mayer wave signal cause large changes in measured hemodynamic signals. Substantial repetition of a stimulus and averaging should reduce the effects of unrelated physiological variations in blood flow in the scalp and brain. However, acquisition of so many repetitions is not always possible with adult or infant subjects. In addition, some responses may not average out, or may intrinsically correlate or become entrained with the stimulus period,208 and can affect both brain blood flow and the rich blood supply to the scalp. Recent evidence has suggested that certain stimuli actually invoke a systemic response possibly through heart rate or blood pressure changes that could manifest as a change in scalp or even baseline cortical blood flow that correlates with the stimulus period.153,211 Several measurement, data analysis, and filtering approaches have been proposed to reduce the influence of hemodynamic activity unrelated to the true cortical response. One approach is to simultaneously acquire measurements between close and more distant pairs of sources and detectors, assuming that the closer detector is predominantly detecting signal from the scalp. This scalp signal can then be approximately removed from the signal from the more distant detector, assuming that a similar superficial contribution is affecting both measurements. A similar approach can be achieved with TD data, where the photons arriving at a detector earlier are likely to have traveled more superficially, and therefore contain information about the scalp signal. The benefit of using TD is that later photons that have traveled more deeply between the same source detector travel through the same region of the scalp as the early photons, so cancellation should be improved relative to the multidistance approach.187,212 An extension of this is reconstructed topography, which creates a depth-resolved map that should be able to spatially separate superficial signals from deeper cortex signals, if suitable measurements are acquired. Spatiotemporal filtering and modeling have also been proposed as a way to remove these systemic signals.213,214 While rapidly becoming more sophisticated, implementation of these data processing techniques can be complex and sometimes subjective. Extensive processing could result in distortion of true cortical signals, particularly when attempting to separate scalp and cortical signals that may both have elements correlating with the stimulus response.211 September/October 2007 051402-17 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Practical measurements. A final difficulty with acquiring functional optical imaging data is reliable placement of sources and detectors on the head. The primary problem is to be able to securely position multiple probes in a pattern such that their interfiber distances are known, and they stay stable during the course of the measurement. Hair on the head is a major challenge, as it can significantly attenuate light entering and leaving the head, and can cause measurement instability over time if the probes shift. Technical advancements have been made in positioning optical fibers and probes onto the head, and relatively stable configurations have been devised.180 However, further work will be required before it becomes possible to routinely, rapidly, and repeatably position optical topography fibers on the head. Reconstructed topography. Direct topography is the simplest way to analyze topographic data, but relies on many assumptions and simplifications. An alternative approach is to reconstruct data onto a 3-D rendering of the head. Although analysis becomes more complex, by creating a depth-resolved image of the measured absorption changes, some of the confounds described earlier can be overcome. Image reconstruction. Instead of requiring a pathlength estimate, reconstruction effectively attempts to recreate the measured data by simulating alterations in the absorption distribution within the head.187,215,216 When the simulated data and the measured data match, the absorption distribution is assumed to be correct. Reconstructing still generally requires an estimate of the baseline absorbing and scattering optical properties of the head, and an accurate model of light propagation, but can help if there are depth-dependent variations in the response, including effects from the scalp. Partial volume effects, particularly where one wavelength of light effectively probes a different region, than a second wavelength, can also be accommodated in reconstructed topography. For reconstructed topography to be effective, it is necessary to measure signals between overlapping sources and detectors, and generally at several source-detector distances see Fig. 5 b . These measurements provide a set of mutually consistent data for which there should be only one possible corresponding differential absorption distribution. This in turn should also lead to improved x-y resolution.175,176,179 Modeling of light propagation in the human head typically uses either Monte Carlo modeling or the diffusion approximation DA to the radiative transfer equation Eq. 6 and a linearizing assumption such as the Born Eq. 7 or Rytov approximations. − · = q0 , M n,m = Jn,m r r , where is the photon density within the tissue, and q0 is the source function. M n,m is the set of changes in measurement from time t = 0 between source n and detector m, and Jn,m r is the Jacobian matrix derived from the Greens function solution to Eq. 6 , representing the change in measurement M n,m corresponding to a change in absorption a at position r within the tissue.215 An example of Jn,m r for a single Journal of Biomedical Optics September/October 2007 051402-18 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… rates.152,181 However, the attenuation difference between a source and detector separated by 2.5 and 4.25 cm can be up to two orders of magnitude.179 For reconstructed topography, these multidistance measurements are vital, and detectors must have excellent dynamic range and noise characteristics to be able to simultaneously distinguish between frequencyencoded signals of such varying amplitudes. Development of faster systems, with better dynamic range or highly advanced and rapid multiplexing schemes is therefore ongoing.176,179 3.2 Tomography The most complex approach to optical brain imaging involves performing a full 3-D reconstruction of the entire head using an approach similar to x-ray computed tomography. This approach has only been demonstrated on small infants, since it requires light to be transmitted across the whole head as shown in Figs. 5 d .145,221–223 However, the benefits are that in addition to functional imaging, this approach can potentially enable imaging of pathologies such as periventricular hemorrhage and ischemia in a perinatal setting, where currently no alternative diagnostic and monitoring tools are available.217 3.2.1 Basic principles Since the use of 3-D tomography for brain imaging is not yet widespread, specific details of one study are given in the next section, and only a brief overview is provided here. For tomographic imaging, it is necessary to acquire a series of measurements or projections through the tissue being imaged. For reconstructed topography, measurements between sources and detectors with different separations allow some depth-resolved information to be deduced. In tomography, the same basic image reconstruction approach is used, but a larger number of projections through the head at different angles are required to calculate the 3-D structure deep within the brain see Figs. 5 c and 5 d . Note that optical tomography has also been used to image the adult breast for the detection or monitoring of breast cancer224 and the human forearm.225 There are three main difficulties with acquiring tomographic data on the head. The first is that light must be detected through greater than 7 cm of tissue, which in general will cause around 10−7 attenuation of the incident light. The second difficulty is that the paths of light through the head are again a function of baseline anatomical structures. Third, to look at pathologies or regions of ischemia for which a baseline e.g., prestimulus state is not available, it can be valuable to attempt to reconstruct this absolute absorbing and scattering structure of the brain, and not just measure changes as with topography . 3-D tomographic imaging of the neonatal had has been demonstrated using TD instruments.145,222,223 TD measurements that use single photon counting have the advantage that very low light levels can be detected by increasing measurement integration time. Background signal and thermal noise will not be temporally correlated, and can often be subtracted from the true signal.226 Therefore, while slow, TD systems can detect usable signals across 15 cm of adult breast, or an entire neonatal infant’s head. TD and FD measurements also have the benefit that they allow the effects of scatterers and absorbJournal of Biomedical Optics ers to be distinguished as illustrated in Fig. 6,215,227 although see Ref. 228 . In fact, using TD data, it is possible to reconstruct images of both absolute baseline absorbing and scattering structures.229 These methods are currently being developed to allow imaging of the baseline optical properties of the infant brain. Differential measurements can be used to reconstruct absorption changes, under the same assumptions as topography: that scattering does not change and the effects of heterogeneous background absorbing and scattering structure on light paths is insignificant or can be incorporated into a model of the head.217,219 Differential imaging of 3-D absorption changes can be achieved with relatively simple linearized reconstruction approaches similar to those used for reconstructed topography described before.217,222 Values of Jn.m r must be calculated based on the 3-D shape of the head, and the precise positions of the sources and detectors see Fig. 5 c . Typically, this is done using finite element models implementing the diffusion approximation, and sometimes incorporating estimates of the baseline absorbing and scattering structures within the head.230,231 For absolute imaging of the baseline absorbing and scattering structure of the brain, a nonlinear algorithm should be used that iteratively seeks a 3-D distribution of absorbing and scattering that can account for the data measured on the head. A highly complex and publicly available software package developed to do both nonlinear reconstruction of absorption and scattering images, as well as forward modeling of Jn.m r using the finite element method is called TOAST timeresolved optical absorption and scattering tomography .215,232,233 The data that are actually fed into these reconstruction algorithms will have a significant effect on image sensitivity to scatter and absorption.217 FD data consists of measurements of amplitude and phase. Measurements are typically made at modulation frequencies between 100 and 200 MHz.234 TD data represents the full impulse response of the tissue, and therefore amplitude and phase at almost any effective modulation frequency can be calculated. Datatypes such as the mean-time or integrated intensity of the signal can also be extracted.226 For images of both absorption and scatter, the image reconstruction algorithm must utilize both amplitude and phase data, or at least two sufficiently different timeresolved parameters. A very accurate model of light propagation and precise calibration of data for systematic errors is also very important, particularly if absolute images of absorption and scatter are required since the data must closely match simulated data corresponding to the correct 3-D distribution of absorption and scatter .217,226 3.3 Examples of Noninvasive Optical Brain Imaging Figure 7 shows a variety of noninvasive optical brain imaging results showing human functional responses to stimulus. Figure 7 a shows results from a study using cw direct topography to examine the hemodynamic response to different kinds of stimuli: finger tapping, median nerve electrical stimulation, and tactile stimulation.153 Stimuli were presented for 10 sec, followed by 18 sec of rest, and were repeated 20 times per subject. The cw system has 16 laser diode sources, eight at 690 nm and eight at 830 nm, frequency encoded between 4 September/October 2007 051402-19 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Fig. 7 Examples of noninvasive optical imaging of functional responses in humans. a Direct functional topography images and time courses acquired for two different motor/somatosensory stimulus paradigms on one subject. A localized HbO2 and HbR response is seen in the opposite contralateral side of the brain to the hand being stimulated. Reproduced with permission from Franceschini et al.153 b Functional topography measurements on adult head during fMRI. Both the spatial and temporal characteristics of the responses can be compared between two modalities. Reproduced with permission from Huppert et al.159 c Black and white visual stimuli were presented to the left and right fields of vision. The hemodynamic response was measured using a dense topographic array of sources and detectors. Data were reconstructed using a linear perturbation model. Reproduced with permission from Culver et al.176 d Time-resolved 3-D tomography data were acquired on a 33-week gestation infant undergoing passive motor stimulus to the left arm. Image series shows successive slices through 3-D head from left ear to right ear. Reproduced with permission from Gibson et al.222 and 7.4 kHz such that all were illuminated while 16 avalanche photodiode detectors collected emerging light. Sources and detectors were configured in two panels, arranged as shown in the photograph in Fig. 7 b where four fibers span the center each delivering both 690 and 830 nm—red arrow , straddled by two rows of four detector fibers—green arrows. The detected signal was digitized at 40 kHz and analyzed to extract signals coming only from the closest sources 3 cm to each detector. Data were then low-pass filtered and downsampled to 10 Hz, and then bandpass filtered between 0.02 and 0.5 Hz to eliminate heart beat pulsations and slow drifts. Data were then block averaged over the 20 trials to reveal the average attenuation changes corresponding to each stimulus. These attenuation values were converted to absorption changes using literature values of the differential pathlength factor DPF at 690 nm 6.51 and 830 nm 5.86 , such that the pathlength is DPFx, where x is the source-detector separation of 3 cm. The resulting absorption changes were then converted to HbO2 and HbR changes using Eq. 4 . These source-detector-specific changes in hemodynamics were then projected onto a map underneath each source-detector pair using an approximate weighted spatial distribution based on the volume sampled. The final map was then interpolated onto a finer grid. The images and time courses in Fig. 7 a show the mean maps and time courses for a single subject. Contralateral responses are those seen in the side of the brain Journal of Biomedical Optics opposite the hand receiving the stimulus where responses would be expected . Figure 7 b shows a similar study, but where optical topography data were acquired during simultaneous fMRI at 3 tesla.159 The subjects performed a finger tapping exercise for 2 sec, with a randomized interstimulus interval between 4 and 20 sec repeated between 27 and 32 times per fMRI scan, with each scan repeated four to six times. fMRI scans with a 2-Hz frame rate were used to measure either BOLD contrast or ASL data, which provides blood flow information. The cw NIRS instrument and data processing approach was similar for this and the study described previously.153 The HbO2, HbR, and HbT time courses for one subject and for each source-detector pair are shown in map form inset lower right of Fig. 7 b . Oil-filled fiducials on the optical fiber heads allowed the optical probe’s location to be coregistered with the MRI scans. The time course of the fMRI BOLD signal is compared to the HbO2, HbR, and HbT temporal responses measured using NIRS. The group average of five subjects is shown, and the work concludes that BOLD most closely temporally correlates with HbR changes. ASL changes were found to correlate most closely with HbT changes data not shown . A second paper continues this analysis to a spatial comparison between the localization of HbR and BOLD responses.235 September/October 2007 051402-20 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… Figure 7 c shows a reconstructed topography result, acquired on the adult head during a 10 sec visual stimulus followed by 35 sec of darkness.176 Data were acquired with a multichannel cw NIRS system consisting of 24 radio frequency encoded light-emitting diode LED sources at 760 and 850 nm, and 28 avalanche photodiode detectors. Frequency encoding allowed simultaneous acquisition at both wavelengths, but LEDs were switched on sequentially allowing detectors with 107 dynamic range to measure signals from sources at three different separations: 13.3, 30.0, and 40.0 mm. Despite this serial source switching, full images were acquired at 10 frames per sec. Data were then reconstructed by inverting the linear problem Eq. 7 and using Jn.m r functions calculated using a finite element model simulation package called NIRFAST,236 which incorporates a simplified baseline structure consisting of a 15-mm scalp/ −1 −1 surrounding the skull layer a = 0.17 cm , s = 7.6 cm −1, −1 . brain of radius= 65 mm a = 0.19 cm s = 11.5 cm Both block averaged and even trial-by-trial data show clear hemodynamic responses to unilateral visual stimulus. Figure 7 d shows 3-D optical tomography images of the functional response to passive motor stimulus of the left arm in the brain of a premature infant.222 The image series represents vertical slices through the 3-D volume of the head from left ear to right ear . A clear volume of absorption increase is seen in the region corresponding to the right motor cortex. Data were acquired using a dual wavelength, 32-channel time-resolved optical tomography system known as MONSTIR.186 The images in Fig. 7 d show the change in absorption at 780 nm, but simultaneously acquired images at 815 nm were also created and converted into maps of the change in HbO2, HbR, and HbT.Four sleeping premature infants were measured, born at a median gestation of 28 weeks and 5 days, and imaged at a median gestational age of 34 weeks 5 days. Optical fibers were positioned on the head using a custom-made foam-lined semi-rigid cap, shaped to allow careful but stable positioning of many fibers over the full cranium.217 Since time-resolved data can be slow to acquire, motor stimulus data were collected one source at a time for up to 12 source positions. Each source was illuminated, and the time-resolved transmitted light signal then collected at up to 32 detector positions for 10 sec while the infant rested. The left arm was then gently moved back and forth by an adult for 15 sec, with transmitted light being measured for the last 10 sec of this exercise to capture the resultant hemodynamic response. The next source position was then illuminated, and the 10 sec rest/15 sec exercise cycle repeated. The resulting time-resolved data were Fourier transformed and the log amplitude and phase components extracted at 100 MHz. The difference in the log amplitude and phase measurements between the active and resting states effectively cancels out errors relating to optical coupling and the system’s impulse response function.217,225,226 The resulting differential measurements were then reconstructed into a 3-D map of absorption changes. In this case, a linearized model similar to Eq. 7 was utilized with careful normalization and scaling, and compensation for errors in optical coupling incorporated into Tikhonov inversion.237 The Jn.m r functions were generated using a finite element model implementing the diffusion equation in a homogenous 3-D volume correspondJournal of Biomedical Optics ing to the contours of the infant’s head. The MONSTIR system and finite element algorithm TOAST have been used for a number of other 3-D optical tomography imaging studies on premature infants.145,223,238 3.4 Methodological Extras 3.4.1 Cytochrome oxidase Significant efforts have been made to topographically measure changes in a third intrinsic brain chromophore called cytochrome oxidase. Cytochrome oxidase is a mitochondrial enzyme that plays an important role in oxygen metabolism in the brain. Cytochrome does not change its concentration, but changes its oxidation state and absorption spectrum during the production of adenosine triphosphate ATP . The concentration of cytochrome oxidase in the brain is lower than that of hemoglobin 2 to 3 micromolar compared to 40 to 60 micromolar , so based on published spectra, the peak absorption in brain of reduced cytochrome C would be 0.06 cm−1 at 550 nm, changing to 0.017 cm−1 when oxidized, and 0.00064 to 0.0002 cm−1 at 800 nm compare to equivalent hemoglobin spectra in Fig. 1 .239 It is therefore very desirable, but quite difficult to measure changes in the redox state of cytochrome oxidase in vivo. Nevertheless, measurement methods have been studied carefully240–242 and NIRS-based cytochrome measurements have been reported in adults.209,243 Exogenous dyes for noninvasive imaging. Exogenous contrast agents can also be used for noninvasive brain imaging. Indocyanine green ICG is an FDA-approved dye commonly used for liver and cardiac function tests.244 ICG has strong NIR absorption and fluorescence, and does not cross the blood brain barrier, and so has also found several optical brain imaging applications. A bolus of ICG can be used to quantify blood flow in the brain, by observing the rate of arrival and accumulation of the dye with NIRS measurements.245 Changes in the dynamics of a dye bolus in patients with perfusion abnormalities such as thrombolitic stroke have also been studied using NIRS.246 ICG has been used to simultaneously determine cerebral blood volume and total circulating blood volume in neonates undergoing blood transfusions.247 An interesting recent study reported that it was possible to also observe the fluorescence of ICG in the brain using noninvasive NIRS.248 This result is important as it supports the possibility that more specifically targeted longwavelength dyes could feasibly be measured noninvasively in the brain. Such molecular probes may have sufficient specificity to allow localization or even diagnosis of diseases such as brain cancer and Alzheimer’s.139,249 Noninvasive neuronal imaging via fast scattering. Since it is not possible to use voltage-sensitive dyes in humans, researchers are highly motivated to isolate a strong source of contrast that can directly report changes in neuronal activity. Observation of the fast scattering signal in exposed rodent cortex was described in Sec. 2.1.4. Noninvasive measurement of this signal is highly challenging, as scattering changes are difficult to detect accurately and separate from the effects of absorption. Recent studies of the fast scattering signal have shown compelling results, and compared findings to EEG and fMRI.161,250,251 However, it is still questioned whether the fast signal is sufficiently robust for routine measurement.104,105 Diffuse correlation spectroscopy. A relatively recent addiSeptember/October 2007 051402-21 Vol. 12 5 Hillman: Optical brain imaging in vivo: techniques… tion to optical topography has been the development of diffuse correlation spectroscopy DCS , which allows noninvasive measurement of cortical blood flow. The theoretical basis of this technique is similar to that of speckle-flow imaging described in Sec. 2.1.2, but DCS measures the temporal autocorrelation of laser light with a long coherence length that has scattered through tissue containing moving red blood cells. For human imaging, DCS acquires signals through single mode fibers placed on the scalp.252–254 Promising results in humans have been reported,255 and it has been shown that it is possible to extend the technique to allow 3-D imaging.256 Summary We have reviewed a wide range of approaches to optical brain imaging, and explored their many important applications to neuroscience. We have provided insights into the fundamental basis of many optical imaging techniques, and highlighted important considerations for their practical implementation. We hope that these details will enable those interested in adopting or developing optical imaging methods to understand the key aspects that should be considered when acquiring measurements or analyzing data. We have sampled the huge body of literature that has successfully used optical brain imaging in all its forms for ground-breaking and fundamentally important research and applications. It is clear that optical brain imaging is a rapidly expanding field, continually evolving to embrace new technologies, contrast mechanisms, and in-vivo applications.

Journal

Journal of Biomedical Optics – SPIE

Published: Sep 1, 2007

Keywords: optical imaging; two-photon microscopy; near-infrared spectroscopy; diffuse optical tomography; neuroimaging; neurovascular coupling

There are no references for this article.