Incorporating Patient-Specific Variability in the Simulation of Realistic Whole-Body $^{18}{\hbox{F-FDG}}$ Distributions for Oncology Applications

2009 ◽  
Vol 97 (12) ◽  
pp. 2026-2038 ◽  
Author(s):  
Amandine Le Maitre ◽  
William Paul Segars ◽  
Simon Marache ◽  
Anthonin Reilhac ◽  
Mathieu Hatt ◽  
...  
Keyword(s):  
2021 ◽  
Vol 8 (1) ◽  
Author(s):  
Jingjie Shang ◽  
Zhiqiang Tan ◽  
Yong Cheng ◽  
Yongjin Tang ◽  
Bin Guo ◽  
...  

Abstract Background Standardized uptake value (SUV) normalized by lean body mass ([LBM] SUL) is recommended as metric by PERCIST 1.0. The James predictive equation (PE) is a frequently used formula for LBM estimation, but may cause substantial error for an individual. The purpose of this study was to introduce a novel and reliable method for estimating LBM by limited-coverage (LC) CT images from PET/CT examinations and test its validity, then to analyse whether SUV normalised by LC-based LBM could change the PERCIST 1.0 response classifications, based on LBM estimated by the James PE. Methods First, 199 patients who received whole-body PET/CT examinations were retrospectively retrieved. A patient-specific LBM equation was developed based on the relationship between LC fat volumes (FVLC) and whole-body fat mass (FMWB). This equation was cross-validated with an independent sample of 97 patients who also received whole-body PET/CT examinations. Its results were compared with the measurement of LBM from whole-body CT (reference standard) and the results of the James PE. Then, 241 patients with solid tumours who underwent PET/CT examinations before and after treatment were retrospectively retrieved. The treatment responses were evaluated according to the PE-based and LC-based PERCIST 1.0. Concordance between them was assessed using Cohen’s κ coefficient and Wilcoxon’s signed-ranks test. The impact of differing LBM algorithms on PERCIST 1.0 classification was evaluated. Results The FVLC were significantly correlated with the FMWB (r=0.977). Furthermore, the results of LBM measurement evaluated with LC images were much closer to the reference standard than those obtained by the James PE. The PE-based and LC-based PERCIST 1.0 classifications were discordant in 27 patients (11.2%; κ = 0.823, P=0.837). These discordant patients’ percentage changes of peak SUL (SULpeak) were all in the interval above or below 10% from the threshold (±30%), accounting for 43.5% (27/62) of total patients in this region. The degree of variability is related to changes in LBM before and after treatment. Conclusions LBM algorithm-dependent variability in PERCIST 1.0 classification is a notable issue. SUV normalised by LC-based LBM could change PERCIST 1.0 response classifications based on LBM estimated by the James PE, especially for patients with a percentage variation of SULpeak close to the threshold.


Author(s):  
D. Keith Walters ◽  
Greg W. Burgreen ◽  
Robert L. Hester ◽  
David S. Thompson ◽  
David M. Lavallee ◽  
...  

Computational fluid dynamics (CFD) simulations were performed for unsteady periodic breathing conditions, using large-scale models of the human lung airway. The computational domain included fully coupled representations of the orotracheal region and large conducting zone up to generation four (G4) obtained from patient-specific CT data, and the small conducting zone (to G16) obtained from a stochastically generated airway tree with statistically realistic geometrical characteristics. A reduced-order geometry was used, in which several airway branches in each generation were truncated, and only select flow paths were retained to G16. The inlet and outlet flow boundaries corresponded to the oronasal opening (superior), the inlet/outlet planes in terminal bronchioles (distal), and the unresolved airway boundaries arising from the truncation procedure (intermediate). The cyclic flow was specified according to the predicted ventilation patterns for a healthy adult male at three different activity levels, supplied by the whole-body modeling software HumMod. The CFD simulations were performed using Ansys FLUENT. The mass flow distribution at the distal boundaries was prescribed using a previously documented methodology, in which the percentage of the total flow for each boundary was first determined from a steady-state simulation with an applied flow rate equal to the average during the inhalation phase of the breathing cycle. The distal pressure boundary conditions for the steady-state simulation were set using a stochastic coupling procedure to ensure physiologically realistic flow conditions. The results show that: 1) physiologically realistic flow is obtained in the model, in terms of cyclic mass conservation and approximately uniform pressure distribution in the distal airways; 2) the predicted alveolar pressure is in good agreement with previously documented values; and 3) the use of reduced-order geometry modeling allows accurate and efficient simulation of large-scale breathing lung flow, provided care is taken to use a physiologically realistic geometry and to properly address the unsteady boundary conditions.


