Improved regional activity quantitation in nuclear medicine using a new approach to correct for tissue partial volume and spillover effects

CT- and MRI-Aided Fluorescence Tomography Reconstructions for Biodistribution Analysis

OBJECTIVES

Optical fluorescence imaging can track the biodistribution of fluorophore-labeled drugs, nanoparticles, and antibodies longitudinally. In hybrid computed tomography-fluorescence tomography (CT-FLT), CT provides the anatomical information to generate scattering and absorption maps supporting a 3-dimensional reconstruction from the raw optical data. However, given the CT's limited soft tissue contrast, fluorescence reconstruction and quantification can be inaccurate and not sufficiently detailed. Magnetic resonance imaging (MRI) can overcome these limitations and extend the options for tissue characterization. Thus, we aimed to establish a hybrid CT-MRI-FLT approach for whole-body imaging and compared it with CT-FLT.

MATERIALS AND METHODS

The MRI-based hybrid imaging approaches were established first by scanning a water and coconut oil-filled phantom, second by quantifying Cy7 concentrations of inserts in dead mice, and finally by analyzing the biodistribution of AF750-labeled immunoglobulins (IgG, IgA) in living SKH1 mice. Magnetic resonance imaging, acquired with a fat-water-separated mDixon sequence, CT, and FLT were co-registered using markers in the mouse holder frame filled with white petrolatum, which was solid, stable, and visible in both modalities.

RESULTS

Computed tomography-MRI fusion was confirmed by comparing the segmentation agreement using Dice scores. Phantom segmentations showed good agreement, after correction for gradient linearity distortion and chemical shift. Organ segmentations in dead and living mice revealed adequate agreement for fusion. Marking the mouse holder frame and the successful CT-MRI fusion enabled MRI-FLT as well as CT-MRI-FLT reconstructions. Fluorescence tomography reconstructions supported by CT, MRI, or CT-MRI were comparable in dead mice with 60 pmol fluorescence inserts at different locations. Although standard CT-FLT reconstruction only considered general values for soft tissue, skin, lung, fat, and bone scattering, MRI's more versatile soft tissue contrast enabled the additional consideration of liver, kidneys, and brain. However, this did not change FLT reconstructions and quantifications significantly, whereas for extending scattering maps, it was important to accurately segment the organs and the entire mouse body. The various FLT reconstructions also provided comparable results for the in vivo biodistribution analyses with fluorescent immunoglobulins. However, MRI additionally enabled the visualization of gallbladder, thyroid, and brain. Furthermore, segmentations of liver, spleen, and kidney were more reliable due to better-defined contours than in CT. Therefore, the improved segmentations enabled better assignment of fluorescence signals and more differentiated conclusions with MRI-FLT.

CONCLUSIONS

Whole-body CT-MRI-FLT was implemented as a novel trimodal imaging approach, which allowed to more accurately assign fluorescence signals, thereby significantly improving pharmacokinetic analyses.

Future trends for patient-specific dosimetry methodology in molecular radiotherapy

Molecular radiotherapy is rapidly expanding, and new radiotherapeutics are emerging. The majority of treatments is still performed using empirical fixed activities and not tailored for individual patients. Molecular radiotherapy dosimetry is often seen as a promising candidate that would allow personalisation of treatments as outcome should ultimately depend on the absorbed doses delivered and not the activities administered. The field of molecular radiotherapy dosimetry has made considerable progress towards the feasibility of routine clinical dosimetry with reasonably accurate absorbed-dose estimates for a range of molecular radiotherapy dosimetry applications. A range of challenges remain with respect to the accurate quantification, assessment of time-integrated activity and absorbed dose estimation. In this review, we summarise a range of technological and methodological advancements, mainly focussed on beta-emitting molecular radiotherapeutics, that aim to improve molecular radiotherapy dosimetry to achieve accurate, reproducible, and streamlined dosimetry. We describe how these new technologies can potentially improve the often time-consuming considered process of dosimetry and provide suggestions as to what further developments might be required.

A Projection-Domain Low-Count Quantitative SPECT Method for α-Particle-Emitting Radiopharmaceutical Therapy

