1.IntroductionToothed whales (odontocetes), such as dolphins and porpoises, use echolocation to hunt prey, navigate, and communicate under water.1 Their brains process information encoded in echoes to extract subtle spatiotemporal cues that allow perception of auditory objects,2 such as prey items, obstacles, or conspecifics. In most mammals, attention is known to modulate what information the brain processes in complex settings with multiple competing sources.3,4 Auditory object forming and attention modulation remains a challenge in humans, especially in crowded environments with multiple sources of sound.5 Understanding how the brain modulates attention in echolocating animals, which rely heavily on sound as a dominant sensor modality in a sensory-challenging and dynamic environment, is therefore offering a unique population to understand this cognitive process. Attempts to measure brain activity in dolphins have been made previously, using modalities such as functional magnetic resonance imaging (fMRI),6,7 measuring hemodynamic responses to neuronal activity, or electroencephalography (EEG),8,9 which measures the electrical field created by coordinated activity of neuronal clusters. fMRI has a high spatial resolution and can therefore distinguish brain activity in neighboring cortical locations but requires aquatic mammals to be removed from the water and transported to a nearby MRI facility. This can limit the variety and effectiveness of conscious attention-based tasks for the animal to perform. Meanwhile, EEG has a high temporal resolution, can be made portable, and can be used underwater if electrodes are insulated from saltwater. However, EEG offers only limited differentiation of brain regions due to a lower spatial resolution. Given the specific limitations of current methodologies, near-infrared spectroscopy (NIRS) may be an additional tool that could improve capacity for measurement of cortical activity in cetaceans. Functional NIRS (fNIRS) is an imaging modality that, like fMRI, measures the hemodynamic response to neuronal activity utilizing the differences in light absorption between oxygenated and deoxygenated blood. Light of specific red to near-infrared wavelengths (typically 700 to 900 nm for human brain measurements) is shown onto the scalp from where it is scattered diffusely into the tissue, reaching the cortex of the brain. A fraction of this diffuse light scatters back to a photo-detector placed on the skin multiple centimeters away from the source. Changes in light intensity can then be correlated to changes in hemoglobin concentrations, using the modified Beer–Lambert’s law.10 NIRS devices can be made portable and water resistant,11–13 and high-density fNIRS has previously been shown to reach fMRI-comparable spatial resolution in human measurements if tomographic reconstruction using multiple spatially separated measurements is performed.14 However, the anatomical and morphological differences of dolphins pose key challenges to the function of NIRS. Challenges such as light absorption in muscle and subcutaneous fat (blubber) tissue above the dolphin brain need to be addressed. The primary challenge of using fNIRS in dolphins is the penetration depth required to reach the brain. Although the skin to cortex distance in human adults is to 20 mm,15 the thicker skin and additional muscle and blubber layer seen in adult bottlenose dolphins increases this distance to 40 mm and more. Approaches to improve the optical penetration depth of NIRS have been reported before. For example, work on breast cancer detection has utilized anatomical reconstruction based on many source and detector combinations (called diffuse optical tomography), benefited from a cylindrical shape of the target tissue16 or compression of the tissue17–19 to overcome penetration depth challenges. In transabdominal fetal pulse-oximetry, large source–detector separations in combination with higher light intensities are utilized, as well as more sensitive detectors and advanced signal processing.20–22 To evaluate the utility of fNIRS for brain imaging and NIRS for physiological measurements in dolphins, we simulated light propagation in a 3D computational dolphin model and determined optimal hardware parameters. This required the creation of an anatomical dolphin model with realistic optical properties of light absorption, scattering, and refraction. Optical properties for terrestrial mammals can be found in the literature, such as Jacques,23 for bone, muscle, fat, skin, and other tissues. However, optical properties for marine mammals, especially with respect to their unique blubber tissue, have not been sufficiently characterized. Here we use optical properties of dolphin blubber and skin measured with frequency-domain NIRS and simulate results of light propagation in a 3D model of a dolphin. We also show expected brain tissue sensitivity (TS) as a function of wavelengths, source–detector distances, and optode placement. Our results indicate that functional brain measurements are challenging and would require specialized equipment. We therefore contrasted the findings in brain TS to muscle TS as an alternative application for NIRS in marine mammal physiology. 2.Simulation of Photon Travel in Dolphin Tissue2.1.Optical Properties of Dolphin TissueWe created a dolphin head model that was divided into four tissue types: “brain,” “bone,” “blubber and skin” (sculp), and “muscle and other tissues.” Optical properties of absorption coefficient (), reduced scattering coefficient (), and the refraction index () were defined for each tissue type. We calculated and as a function of wavelength using literature values reported by Jacques24 for brain, bone, and muscle tissues. For the marine mammal–specific blubber tissue, optical properties were measured using a combination of frequency domain NIRS (FD-NIRS) and broadband NIRS. For the optical properties of blubber and skin, we measured bottlenose dolphins (Tursiops truncatus) at the Fundación Oceanogràfic de la Comunitat Valenciana, Valencia, Spain; the Siegfried and Roy’s Secret Garden and Dolphin Habitat at the Mirage Las Vegas Hotel & Casino, Las Vegas, Nevada, United States; and Dolphin Quest Oahu, Honolulu, Hawaii, United States. All research procedures were approved by the Oceanogràfic Animal Care & Welfare Committee at the Fundación Oceanogràfic Valencia, Spain (OCE-10-20), and the Bureau of Medicine (BUMED, NRD1170). The animals were trained to either slide onto a ledge or remain stationary inside the water where they rested voluntarily for up to 3 min, during which optical properties were measured. The animals had free access to the pool at any time and no restraints were used to hold animals in place, thus animals could leave the measurement site at any time. No sign of discomfort was noticed in the animal’s behavior while they were exposed to the measurement devices. An FD-NIRS system (Imagent 2.0, ISS, Inc., Champaign, Illinois, United States) was used to measure and at seven different wavelengths in the near-infrared range between 730 and 830 nm. An optical probe with source and detector fibers spaced at 9, 18, 27, and 36 mm was placed on the dolphin’s skin between the pectoral fins on the chest. Light intensity was measured with optical fibers connected to four individual photo-multiplier tubes. Given the source–detector distances, the penetration depth of the light was not expected to exceed the deep blubber layer. Therefore, the measured optical properties are assumed to represent only a combination of blubber and skin pigmentation. The light was modulated at 130 MHz to enable the calculation of phase delays between the source and detected light, alternating current (AC) component, and direct current (DC) component. Using the multi-distance approach for optical property calculation,25,26 we determined the average and across the measurement period. We recorded a total of 117 measurements on 18 animals, in which only the 51 chest measurements were considered for this work due to a greater similarity in pigmentation across animals in this location. Every measurement was inspected for signal quality by ensuring that the modulation percentage was , no motion artifacts were observed, and no detector was saturated. Individual time frames and data channels showing data quality deficiencies were removed. In the end, a total of 41 measurements across 10 animals (2 Valencia + 2 Las Vegas + 6 Hawaii) were included in this study. Individual time averages of and were calculated for each wavelength and measurement. The median across all measurements was calculated to estimate the skin and blubber optical properties. The median excluded any outlier values exceeding the interquartile range by more than 1.5 times the interquartile range. In addition, for a subset of animals, broadband absorbance spectra were measured with a photo-spectrometer (FLAME-T-VIS-NIR-ES, Ocean Insight, Inc., Orlando, Florida, USA) and a tungsten halogen light source (HL-2000-HP-FHSA, Ocean Insight, Inc., Orlando, Florida, United States). Source and detector fibers were placed on the skin close to the location of the FD-NIRS measurements with a source–detector distance of 15 to 18 mm. With integration times of 4 to 20 s, depending on the environmental conditions and absorbance of the dolphin tissue, the received light was accumulated and saved over a wavelength range of 350 to 1000 nm. A total of 51 trials (32 on the chest) on a total of 6 animals were recorded. The measurement duration of a trial was 2 min, comparable to FD-NIRS measurements, however, the sampling rate depended on the integration time and the number of samples is therefore inconsistent. After rejection of saturated samples, often due to motion or influence of ambient light occurring within the long integration times, all samples were smoothed using a boxcar of 15 indices, corresponding to about 4.2 nm wavelength intervals, and averaged. The average absorbance across all remaining 18 measurements was calculated. To obtain absolute , the reflectance in the wavelength range of 700 to 900 nm was fitted to the FD-NIRS-measured of the blubber tissue. We used a least square fitting algorithm (“lsqcurvefit,” MATLAB 2022a, The MathWorks, Inc., Natick, Massachusetts, United States) to estimate oxygenated and deoxygenated hemoglobin concentrations, water (), melanin (), and fat () based on the spectrogram between 700 and 900 nm. For hemoglobin concentrations, the weights represented the blood volume fraction () relative to whole blood and the tissue saturation with oxygenated hemoglobin () in accordance with the nomenclature in Ref. 23. Hemoglobin concentrations in whole blood were assumed to be . A lower bound of 20% fat and 10% water was set to prevent fitting to local minima not representing blubber histology. The scattering coefficient follows a power law , where is the wavelengths of light and and are coefficients that were fitted to match the measured by the FD-NIRS system. The optical properties for brain, bone, and muscle tissue were taken from the literature. Specifically, and were taken from Ref. 24. Using the tissue specific concentration values from Ref. 24, the overall tissue absorption can be calculated through weighted summation of .27,28 Thus the total tissue specific absorption coefficient is defined as The scattering coefficient over wavelength was reconstructed using the power law coefficients provided by Jacques.24 Finally, were taken from the literature with brain24 and muscle29 , blubber30 with a higher lipid concentration was set to , and bone31 . 2.2.Finite-Element ModelA 3D head model was created based on a computed-tomographic (CT) scan [Fig. 1(a)] of a stranded striped dolphin (Stenella coeruleoalba), provided by the Fundación Oceanogràfic de la Comunitat Valenciana, Valencia, Spain. The 3D CT scan contained the full circumference of the dolphin from the tip of the rostrum to the third thoracic vertebra. Each voxel had a size of . The CT scan was segmented by hand using the image processing and segmentation program 3D Slicer.32,33 After segmentation, a 3D image () was saved containing the segmented tissues separated into five categories: (1) brain, (2) bone, (3) blubber and skin, and (4) muscle and other tissues, while the background was given a default value of 0—air [Fig. 1(b)]. No voxel of the CT scan was left unassigned, and no air was assigned to a voxel inside the volume. MATLAB R2022a (The MathWorks, Inc., Natick, Massachusetts, United States) was used for processing, and the image was loaded using Fieldtrip34 function “ft_read_mri.”34 The segmented image was converted into a tetrahedral mesh using the Computational Geometry Algorithms Library’s 3D triangulation algorithm.35,36 The generated 82,512 nodes created 467,648 tetrahedrons [Figs. 1(c) and 1(d)]. The average tetrahedral volume was determined to be (mean ± standard deviation). Every node inside the generated mesh was assigned , diffusion coefficient , and . 2.3.Photon SimulationThe diffusion equation describing the propagation of light in scattering media was solved at every node in the mesh to calculate the fluence rate of light at every point, using the NIRFASTer toolbox.37–39 The light sources and detectors were placed on the animal by finding the edge between dolphin skin and surrounding air in every labeled slice ( images along the axis). Along the axis, 11 equidistant slices were chosen for optode placement, as seen in Fig. 1(d). Along the skin, 21 radially equidistant optodes were placed on every slice. After removing optodes to ensure that no pair fell within 10 mm of one another, a total of 222 optodes were alternated in sources and detectors, leading to an equal number of 111 sources and 111 detectors. Every source–detector pair within 50 to 120 mm of each other created a data channel. The resulting 2746 channels created a dense grid around the head. Using the NIRFASTer toolbox,37–39 the forward model was calculated for continuous wave (CW) NIRS. The light projection was calculated for every source optode and the detector amplitude for every source–detector combination, or channel, given by the fluence rate at the detector position on the surface of the mesh assuming that only one source is active at a time. All fluence rates and amplitudes were computed relative to the input intensity of . The optical properties of all tissues were reassigned, and the calculation was repeated for all wavelengths in the target range of 600 to 1700 nm, in 50 nm intervals. After the forward modeling, the Jacobian matrix was derived.40 The Jacobian (mm) is a sensitivity matrix that quantifies the fluence rate change in the detector given an absorption change (of ) in every node in the mesh for a specific source and detector combination, i.e., Thus is the partial derivative of the fluence rate for a change in absorption coefficient . The influence of the change in scattering is assumed to be negligible under the assumptions of a CW illumination. 2.4.Quantitative Comparison and StatisticsTo determine optimal source and detector positions, distances, and wavelengths, the following sensitivity measurements were calculated for wavelengths ranging from 600 to 1700 nm, in steps of 50 nm, varying the optical properties of the modeled tissue. To calculate the sensitivity to brain or muscle tissue, we utilized separate metrics: (1) Tissue sensitivity (TS) was calculated as the difference between simulated amplitudes before and after a local increase in oxygenated hemoglobin inside the target tissue. (2) Depth sensitivity based on the channel-specific sensitivity matrix . (3) Depth sensitivity based on a flatfield change in absorption, set to 1% across the entire mesh. This analysis utilizes back projection to find the location of absorption coefficient change in the mesh and is performed within groups of channels with similar source–detector distance. (4) Light power inside a simulated large detector covering a . 2.4.1.Tissue sensitivityFor the TS measurement, we ran a separate forward model, in which the hemoglobin content of the target tissue, here the brain, was increased by 10%, i.e., . We calculated the TS as the difference between the two forward models, such that where is the fluence rate measured at the detector for every channel before the brain activation, and is the fluence rate during activation. To ensure no channel has a uniquely and unrealistically high sensitivity to the target tissue and reduce the effect of numerical inaccuracies and poor mesh quality, we averaged all following sensitivity and intensity metrics across the highest 1% TS.2.4.2.Depth sensitivity: JacobianThe Jacobian or sensitivity matrix shows the influence of brain tissue to measured light intensity. We therefore looked at the highest absolute value of the Jacobian inside the brain sub-mesh, relative to the highest absolute value of the Jacobian in the entire mesh. This metric indicates the total sensitivity of a specific channel to the target tissue. 2.4.3.Depth sensitivity: flatfieldThe flatfield analysis of depth sensitivity was based on work published by Dehghani et al.38 A flatfield, or uniform, increase of 1% in absorption is used to probe the influence of every point in the mesh on a given cluster of channels by where is the change in absorption and is the resulting change in fluence rate at the detector. Note that this inverse model is a tomographic approach to determine the depth in which the highest absorption change occurs. We compared source–detector distances by binning all channels into seven 10-mm-wide groups of source–detector distances between 50 and 120 mm. Using the back projection of the light propagation model, which calculates for every node in the mesh based on measured detector intensities (here calculated as a flatfield), we answer the question of where in the mesh a change in absorption occurred. The back projection requires an inversion of , defined according to Dehghani et al.38 and Culver et al.40 using the Moore–Penrose pseudoinverse41 as where is the Tikhonov regularization parameter applied to the identity matrix . Further,The regularization of the inversion algorithm is done in two parts, is a depth-dependent spatial regularization defined as with and is a spatial regularization defined as with2.4.4.Large detector placementTo determine if a detector can measure a subtle change in absorption, given a large source–detector-distance, we simulated placing large detectors with a total area of onto the mesh’s surface. We then compared the fluence rate measured in these detectors before and during a simulated activation by increasing the hemoglobin content inside the brain by 10%. The detector size was simulated by placing a square grid of simulated sub-detectors () around the original detector location. The distance between sub-detectors was set to 0.2 mm. All sub-detectors amplitudes were calculated as described in 2.3. The total large detector fluence rate, or power, summed over the detector area was then approximated with A detector must achieve a high dynamic range to not saturate at the maximum expected power , while also being sensitive to the subtle difference in fluence rate between the activated brain and the brain at rest . A small range between and is therefore desirable for real-life detector selection. 2.5.Juvenile 3D ModelAll analysis was repeated for a “juvenile” animal model, in which all dimensions of the model were shrunk by 25%. Although the number of nodes and elements in the mesh is unchanged, the average tetrahedral volume decreased to (mean ± standard deviation). The same amount of optodes on a smaller mesh resulted in a greater overlap within the defined 50 to 120 mm source–detector separation, totaling 4214 channels. All data post-processing and statistical analysis were repeated for the reduced mesh size and compared to the normal sized “adult” model. 3.ResultsOptical properties measured from the chest of 10 voluntarily participating Atlantic bottlenose dolphins in three different professional care facilities were averaged to obtain a general blubber and skin [Fig. 2(a)] and [Fig. 2(b)]. The measured absorption properties of the blubber tissue were lower than highly perfused brain and muscle tissues, but at some wavelengths higher than bone optical properties taken from the literature [Fig. 2(c), black line]. The scattering properties measured with the FD-NIRS device were, overall, decreased with wavelength, and when compared to other tissues [Fig. 2(d)] were overall low. The fitting results for the measured blubber and skin layer are shown in Fig. 2(e), where the black line shows the broadband photo-spectrometer measured absorbance after it had been fitted to the FD-NIRS measured (red asterisk). The least squares fit of tissue components then generated the red line. The weights for the fitting were determined to be , , , , and . Similarly, the blue asterisks show the measured , which were used to fit the blue line according to the power law. The scaling factor and exponent to the scattering fit were determined to be and . At a typical NIRS operation wavelengths of 850 nm, the muscles and brain are the most absorbing tissues [Fig. 2(f)], whereas the bone tissue is the most scattering component [Fig. 2(g)]. To fully understand the requirements for dolphin fNIRS measurements, we observed how light travels inside the dolphin’s body. We found that given a light source fluence rate of , the light fluence rate decreased rapidly [see Figs. 3(a) and 3(b)], with a single light source losing as much as eight orders of magnitude in fluence rate before reaching the brain cortex. A circular arrangement of light sources can evenly illuminate the brain cortex for a full coverage of from the midsagittal plane [Fig. 3(b)]. To find the depth and tissue of most influence to the signal, we calculated the flatfield metric. At 60 to 70 mm source–detector separation, blubber tissue was the most sensitive area [Fig. 3(c)]. With a wider source–detector separation (e.g., 110 to 120 mm), we see the strongest influence on the measured light intensities originate from the border of blubber to muscle tissue layer [Fig. 3(d)]. To see how strong the contribution of the brain tissue is, we looked at the sensitivity matrix (). Figures 3(e) and 3(f) illustrate the sensitivity as a percentage of the strongest sensitivity in the entire mesh. At a 60 to 70 mm optode separation, we achieve a sensitivity of . A larger source–detector distance can increase the sensitivity to the brain to about 0.01% [Fig. 3(f)]. The sensitivity increases when simulating a juvenile dolphin—shrinking the mesh size by 25%. Fluence rates in the brain are [Figs. 4(a) and 4(b)], with an increased sensitivity to deeper layers of the muscle tissue, where the maximum sensitivity borders on the skull [Fig. 4(d)]. The Jacobian confirms these improvements, with a sensitivity of 0.1% in the shorter source–detector distance [Fig. 4(e)] and a sensitivity of with long source–detector distances [Fig. 4(f)]. We compared the adult and juvenile models for different quantitative measurements. The most sensitive 1% of channels was found based on the TS metric. When comparing a mean across the 1% most sensitive channels to brain tissue based either on the maximum inside the brain mesh or the maximum TS index, we found a more than two orders of magnitude increase in sensitivity [Figs. 5(a) and 5(b)] for the juvenile model compared to the adult model. The most sensitive wavelength is between 700 and 900 nm in the adult model and 1400 and 1600 nm in the juvenile model. Because this sensitivity metric does not consider the total amount of light reaching the detector, we also looked at the detector amplitudes of the most sensitive channel. The total measured power in a detector was calculated to be with respect to a light source fluence rate of at 1100 nm and slightly lower at 700 to 900 nm. The power decreases strongly as wavelength increases beyond 1100 nm [Fig. 5(c)]. The detector power difference between the activated brain and the brain at rest showed an intensity difference of at 800 to 900 nm and 1050 to 1100 nm. The juvenile model, while receiving less absolute light intensities, showed a of at 800 to 900 nm, and thus an order of magnitude improvement [Fig. 5(d)]. Both juvenile and adult meshes showed better sensitivity at the longest source–detector distance [Fig. 5(e)]. Another practical consideration for selection of a detector capable of measuring the brain activity is the dynamic range of intensities it must cover. The juvenile dolphin model had greater brain sensitivity with slightly lower absolute fluence rates and thus requires a smaller dynamic range. We estimated that the required dynamic range was about for wavelengths in the ranges of 800 to 900 nm and 1050 to 1100 nm in an adult dolphin model; it was about in the juvenile model [Fig. 5(f)]. The best source–detector distance can be found from the above results by considering one of the most sensitive wavelengths, 850 nm. For this wavelength, we plotted the sensitivity and optical power metrics as a function of optode separation. In both the adult and juvenile models, we found that the largest source–detector separations (110 to 120 mm) yield the highest sensitivity to the brain tissue; however, these source–detector distances show lower changes (see Fig. 6). Finally, to find the optimal location for optode placement, we plotted the maximum Jacobian inside the brain for every channel at the position halfway between the source and detector. The location most sensitive to changes in brain signal was located between 45 and 90 deg from the midsagittal plane, about 30 to 50 mm behind the blowhole (see Fig. 7). Interpretation of the data with respect to brain sensitivity has led to discussions on sensitivities of other tissues and the utility of NIRS as a measurement of physiological parameters, such as tissue saturation with oxygenated hemoglobin or heart rate dynamics. We thus investigated the sensitivity to muscle tissue by changing the hemoglobin concentration inside the muscle tissue by 10%, comparable to the brain activation in Fig. 5. A much higher sensitivity overall and dependency as a function of wavelength was found in the TS index [Fig. 8(b)]. Despite comparable absolute power measurements in the simulated detector [Fig. 8(c)], given similar source–detector separations, the change in power as a result of the total hemoglobin concentration increase was much greater than previously observed in the brain [Fig. 8(d)], leading to a greatly reduced dynamic range required to measure these signal [Fig. 8(f)]. 4.DiscussionThe main goal of this work was to determine the sensitivity of NIRS to functional activation in a dolphin’s brain within practical limits, to explore the potential utility of fNIRS to measure brain activity in medium sized cetaceans. We measured optical properties of dolphin blubber and combined it with literature values of optical properties from muscle, bone, and brain tissue. We simulated light travel through a 3D model of a striped dolphin and extracted brain sensitivity metrics. Overall, we found that the subcutaneous blubber layer and muscle layer absorption and the increased distance to the brain as compared to human fNIRS applications, makes functional brain measurements in dolphins challenging. We therefore investigated the influence of animal size, by simulating a juvenile animal through 3D mesh size reduction. We further laid out sensitivity metrics for muscle tissue to discuss the use of NIRS as a tool for animal healthcare and research. These findings are further discussed below. We estimated that, as a percentage of the maximum value in the Jacobian matrix, sensitivity to brain tissue in adult bottlenose dolphin falls between 0.01% and 0.001% for source–detector separations of 110 to 120 mm [see Figs. 3(e) and 3(f)]. In contrast, for a human model,38 a brain sensitivity of can be achieved with source–detector separations of . Given anatomical constraints in adult dolphins, a source–detector separation of up to 120 mm cannot achieve a signal quality comparable to human models. However, we have shown some sensitivity to signals deriving from the brain (Figs. 3 and 5). In a setup with external stimuli or cognitive tasks, data quality like that achieved in human research would require a larger number of trials to increase the signal-to-noise ratio and average out physiological artifacts. The sensitivity is improved by an order of magnitude for a juvenile dolphin or a smaller dolphin species (simulated here with a 25% reduction in size). For such smaller animals, we can achieve comparable results to human fNIRS with a sensitivity close to or exceeding 1% at source–detector distances (Fig. 4). These sensitivity percentages assume a wavelength of light of 850 nm, which is one of the most sensitive wavelengths to brain tissue [Fig. 5(b)]. It is also a wavelength used in many commercially available NIRS devices for use in humans. Yet it is not the only wavelength that is suitable for fNIRS applications in dolphins. High sensitivity to brain tissue, as measured by the TS index, was also found for wavelengths between 700 and 900 nm, 1050 and 1100 nm, and 1400 and 1550 nm [Fig. 5(b)], which is reflected in the sensitivity matrix [Fig. 5(a)]. Despite high sensitivity in the longer wavelength range, we see a steep decline in measured detector amplitudes for wavelength . This is likely due to the increase in water absorption, which is a major contributor in most tissues. Combining the TS to the brain with the returned light power, optimal wavelengths for imaging the brain are between 700 and 900 nm or 1050 and 1100 nm. The gap at 950 and 1000 nm is likely attributed to high fat absorption (see Fig. 2). Selection of a suitable detector depends both on the minimum required light sensitivity to pick up brain activations and the largest expected light intensity, which together determine the required dynamic range of the detector. Increasing a sensor’s dynamic range tends to increase its cost, which is an important practical limitation. We found that the wavelengths that maximized the light received and achieving the highest brain TS also required a smaller dynamic range. The smallest dynamic range required in the adult model was ; the juvenile model required a dynamic range two orders of magnitude smaller, i.e., [see Fig. 5(f)]. This dynamic range was based on results from the 1% of channels that were most sensitive to signals originating from the brain according to the TS [Fig. 5(a)], with source–detector separations between 110 and 120 mm. Large source–detector separations are necessary for enough light to reach the brain but lead to low received light intensity and require sensors to have a large dynamic ranges. However, the optimal compromise across these constraints depends on the size of the animal as the dynamic range decreases by a factor of 100 and the light power difference to brain activation increases by a factor of 10 when imaging a juvenile dolphin [Fig. 5(d)]. The optimal source–detector location for fNIRS depends jointly on the amount of light received in the detector and the sensitivity of the received light to signals from the brain (Fig. 6). Selection of the optode separation has a tradeoff, and the optimal parameters vary with the size of the animal. For our adult dolphin model, a source–detector separation of is within one order of magnitude of that of the most sensitive channel and can create a trade-off between received light intensity and source–detector separation. This distance will increase for larger animals and can be reduced for smaller ones. In addition to wavelength and source–detector separation, we also explored how optode placement affected sensitivity. The Jacobian matrix can tell us about the sensitivity at a specific location, i.e., the contribution of a specific node inside the mesh to the measured light intensity in any channel. The highest sensitivity to brain tissue, according to the largest Jacobian values, was found in channels covering an area 30 to 50 mm behind the blowhole at 45 to 90 deg from the midsagittal line (Fig. 7). This location refers to the midpoint of the source–detector distance. Although the overall sensitivity changed with wavelength, the location of best sensitivity remained comparable. Although brain fNIRS measurements appeared challenging with today’s commercial NIRS devices, physiological measurements of muscle and blubber tissue for behavioral research and veterinary purposes appear feasible. We have previously demonstrated the ability to measure cardiac pulsation and respiration from a bottlenose dolphin’s chest13 at separations . Our simulation results confirm these findings. A source–detector distance of 60 mm shows a high sensitivity to the muscle layer underneath the blubber in both adult and juvenile dolphins, as shown by the signal percentage of the flatfield images [see Figs. 3(c) and 4(c), respectively]. The muscle sensitivity is increased greatly when comparing brain [Fig. 3(b)] to muscle tissue [Fig. 8(b)] within the same model. In fact, a 10% total hemoglobin concentration increase in the brain reveals a sensitivity more than three orders of magnitude lower than in muscle tissue. This trend is reflected in the dynamic range between activated and rested tissues. Our model therefore suggests that muscle tissue measurements in adult dolphins are possible with commercially available NIRS devices, whereas brain fNIRS measurements will require highly sensitive, specialized equipment. The models and assumptions used here represent our best good faith estimates of tissue properties and anatomical structures in dolphins; however, some caveats should be highlighted. To simulate how sensitive a NIRS measurement on the surface would be to hemodynamic signals originating from deeper tissues, we first needed to create a finite-element model and assign optical properties to different layers. The model created here is based on the CT scan [see Fig. 1(a)] of a stranded, deceased dolphin. In creating the segmentation, some tissue had to be interpolated. The segmentation was compared to information about dissected animals from the literature to ensure a realistic allocation of tissue groups, which at times required the interpolation of smaller air pockets created by the decay of connective tissues. These pockets were interpolated with surrounding muscle tissue. The anatomical structure underlying all simulations is only based on a single individual. Subject variations as a function of age, sex, season, and location are beyond the scope of this work; however, they may play an important role in the interpretation of this data. For example, the thickness of the blubber layers can vary between summer and winter season. Additionally, a juvenile animal might have anatomical differences, which are not captured by simple miniaturization of an adult 3D model. Each node of the mesh generated from the segmentation is a point in the finite-element model to which we assigned light absorption, scattering, and refraction properties. The assumed values of these parameters were taken from a review by Jacque,24 most of which came from terrestrial mammal tissue measurements. These optical properties were used for the brain, bone, and muscle layers. Light-absorbing tissue in terrestrial mammals may have a different composition from that in marine mammals. For instance, proteins like myoglobin tend to be more common in the muscle tissue of deep diving marine mammals than in terrestrial mammals. Finally, we limited the source–detector separation to 120 mm. Further increasing the source–detector separation is theoretically possible to increase the brain sensitivity, but the returned light intensity would be further decreased. Additionally, a larger source–detector separation will decrease the spatial resolution of fNIRS and limit the purpose for observation of specific cortical activation. Although we have shown that brain imaging in dolphins with NIRS is in principle possible at 120 mm source–detector distance, further evaluation of practicality will be needed. 5.ConclusionWe have shown that with our assumptions, achieving brain fNIRS measurements in adult bottlenose dolphins will be challenging. Signal strengths of brain activations are close to the limit of what even a specialized maximum light power, maximum sensitivity NIRS device could theoretically achieve. Nonetheless, fNIRS might be possible in specific locations at source–detector distances and at wavelength from 700 to 900 nm or 1100 to 1150 nm. We also found that smaller dolphins, such as juvenile animals or potentially a smaller species, can greatly improve the received signal from the brain. Our simulations also showed that in both adult and juvenile models, sensitivity to muscle tissue layers dominated measurements with shorter source–detector separations of 50 mm and likely below. Still, development of NIRS to study cerebral physiology in dolphins could yield a valuable tool in animal care and help unravel adaptations that allow these marine mammals to perform extreme breath-hold diving. Although muscle perfusion and oxygenation can be measured using off-the-shelf NIRS devices, functional brain measurements in dolphins will require development of specialized tools that deliver and sense greater amounts of light than in human fNIRS to support the large source–detector separations necessary to reach the brain. DisclosuresThe authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. AF was employed without salary by the company Global Diving Research SL. Code, Data, and Materials AvailabilityThe work presented here utilizes the NIRFASTer toolbox, which is publicly available on GitHub: https://github.com/nirfaster/NIRFASTer.git. The 3D model, optical properties, and code required to reproduce the analysis discussed in this article will be made available under doi: 10.24433/CO.6763388.v1. AcknowledgmentsWe would like to thank all trainers and staff at the Fundación Oceanogràfic de la Comunitat Valenciana, Valencia, Spain, and the Siegfried & Roy Secret Garden and Dolphin Habitat at the Mirage Hotel and Casino in Las Vegas, Nevada, United States, for their contribution to the data collection in dolphins. We further would like to acknowledge Daniel García Párraga of the Fundación Oceanogràfic de la Comunitat Valenciana for providing a bottlenose dolphin CT scan. The authors also would like to thank Hamid Dehghani for advice on analyzing and summarizing model results. ReferencesB. Würsig, Ethology and Behavioral Ecology of Odontocetes, 483
–494 1st ed.Springer Nature Switzerland AG, Cham, Switzerland
(2019). Google Scholar
M. W. Muller et al.,
“Time-frequency analysis and modeling of the backscatter of categorized dolphin echolocation clicks for target discrimination,”
J. Acoust. Soc. Am., 124 657
–666 https://doi.org/10.1121/1.2932060 JASMAN 0001-4966
(2008).
Google Scholar
R. Desimone and J. Duncan,
“Neural mechanisms of selective visual attention,”
Annu. Rev. Neurosci., 18 193
–222 https://doi.org/10.1146/annurev.ne.18.030195.001205 ARNSD5 0147-006X
(1995).
Google Scholar
B. G. Shinn-Cunningham,
“Object-based auditory and visual attention,”
Trends Cogn. Sci., 12 182
–186 https://doi.org/10.1016/j.tics.2008.02.003 TCSCFK 1364-6613
(2008).
Google Scholar
H. An et al.,
“Springer Handbook of Auditory Research the Auditory System at the Cocktail Party,”
http://download.springer.com/static/pdf/730/bok%253A978-3-319-51662-2.pdf?originUrl=http%3A%2F%2Flink.springer.com%2Fbook%2F10.1007%2F978-3-319-51662-2&token2=exp=1497617002~acl=%2Fstatic%2Fpdf%2F730%2Fbok%25253A978-3-319-51662-2.pdf%3ForiginUrl%3Dhttp%25
(2017).
Google Scholar
D. S. Houser et al.,
“Structural and functional imaging of bottlenose dolphin (Tursiops truncatus) cranial anatomy,”
J. Exp. Biol., 207 3657
–3665 https://doi.org/10.1242/jeb.01207 JEBIAM 0022-0949
(2004).
Google Scholar
S. Ridgway et al.,
“Functional imaging of dolphin brain metabolism and blood flow,”
J. Exp. Biol., 209 2902
–2910 https://doi.org/10.1242/jeb.02348 JEBIAM 0022-0949
(2006).
Google Scholar
Y. Yu et al.,
“A portable waterproof EEG acquisition device for dolphins,”
Sensors, 21 3336 https://doi.org/10.3390/s21103336 SNSRES 0746-9462
(2021).
Google Scholar
M. D. Schalles et al.,
“Auditory oddball responses in Tursiops truncatus,”
JASA Express Lett., 1 081202 https://doi.org/10.1121/10.0005991
(2021).
Google Scholar
A. Sassaroli and S. Fantini,
“Comment on the modified Beer–Lambert law for scattering media,”
Phys. Med. Biol., 49 N255
–N257 https://doi.org/10.1088/0031-9155/49/14/N07 PHMBA7 0031-9155
(2004).
Google Scholar
J. C. McKnight et al.,
“Shining new light on mammalian diving physiology using wearable near-infrared spectroscopy,”
PLoS Biol., 17 1
–20 https://doi.org/10.1371/journal.pbio.3000306
(2019).
Google Scholar
J. C. McKnight et al.,
“Shining new light on sensory brain activation and physiological measurement in seals using wearable optical technology,”
Philos. Trans. R. Soc. B, 376 20200224 https://doi.org/10.1098/rstb.2020.0224 PTRBAE 0962-8436
(2021).
Google Scholar
A. Ruesch et al.,
“Near-infrared spectroscopy as a tool for marine mammal research and care,”
Front. Physiol., 12 816701 https://doi.org/10.3389/fphys.2021.816701 FROPBK 0301-536X
(2022).
Google Scholar
A. T. Eggebrecht et al.,
“A quantitative spatial comparison of high-density diffuse optical tomography and fMRI cortical mapping,”
Neuroimage, 61 1120
–1128 https://doi.org/10.1016/j.neuroimage.2012.01.124 NEIMEF 1053-8119
(2012).
Google Scholar
M. G. Stokes et al.,
“Simple metric for scaling motor threshold based on scalp-cortex distance: application to studies using transcranial magnetic stimulation,”
J. Neurophysiol., 94 4520
–4527 https://doi.org/10.1152/jn.00067.2005 JONEA4 0022-3077
(2005).
Google Scholar
B. Pogue et al.,
“Instrumentation and design of a frequency-domain diffuse optical tomography imager for breast cancer detection,”
Opt. Express, 1 391 https://doi.org/10.1364/OE.1.000391 OPEXFF 1094-4087
(1997).
Google Scholar
J. M. Cochran et al.,
“Hybrid time-domain and continuous-wave diffuse optical tomography instrument with concurrent, clinical magnetic resonance imaging for breast cancer imaging,”
J. Biomed. Opt., 24 051409 https://doi.org/10.1117/1.JBO.24.5.051409 JBOPFO 1083-3668
(2019).
Google Scholar
E. Y. Chae et al.,
“Development of digital breast tomosynthesis and diffuse optical tomography fusion imaging for breast cancer detection,”
Sci. Rep., 10 13127 https://doi.org/10.1038/s41598-020-70103-0 SRCEC3 2045-2322
(2020).
Google Scholar
C. M. Robbins et al.,
“Two-layer spatial frequency domain imaging of compression-induced hemodynamic changes in breast tissue,”
J. Biomed. Opt., 26 056005 https://doi.org/10.1117/1.JBO.26.5.056005 JBOPFO 1083-3668
(2021).
Google Scholar
A. Zourabian et al.,
“Trans-abdominal monitoring of fetal arterial blood oxygenation using pulse oximetry,”
https://www.spiedigitallibrary.org/terms-of-use
(2000).
Google Scholar
K. B. Gan, E. Zahedi and M. A. Mohd Ali,
“Transabdominal fetal heart rate detection using NIR photopleythysmography: instrumentation and clinical results,”
IEEE Trans. Biomed. Eng., 56 2075
–2082 https://doi.org/10.1109/TBME.2009.2021578 IEBEAX 0018-9294
(2009).
Google Scholar
D. D. Fong et al.,
“Design and in vivo evaluation of a non-invasive transabdominal fetal pulse oximeter,”
IEEE Trans. Biomed. Eng., 68 256
–266 https://doi.org/10.1109/TBME.2020.3000977 IEBEAX 0018-9294
(2021).
Google Scholar
S. L. Jacques,
“Optical properties of biological tissues: a review,”
Phys. Med. Biol., 58 R37
–R61 https://doi.org/10.1088/0031-9155/58/11/R37 PHMBA7 0031-9155
(2013).
Google Scholar
S. L. Jacques,
“Erratum: optical properties of biological tissues: a review,”
Phys. Med. Biol., 58 5007
–5008 https://doi.org/10.1088/0031-9155/58/14/5007 PHMBA7 0031-9155
(2013).
Google Scholar
S. Fantini et al.,
“Frequency-domain multichannel optical detector for noninvasive tissue spectroscopy and oximetry,”
Opt. Eng., 34 32 https://doi.org/10.1117/12.183988
(1995).
Google Scholar
S. Fantini et al.,
“Non-invasive optical monitoring of the newborn piglet brain using continuous-wave and frequency-domain spectroscopy,”
Phys. Med. Biol., 44 1543
–1563 https://doi.org/10.1088/0031-9155/44/6/308 PHMBA7 0031-9155
(1999).
Google Scholar
N. Bosschaart et al.,
“A literature review and novel theoretical approach on the optical properties of whole blood,”
Lasers Med. Sci., 29 453
–479 https://doi.org/10.1007/s10103-013-1446-7
(2014).
Google Scholar
S. Prahl,
“Optical absorption of hemoglobin,”
https://omlc.org/spectra/hemoglobin/index.html
(1999).
Google Scholar
D. K. Tuchina et al.,
“Quantification of glucose and glycerol diffusion in myocardium,”
J. Innov. Opt. Health Sci., 8 1541006 https://doi.org/10.1142/S1793545815410060
(2015).
Google Scholar
F. P. Bolin et al.,
“Refractive index of some mammalian tissues using a fiber optic cladding method,”
Appl. Opt., 28
(12), 2297
–2303 https://doi.org/10.1364/AO.28.002297 APOPAI 0003-6935
(1989).
Google Scholar
A. Ascenzi and C. Fabry,
“Technique for dissection and measurement of refractive index of osteones,”
http://rupress.org/jcb/article-pdf/6/1/139/1071808/139.pdf
().
Google Scholar
A. Fedorov et al.,
“3D Slicer as an image computing platform for the quantitative imaging network,”
Magn. Reson. Imaging, 30 1323
–1341 https://doi.org/10.1016/j.mri.2012.05.001 MRIMDQ 0730-725X
(2012).
Google Scholar
A. Fedorov et al.,
“3D Slicer as an image computing platform for the quantitative imaging network,”
Magn. Reson. Imaging, https://www.slicer.org/
(2022).
Google Scholar
R. Oostenveld et al.,
“FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data,”
Comput. Intell. Neurosci., 2011 1
–9 https://doi.org/10.1155/2011/156869
(2011).
Google Scholar
, “CGAL User and Reference Manual,”
https://doc.cgal.org/5.4/Manual/packages.html
(2022).
Google Scholar
P. Alliez et al.,
“3D mesh generation,”
https://doc.cgal.org/5.4/Manual/packages.html#PkgMesh3
(2022).
Google Scholar
H. Dehghani et al.,
“Near infrared optical tomography using NIRFAST: algorithm for numerical model and image reconstruction,”
Commun. Numer. Methods Eng., 25 711
–732 https://doi.org/10.1002/cnm.1162 CANMER 0748-8025
(2009).
Google Scholar
H. Dehghani et al.,
“Depth sensitivity and image reconstruction analysis of dense imaging arrays for mapping brain function with diffuse optical tomography,”
Appl. Opt., 48 D137
–D143 https://doi.org/10.1364/AO.48.00D137 APOPAI 0003-6935
(2009).
Google Scholar
M. Jermyn et al.,
“Fast segmentation and high-quality three-dimensional volume mesh creation from medical images for diffuse optical tomography,”
J. Biomed. Opt., 18 086007 https://doi.org/10.1117/1.JBO.18.8.086007 JBOPFO 1083-3668
(2013).
Google Scholar
J. P. Culver et al.,
“Diffuse optical tomography of cerebral blood flow, oxygenation, and metabolism in rat during focal ischemia,”
J. Cereb. Blood Flow Metab., 23 911
–924 https://doi.org/10.1097/01.WCB.0000076703.71231.BB
(2003).
Google Scholar
R. Penrose,
“A generalized inverse for matrices,”
Math. Proc. Cambridge Philos. Soc., 51 406
–413 https://doi.org/10.1017/S0305004100030401 MPCPCO 0305-0041
(1955).
Google Scholar
BiographyAlexander Ruesch received his PhD from Carnegie Mellon University (CMU), where he adapted instrumentation and developed computational models to predict intracranial pressure and cerebral autoregulation impairment in humans and non-human primates. He is a postdoctoral research fellow at CMU, where he explores the application of non-invasive, diffuse optical instruments to study diving physiology, and cognition in marine mammals and human free divers. Deepshikha Acharya acquired her PhD from CMU, where she was co-mentored by Prof. Jana M. Kainerstorfer and Prof. Matthew A. Smith. Her doctoral work was focused on finding non-invasive neurovascular biomarkers for distribution in cerebral blood flow regulation. She is a postdoctoral researcher working with Prof. Gemma Bale at the University of Cambridge interested in the effects of aging and dementia on cerebral function and metabolism. Eli Bulger received his BASc degree in nanotechnology engineering from the University of Waterloo in 2019. He is a PhD student in biomedical engineering at CMU, co-advised by Barbara Shinn-Cunningham and Jana Kainerstorfer. His research interests include methodology and applications of cognitive neuroscience, neurovascular coupling, and non-invasive source imaging. Jiaming Cao received his PhD from the Department of Biomedical Engineering at CMU, advised by Dr. Jana Kainerstorfer and Dr. Pulkit Grover. He is a research fellow at the School of Computer Science of the University of Birmingham in Dr. Hamid Dehghani’s Group. His work mainly focuses on computational modeling and improving data processing algorithms in NIRS. J. Christopher McKnight is a research fellow in the Sea Mammal Research Unit, University of St Andrews, United Kingdom. His work mainly focuses on the diving behavior and diving physiology of marine mammals and humans, and the development of animal-borne technology. Mercy Manley began in veterinary medicine as a veterinary assistant in a veterinary general practice in 2002. Earning her RVT in 2004, she continued in research at Stanford Medical School, Roche Pharmaceuticals, and MD Anderson Cancer Center. She worked in veterinary referral practice, focusing on internal medicine and oncology. She was the research coordinator and veterinary nurse for the Secret Garden and Dolphin Habitat, working with dolphins and large cats. She is currently teaching in a veterinary nursing program. Andreas Fahlman is a comparative physiologist and studies the physiological limitations for survival in breath-hold diving vertebrates. Of special focus is the cardiorespiratory coupling seen in marine mammals that allow them to dive to extreme depths for prolonged durations. Barbara G. Shinn-Cunningham is director and Cowan Professor of Auditory Neuroscience at the Neuroscience Institute of CMU. She has won mentorship awards from the Society for Neuroscience and the Acoustical Society of America (ASA), the ASA Helmholtz-Rayleigh Interdisciplinary Silver Medal, a Sloan Foundation Fellowship, and a Vannevar Bush Fellowship. She is a fellow of ASA and the American Institute for Medical and Biological Engineers. Her work explores auditory attention, spatial perception, and scene analysis. Jana M. Kainerstorfer earned her PhD in physics from the University of Vienna in Austria in partnership with the National Institutes of Health and worked as a postdoctoral fellow at Tufts University. She is an associate professor of biomedical engineering at CMU. Her lab at CMU focuses on diffuse optical imaging technology development for monitoring cerebral and tissue health, with an emphasis on clinical translation. She serves on program committees at national and international conferences, is an associate editor for the Journal of Biomedical Optics, and is a senior member of the Optical Society of America. |
Brain
Near infrared spectroscopy
Brain tissue
Muscles
Tissues
Animals
Optical properties