PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(9): e0184743.
Published online 2017 September 13. doi:  10.1371/journal.pone.0184743
PMCID: PMC5597218

Simultaneous trimodal PET-MR-EEG imaging: Do EEG caps generate artefacts in PET images?

Ravichandran Rajkumar, Conceptualization, Formal analysis, Investigation, Project administration, Writing – original draft, Writing – review & editing,1,2,3 Elena Rota Kops, Conceptualization, Investigation, Project administration, Software, Writing – review & editing,1 Jörg Mauler, Investigation, Validation,1 Lutz Tellmann, Investigation, Methodology,1 Christoph Lerche, Project administration, Resources, Validation, Writing – review & editing,1 Hans Herzog, Supervision, Writing – review & editing,1 N. Jon Shah, Funding acquisition, Resources, Supervision, Writing – review & editing,1,3,4,5 and Irene Neuner, Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing1,2,3,*
Kewei Chen, Editor

Abstract

Trimodal simultaneous acquisition of positron emission tomography (PET), magnetic resonance imaging (MRI), and electroencephalography (EEG) has become feasible due to the development of hybrid PET-MR scanners. To capture the temporal dynamics of neuronal activation on a millisecond-by-millisecond basis, an EEG system is appended to the quantitative high resolution PET-MR imaging modality already established in our institute. One of the major difficulties associated with the development of simultaneous trimodal acquisition is that the components traditionally used in each modality can cause interferences in its counterpart. The mutual interferences of MRI components and PET components on PET and MR images, and the influence of EEG electrodes on functional MRI images have been studied and reported on. Building on this, this study aims to investigate the influence of the EEG cap on the quality and quantification of PET images acquired during simultaneous PET-MR measurements. A preliminary transmission scan study on the ECAT HR+ scanner, using an Iida phantom, showed visible attenuation effect due to the EEG cap. The BrainPET-MR emission images of the Iida phantom with [18F]Fluordeoxyglucose, as well as of human subjects with the EEG cap, did not show significant effects of the EEG cap, even though the applied attenuation correction did not take into account the attenuation of the EEG cap itself.

Introduction

The modalities of positron emission tomography (PET), magnetic resonance imaging (MRI), and electroencephalography (EEG) were combined into a single trimodal approach in order to best utilise their complementary features. The trimodal imaging set-up yields high spatial resolution MR images, high temporal resolution EEG signals, and metabolic PET images simultaneously in a single scanning session [1]. The feasibility of measuring PET-MR-EEG data simultaneously in single session has already been reported [24] and the benefits of such a set-up are numerous. For example, the use of a PET-MR-EEG set-up can advance the understanding of brain structure and function and can be used to identify new biomarkers for neurological and psychiatric disorders, as well as improving accuracy in early diagnosis [5]. A potential in vivo marker for neocortical neuronal loss in Alzheimer’s disease was found based on evidence from PET, MRI, and EEG [6]. Furthermore, this setup is especially desirable for pharmacological challenge studies where a sequential design would add confounding factors [7]. A PET-MR-EEG set-up could also be used to improve pre-surgical evaluation and post-surgical monitoring in epilepsy patients in order to identify and localise the epileptogenic zone more quickly [8,9]. Further possibilities relating to the potential application of multimodal neuroimaging are well reviewed and discussed in the literature [1012]. Interpretations and quantitative results from trimodal studies should, however, consider mutual interferences between the three modalities. In the literature, studies on the mutual interferences of PET and MR components in hybrid PET-MR scanners on MR and PET images [1316], on the successful removal of MR-related artefacts from EEG signals [1719], as well as on the influence of EEG electrodes on functional MRI (fMRI) blood oxygenation level dependent (BOLD) signal [20] have already been reported on. Simultaneous PET-EEG measurements were performed by a number of groups using PET only scanners [2124] and applying 68Ge-transmission based attenuation correction. It is not necessary to investigate and analyse the attenuation effect of an EEG cap while using PET only scanners for simultaneous PET-EEG measurements because the attenuation map generated during transmission scans will automatically incorporate the attenuation details of the EEG cap. Thus, final reconstructed images from PET only scanners do not show any attenuation effect due to the components of the EEG cap. This is, however, not the case in hybrid PET-MR scanners, where acquisition of transmission images is not possible. In the hybrid BrainPET-MR scanner, the attenuation correction is performed by using an MR based attenuation correction method [25]. To our knowledge, an explorative study investigating and quantifying the effect of the components of an MR-compatible EEG cap during a simultaneous trimodal study on PET images has not yet been performed.

