Virtual cylindrical PET for efficient DOI image reconstruction with sub-millimetre resolution.

OBJECTIVE: Image reconstruction in high resolution, narrow bore PET scanners with depth of interaction (DOI) capability presents a substantial computational challenge due to the very high sampling in detector and image space. The aim of this study is to evaluate the use of a virtual cylinder in reducing the number of lines of response (LOR) for DOI-based reconstruction in high resolution PET systems while maintaining uniform sub-millimetre spatial resolution. Approach: Virtual geometry was investigated using the awake animal mousePET as a high resolution test case. Using GATE, we simulated the physical scanner and three virtual cylinder implementations with detector size 0.74 mm, 0.47 mm and 0.36 mm (vPET1, vPET2 and vPET3, respectively). The virtual cylinder condenses physical LORs stemming from various crystal pairs and DOI combinations, and which intersect a single virtual detector pair, into a single virtual LOR. Quantitative comparisons of the point spread function (PSF) at various positions within the field of view (FOV) were compared for reconstructions based on the vPET implementations and the physical scanner. We also assessed the impact of the anisotropic PSFs by reconstructing images of a micro Derenzo phantom. Main results: All virtual cylinder implementations achieved LOR data compression of at least 50% for DOI PET reconstruction. PSF anisotropy in radial and tangential profiles was chiefly influenced by DOI resolution and only marginally by virtual detector size. Spatial degradation introduced by virtual cylinders was most prominent in the axial profile. All virtual cylinders achieved sub-millimetre volumetric resolution across the FOV when 6-bin DOI reconstructions (3.3 mm DOI resolution) were performed. Using vPET2 with 6 DOI bins yielded nearly identical reconstructions to the non-virtual case in the transaxial plane, with an LOR compression ratio of 86%. Resolution modelling significantly reduced the effects of the asymmetric PSF arising from the non-cylindrical geometry of mousePET. Significance: Narrow bore and high resolution PET scanners require detectors with DOI capability, leading to computationally demanding reconstructions due to the large number of LORs. In this study, we show that DOI PET reconstruction with 50-86% LOR compression is possible using virtual cylinders while maintaining sub-millimetre spatial resolution throughout the FOV. The methodology and analysis can be extended to other scanners with DOI capability intended for high resolution PET imaging. .

Optoacoustic biomarkers of lipids, hemorrhage and inflammation in carotid atherosclerosis.

Imaging plays a critical role in exploring the pathophysiology and enabling the diagnostics and therapy assessment in carotid artery disease. Ultrasonography, computed tomography, magnetic resonance imaging and nuclear medicine techniques have been used to extract of known characteristics of plaque vulnerability, such as inflammation, intraplaque hemorrhage and high lipid content. Despite the plethora of available techniques, there is still a need for new modalities to better characterize the plaque and provide novel biomarkers that might help to detect the vulnerable plaque early enough and before a stroke occurs. Optoacoustics, by providing a multiscale characterization of the morphology and pathophysiology of the plaque could offer such an option. By visualizing endogenous (e.g., hemoglobin, lipids) and exogenous (e.g., injected dyes) chromophores, optoacoustic technologies have shown great capability in imaging lipids, hemoglobin and inflammation in different applications and settings. Herein, we provide an overview of the main optoacoustic systems and scales of detail that enable imaging of carotid plaques in vitro, in small animals and humans. Finally, we discuss the limitations of this novel set of techniques while investigating their potential to enable a deeper understanding of carotid plaque pathophysiology and possibly improve the diagnostics in future patients with carotid artery disease.

MRI tractography-guided PET image reconstruction regularisation using connectome-based nonlocal means filtering.

Positron emission tomography (PET) molecular biomarkers and diffusion magnetic resonance imaging (dMRI) derived information show associations and highly complementary information in a number of neurodegenerative conditions, including Alzheimer's disease. Diffusion MRI provides valuable information about the microstructure and structural connectivity (SC) of the brain which could guide and improve the PET image reconstruction when such associations exist. However, this potental has not been previously explored. In the present study, we propose a CONNectome-based non-local means one-atep late maximuma posteriori(CONN-NLM-OSLMAP) method, which allows diffusion MRI-derived connectivity information to be incorporated into the PET iterative image reconstruction process, thus regularising the estimated PET images. The proposed method was evaluated using a realistic tau-PET/MRI simulated phantom, demonstrating more effective noise reduction and lesion contrast improvement, as well as the lowest overall bias compared with both a median filter applied as an alternative regulariser and CONNectome-based non-local means as a post-reconstruction filter. By adding complementary SC information from diffusion MRI, the proposed regularisation method offers more useful and targeted denoising and regularisation, demonstrating the feasibility and effectiveness of integrating connectivity information into PET image reconstruction.