Author(s):  
Mao Li ◽  
Karol Miller ◽  
Grand Joldes ◽  
Ron Kikinis ◽  
Adam Wittek

2014 ◽  
Vol 136 (10) ◽  
Author(s):  
D. Keith Walters ◽  
Greg W. Burgreen ◽  
Robert L. Hester ◽  
David S. Thompson ◽  
David M. Lavallee ◽  
...  

Computational fluid dynamics (CFD) simulations were performed using large-scale models of the human lung airway and unsteady periodic breathing conditions. The computational domain included fully coupled representations of the orotracheal region and large conducting zone up to generation four (G4) obtained from patient-specific CT data, and the small conducting zone (to the 16th generation) obtained from a stochastically generated airway tree with statistically realistic morphological characteristics. A reduced-geometry airway model was used, in which several airway branches in each generation were truncated, and only select flow paths were retained to the 16th generation. The inlet and outlet flow boundaries corresponded to the oral opening, the physical inlet/outlet boundaries at the terminal bronchioles, and the unresolved airway boundaries created from the truncation procedure. The total flow rate was specified according to the expected ventilation pattern for a healthy adult male, which was supplied by the whole-body modeling software HumMod. The unsteady mass flow distribution at the distal boundaries was prescribed based on a preliminary steady-state simulation with an applied flow rate equal to the average flow rate during the inhalation phase of the breathing cycle. In contrast to existing studies, this approach allows fully coupled simulation of the entire conducting zone, with no need to specify distal mass flow or pressure boundary conditions a priori, and without the use of impedance or one-dimensional (1D) flow models downstream of the truncated boundaries. The results show that: (1) physiologically realistic flow is obtained in the model, in terms of cyclic mass conservation and approximately uniform pressure distribution in the distal airways; (2) the predicted alveolar pressure is in good agreement with correlated experimental data; and (3) the use of reduced-order geometry modeling allows accurate and efficient simulation of large-scale breathing lung flow, provided care is taken to use a physiologically realistic geometry and to properly address the unsteady boundary conditions.


2020 ◽  
Vol 59 (05) ◽  
pp. 365-374
Author(s):  
Theresa Ida Götz ◽  
Elmar Wolfgang Lang ◽  
Olaf Prante ◽  
Michael Cordes ◽  
Torsten Kuwert ◽  
...  

Abstract Objective Patients with advanced prostate cancer are suitable candidates for [177Lu]PSMA-617 therapy. Integrated SPECT/CT systems have the potential to improve the accuracy of patient-specific tumor dosimetry. We present a novel patient-specific Monte Carlo based voxel-wise dosimetry approach to determine organ and total tumor doses (TTD). Methods 13 patients with histologically confirmed metastasized castration-resistant prostate cancer were treated with a total of 18 cycles of [177Lu]PSMA-617 therapy. In each patient, dosimetry was performed after the first cycle of [177Lu]PSMA-617 therapy. Regions of interest were defined manually on the SPECT/CT images for the kidneys, spleen and all 295 PSMA-positive tumor lesions in the field of view. The absorbed dose to normal organs and to all tumor lesions were calculated by a three dimensional dosimetry method based on Monte Carlo Simulations. Results The average dose values yielded the following results: 2.59 ± 0.63 Gy (1.67–3.92 Gy) for the kidneys, 0.79 ± 0.46 Gy (0.31–1.90 Gy) for the spleen and 11.00 ± 11.97 Gy (1.28–49.10 Gy) for all tracer-positive tumor lesions. A trend towards higher TTD was observed in patients with Gleason Scores > 8 compared to Gleason Scores ≤ 8 and in lymph node metastases compared to bone metastases. A significant correlation was determined between the serum-PSA level before RLT and the TTD (r = –0.57, p < 0.05), as well as between the TTD with the percentage change of serum-PSA levels before and after therapy was observed (r = –0.57, p < 0.05). Patients with higher total tumor volumes of PSMA-positive lesions demonstrated significantly lower kidney average dose values (r = –0.58, p < 0.05). Conclusion The presented novel Monte Carlo based voxel-wise dosimetry calculates a patient specific whole-body dose distribution, thus taking into account individual anatomies and tissue compositions showing promising results for the estimation of radiation doses of normal organs and PSMA-positive tumor lesions.


Author(s):  
Wenli Liu ◽  
Hongkai Wang ◽  
Pu Zhang ◽  
Chengwei Li ◽  
Jie Sun ◽  
...  