The MR compatible EEG cap, used during trimodal studies for detecting brain activity signals from the scalp, contains tiny silver/silver chloride (Ag-AgCl) metal electrodes embedded in a plastic housing made of polycarbonate and acrylate. The 511 keV photons emitted by the PET radiotracer may be attenuated by the tiny metal electrodes, resistors, cables and the plastic housing in the EEG cap, leading to possible artefacts in the PET images. A previous study [26] relating to the impact of metal artefacts due to EEG electrodes in brain positron emission tomography / computed tomography (PET/CT) imaging, reported that EEG electrodes gave rise to local hot spots and to a significant quantification bias in the PET/CT brain images. These errors, however, did not lead to a modification of the visual interpretation of the brain images or the subsequent diagnoses. However, it should be taken into account that EEG Genuine Grass Precious Gold Disc Electrodes with a disc diameter of 10 mm, supplied by Astro-Med, Inc. (West Warwick, RI, USA), were used in this study, and such disc electrodes are not used in our trimodal studies [2,3]. Furthermore, the attenuation correction was based on CT images and is more sensitive to attenuating material in the PET/CT field of view (FOV). Hence, the study results obtained from PET/CT scanners with a different type of EEG electrodes are not applicable to trimodal PET-MR-EEG. The study presented here, however, is designed to investigate and quantify the influence of EEG caps on PET images measured using hybrid PET-MR scanners, for which the attenuation maps do not include any attenuation information related to the EEG cap.

Materials and method

A MR compatible BrainCap MR EEG cap (EASYCAP GmbH, Herrsching, Germany) was used to investigate the attenuation effect. The components in the 32 channel BrainCap MR EEG cap and the materials these components are made of are shown in Table 1. In addition to the major components shown in Table 1, the BrainCap MR EEG cap also contains chip resistors (~ 2 mm length, ~ 1 mm breadth) and minute soldering (~ 1 mm length) of copper wire connections between chip resistor and electrode and between chip resistor and thin cables.

Table 1
BrainCap MR components and materials.

Measurements using the Iida PET phantom [27] (Fig 1) and the BrainCap MR EEG cap were performed with the ECAT EXACT HR+ scanner (CTI/Siemens, Knoxville, TN, USA; 32 rings, axial field of view 15.5 cm [28]) and with the 3T-BrainPET-MR scanner (Siemens, Erlangen, Germany; 72 rings, axial field of view 19.2 cm [29,30]).

Fig 1
Iida phantom (upper row) without (left) and with (right) MR compatible BrainCap MR EEG cap.

Phantom data acquisition and reconstruction

The BrainCap MR EEG cap was placed on the Iida PET phantom (Fig 1) and a 60 minute transmission scan was performed in the ECAT HR+ scanner. The same measurement was repeated without the BrainCap MR EEG cap. The transmission data were iteratively reconstructed (OP-OSEM, 16 subsets, 6 iterations) yielding attenuation maps with and without the EEG cap (AMw and AMwo).

PET emission measurements with simultaneous MR measurements were performed in the 3T-BrainPET-MR scanner on the Iida phantom filled with about 550 ml of radioactive solution (the radioactive solution was prepared by adding 376 MBq [18F]Fluorodeoxyglucose ([18F]FDG) to one litre water). The PET data of the Iida phantom were also acquired in two sessions, one with the BrainCap MR EEG cap and one without the EEG cap. The data were acquired in list mode and the duration of each session was 30 minutes. The MR protocol incorporated in each session during simultaneous PET-MR measurement and their parameter are as follows:

  • magnetisation-prepared rapid acquisition gradient-echo (MP-RAGE) sequence (T1 weighted anatomical imaging sequence, repetition time (TR) = 2250 ms, echo time (TE) = 3.03 ms, field-of-view (FOV) 256 × 256 × 176 mm3, voxel size = 1 x 1 x 1 mm3, flip angle (FA) = 9°, 176 sagittal slices and a GRAPPA factor of 2 with 70 auto-calibration signal lines)
  • ultra short echo time (UTE) sequence (T2* weighted sequence for attenuation correction), TR = 200 ms, TE1 = 0.07 ms, TE2 = 2.46 ms, FOV = 320 x 320 x 320 mm3, voxel size = 1.67 x 1.67 x 1.67 mm3, FA = 15°)