Growth of chemical gardens in gaseous acidic atmospheres.

The generation of chemobrionic architectures through slow injection of aqueous silicate solution in gaseous TiCl4 is demonstrated. The tubes were characterized by XRD, SEM and wet chemistry control experiments, and their mechanism of formation was unraveled. These structures serve as laboratory models for calthemites or soda straws.

A chemogenomic approach is required for effective treatment of amyotrophic lateral sclerosis.

ALS is a fatal untreatable disease involving degeneration of motor neurons. Μultiple causative genes encoding proteins with versatile functions have been identified indicating that diverse biological pathways lead to ALS. Chemical entities still represent a promising choice to delay ALS progression, attenuate symptoms and/or increase life expectancy, but also gene-based and stem cell-based therapies are in the process of development, and some are tested in clinical trials. Various compounds proved effective in transgenic models overexpressing distinct ALS causative genes unfortunately though, they showed no efficacy in clinical trials. Notably, while animal models provide a uniform genetic background for preclinical testing, ALS patients are not stratified, and the distinct genetic forms of ALS are treated as one group, which could explain the observed discrepancies between treating genetically homogeneous mice and quite heterogeneous patient cohorts. We suggest that chemical entity-genotype correlation should be exploited to guide patient stratification for pharmacotherapy, that is administered drugs should be selected based on the ALS genetic background.

Plausible Emergence of Biochemistry in Enceladus Based on Chemobrionics.

Saturn's satellite Enceladus is proposed to have a soda-type subsurface ocean with temperature able to support life and an iron ore-based core. Here, it was demonstrated that ocean chemistry related to Enceladus can support the development of Fe-based hydrothermal vents, one of the places suggested to be the cradle of life. The Fe-based chemical gardens were characterized with Fourier-transform (FT)IR spectroscopy and XRD. The developed chemobrionic structures catalyzed the condensation polymerization of simple organic prebiotic molecules to kerogens. Further, they could passively catalyze the condensation of the prebiotic molecule formamide to larger polymers, suggesting that elementary biochemical precursors could have emerged in Enceladus.

MLC tracking for lung SABR is feasible, efficient and delivers high-precision target dose and lower normal tissue dose.

BACKGROUND AND PURPOSE: The purpose of this work is to present the clinical experience from the first-in-human trial of real-time tumor targeting via MLC tracking for stereotactic ablative body radiotherapy (SABR) of lung lesions. METHODS AND MATERIALS: Seventeen patients with stage 1 non-small cell lung cancer (NSCLC) or lung metastases were included in a study of electromagnetic transponder-guided MLC tracking for SABR (NCT02514512). Patients had electromagnetic transponders inserted near the tumor. An MLC tracking SABR plan was generated with planning target volume (PTV) expanded 5 mm from the end-exhale gross tumor volume (GTV). A clinically approved comparator plan was generated with PTV expanded 5 mm from a 4DCT-derived internal target volume (ITV). Treatment was delivered using a standard linear accelerator to continuously adapt the MLC based on transponder motion. Treated volumes and reconstructed delivered dose were compared between MLC tracking and comparator ITV-based treatment. RESULTS: All seventeen patients were successfully treated with MLC tracking (70 successful fractions). MLC tracking treatment delivery time averaged 8 minutes. The time from the start of CBCT to the end of treatment averaged 22 minutes. The MLC tracking PTV for 16/17 patients was smaller than the ITV-based PTV (range -1.6% to 44% reduction, or -0.6 to 18 cc). Reductions in mean lung dose (27 cGy) and V20Gy (50 cc) were statistically significant (p < 0.02). Reconstruction of treatment doses confirmed a statistically significant improvement in delivered GTV D98% (p < 0.05) from planned dose compared with the ITV-based plans. CONCLUSION: The first treatments with lung MLC tracking have been successfully performed in seventeen SABR patients. MLC tracking for lung SABR is feasible, efficient and delivers high-precision target dose and lower normal tissue dose.

