In vivo real-time positron emission particle tracking (PEPT) and single particle PET.

Positron emission particle tracking (PEPT) enables 3D localization and tracking of single positron-emitting radiolabelled particles with high spatiotemporal resolution. The translation of PEPT to the biomedical imaging field has been limited due to the lack of methods to radiolabel biocompatible particles with sufficient specific activity and protocols to isolate a single particle in the sub-micrometre size range, below the threshold for capillary embolization. Here we report two key developments: the synthesis and 68Ga-radiolabelling of homogeneous silica particles of 950 nm diameter with unprecedented specific activities (2.1 ± 1.4 kBq per particle), and the isolation and manipulation of a single particle. We have combined these developments to perform in vivo PEPT and dynamic positron emission tomography (PET) imaging of a single radiolabelled sub-micrometre size particle using a pre-clinical positron emission tomography/computed tomography scanner. This work opens possibilities for quantitative assessment of haemodynamics in vivo in real time, at the whole-body level using minimal amounts of injected radioactive dose and material.

Evaluation of manual and automated approaches for segmentation and extraction of quantitative indices from [18F]FDG PET-CT images.

Utilisation of whole organ volumes to extract anatomical and functional information from computed tomography (CT) and positron emission tomography (PET) images may provide key information for the treatment and follow-up of cancer patients. However, manual organ segmentation, is laborious and time-consuming. In this study, a CT-based deep learning method and a multi-atlas method were evaluated for segmenting the liver and spleen on CT images to extract quantitative tracer information from Fluorine-18 fluorodeoxyglucose ([18F]FDG) PET images of 50 patients with advanced Hodgkin lymphoma (HL). Manual segmentation was used as the reference method. The two automatic methods were also compared with a manually defined volume of interest (VOI) within the organ, a technique commonly performed in clinical settings. Both automatic methods provided accurate CT segmentations, with the deep learning method outperforming the multi-atlas with a DICE coefficient of 0.93 ± 0.03 (mean ± standard deviation) in liver and 0.87 ± 0.17 in spleen compared to 0.87 ± 0.05 (liver) and 0.78 ± 0.11 (spleen) for the multi-atlas. Similarly, a mean relative error of -3.2% for the liver and -3.4% for the spleen across patients was found for the mean standardized uptake value (SUVmean) using the deep learning regions while the corresponding errors for the multi-atlas method were -4.7% and -9.2%, respectively. For the maximum SUV (SUVmax), both methods resulted in higher than 20% overestimation due to the extension of organ boundaries to include neighbouring, high-uptake regions. The conservative VOI method which did not extend into neighbouring tissues, provided a more accurate SUVmaxestimate. In conclusion, the automatic, and particularly the deep learning method could be used to rapidly extract information of the SUVmeanwithin the liver and spleen. However, activity from neighbouring organs and lesions can lead to high biases in SUVmaxand current practices of manually defining a volume of interest in the organ should be considered instead.

Improved 3D tumour definition and quantification of uptake in simulated lung tumours using deep learning.

Objective.In clinical positron emission tomography (PET) imaging, quantification of radiotracer uptake in tumours is often performed using semi-quantitative measurements such as the standardised uptake value (SUV). For small objects, the accuracy of SUV estimates is limited by the noise properties of PET images and the partial volume effect. There is need for methods that provide more accurate and reproducible quantification of radiotracer uptake.Approach.In this work, we present a deep learning approach with the aim of improving quantification of lung tumour radiotracer uptake and tumour shape definition. A set of simulated tumours, assigned with 'ground truth' radiotracer distributions, are used to generate realistic PET raw data which are then reconstructed into PET images. In this work, the ground truth images are generated by placing simulated tumours characterised by different sizes and activity distributions in the left lung of an anthropomorphic phantom. These images are then used as input to an analytical simulator to simulate realistic raw PET data. The PET images reconstructed from the simulated raw data and the corresponding ground truth images are used to train a 3D convolutional neural network.Results.When tested on an unseen set of reconstructed PET phantom images, the network yields improved estimates of the corresponding ground truth. The same network is then applied to reconstructed PET data generated with different point spread functions. Overall the network is able to recover better defined tumour shapes and improved estimates of tumour maximum and median activities.Significance.Our results suggest that the proposed approach, trained on data simulated with one scanner geometry, has the potential to restore PET data acquired with different scanners.