The image data were corrected for random coincidences, scatter coincidences, attenuation, and dead time prior to reconstruction. The frame length for the reconstruction of each session was estimated in such a way that the total number of prompt counts within the chosen frame length is equal for both sessions. This was done in order to eliminate the effect of change in radioactivity level between sessions due to Fluorine-18 (18F) physical decay. Since FDG in the phantom does not undergo any metabolic changes, such a correction approach is feasible. The estimated frame lengths are 1095 s for the first session with EEG cap and 1300 s for the second session without EEG cap. The start offset was 100 s for both sessions. Reconstruction was done using the OP-OSEM reconstruction algorithm provided by the manufacturer (2 subsets, 32 iterations); finally, a 3D Gaussian filtering with 3 mm kernel was performed. The attenuation map AMwo provided the attenuation correction for both sessions’ data, simulating the real situation of human subject scans, for which attenuation maps of EEG caps cannot be produced in routine study protocols. The acquisition and application of an individual attenuation mask for each volunteer would require additional measurement time in the PET-MR scanner. This would interfere with the optimal timing for bolus infusion in PET-MR protocols and the application of radiotracers [31], thus increasing the radiation exposure inadequately. In addition, due to head motion creating mutual small registration mismatches, an individual mask would not guarantee a more precise PET measurement. Thus, the PET emission images were reconstructed using AMwo. The exact position of AMwo for the BrainPET images was obtained by matching the AMwo to the anatomical MR images of the Iida phantom, obtained during the simultaneous PET measurement.

Human data acquisition and reconstruction

A similar experiment with two sessions was performed on eight healthy human subjects (age: 28 ± 3.6 years) injected with 202.8 ± 35.1 MBq of [18F]FDG at the start of the trimodal experiment and while the subjects were lying inside the BrainPET-MR scanner. The first dynamic emission scan with an EEG cap was started simultaneously with the injection of [18F]FDG for a duration of 60 minutes (measurement with cap—MwC); a predefined MR protocol, adapted for measurement with EEG cap, started at the same time. The sequences incorporated in the MR protocol and their parameters are as follows:

  • MP-RAGE (same parameters as the mesurement done on the phantom)
  • UTE (same parameters as the mesurement done on the phantom)
  • spectroscopy (Point Resolved Spectroscopy (PRESS) sequence, TR = 2500 ms, TE1 = 14 ms, TE2 = 105 ms, number of averages (NA) = 128, voxel size = 25 × 25 × 25 mm3, RF pulse centred at 2.4 ppm, 16 step phase cycling)
  • echo planar imaging (T2*-weighted EPI sequence, TR = 2200 ms, TE = 30 ms, FOV = 200 x 200 x 108 mm3, voxel size = 3.125 x 3.125 x 3.0 mm3, FA = 80°, number of slices = 36, number of volumes = 165, eyes closed resting state measurement)

EEG data were also recorded during this simultaneous measurement. After MwC, the subject was taken out from the BrainPET-MR scanner and the EEG cap was removed. After an interval of 17:08 ± 3:56 minutes from the end of MwC, a second emission scan over 20 minutes without the EEG cap was performed (measurement without cap—MwoC) with a predefined MR protocol. The MR protocol during MwoC comprised diffusion tensor imaging (DTI) using a standard double-refocused spin-echo EPI sequence with bipolar gradient pulses (TR = 9100 ms, TE = 87 ms, FOV = 240 x 240 x 136 mm3, voxel-size = 1.9 × 1.9 × 1.9 mm3, flip angle = 90°, b-values = 0, 1000 s mm-2, number of slices = 72, number of field gradient directions = 30, NA = 4). The protocols during MwC and MwoC did not include any task. The measurements on the human subjects were part of an ongoing study approved by the Ethics Committee of the Medical Faculty of the RWTH Aachen University, Germany. No additional activity was injected for the purpose of the report presented here. Written informed consent was obtained from all subjects and the study was conducted according to the principles expressed in the Declaration of Helsinki.

The human PET emission data were reconstructed with the same protocol as used for the phantom data (OP-OSEM reconstruction provided by the manufacturer with 2 subsets and 32 iterations). The count rate reduction in PET caused by MR sequences such as EPI and DTI sequences [14] can be, and was indeed corrected by interpolation. The attenuation corrections were performed with template-based attenuation maps [25], which do not contain any attenuation information from the EEG cap. The template-based attenuation map was obtained from a template attenuation map with a corresponding MR template. The former was built as an average image from eight available transmission images performed with the ECAT Exact HR+ scanner (reconstruction: OP-OSEM, 16 subsets, 6 iterations). The MR template was generated with the corresponding and co-registered MR images. The idea of the template-based attenuation correction is based upon the possibility of adapting the MR template to the patient MR image with nonlinear registration (SPM8; Wellcome Department of Neurology, London, UK). The resulting nonlinear transformations are applied to the template attenuation map finally adapting it to the PET data of the subject investigated. Further processing steps were done using the PMOD software package (version 3.4, Zurich, Switzerland), the MPI tool (version 6.48, Advanced Tomo Vision, Kerpen, Germany), the Matlab software package (version 8.5 (R2015a)) and the SPM tool box (version 12).

Phantom data pre-processing

