Introduction

In the past two decades, functional brain imaging has become a major tool to investigate the neural bases of cognition. The field of cognitive neuroscience has witnessed rapid growth following the introduction of water-labelled positron emission tomography (15O-PET) and, more recently, functional magnetic resonance imaging (fMRI). Both techniques have major advantages, among them providing a true measure of brain metabolism (PET) and remarkably good spatial and temporal resolutions (fMRI). Furthermore, these techniques allow multiple repeats of various conditions, a fact that increases the accuracy of the data collected.

However, despite these major advantages, both PET and fMRI require the cognitive task to be performed inside the scanner. The several drawbacks of this situation (uncomfortable supine position, noisy environment, limited “ecological” assessment), may affect the participant’s behaviour by causing interference with the ongoing cognitive activity [1]. Additionally, participants need to remain in the scanner throughout the whole experiment, which may be hard for individuals who, for whatever reason, have difficulty staying in a motion-free decubitus position for a long time.

By contrast, single photon emission computerized tomography (SPECT) using 99mTc-HMPAO (d,l-hexamethyl propyleneamine oxime) or 99mTc-ECD (L,L-ethylcysteinate dimer) allows dissociation of task administration and image acquisition. Cognitive activity occurs outside the scanner, and the time for which participants are required to remain in the scanner is substantially reduced, corresponding only to the duration of image acquisition. Hence, contrary to PET or fMRI, SPECT provides a means to measure brain regional metabolism associated with cognitive processes while keeping conditions close to standard laboratory ones.

Only few studies have taken advantage of SPECT for its use with healthy young adults [2-10], although they have provided encouraging support for the feasibility of activation protocols with 99mTc-ECD or 133Xe. To our knowledge, however, only one cognitive activation study using 99mTc-HMPAO has been reported with healthy young subjects [11]. Hence, there is a need to validate the use of this particular radiotracer in experiments devoted to investigate brain-cognition relations. Moreover, and given that HMPAO uptake is closely related to neuro-energetic mechanisms in the brain [12], it is undoubtedly an appropriate candidate to trace differential metabolic changes associated with cognitive activity. Another appealing advantage of 99mTc-HMPAO SPECT is its application in the diagnosis of neurodegenerative diseases [13]. There is evidence showing that 99mTc-HMPAO SPECT allows to discriminate between Alzheimer’s disease (AD) on the one hand, and frontotemporal dementia [14] or Lewy bodies disease [15, 16] on the other hand. More, and of greater interest to the present purpose, there is very recent evidence showing that 99mTc-HMPAO SPECT allows to discriminate early dysfunction associated with AD, most particularly when the brain is under cognitive “stress”, as reached using cognitive activation [17, 18]. Thus, not only there is a need to validate the use of 99mTc-HMPAO SPECT for activation studies in healthy young populations, but also to offer a protocol easily applicable in clinical settings for investigations of various states, among which normal and pathological aging.

Another challenge of the present study was to considerably reduce the recommended 555 MBq [7, 8] or 370 MBq [4, 5] of radioactivity to 270 MBq per scan, yet preserving the high quality of the images. We used a 1 to 1 dose ratio across scans as commonly applied in 99mTc SPECT activation studies [4, 5, 7, 8], although other methods have been proposed, such as increasing the second dose by a factor of approximately 2.0 [19, 20]. Notably, while the latter is meant to increase the signal to noise ratio, lowered by the pile-up effect inherent to the split-dose protocol, the former allows a larger reduction of the subject’s exposure to radioactivity. Moreover, images can be post-processed using proportional scaling to reliably grasp between-conditions differential activities [21], despite the pile-up effect (see the post-processing and statistical analyses section). For all these reasons, and because of particular radioprotection concerns, we favoured a 1 to 1 ratio procedure.

Finally, we addressed the feasibility of the activation protocol by convergent validity of the results. We used a cognitive activation task, the n-back, for which a large ensemble of comparable results has been already reported in the neuroimaging literature [22-26]. This task is currently one of the paradigms for which the most precise knowledge of the underlying functional architecture is available [27]. Replicating the findings reported both with 15O-PET and fMRI would support the validity and the feasibility of cognitive activation studies with low-dose 99mTc-HMPAO SPECT.

Materials and methods

Subjects

