Linear scaling computation of forces for the domain-decomposition linear Poisson-Boltzmann method.

The Linearized Poisson-Boltzmann (LPB) equation is a popular and widely accepted model for accounting solvent effects in computational (bio-) chemistry. In the present article, we derive the analytical forces using the domain-decomposition-based LPB-method with a van-der Waals or solvent-accessible surface. We present an efficient strategy to compute the forces and its implementation, allowing linear scaling of the method with respect to the number of atoms using the fast multipole method. Numerical tests illustrate the accuracy of the computation of the analytical forces and compare the efficiency with other available methods.

Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution.

The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) imaging without a separate transmission scan is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Based on this premise, we investigated whether the SPECT scattered-photon data acquired in list-mode (LM) format and including the energy information can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramer-Rao bound (CRB) for the task of jointly estimating the activity and attenuation distribution using only the SPECT emission data. In the process, a path-based formalism to process the LM SPECT emission data, including the scattered-photon data, is proposed. The Fisher information method was implemented on NVIDIA graphics processing units (GPU) for acceleration. The method was applied to analyze the information content of SPECT LM emission data, which contains up to first-order scattered events, in a simulated SPECT system with parameters modeling a clinical system using realistic computational studies with 2-D digital synthetic and anthropomorphic phantoms. The method was also applied to LM data containing up to second-order scatter for a synthetic phantom. Experiments with anthropomorphic phantoms simulated myocardial perfusion and dopamine transporter (DaT)-Scan SPECT studies. The results show that the CRB obtained for the attenuation and activity coefficients was typically much lower than the true value of these coefficients. An increase in the number of detected photons yielded lower CRB for both the attenuation and activity coefficients. Further, we observed that systems with better energy resolution yielded a lower CRB for the attenuation coefficient. Overall, the results provide evidence that LM SPECT emission data, including the scattered photons, contains information to jointly estimate the activity and attenuation coefficients.

Cariostatic Efficacy of Aqueous and Ethanolic Extracts of Liquorice in the Schoolchildren: In Vivo Comparative Study.

AIM: The present study was done to evaluate the in vivo cariostatic efficacy of aqueous and ethanolic extracts of liquorice to ascertain whether it could be developed into a caries-preventive regimen basically targeted for use in the pediatric population. MATERIALS AND METHODS: Thirty schoolchildren of 6-12-year-old were selected for the study. Powder of Glycyrrhiza glabra is used to prepare the gel with various concentration of aqueous and ethanolic liquorice gel. The preweighed dose was delivered through the vials. The drug concentrations were based on their respective minimum bactericidal concentration (MBC) values against Streptococcus mutans, which were calculated earlier. And it is divided into three groups, i.e., group I: aqueous liquorice extract 1.75 g/10 mL saline, group II: ethanolic liquorice extract 350 mg/10 mL, and group III: hexidine (0.2% chlorhexidine, CHX). For statistical analysis, Tukey's post hoc with one-way analysis of variance (ANOVA) and t test were applied. RESULTS: It was found that hexigel has a potential antibacterial activity against S. mutans, with minimum inhibitory concentration (MIC) of 3.14 ± 2.02. Ethanolic liquorice shows MIC of 2.15 ± 0.91 and aqueous liquorice shows MIC of 1.30 ± 1.08. Tukey's post hoc test showed statistically nonsignificant difference between hexigel and ethanolic liquorice against S. mutans. CONCLUSION: On conclusion, the present study found that hexigel was better than both the ethanolic and aqueous solutions of liquorice. And ethanolic liquorice was found to be better than aqueous solution, but it was not statistically significant, which could be due to the small sample size. CLINICAL SIGNIFICANCE: Dental caries is one of the most common infectious microbial diseases. Various steps have been taken to prevent dental caries, fluoride being the most common among them. Nowadays, G. glabra, commonly known as liquorice (mulethi), is one such medicinal plant used by various cultures for thousands of years to relieve coughs, sore throats, and gastric inflammation. This drug in our study demonstrated inhibitory effect on the growth of S. mutans.

A maximum-likelihood method to estimate a single ADC value of lesions using diffusion MRI.

