30 years ago, Jöbsis demonstrated that it was possible to measure blood and tissue oxygenation changes in the brain of a cat using near-infrared (NIR) light.1 Since then, optical brain imaging has grown into a rich and diverse field, encompassing both basic science and clinical applications that exploit the unique properties of light to image the brain.
There are significant benefits to using light to image living tissue. Most importantly, light can provide exquisite sensitivity to functional changes, either via intrinsic changes in absorption, fluorescence, or scatter, or through the use of extrinsic contrast. Oxy- and deoxyhemoglobin ( and ) are the most widely measured intrinsic chromophores, alongside cytochromes and metabolites, which have distinctive absorption or fluorescence characteristics. Absorbing and fluorescent dyes, as well as cutting-edge transgenic methods, can also provide highly specific optical contrast enhancement, often capable of actively reporting functional parameters such as changes in membrane potential or ion concentrations. 2, 3, 4, 5 A continuing rich supply of new optical contrast agents is also anticipated, particularly since unlike contrast agents for most 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, 7, 8, 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 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-fluoro-2-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. 1a ].
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 imaging of the brain’s surface. Extrinsic dyes and visible wavelengths can also be used [see Fig. 1b].
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, 21, 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.
Macroscopic Optical Imaging of Exposed Cortex in Two Dimensions
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 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 administration 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 better signal to noise (where noise includes normal background physiological variations).
Basic principles: two-dimensional camera-based 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 ( and ) as shown in Fig. 1a. 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 to , so changes in the oxygenation state of blood correspond to changes in the relative concentrations of and . Consequently, variations in the amount and oxygenation state of hemoglobin modulate the absorption properties of the brain. By simply shining light into the exposed cortex and taking pictures with a camera, the hemodynamic response to stimulus can be readliy observed. Further, since and have different absorption spectra [see Figs. 1a and 1b], measurements with different wavelengths of light can produce images that are preferentially sensitive to changes in the concentration of either or . At isobestic points, where the and absorption is the same (e.g., , 530, 570, and ), changes in total hemoglobin concentration 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 .27 By adding a second wavelength (often ), researchers were able to investigate oximetric changes in addition to changes in .28, 29 Since absorption at is low (around 16% of the absorption of ), changes in absorption at are largely attributed to, and interpreted as, changes in . However, it should be noted that even though absorption at is lower than that of , changes can interfere with and distort inference of concentration changes from measurements at .
To calculate the actual changes in concentration of , , and , 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 source, as shown in Fig. 2.30, 31 Full spectra of the reflectance changes along a single line on the cortex have also been used to provide a more detailed spectral picture of the changes occurring during functional activation.32, 33
Calculation of hemoglobin concentrations from such measurements is not trivial, as the effects of scattering must be considered, and then data must be combined with the known spectra for and to deduce the concentration changes corresponding to the measurements. In this case, we assume that scattering in the brain does not change significantly during the hemodynamic response. (Note that the following approach is also applied to noninvasive clinical optical topography, as is discussed in Sec. 2). The basis of this approach is the modified Beer-Lambert law,34 which states that:is the detected intensity, is the incident intensity, is the mean pathlength traveled by light in the tissue (a function of baseline scattering and absorbing properties), is a geometry-dependent factor, and is the absorption coefficient of the tissue as given by Beer’s law: chromophores in the tissue, where is the wavelength-dependent specific absorption coefficient of a particular chromophore (e.g., the absorption spectrum) and 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 apart. Thus, the true pathlength is .
may be spatially varying and difficult to accurately determine. is also very difficult to quantify. Therefore, it is common to consider changes in detected signal relative to an initial time point or a reference state:and ) are acquired, it is possible to use these differential measurements to determine the changes in and concentration occurring relative to time using: and can be calculated using a linear model:
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. 1a], each measurement wavelength will experience different levels of scattering. The mean path of light through tissue 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 in almost all cases. Therefore, as given by Eqs. 4, 5, accurate estimates of these wavelength-dependent pathlength factors are required even to calculate differential changes in hemoglobin concentration ( and ).
Note, however, that if is measured at an isobestic point, and does not change over time, then will be directly proportional to .
Estimation of this wavelength-dependent pathlength 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 can directly impact the correct calculation of , , and changes in tissue.35 Note also that if only values of and are derived, that it is not possible to infer changes in absolute oxygen saturation without first estimating the baseline absolute concentrations of and .
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 , , and .
A typical hemodynamic response in the somatosensory cortex of a rat during a four-second forepaw stimulus ( pulses at ) 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 is generally seen, along with an increase in and a decrease in . 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 , 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, 44, 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-vivo 2, 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 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 after the stimulus begins. The neuronal response is much faster. Measured using VSDs and CaSDs, neuronal responses typically peak within 25 to , and are usually followed by an inhibitory or refractory period of negative signal, which returns to baseline at around 150 to 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 forepaw stimulus ( pulses at ). 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 per frame).
Voltage sensitive dyes typically have quite low contrast (typical ). 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. 1a]. 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 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 decrease in absorption or fluorescence is observed in response to an increase in membrane potential.
CaSDs have much higher contrast than VSDs 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 and absorption). The challenge with 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 to .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 speckle-flow imaging. While hyperspectral and intrinsic imaging provide measurements of the changes in and 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 ). The camera then simply acquires rapid and high-resolution 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., or ), 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.
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 3a 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 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 3b 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 , emission: 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. 3b shows intrinsic signal ocular dominance maps ( 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 images 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. 3c. 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, 68, 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 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 .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 veins.73 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 , , and 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 . These calculations are based on models of the hemodynamic response to stimulus. 32, 74, 75, 76, 77, 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. 3d, where an MRI-compatible fiber optic endoscope was configured to allow synchronized fMRI and functional optical imaging at two wavelengths (570 and ) 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 , 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.
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. 3e, 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
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, 98, 99, 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 and emission is between 500 and . 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).
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 light occurred within of stimulation and correspond well with electrophysiological recordings. A 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 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 ( and ).
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 depth-resolved 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 . 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, 110, 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 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.
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 with 100 to resolution. A recent result is shown in Fig. 3f, where LOT was used to image the exposed-cortex hemodynamic response to rat forepaw stimulus at 473 and . 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 3f shows two cross sectional images through the 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 cross-validated 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 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 depth-resolved electrode arrays, although demonstration of a gradient index lens-based system also recently showed promising results imaging VSDs in mouse cortex to depths of .114
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 1990,115 and has since become an invaluable tool for imaging intact biological specimens at very high resolution, to depths of up to . 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.
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 , 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. 4a . This high photon flux is achievable if a pulsed laser is used, which delivers very high energy, rapid pulses (typically a :sapphire laser, tunable between and , with pulse repetition, pulse width, mean power , peak power ). 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. 1b], the laser beam’s focus can be maintained to depths of , 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 signal, when reformed according to the beam scanning pattern, produces a high resolution image of the fluorescent structures within tissue. Where traditional microscopy required tissues to be stained and sectioned, two-photon microscopy therefore allows high-resolution 3-D imaging of intact, functioning tissue.118
Techniques have been developed to allow in-vivo imaging of the exposed cortex of animals, and the use of two-photon microscopy for neuroscience research is now widespread.119 After exposing the cortex as described in Sec. 2.1.1, a typical approach is to then apply or inject topical dyes and subsequently cover the exposed brain with agarose in artificial cerebrospinal fluid (CSF).120 A small glass cover slip placed over the opening in the skull and sealed in place using acrylic cement can be used to effectively reseal the skull, allowing intracranial pressure to re-equalize, preventing contamination or dissolution of ambient gases and controlling brain movement (which would otherwise be substantially affected by heart rate and breathing). Measurements would typically proceed while the animal continued to be anesthetized, although some researchers have implanted sterile windows or artificial duras that can stay in place for up to a year, allowing longitudinal studies.14, 24, 121
Apart from access to the brain, two difficulties with implementing two-photon brain imaging remain. First, the cost of the laser required for imaging is currently high ($100,000 to $200,000 USD). The second difficulty is that, while commercial two-photon microscopy systems exist, to date most are based on modified confocal microscopy systems, and are also very expensive. This means that some of the benefits of two-photon microscopy are not fully exploited, and modification of the systems to allow live animal imaging is difficult or impossible. As a result, many researchers have instead built their own two-photon microscopy systems, optimized for in-vivo imaging. 122, 123, 124, 125, 126 Table 1 summarizes the key attributes of in-vivo two-photon microscopy systems that are specifically designed for in-vivo imaging.
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.
A recent example of in-vivo two-photon microscopy of brain function was the use of calcium-sensitive 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 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. 3b] were found to occur even on the scale of discrete individual neurons.
The advent of functional two-photon 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, 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. 4c. 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 4d (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 and 550-long-pass dichroic and emission filters to separate the two emission bands onto different photomultiplier detectors. Figure 4d (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 4f 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 ( , pulses at for ), followed by two-photon imaging of individual neurons at a depth of 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 , 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 , with spectrally separated detectors for rhod-2 AM and fluorescein. A second laser was then focused onto an endfoot and pulsed to cause photolysis of the DMNP-EDTA, 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 dextran-conjugated 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-hydroxy-benzothiazole, 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 4e 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
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), 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 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 voltage-sensitive 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
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 , , and . Optical imaging is being widely adopted for studies of the cortical hemodynamic response to a wide range of stimuli. 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 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, 164, 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.
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, 167, 168, 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. 5a and 5b ], and 3-D tomographic imaging of the brains of infants [as shown in Figs. 5c and 5d]. The basic principles of these approaches are described next, followed by specific examples of their in-vivo neuroimaging applications.
Noninvasive topographic optical measurements of the human head can interrogate hemodynamic changes in the cortical surface of the brain [Fig. 1b]. 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 5e 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, 148, 149, 150 and the motor and somatosensory cortices between the crown and the ear (in response to movements, tactile stimulation, or pain). 151, 152, 153, 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, 164, 165
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 simplest 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. 1b, 5b, 5f]. Despite significant attenuation, variations in absorption due to changes in local , , and 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 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 depth-resolved reconstruction of functional activation.174, 175, 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 ), such that multiple wavelengths177, 178 or multiple sources 149, 179, 180, 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 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, 185, 186, 187, 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 .
Calculation of and HbR concentrations. The change in intensity measured between a topographic source and detector (the attenuation) will generally reflect a change in the absorption of underlying tissue [Eq. 1]. The basic approach for deriving changes in and concentration from changes in attenuation was outlined in Sec. 1.1.2 [Eqs. 1, 5]. Determining and 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 [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, 5 are to be used for “direct” topography, an estimate of for each wavelength and measurement pair is required.
FD and TD systems can be used to measure directly, and as an adjunct to cw imaging can provide 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 measurements (or generalized values of DPF) are to be used, it should be noted that 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 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. 1b], and even the effects of isotropic scattering in the brain. 88, 192, 193, 194, 195, 196, 197, 198, 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 .
In general, for direct topography, simplified models and scaling factors are applied to account for 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 and its 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 , heart rate , and the “vasomotor” or Mayer wave signal 207, 208, 209, 210 can 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
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.
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. 5b]. 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 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.is the photon density within the tissue, and is the source function. is the set of changes in measurement (from time ) between source and detector , and is the Jacobian matrix derived from the Greens function solution to Eq. 6, representing the change in measurement corresponding to a change in absorption at position within the tissue.215 An example of for a single source-detector pair is shown in Fig. 5f. This distribution was calculated using a Monte Carlo model of light propagation in a realistic 3-D model of the adult head, based on a segmented MRI scan.202
Analytical solutions to the DA for simple shapes such as slabs and semi-infinite planes exist.197 More complex approaches include applying the DA via finite element or finite difference methods, where heterogeneous structures can be incorporated into the simulation. Monte Carlo modeling can also be used to simulate light propagation in the head, although for large tissue volumes, it can be prohibitively slow. Once the model has been used to generate the Jacobian (otherwise known as the spatial sensitivity matrix, weight matrix, or photon measurement density function), Eq. 7 can be solved for using a linear inversion technique such as Tikhonov regularization, or more complex nonlinear iterative algorithms.217
Plots of [Fig. 5f] reveal the varying paths of light traveling from source to detector, and show how only some of the detected light has sampled the cortex. The areas of highest sensitivity are directly below the source and detector, revealing why changes in coupling and the scalp have such a significant effect on signal amplitude. The banana shape of the light path is characteristic of diffuse imaging, and is maintained as source-detector distances get larger or smaller, such that the center region gets deeper or shallower, respectively.
Implementation of reconstructed topography can be achieved with various levels of complexity. The simplest approach assumes that changes in absorption are occurring in a homogenous volume of known optical properties. This is an approximation, since the baseline structures of the skull and brain will affect the shape of [as shown in Fig. 5f] in a way that cannot be accounted for in a homogenous model.198, 218, 219 The effects of CSF have also been widely considered, since the diffusion approximation is only valid for scattering regions where , which is not likely to be true for large thicknesses of CSF.196, 199 Several papers further suggested that light piping could occur in the CSF layer, significantly reducing penetration of light into the cortex. However, more recent publications have questioned the magnitude of the effect once the complex 3-D nature of the CSF layer, trabeculi, vessels, and refractive index differences are taken into account.193, 194, 217 By incorporating a priori structural information into the image reconstruction scheme, constraints can be applied to help improve the localization of activation to the cortex, although care must be taken to account for changes in the data that may originate from extra-cerebral contributions.202
Data acquisition. As NIR light sources and detectors become increasingly affordable, it has become possible to construct NIRS instruments with large numbers of multiwavelength source and detector channels. While larger arrays allow simultaneous mapping of several regions of the cortex simultaneously, they also present significant challenges for fast data acquisition.180 Since the human heart beat causes substantial variations in NIRS signals,220 aliasing can be a significant problem unless measurements from all source-detectors pairs can be acquired at a rate of at least . Systems that utilize frequency encoding to allow all sources to be illuminated simultaneously can acquire data in parallel at high frame rates.152, 181 However, the attenuation difference between a source and detector separated by 2.5 and 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 frequency-encoded 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
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. 5d]. 145, 221, 222, 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
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. 5c and 5d]. 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 of tissue, which in general will cause around 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 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 absorbers 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 must be calculated based on the 3-D shape of the head, and the precise positions of the sources and detectors [see Fig. 5c]. 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 using the finite element method is called TOAST (time-resolved 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 .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 time-resolved 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
Examples of Noninvasive Optical Brain Imaging
Figure 7 shows a variety of noninvasive optical brain imaging results showing human functional responses to stimulus. Figure 7a 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 , followed by of rest, and were repeated 20 times per subject. The cw system has 16 laser diode sources, eight at and eight at , frequency encoded between 4 and 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. 7b where four fibers span the center (each delivering both 690 and —red arrow), straddled by two rows of four detector fibers—green arrows. The detected signal was digitized at and analyzed to extract signals coming only from the closest sources to each detector. Data were then low-pass filtered and downsampled to , and then bandpass filtered between 0.02 and 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 (6.51) and (5.86), such that the pathlength is , where is the source-detector separation of . The resulting absorption changes were then converted to and 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. 7a show the mean maps and time courses for a single subject. Contralateral responses are those seen in the side of the brain opposite the hand receiving the stimulus (where responses would be expected).
Summary of optimal two-photon microscope characteristics for in-vivo brain imaging.
|Feature||Custom-built in-vivo two-photon microscopy systems||Commercial two-photon microscopy systems|
|Form factor||Almost always upright.||Often inverted.|
|Imaging platform||Fixed stage (objective moves up and down, not stage). Motorized stage. Very open clearance around objective for animal and monitoring equipment.||Require translation of stage, poor clearance between stage and objective. Motorized stage movement very expensive.|
|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.||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.|
|Numerical aperture||Very high numerical aperture (but long working distance) allows good resolution, high efficiency, and ability to image through exposed-cortex preparation [see Fig. 4b].||Similar objectives can be chosen for commercial systems if compatible.|
|Scan rate||Typically designed to allow very high frame or line-scan rates to image dynamics of functional activity.||Frame rate well optimized, but limited adaptability for patterns and line scans.|
|Software||Custom software allows wide variety of scan patterns and routines for a wide range of experiments.||Often proprietary and difficult to modify. Fixed scan patterns and routines.|
|Multimodality||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.||Very difficult to modify or interrupt light path.|
Figure 7b shows a similar study, but where optical topography data were acquired during simultaneous fMRI at .159 The subjects performed a finger tapping exercise for , with a randomized interstimulus interval between 4 and repeated between 27 and 32 times per fMRI scan, with each scan repeated four to six times. fMRI scans with a 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 , , and time courses for one subject and for each source-detector pair are shown in map form [inset lower right of Fig. 7b]. 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 , , and temporal responses measured using NIRS. The group average of five subjects is shown, and the work concludes that BOLD most closely temporally correlates with changes. ASL changes were found to correlate most closely with changes (data not shown). A second paper continues this analysis to a spatial comparison between the localization of and BOLD responses.235
Figure 7c shows a reconstructed topography result, acquired on the adult head during a visual stimulus followed by 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 , and 28 avalanche photodiode detectors. Frequency encoding allowed simultaneous acquisition at both wavelengths, but LEDs were switched on sequentially allowing detectors with dynamic range to measure signals from sources at three different separations: 13.3, 30.0, and . 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 functions calculated using a finite element model simulation package called NIRFAST,236 which incorporates a simplified baseline structure consisting of a scalp/skull layer ( , ) surrounding the brain of ( , ). Both block averaged and even trial-by-trial data show clear hemodynamic responses to unilateral visual stimulus.
Figure 7d 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. 7d show the change in absorption at , but simultaneously acquired images at were also created and converted into maps of the change in , , and .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 while the infant rested. The left arm was then gently moved back and forth by an adult for , with transmitted light being measured for the last of this exercise to capture the resultant hemodynamic response. The next source position was then illuminated, and the rest/ exercise cycle repeated.
The resulting time-resolved data were Fourier transformed and the log(amplitude) and phase components extracted at . 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 functions were generated using a finite element model implementing the diffusion equation in a homogenous 3-D volume corresponding 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
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 at , changing to when oxidized, and 0.00064 to at (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, 241, 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 long-wavelength 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 addition 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, 253, 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
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.
The author acknowledges the help of Matthew B. Bouchard and Sean A. Burgess in proof reading this manuscript, as well as the help and support of many at the Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital. Funding was provided by NIH/NINDS grant R21NS053684, NIH/NHLBI grant R01HL071635, and the Wallace Coulter Foundation.