In order to evaluate the attenuation caused by the electrodes of the EEG cap on the ECAT HR+ attenuation image, the relative difference image in percentage was calculated between AMw and AMwo using the formula:

AMwAMwoAMwo100
(1)

The attenuation assessment was done by defining six different VOIs at the position of electrodes on the relative difference attenuation image in percentage. The VOIs were drawn using the tools available in the PMOD software package; the shape and size of the VOIs were drawn to exactly match the attenuation spots on the difference image.

The reconstructed and smoothed FDG-PET Iida phantom image without BrainCap MR EEG cap (PIwoC) was rigidly matched to the FDG-PET Iida phantom image with BrainCap MR EEG cap (PIwC) using the PMOD software package. In order to reveal qualitative differences between PIwoC and PIwC, two professionally experienced physicians were asked to compare the reconstructed emission images of the Iida phantom; the comparisons in each category were blinded.

In order to be able to compare the two images for quantitative assessment, PIwoC and PIwC were normalised to Z-score images relative to the global mean and standard deviation. A relative difference image between Z-score normalised PIwoC and PIwc in percentage was created using the formula:

PIwCPIwoCPIwoC100
(2)

A whole phantom binary mask, covering only the grey matter (GM) area in the phantom, was created by thresholding PIwoC using the SPM software package for further processing and calculations. Eight other spherical volumes of interest (VOI) of 9 mm radius were created using the FSL software package. These eight VOIs were directly positioned below the electrodes; the position of the electrodes was determined by rigidly matching the transmission scan images to the phantom emission images. In order to further restrict the calculations to the GM regions of the phantom, only the VOI voxels within the GM region were considered (GM correction). This was achieved by using the whole phantom binary GM mask created previously. In order to quantify the attenuation effect of the components of the EEG cap, the statistical parameters such as mean, median, standard deviation (SD) and interquartile range (IQR) were calculated from the relative difference image in percentage for the whole phantom but only in relation to the grey matter region and for eight other spherical GM corrected VOIs created before. The statistical parameters were calculated using the Matlab software package.

FDG-PET human data pre-processing

FDG-PET emission data acquired during MwC were reconstructed into 18 different time frames and an average over the last four frames, with a frame length of 5 minutes each—corresponding to 40–60 minutes after [18F]FDG injection, was considered for further analysis (FDG-PET image with cap—IwC). Similarly, FDG-PET emission data acquired during MwoC were reconstructed into 4 frames with a frame length of 5 minutes each and averaged into one image (FDG-PET image without cap—IwoC). Before averaging, both IwC and IwoC were smoothed with a 3 mm Gaussian filter and motion corrected; IwoC was decay corrected back to IwC scan time and both normalised to standard uptake value (SUV) images, with respect to body weight and total injected activity of the individual subject. Since FDG uptake is very high in the grey matter of the brain, further processing and computations are done only on the grey matter regions. In order to accomplish this, individual grey matter masks (GM mask) of each subject were created from the IwC co-registered MR structural image (MPRAGE) using the brain extraction tool (BET) [32]; the extracted GM mask then was binarised at a threshold of 10%.

The IwoC was rigidly matched to IwC using PMOD software package and the relative difference image in percentage between IwC and IwoC was calculated using the formula:

IwCIwoCIwoC100
(3)

The relative difference image in percentage was corrected for grey matter using the binarised GM mask created before. The processing steps are illustrated in Fig 2.

Fig 2
Flow chart showing preprocessing steps for FDG PET human data.

Results

FDG-PET phantom and human subject data were pre-processed as explained in the methods section and the results for phantom and human subject data are presented below accordingly.

Phantom results

The ECAT HR+ relative difference attenuation image in percentage between AMw and AMwo of the Iida phantom is shown in Fig 3. The attenuation assessment at six different VOIs showed a group average relative difference of about 41.6 ± 24.5% directly on the accessed electrode positions.

Fig 3
Attenuation maps of the Iida phantom with EEG cap (AMw, upper row) and without EEG cap (AMwo, middle row); relative difference image in percentage (bottom row).

When comparing the reconstructed emission images of the Iida phantom measured in the 3T-BrainPET-MR scanner, the professionally experienced physicians could not detect any differences between the emission images obtained with and without the EEG cap. In order to underpin this, a relative difference image in percentage between the emission images of the Iida phantom with and without an EEG cap was created (Fig 4) and visually inspected. The difference image showed no artefacts due to the presence of the EEG cap components.

Fig 4
Emission images of the Iida phantom with EEG cap (a), and without EEG cap (b), both reconstructed using the attenuation map without EEG cap; relative difference image in percentage (c).

