Functional optical topography (OT) measures the changes in oxygenated and deoxygenated hemoglobin (HbO(2), HHb) across multiple brain sites which occur in response to neuronal activation of the cerebral cortex. However, identification of areas of cortical activation is a complex task due to intrinsic physiological noise and systemic interference and careful statistical analysis is therefore required. A total of 10 young healthy adults were studied. The activation paradigm comprised of anagrams followed by finger tapping. 12 channels of the OT system were positioned over the frontal cortex and 12 channels over the motor cortex while the systemic physiology (mean blood pressure (MBP), heart rate (HR), scalp flux) was simultaneously monitored. Analysis was done using the functional Optical Signal Analysis (fOSA) software and Statistical Parametric Mapping (SPM), where we utilized two approaches: (i) using only HbO(2) as a regressor in the general linear model (GLM) and (ii) using all of the explanatory variables (HbO(2), MBP, HR and scalp flux) as regressors. Group analysis using SPM showed significant correlation in a large number of OT channels between HbO(2) and systemic regressors; however no differences in activation areas were seen between the two approaches

Elwell, CE

Koh, PH

Stubbs, C

Tachtsidis, I

English

PubMed

UCL Discovery

Functional optical topography analysis using Statistical Parametric Mapping (SPM) methodology with and without physiological confounds Ilias Tachtsidis1, Peck H. Koh1, Charlotte Stubbs1, Clare E. Elwell1 1Biomedical Optics Research Laboratory, Department of Medical Physics and Bioengineering, University College London, Gower Street, London WC1E 6BT Abstract   Functional optical topography (OT) measures the changes in oxyge nated  and  deoxygenated  hemoglobin  (HbO2,  HHb)  across  multiple  brain  sites which occur in response to neuronal activation of the cerebral cortex. However, identification of areas of cortical activation is a complex task due to intrinsic phy siological noise and systemic interference and careful statistical analysis is there fore required. A total of 10 young healthy adults were studied. The activation pa radigm comprised of anagrams followed by finger tapping. 12 channels of the OT system were positioned over the frontal cortex and 12 channels over the motor cortex  while  the  systemic  physiology  (mean  blood  pressure  (MBP),  heart  rate (HR), scalp  flux)  was  simultaneously  monitored. Analysis  was  done  using  the functional  Optical  Signal  Analysis  (fOSA)  software  and  Statistical  Parametric Mapping (SPM), where we utilized two approaches: (i) using only HbO2 as a re gressor in the general linear model (GLM) and (ii) using all of the explanatory va riables (HbO2, MBP, HR and scalp flux) as regressors. Group analysis using SPM showed significant correlation in a large number of OT channels between HbO2 and systemic regressors; however no differences in activation areas were seen be tween the two approaches. 1 Introduction Optical topography (OT) techniques to map functional brain activation rely on making  simultaneous  near infrared  spectroscopy  (NIRS)  measurements  of changes in oxygenated and deoxygenated hemoglobin (HbO2, HHb) concentration at multiple brain sites. The generated spatial maps of the hemoglobin concentra tion changes correspond to specific regions of the cerebral cortex. However, iden tification of areas of cortical activation require statistical analysis with researchers using classical statistical approaches such as “Student’s t test” and/or more com plex tools such as Principal Component Analysis, Independent Component Analy sis, and more recently Statistical Parametric Mapping (SPM) [1 5]. Interpretation of the optical functional activation remains a complex task due to the intrinsic physiological noise and systemic interference [6 8]. The aim of this study is to in vestigate the use of SPM analysis for interpretation of functional optical topogra phy hemoglobin changes during frontal lobe anagram activation and motor cortex finger tapping activation with and without accounting for physiological and sys temic changes. 2 Methods We studied 10 subjects (5 male and 5 female), all had English as their first lan guage and were right handed (this study was approved by the UCL Research Eth ics Committee). The age range was between 19 and 27 years with a mean age of 22 years.  NIRS measurements were conducted with the ETG 100 Optical Topography System (Hitachi Medical Co., Japan) using two 12 channel arrays. Each optode ar ray consisted of 5 source optodes (each delivering light at 780 and 830 nm) and 4 detector  optodes.  The  source detector interoptode spacing was  30mm and data were acquired at 10Hz. The optodes were placed over the subject’s left frontal cor tex and positioned according to the International 10 20 system of electrode place ment such that channels 1 12 were centered approximately over the frontopolar region (Fp) and channels 13 24 were centered approximately over the left primary motor cortex (C3) (see Fig. 1). The Cartesian coordinates of the NIRS optodes placement were measured   in relation to known fiducial landmarks (nasion, inion, left and right ear and top of the head), using the Polhemus tracking device (Isotrak II, USA) and the inter subject variability were calculated. For this study the x, y and z coordinates of each of the optodes did not vary between subjects by more than ±2cm of the mean position of each probe. A Portapres® system (TNO Institute of Applied Physics) was used to conti nuously and non invasively measure mean blood pressure (MBP) and heart rate (HR) from the finger. A laser Doppler probe (FloLab, Moore Instruments) was placed over the forehead to monitor the changes in scalp blood flow (flux). Data were recorded during 1 minute of the subject at rest (baseline), followed by 45 seconds of the subject solving anagrams (a mixture of 4 letter, 6 letter and 8 letter anagrams), followed by 30 seconds of rest, followed by 30 seconds of fin ger tapping. In this study solving an anagram was defined as producing one cohe rent word using only the letters from another word (e.g. icon – coin). Each ana gram solving and finger tapping period was repeated a total of eight times, with the study ending after a 60 second rest period (total study time 19.5 minutes). Dur 3 ing this functional paradigm the participants were asked to be still and silent; and no response to anagram solving was requested. All data were subjected to an identical preprocessing procedure using the func tional Optical Signal Analysis program (fOSA) [5] (University College London, UK) to convert the relative changes in light intensities to concentration changes in HbO2 and HHb using a differential pathlength factor correction of 6.26. All the signals including MBP, HR and flux, were linear detrended, then decimated from 10Hz to 1Hz and smoothed using a 3 second moving average. Before performing the statistical analysis using the SPM, image processing of the HbO2 data was performed using fOSA as described by Koh et al. [5]. In brief, it involves the spatial conversion of the OT maps of each measurement channel to its corresponding  pixel location for each of the 24 channels. A  series  of two dimensional SPM OT maps were generated with each SPM OT map representing a 1 second time point. To improve on the sensitivity and facilitate inter subject comparisons, an interpolation process using nearest neighbors is used to produce a more realistic data grid.  In our SPM analysis we utilized two approaches to model the OT data; in the first approach we used only the HbO2 signal in the General Linear Model (GLM) and in the second approach we incorporated the physiological measures (MBP, HR and flux) in addition to the HbO2 signal as additional regressors in our statis tical model. To treat the variability of the haemodynamic response (HbO2) arising from different events between different brain pixels, we modeled both latency and dispersion  derivatives as  additional  regressors  to  the  canonical  Haemodynamic Response Function (HRF) (see Fig. 2) in the two approaches. The design matrix is shown in Figure 3. The resulting statistical maps were family wise error corrected to account for the spatial correlation between neighboring pixels in OT measurements. An F contrast was defined in the SPM test for each subject to infer on any interaction (canonical HRF, its temporal and dispersion derivatives) in each of the two tasks. A group t statistic was then conducted in SPM to infer on its population level. 3 Results Figure 4 shows the SPM group t results, without physiological confounds, for (a) the anagram task and (b) the finger tapping for the HbO2 signal. Figure 5 shows the SPM group t results, with physiological confounds, for (a) the anagram and (b) the finger tapping task for the HbO2 signal. Darker pixels correspond to higher t values, with positive t scores for HbO2 signal indicating a significant task related increase in HbO2. Figure 5(c) shows the SPM group t results for the physiological confounds; in this map darker pixels correspond to significant correlations of one or a combination of the regressors MBP, HR and flux, with the HbO2 haemody namic response.   Results from the group level analysis in the current study have demonstrated that the most significant channels (dark regions) during anagram activation are be tween 19, 21, 22 and 24 which are on the lower site of the temporal lobe; and dur ing finger tapping are channels 2 and 18 which are on the upper part of the tem poral lobe. 4 Discussion There were no differences in the activation responses in the SPM group analysis with and without physiological confounds in this functional paradigm. Group level analysis on the systemic regressors shows highly significant correlation of the sys temic variables in most channels in the study; however this did not make a differ ence in the interpretation of the task related HbO2 response.  Interpretation of the optical data during functional activation remains a com plex task as there are many confounds due to different degrees of physiological noise and systemic interference both at the intra and inter subject level. A unified approach of statistical investigation for inference analysis of functional activation with optical topography is needed. We have previously compared “Student’s t test” with SPM methods and found the SPM to be a better analysis tool [8]. Here we detail the SPM method for OT data. However further work is needed to refine the Optical SPM analysis, especially for the development of an optical response function similar to the haemodynamic response function (used in this study). Acknowledgments    The  authors  would  like  to  acknowledge  the  EPSRC  (Grant  No EP/D060982/1). References 1. Hirosaka R, Katura T, Kawaguchi H et al (2004) Noisy time delayed decorrelation and its ap plication to extraction of neural activity from single optical recordings in guinea pigs. Physica D 320 332 2. Yoshida T, Sakagami M, Yamazaki K et al (2007) Extraction of neural activity from in vivo optical  recordings  using  multiple  independent  component  analysis.  IEEJ  Trans  EIS 127(10):1642 1650 3. Zhang X, Toronov V, Webb A (2005) Simultaneous integrated diffuse optical tomography and functional magnetic resonance imaging of the human brain. Opt Express 13(14): 5513 5521 4. Schiessl I, Stetter M, Mayhew JEW et al (2000) Blind signal separation from optical imaging recordings with extended spatial decorrelation. IEEE Trans on Biomed Eng 47(5):573 577 5. Koh PH, Glaser DE, Flandin G et al (2007) Functional optical signal analysis (fOSA): a soft ware tool for NIRS data processing incorporating statistical parametric mapping (SPM). JBO 12(6) 5 6. Tachtsidis I, Leung TS, Devoto L et al (2008) Measurement of frontal lobe functional activa tion and related systemic effects: a near infrared spectroscopy investigation. Adv Exp Med Biol 614:397 403 7. Tachtsidis I, Leung TS, Tisdall MM et al (2008) Investigation of frontal cortex, motor cortex and systemic haemodynamic changes during anagram solving. Adv Exp Med Biol 614: 21 28 8. Tachtsidis I, Leung TS, Chopra A et al (2009) False positives in functional near infrared topo graphy. Adv Exp Med Biol In Press 9. Friston KJ, Holmes AP, Poline JB et al (1995) Analysis of fMRI time series revisited. Neu roImage 2(1):45 53  Fig. 1 A picture illustrating the approximate positions of the OT light sources, detectors and locations of corres ponding measuring positions/channels co registered with brain  MRI  using  the  Pholemus  information  (numbers represent OT channels).          Fig. 2 Canonical HRF (2 gamma functions) plus Multiva riate Taylor expansion in time (Temporal Derivative) and width (Dispersion Derivative) [9].        Fig. 3 Design matrix for the oxyhaemoglobin model for the anagrams (first three columns), finger tapping (next three columns) tasks followed by the three re gressors (physiological confounds) and a constant er ror term for the design matrix. Each row represents a time slice of the topographic brain map.  Fig. 4 The SPM group t results without physiological confounds with a threshold of t=4.29 and significance of p≤0.05; darker pixels correspond to higher significant t values.  Fig. 5 The  SPM group t results with physiological  confounds  and  a  thre shold  of  t=4.29  and  significance  of p≤0.05;  darker  pixels  correspond  to higher significant t values.  

Analysis of fMRI time series revisited.

Blind signal separation from optical imaging recordings with extended spatial decorrelation.

Chopra A et al (2009) False positives in functional near infrared topo graphy. Adv Exp Med Biol In

Extraction of neural activity from in vivo optical recordings using multiple independent component analysis.

Functional optical signal analysis (fOSA): a soft ware tool for NIRS data processing incorporating statistical parametric mapping (SPM).

Measurement of frontal lobe functional activa tion and related systemic effects: a near infrared spectroscopy investigation.

Noisy time delayed decorrelation and its ap plication to extraction of neural activity from single optical recordings in guinea pigs.

Tisdall MM et al (2008) Investigation of frontal cortex, motor cortex and systemic haemodynamic changes during anagram solving. Adv Exp Med Biol 614:

Functional Optical Topography Analysis Using Statistical Parametric Mapping (SPM) Methodology with and without Physiological Confounds

Functional Optical Topography Analysis Using Statistical Parametric Mapping (SPM) Methodology with and without Physiological Confounds

Abstract

Similar works

Full text

Available Versions

UCL Discovery