Production of copper-64 using a hospital cyclotron: targetry, purification and quality analysis.

OBJECTIVES: To construct and evaluate a 64Cu production system that minimises the amount of costly 64Ni, radionuclidic impurities and nonradioactive metal contamination and maximises radiochemical and radionuclidic purity and molar activity; and to report analytical and quality control methods that can be used within typical PET radiochemistry production facilities to measure metal ion concentrations and radiometal molar activities. METHODS: Low volume was ensured by dissolving the irradiated nickel in a low volume of hydrochloric acid (<1 mL) using the concave gold target backing as a reaction vessel in a custom-built target holder. Removal of contaminating 55Co and nonradioactive trace metals was ensured by adding an intermediate hydrochloric acid concentration step during the conventional ion-exchange elution process. The radionuclidic purity of the product was determined by half-life measurements, gamma spectroscopy and ion radiochromatography. Trace metal contamination and molar activity were determined by ion chromatography. RESULTS AND CONCLUSIONS: On a small scale, suitable for preclinical research, the process produced typically 3.2 GBq 64Cu in 2 mL solution from 9.4 ± 2.1 mg nickel-64 electroplated onto a gold target backing. The product had high molar activity (121.5 GBq/µmol), was free of trace metal contamination detectable by ion chromatography and has been used for many preclinical and clinical PET imaging applications.

Standardisation of conventional and advanced iterative reconstruction methods for Gallium-68 multi-centre PET-CT trials.

PURPOSE: To assess the applicability of the Fluorine-18 performance specifications defined by EANM Research Ltd (EARL), in Gallium-68 multi-centre PET-CT trials using conventional (ordered subset expectation maximisation, OSEM) and advanced iterative reconstructions which include the systems' point spread function (PSF) and a Bayesian penalised likelihood algorithm (BPL) commercially known as Q.CLEAR. The possibility of standardising the two advanced reconstruction methods was examined. METHODS: The NEMA image quality phantom was filled with Gallium-68 and scanned on a GE PET-CT system. PSF and BPL with varying post-reconstruction Gaussian filter width (2-6.4 mm) and penalisation factor (200-1200), respectively, were applied. The average peak-to-valley ratio from six profiles across each sphere was estimated to inspect any edge artefacts. Image noise was assessed using background variability and image roughness. Six GE and Siemens PET-CT scanners provided Gallium-68 images of the NEMA phantom using both conventional and advanced reconstructions from which the maximum, mean and peak recoveries were drawn. Fourteen patients underwent 68Ga-PSMA PET-CT imaging. BPL (200-1200) reconstructions of the data were compared against PSF smoothed with a 6.4-mm Gaussian filter. RESULTS: A Gaussian filter width of approximately 6 mm for PSF and a penalisation factor of 800 for BPL were needed to suppress the edge artefacts. In addition, those reconstructions provided the closest agreement between the two advanced iterative reconstructions and low noise levels with the background variability and the image roughness being lower than 7.5% and 11.5%, respectively. The recoveries for all methods generally performed at the lower limits of the EARL specifications, especially for the 13- and 10-mm spheres for which up to 27% (conventional) and 41% (advanced reconstructions) lower limits are suggested. The lesion standardised uptake values from the clinical data were significantly different between BPL and PSF smoothed with a Gaussian filter of 6.4 mm wide for all penalisation factors except for 800 and 1000. CONCLUSION: It is possible to standardise the advanced reconstruction methods with the reconstruction parameters being also sufficient for minimising the edge artefacts and noise in the images. For both conventional and advanced reconstructions, Gallium-68 specific recovery coefficient limits were required, especially for the smallest phantom spheres.

Evaluation of [13N]ammonia positron emission tomography as a potential method for quantifying glutamine synthetase activity in the human brain.