In order to gain a visual impression of the position of the electrodes in relation to the PET emission images of the Iida phantom, the PET emission images were co-registered to the transmission image and all the slices were visually inspected for mutual artefacts due to EEG cap components. Fig 5 shows three different axial slices from the emission image of the Iida phantom with EEG cap, together with the corresponding axial slices of transmission and emission relative difference images in percentage.

Fig 5
Axial slices at three different axial positions of the emission image of the Iida phantom with EEG cap (upper row), corresponding axial slices of the relative difference images from Fig 3 (transmission image, middle row), and from Fig 4 (emission image, ...

The statistical parameters were calculated from the relative difference image in percentage for the whole phantom (only grey matter region) and for the eight spherical VOIs created previously; they are summarised in S1 Table and a visual representation is shown in the bar plot (Fig 6).

Fig 6
Bar plot showing the mean (green bar) and standard deviation (black lines) of the relative difference emission image of the phantom for whole phantom grey matter and eight spherical VOIs under the position of the electrodes.

FDG—PET human data results

As stated for the phantom data, the visual inspection of the relative difference image in percentage between IwC and IwoC (both reconstructed with AMwo) of the human subjects also did not show any visible artefacts due to the EEG cap components. Emission PET images with and without EEG cap and the relative difference image of a representative subject is shown in Fig 7. Relative difference images in percentage of the remaining seven subjects are shown in S1 Fig. The statistical parameters calculated from human emission relative difference images in percentage are given in the S2 Table. The histogram plot of the voxel values in the human relative difference images is shown in S2 Fig.

Fig 7
FDG PET SUV image of a human subject with EEG cap (top row) and without EEG cap (middle row), both reconstructed with an AC map without EEG cap; relative difference image in percentage between the above FDG PET images (bottom row).

Discussion and conclusions

The transmission scan performed with the ECAT HR+ scanner on the Iida phantom with the EEG cap clearly shows attenuation on the surface of the phantom scalp caused by the components of EEG cap (Fig 3). The presence of the tiny metal electrodes and chip resistors, plastic housing and other minute soldering materials causing attenuation, generate a higher attenuation correction factor (ACF) for the affected line of response (LOR), consequently leading to an overcorrection along these LORs. The result is a possible overcorrection in the areas directly below the electrodes. In contrast, Figs Figs44 and and55 clearly illustrate that PET emission images of the phantom showed no significantly visible artefacts due to the presence of EEG cap. The quantitative assessment on the relative difference in percentage of the phantom PET emission images showed the mean values of 0.45(±11.21) % and maximum of -2.36 (±11.04) % on the whole phantom grey matter region and over the eight spherical VOIs, respectively. However, it cannot be concluded that the relative differences shown are only due to the presence of EEG cap. Other factors contributing to this relative difference may be the influence of count rate, statistical uncertainties arising from both the measurement and the reconstruction of PET images, and mismatches due to the rigid matching approach. Furthermore, the emission phantom measurement with and without EEG cap was performed at two different sessions. The counts recorded in a single detector element during two different PET measurements with identical measurement time and phantom radioactivity differ due the Poisson noise of the radioactive events. This noise propagates via the image reconstruction into single image pixels. Therefore, even if the mean values in a VOI are exactly zero between two PET studies, the standard deviation (SD) (and interquartile range (IQR)) will have relatively large values, as tabulated in S1 Table. Only a PET measurement carried out over a very long time, with its better count statistics and optimum reconstruction parameters will result in low SD. Thus, we conclude that even though there might be small artefacts due to the presence of the EEG cap in the emission images, these artefacts do not stand out from the general noise, which includes the factors discussed above. Furthermore, the results from the emission scans of the human subjects show that these additional attenuation effects, caused by the presence of the EEG electrodes, do not lead to any significant visible artefacts. Quantitative assessment on the relative difference images of the human subject data was not performed due to the fact that the metabolism of the FDG tracer changes during the acquisition time. It is well known and illustrated via the whole brain grey matter activity curve, shown in S3 Fig, that FDG accumulates in grey matter over time. Consequently, applying Eq 3 would result in images showing a negative activity distribution. These relative difference images can be observed in Fig 7 and S1 Fig. Even though a quantitative assessment of the human subject brain at different scan times is performed, as done for the phantom emission data, it is not possible to interpret this result because the difference image will show significant changes due to metabolic activity and accumulation of FDG. Factors, such as metabolic activity of the human brain, accumulation of FDG over time, and similar factors, as discussed for the phantom emission measurement in the above paragraph, contributed to the negative relative differences and the large standard deviation values tabulated in S2 Table.

In summary, there may be several factors responsible for the fact that the emission scan images of phantom and human subjects obtained in the PET-MR scanner using an EEG cap did not show any observable variations:

  • The EEG electrodes are housed in plastic made of polycarbonate and acrylate, which has attenuation coefficients similar to water (attenuation coefficient calculation is shown in S1 Appendix)—and thus brain tissue.
  • The attenuation factor along a line of response, which crosses an EEG electrode, results from the integral through all attenuating material along this line. The attenuation caused by the brain tissue and the bone along this line is much higher than that caused by the EEG electrodes, so their contribution can be regarded as negligible.
  • The diameter of the Ag/AgCl EEG electrode is 2 mm and the thickness varies between 1.1 and 1.4 mm. The total area of the 32 Ag/AgCl electrodes only covers about 0.1% of the total area of the head surface (circumference of average human head was considered as 58 cm [33] for this computation. Calculation of percentage of area covered by the electrodes on the surface of the human head is shown in S1 Appendix) so that the percentage of annihilation photons hitting the EEG electrodes on their path from the brain to the PET detectors crystals is very low. Thus, the influence of the lines of responses being attenuated by the EEG electrodes and used as input for image reconstruction can be regarded as negligible.
  • The electrode cables used in the EEG cap are very thin and most of the cable bundle is placed out of the field of view of PET detectors. Since the transmission study did not reveal any effect due to the cables, it follows that no noticeable effect should be expected in the emission scan.

In summary, we conclude that the attenuation due to EEG caps at the location of the electrodes is visible and quantifiable in ECAT HR+ transmission scan data; however, its influence is not visible in the reconstructed emission scans from BrainPET-MR scanner. The effects are inconsequential, given that the EEG cap components, such as plastic electrode housing and fabric material, have very low attenuation coefficients and the metal electrodes and chip resistors used are extremely thin.

Limitations

The results presented in this study are based on a 32 channel BrainCap MR EEG cap, however it may be of interest, especially with the background of source localisation, to study and compare the influence of a larger number of electrodes (64 or 128 channels) on PET emission data. Furthermore, the PET emission images of the human subjects were reconstructed without attenuation information about the EEG cap. The present ethical and radiation protection approvals do not include additional CT scans for healthy human subjects in order to create corresponding attenuation maps. Future studies will be planned and performed to investigate any mutual quantification bias with the attenuation maps containing the electrode information.

Supporting information

S1 Fig

FDG PET relative difference images in percentage between SUV images of human subjects with and without EEG cap.

The subjects 2 to 8 are shown.

(TIF)

S2 Fig

Histogram plots of the human emission relative difference images in percentage (only the whole grey matter region is considered in the plots.)

(TIF)

S3 Fig

Activity curves of all 8 subjects showing SUV in whole brain grey matter.

(TIF)

S1 Table

Statistical parameters calculated from VOIs drawn on the phantom emission relative difference image.

(DOCX)

S2 Table

Statistical parameters calculated from the whole grey matter region of the human emission relative difference images.

(DOCX)

S1 Appendix

Calculation of attenuation coefficient for polycarbonate and acrylate and calculation of percentage of area covered by the electrodes on the surface of the human head.

(PDF)

Acknowledgments

This study is part of the doctoral thesis (Dr. rer. medic.) of Ravichandran Rajkumar at the Medical Faculty of the RWTH Aachen University, Germany. We gratefully thank Dr. Jorge Arrubla for assistance in trimodal human data acquisition. We thank Andrea Muren, Cornelia Frey, Silke Frensch, and Suzanne Schaden for their technical assistance. We appreciate Mrs. Claire Rick for proof reading the manuscript. This study was in part supported by the EU FP7 funded project TRIMAGE (Nr. 602621).

Funding Statement

This trimodal imaging approach was supported by the EU-FP7 funded project TRIMAGE (Nr. 602621). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

According to our ethical approval, granted by the Ethics Committee of the Medical Faculty of RWTH Aachen University to the Forschungszentrum Juelich, data and results presented in this manuscript will be available upon request via Prof. Dr. Günther Schmalzing, Ethics Committee of the Medical Faculty, RWTH Aachen, Pauwelsstrasse 30, 52074 Aachen, Germany (ed.nehcaaku@gnizlamhcsg, ed.nehcaaku@nehcaake).

References

1. Shah NJ, Oros-Peusquens A-MM, Arrubla J, Zhang K, Warbrick T, Mauler J, et al. Advances in multimodal neuroimaging: Hybrid MR-PET and MR-PET-EEG at 3 T and 9.4 T. J Magn Reson. 2013;229: 101–15. doi: 10.1016/j.jmr.2012.11.027 [PubMed]
2. Neuner I, Mauler J, Arrubla J, Rota Kops E, Tellmann L, Scheins J, et al. Simultaneous trimodal MR-PET-EEG imaging for the investigation of resting state networks in humans. EJNMMI Physics. 2015;2: A71 doi: 10.1186/2197-7364-2-S1-A71 [PMC free article] [PubMed]
3. Neuner I, Wyss C, Rajkumar R, Mauler J, Boers F, Tellmann L, et al. Simultaneous trimodal imaging approach using [11C]—Flumazenil for MR-PET-EEG: feasibility of set-up. PSMR 2016—5th Conference on PET/MR and SPECT/MR; 2016 May 23–25; Cologne, Germany.
4. Wielaard J, Hammers A, Stirling J, Mackewn J, Mcginnity CJ, Schleyer P, et al. Simultaneous EEG-PET-fMRI/ASL recordings during a simulated seizure. International Conference on Basic and Clinical Multimodal Imaging, 2015, Sep 01–05; Utrecht, Netherlands.
5. Polikar R, Tilley C, Hillis B, Clark CM. Multimodal EEG, MRI and PET data fusion for Alzheimer’s disease diagnosis. 32nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, 2010 Aug 31–Sep 4, Buenos Aires. 2010; 6058–6061. 10.1109/IEMBS.2010.5627621 [PubMed]
6. Hampel H, Teipel SJ, Alexander GE, Pogarell O, Rapoport SI, Möller HJ. In vivo imaging of region and cell type specific neocortical neurodegeneration in Alzheimer’s disease: Perspectives of MRI derived corpus callosum measurement for mapping disease progression and effects of therapy. Evidence from studies with MRI, EEG and. J Neural Transm. 2002;109: 837–855. doi: 10.1007/s007020200069 [PubMed]
7. Laruelle M, Iyer RN, Al-Tikriti MS, Zea-Ponce Y, Malison R, Zoghbi SS, et al. Microdialysis and SPECT measurements of amphetamine-induced dopamine release in nonhuman primates. Synapse. 1997;25: 1–14. doi: 10.1002/(SICI)1098-2396(199701)25:1<1::AID-SYN1>3.0.CO;2-H [PubMed]
8. Gosseries O, Pistoia F, Charland-Verville V, Carolei A, Sacco S, Laureys S. The role of neuroimaging techniques in establishing diagnosis, prognosis and therapy in disorders of consciousness. Open Neuroimag J. 2016;10: 52–68. doi: 10.2174/1874440001610010052 [PMC free article] [PubMed]
9. Struck A, Hall L, Floberg J, Scott P, Douglas D. Surgical decision making in Temporal Lobe Epilepsy (TLE): a comparison of 18Fluorodeoxyglucose (FDG) Positron Emission Tomography (PET), MRI, and EEG. Epilepsy Behav. 2011;22: 293–297. doi: 10.1016/j.yebeh.2011.06.022 [PMC free article] [PubMed]
10. Uludağ K, Roebroeck A. General overview on the merits of multimodal neuroimaging data fusion. NeuroImage. 2014. pp. 3–10. doi: 10.1016/j.neuroimage.2014.05.018 [PubMed]
11. Liu S, Cai W, Liu S, Zhang F, Fulham M, Feng D, et al. Multimodal neuroimaging computing: a review of the applications in neuropsychiatric disorders. Brain Inform. Springer Berlin Heidelberg; 2015;2: 167–180. doi: 10.1007/s40708-015-0019-x [PMC free article] [PubMed]
12. O’Halloran R, Kopell BH, Sprooten E, Goodman WK, Frangou S. Multimodal neuroimaging-informed clinical applications in neuropsychiatric disorders. Front Psychiatry. 2016;7. [PMC free article] [PubMed]
13. Herzog H. PET / MRI: Challenges, solutions and perspectives. ZMedPhys. 2012;22: 281–298. [PubMed]
14. Weirich C, Brenner D, Scheins J, Besançon É, Tellmann L, Herzog H, et al. Analysis and correction of count rate reduction during simultaneous MR-PET measurements with the BrainPET scanner. IEEE Trans Med Imaging. 2012;31: 1372–1380. doi: 10.1109/TMI.2012.2188903 [PubMed]
15. Scheins JJ, Herzog H, Shah NJ. Fully-3D PET image reconstruction using scanner-independent, adaptive projection data and highly rotation-symmetric voxel assemblies. IEEE Trans Med Imaging. 2011;30: 879–892. doi: 10.1109/TMI.2011.2109732 [PubMed]
16. Tellmann L, Quick HH, Bockisch A, Herzog H, Beyer T. The effect of MR surface coils on PET quantification in whole-body PET / MR: Results from a pseudo-PET / MR phantom study. Med Phys. 2011;38: 2795–2805. doi: 10.1118/1.3582699 [PubMed]
17. Allen PJ, Polizzi G, Krakow K, Fish DR, Lemieux L. Identification of EEG events in the MR scanner: the problem of pulse artifact and a method for its subtraction. NeuroImage. 1998;8: 229–239. doi: 10.1006/nimg.1998.0361 [PubMed]
18. Allen PJ, Josephs O, Turner R. A method for removing imaging artifact from continuous EEG recorded during functional MRI. NeuroImage. 2000;12: 230–239. doi: 10.1006/nimg.2000.0599 [PubMed]
19. Neuner I, Arrubla J, Felder J, Shah NJ. Simultaneous EEG-fMRI acquisition at low, high and ultra-high magnetic fields up to 9.4T: Perspectives and challenges. NeuroImage. Elsevier B.V.; 2013; doi: 10.1016/j.neuroimage.2013.06.048 [PubMed]
20. Bonmassar G, Hadjikhani N, Ives JR, Hinton D, Belliveau JW. Influence of EEG Electrodes on the BOLD fMRI Signal. Hum Brain Mapp. 2001;14: 108–115. [PubMed]
21. Schreckenberger M, Lange-Asschenfeld C, Lochmann M, Mann K, Siessmeier T, Buchholz HG, et al. The thalamus as the generator and modulator of EEG alpha rhythm: A combined PET/EEG study with lorazepam challenge in humans. NeuroImage. 2004;22: 637–644. doi: 10.1016/j.neuroimage.2004.01.047 [PubMed]
22. Oakes TR, Pizzagalli DA, Hendrick AM, Horras KA, Larson CL, Abercrombie HC, et al. Functional coupling of simultaneous electrical and metabolic activity in the human brain. Hum Brain Mapp. 2004;21: 257–270. doi: 10.1002/hbm.20004 [PubMed]
23. Pizzagalli DA, Oakes TR, Davidson RJ. Coupling of theta activity and glucose metabolism in the human rostral anterior cingulate cortex: An EEG/PET study of normal and depressed subjects. Psychophysiology. 2003;40: 939–949. doi: 10.1111/1469-8986.00112 [PubMed]
24. Lindgren KA, Larson CL, Schaefer SM, Abercrombie HC, Ward RT, Oakes TR, et al. Thalamic metabolic rate predicts EEG alpha power in healthy control subjects but not in depressed patients. Biol Psychiatry. 1999;45: 943–952. doi: 10.1016/S0006-3223(98)00350-3 [PubMed]
25. Rota Kops E, Hautzel H, Herzog H, Antoch G, Shah NJ. Comparison template-based versus CT-based attenuation correction for hybrid MR/PET scanners. IEEE Trans Nucl Sci. 2015;62: 2115–2121. doi: 10.1109/TNS.2015.2452574
26. Lemmens C, Montandon M, Nuyts J. Impact of metal artefacts due to EEG electrodes in brain PET / CT imaging. Phys Med Biol. 2008;53: 4417–4429. doi: 10.1088/0031-9155/53/16/013 [PubMed]
27. Iida H, Hori Y, Ishida K, Imabayashi E, Matsuda H, Takahashi M, et al. Three-dimensional brain phantom containing bone and grey matter structures with a realistic head contour. Ann Nucl Med. 2013; 25–36. doi: 10.1007/s12149-012-0655-7 [PMC free article] [PubMed]
28. Brix G, Zaers J, Adam L-E, Bellemann ME, Ostertag H, Trojan H, et al. Performance evaluation of a whole-body PET scanner Using the NEMA protocol. J Nucl Med. 1997;Vol.38 No.: 1614–1623. [PubMed]
29. Schlemmer H-PW, Pichler BJ, Schmand M, Burbar Z, Michel C, Ladebeck R, et al. Simultaneous MR/PET imaging of the human brain: feasibility study. Radiology. 2008;248: 1028–1035. doi: 10.1148/radiol.2483071927 [PubMed]
30. Herzog H, Langen KJ, Weirich C, Rota Kops E, Kaffanke J, Tellmann L, et al. High resolution BrainPET combined with simultaneous MRI. Nuklearmedizin. 2011;50: 74–82. doi: 10.3413/Nukmed-0347-10-09 [PubMed]
31. Villien M, Wey H-Y, Mandeville JB, Catana C, Polimeni JR, Sander CY, et al. Dynamic functional imaging of brain glucose utilization using fPET-FDG. NeuroImage. Elsevier Inc.; 2014;100: 192–9. doi: 10.1016/j.neuroimage.2014.06.025 [PMC free article] [PubMed]
32. Smith SM. Fast robust automated brain extraction. Hum Brain Mapp. 2002;17: 143–155. doi: 10.1002/hbm.10062 [PubMed]
33. Bushby KM, Cole T, Matthews JN, Goodship JA. Centiles for adult head circumference. Arch Dis Child. 1992;67: 1286–1287. doi: 10.1136/adc.67.10.1286 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science