The population consisted of 12 healthy young adults (10 women, 2 men), with an average age of 26.42 years (SD = 2.57, range 23 to 31) and an educational level of 20.08 years (SD = 1.98, range 17 to 24). All participants were right-handed (as assessed by the Oldfield-Edinburgh Inventory [28]) and fluent in French. None reported a history of medical, neurological or psychiatric disease. The participants were recruited in the university community by means of advertisements, and all signed written informed consent for participation. The experiment was approved by the ethics committees of the Geneva University Hospitals and of the Faculty of Psychology and Educational Sciences at the University of Geneva.

Neuropsychological procedure

The verbal n-back task [26] was administered to the participants. The subjects were visually presented with series of letters (upper and lower case consonants), displayed one after the other at the centre of the computer screen. The 0-back condition, which served as a reference, consisted in deciding whether each stimulus corresponded to a specific target (letter “X”). The 2-back condition consisted in deciding whether each stimulus corresponded to the antepenultimate one in the sequence. In each condition, 6 blocks of 36 trials were built, for a total of 216 trials. The trial sequence was made up of a 500 ms stimulus display followed by a 2,500 ms fixation point (see Fig. 1). One third of the items were selected as targets, randomly arranged in each block. Within each condition, block and trial order were maintained constant across participants. The 0-back was always provided before the 2-back. The task was computerized and built using E-Prime Software (Psychology Software Tools, Pittsburgh, PA, USA). The script was run on a Dell (Round Rock, TX, USA) personal computer provided with an E-Prime Serial Response Box device and a 15-in. monitor with 640 × 480 pixels resolution. For each trial, time and accuracy of responses were recorded online by the computer.

Fig. 1
figure 1

Schematic representation of the trial sequences used in the 0-back condition (left panel) and in the 2-back condition (right panel)

Nuclear medicine procedure

Each of the two experimental conditions of the cognitive task was associated with a single full brain scan recording. A split-dose technique was applied to the participants using 270 MBq (7.3 mCi) of 99mTc-HMPAO per scan. 99mTc-HMPAO was prepared according to the recommendations of the manufacturer (Ceretec, Amersham, UK) and intravenously injected no longer than 20 min after preparation. Both tasks were performed in the scanning room but outside the scanner apparatus, providing settings close to laboratory conditions. The participants were invited to sit comfortably in front of the computer, with the left arm lying onto a small table. An intravenous butterfly catheter was placed in the elbow fold in preparation for subsequent radiopharmaceutical intravenous injection. The task instructions for the 0-back condition were provided, and the six experimental blocks were administered. The tracer was injected at the beginning of the third block. Due to the rapid first-pass extraction of 99mTc-HMPAO [29], brain images reflect the average activity that was labelled during the third and fourth blocks. At the end of the task, the participants were required to rest during the washout period (at least 20 min, see Table 1) before being placed in the scanner in a decubitus position. Image acquisition lasted 30 min. The same procedure was applied for the acquisition of the second brain image associated with the 2-back condition. The interval between the end of the first scan and the acquisition of the second one was approximately 50 min. It is worth pointing out that inherent to the split-dose single day protocol is contamination of the second scan by the remaining activity of the first one [5, 6]. However, it has been argued that the pile-up does not prevent the conduction of neurofunctional activation studies [4, 5, 21], particularly if post-processing involves proportional scaling procedures. Indeed, the latter allows statistical control for global variations in the signal across scans and conditions [21].

Table 1 Quality control values, by scan and participant

Quality control procedure

For each scan, quality controls concerned the dose preparation, the effective dose administered to the subjects, the injection/acquisition delay and the total counts obtained in the reconstructed images. For each dose administered, the syringe radioactive content, in MBq, was measured by means of standard laboratory procedures. Pre- and post-injection counts were subsequently acquired by recording the radioactivity of the source using the Toshiba-GCA gamma camera, calibrated on 140 KeV. Acquisition was performed in a fixed mode for 5 min. A region of interest (ROI) was manually drawn around radioactive foci. To complete quality control measures, the time between injection and acquisition was recorded.

Acquisition protocol

For both scans, SPECT acquisition was performed at least 20 min after i.v. injection using a triple head Toshiba GCA-9300 gamma camera (Toshiba Medical Systems, Tokyo, Japan) with a fan beam high resolution collimation. Data were acquired on a 128 × 128 matrix, with a 360° rotation angle and a step-and-shoot mode (6° step angle and 90 s rotation time), providing 60 projections.