PURPOSE: The conversion of synaptic glutamate to glutamine in astrocytes by glutamine synthetase (GS) is critical to maintaining healthy brain activity and may be disrupted in several brain disorders. As the GS catalysed conversion of glutamate to glutamine requires ammonia, we evaluated whether [13N]ammonia positron emission tomography (PET) could reliability quantify GS activity in humans. METHODS: In this test-retest study, eight healthy volunteers each received two dynamic [13N]ammonia PET scans on the morning and afternoon of the same day. Each [13N]ammonia scan was preceded by a [15O]water PET scan to account for effects of cerebral blood flow (CBF). RESULTS: Concentrations of radioactive metabolites in arterial blood were available for both sessions in five of the eight subjects. Our results demonstrated that kinetic modelling was unable to reliably distinguish estimates of the kinetic rate constant k3 (related to GS activity) from K1 (related to [13N]ammonia brain uptake), and indicated a non-negligible back-flux of [13N] to blood (k2). Model selection favoured a reversible one-tissue compartmental model, and [13N]ammonia K1 correlated reliably (r2 = 0.72-0.92) with [15O]water CBF. CONCLUSION: The [13N]ammonia PET method was unable to reliably estimate GS activity in the human brain but may provide an alternative index of CBF.

Weighted Manifold Alignment using Wave Kernel Signatures for Aligning Medical Image Datasets.

Manifold alignment (MA) is a technique to map many high-dimensional datasets to one shared low-dimensional space. Here we develop a pipeline for using MA to reconstruct high-resolution medical images. We present two key contributions. First, we develop a novel MA scheme in which each high-dimensional dataset can be differently weighted preventing noisier or less informative data from corrupting the aligned embedding. We find that this generalisation improves performance in our experiments in both supervised and unsupervised MA problems. Second, we use the wave kernel signature as a graph descriptor for the unsupervised MA case finding that it significantly outperforms the current state-of-the-art methods and provides higher quality reconstructed magnetic resonance volumes than existing methods.

Assessing the feasibility of intranasal radiotracer administration for in brain PET imaging.

INTRODUCTION: The development of clinically useful tracers for PET imaging is enormously challenging and expensive. The intranasal (IN) route of administration is purported to be a viable route for delivering drugs to the brain but has, as yet, not been investigated for the delivery of PET tracers. If the intranasal (IN) pathway presents a viable option, it extends the PET imaging field by increasing the number of tracers available for human use. Here we report the results of a rodent study testing the feasibility of the IN route to administer radiotracers for brain PET imaging. METHODS: We used two different, well characterised, brain penetrant radiotracers, [18F]fluorodeoxyglucose ([18F]FDG) and [18F]fallypride, and aimed to evaluate the pharmacokinetics after administration of the tracers via the intranasal route, and contrast this to intravenous administration. Image acquisition was carried out after tracer administration and arterial blood samples were collected at different time intervals, centrifuged to extract plasma and gamma counted. We hypothesised that [brain region]:[plasma] ratios would be higher via the intranasal route as there are two inputs, one directly from the nose to the brain, and another from the peripheral circulation. To assess the feasibility of using this approach clinically, we used these data to estimate radiation dosimetry in humans. RESULTS: Contrary to our hypothesis, in case of both radiotracers, we did not see a higher ratio in the expected brain regions, except in the olfactory bulb, that is closest to the nose. It appears that the radiotracers move into the olfactory bulb region, but then do not progress further into other brain regions. Moreover, as the nasal cavity has a small surface area, the extrapolated dosimetry estimations for intranasal human imaging showed an unacceptably high level (15 mSv/MBq) of cumulative skin radiation exposure. CONCLUSIONS: Therefore, although an attractive route for brain permeation, we conclude that the intranasal route would present difficulties due to non-specific signal and radiation dosimetry considerations for brain PET imaging.

The impact of hypoglycaemia awareness status on regional brain responses to acute hypoglycaemia in men with type 1 diabetes.