A Closed Chemobrionic System as a Biochemical Delivery Platform.

Inorganic cells bearing calcium silicate membranes were prepared and resembled closed chemical gardens. It was demonstrated that these inorganic cells can successfully be loaded with natural products, proteins and plasmid DNA, and their cargo can be released in a controlled manner. These cells demonstrated the ability of chemical gardens to act as platforms for the sustained delivery of biomolecules and are expected to introduce chemical gardens in the field of biosciences.

Open-field PET: Simultaneous brain functional imaging and behavioural response measurements in freely moving small animals.

A comprehensive understanding of how the brain responds to a changing environment requires techniques capable of recording functional outputs at the whole-brain level in response to external stimuli. Positron emission tomography (PET) is an exquisitely sensitive technique for imaging brain function but the need for anaesthesia to avoid motion artefacts precludes concurrent behavioural response studies. Here, we report a technique that combines motion-compensated PET with a robotically-controlled animal enclosure to enable simultaneous brain imaging and behavioural recordings in unrestrained small animals. The technique was used to measure in vivo displacement of [11C]raclopride from dopamine D2 receptors (D2R) concurrently with changes in the behaviour of awake, freely moving rats following administration of unlabelled raclopride or amphetamine. The timing and magnitude of [11C]raclopride displacement from D2R were reliably estimated and, in the case of amphetamine, these changes coincided with a marked increase in stereotyped behaviours and hyper-locomotion. The technique, therefore, allows simultaneous measurement of changes in brain function and behavioural responses to external stimuli in conscious unrestrained animals, giving rise to important applications in behavioural neuroscience.

Direct Estimation of Voxel-Wise Neurotransmitter Response Maps From Dynamic PET Data.

Computational methods, such as the linear parametric neurotransmitter PET (lp-ntPET) method, have been developed to characterize the transient changes in radiotracer kinetics in the target tissue during endogenous neurotransmitter release. In this paper, we describe and evaluate a parametric reconstruction algorithm that uses an expectation maximization framework, along with the lp-ntPET model, to estimate the endogenous neurotransmitter response to stimuli directly from the measured PET data. Computer simulations showed that the proposed direct reconstruction method offers improved accuracy and precision for the estimated timing parameters of the neurotransmitter response at the voxel level ( td=1±2 min, for activation onset bias and standard deviation) compared with conventional post reconstruction modeling ( td=4±7 min). In addition, we applied the proposed direct parameter estimation methodology to a [11C]raclopride displacement study of an awake rat and generated parametric maps illustrating the magnitude of ligand displacement from striatum. Although the estimated parametric maps of activation magnitude obtained from both direct and post reconstruction methodologies suffered from false positive activations, the proposed direct reconstruction framework offered more reliable parametric maps when the activation onset parameter was constrained.

Determining Glucose Metabolism Kinetics Using 18F-FDG Micro-PET/CT.

This paper describes the use of 18F-FDG and micro-PET/CT imaging to determine in vivo glucose metabolism kinetics in mice (and is transferable to rats). Impaired uptake and metabolism of glucose in multiple organ systems due to insulin resistance is a hallmark of type 2 diabetes. The ability of this technique to extract an image-derived input function from the vena cava using an iterative deconvolution method eliminates the requirement of the collection of arterial blood samples. Fitting of tissue and vena cava time activity curves to a two-tissue, three compartment model permits the estimation of kinetic micro-parameters related to the 18F-FDG uptake from the plasma to the intracellular space, the rate of transport from intracellular space to plasma and the rate of 18F-FDG phosphorylation. This methodology allows for multiple measures of glucose uptake and metabolism kinetics in the context of longitudinal studies and also provides insights into the efficacy of therapeutic interventions.

Motion compensation using origin ensembles in awake small animal positron emission tomography.