Single-photon emission-computed tomography (SPECT) provides a mechanism to estimate regional isotope uptake in lesions and at-risk organs after administration of α-particle-emitting radiopharmaceutical therapies (α-RPTs). However, this estimation task is challenging due to the complex emission spectra, the very low number of detected counts (~20 times lower than in conventional SPECT), the impact of stray-radiation-related noise at these low counts, and the multiple image-degrading processes in SPECT. The conventional reconstruction-based quantification methods are observed to be erroneous for α-RPT SPECT. To address these challenges, we developed a low-count quantitative SPECT (LC-QSPECT) method that directly estimates the regional activity uptake from the projection data (obviating the reconstruction step), compensates for stray-radiation-related noise, and accounts for the radioisotope and SPECT physics, including the isotope spectra, scatter, attenuation, and collimator-detector response, using a Monte Carlo-based approach. The method was validated in the context of 3-D SPECT with 223Ra, a commonly used radionuclide for α-RPT. Validation was performed using both realistic simulation studies, including a virtual clinical trial, and synthetic and 3-D-printed anthropomorphic physical-phantom studies. Across all studies, the LC-QSPECT method yielded reliable regional-uptake estimates and outperformed the conventional ordered subset expectation-maximization (OSEM)-based reconstruction and geometric transfer matrix (GTM)-based post-reconstruction partial-volume compensation methods. Furthermore, the method yielded reliable uptake across different lesion sizes, contrasts, and different levels of intralesion heterogeneity. Additionally, the variance of the estimated uptake approached the Cramér-Rao bound-defined theoretical limit. In conclusion, the proposed LC-QSPECT method demonstrated the ability to perform reliable quantification for α-RPT SPECT.

A Narrative Review of 99mTc-Aprotinin in the Diagnosis of Cardiac Amyloidosis and a New Life for an Unfairly Abandoned Drug

Several studies investigated the use of 99mTc-labelled Aprotinin as an amyloid seeker some years ago. In vitro tests showed high binding affinity for several types of amyloid fibrils accompanied by an excellent specificity. Initial human studies demonstrated good accuracy in detecting cardiac involvement. Scintigraphy results were confirmed in a group of 28 endomyocardial biopsies. Unfortunately, clinical studies were halted because of a temporary suspension of the vector protein (Trasylol) and public health concerns over prion contamination of the bovine origin compound. To obviate these limitations, efforts have been made to label a recombinant Aprotinin with 99mTc, which exhibits the same affinity for h-insulin fibrils. With the aim of developing a PET tracer, the same recombinant protein was labeled with Gallium. The introduction of a bifunctional chelator did not affect fibril affinity. Finally, a synthetic peptidic fragment, the cyclic 30-51 SS, was synthetized. After direct technetium labeling, an impressive increase in affinity was demonstrated. This peptide appears to be a potential candidate for Gallium labeling through a bifunctional chelator for PET imaging.

Antiproliferative and antioxidant properties of nematocysts crude venom from jellyfish Acromitus flagellatus against human cancer cell lines

This study aimed to investigate the antiproliferative and antioxidant properties of crude venom from the nematocyst of Jellyfish Acromitus flagellates on human lung cancer (A549) and liver cancer (HepG2) cell lines. The prepared crude venom was subjected to analyses of the biochemical constituents, protein profiles, antioxidant and anticancer activities by standard methods. The extracted venom was pale-yellow in color and viscous/sticky. The biochemical composition such as, protein (1.547 mg/ml), lipid (0.039 mg/ml) and carbohydrate (0.028 mg/ml) was estimated. Protein profiles were determined by SDS PAGE, the result revealed that the molecular weight range from 205 - 3.5 kDa. The free radical scavenging activity was analyzed by the reducing potential (56.36%), DPPH (72.47%), hydroxyl (68.50%), superoxide anion (65.75%), and nitric oxide (33.04%). The cell viability was observed by using different concentrations (20 to 100 µg/ml) of crude venom on A549 and HepG2 cancer cell lines and the IC50 values were recorded in (60 μg/ml and 40 μg/ml) respectively, while it had none cytotoxic effects on Vero cell line up to the concentration of 90 μg/ml. These results suggest that crude venom from nematocyst of A. flagellatus possesses anti-cancer activity and able to develop novel drugs on marine-derived compounds.