AIMS/HYPOTHESIS: Impaired awareness of hypoglycaemia (IAH) in type 1 diabetes increases the risk of severe hypoglycaemia sixfold and can be resistant to intervention. We explored the impact of IAH on central responses to hypoglycaemia to investigate the mechanisms underlying barriers to therapeutic intervention. METHODS: We conducted [15O]water positron emission tomography studies of regional brain perfusion during euglycaemia (target 5 mmol/l), hypoglycaemia (achieved level, 2.4 mmol/l) and recovery (target 5 mmol/l) in 17 men with type 1 diabetes: eight with IAH, and nine with intact hypoglycaemia awareness (HA). RESULTS: Hypoglycaemia with HA was associated with increased activation in brain regions including the thalamus, insula, globus pallidus (GP), anterior cingulate cortex (ACC), orbital cortex, dorsolateral frontal (DLF) cortex, angular gyrus and amygdala; deactivation occurred in the temporal and parahippocampal regions. IAH was associated with reduced catecholamine and symptom responses to hypoglycaemia vs HA (incremental AUC: autonomic scores, 26.2 ± 35.5 vs 422.7 ± 237.1; neuroglycopenic scores, 34.8 ± 88.8 vs 478.9 ± 311.1; both p < 0.002). There were subtle differences (p < 0.005, k ≥ 50 voxels) in brain activation at hypoglycaemia, including early differences in the right central operculum, bilateral medial orbital (MO) cortex, and left posterior DLF cortex, with additional differences in the ACC, right GP and post- and pre-central gyri in established hypoglycaemia, and lack of deactivation in temporal regions in established hypoglycaemia. CONCLUSIONS/INTERPRETATION: Differences in activation in the post- and pre-central gyri may be expected in people with reduced subjective responses to hypoglycaemia. Alterations in the activity of regions involved in the drive to eat (operculum), emotional salience (MO cortex), aversion (GP) and recall (temporal) suggest differences in the perceived importance and urgency of responses to hypoglycaemia in IAH compared with HA, which may be key to the persistence of the condition.

The effect of high count rates on cardiac perfusion quantification in a simultaneous PET-MR system using a cardiac perfusion phantom.

BACKGROUND: PET-MRI is under investigation as a new strategy for quantitative myocardial perfusion imaging. Consideration is required as to the maximum scanner count rate in order to limit dead-time losses resulting from administered activity in the scanner field of view during the first pass of the radiotracer. RESULTS: We performed a decaying-source experiment to investigate the high count-rate performance of a PET-MR system (Siemens mMR) over the expected range of activities during a clinical study. We also performed imaging of a cardiac perfusion phantom, which provides an experimental simulation of clinical transit of a simultaneous radiotracer (phantom injected activities range 252 to 997 MBq) and gadolinium-based contrast agent (GBCA). Time-activity and time-intensity curves of the aorta and myocardium compartments from PET and MR images were determined, and quantification of perfusion was then performed using a standard cardiac kinetic model. The decaying-source experiment showed a maximum noise equivalent count rate (NECRmax) of 286 kcps at a singles rate of 47.1 Mcps. NECR was maintained within 5% (NECR95%) of the NECRmax with a singles rate of 34.1 Mcps, corresponding to 310 MBq in the phantom. Count-rate performance was degraded above the singles rate of 64.9 Mcps due to the number of detection events impacting the quantitative accuracy of reconstructed images. A 10% bias in image activity concentration was observed between singles rates of 78.2 and 82.9 Mcps. Perfusion phantom experiments showed that image-based activity concentration and quantified values of perfusion were affected by count losses when the total singles rate was greater than 64.9 Mcps. This occurred during the peak arterial input function (AIF) phase of imaging for injected activities to the phantom of 600 MBq and greater. CONCLUSIONS: Care should be taken to avoid high count-rate losses in simultaneous PET-MRI studies. Based on our results in phantoms, bias in reconstructed images should be avoided by adhering to a singles rate lower than 64.9 Mcps on the mMR system. Quantification of perfusion values using singles rates higher than 64.9 Mcps on this system may be compromised and should be avoided.

18F-Tetrafluoroborate, a PET Probe for Imaging Sodium/Iodide Symporter Expression: Whole-Body Biodistribution, Safety, and Radiation Dosimetry in Thyroid Cancer Patients.