In emission tomographic imaging, the stochastic origin ensembles algorithm provides unique information regarding the detected counts given the measured data. Precision in both voxel and region-wise parameters may be determined for a single data set based on the posterior distribution of the count density allowing uncertainty estimates to be allocated to quantitative measures. Uncertainty estimates are of particular importance in awake animal neurological and behavioral studies for which head motion, unique for each acquired data set, perturbs the measured data. Motion compensation can be conducted when rigid head pose is measured during the scan. However, errors in pose measurements used for compensation can degrade the data and hence quantitative outcomes. In this investigation motion compensation and detector resolution models were incorporated into the basic origin ensembles algorithm and an efficient approach to computation was developed. The approach was validated against maximum liklihood-expectation maximisation and tested using simulated data. The resultant algorithm was then used to analyse quantitative uncertainty in regional activity estimates arising from changes in pose measurement precision. Finally, the posterior covariance acquired from a single data set was used to describe correlations between regions of interest providing information about pose measurement precision that may be useful in system analysis and design. The investigation demonstrates the use of origin ensembles as a powerful framework for evaluating statistical uncertainty of voxel and regional estimates. While in this investigation rigid motion was considered in the context of awake animal PET, the extension to arbitrary motion may provide clinical utility where respiratory or cardiac motion perturb the measured data.

List-mode image reconstruction for positron emission tomography using tetrahedral voxels.

Image space decomposition based on tetrahedral voxels are interesting candidates for use in emission tomography. Tetrahedral voxels provide many of the advantages of point clouds with irregular spacing, such as being intrinsically multi-resolution, yet they also serve as a volumetric partition of the image space and so are comparable to more standard cubic voxels. Additionally, non-rigid displacement fields can be applied to the tetrahedral mesh in a straight-forward manner. So far studies incorporating tetrahedral decomposition of the image space have concentrated on pre-calculated, node-based, system matrix elements which reduces the flexibility of the tetrahedral approach and the capacity to accurately define regions of interest. Here, a list-mode on-the-fly calculation of the system matrix elements is described using a tetrahedral decomposition of the image space and volumetric elements-voxels. The algorithm is demonstrated in the context of awake animal PET which may require both rigid and non-rigid motion compensation, as well as quantification within small regions of the brain. This approach allows accurate, event based, motion compensation including non-rigid deformations.

4D PET iterative deconvolution with spatiotemporal regularization for quantitative dynamic PET imaging.

Quantitative measurements in dynamic PET imaging are usually limited by the poor counting statistics particularly in short dynamic frames and by the low spatial resolution of the detection system, resulting in partial volume effects (PVEs). In this work, we present a fast and easy to implement method for the restoration of dynamic PET images that have suffered from both PVE and noise degradation. It is based on a weighted least squares iterative deconvolution approach of the dynamic PET image with spatial and temporal regularization. Using simulated dynamic [(11)C] Raclopride PET data with controlled biological variations in the striata between scans, we showed that the restoration method provides images which exhibit less noise and better contrast between emitting structures than the original images. In addition, the method is able to recover the true time activity curve in the striata region with an error below 3% while it was underestimated by more than 20% without correction. As a result, the method improves the accuracy and reduces the variability of the kinetic parameter estimates calculated from the corrected images. More importantly it increases the accuracy (from less than 66% to more than 95%) of measured biological variations as well as their statistical detectivity.

Full field spatially-variant image-based resolution modelling reconstruction for the HRRT.

Accurate characterisation of the scanner's point spread function across the entire field of view (FOV) is crucial in order to account for spatially dependent factors that degrade the resolution of the reconstructed images. The HRRT users' community resolution modelling reconstruction software includes a shift-invariant resolution kernel, which leads to transaxially non-uniform resolution in the reconstructed images. Unlike previous work to date in this field, this work is the first to model the spatially variant resolution across the entire FOV of the HRRT, which is the highest resolution human brain PET scanner in the world. In this paper we developed a spatially variant image-based resolution modelling reconstruction dedicated to the HRRT, using an experimentally measured shift-variant resolution kernel. Previously, the system response was measured and characterised in detail across the entire FOV of the HRRT, using a printed point source array. The newly developed resolution modelling reconstruction was applied on measured phantom, as well as clinical data and was compared against the HRRT users' community resolution modelling reconstruction, which is currently in use. Results demonstrated improvements both in contrast and resolution recovery, particularly for regions close to the edges of the FOV, with almost uniform resolution recovery across the entire transverse FOV. In addition, because the newly measured resolution kernel is slightly broader with wider tails, compared to the deliberately conservative kernel employed in the HRRT users' community software, the reconstructed images appear to have not only improved contrast recovery (up to 20% for small regions), but also better noise characteristics.