PURPOSE: Design a statistically rigorous procedure to estimate a single apparent diffusion coefficient (ADC) of lesion from the mean lesion signal intensity in diffusion MRI. THEORY AND METHODS: A rigorous maximum-likelihood technique that incorporated the statistics of the mean lesion intensity and accounted for lesion heterogeneity was derived to estimate the ADC value. Performance evaluation included comparison with the conventionally used linear-regression and a statistically rigorous state-of-the-art ADC-map technique using realistic and clinically relevant simulation studies conducted with assistance of patient data for homogeneous and heterogeneous lesion models. RESULTS: The proposed technique outperformed the linear-regression and ADC-map approaches over a large spectrum of signal-to-noise ratio, ADC, lesion size, image-misalignment parameters, including at no image misalignment, and different amounts of lesion heterogeneity. The method was also superior at different sets of b values and in studies from specific patient-image-derived data. The technique took less than a second to execute. CONCLUSIONS: A rigorous, computationally fast, easy-to-implement, and convenient-to-use maximum-likelihood technique was proposed to estimate a single ADC value of the lesion. Results provide strong evidence in support of the method. Magn Reson Med 76:1919-1931, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Impact of the Fano Factor on Position and Energy Estimation in Scintillation Detectors.

The Fano factor for an integer-valued random variable is defined as the ratio of its variance to its mean. Light from various scintillation crystals have been reported to have Fano factors from sub-Poisson (Fano factor < 1) to super-Poisson (Fano factor > 1). For a given mean, a smaller Fano factor implies a smaller variance and thus less noise. We investigated if lower noise in the scintillation light will result in better spatial and energy resolutions. The impact of Fano factor on the estimation of position of interaction and energy deposited in simple gamma-camera geometries is estimated by two methods - calculating the Cramér-Rao bound and estimating the variance of a maximum likelihood estimator. The methods are consistent with each other and indicate that when estimating the position of interaction and energy deposited by a gamma-ray photon, the Fano factor of a scintillator does not affect the spatial resolution. A smaller Fano factor results in a better energy resolution.

How accurately can quantitative imaging methods be ranked without ground truth: An upper bound on no-gold-standard evaluation.

Objective evaluation of quantitative imaging (QI) methods with patient data, while important, is typically hindered by the lack of gold standards. To address this challenge, no-gold-standard evaluation (NGSE) techniques have been proposed. These techniques have demonstrated efficacy in accurately ranking QI methods without access to gold standards. The development of NGSE methods has raised an important question: how accurately can QI methods be ranked without ground truth. To answer this question, we propose a Cramer-Rao bound (CRB)-based framework that quantifies the upper bound in ranking QI methods without any ground truth. We present the application of this framework in guiding the use of a well-known NGSE technique, namely the regression-without-truth (RWT) technique. Our results show the utility of this framework in quantifying the performance of this NGSE technique for different patient numbers. These results provide motivation towards studying other applications of this upper bound.

Deep learning generation of preclinical positron emission tomography (PET) images from low-count PET with task-based performance assessment.