The accurate estimation of patient’s exposure to the radiofrequency (RF) electromagnetic field of magnetic resonance imaging (MRI) significantly depends on a precise individual anatomical model. In the study, we investigated the applicability of an efficient whole-body individual modelling method for the assessment of MRI RF exposure. The individual modelling method included a deformable human model and tissue simplification techniques. Besides its remarkable efficiency, this approach utilized only a low specific absorption rate (SAR) sequence or even no MRI scan to generate the whole-body individual model. Therefore, it substantially reduced the risk of RF exposure. The dosimetric difference of the individual modelling method was evaluated using the manually segmented human models. In addition, stochastic dosimetry using a surrogate model by polynomial chaos presented SAR variability due to body misalignment and tilt in the coil, which were frequently occurred in the practical scan. In conclusion, the dosimetric equivalence of the individual models was validated by both deterministic and stochastic dosimetry. The proposed individual modelling method allowed the physicians to quantify the patient-specific SAR while the statistical results enabled them to comprehensively weigh over the exposure risk and get the benefit of imaging enhancement by using the high-intensity scanners or the high-SAR sequences.


2020 ◽  
Author(s):  
Srinivas Chilukuri ◽  
Sham Sundar ◽  
Rajesh Thiyagarajan ◽  
Jose Easow ◽  
Mayur Sawant ◽  
...  

Abstract Objective To standardize the technique and resources for total marrow and lymphoid irradiation (TMLI) as part of the conditioning regimen before allogenic bone marrow transplantation (ABMT) using helical tomotherapy.Methods We used this technique in our first 5 patients requiring TMLI. Patients were immobilized using a mask and a whole-body vacuum cushion. CT scanning was performed in head first supine (HFS) and feet first supine (FFS) orientations with an overlap at mid-thigh. Target consisted of the entire skeleton, spleen, sanctuary sites and major lymphatics whereas lungs, kidneys, aero-digestive tract, bowel, parotids, heart and liver were defined as organs at risk (OAR). Treatment was performed in two parts based on 2 different plans generated in HFS and FFS orientations with an overlap at the mid thigh. Patients along with the immobilization device were manually rotated by 180° to change the orientation after the delivery of HFS plan. The dose at the junction was contributed by a complementary dose gradient from each of the plans. Plan was to deliver 95% of 12Gy to 98% of CTV with dose heterogeneity < 10% and pre-specified OAR doe constraints. Megavoltage-CT was used for position verification before each fraction. Patient specific quality assurance and an in-vivo film dosimetry to verify junction dose were performed in all patients.Results Treatment was delivered in two daily fractions of 2Gy each for 3 days with at least 8-hours gap between each fraction. The target coverage goals were met in all the patients. The average person-hours per patient were 16.5, 21.5 and 25.75 for radiation oncologist, radiation therapist and medical physicist respectively. Average in-room time per patient was 9.25 hours with an average beam-on time of 3.32 hours for all the six fractions. Conclusion This report comprehensively describes technique and resource requirements for TMLI and would serve as a practical guide for departments keen to start this service. Despite being time and labor intensive, it can be implemented safely and robustly. We will be using this methodology in a prospective phase II trial to study safety and feasibility of dose escalated TMLI as part of conditioning regimen before ABMT.


Author(s):  
Werner Backfrieder

"Patient specific dosimetry established during the last decade in modern radio-therapy. Usually, tracer kinetics in main compartments of observed metabolism is assessed from anterior and posterior whole body scans. The effective doses for each organ, derived by the MIRD scheme, provide evidence for following radiotherapeutic treatment and helps to meet vital dose limits for critical organs, e.g. kidneys. The calculation of individual dose in a three-dimensional context leads to more accurate dose estimates, as was proven by intensive research, but is still on the cusp to clinical application. In this work, a statistical approach, based on multi-modal image and feature data, is presented, to overcome manual segmentation, the most time consuming step, in 3D based dose calculation. 3D data volumes from a hybrid SPECT study, comprising SPECT and CT data, covering main compartments of metabolism, build the image features of a Gaussian classifier. From prior segmentations organspecific membership maps are derived, and substituted as additional feature into the segmentation procedure. Centroids, eccentricity and principal axes of organ models are registered to a rough thresholded image of the SPECT component, and define membership coefficients of the voxels. The new approach yields accurate results, even with real patient data. The new method needs minimal user interaction during selection of some sample regions, thus showing high potential for implementation in a clinical workflow."


Blood ◽  
2016 ◽  
Vol 128 (22) ◽  
pp. 662-662
Author(s):  
Issa F. Khouri ◽  
Bill Erwin ◽  
Alison M Gulbis ◽  
Francesco Turturro ◽  
S Cheenu Kappadath ◽  
...  