Impact of extraneous mispositioned events on motion-corrected brain SPECT images of freely moving animals.

PURPOSE: Single photon emission computed tomography (SPECT) brain imaging of freely moving small animals would allow a wide range of important neurological processes and behaviors to be studied, which are normally inhibited by anesthetic drugs or precluded due to the animal being restrained. While rigid body motion of the head can be tracked and accounted for in the reconstruction, activity in the torso may confound brain measurements, especially since motion of the torso is more complex (i.e., nonrigid) and not well correlated with that of the head. The authors investigated the impact of mispositioned events and attenuation due to the torso on the accuracy of motion corrected brain images of freely moving mice. METHODS: Monte Carlo simulations of a realistic voxelized mouse phantom and a dual compartment phantom were performed. Each phantom comprised a target and an extraneous compartment which were able to move independently of each other. Motion correction was performed based on the known motion of the target compartment only. Two SPECT camera geometries were investigated: a rotating single head detector and a stationary full ring detector. The effects of motion, detector geometry, and energy of the emitted photons (hence, attenuation) on bias and noise in reconstructed brain regions were evaluated. RESULTS: The authors observed two main sources of bias: (a) motion-related inconsistencies in the projection data and (b) the mismatch between attenuation and emission. Both effects are caused by the assumption that the orientation of the torso is difficult to track and model, and therefore cannot be conveniently corrected for. The motion induced bias in some regions was up to 12% when no attenuation effects were considered, while it reached 40% when also combined with attenuation related inconsistencies. The detector geometry (i.e., rotating vs full ring) has a big impact on the accuracy of the reconstructed images, with the full ring detector being more advantageous. CONCLUSIONS: Motion-induced inconsistencies in the projection data and attenuation/emission mismatch are the two main causes of bias in reconstructed brain images when there is complex motion. It appears that these two factors have a synergistic effect on the qualitative and quantitative accuracy of the reconstructed images.

Evaluation of a direct 4D reconstruction method using generalised linear least squares for estimating nonlinear micro-parametric maps.

OBJECTIVE: Estimation of nonlinear micro-parameters is a computationally demanding and fairly challenging process, since it involves the use of rather slow iterative nonlinear fitting algorithms and it often results in very noisy voxel-wise parametric maps. Direct reconstruction algorithms can provide parametric maps with reduced variance, but usually the overall reconstruction is impractically time consuming with common nonlinear fitting algorithms. METHODS: In this work we employed a recently proposed direct parametric image reconstruction algorithm to estimate the parametric maps of all micro-parameters of a two-tissue compartment model, used to describe the kinetics of [[Formula: see text]F]FDG. The algorithm decouples the tomographic and the kinetic modelling problems, allowing the use of previously developed post-reconstruction methods, such as the generalised linear least squares (GLLS) algorithm. RESULTS: Results on both clinical and simulated data showed that the proposed direct reconstruction method provides considerable quantitative and qualitative improvements for all micro-parameters compared to the conventional post-reconstruction fitting method. Additionally, region-wise comparison of all parametric maps against the well-established filtered back projection followed by post-reconstruction non-linear fitting, as well as the direct Patlak method, showed substantial quantitative agreement in all regions. CONCLUSIONS: The proposed direct parametric reconstruction algorithm is a promising approach towards the estimation of all individual microparameters of any compartment model. In addition, due to the linearised nature of the GLLS algorithm, the fitting step can be very efficiently implemented and, therefore, it does not considerably affect the overall reconstruction time.

Markerless motion tracking of awake animals in positron emission tomography.