Given that our aim was to develop a protocol that could be easily implemented in clinical settings, and despite the recent development of alternative methods [19], image processing was conducted by means of methods currently available in most commercially distributed processing software. These methods are routinely performed, and have been discussed in the literature for SPECT activation protocols [4, 5, 7, 8]. Thus, scatter correction was performed using the TEW method [30], with 20% set around the 140 KeV peak, and an additional 7% in the upper and lower sides of the main window. Images were reconstructed in a 128 × 128 matrix. Attenuation correction was done by filtered back projection (Sorenson method) using a Butterworth filter of 8.0.16 for the main window and of 8.0.09 for the sub-windows. The attenuation correction map was automatically drawn with a contour set at 5%, and an attenuation coefficient set to 0.09 cm−1. Image reconstructions used a Shepp & Logan filter in the axial, coronal and sagittal planes with a slice thickness of 1 pixel. Voxel size was 3.44 mm3. Examples of reconstructed images are provided in Fig. 2.

Fig. 2
figure 2

Example of raw SPECT images of the 0-back and 2-back conditions in the coronal and axial planes (images are displayed according to the radiological convention)

Post-processing and statistical analyses

The original Toshiba raw data files were converted into Analyze format using MedX software (Sensor Systems, Sterling, VA, USA). By means of SPM2 (Department of Cognitive Neurology, Institute of Neurology, London, UK) implemented in MatLab 6.5 (Math-Works, Sherborn, MA, USA), the different images of a participant were realigned, creating a mean volume of resliced scans. Images were further normalized into the Talairach and Tournoux stereotaxic space [31], using the mean images to obtain the parameters’ values and applying to the resliced images a trilinear interpolation with eight non-linear iterations. Co-registration, realignment and normalization were conducted to match each scan onto a common reference space, hence minimizing individual differences in anatomy and allowing subsequent voxel-wise analyses. Finally, images were smoothed using a Gaussian filter of 10 × 10 × 12 mm3 at full width half maximum (FWHM). Subtraction method was applied to count densities on a voxel-by-voxel basis to determine brain areas of significantly increased or decreased activity. Differences in global activity between scans were removed by scaling the activity in each voxel proportional to the averaged global activity [21]. Grand mean values were computed using a grey matter threshold of 0.5. This proportional scaling procedure not only allows statistical control for differences in global scan activity and cross-scan activity accumulation, but also increases the sensitivity in detecting small variation in the signal associated with the experimental conditions. Finally, to obtain a more readily comprehensible unit of activity, the corrected values were adjusted around a grand mean of 50 ml/min/100 g, corresponding to the approximated average blood flow as measured by 15O-PET. The effect of task conditions (2-back vs 0-back) on the count density at each voxel was estimated using a general linear model, wherein the changes in global counts are considered as a covariate. Linear contrasts of the means across conditions were used on a voxel-by-voxel basis using t statistics. The resulting sets of t values constituted statistical parametric maps SPM(t), that were subsequently transformed to unit normal distribution SPM(Z) maps.

Results

Quality control

Quality control values are reported in Table 1. Of particular interest is the low inter-subject variability in (a) the radioactivity values collected by means of laboratory procedures, (b) the percentage reduction in counts post-injection and (c) the time interval between injection and acquisition. Furthermore, these values measured for scan 1 and scan 2 were not significantly different, as assessed by paired-samples t tests conducted on each variable separately. All analyses yield p values larger than 0.15. Finally, it is important to stress that all images acquired were of good quality (see Fig. 2 as an example), and all were taken into account for further analysis.

Imaging data

