Using a multimodal near-infrared spectroscopy and MRI to quantify gray matter metabolic rate for oxygen: A hypothermia validation study

Non-invasive quantitative imaging of cerebral oxygen metabolism (CMRO2) in small animal models is crucial to understand the role of oxidative metabolism in healthy and diseased brains. In this study, we developed a multimodal method combining near-infrared spectroscopy (NIRS) and MRI to non-invasively study oxygen delivery and consumption in the cortex of mouse and rat models. The term CASNIRS is proposed to the technique that measures CMRO2 with ASL and NIRS. To determine the reliability of this method, CMRO2 values were compared with reported values measured with other techniques. Also, the sensitivity of the CASNIRS technique to detect changes in CMRO2 in the cortex of the animals was assessed by applying a reduction in core temperature, which is known to reduce CMRO2. Cerebral blood flow (CBF) and CMRO2 were measured in five mice and five rats at a core temperature of 37 °C followed by another measurement at 33 °C. CMRO2 was 7.8 ± 1.8 and 3.7 ± 0.9 (ml/100 g/min, mean ± SD) in mice and rats respectively. These values are in good agreement with reported values measured by 15O PET, 17O NMR, and BOLD fMRI. In hypothermia, we detected a significant decrease of 37% and 32% in CMRO2 in the cortex of mice and rats, respectively. Q10 was calculated to be 3.2 in mice and 2.7 in rats. In this study we showed that it is possible to assess absolute values of metabolic correlates such as CMRO2, CBF and oxygen extraction fraction (OEF) noninvasively in living brain of mice and rats by combining NIRS with MRI. This will open new possibilities for studying brain metabolism in patients as well as the many mouse/rat models of brain disorders.