Noninvasive functional imaging of awake, unrestrained small animals using motion-compensation removes the need for anesthetics and enables an animal's behavioral response to stimuli or administered drugs to be studied concurrently with imaging. While the feasibility of motion-compensated radiotracer imaging of awake rodents using marker-based optical motion tracking has been shown, markerless motion tracking would avoid the risk of marker detachment, streamline the experimental workflow, and potentially provide more accurate pose estimates over a greater range of motion. We have developed a stereoscopic tracking system which relies on native features on the head to estimate motion. Features are detected and matched across multiple camera views to accumulate a database of head landmarks and pose is estimated based on 3D-2D registration of the landmarks to features in each image. Pose estimates of a taxidermal rat head phantom undergoing realistic rat head motion via robot control had a root mean square error of 0.15 and 1.8 mm using markerless and marker-based motion tracking, respectively. Markerless motion tracking also led to an appreciable reduction in motion artifacts in motion-compensated positron emission tomography imaging of a live, unanesthetized rat. The results suggest that further improvements in live subjects are likely if nonrigid features are discriminated robustly and excluded from the pose estimation process.

Isotope specific resolution recovery image reconstruction in high resolution PET imaging.

PURPOSE: Measuring and incorporating a scanner-specific point spread function (PSF) within image reconstruction has been shown to improve spatial resolution in PET. However, due to the short half-life of clinically used isotopes, other long-lived isotopes not used in clinical practice are used to perform the PSF measurements. As such, non-optimal PSF models that do not correspond to those needed for the data to be reconstructed are used within resolution modeling (RM) image reconstruction, usually underestimating the true PSF owing to the difference in positron range. In high resolution brain and preclinical imaging, this effect is of particular importance since the PSFs become more positron range limited and isotope-specific PSFs can help maximize the performance benefit from using resolution recovery image reconstruction algorithms. METHODS: In this work, the authors used a printing technique to simultaneously measure multiple point sources on the High Resolution Research Tomograph (HRRT), and the authors demonstrated the feasibility of deriving isotope-dependent system matrices from fluorine-18 and carbon-11 point sources. Furthermore, the authors evaluated the impact of incorporating them within RM image reconstruction, using carbon-11 phantom and clinical datasets on the HRRT. RESULTS: The results obtained using these two isotopes illustrate that even small differences in positron range can result in different PSF maps, leading to further improvements in contrast recovery when used in image reconstruction. The difference is more pronounced in the centre of the field-of-view where the full width at half maximum (FWHM) from the positron range has a larger contribution to the overall FWHM compared to the edge where the parallax error dominates the overall FWHM. CONCLUSIONS: Based on the proposed methodology, measured isotope-specific and spatially variant PSFs can be reliably derived and used for improved spatial resolution and variance performance in resolution recovery image reconstruction. The benefits are expected to be more substantial for more energetic positron emitting isotopes such as Oxygen-15 and Rubidium-82.

Simulation-based optimisation of the PET data processing for partial saturation approach protocols.

Positron emission tomography (PET) with [(11)C]Raclopride is an important tool for studying dopamine D2 receptor expression in vivo. [(11)C]Raclopride PET binding experiments conducted using the Partial Saturation Approach (PSA) allow the estimation of receptor density (B(avail)) and the in vivo affinity appK(D). The PSA is a simple, single injection, single scan experimental protocol that does not require blood sampling, making it ideal for use in longitudinal studies. In this work, we generated a complete Monte Carlo simulated PET study involving two groups of scans, in between which a biological phenomenon was inferred (a 30% decrease of B(avail)), and used it in order to design an optimal data processing chain for the parameter estimation from PSA data. The impact of spatial smoothing, noise removal and image resolution recovery technique on the statistical detection was investigated in depth. We found that image resolution recovery using iterative deconvolution of the image with the system point spread function associated with temporal data denoising greatly improves the accuracy and the statistical reliability of detecting the imposed phenomenon. Before optimisation, the inferred B(avail) variation between the two groups was underestimated by 42% and detected in 66% of cases, while a false decrease of appK(D) by 13% was detected in more than 11% of cases. After optimisation, the calculated B(avail) variation was underestimated by only 3.7% and detected in 89% of cases, while a false slight increase of appK(D) by 3.7% was detected in only 2% of cases. We found during this investigation that it was essential to adjust a factor that accounts for difference in magnitude between the non-displaceable ligand concentrations measured in the target and in the reference regions, for different data processing pathways as this ratio was affected by different image resolutions.