Comparison of Correction Techniques for the Spill in Effect in Emission Tomography

In positron emission tomography (PET) imaging, accurate clinical assessment is often affected by the partial volume effect (PVE) leading to overestimation (spill-in) or underestimation (spill-out) of activity in various small regions. The spill-in correction, in particular, can be very challenging when the target region is close to a hot background region. Therefore, this study evaluates and compares the performance of various recently developed spill-in correction techniques, namely: background correction (BC), local projection (LP), and hybrid kernelized (HKEM) methods. We used a simulated digital phantom and [¹⁸F]-NaF PET data of three patients with abdominal aortic aneurysms (AAA) acquired with Siemens Biograph mMR™ and mCT™ scanners respectively. Region of Interest (ROI) analysis was performed and the extracted SUV _mean , SUV _max and target-to-background ratio (TBR) scores were compared. Results showed substantial spill-in effects from hot regions to targeted regions, which are more prominent in small structures. The phantom experiment demonstrated the feasibility of spill-in correction with all methods. For the patient data, large differences in SUV _mean , SUV _max and TBR _max scores were observed between the ROIs drawn over the entire aneurysm and ROIs excluding some regions close to the bone. Overall, BC yielded the best performance in spill-in correction in both phantom and patient studies.

AI-based applications in hybrid imaging: how to build smart and truly multi-parametric decision models for radiomics

INTRODUCTION

The quantitative imaging features (radiomics) that can be obtained from the different modalities of current-generation hybrid imaging can give complementary information with regard to the tumour environment, as they measure different morphologic and functional imaging properties. These multi-parametric image descriptors can be combined with artificial intelligence applications into predictive models. It is now the time for hybrid PET/CT and PET/MRI to take the advantage offered by radiomics to assess the added clinical benefit of using multi-parametric models for the personalized diagnosis and prognosis of different disease phenotypes.

OBJECTIVE

The aim of the paper is to provide an overview of current challenges and available solutions to translate radiomics into hybrid PET-CT and PET-MRI imaging for a smart and truly multi-parametric decision model.

Partial volume correction for improved PET quantification in 18F-NaF imaging of atherosclerotic plaques

BACKGROUND

Accurate quantification of plaque imaging using 18F-NaF PET requires partial volume correction (PVC).

METHODS

PVC of PET data was implemented by the use of a local projection (LP) method. LP-based PVC was evaluated with an image quality (NEMA) and with a thorax phantom with "plaque-type" lesions of 18-36 mL. The validated PVC method was then applied to a cohort of 17 patients, each with at least one plaque in the carotid or ascending aortic arteries. In total, 51 calcified (HU > 110) and 16 non-calcified plaque lesions (HU < 110) were analyzed. The lesion-to-background ratio (LBR) and the relative change of LBR (ΔLBR) were measured on PET.

RESULTS

Following PVC, LBR of the spheres (NEMA phantom) was within 10% of the original values. LBR of the thoracic lesions increased by 155% to 440% when the LP-PVC method was applied to the PET images. In patients, PVC increased the LBR in both calcified [mean = 78% (-8% to 227%)] and non-calcified plaques [mean = 41%, (-9%-104%)].

CONCLUSIONS

PVC helps to improve LBR of plaque-type lesions in both phantom studies and clinical patients. Better results were obtained when the PVC method was applied to images reconstructed with point spread function modeling.

Impact of motion compensation and partial volume correction for 18F-NaF PET/CT imaging of coronary plaque