We report the safety, biodistribution, and internal radiation dosimetry, in humans with thyroid cancer, of 18F-tetrafluoroborate (18F-TFB), a novel PET radioligand for imaging the human sodium/iodide symporter (hNIS). Methods: Serial whole-body PET scans of 5 subjects with recently diagnosed thyroid cancer were acquired before surgery for up to 4 h after injection of 184 ± 15 MBq of 18F-TFB. Activity was determined in whole blood, plasma, and urine. Mean organ-absorbed doses and effective doses were calculated via quantitative image analysis and using OLINDA/EXM software. Results: Images showed a high uptake of 18F-TFB in known areas of high hNIS expression (thyroid, salivary glands, and stomach). Excretion was predominantly renal. No adverse effects in relation to safety of the radiopharmaceutical were observed. The effective dose was 0.0326 ± 0.0018 mSv/MBq. The critical tissues/organs receiving the highest mean sex-averaged absorbed doses were the thyroid (0.135 ± 0.079 mSv/MBq), stomach (0.069 ± 0.022 mSv/MBq), and salivary glands (parotids, 0.031 ± 0.011 mSv/MBq; submandibular, 0.061 ± 0.031 mSv/MBq). Other organs of interest were the bladder (0.102 ± 0.046 mSv/MBq) and kidneys (0.029 ± 0.009 mSv/MBq). Conclusion: Imaging using 18F-TFB imparts a radiation exposure similar in magnitude to many other 18F-labeled radiotracers. 18F-TFB shows a biodistribution similar to 99mTc-pertechnetate, a known nonorganified hNIS tracer, and is pharmacologically and radiobiologically safe in humans. Phase 2 trials for 18F-TFB as an hNIS imaging agent are warranted.

High-Resolution Self-Gated Dynamic Abdominal MRI Using Manifold Alignment.

We present a novel retrospective self-gating method based on manifold alignment (MA), which enables reconstruction of free breathing, high spatial, and temporal resolution abdominal magnetic resonance imaging sequences. Based on a radial golden-angle acquisition trajectory, our method enables a multidimensional self-gating signal to be extracted from the k -space data for more accurate motion representation. The k -space radial profiles are evenly divided into a number of overlapping groups based on their radial angles. MA is then used to simultaneously learn and align the low dimensional manifolds of all groups, and embed them into a common manifold. In the manifold, k -space profiles that represent similar respiratory positions are close to each other. Image reconstruction is performed by combining radial profiles with evenly distributed angles that are close in the manifold. Our method was evaluated on both 2-D and 3-D synthetic and in vivo data sets. On the synthetic data sets, our method achieved high correlation with the ground truth in terms of image intensity and virtual navigator values. Using the in vivo data, compared with a state-of-the-art approach based on the center of k -space gating, our method was able to make use of much richer profile data for self-gating, resulting in statistically significantly better quantitative measurements in terms of organ sharpness and image gradient entropy.

Characterization of attenuation and respiratory motion artifacts and their influence on SPECT MP image evaluation using a dynamic phantom assembly with variable cardiac defects.

BACKGROUND: A phantom assembly that simulates the respiratory motion of the heart was used to investigate artifacts and their impact on defect detection. METHODS: SPECT/CT images were acquired for phantoms with and without small and large cardiac defects during normal and deep breathing, and also at four static respiratory phases. Acquisitions were reconstructed with and without AC, and with misalignment of transmission and emission scans. A quantitative analysis was performed to assess artifacts. Two physicians reported on defect presence or absence and their results were evaluated. RESULTS: All large defects were correctly reported. Attenuation reduced uptake in the basal LV walls, increasing FN physicians' reports for small defects. Respiratory motion reduced uptake mainly in the anterior and inferior walls increasing FP and FN reports on images without and with small defects, respectively. Artifacts due to misalignment between CT and SPECT scans in normal breathing phantoms did not influence the physicians' reports. CONCLUSIONS: Attenuation and respiratory motion correction should be applied to reduce artifacts before reporting on small defects in deep breathing conditions. Artifacts due to misalignment between CT and SPECT scans do not affect defect detection in normal breathing when the LV is co-registered in SPECT and CT images prior to AC.

