Search

Simple Summary: Radiomics analysis of positron emission tomography (PET) images can provide objective measurements of tumor heterogeneity and spatial patterns. However, the relatively low resolution, high noise, and limited longitudinal data availability make it difficult to systematically investigate the relationship between the microscopic tumor phenotypes and corresponding PET radiomics signatures. To address this challenge, we use a multiscale, stochastic mathematical model of tumor growth to generate cross-sections of tumors in vascularized normal tissue on a microscopic level. By varying the biological parameters of the model, distinct tumor phenotypes are obtained, and their corresponding PET images are generated. The simulated data are then used to find the optimal combination of PET radiomics features that can reliably distinguish visually similar tumor phenotypes. In addition, we study the longitudinal changes in the discriminative power of radiomics features with tumor growth from a single cell to approximately one million cells. Cancers can manifest large variations in tumor phenotypes due to genetic and microenvironmental factors, which has motivated the development of quantitative radiomics-based image analysis with the aim to robustly classify tumor phenotypes in vivo. Positron emission tomography (PET) imaging can be particularly helpful in elucidating the metabolic profiles of tumors. However, the relatively low resolution, high noise, and limited PET data availability make it difficult to study the relationship between the microenvironment properties and metabolic tumor phenotype as seen on the images. Most of previously proposed digital PET phantoms of tumors are static, have an over-simplified morphology, and lack the link to cellular biology that ultimately governs the tumor evolution. In this work, we propose a novel method to investigate the relationship between microscopic tumor parameters and PET image characteristics based on the computational simulation of tumor growth. We use a hybrid, multiscale, stochastic mathematical model of cellular metabolism and proliferation to generate simulated cross-sections of tumors in vascularized normal tissue on a microscopic level. The generated longitudinal tumor growth sequences are converted to PET images with realistic resolution and noise. By changing the biological parameters of the model, such as the blood vessel density and conditions for necrosis, distinct tumor phenotypes can be obtained. The simulated cellular maps were compared to real histology slides of SiHa and WiDr xenografts imaged with Hoechst 33342 and pimonidazole. As an example application of the proposed method, we simulated six tumor phenotypes that contain various amounts of hypoxic and necrotic regions induced by a lack of oxygen and glucose, including phenotypes that are distinct on the microscopic level but visually similar in PET images. We computed 22 standardized Haralick texture features for each phenotype, and identified the features that could best discriminate the phenotypes with varying image noise levels. We demonstrated that "cluster shade" and "difference entropy" are the most effective and noise-resilient features for microscopic phenotype discrimination. Longitudinal analysis of the simulated tumor growth showed that radiomics analysis can be beneficial even in small lesions with a diameter of 3.5–4 resolution units, corresponding to 8.7–10.0 mm in modern PET scanners. Certain radiomics features were shown to change non-monotonically with tumor growth, which has implications for feature selection for tracking disease progression and therapy response. [ABSTRACT FROM AUTHOR]

Background: Prostate‐specific membrane antigen (PSMA) PET imaging represents a valuable source of information reflecting disease stage, response rate, and treatment optimization options, particularly with PSMA radioligand therapy. Quantification of radiopharmaceutical uptake in healthy organs from PSMA images has the potential to minimize toxicity by extrapolation of the radiation dose delivery towards personalization of therapy. However, segmentation and quantification of uptake in organs requires labor‐intensive organ delineations that are often not feasible in the clinic nor scalable for large clinical trials. Purpose: In this work we develop and test the PSMA Healthy organ segmentation network (PSMA‐Hornet), a fully‐automated deep neural net for simultaneous segmentation of 14 healthy organs representing the normal biodistribution of [18F]DCFPyL on PET/CT images. We also propose a modified U‐net architecture, a self‐supervised pre‐training method for PET/CT images, a multi‐target Dice loss, and multi‐target batch balancing to effectively train PSMA‐Hornet and similar networks. Methods: The study used manually‐segmented [18F]DCFPyL PET/CT images from 100 subjects, and 526 similar images without segmentations. The unsegmented images were used for self‐supervised model pretraining. For supervised training, Monte‐Carlo cross‐validation was used to evaluate the network performance, with 85 subjects in each trial reserved for model training, 5 for validation, and 10 for testing. Image segmentation and quantification metrics were evaluated on the test folds with respect to manual segmentations by a nuclear medicine physician, and compared to inter‐rater agreement. The model's segmentation performance was also evaluated on a separate set of 19 images with high tumor load. Results: With our best model, the lowest mean Dice coefficient on the test set was 0.826 for the sublingual gland, and the highest was 0.964 for liver. The highest mean error in tracer uptake quantification was 13.9% in the sublingual gland. Self‐supervised pretraining improved training convergence, train‐to‐test generalization, and segmentation quality. In addition, we found that a multi‐target network produced significantly higher segmentation accuracy than single‐organ networks. Conclusions: The developed network can be used to automatically obtain high‐quality organ segmentations for PSMA image analysis tasks. It can be used to reproducibly extract imaging data, and holds promise for clinical applications such as personalized radiation dosimetry and improved radioligand therapy. [ABSTRACT FROM AUTHOR]

Purpose: Respiratory motion during positron emission tomography (PET) scans can be a major detriment to image quality in oncological imaging. The impact of motion on lesion quantification and detectability can be assessed using phantoms with realistic anatomy representation and motion modeling. In this work, we develop an anthropomorphic phantom for PET imaging that combines anatomic fidelity and a realistic breathing mechanism with deformable lungs. Methods: We start from a previously developed anatomically accurate but static phantom of a human torso, and add elastic lungs with a highly controllable actuation mechanism which replicates the physics of breathing. The space outside the lungs is filled with a radioactive water solution. To maintain anatomical accuracy and realistic gamma ray attenuation in the torso, all motion mechanisms and actuators are positioned outside of the phantom compartment. The actuation mechanism can produce custom respiratory waveforms with breathing rates up to 25 breaths per minute and tidal volumes up to 1200 mL. Results: Several tests were performed to validate the performance of the phantom assembly, in which the phantom was filled with water and given respiratory waveforms to execute. All parts demonstrated expected performance. Force requirements were not exceeded and no leaks were detected, although continued use of the phantom is required to evaluate wear. The motion of the lungs was determined to be within a reasonable realistic range. Conclusions: The full mechanical design is described in this paper, as well as a software application with graphical user interface which was developed to plan and visualize respiratory patterns. Both are available online as open source files. The developed phantom will facilitate future work in evaluating the impact of respiratory motion on lesion quantification and detectability in clinical practice. [ABSTRACT FROM AUTHOR]

Current methods using a single PET scan to detect voxel-level transient dopamine release—using F-test (significance) and cluster size thresholding—have limited detection sensitivity for clusters of release small in size and/or having low release levels. Specifically, simulations show that voxels with release near the peripheries of such clusters are often rejected—becoming false negatives and ultimately distorting the F-distribution of rejected voxels. We suggest a Monte Carlo method that incorporates these two observations into a cost function, allowing erroneously rejected voxels to be accepted under specified criteria. In simulations, the proposed method improves detection sensitivity by up to 50% while preserving the cluster size threshold, or up to 180% when optimizing for sensitivity. A further parametric-based voxelwise thresholding is then suggested to better estimate the release dynamics in detected clusters. We apply the Monte Carlo method to a pilot scan from a human gambling study, where additional parametrically unique clusters are detected as compared to the current best methods—results consistent with our simulations. [ABSTRACT FROM AUTHOR]