Inferential statistics on brain imaging data were conducted using SPM2 (Department of Cognitive Neurology, Institute of Neurology, London, UK). A multi-subjects, condition, and covariates design was applied, with two conditions (0-back, 2-back), but no behavioural covariates. The analysis was performed without correction for multiple comparisons, using a p value threshold of p = 0.001 and a cluster of at least 20 voxels for the peaks to be considered. Detailed results are reported in Table 2. Given that the sample included two men and ten women, and that controversies exist in the literature with respect to gender differences [32, 33], we individually compared each man’s pattern of activity to the average pattern found in the group of women using the method proposed by Signorini [34]. No significant differences were found. Therefore, we pooled all scans together to increase statistical power and report the findings relative to the complete sample of 12 participants. Of direct interest to the present purpose is the increased activity specifically associated with the 2-back condition [for the corresponding SPM(Z) maps, see Fig. 3]. The peaks occur in an ensemble of areas involving, on the left side, the frontopolar (BA10) and dorsolateral (BA46) regions of the prefrontal cortex, the frontal eye field (BA8), and the left inferior parietal lobule (BA40). Additionally to these left-sided foci, areas of the right hemisphere demonstrated specific activations during the 2-back task; these include the orbitofrontal (BA11) and frontopolar (BA10) cortices as well as the somatosensory and motor primary regions (BA3/4) and the superior parietal lobule (BA7). Finally, increased bilateral activity was found in the visual association areas, mainly involving the cuneus (BA18).

Table 2 Peak coordinates (x, y, z) of significant condition effects, with associated Z value and anatomical areas
Fig. 3
figure 3

SPM(Z) maps for the 2-back minus 0-back contrast, displayed according to the neurological convention

Behavioural data

Analysis of the behavioural data was performed on the median response times, computed only on correct responses and on the proportion of correct responses. Scores were averaged for blocks 1 and 2, blocks 3 and 4 and blocks 5 and 6, corresponding to the events occurring before injection and tracer uptake, during uptake and after uptake, respectively. Descriptive statistics by Period (before, during, after) and Condition (0-back, 2-back) are graphically displayed in Fig. 4.

Fig. 4
figure 4

Behavioural results: median response times (RT, bars) and proportion of correct responses (Pcor, lines) before, during and after uptake, for the 0-back and the 2-back conditions

Inferential statistics were conducted using a 2 × 3 ANOVA design, with the Condition (0-back, 2-back) and the Period (before, during, after uptake) as repeated measures. With respect to the median response times, results revealed a significant main effect of Condition F(1,11) = 22.11, p < .001, η p 2 = 0.67, but the main effect of Period, F(2,22) = 0.64, p = 0.54, and the Condition × Period interaction, F(2,22) = 0.80, p = 0.46, were not significant. As concerns the proportion of correct responses, a similar pattern was observed with a significant main effect of Condition, F(1,11) = 7.28, p < .05, η p 2 = 0.40, but no significant main effect of Period, F(2,22) = 2.04, p = 0.15, and no significant Condition × Period interaction, F(2,22) = 1.10, p = 0.35. Altogether, these results are in line with previous findings reported in the literature [35, 36] and support the hypothesis of an increase in cognitive demand associated with the 2-back condition, as compared to the 0-back. Indeed, in the former, participants take more time to respond and demonstrate lower accuracy. It is also worth pointing out that performance remains stable across the three time periods, with similar patterns in both tasks.

Discussion

Findings from the present experiment using a 99mTc-HMPAO SPECT split-dose activation protocol are rather encouraging with respect to the two goals pursued. First, from a technical and methodological point of view, these results were supportive of the use of 99mTc-HMPAO as a radiotracer labelling differences in metabolism during cognitive activation and of the possibility of collecting good quality brain images despite relatively low doses of radiopharmaceuticals. Second, results regarding the localizations of brain activations during a verbal working memory task were in good agreement with those reported in the literature. Convergent validity was adequately demonstrated, hence providing additional support for the feasibility of cognitive activation studies with 99mTc-HMPAO SPECT.

Concerning the behavioural data, performance was significantly lower in the 2-back condition than in the 0-back one, which replicated previous findings [35, 36] and supported the assumption that the 2-back was more demanding in processing resources from the executive component of working memory [23, 24]. Involvement of working memory was further supported by brain imaging results, which revealed that the 2-back condition was associated with increased activity in left prefrontal areas (extending to the right hemisphere), the left parietal lobule (BA40) and bilateral occipital areas. As pointed out by Owen et al. [27] in a recent meta-analysis of functional imaging studies using the n-back task, the fronto-parietal network involvement in the 2-back task is currently largely undisputed. The fact that the major brain areas reported in the literature using 15O-PET or fMRI were also found in the present study strongly supports the validity of the measures acquired with 99mTc-HMPAO SPECT, and provides good evidence for a reliable application of the present protocol for brain imaging experiments.