Recent studies have suggested that ¹⁸F-NaF-PET enables visualization and quantification of plaque micro-calcification in the coronary tree. However, PET imaging of plaque calcification in the coronary arteries is challenging because of the respiratory and cardiac motion as well as partial volume effects. The objective of this work is to implement an image reconstruction framework, which incorporates compensation for respiratory as well as cardiac motion (MoCo) and partial volume correction (PVC), for cardiac ¹⁸F-NaF PET imaging in PET/CT. We evaluated the effect of MoCo and PVC on the quantification of vulnerable plaques in the coronary arteries. Realistic simulations (Biograph TPTV, Biograph mCT) and phantom acquisitions (Biograph mCT) were used for these evaluations. Different uptake values in the calcified plaques were evaluated in the simulations, while three 'plaque-type' lesions of 36, 31 and 18 mm³ were included in the phantom experiments. After validation, the MoCo and PVC methods were applied in four pilot NaF-PET patient studies. In all cases, the MoCo-based image reconstruction was performed using the STIR software. The PVC was obtained from a local projection (LP) method, previously evaluated in preclinical and clinical PET. The results obtained show a significant increase of the measured lesion-to-background ratios (LBR) in the MoCo  +  PVC images. These ratios were further enhanced when using directly the tissue-activities from the LP method, making this approach more suitable for the quantitative evaluation of coronary plaques. When using the LP method on the MoCo images, LBR increased between 200% and 1119% in the simulated data, between 212% and 614% in the phantom experiments and between 46% and 373% in the plaques with positive uptake observed in the pilot patients. In conclusion, we have built and validated a STIR framework incorporating MoCo and PVC for ¹⁸F-NaF PET imaging of coronary plaques. First results indicate an improved quantification of plaque-type lesions.

Noise suppressed partial volume correction for cardiac SPECT/CT

PURPOSE

Partial volume correction (PVC) methods typically improve quantification at the expense of increased image noise and reduced reproducibility. In this study, the authors developed a novel voxel-based PVC method that incorporates anatomical knowledge to improve quantification while suppressing noise for cardiac SPECT/CT imaging.

METHODS

In the proposed method, the SPECT images were first reconstructed using anatomical-based maximum a posteriori (AMAP) with Bowsher's prior to penalize noise while preserving boundaries. A sequential voxel-by-voxel PVC approach (Yang's method) was then applied on the AMAP reconstruction using a template response. This template response was obtained by forward projecting a template derived from a contrast-enhanced CT image, and then reconstructed using AMAP to model the partial volume effects (PVEs) introduced by both the system resolution and the smoothing applied during reconstruction. To evaluate the proposed noise suppressed PVC (NS-PVC), the authors first simulated two types of cardiac SPECT studies: a (99m)Tc-tetrofosmin myocardial perfusion scan and a (99m)Tc-labeled red blood cell (RBC) scan on a dedicated cardiac multiple pinhole SPECT/CT at both high and low count levels. The authors then applied the proposed method on a canine equilibrium blood pool study following injection with (99m)Tc-RBCs at different count levels by rebinning the list-mode data into shorter acquisitions. The proposed method was compared to MLEM reconstruction without PVC, two conventional PVC methods, including Yang's method and multitarget correction (MTC) applied on the MLEM reconstruction, and AMAP reconstruction without PVC.

RESULTS

The results showed that the Yang's method improved quantification, however, yielded increased noise and reduced reproducibility in the regions with higher activity. MTC corrected for PVE on high count data with amplified noise, although yielded the worst performance among all the methods tested on low-count data. AMAP effectively suppressed noise and reduced the spill-in effect in the low activity regions. However it was unable to reduce the spill-out effect in high activity regions. NS-PVC yielded superior performance in terms of both quantitative assessment and visual image quality while improving reproducibility.

CONCLUSIONS

The results suggest that NS-PVC may be a promising PVC algorithm for application in low-dose protocols, and in gated and dynamic cardiac studies with low counts.

Development of computational small animal models and their applications in preclinical imaging and therapy research