BACKGROUND: Preclinical low-count positron emission tomography (LC-PET) imaging offers numerous advantages such as facilitating imaging logistics, enabling longitudinal studies of long- and short-lived isotopes as well as increasing scanner throughput. However, LC-PET is characterized by reduced photon-count levels resulting in low signal-to-noise ratio (SNR), segmentation difficulties, and quantification uncertainties. PURPOSE: We developed and evaluated a novel deep-learning (DL) architecture-Attention based Residual-Dilated Net (ARD-Net)-to generate standard-count PET (SC-PET) images from LC-PET images. The performance of the ARD-Net framework was evaluated for numerous low count realizations using fidelity-based qualitative metrics, task-based segmentation, and quantitative metrics. METHOD: Patient Derived tumor Xenograft (PDX) with tumors implanted in the mammary fat-pad were subjected to preclinical [18F]-Fluorodeoxyglucose (FDG)-PET/CT imaging. SC-PET images were derived from a 10 min static FDG-PET acquisition, 50 min post administration of FDG, and were resampled to generate four distinct LC-PET realizations corresponding to 10%, 5%, 1.6%, and 0.8% of SC-PET count-level. ARD-Net was trained and optimized using 48 preclinical FDG-PET datasets, while 16 datasets were utilized to assess performance. Further, the performance of ARD-Net was benchmarked against two leading DL-based methods (Residual UNet, RU-Net; and Dilated Network, D-Net) and non-DL methods (Non-Local Means, NLM; and Block Matching 3D Filtering, BM3D). The performance of the framework was evaluated using traditional fidelity-based image quality metrics such as Structural Similarity Index Metric (SSIM) and Normalized Root Mean Square Error (NRMSE), as well as human observer-based tumor segmentation performance (Dice Score and volume bias) and quantitative analysis of Standardized Uptake Value (SUV) measurements. Additionally, radiomics-derived features were utilized as a measure of quality assurance (QA) in comparison to true SC-PET. Finally, a performance ensemble score (EPS) was developed by integrating fidelity-based and task-based metrics. Concordance Correlation Coefficient (CCC) was utilized to determine concordance between measures. The non-parametric Friedman Test with Bonferroni correction was used to compare the performance of ARD-Net against benchmarked methods with significance at adjusted p-value ≤0.01. RESULTS: ARD-Net-generated SC-PET images exhibited significantly better (p ≤ 0.01 post Bonferroni correction) overall image fidelity scores in terms of SSIM and NRMSE at majority of photon-count levels compared to benchmarked DL and non-DL methods. In terms of task-based quantitative accuracy evaluated by SUVMean and SUVPeak, ARD-Net exhibited less than 5% median absolute bias for SUVMean compared to true SC-PET and lower degree of variability compared to benchmarked DL and non-DL based methods in generating SC-PET. Additionally, ARD-Net-generated SC-PET images displayed higher degree of concordance to SC-PET images in terms of radiomics features compared to non-DL and other DL approaches. Finally, the ensemble score suggested that ARD-Net exhibited significantly superior performance compared to benchmarked algorithms (p ≤ 0.01 post Bonferroni correction). CONCLUSION: ARD-Net provides a robust framework to generate SC-PET from LC-PET images. ARD-Net generated SC-PET images exhibited superior performance compared other DL and non-DL approaches in terms of image-fidelity based metrics, task-based segmentation metrics, and minimal bias in terms of task-based quantification performance for preclinical PET imaging.

WIN-PDQ: A Wiener-estimator-based projection-domain quantitative SPECT method that accounts for intra-regional uptake heterogeneity.

SPECT can enable the quantification of activity uptake in lesions and at-risk organs in {\alpha}-particle-emitting radiopharmaceutical therapies ({\alpha}-RPTs). But this quantification is challenged by the low photon counts, complicated isotope physics, and the image-degrading effects in {\alpha}-RPT SPECT. Thus, strategies to optimize the SPECT system and protocol designs for the task of regional uptake quantification are needed. Objectively performing this task-based optimization requires a reliable (accurate and precise) regional uptake quantification method. Conventional reconstruction-based quantification (RBQ) methods have been observed to be erroneous for {\alpha}-RPT SPECT. Projection-domain quantification methods, which estimate regional uptake directly from SPECT projections, have demonstrated potential in providing reliable regional uptake estimates, but these methods assume constant uptake within the regions, an assumption that may not hold. To address these challenges, we propose WIN-PDQ, a Wiener-estimator-based projection-domain quantitative SPECT method. The method accounts for the heterogeneity within the regions of interest while estimating mean uptake. An early-stage evaluation of the method was conducted using 3D Monte Carlo-simulated SPECT of anthropomorphic phantoms with radium-223 uptake and lumpy-model-based intra-regional uptake heterogeneity. In this evaluation with phantoms of varying mean regional uptake and intra-regional uptake heterogeneity, the WIN-PDQ method yielded ensemble unbiased estimates and significantly outperformed both reconstruction-based and previously proposed projection-domain quantification methods. In conclusion, based on these preliminary findings, the proposed method is showing potential for estimating mean regional uptake in {\alpha}-RPTs and towards enabling the objective task-based optimization of SPECT system and protocol designs.