Two major issues need further discussion: the advantages vs. limitations of SPECT for cognitive activation studies, and the choice of HMPAO over ECD as ligand to assess brain energy metabolism. Although widely applied in clinical routine, SPECT has been rarely used for activation studies in the field of cognitive neurosciences compared to fMRI, or even 15O-PET. Both latter techniques benefit from almost real-time imaging of the brain, allow repetitions of intra-individual scan acquisitions and offer a precise spatial resolution. However, they require a fairly rigid experimental set-up, and the participant has to perform the cognitive test within the scanner and remain in a decubitus position throughout the whole experiment. Such technical and environmental constraints not only restrict comfort for the participant but also prevent assessment of individuals who cannot remain still in the scanner for a long time. Even more important is the non-negligible disruptive impact of the fMRI noisy environment on the cognitive processing under investigation [1]. This effect may possibly compromise between-group comparisons, given that various pathological states, or even advancing age, are related to differences in the susceptibility to interfering sources of information [37]. Hence, with respect to this particular issue, SPECT undoubtedly demonstrates large advantages by allowing cognitive assessment to be conducted outside the scanner. Another major advantage of the protocol we propose is its potential application in clinical populations, more specifically in cases of neurodegenerative diseases. Indeed, there is increasing evidence that 99mTc-HMPAO SPECT permits discrimination of functional specificities across different types of dementias [14, 15] even before abnormalities are reported with structural imaging. Of additional interest is very recent evidence demonstrating that early dysfunction found in Alzheimer’s disease is best apprehended during cognitive activation as compared to rest [17, 18]. Hence, although 99mTc-HMPAO SPECT might not be the most appropriate technique to investigate the neurofunctional bases of cognition in young healthy populations, its application might become particularly relevant for the study of normal and pathological aging, and the early detection of functional deficits associated with neurodegenerative diseases.

SPECT does have several limitations. Among these are the small number of intra-individual brain recordings and the limited ability, due to its moderate resolution, to detect signal differences between individuals. Although these limitations need to be recognized, encouraging findings already support the use of SPECT for activation studies [4, 5, 7, 8], which should not be dismissed. Moreover, previous studies have demonstrated that the moderate resolution of SPECT does not prevent applying voxel-based analyses to SPECT images, and that approaches such as statistical parametric mapping (SPM) can adequately be employed [2, 21], even with rather small samples [38]. The findings from the present study also confirm the feasibility of activation protocols with SPECT, despite intrinsic limitations.

The present work is also one of the first contributions providing a clear demonstration that SPECT can be successfully performed with 99mTc-HMPAO in activation experiments. To date, at least in healthy young populations, most comparable studies exclusively used 99mTc-ECD, probably because ECD has long been considered to reflect regional cerebral blood flow (rCBF) better than 99mTc-HMPAO. Recent work, however, has pointed out a cardinal linkage between HMPAO retention dynamics and neuro-energetic mechanisms in the brain [13, 39]. More specifically, HMPAO is prominently retained in astrocytes [12], a cell type for which a central role in coupling neuronal activity to glucose utilization [40] and blood flow [41] has been demonstrated. Targeting this particular metabolic dynamics supports favoring HMPAO over ECD to precisely trace metabolic changes related to cognitive activity, and particularly those associated with memory [42, 43]. In this vein, and as already underlined earlier, the use of HMPAO may be very helpful in tracing dysfunctions associated with pathological processes that affect glial cells, such as alterations of astrocyte function in degenerative diseases [44, 45], or changes in microglia occurring in the course of normal aging [46].

The present study provides substantial support for the feasibility of a split-dose 99mTc-HMPAO SPECT activation protocol and reveals its ability to identify brain regions involved in cognitive activity. Our findings further warrant the use of HMPAO as a radioligand. Considering the emerging role of neuron-astrocyte interactions in coupling synaptic activity with metabolic and vascular responses [47], HMPAO may be better suited than ECD to capture neuro-energetic modulations associated with cognitive activity. Thus, 99mTc-HMPAO SPECT activation protocols may attract increasing interest not only in the study of neurofunctional bases of cognition, but also, and most importantly, in further investigations of neuro-energetic dysfunctions in normal and pathological aging. Finally, the results broadly support performing brain-imaging experiments in settings close to standard clinical ones and associating low irradiation with high-quality images. In conclusion, this empirical validation study has shown that functional imaging studies of cognition using 99mTc-HMPAO SPECT are both technically and methodologically feasible.