Introduction
The brain is highly aerobic, using approximately 20% of the body's total oxygen uptake. The gray matter (GM) consumes 3-4 times more oxygen than the white matter (Dunning and Wolff, 1937;Pantano et al., 1984) and has a higher mitochondrial density (Gerstl et al., 1969;Santuy et al., 2018). This requirement for oxygen means that changes in the rate of oxygen uptake, i.e., the cerebral metabolic rate of oxygen (CMRO 2 ), should be a sensitive index of injury and disease. In addition, CMRO 2 is usually tightly coupled to both metabolic rate in the brain and to cerebral blood flow (CBF). A method of monitoring all parameters simultaneously would be of great benefit for studying the regulation of oxygen delivery as well as for understanding how oxygen uptake (related to mitochondrial function) is impacted by disease.
In most neurological diseases, including multiple sclerosis (Lin and Beal, 2006;Sun et al., 1998), Alzheimer's Disease (Ni et al., 2018) and Parkinson's disease (Borghammer et al., 2012), abnormal oxidative metabolism is likely to occur. Consequently, both CMRO 2 and CBF may be affected. These physiological parameters can be detected currently by neuroimaging modalities such as positron emission tomography (PET) (Mintun et al., 1984;Ter-Pogossian et al., 1970), functional magnetic resonance imaging (fMRI) (Hyder et al., 2001;Sicard and Duong, 2005), and 17 O-NMR (Zhu et al., 2005). However, these techniques involve the use of either O 2 isotopes as exogenous tracers, or complex biophysical models. In addition, technical procedures, such as intubation, and placing arterial and venous catheters, are required in these techniques and present more challenges in small animals (Cui et al., 2013). Since many of the neurological disease models are mouse or rat models, a non-invasive and relatively simple means of accurately measuring oxidative metabolism and hemodynamics in these models, would provide an important tool for exploring cerebral pathophysiology. In this paper, we propose a novel method for quantification of CMRO 2 , that is tailored specifically to rodent models, but also has potential for use in patients.
Based on the Fick principle, CMRO 2 can be derived from the product of CBF, hemoglobin concentration in blood, and the oxygen extraction fraction (OEF) which is equal to the change in oxygen content of blood between the arterial and venous systems. To obtain these values, we need to quantify CBF, arterial oxy-hemoglobin saturation (S a O 2 ), tissue microvascular oxyhemoglobin saturation (S t O 2 ) and total Hb (tHb) concentration in the brain. In our multimodal system, CBF is quantified using arterial spin labeling (ASL) MRI. S a O 2 is quantified with a pulse oximeter. Using broadband near-infrared spectroscopy (NIRS), we quantify tHb and S t O 2 , which can then be used to estimate venous oxygen content (see section 2.5.2).
NIRS is widely used in experimental and clinical applications to measure changes in oxygenation and hemodynamics. We showed that with a custom-built broadband NIRS system and an adapted processing pipeline based on the second differential method, it is possible to quantify the absolute concentration of dHb at any given time point in the cortex of rat models (Zhang et al., 2010). By applying an anoxia pulse method (Cooper et al., 1998), we can quantify tHb concentration and S t O 2 . Due to the "banana" shaped path of light in the brain (Chuang et al., 2013;Johnson, 2017;Mansouri et al., 2010), we will be measuring largely from the animal GM.
The objective of the current study is to validate our proposed CMRO 2 measurement technique by manipulating metabolic rate in commonly used rodent species (mouse and rat). For this, we will be applying mild hypothermia, as it is known to reduce brain CMRO 2 (Klementavicius et al., 1996;Zhu et al., 2005Zhu et al., , 2007. Additionally, it has been shown that hypothermia, and also isoflurane anesthesia, induce burst suppression (BS) in the cortex, which is an EEG state associated with brain inactivation (Ching et al., 2012). This EEG pattern was found to have neuroprotective effects, achieved through the reduction of brain metabolism (Doyle and Matta, 1999). Therefore, our experimental model was to reversibly change the core temperature of the animals from 37 C to 33 C. We show that there is reasonable agreement with literature values (when measured at 37 C), and that we can reproducibly detect a decline in CMRO 2 with hypothermia.
We propose the term CASNIRS to represent our technique, measuring CMRO 2 using ASL and NIRS. With CASNIRS, we show that one can quantify tissue oxygenation, CBF, OEF and CMRO 2 in the GM of rodent models of neurological disorders. This method could be translated to patients as well.

Study design
The animals were anesthetized with 5% isoflurane added to a gas mixture of 70% N 2 and 30% O 2 and maintained at 2% isoflurane. Hair from the mice and rats head was removed using depilatory cream (Nair; Church & Dwight Co. Inc., USA) to reduce the scattering effect and improve the coupling between the optic fibers and the head. The fibers were placed on top of the head, with the posterior edge of the prism placed on the line between the two external auditory meatus (within 1 mm of the interaural line). The fibers were secured by covering the area with masking tape connected to the outside of the MRI coil cradle. Fibers were spaced with a black rubber separator to keep the optimal distance constant between the source and detector fibers (4 mm for mouse and 7 mm for rat) and to help reduce direct illumination between the source and detector. These separation distances, which were calculated by simulations for light propagation modelling (unpublished data), ensure that the optical signal emanates mainly from the GM of the cerebral cortex (Johnson, 2017). Heart rate, breathing rate and S a O 2 were monitored by the MouseOx MRI-compatible pulse oximeter (Starr Life Sciences, USA) on the shaved left thigh of the animal. Core temperature was monitored by a rectal thermometer. The animal was placed on a heating pad in the MRI coil together with the aforementioned instruments and positioned for optimal image quality. Animals were monitored for 30 min to ensure a stable core temperature of 37 C was reached before data acquisition. After 30 min, ASL acquisition was initiated and required approximately 14 min. Core temperature was then reduced to 33 C (over approximately 8 min). When the temperature was stable (8 min later), ASL was repeated. Core temperature was returned to 37 C (approximately 5 min to return core temperature and 8 min to allow the physiology to stabilize). An anoxia pulse was given for 50s followed by 5 min of 70%N 2 /30%O 2 . The animal was removed from the MRI and a blood sample (150 μl) was taken from the tail and analysed with the blood gas analyser (Stat Profile pHOx Ultra, Nova Biomedical Corporation, USA) to quantify the hemoglobin concentration in large vessels. NIRS data were averaged over the time of the ASL for use in the CMRO 2 calculations.

Animals
5 Male Wistar rats (150-200 g, 6 weeks old from Charles River, Qu ebec, Canada), and 5 male C57BL/6J mice (21-23 g, 5-6 weeks old from Charles River, Qu ebec, Canada) were housed and maintained in the University of Calgary Animal Care facility with a 12 h light and dark cycle with access to water and food pellets ad libitum. Animal protocols were approved by the Animal Care Committee of the University of Calgary and conformed to the guidelines established by the Canadian Council of Animal Care (CCAC).

NIRS system
A custom-built continuous wave (CW) NIRS system, consisting of a broadband white light source and a spectrograph, was employed to measure the attenuation spectrum (255 wavelengths) over the range of 705-960 nm. A single MR-compatible optic fiber with a 1000 μm core diameter and 0.39 NA (FT1000EMT, Thorlabs Inc, USA) was used to deliver the light provided by the 100-W broadband quartz-tungstenhalogen light bulb of the fiber optic illuminator (model 77501, Oriel Instruments Inc., USA) to the animal brain. A focal lens (AC127-019-B-ML, Thorlabs Inc, USA) was used to focus the output light onto the fiber optic. The light intensity was moderated and optimized by an iris and a shutter to ensure that the power deposition did not damage the biological tissue. A GRIN lens with a 1.80 mm diameter and 0.55 NA (#64-525, Edmund Optics, USA) and a 90 prism with a leg length of 2 mm (#45-524,TECHSPEC N-SF11, Edmund Optics, USA) were glued with a cyanoacrylate containing adhesive (Krazy Glue, Elmer's Products, Inc, USA) to the end of the fiber to collimate the light and direct it into the tissue (Johnson, 2017). A schematic of the experimental setup is shown in Fig. 1. Light coupling to the skin was enhanced by covering the prisms with a thin layer of glycerol (Bashkatov et al., 2002). A similar optic fiber was used to transport the NIR light from the mouse brain to the spectrograph (Shamrock 303i, Andor Technology Inc, Northern Ireland). Input NIR light was spread along wavelengths, using a diffraction grating, to form the NIR spectrum. Collected spectra were projected onto the CCD camera (Andor iDus) mounted on the spectrograph for digitization. The CCD camera was cooled to -40 C to significantly reduce the dark noise contribution.

Magnetic resonance imaging (MRI)
A 9.4 T horizontal bore MRI (Bruker Avance console, Bruker Biospin GmbH, Rheinstetten, Germany) with a 35 mm quadrature volume coil was used to non-invasively quantify absolute CBF with an ASL sequence. An axial slice was acquired around the bregma, where the optic fibers were located, using a CASL-HASTE sequence (Lei et al., 2001) with the following parameters: TR ¼ 3000 ms, TE ¼ 13.3 ms, FOV ¼ 30 Â 30 mm, matrix size ¼ 128 Â 128 pixels, slice thickness ¼ 1 mm, 16 averages. Four perfusion images were collected per measurement: 2 control images and 2 tagged images, to correct for magnetization transfer. Following these images, a T 1 map was obtained in the same location using a RARE-VTR sequence where TE ¼ 10 ms, TR ¼ 100, 500, 1000, 3000, and 7500 ms. Together, the four perfusion images and the T 1 map were collected over a period of 14 min.
2.5. Data analysis 2.5.1. CBF quantification A perfusion map was generated from the four perfusion images and the T 1 map, using a custom Matlab script (Johnson, 2017). In this map, CBF was calculated on a voxel-by-voxel basis using the following equation (Buxton, 2005;Pekar et al., 1996): where ƒ is the tissue perfusion in ml/100 g/min, λ ¼ 0:9 is the bloodbrain partition coefficient (Herscovitch and Raichle, 1985;Leithner et al., 2010), M c is the average signal of the two control images, M L is the average signal of the two tagged images, α ¼ 0:675 AE 0:044 is the spin-labelling efficiency (Johnson, 2017), and T ' 1 is the measured T 1 value in each voxel. Individual perfusion maps were calculated for each animal. An ROI was created using a tagged image, which has more contrast, and which encompassed a region that is estimated to be similar to the region of sensitivity for NIRS. This is largely in cortical GM (Fig. 2). This ROI was then used on the perfusion and T 1 maps to obtain a mean AE SD (ml/100 g/min) for perfusion and T 1 (ms) values for each condition.

dHb and S t O 2 quantification
A Matlab algorithm (Zhang et al., 2010) based on the modified Beer-Lambert law was used to convert the attenuation light, measured by NIRS, into concentration of the chromophore dHb. The second differential spectral analysis method (Matcher et al., 1994) was employed in the algorithm to remove the scattering effect from the measured absorbance spectra. The well-defined spectral absorption feature of water at 800-850 nm and the cerebral water content, which is assumed to be 80% in adult rodents (Reinoso et al., 1997), were used to estimate the pathlength that light travels through the tissue. The wavelengths analysed for dHb were 720-810 nm. The algorithm applied multilinear regression to fit the measured attenuation spectra to the specific extinction coefficients of the pure chromophores, downloaded from the University College of London medical physics website (UCL, 2005). The concentration of tHb in the cerebral tissue was determined using the anoxia pulse method (Cooper et al., 1998;Zhang et al., 2010), in which it is assumed that during an anoxic pulse the total hemoglobin is approximated by the total [dHb] ([tHb] ¼ [dHb]). The anoxic pulse was obtained by switching the inspired gas to 0% O 2 and 100% N 2 for 50 s. This period was sufficient to reach steady state of dHb while minimizing the effect on other physiological parameters (Cooper et al., 1998;Zhang et al., 2010). Calculated dHb and tHb concentrations were used to determine the oxygen saturation in tissue (S t O 2 ) with the following equation:

CMRO 2 quantification
Based on the modified Fick principle (Tichauer et al., 2006b), where k ¼ 1.39 (ml(O 2 )/g(Hb)) is a factor describing the amount of O 2 bound to Hb when completely saturated (Brown et al., 2003), and [Hb] is the concentration of hemoglobin (g/l) in large vessels. Since OEF is equal to the arteriovenous oxygen difference, Eq.
(3) can be rewritten as the following:   (Phelps et al., 1979;Tichauer et al., 2006b), an estimate of S v O 2 can be expressed as: By substituting Eq. (5) in Eq. (4), CMRO 2 can be quantified by combining CBF values with NIRS and pulse oximeter data, using the following equation (Eq. (6)): CBF is the average value of cerebral blood flow (ml/100 g/min) calculated from the ROI in the cortex (Fig. 2). S a O 2 is the arterial blood oxygen saturation, measured by the MouseOx MRI-compatible pulse oximeter on the thigh of the animal, and S t O 2 is the tissue oxygen saturation calculated from [dHb] and [tHb] obtained with the NIRS system (Eq. (2)).
[Hb] is the concentration of hemoglobin (g/l) in large vessels, measured from the blood sample with the blood gas analyser. R is the ratio of small to large vessel hematocrit set at 0.61 (Bereczki et al., 1993).

Statistical analysis
All the measured parameters were compared during normothermia and hypothermia using a paired t-test, where p < 0.05 was used as the threshold for significance.

Data availability
The raw data used to calculate CMRO 2 are available in PRISM Dataverse: University of Calgary's Data Repository. More data may be available upon direct request.

Results
A representative mouse perfusion map is shown in Fig. 2. Perfusion images and a T 1 map were acquired first when the core temperature of the animal was stable at 37 C, and again after stabilizing the animal temperature at 33 C. Perfusion maps showed bilateral reduction in CBF in the cortex of the animals under hypothermic conditions (Fig. 2,  Table 1). A decrease in T 1 values was observed as well (Table 1). Mean values of T 1 and CBF were quantified from the same ROI and are presented per each mouse (Fig. 3A) and rat (Fig. 3B) during normothermia and hypothermia.
To calculate CMRO 2 , CBF values were combined with NIRS and pulse oximeter data, using equation (6), which is based on the modified Fick principle (Tichauer et al., 2006b). Fig. 4A and B shows representative data extracted from the recorded NIR spectra over the course of 90 min. The concentration of dHb (μM) (Fig. 4A) was quantified from the raw NIR spectra, using the developed processing algorithm based on the second-differential method. Blue arrows in Fig. 4A indicate the initiation of ASL measurements, anoxic pulse, or changing core temperature. Fig. 4B shows the S t O 2 calculated using Equation (2). Fig. 4C displays a time-series of the arterial oxygen saturation (S a O 2 ) measured by the MouseOx pulse oximeter from the left thigh of the animal concurrently to the NIRS data collection. Upon applying the anoxic pulse, indicated by the 5th arrow (Fig. 4A), the concentration of dHb started to increase sharply and remained steady when it reached the tHb concentration, while the tissue oxygen (Fig. 4B) and arterial oxygen (Fig. 4C) saturation dropped immediately. Following the cessation of the anoxic pulse, which lasted 50 s, the values of dHb, S t O 2, and S a O 2 recovered to their baseline values. Fig. 5 shows S t O 2 and S a O 2 in each of the 10 animals. S t O 2 increased significantly (p < 0.001) in mice under hypothermic conditions and trended toward an increase in rats (p ¼ 0.1), while S a O 2 decreased in mice (p ¼ 0.04) and no significant change was found in rat models. OEF and CMRO 2 were quantified per each mouse (Fig. 6A) and each rat (Fig. 6B) during normothermia (37 C) and mild hypothermia (33 C). Table 1 summarizes the average values of these parameters in each Table 1 Quantification of T 1 , CBF, S t O 2, S a O 2 , OEF and CMRO 2 in mouse (n ¼ 5) and rat (n ¼ 5) models during normothermia (37 C) and hypothermia (33 C). (mean AE SD; paired t-test).  condition. In mouse models, a significant drop of 14% in CBF (p ¼ 0.01), 24% in OEF (p < 0.01) and 37% in CMRO 2 (p < 0.01) were observed. In rat models, CMRO 2 declined significantly by 32% (p < 0.01), although the drop in CBF and OEF was not significant.

Discussion
High-field MRI was combined with broadband NIRS and a pulse oximeter to simultaneously measure CBF, OEF, and CMRO 2 in the cortex of mouse and rat models. The term CASNIRS is proposed here to represent Fig. 3. Quantification of T 1 and CBF in each of the 5 mice (A) and rats (B) during normothermia (37 C) and hypothermia (33 C). Each symbol represents a different animal. Statistical analysis was performed comparing normothermia and hypothermia (** -p 0.01) using paired t-test (Table 1). Fig. 4. Representative data extracted from the NIR spectra and the MouseOx pulse oximeter collected from one mouse: A) The concentration of dHb (μM) quantified from the raw NIR spectra, B) The tissue O 2 saturation, S t O 2 (%), calculated with the dHb and tHb from (A), C) The arterial O 2 saturation, S a O 2 (%), measured by the MouseOx pulse oximeter on the thigh. Each arrow indicates a different challenge or measurement: 1) First ASL measurement at 37 C and 30% O 2 , 2) Core temperature is stable at 33 C, 3) Second ASL measurement at 33 C and 30% O 2 , 4) Core temperature is stable at 37 C, 5) Initiation of an anoxic pulse which lasts for 50 s. the technique measuring CMRO 2 with ASL and NIRS. One type of validation for the CASNIRS method is to compare the absolute values of our data to reported values for rodent CMRO 2 .
The most common methods currently available to determine CBF and CMRO 2 include 15 O PET (Frackowiak et al., 1980) and 17 O MR (Cui et al., 2013;Zhu et al., 2005). These methods involve inhalation of oxygen isotopes and tracking their consumption through their oxidative pathways or determining the washout rate of 17 O-labeled water tracer introduced by a bolus injection. 15 O PET requires a radioactive and unstable oxygen isotope with a half-life of about 2 min (Frackowiak et al., 1980). 17 O MR provides a direct measure for regional CMRO 2 in vivo, but it is still challenging to perform in small animals, like mouse models, due to the requirement for intubation. Although there are no direct measures for CMRO 2 conducted by NIRS techniques, some NIRS approaches have been used to indirectly determine CMRO 2 in piglets and humans. Previous studies (Tichauer et al., 2006a;Verdecchia et al., 2013) have measured CMRO 2 in newborn piglets by combining two near-infrared techniques: diffuse correlation spectroscopy (DCS) and time resolved (TR) NIRS. However, CBF was determined in an invasive way which involved the injection of indocyanine green (ICG) as an intravascular tracer. Other studies (Elwell et al., 2005;Yoxall and Weindling, 1998) were conducted on full and pre-term infants, using spatially resolved spectroscopy to measure absolute mean cerebral oxygen saturation, and CBF by delivering a bolus of HbO 2 to the brain, and then deriving CMRO 2 values from these measurements. To our knowledge, no similar studies have been conducted on small animals like rats and mice. Table 2 summarizes some values of CMRO 2 from the literature, quantified by currently existing methods in mouse and rat models under different anesthesia conditions. The CMRO 2 values for mice in Table 2 are fairly similar, with a mean AE SD of 6.4 AE 0.7 ml/100 g/min. The rat data ranged from 3.6 AE 0.1 to 7.5 AE 0.4 ml/100 g/min. These values are broadly consistent with the normothermic values measured with the multimodal CASNIRS technique (7.8 AE 1.8 and 3.7 AE 0.9 ml/100 g/min for mice and rats, respectively), since they fall in the same range.
Several factors can contribute to the differences between our measured values and the ones reported in literature. This wide variability can be due to the different types of anesthesia used in different studies, as anesthesia is known to play a crucial role in altering CBF values (Nakao et al., 2001). Another reason for this difference could be the different values of core temperature in different studies. Moreover, the values presented in this study are averaged over one cortical region only. However, this limitation can be overcome by applying multifiber imaging that covers the whole brain, as used in human fNIRS, instead of the two-fibers system described here.
In order to provide another level of validation for the CASNIRS method, a mild hypothermic challenge was applied. Mild hypothermia is known to reduce GM CMRO 2 (Klementavicius et al., 1996;Zhu et al., 2005Zhu et al., , 2007. Using CASNIRS, we detected a 37% and 32% decline in CMRO 2 in the brain of mice and rats, respectively, when exposed to hypothermia. Using a literature value of 4.5 ml O 2 /100 g/min at 37 C and 2.6 ml O 2 /100 g/min at 32 C, obtained in a rat model by 17 O MRSI (Zhu et al., 2007), there was a decrease of 42% in CMRO 2 when the core temperature dropped by 5 C. This extrapolates to a 33.7% decline in CMRO 2 with a 4 C (37-33 C) decline in core temperate, which is a decrease in CMRO 2 that is comparable to that found in our study.
An average of 44% decrease in CMRO 2 has been shown previously in rat brains using direct 17 O NMR detection (Zhu et al., 2003(Zhu et al., , 2007. Another study has found a reduction in CMRO 2 of 15-20% in mild hypothermia (Klementavicius et al., 1996). CMRO 2 decreased by more than 50% when mild hypothermia (37 C-34 C) was applied to control elevated intracranial pressure in patients with head injury (Shiozaki et al., 1993). A 40% decrease in CMRO 2 was reported as well in monkeys under pentobarbital anesthesia during hypothermia (33 C) (Bering, 1961).
Other studies have found a linear relationship between brain temperature and cerebral metabolism, reporting a change of 6-8% in metabolism per degree C temperature (Lanier, 1995). This corroborates with our data, since a change of 4 C would change the metabolism by 24-32%.
Finally, Q10 which is an estimate for the expected change in metabolism for a 10 C change in temperature is calculated here. According to the formula: Q10 ¼ CMRO2 CMRO 0 2 ! 10 C=ðTÀT0Þ (Michenfelder and Milde, 1991), and to the values measured by CASNIRS, Q10 is 3.2 for mice and 2.7 for rats. These two Q10 values fall within the range of the reported values (2.0-4.4) (Klementavicius et al., 1996), and contribute to the validation of our technique. This wide variability in the reported values may be explained by the fact that total CMRO 2 is comprised of basal and functional CMRO 2 .
The hypothermia intervention showed that a change in CMRO 2 is detected with CASNIRS. The decline in CMRO 2 was very close to published data. These data support the conclusion that this method has merit for quantifying CMRO 2 in rodent models.
When a decrease in CMRO 2 occurs, the hemodynamic response could be to reduce either CBF, OEF, or both. It appears that there is inter-subject variability in this response. The underlying reason is not understood, but different patterns of burst suppression induced by hypothermia may be a contributing factor. Figs. 3 and 6 demonstrate that some animals show different patterns of response, especially in the rat data. While it is too early to say whether this is an interspecies difference, there is a definite need for looking into the change in CMRO 2 as a more robust marker, rather than the change in either CBF or OEF. It has been previously suggested that quantification of CMRO 2 provides unique and additional data over quantification of either CBF or OEF alone (Tichauer et al., 2006b). Using NIRS in a pig hypoxia model, Tichauer et al. also showed that individual subjects showed variable responses in terms of OEF and CBF. This supports the conclusion that such variability occurs across species and is also likely to occur in humans.
An additional parameter that this technique can measure is S t O 2 in regional micro-vessels. This parameter is crucial as it represents the tissue Fig. 6. Quantification of OEF and CMRO 2 in each of the 5 mice (A) and rats (B) during normothermia (37 C) and hypothermia (33 C). Each symbol represents a different animal. Statistical analysis was performed comparing normothermia and hypothermia (* -p 0.05, ** -p 0.01) using paired t-test (Table 1).  (2006) oxygenation status. During normothermia, the S t O 2 value in mice was 0.72 and 0.71 in rats. In humans, S t O 2 was found to be 0.63 (Yang and Dunn, 2015). Any significant changes from these values would help determine if the tissue is hypoxic or hyperoxic. In the mild hypothermia condition, there was a significant increase (5.5%, p < 0.001) in the S t O 2 of mice while in rats the value trended toward an increase (2.9%, p ¼ 0.1). This means that mild hypothermia didn't create a hypoxic environment. In fact, there was a trend to hyperoxygenation indicating that the oxygen demand was reduced relative to the change reduction in oxygen supply (32-37% reduction in CMRO 2 vs. 14-20% drop in CBF) so oxygen which was not consumed, remained in the microvessels. This is supported by the significant decline in OEF in mice. In rats, the variability of the change in OEF resulted in no significant difference in groups, but 2 of 5 subjects showed a marked decline in OEF with hypothermia. A decrease in T 1 values was also observed as the temperature was dropped. This may be expected as the T 1 in ex vivo brain also declines with a reduction in temperature (Birkl et al., 2013).
The CASNIRS method has significant advantages for measuring oxygenation and CMRO 2 . It is non-invasive, allowing for repeated measures, does not use isotopes or ionizing radiation, and can be undertaken in any current MRI system on animals or humans.

Conclusions
The study was designed to determine the feasibility of a novel NIR-MRI technique to simultaneously assess absolute values of metabolic correlates including CMRO 2 , CBF, tHb, S t O 2 and OEF noninvasively in living mouse and rat cortex.
Using the CASNIRS multimodal technique we were able to noninvasively quantify absolute values and detect changes in CMRO 2 in the cerebral cortex of mouse and rat models under normothermic and hypothermic conditions, with a Q10 of 3.2 for mice and 2.7 for rats.
This system provides a new tool to study cerebral physiology and pathophysiology. It is relatively simple to implement within the MRI and could be useful for future research focusing on changes in oxidative metabolism in neurological diseases. This unique technique will open new possibilities for studying brain metabolism using the many mouse and rat models of brain disease available.