Joint regional uptake quantification of Thorium-227 and Radium-223 using a multiple-energy-window projection-domain quantitative SPECT method.

Thorium-227-based alpha-particle radiopharmaceutical therapies ({\alpha}-RPTs) are being investigated in several clinical and pre-clinical studies. After administration, Thorium-227 decays to Radium-223, another alpha-particle-emitting isotope, which redistributes within the patient. Reliable dose quantification of both Thorium-227 and Radium-223 is clinically important, and SPECT may perform this quantification as these isotopes also emit X- and gamma-ray photons. However, reliable quantification is challenged by the orders-of-magnitude lower activity compared to conventional SPECT, resulting in a very low number of detected counts, the presence of multiple photopeaks, substantial overlap in the emission spectra of these isotopes, and the image-degrading effects in SPECT. To address these issues, we propose a multiple-energy-window projection-domain quantification (MEW-PDQ) method that jointly estimates the regional activity uptake of both Thorium-227 and Radium-223 directly using the SPECT projection from multiple energy windows. We evaluated the method with realistic simulation studies using anthropomorphic digital phantoms, in the context of imaging patients with bone metastases of prostate cancer and treated with Thorium-227-based {\alpha}-RPTs. The proposed method yielded reliable (accurate and precise) regional uptake estimates of both isotopes and outperformed state-of-the-art methods across different lesion sizes and contrasts, in a virtual imaging trial, as well as with moderate levels of intra-regional heterogeneous uptake and with moderate inaccuracies in the definitions of the support of various regions. Additionally, we demonstrated the effectiveness of using multiple energy windows and the variance of the estimated uptake using the proposed method approached the Cram\'er-Rao-lower-bound-defined theoretical limit.

Joint regional uptake quantification of thorium-227 and radium-223 using a multiple-energy-window projection-domain quantitative SPECT method.

Thorium-227 (227Th)-based α-particle radiopharmaceutical therapies (α-RPTs) are currently being investigated in several clinical and pre-clinical studies. After administration, 227Th decays to 223Ra, another α-particle-emitting isotope, which redistributes within the patient. Reliable dose quantification of both 227Th and 223Ra is clinically important, and SPECT may perform this quantification as these isotopes also emit X- and γ-ray photons. However, reliable quantification is challenging for several reasons: the orders-of-magnitude lower activity compared to conventional SPECT, resulting in a very low number of detected counts, the presence of multiple photopeaks, substantial overlap in the emission spectra of these isotopes, and the image-degrading effects in SPECT. To address these issues, we propose a multiple-energy-window projection-domain quantification (MEW-PDQ) method that jointly estimates the regional activity uptake of both 227Th and 223Ra directly using the SPECT projection data from multiple energy windows. We evaluated the method with realistic simulation studies conducted with anthropomorphic digital phantoms, including a virtual imaging trial, in the context of imaging patients with bone metastases of prostate cancer who were treated with 227Th-based α-RPTs. The proposed method yielded reliable (accurate and precise) regional uptake estimates of both isotopes and outperformed state-of-the-art methods across different lesion sizes and contrasts, as well as in the virtual imaging trial. This reliable performance was also observed with moderate levels of intra-regional heterogeneous uptake as well as when there were moderate inaccuracies in the definitions of the support of various regions. Additionally, we demonstrated the effectiveness of using multiple energy windows and the variance of the estimated uptake using the proposed method approached the Cramér-Rao-lower-bound-defined theoretical limit. These results provide strong evidence in support of this method for reliable uptake quantification in 227Th-based α-RPTs.

ISIT-QA: In Silico Imaging Trial to Evaluate a Low-Count Quantitative SPECT Method Across Multiple Scanner-Collimator Configurations for 223Ra-Based Radiopharmaceutical Therapies.