The development of multimodality preclinical imaging techniques and the rapid growth of realistic computer simulation tools have promoted the construction and application of computational laboratory animal models in preclinical research. Since the early 1990s, over 120 realistic computational animal models have been reported in the literature and used as surrogates to characterize the anatomy of actual animals for the simulation of preclinical studies involving the use of bioluminescence tomography, fluorescence molecular tomography, positron emission tomography, single-photon emission computed tomography, microcomputed tomography, magnetic resonance imaging, and optical imaging. Other applications include electromagnetic field simulation, ionizing and nonionizing radiation dosimetry, and the development and evaluation of new methodologies for multimodality image coregistration, segmentation, and reconstruction of small animal images. This paper provides a comprehensive review of the history and fundamental technologies used for the development of computational small animal models with a particular focus on their application in preclinical imaging as well as nonionizing and ionizing radiation dosimetry calculations. An overview of the overall process involved in the design of these models, including the fundamental elements used for the construction of different types of computational models, the identification of original anatomical data, the simulation tools used for solving various computational problems, and the applications of computational animal models in preclinical research. The authors also analyze the characteristics of categories of computational models (stylized, voxel-based, and boundary representation) and discuss the technical challenges faced at the present time as well as research needs in the future.

Improved quantification for local regions of interest in preclinical PET imaging

In Positron Emission Tomography, there are several causes of quantitative inaccuracy, such as partial volume or spillover effects. The impact of these effects is greater when using radionuclides that have a large positron range, e.g. (68)Ga and (124)I, which have been increasingly used in the clinic. We have implemented and evaluated a local projection algorithm (LPA), originally evaluated for SPECT, to compensate for both partial-volume and spillover effects in PET. This method is based on the use of a high-resolution CT or MR image, co-registered with a PET image, which permits a high-resolution segmentation of a few tissues within a volume of interest (VOI) centered on a region within which tissue-activity values need to be estimated. The additional boundary information is used to obtain improved activity estimates for each tissue within the VOI, by solving a simple inversion problem. We implemented this algorithm for the preclinical Argus PET/CT scanner and assessed its performance using the radionuclides (18)F, (68)Ga and (124)I. We also evaluated and compared the results obtained when it was applied during the iterative reconstruction, as well as after the reconstruction as a postprocessing procedure. In addition, we studied how LPA can help to reduce the 'spillover contamination', which causes inaccurate quantification of lesions in the immediate neighborhood of large, 'hot' sources. Quantification was significantly improved by using LPA, which provided more accurate ratios of lesion-to-background activity concentration for hot and cold regions. For (18)F, the contrast was improved from 3.0 to 4.0 in hot lesions (when the true ratio was 4.0) and from 0.16 to 0.06 in cold lesions (true ratio  =  0.0), when using the LPA postprocessing. Furthermore, activity values estimated within the VOI using LPA during reconstruction were slightly more accurate than those obtained by post-processing, while also visually improving the image contrast and uniformity within the VOI.

Comparison of the scanning linear estimator (SLE) and ROI methods for quantitative SPECT imaging

In quantitative emission tomography, tumor activity is typically estimated from calculations on a region of interest (ROI) identified in the reconstructed slices. In these calculations, unpredictable bias arising from the null functions of the imaging system affects ROI estimates. The magnitude of this bias depends upon the tumor size and location. In prior work it has been shown that the scanning linear estimator (SLE), which operates on the raw projection data, is an unbiased estimator of activity when the size and location of the tumor are known. In this work, we performed analytic simulation of SPECT imaging with a parallel-hole medium-energy collimator. Distance-dependent system spatial resolution and non-uniform attenuation were included in the imaging simulation. We compared the task of activity estimation by the ROI and SLE methods for a range of tumor sizes (diameter: 1-3 cm) and activities (contrast ratio: 1-10) added to uniform and non-uniform liver backgrounds. Using the correct value for the tumor shape and location is an idealized approximation to how task estimation would occur clinically. Thus we determined how perturbing this idealized prior knowledge impacted the performance of both techniques. To implement the SLE for the non-uniform background, we used a novel iterative algorithm for pre-whitening stationary noise within a compact region. Estimation task performance was compared using the ensemble mean-squared error (EMSE) as the criterion. The SLE method performed substantially better than the ROI method (i.e. EMSE(SLE) was 23-174 times lower) when the background is uniform and tumor location and size are known accurately. The variance of the SLE increased when a non-uniform liver texture was introduced but the EMSE(SLE) continued to be 5-20 times lower than the ROI method. In summary, SLE outperformed ROI under almost all conditions that we tested.

Partial volume effect estimation and correction in the aortic vascular wall in PET imaging