Impact of respiratory motion correction on SPECT myocardial perfusion imaging using a mechanically moving phantom assembly with variable cardiac defects.

BACKGROUND: The aim of this study was to determine the impact of respiratory motion correction on SPECT MPI and on defect detection using a phantom assembly. METHODS: SPECT/CT data were acquired using an anthropomorphic phantom with inflatable lungs and with an ECG beating and moving cardiac compartment. The heart motion followed the respiratory pattern in the cranio-caudal direction to simulate normal or deep breathing. Small or large transmural defects were inserted into the myocardial wall of the left ventricle. SPECT/CT images were acquired for each of the four respiratory phases, from exhale to inhale. A respiratory motion correction was applied using an image-based method with transformation parameters derived from the SPECT data by a non-rigid registration algorithm. A report on defect detection from two physicians and a quantitative analysis on MPI data were performed before and after applying motion correction. RESULTS: Respiratory motion correction eliminated artifacts present in the images, resulting in a uniform uptake and reduction of motion blurring, especially in the inferior and anterior regions of the LV myocardial walls. The physicians' report after motion correction showed that images were corrected for motion. CONCLUSIONS: A combination of motion correction with attenuation correction reduces artifacts in SPECT MPI. AC-SPECT images with and without motion correction should be simultaneously inspected to report on small defects.

Feasibility of simultaneous PET-MR perfusion using a novel cardiac perfusion phantom.

BACKGROUND: PET-MR scanners are beginning to be employed for quantitative myocardial perfusion imaging. In order to examine simultaneous perfusion calculations, this work describes a feasibility study of simultaneous PET-MR of gadolinium-based contrast agent (GBCA) and PET radiotracer in a novel cardiac perfusion phantom. RESULTS: [18F]F- and GBCA were injected simultaneously into a cardiac phantom using a range of ground-truth myocardial perfusion rates of 1 to 5 ml/g/min. PET quantification of K1 (ml/g/min) was performed using a single tissue compartment model. MR perfusion was calculated using a model-independent signal deconvolution technique. PET and MR signal traces from the phantom aorta and myocardial sections show true simultaneous PET and MR arterial input functions (AIF) and myocardial uptake respectively at each perfusion rate. Calculation of perfusion parameters showed both K1 and h(t = 0) (PET and MR perfusion parameters respectively) to be linearly related with the ground truth perfusion rate (PT ), and also linearly related to each other (R2 = 0.99). The highest difference in perfusion values between K1 and PT was 16% at 1 ml/g/min, and the mean difference for all other perfusion rates was <3%. CONCLUSIONS: The perfusion phantom allows accurate and reproducible simulation of the myocardial kinetics for simultaneous PET-MR imaging, and may find use in protocol design and development of PET-MR based quantification techniques and direct comparison of quantification of the two modalities.

Differences in Regional Brain Responses to Food Ingestion After Roux-en-Y Gastric Bypass and the Role of Gut Peptides: A Neuroimaging Study.

OBJECTIVE: Improved appetite control, possibly mediated by exaggerated gut peptide responses to eating, may contribute to weight loss after Roux-en-Y gastric bypass (RYGB). This study compared brain responses to food ingestion between post-RYGB (RYGB), normal weight (NW), and obese (Ob) unoperated subjects and explored the role of gut peptide responses in RYGB. RESEARCH DESIGN AND METHODS: Neuroimaging with [(18)F]-fluorodeoxyglucose (FDG) positron emission tomography was performed in 12 NW, 21 Ob, and 9 RYGB (18 ± 13 months postsurgery) subjects after an overnight fast, once FED (400 kcal mixed meal), and once FASTED, in random order. RYGB subjects repeated the studies with somatostatin infusion and basal insulin replacement. Fullness, sickness, and postscan ad libitum meal consumption were measured. Regional brain FDG uptake was compared using statistical parametric mapping. RESULTS: RYGB subjects had higher overall fullness and food-induced sickness and lower ad libitum consumption. Brain responses to eating differed in the hypothalamus and pituitary (exaggerated activation in RYGB), left medial orbital cortex (OC) (activation in RYGB, deactivation in NW), right dorsolateral frontal cortex (deactivation in RYGB and NW, absent in Ob), and regions mapping to the default mode network (exaggerated deactivation in RYGB). Somatostatin in RYGB reduced postprandial gut peptide responses, sickness, and medial OC activation. CONCLUSIONS: RYGB induces weight loss by augmenting normal brain responses to eating in energy balance regions, restoring lost inhibitory control, and altering hedonic responses. Altered postprandial gut peptide responses primarily mediate changes in food-induced sickness and OC responses, likely to associate with food avoidance.