Personalized dose-based treatment planning requires accurate and reproducible noninvasive measurements to ensure safety and effectiveness. Dose estimation using SPECT is possible but challenging for alpha (α)-particle-emitting radiopharmaceutical therapy (α-RPT) because of complex γ-emission spectra, extremely low counts, and various image-degrading artifacts across a plethora of scanner-collimator configurations. Through the incorporation of physics-based considerations and skipping of the potentially lossy voxel-based reconstruction step, a recently developed projection-domain low-count quantitative SPECT (LC-QSPECT) method has the potential to provide reproducible, accurate, and precise activity concentration and dose measures across multiple scanners, as is typically the case in multicenter settings. To assess this potential, we conducted an in silico imaging trial to evaluate the LC-QSPECT method for a 223Ra-based α-RPT, with the trial recapitulating patient and imaging system variabilities. Methods: A virtual imaging trial titled In Silico Imaging Trial for Quantitation Accuracy (ISIT-QA) was designed with the objectives of evaluating the performance of the LC-QSPECT method across multiple scanner-collimator configurations and comparing performance with a conventional reconstruction-based quantification method. In this trial, we simulated 280 realistic virtual patients with bone-metastatic castration-resistant prostate cancer treated with 223Ra-based α-RPT. The trial was conducted with 9 simulated SPECT scanner-collimator configurations. The primary objective of this trial was to evaluate the reproducibility of dose estimates across multiple scanner-collimator configurations using LC-QSPECT by calculating the intraclass correlation coefficient. Additionally, we compared the reproducibility and evaluated the accuracy of both considered quantification methods across multiple scanner-collimator configurations. Finally, the repeatability of the methods was evaluated in a test-retest study. Results: In this trial, data from 268 223RaCl2 treated virtual prostate cancer patients, with a total of 2,903 lesions, were used to evaluate LC-QSPECT. LC-QSPECT provided dose estimates with good reproducibility across the 9 scanner-collimator configurations (intraclass correlation coefficient > 0.75) and high accuracy (ensemble average values of recovery coefficients ranged from 1.00 to 1.02). Compared with conventional reconstruction-based quantification, LC-QSPECT yielded significantly improved reproducibility across scanner-collimator configurations, accuracy, and test-retest repeatability ([Formula: see text] Conclusion: LC-QSPECT provides reproducible, accurate, and repeatable dose estimations in 223Ra-based α-RPT as evaluated in ISIT-QA. These findings provide a strong impetus for multicenter clinical evaluations of LC-QSPECT in dose quantification for α-RPTs.

DEMIST: A Deep-Learning-Based Detection-Task-Specific Denoising Approach for Myocardial Perfusion SPECT.

There is an important need for methods to process myocardial perfusion imaging (MPI) single-photon emission computed tomography (SPECT) images acquired at lower-radiation dose and/or acquisition time such that the processed images improve observer performance on the clinical task of detecting perfusion defects compared to low-dose images. To address this need, we build upon concepts from model-observer theory and our understanding of the human visual system to propose a detection task-specific deep-learning-based approach for denoising MPI SPECT images (DEMIST). The approach, while performing denoising, is designed to preserve features that influence observer performance on detection tasks. We objectively evaluated DEMIST on the task of detecting perfusion defects using a retrospective study with anonymized clinical data in patients who underwent MPI studies across two scanners (N = 338). The evaluation was performed at low-dose levels of 6.25%, 12.5%, and 25% and using an anthropomorphic channelized Hotelling observer. Performance was quantified using area under the receiver operating characteristics curve (AUC). Images denoised with DEMIST yielded significantly higher AUC compared to corresponding low-dose images and images denoised with a commonly used task-agnostic deep learning-based denoising method. Similar results were observed with stratified analysis based on patient sex and defect type. Additionally, DEMIST improved visual fidelity of the low-dose images as quantified using root mean squared error and structural similarity index metric. A mathematical analysis revealed that DEMIST preserved features that assist in detection tasks while improving the noise properties, resulting in improved observer performance. The results provide strong evidence for further clinical evaluation of DEMIST to denoise low-count images in MPI SPECT.

Can patient-specific acquisition protocol improve performance on defect detection task in myocardial perfusion SPECT?

Myocardial perfusion imaging using single-photon emission computed tomography (SPECT), or myocardial perfusion SPECT (MPS) is a widely used clinical imaging modality for the diagnosis of coronary artery disease. Current clinical protocols for acquiring and reconstructing MPS images are similar for most patients. However, for patients with outlier anatomical characteristics, such as large breasts, images acquired using conventional protocols are often sub-optimal in quality, leading to degraded diagnostic accuracy. Solutions to improve image quality for these patients outside of increased dose or total acquisition time remain challenging. Thus, there is an important need for new methodologies to improve image quality for such patients. One approach to improving this performance is adapting the image acquisition protocol specific to each patient. For this study, we first designed and implemented a personalized patient-specific protocol-optimization strategy, which we term precision SPECT (PRESPECT). This strategy integrates ideal observer theory with the constraints of tomographic reconstruction to optimize the acquisition time for each projection view, such that MPS defect detection performance is maximized. We performed a clinically realistic simulation study on patients with outlier anatomies on the task of detecting perfusion defects on various realizations of low-dose scans by an anthropomorphic channelized Hotelling observer. Our results show that using PRESPECT led to improved performance on the defect detection task for the considered patients. These results provide evidence that personalization of MPS acquisition protocol has the potential to improve defect detection performance, motivating further research to design optimal patient-specific acquisition and reconstruction protocols for MPS, as well as developing similar approaches for other medical imaging modalities.

Adverse events following immunization of ChAdOx1 nCoV-19 vaccine among healthcare workers of a medicine-teaching institution of North India.

Objective: This study sought to assess the prevalence of adverse events following immunization (AEFI) and factors associated with AEFI of the ChAdOx1 nCoV-19 vaccine (Covishield) among healthcare workers (HCW) of a medicine-teaching institution of North India. Materials and Methods: A cross-sectional study was conducted in the months of June and July 2021 among HCW (N = 203) of 18 years and above, vaccinated with at least the first dose of Covishield. A semi-structured, prevalidated, and pretested questionnaire was used to collect information through an interview schedule. The questionnaire was divided into five sections: the sociodemographic profile, behavioral characteristics, past medical history, COVID-19 awareness, and past infection and COVID-19 vaccine related information. Chi-squared test was applied to check the association of different factors with AEFI. Results: In our study, 73.89% of participants suffered from at least one AEFI after the first dose of the vaccine, while 48.66% had at least one AEFI after the second dose. Females reported significantly high AEFI for both doses (P = 0.001, 0.000). We found a significant association between the occurrence of AEFI and occupation (first dose P = 0.015), substance abuse (first dose P = 0.002), diet (first dose P = 0.016), and allergy (first dose P = 0.027). Other significant findings were headaches among HCW ≥40 years of age (dose P = 0.034) and systemic AEFI in participants with comorbidity (first dose P = 0.020). Conclusion: More AEFI were reported after the first dose as compared to the second dose. AEFI were more among females after both the doses. Occupation, substance use, diet, and history of allergy were significantly associated with AEFI.

Need for Objective Task-Based Evaluation of Image Segmentation Algorithms for Quantitative PET: A Study with ACRIN 6668/RTOG 0235 Multicenter Clinical Trial Data.

Reliable performance of PET segmentation algorithms on clinically relevant tasks is required for their clinical translation. However, these algorithms are typically evaluated using figures of merit (FoMs) that are not explicitly designed to correlate with clinical task performance. Such FoMs include the Dice similarity coefficient (DSC), the Jaccard similarity coefficient (JSC), and the Hausdorff distance (HD). The objective of this study was to investigate whether evaluating PET segmentation algorithms using these task-agnostic FoMs yields interpretations consistent with evaluation on clinically relevant quantitative tasks. Methods: We conducted a retrospective study to assess the concordance in the evaluation of segmentation algorithms using the DSC, JSC, and HD and on the tasks of estimating the metabolic tumor volume (MTV) and total lesion glycolysis (TLG) of primary tumors from PET images of patients with non-small cell lung cancer. The PET images were collected from the American College of Radiology Imaging Network 6668/Radiation Therapy Oncology Group 0235 multicenter clinical trial data. The study was conducted in 2 contexts: (1) evaluating conventional segmentation algorithms, namely those based on thresholding (SUVmax40% and SUVmax50%), boundary detection (Snakes), and stochastic modeling (Markov random field-Gaussian mixture model); (2) evaluating the impact of network depth and loss function on the performance of a state-of-the-art U-net-based segmentation algorithm. Results: Evaluation of conventional segmentation algorithms based on the DSC, JSC, and HD showed that SUVmax40% significantly outperformed SUVmax50%. However, SUVmax40% yielded lower accuracy on the tasks of estimating MTV and TLG, with a 51% and 54% increase, respectively, in the ensemble normalized bias. Similarly, the Markov random field-Gaussian mixture model significantly outperformed Snakes on the basis of the task-agnostic FoMs but yielded a 24% increased bias in estimated MTV. For the U-net-based algorithm, our evaluation showed that although the network depth did not significantly alter the DSC, JSC, and HD values, a deeper network yielded substantially higher accuracy in the estimated MTV and TLG, with a decreased bias of 91% and 87%, respectively. Additionally, whereas there was no significant difference in the DSC, JSC, and HD values for different loss functions, up to a 73% and 58% difference in the bias of the estimated MTV and TLG, respectively, existed. Conclusion: Evaluation of PET segmentation algorithms using task-agnostic FoMs could yield findings discordant with evaluation on clinically relevant quantitative tasks. This study emphasizes the need for objective task-based evaluation of image segmentation algorithms for quantitative PET.

A quality assurance framework for routine monitoring of deep learning cardiac substructure computed tomography segmentation models in radiotherapy.

BACKGROUND: For autosegmentation models, the data used to train the model (e.g., public datasets and/or vendor-collected data) and the data on which the model is deployed in the clinic are typically not the same, potentially impacting the performance of these models by a process called domain shift. Tools to routinely monitor and predict segmentation performance are needed for quality assurance. Here, we develop an approach to perform such monitoring and performance prediction for cardiac substructure segmentation. PURPOSE: To develop a quality assurance (QA) framework for routine or continuous monitoring of domain shift and the performance of cardiac substructure autosegmentation algorithms. METHODS: A benchmark dataset consisting of computed tomography (CT) images along with manual cardiac substructure delineations of 241 breast cancer radiotherapy patients were collected, including one "normal" image domain of clean images and five "abnormal" domains containing images with artifact (metal, contrast), pathology, or quality variations due to scanner protocol differences (field of view, noise, reconstruction kernel, and slice thickness). The QA framework consisted of an image domain shift detector which operated on the input CT images and a shape quality detector on the output of an autosegmentation model, and a regression model for predicting autosegmentation model performance. The image domain shift detector was composed of a trained denoising autoencoder (DAE) and two hand-engineered image quality features to detect normal versus abnormal domains in the input CT images. The shape quality detector was a variational autoencoder (VAE) trained to estimate the shape quality of the auto-segmentation results. The output from the image domain shift and shape quality detectors was used to train a regression model to predict the per-patient segmentation accuracy, measured by Dice coefficient similarity (DSC) to physician contours. Different regression techniques were investigated including linear regression, Bagging, Gaussian process regression, random forest, and gradient boost regression. Of the 241 patients, 60 were used to train the autosegmentation models, 120 for training the QA framework, and the remaining 61 for testing the QA framework. A total of 19 autosegmentation models were used to evaluate QA framework performance, including 18 convolutional neural network (CNN)-based and one transformer-based model. RESULTS: When tested on the benchmark dataset, all abnormal domains resulted in a significant DSC decrease relative to the normal domain for CNN models ( p < 0.001 $p < 0.001$ ), but only for some domains for the transformer model. No significant relationship was found between the performance of an autosegmentation model and scanner protocol parameters ( p = 0.42 $p = 0.42$ ) except noise ( p = 0.01 $p = 0.01$ ). CNN-based autosegmentation models demonstrated a decreased DSC ranging from 0.07 to 0.41 with added noise, while the transformer-based model was not significantly affected (ANOVA, p = 0.99 $p=0.99$ ). For the QA framework, linear regression models with bootstrap aggregation resulted in the highest mean absolute error (MAE) of 0.041 ± 0.002 $0.041 \pm 0.002$ , in predicted DSC (relative to true DSC between autosegmentation and physician). MAE was lowest when combining both input (image) detectors and output (shape) detectors compared to output detectors alone. CONCLUSIONS: A QA framework was able to predict cardiac substructure autosegmentation model performance for clinically anticipated "abnormal" domain shifts.

Beyond Average: α-Particle Distribution and Dose Heterogeneity in Bone Metastatic Prostate Cancer.

α-particle emitters are emerging as a potent modality for disseminated cancer therapy because of their high linear energy transfer and localized absorbed dose profile. Despite great interest and pharmaceutical development, there is scant information on the distribution of these agents at the scale of the α-particle pathlength. We sought to determine the distribution of clinically approved [223Ra]RaCl2 in bone metastatic castration-resistant prostate cancer at this resolution, for the first time to our knowledge, to inform activity distribution and dose at the near-cell scale. Methods: Biopsy specimens and blood were collected from 7 patients 24 h after administration. 223Ra activity in each sample was recorded, and the microstructure of biopsy specimens was analyzed by micro-CT. Quantitative autoradiography and histopathology were segmented and registered with an automated procedure. Activity distributions by tissue compartment and dosimetry calculations based on the MIRD formalism were performed. Results: We revealed the activity distribution differences across and within patient samples at the macro- and microscopic scales. Microdistribution analysis confirmed localized high-activity regions in a background of low-activity tissue. We evaluated heterogeneous α-particle emission distribution concentrated at bone-tissue interfaces and calculated spatially nonuniform absorbed-dose profiles. Conclusion: Primary patient data of radiopharmaceutical therapy distribution at the small scale revealed that 223Ra uptake is nonuniform. Dose estimates present both opportunities and challenges to enhance patient outcomes and are a first step toward personalized treatment approaches and improved understanding of α-particle radiopharmaceutical therapies.

Simulating Silicon Photomultiplier Response to Scintillation Light.

The response of a Silicon Photomultiplier (SiPM) to optical signals is affected by many factors including photon-detection efficiency, recovery time, gain, optical crosstalk, afterpulsing, dark count, and detector dead time. Many of these parameters vary with overvoltage and temperature. When used to detect scintillation light, there is a complicated non-linear relationship between the incident light and the response of the SiPM. In this paper, we propose a combined discrete-time discrete-event Monte Carlo (MC) model to simulate SiPM response to scintillation light pulses. Our MC model accounts for all relevant aspects of the SiPM response, some of which were not accounted for in the previous models. We also derive and validate analytic expressions for the single-photoelectron response of the SiPM and the voltage drop across the quenching resistance in the SiPM microcell. These analytic expressions consider the effect of all the circuit elements in the SiPM and accurately simulate the time-variation in overvoltage across the microcells of the SiPM. Consequently, our MC model is able to incorporate the variation of the different SiPM parameters with varying overvoltage. The MC model is compared with measurements on SiPM-based scintillation detectors and with some cases for which the response is known a priori. The model is also used to study the variation in SiPM behavior with SiPM-circuit parameter variations and to predict the response of a SiPM-based detector to various scintillators.

Need for objective task-based evaluation of AI-based segmentation methods for quantitative PET.

Artificial intelligence (AI)-based methods are showing substantial promise in segmenting oncologic positron emission tomography (PET) images. For clinical translation of these methods, assessing their performance on clinically relevant tasks is important. However, these methods are typically evaluated using metrics that may not correlate with the task performance. One such widely used metric is the Dice score, a figure of merit that measures the spatial overlap between the estimated segmentation and a reference standard (e.g., manual segmentation). In this work, we investigated whether evaluating AI-based segmentation methods using Dice scores yields a similar interpretation as evaluation on the clinical tasks of quantifying metabolic tumor volume (MTV) and total lesion glycolysis (TLG) of primary tumor from PET images of patients with non-small cell lung cancer. The investigation was conducted via a retrospective analysis with the ECOG-ACRIN 6668/RTOG 0235 multi-center clinical trial data. Specifically, we evaluated different structures of a commonly used AI-based segmentation method using both Dice scores and the accuracy in quantifying MTV/TLG. Our results show that evaluation using Dice scores can lead to findings that are inconsistent with evaluation using the task-based figure of merit. Thus, our study motivates the need for objective task-based evaluation of AI-based segmentation methods for quantitative PET.