Longitudinal progression of sporadic Parkinson's disease: a multi-tracer positron emission tomography study.

Parkinson's disease is a heterogeneous disorder with multiple factors contributing to disease initiation and progression. Using serial, multi-tracer positron emission tomography imaging, we studied a cohort of 78 subjects with sporadic Parkinson's disease to understand the disease course better. Subjects were scanned with radiotracers of presynaptic dopaminergic integrity at baseline and again after 4 and 8 years of follow-up. Non-linear multivariate regression analyses, using random effects, of the form BP(ND)(t) or K(occ)(t) = a*e((-)(bt)(-d)(A) + c, where BP(ND) = tracer binding potential (nondispaceable), K(OCC) = tracer uptake constant a, b, c and d are regression parameters, t is the symptom duration and A is the age at onset, were utilized to model the longitudinal progression of radiotracer binding/uptake. We found that the initial tracer binding/uptake was significantly different in anterior versus posterior striatal subregions, indicating that the degree of denervation at disease onset was different between regions. However, the relative rate of decline in tracer binding/uptake was similar between the striatal subregions. While an antero-posterior gradient of severity was maintained for dopamine synthesis, storage and reuptake, the asymmetry between the more and less affected striatum became less prominent over the disease course. Our study suggests that the mechanisms underlying Parkinson's disease initiation and progression are probably different. Whereas factors responsible for disease initiation affect striatal subregions differently, those factors contributing to disease progression affect all striatal subregions to a similar degree and may therefore reflect non-specific mechanisms such as oxidative stress, inflammation or excitotoxicity.

Expectation and dopamine release: mechanism of the placebo effect in Parkinson's disease.

The power of placebos has long been recognized for improving numerous medical conditions such as Parkinson's disease (PD). Little is known, however, about the mechanism underlying the placebo effect. Using the ability of endogenous dopamine to compete for [11C]raclopride binding as measured by positron emission tomography, we provide in vivo evidence for substantial release of endogenous dopamine in the striatum of PD patients in response to placebo. Our findings indicate that the placebo effect in PD is powerful and is mediated through activation of the damaged nigrostriatal dopamine system.

Quantitative SPECT imaging and biodistribution point to molecular weight independent tumor uptake for some long-circulating polymer nanocarriers.

Polymeric nanocarriers are promising entities for cancer diagnosis and therapy. The aim of such nanocarriers is to selectively accumulate in cancerous tissue that is difficult to visualize or treat. The passive accumulation of a nanocarrier in a tumor through extravasation is often attributed to the enhanced permeation and retention (EPR) effect and the size and shape of the nanocarrier. However, the tumor microenvironment is very heterogeneous and the intratumoral pressure is usually high, leading to different opinions about how the EPR of nanocarriers through the irregular vasculature of a tumor leads to accumulation. In order to investigate this topic, we studied methods for the determination of pharmacokinetic parameters, biodistribution and the tumor uptake of nanocarriers. More specifically, we used non-invasive quantitative Single-Photon Emission Computed Tomography/Computed Tomography (qSPECT/CT) imaging of hyperbranched polyglycerols (HPGs) to explore the specific biodistribution and tumor uptake of six model nanocarriers in Rag2m mice. We were interested to see if a distinct molecular weight (MW) of nanocarriers (HPG 25, 50, 100, 200, 300, 500 kDa) is favoured by the tumor. To trace the model nanocarriers, HPGs were covalently linked to the strong chelator desferrioxamine (DFO), and radiolabeled with the gamma emitter 67Ga (EC = 100%, E γ = 185 keV (21.4%), 300 keV (16.6%), half-life = 3.26 d). Without the need for blood collection, but instead using qSPECT/CT imaging inside the heart, the blood circulation half-lives of the 67Ga labeled HPGs were determined and increased from 9.9 ± 2.9 to 47.8 ± 7.9 hours with increasing polymer MW. Total tumor accumulation correlated positively with the circulation time of the HPGs. Comparing the tumor-to-blood ratio dynamically revealed how blood and tumor concentrations of the nanocarrier change over time and when equilibrium is reached. The time of equilibrium is size-dependent and increases with molecular weight. Furthermore, the data indicate that for larger MWs, nanocarrier uptake and retention by the tumor is size independent. Further studies are necessary to advance our understanding of the interplay between MW and nanoparticle accumulation in tumors.

Evaluation of 209At as a theranostic isotope for 209At-radiopharmaceutical development using high-energy SPECT.

The development of alpha-emitting radiopharmaceuticals using 211At requires quantitative determination of the time-dependent nature of the 211At biodistribution. However, imaging-based methods for acquiring this information with 211At have not found wide-spread use because of its low abundance of decay emissions suitable for external detection. In this publication we demonstrate the theranostic abilities of the 211At/209At isotope pair and present the first-ever 209At SPECT images. The VECTor microSPECT/PET/CT scanner was used to image 209At with a collimator suitable for the 511 keV annihilation photons of PET isotopes. Data from distinct photopeaks of the 209At energy spectrum (195 keV (22.6%), 239 keV (12.4 %), 545 keV (91.0 %), a combined 782/790 keV peak (147 %), and 209Po x-rays (139.0 %)) were independently evaluated for use in image reconstructions using Monte Carlo (GATE) simulations and phantom studies. 209At-imaging in vivo was demonstrated in a healthy mouse injected with 10 MBq of free [209At]astatide. Image-based measurements of 209At uptake in organs of interest-acquired in 5 min intervals-were compared to ex vivo gamma counter measurements of the same organs. Simulated and measured data indicated that-due to the large amount of scatter from high energy (>750 keV) gammas-reconstructed images using the x-ray peak outperformed those obtained from other peaks in terms of image uniformity and spatial resolution, determined to be <0.85 mm. 209At imaging using the x-ray peak revealed a biodistribution that matched the known distribution of free astatide, and in vivo image-based measurements of 209At uptake in organs of interest matched ex vivo measurements within 10%. We have acquired the first 209At SPECT images and demonstrated the ability of quantitative SPECT imaging with 209At to accurately determine astatine biodistributions with high spatial and temporal resolution.

MR-compatibility of a high-resolution small animal PET insert operating inside a 7 T MRI.

A full-ring PET insert consisting of 16 PET detector modules was designed and constructed to fit within the 114 mm diameter gradient bore of a Bruker 7 T MRI. The individual detector modules contain two silicon photomultiplier (SiPM) arrays, dual-layer offset LYSO crystal arrays, and high-definition multimedia interface (HDMI) cables for both signal and power transmission. Several different RF shielding configurations were assessed prior to construction of a fully assembled PET insert using a combination of carbon fibre and copper foil for RF shielding. MR-compatibility measurements included field mapping of the static magnetic field (B 0) and the time-varying excitation field (B 1) as well as acquisitions with multiple pulse sequences: spin echo (SE), rapid imaging with refocused echoes (RARE), fast low angle shot (FLASH) gradient echo, and echo planar imaging (EPI). B 0 field maps revealed a small degradation in the mean homogeneity (+0.1 ppm) when the PET insert was installed and operating. No significant change was observed in the B 1 field maps or the image homogeneity of various MR images, with a 9% decrease in the signal-to-noise ratio (SNR) observed only in EPI images acquired with the PET insert installed and operating. PET detector flood histograms, photopeak amplitudes, and energy resolutions were unchanged in individual PET detector modules when acquired during MRI operation. There was a small baseline shift on the PET detector signals due to the switching amplifiers used to power MRI gradient pulses. This baseline shift was observable when measured with an oscilloscope and varied as a function of the gradient duty cycle, but had no noticeable effect on the performance of the PET detector modules. Compact front-end electronics and effective RF shielding led to minimal cross-interference between the PET and MRI systems. Both PET detector and MRI performance was excellent, whether operating as a standalone system or a hybrid PET/MRI.

Combined PET/MRI: from Status Quo to Status Go. Summary Report of the Fifth International Workshop on PET/MR Imaging; February 15-19, 2016; Tübingen, Germany.

This article provides a collaborative perspective of the discussions and conclusions from the fifth international workshop of combined positron emission tomorgraphy (PET)/magnetic resonance imaging (MRI) that was held in Tübingen, Germany, from February 15 to 19, 2016. Specifically, we summarise the second part of the workshop made up of invited presentations from active researchers in the field of PET/MRI and associated fields augmented by round table discussions and dialogue boards with specific topics. This year, this included practical advice as to possible approaches to moving PET/MRI into clinical routine, the use of PET/MRI in brain receptor imaging, in assessing cardiovascular diseases, cancer, infection, and inflammatory diseases. To address perceived challenges still remaining to innovatively integrate PET and MRI system technologies, a dedicated round table session brought together key representatives from industry and academia who were engaged with either the conceptualisation or early adoption of hybrid PET/MRI systems. Discussions during the workshop highlighted that emerging unique applications of PET/MRI such as the ability to provide multi-parametric quantitative and visual information which will enable not only overall disease detection but also disease characterisation would eventually be regarded as compelling arguments for the adoption of PET/MR. However, as indicated by previous workshops, evidence in favour of this observation is only growing slowly, mainly due to the ongoing inability to pool data cohorts from independent trials as well as different systems and sites. The participants emphasised that moving from status quo to status go entails the need to adopt standardised imaging procedures and the readiness to act together prospectively across multiple PET/MRI sites and vendors.

Premature termination mutations in exon 3 of the SMN1 gene are associated with exon skipping and a relatively mild SMA phenotype.

Autosomal recessive spinal muscular atrophy (SMA) is a common motor neuron disease caused by absence or mutation in the survival motor neuron (SMN1) gene. SNM1 and a nearly identical copy, SMN2, encode identical proteins, but SMN2 only produces a little full length protein due to alternative splicing. The level of functional SMN protein and the number of SMN2 genes correlate with the clinical phenotype ranging from severe to very mild. Here, we report on premature termination mutations in SMN1 exon 3 (425del5 and W102X) which induce skipping of the mutated exon. The novel nonsense mutation W102X was detected in two patients with a relatively mild phenotype who had only two copies of the SMN2 gene, a number that has previously been found associated with the severe form of SMA. We show that the shortened transcripts are translated into predicted in frame protein isoforms. Aminoglycoside treatment suppressed the nonsense mutation in cultured cells and abolished exon skipping. Fibroblasts from both patients show a high number of nuclear structures containing SMN protein (gems). These findings suggest that the protein isoform lacking the exon 3 encoded region contributes to the formation of the nuclear protein complex which may account for the milder clinical phenotype.

A reversible tracer analysis approach to the study of effective dopamine turnover.

Changes in dopamine turnover resulting from disease states such as Parkinson's disease may be reflected in corresponding changes in the kinetics of the positron emission tomographic tracer [(18)F]fluorodopa. The authors had previously refined the conventional irreversible-tracer graphical approach to determine both the uptake rate constant K(i) and the rate constant kloss that describes the slow loss of the trapped kinetic component. Because these parameters change in the opposite sense with disease, their ratios may be more powerfully discriminating than either one alone. The ratio k(loss)/K(i) is indicative of effective dopamine turnover. Its inverse, K(i)/k(loss), can be interpreted as the effective distribution volume (EDV) of the specific uptake compartment referred to the fluorodopa concentration in plasma. Here the authors present a new approach to the estimation of EDV based on reversible-tracer graphical methods. When implemented with a plasma input function, the method evaluates EDV directly. When implemented with a tissue input function, the outcome is proportional to the ratio of the distribution volumes of the specific uptake and precursor compartments. Comparison of the new and previous approaches strongly validates this alternative approach to the study of effective dopamine turnover.

Analysis of four dopaminergic tracers kinetics using two different tissue input function methods.

The integrity of the dopaminergic system can be studied using positron emission tomography. The presynaptic tracers [11C]-methylphenidate and [11C]dihydrotetrabenazine (DTBZ) are used to investigate the dopamine transporter availability, the dopamine vesicular transporter integrity; the postsynaptic tracers [11C]-raclopride and [11C]-Schering 23990 (SCH) are used to probe the D2 and D1 receptors. These are reversible tracers, where the binding potential (BP) = Bmax/Kd often is used to quantify the amount of their specific binding to the sites of interest. The simplified tissue input methods to calculate BP are attractive, since they do not require a blood input function. The suitability and performance of two such methods were evaluated: the Logan graphical tissue method, and the Lammertsma reference tissue method (RTM). The BP estimates obtained with the two methods were nearly identical in most cases, with similar reliability and reproducibility indicating that all four tracers satisfy the assumptions required by each method. The correlations among the fitted parameters obtained from the RTM were estimated and were found not to introduce noticeable bias in the RTM BP and R1 estimates. R1 showed low intersubject and intrasubject variability. The k2 estimate showed good reliability for SCH with cerebellar input function and DTBZ with occipital input function.

Effect of dopamine loss and the metabolite 3-O-methyl-[18F]fluoro-dopa on the relation between the 18F-fluorodopa tissue input uptake rate constant Kocc and the [18F]fluorodopa plasma input uptake rate constant Ki.

Parkinson disease is characterized by the loss of dopaminergic neurons, thus decreasing the system's ability to produce and store dopamine (DA). Such ability is often investigated using 18F-fluorodopa (FD) positron emission tomography. A commonly used model to investigate the DA synthesis and storage rate is the modified Patlak graphical approach. This approach allows for both plasma and tissue input functions, yielding the respective uptake rate constants K(i) and K(occ). This method requires the presence of an irreversible compartment and the absence of any nontrapped tracer metabolite. In the case of K(occ), this last assumption is violated by the presence of the FD metabolite 3-O-methyl-[18F]fluoro-dopa (3OMFD), which makes the K(occ) evaluation susceptible to a downward bias. It was found that both K(i) and K(occ) are influenced by DA loss and thus are not pure measures of DA synthesis and storage. In the case of K(occ), the presence of 3OMFD exacerbates the effect of DA egress, thus introducing a disease-dependent bias in the K(occ) determination. These findings imply that K(i) and K(occ) provide different assessments of disease severity and that, as disease progresses, K(i) and especially K(occ) become more related to DA storage capacity and less to the DA synthesis rate.

Apomorphine-induced changes in synaptic dopamine levels: positron emission tomography evidence for presynaptic inhibition.

The authors developed a novel positron emission tomography method to estimate changes in the synaptic level of dopamine ([DA]) induced by direct dopamine agonists (for example, apomorphine) in patients with Parkinson disease. The method is based on the typical asymmetry of the nigrostriatal lesion that often occurs in Parkinson disease. Using the between-side difference (ipsilateral (I) and contralateral (C) putamen to the more affected body side) of the inverse of the putamen [11C]raclopride binding potential (BP), the authors obtained [equation: see text] at baseline (that is, before apomorphine administration) and [equation: see text] after apomorphine administration (assuming the concentration of apomorphine is equal in both putamina). The between-side difference in the estimated synaptic concentration of dopamine (diff[DA]) should remain constant unless apomorphine affects dopamine release differently between the two sides. The authors found that apomorphine given subcutaneously at doses of 0.03 and 0.06 mg/kg induced significant changes in their estimate of diff[DA] (P < 0.05). Such changes were more pronounced when only patients with a stable response to levodopa were considered (P < 0.01). These findings provide in vivo evidence that direct dopamine agonists can inhibit the release of endogenous dopamine. The authors propose that this effect is mainly mediated by the activation of presynaptic D2/D3 dopamine receptors.

PET studies of parkinsonism associated with mutation in the alpha-synuclein gene.

OBJECTIVE: To assess the pattern of dopaminergic abnormalities in a Greek-American kindred (family H) with autosomal dominantly inherited, levodopa-responsive parkinsonism caused by a mutation of the gene encoding alpha-synuclein. BACKGROUND: Mutations of alpha-synuclein have been associated recently with dominantly inherited, levodopa-responsive parkinsonism. The pattern of dopamine deficiency and status of postsynaptic dopamine receptors in this condition have not been reported previously. The authors followed a large, six-generation family in whom the affected members carry the recently reported G209A mutation in the gene encoding alpha-synuclein. METHODS: The authors studied four affected and two clinically unaffected gene-negative members of family H using [18F]-6-fluoro-L-dopa (FD) and [11C]-raclopride (RAC) PET to assess presynaptic dopaminergic function and dopamine D2 receptors. The results were compared with normal subjects and patients with sporadic, idiopathic PD (IP). RESULTS: In affected individuals, FD uptake was reduced in both the caudate and the putamen, but the putamen was affected more severely than the caudate, as seen in IP. RAC binding was within the normal range, but the ratio of RAC binding in the putamen to that in the caudate was increased in affected members of family H. This pattern is similar to that seen in IP. CONCLUSIONS: PET of the nigrostriatal system in parkinsonism associated with a mutation in the ac-synuclein gene indicates that it results in a pattern of dopamine deficiency, with preserved D2 binding, indistinguishable from IP.

Bilateral human fetal striatal transplantation in Huntington's disease.

BACKGROUND: Transplanted striatal cells have been demonstrated to survive, grow, establish afferent and efferent connections, and improve behavioral signs in animal models of Huntington's disease (HD). OBJECTIVE: To evaluate feasibility and safety and to provide preliminary information regarding the efficacy of bilateral human fetal striatal transplantation in HD. METHODS: Seven symptomatic patients with genetically confirmed HD underwent bilateral stereotactic transplantation of two to eight fetal striata per side in two staged procedures. Tissue was dissected from the lateral half of the lateral ventricular eminence of donors 8 to 9 weeks postconception. Subjects received cyclosporine for 6 months. RESULTS: Three subjects developed subdural hemorrhages (SDHs) and two required surgical drainage. One subject died 18 months after surgery from probable cardiac arrhythmia secondary to severe atherosclerotic cardiac disease. Autopsy demonstrated clearly demarcated grafts of typical developing striatal morphology, with host-derived dopaminergic fibers extending into the grafts and no evidence of immune rejection. Other adverse events were generally mild and transient. Mean Unified HD Rating Scale (UHDRS) motor scores were 32.9 plus minus 6.2 at baseline and 29.7 plus minus 7.5 12 months after surgery (p = 0.24). Post-hoc analysis, excluding one subject who experienced cognitive and motor deterioration after the development of symptomatic bilateral SDHs, found that UHDRS motor scores were 33.8 plus minus 6.2 at baseline and 27.5 plus minus 5.2 at 12 months (p = 0.03). CONCLUSIONS: Transplantation of human fetal striatal cells is feasible and survival of transplanted cells was demonstrated. Patients with moderately advanced HD are at risk for SDH after transplantation surgery.

Quantitative comparison of three- and two-dimensional PET with human brain studies.

UNLABELLED: The aim of this study was to test the quantitation accuracy of three-dimensional PET in brain scanning. METHODS: Three-dimensional data from 11 human subjects were tested using 11C-dihydrotetrabenazine, 11C-Schering 23390 and 18F-FDG as tracers. Two-dimensional scans were performed on the same subjects and the distribution volume, distribution volume ratio and local metabolic rate of glucose (LMRGlu) values obtained from these were used as reference. Three-dimensional data were processed as follows: iterative convolution subtraction scatter correction, detector normalization including radial and axial geometric factors, attenuation correction extracted from a two-dimensional transmission scan, Kinahan-Rogers reconstruction and region-of-interest-based sensitivity calibration. RESULTS: No major systematic differences between the two methods were found. The agreement between the two-dimensional and three-dimensional data was within 5%. Although statistical analysis generally did not show this difference to be significant, reliability analysis indicated that comparing two-dimensional and three-dimensional data might introduce some inaccuracies. CONCLUSION: Three-dimensional PET yields quantitatively valid results for brain scanning.

Recovery of the human striatal signal in a slice oriented positron emission tomograph.

The human striatum is small enough for partial volume effects to be important when imaged in positron tomographs with slice widths 10 mm or greater. The combination of interslice distance and slice width in such tomographs results in an axial undersampling of the striatal activity which introduces the additional problem of variation of axial recovery as a function of position of the striatum along the tomograph axis. Using striatal phantoms, we have developed a method that corrects the recovered striatal signal to a maximum value equivalent to that measured when the object is centered with respect to a slice. This makes the recovery independent of the axial position of the striatum. The method also provides an estimate of the total striatal activity by integrating the axial image intensity distribution along the tomograph axis. The method is able to detect and correct for relative axial tilt of the left and right striatum. We applied it to 26 human [18F]-6-L-fluorodopa scans and obtained an average uptake rate constant k value of 0.25 +/- 0.05 ml/min/striatum and a left to right k value percentage asymmetry of 0.1% +/- 6.3%.

Graphical analysis of 6-fluoro-L-dopa trapping: effect of inhibition of catechol-O-methyltransferase.

UNLABELLED: Graphical methods to analyze tracer time-course data allow reliable quantitation of the rate of incorporation of tracer from plasma into a "trapped" kinetic component, even when the details of the kinetic model are unknown. Applications of the method over long time periods often expose the slow reversibility of the trapping process. In the extended graphical method, both trapping rate and a presumed first-order loss rate constant are estimated simultaneously from the time-course data. METHODS: We applied the extended graphical method to 6-fluoro-L-dopa (6-FD), simultaneously estimating the rate of uptake (Ki) and the rate constant for loss from the trapped component (K(loss)) in a single fitting procedure. We applied this approach to study the effects of two catechol-O-methyl-transferase inhibitors on the kinetics of 6-FD in cynomolgus monkeys. RESULTS: Inhibition of peripheral O-methylation with either inhibitor, confirmed by high-performance liquid chromatography analysis of labeled compounds in arterial plasma, had no significant effect on Ki, in agreement with previously reported studies. In contrast, tolcapone, a catechol-O-methyl-transferase inhibitor, having central effects in addition to peripheral effects at the dosage used, decreased K(loss) by 40% from control values (p < 0.002), whereas nitecapone, which has no known central activity, had no significant effect. CONCLUSION: This method provides insight into the neurochemical basis for the kinetic behavior of 6-FD in both health and disease and may be used to define the action of centrally active drugs that influence the metabolism of dopamine.

Normalization for 3D PET with a low-scatter planar source and measured geometric factors.

For 3D PET normalization methods, a balance must be struck between statistical accuracy and individual detector or line-of-response (LOR) fidelity. Methods with potentially the best LOR accuracy tend to be statistically poor, while techniques to improve the statistical quality tend to reduce the individual detector fidelity. We have developed and implemented a 3D PET normalization method for our ECAT 953B scanner (Siemens/CTI) that determines the detector normalization factors (NFs) as a product of a four-dimensional matrix of measured geometric factors (GFs) and single detector efficiency factors (epsilon). The effects of various alterations to the algorithm on the accuracy of the normalization have been examined through the impact on reconstructed images. An accurate set of GFs is crucial, as inaccurate NFs can result if LORs with similar but not identical geometric symmetries are grouped together. The general method can be extended to other tomographs, although the dimensionality of a GF matrix may be scanner-specific; the key is to determine the optimal number of dimensions in the GF matrix. The GFs for our scanner are specified by: (i) the two detector rings for each LOR; (ii) the radial distance of the LOR from the tomograph centre; and (iii) the positions within the detector block of the two crystals defining the LOR. Some residual radial non-uniformities are present in all the NF variations we examined. For the NF method presented here, the radial non-uniformities are attributed to the interaction between object-dependent scatter and normalization. Results indicate that this non-uniformity is detectable for scans with as few as 13 million total counts.

A phantom study evaluating the quantitative aspect of 3D PET imaging of the brain.

Phantom studies are used to develop a reliable quantitative data processing protocol for 3D PET brain scanning for conditions typically encountered in FDG and neuroreceptor brain imaging. These protocols often span several half-lives of the injected radiotracer thus resulting in a greatly varying statistical content of the acquired data over the study duration. Detector normalization, scatter correction and their interplay over a wide range of statistical content of acquired data were evaluated. Overall sensitivity calibration factors were determined after all other quantification corrections were applied to the data. The result is an optimum data processing protocol that includes an iterative convolution subtraction scatter correction method, a normalization procedure that takes into account the geometric properties of the scanner and a region of interest based calibration procedure, applied in this order. This protocol yields a 3D PET quantification accuracy within approximately 3% of independently measured concentration values for scanning conditions that include variation in the number of acquired counts from one million to several hundred millions and variation in size and shape from a 20 cm diameter phantom to a tapered phantom with minimum cross section of 3.7 x 14.5 cm2. This performance is comparable with that of the 2D acquisition mode.

The performance of the single-slice rebinning technique for imaging the human striatum as evaluated by phantom studies.

The single-slice rebinning (SSRB) technique has been implemented on a ring-based tomograph with removable septa and its performance evaluated with phantom studies. The studies involved the determination of the plane efficiency profile and the mapping of the point spread function. The image quality obtained using the SSRB technique has been compared to that obtained using the standard 3D reconstruction method in terms of the following figures of merit: recovery as a function of source axial position and slice averaging, contrast recovery and contrast to background noise ratio as a function of source axial position. The effect of the region of interest (ROI) size on such a comparison has been evaluated. Finally the effects of different contrast conditions and ROI position on a striatal phantom image have been thoroughly explored.

Statistical list-mode image reconstruction for the high resolution research tomograph.

We have investigated statistical list-mode reconstruction applicable to a depth-encoding high resolution research tomograph. An image non-negativity constraint has been employed in the reconstructions and is shown to effectively remove the overestimation bias introduced by the sinogram non-negativity constraint. We have furthermore implemented a convergent subsetized (CS) list-mode reconstruction algorithm, based on previous work (Hsiao et al 2002 Conf. Rec. SPIE Med. Imaging 4684 10-19; Hsiao et al 2002 Conf. Rec. IEEE Int. Symp. Biomed. Imaging 409-12) on convergent histogram OSEM reconstruction. We have demonstrated that the first step of the convergent algorithm is exactly equivalent (unlike the histogram-mode case) to the regular subsetized list-mode EM algorithm, while the second and final step takes the form of additive updates in image space. We have shown that in terms of contrast, noise as well as FWHM width behaviour, the CS algorithm is robust and does not result in limit cycles. A hybrid algorithm based on the ordinary and the convergent algorithms is also proposed, and is shown to combine the advantages of the two algorithms (i.e. it is able to reach a higher image quality in fewer iterations while maintaining the convergent behaviour), making the hybrid approach a good alternative to the ordinary subsetized list-mode EM algorithm.