Development of an MRI-compatible digital SiPM detector stack for simultaneous PET/MRI.

Advances in solid-state photon detectors paved the way to combine positron emission tomography (PET) and magnetic resonance imaging (MRI) into highly integrated, truly simultaneous, hybrid imaging systems. Based on the most recent digital SiPM technology, we developed an MRI-compatible PET detector stack, intended as a building block for next generation simultaneous PET/MRI systems. Our detector stack comprises an array of 8 × 8 digital SiPM channels with 4 mm pitch using Philips Digital Photon Counting DPC 3200-22 devices, an FPGA for data acquisition, a supply voltage control system and a cooling infrastructure. This is the first detector design that allows the operation of digital SiPMs simultaneously inside an MRI system. We tested and optimized the MRI-compatibility of our detector stack on a laboratory test bench as well as in combination with a Philips Achieva 3 T MRI system. Our design clearly reduces distortions of the static magnetic field compared to a conventional design. The MRI static magnetic field causes weak and directional drift effects on voltage regulators, but has no direct impact on detector performance. MRI gradient switching initially degraded energy and timing resolution. Both distortions could be ascribed to voltage variations induced on the bias and the FPGA core voltage supply respectively. Based on these findings, we improved our detector design and our final design shows virtually no energy or timing degradations, even during heavy and continuous MRI gradient switching. In particular, we found no evidence that the performance of the DPC 3200-22 digital SiPM itself is degraded by the MRI system.

A method for accurate spatial registration of PET images and histopathology slices.

BACKGROUND: Accurate alignment between histopathology slices and positron emission tomography (PET) images is important for radiopharmaceutical validation studies. Limited data is available on the registration accuracy that can be achieved between PET and histopathology slices acquired under routine pathology conditions where slices may be non-parallel, non-contiguously cut and of standard block size. The purpose of this study was to demonstrate a method for aligning PET images and histopathology slices acquired from patients with laryngeal cancer and to assess the registration accuracy obtained under these conditions. METHODS: Six subjects with laryngeal cancer underwent a (64)Cu-copper-II-diacetyl-bis(N4-methylthiosemicarbazone) ((64)Cu-ATSM) PET computed tomography (CT) scan prior to total laryngectomy. Sea urchin spines were inserted into the pathology specimen to act as fiducial markers. The specimen was fixed in formalin, as per standard histopathology operating procedures, and was then CT scanned and cut into millimetre-thick tissue slices. A subset of the tissue slices that included both tumour and fiducial markers was taken and embedded in paraffin blocks. Subsequently, microtome sectioning and haematoxylin and eosin staining were performed to produce 5-µm-thick tissue sections for microscopic digitisation. A series of rigid registration procedures was performed between the different imaging modalities (PET; in vivo CT-i.e. the CT component of the PET-CT; ex vivo CT; histology slices) with the ex vivo CT serving as the reference image. In vivo and ex vivo CTs were registered using landmark-based registration. Histopathology and ex vivo CT images were aligned using the sea urchin spines with additional anatomical landmarks where available. Registration errors were estimated using a leave-one-out strategy for in vivo to ex vivo CT and were estimated from the RMS landmark accuracy for histopathology to ex vivo CT. RESULTS: The mean ± SD accuracy for registration of the in vivo to ex vivo CT images was 2.66 ± 0.66 mm, and the accuracy for registration of histopathology to ex vivo CT was 0.86 ± 0.41 mm. Estimating the PET to in vivo CT registration accuracy to equal the PET-CT alignment accuracy of 1 mm resulted in an overall average registration error between PET and histopathology slices of 3.0 ± 0.7 mm. CONCLUSIONS: We have developed a registration method to align PET images and histopathology slices with an accuracy comparable to the spatial resolution of the PET images.