Abstract Background: Nonmyeloablative allogeneic transplantation has the potential to induce long-term remissions in patients with relapsed lymphoma. However, a non-intense conditioning regimen enhances the risk of early relapse. Anti-CD20 antibody radioimmunotherapy (90YIT) delivers radiation dose not only to the tumor cells that bind the antibody but also to inaccessible neighboring cells as a result of the cross-fire effect. Thus, we hypothesized that the addition of escalated 90YIT dose to the recently published bendamustine+fludarabine conditioning regimen (Khouri et al. Blood 2014) would facilitate early cytoreduction in such patients and promote improved long-term disease control by the allogeneic graft. Organ doses from a 90YIT weight-based activity prescription (mCi/kg) vary considerably, which justifies a dosimetry-based strategy for mCi/kg escalation. Methods and patients: On days -22 and -14, rituximab was given at 250 mg/m2 preceding 111In ibritumumab and 90YIT administration, respectively. Organ dosimetric assessment was performed based on serial 111In ibritumumab whole body scanning (0, 4, 24, 72 and 144 hours) , to select from among five 90YIT mCi/kg prescriptions (0.5, 0.75, 1, 1.25 or 1.5) that would result in an estimated 10 - 12 Gy dose to the liver, lungs or kidneys. Organ dose was corrected for patient-specific mass, based on a CT volume estimate times 1.03 g/cc for liver and kidneys, and a variable specific gravity for lungs (Simon, J Clin Monit Comput, 2000). Bendamustine 130 mg/m2 plus 30 mg/m2 of fludarabine IV were given daily on days -5 to - 3 prior to transplantation. Tacrolimus and mini-methotrexate (Mycophenolate mofetil in case of cord blood transplantation) were used for GVHD prophylaxis. In addition, thymoglobulin 1 mg/kg IV was given on days -2, and -1 in patients receiving an unrelated donor transplant. Results: Twenty patients were studied. The median age was 58 years (range, 37-71). Lymphoma histologies included: indolent (n=8, 40%), diffuse large cell (n=6; 30%), double-hit (n=2; 10%) and mantle cell (n= 4, 20%). The median number of prior chemotherapies received was 4 (range, 2-7). At study entry, 8 patients (40%) were in complete remission following salvage therapy, 7 (35%) were in partial response, and 5 (25%) had refractory disease. Six of 16 (37.5%) patients tested were PET+. Dosimetry: The most exposed organ was either liver (16 patients) or lungs (4 patients). The distribution among the five 90YIT mCi/kg prescriptions (smallest to largest) was 2, 4, 12, 1 and 1, with a mean of 0.94 ± 0.23 mCi/kg. If all twenty patients were treated at 1 mCi/kg (the most common prescription), the 20 Gy limit employed for 90YIT clinical trials prior to approval would have been exceeded in only one patient for the liver (22.9 Gy) or lungs (20.9 Gy). The maximum liver and lung doses at 0.75 mCi/kg would have been 17.2 and 15.7 Gy, respectively. Transplant outcomes: Fifteen patients (75%) received their transplants from unrelated donors (including 1 mismatched and 2 cord blood), and only 5 (25%) from HLA-compatible siblings. The median number of CD34+ cells infused was 6.2 × 106/kg. Neutrophil counts recovered to > 0.5 × 109/L after a median of 12 days (range, 0-24 days). Platelet counts recovered to > 20 × 109/L after a median of 19 days (range, 9-30 days). By day 30, median donor myeloid and T-cells were 100% (range, 98-100). The cumulative incidence of acute grade 2-4 GVHD and chronic extensive GVHD were 25% (5% for acute grade 3-4) and 32%, respectively. Treatment-related mortality (TRM) rates at day 100 and 1 year after transplantation were 0% and 10%, respectively. The 2 cord blood transplants engrafted with 100% donor cells and none had GVHD. With a median follow-up duration of 14 months (range, 3-34 months), the overall survival and progression-free survival rates were 85% and 70%, respectively. No significant difference in survival or TRM could be detected by age, donor type, histology, disease status, PET status or number of prior therapies. Conclusions: Our results indicate that dose-intense 90YIT combined with fludarabine and bendamustine is a well-tolerated nonmyeloablative allogeneic conditioning for lymphoid malignancies, with promising results of engraftment, GVHD and survival. Our stratified 90YIT prescription results suggest that future studies with a fixed dose of 1 mCi/kg level without dosimetry would have an acceptable radiation risk to vital organs in this setting. Disclosures Jabbour: ARIAD: Consultancy, Research Funding; Pfizer: Consultancy, Research Funding; Novartis: Research Funding; BMS: Consultancy.


Sign in / Sign up

Export Citation Format

Share Document