We evaluated the impact of partial volume effect (PVE) in the assessment of arterial diseases with (18)FDG PET. An anthropomorphic digital phantom enabling the modeling of aorta related diseases like atherosclerosis and arteritis was used. Based on this phantom, we performed GATE Monte Carlo simulations to produce realistic PET images with a known organ segmentation and ground truth activity values. Images corresponding to 15 different activity-concentration ratios between the aortic wall and the blood and to 7 different wall thicknesses were generated. Using the PET images, we compared the theoretical wall-to-blood activity-concentration ratios (WBRs) with the measured WBRs obtained with five measurement methods: (1) measurement made by a physician (Expert), (2) automated measurement supposed to mimic the physician measurements (Max), (3) simple correction based on a recovery coefficient (Max-RC), (4) measurement based on an ideal VOI segmentation (Mean-VOI) and (5) measurement corrected for PVE using an ideal geometric transfer matrix (GTM) method. We found that Mean-VOI WBRs values were strongly affected by PVE. WBRs obtained by the physician measurement, by the Max method and by the Max-RC method were more accurate than WBRs obtained with the Mean-VOI approach. However Expert, Max and Max-RC WBRs strongly depended on the wall thickness. Only the GTM corrected WBRs did not depend on the wall thickness. Using the GTM method, we obtained more reproducible ratio values that could be compared across wall thickness. Yet, the feasibility of the implementation of a GTM-like method on real data remains to be studied.

Least-squares dual characterization for ROI assessment in emission tomography

Our aim is to describe an original method for estimating the statistical properties of regions of interest (ROIs) in emission tomography. Drawn upon the works of Louis on the approximate inverse, we propose a dual formulation of the ROI estimation problem to derive the ROI activity and variance directly from the measured data without any image reconstruction. The method requires the definition of an ROI characteristic function that can be extracted from a co-registered morphological image. This characteristic function can be smoothed to optimize the resolution-variance tradeoff. An iterative procedure is detailed for the solution of the dual problem in the least-squares sense (least-squares dual (LSD) characterization), and a linear extrapolation scheme is described to compensate for sampling partial volume effect and reduce the estimation bias (LSD-ex). LSD and LSD-ex are compared with classical ROI estimation using pixel summation after image reconstruction and with Huesman's method. For this comparison, we used Monte Carlo simulations (GATE simulation tool) of 2D PET data of a Hoffman brain phantom containing three small uniform high-contrast ROIs and a large non-uniform low-contrast ROI. Our results show that the performances of LSD characterization are at least as good as those of the classical methods in terms of root mean square (RMS) error. For the three small tumor regions, LSD-ex allows a reduction in the estimation bias by up to 14%, resulting in a reduction in the RMS error of up to 8.5%, compared with the optimal classical estimation. For the large non-specific region, LSD using appropriate smoothing could intuitively and efficiently handle the resolution-variance tradeoff.

Evaluation of a method for projection-based tissue-activity estimation within small volumes of interest

A new method of compensating for tissue-fraction and count-spillover effects, which require tissue segmentation only within a small volume surrounding the primary lesion of interest, was evaluated for SPECT imaging. Tissue-activity concentration estimates are obtained by fitting the measured projection data to a statistical model of the segmented tissue projections. Multiple realizations of two simulated human-torso phantoms, each containing 20 spherical 'tumours', 1.6 cm in diameter, with tumour-to-background ratios of 8:1 and 4:1, were simulated. Estimates of tumour- and background-activity concentration values for homogeneous as well as inhomogeneous tissue activities were compared to the standard uptake value (SUV) metrics on the basis of accuracy and precision. For perfectly registered, high-contrast, superficial lesions in a homogeneous background without scatter, the method yielded accurate (<0.4% bias) and precise (<6.1%) recovery of the simulated activity values, significantly outperforming the SUV metrics. Tissue inhomogeneities, greater tumour depth and lower contrast ratios degraded precision (up to 11.7%), but the estimates remained almost unbiased. The method was comparable in accuracy but more precise than a well-established matrix inversion approach, even when errors in tumour size and position were introduced to simulate moderate inaccuracies in segmentation and image registration. Photon scatter in the object did not significantly affect the accuracy or precision of the estimates.