Respiratory motion correction of PET using MR-constrained PET-PET registration.

BACKGROUND: Respiratory motion in positron emission tomography (PET) is an unavoidable source of error in the measurement of tracer uptake, lesion position and lesion size. The introduction of PET-MR dual modality scanners opens a new avenue for addressing this issue. Motion models offer a way to estimate motion using a reduced number of parameters. This can be beneficial for estimating motion from PET, which can otherwise be difficult due to the high level of noise of the data. METHOD: We propose a novel technique that makes use of a respiratory motion model, formed from initial MR scan data. The motion model is used to constrain PET-PET registrations between a reference PET gate and the gates to be corrected. For evaluation, PET with added FDG-avid lesions was simulated from real, segmented, ultrashort echo time MR data obtained from four volunteers. Respiratory motion was included in the simulations using motion fields derived from real dynamic 3D MR volumes obtained from the same volunteers. RESULTS: Performance was compared to an MR-derived motion model driven method (which requires constant use of the MR scanner) and to unconstrained PET-PET registration of the PET gates. Without motion correction, a median drop in uncorrected lesion [Formula: see text] intensity to [Formula: see text] and an increase in median head-foot lesion width, specified by a minimum bounding box, to [Formula: see text] was observed relative to the corresponding measures in motion-free simulations. The proposed method corrected these values to [Formula: see text] ([Formula: see text]) and [Formula: see text] ([Formula: see text]) respectively, with notably improved performance close to the diaphragm and in the liver. Median lesion displacement across all lesions was observed to be [Formula: see text] without motion correction, which was reduced to [Formula: see text] ([Formula: see text]) with motion correction. DISCUSSION: This paper presents a novel technique for respiratory motion correction of PET data in PET-MR imaging. After an initial 30 second MR scan, the proposed technique does not require use of the MR scanner for motion correction purposes, making it suitable for MR-intensive studies or sequential PET-MR. The accuracy of the proposed technique was similar to both comparative methods, but robustness was improved compared to the PET-PET technique, particularly in regions with higher noise such as the liver.

False Discovery Rates in PET and CT Studies with Texture Features: A Systematic Review.

PURPOSE: A number of recent publications have proposed that a family of image-derived indices, called texture features, can predict clinical outcome in patients with cancer. However, the investigation of multiple indices on a single data set can lead to significant inflation of type-I errors. We report a systematic review of the type-I error inflation in such studies and review the evidence regarding associations between patient outcome and texture features derived from positron emission tomography (PET) or computed tomography (CT) images. METHODS: For study identification PubMed and Scopus were searched (1/2000-9/2013) using combinations of the keywords texture, prognostic, predictive and cancer. Studies were divided into three categories according to the sources of the type-I error inflation and the use or not of an independent validation dataset. For each study, the true type-I error probability and the adjusted level of significance were estimated using the optimum cut-off approach correction, and the Benjamini-Hochberg method. To demonstrate explicitly the variable selection bias in these studies, we re-analyzed data from one of the published studies, but using 100 random variables substituted for the original image-derived indices. The significance of the random variables as potential predictors of outcome was examined using the analysis methods used in the identified studies. RESULTS: Fifteen studies were identified. After applying appropriate statistical corrections, an average type-I error probability of 76% (range: 34-99%) was estimated with the majority of published results not reaching statistical significance. Only 3/15 studies used a validation dataset. For the 100 random variables examined, 10% proved to be significant predictors of survival when subjected to ROC and multiple hypothesis testing analysis. CONCLUSIONS: We found insufficient evidence to support a relationship between PET or CT texture features and patient survival. Further fit for purpose validation of these image-derived biomarkers should be supported by appropriate biological and statistical evidence before their association with patient outcome is investigated in prospective studies.