Comparison between model-predicted tumor oxygenation dynamics and vascular-/flow-related Doppler indices

PURPOSE

Mathematical modeling is a powerful and flexible method to investigate complex phenomena. It discloses the possibility of reproducing expensive as well as invasive experiments in a safe environment with limited costs. This makes it suitable to mimic tumor evolution and response to radiotherapy although the reliability of the results remains an issue. Complexity reduction is therefore a critical aspect in order to be able to compare model outcomes to clinical data. Among the factors affecting treatment efficacy, tumor oxygenation is known to play a key role in radiotherapy response. In this work, we aim at relating the oxygenation dynamics, predicted by a macroscale model trained on tumor volumetric data of uterine cervical cancer patients, to vascularization and blood flux indices assessed on Ultrasound Doppler images.

METHODS

We propose a macroscale model of tumor evolution based on three dynamics, namely active portion, necrotic portion, and oxygenation. The model parameters were assessed on the volume size of seven cervical cancer patients administered with 28 fractions of intensity modulated radiation therapy (IMRT) (1.8 Gy/fraction). For each patient, five Doppler ultrasound tests were acquired before, during, and after the treatment. The lesion was manually contoured by an expert physician using 4D View^® (General Electric Company - Fairfield, Connecticut, United States), which automatically provided the overall tumor volume size along with three vascularization and/or blood flow indices. Volume data only were fed to the model for training purpose, while the predicted oxygenation was compared a posteriori to the measured Doppler indices.

RESULTS

The model was able to fit the tumor volume evolution within 8% error (range: 3-8%). A strong correlation between the intrapatient longitudinal indices from Doppler measurements and oxygen predicted by the model (about 90% or above) was found in three cases. Two patients showed an average correlation value (50-70%) and the remaining two presented poor correlations. The latter patients were the ones featuring the smallest tumor reduction throughout the treatment, typical of hypoxic conditions. Moreover, the average oxygenation value predicted by the model was close to the average vascularization-flow index (average difference: 7%).

CONCLUSIONS

The results suggest that the modeled relation between tumor evolution and oxygen dynamics was reasonable enough to provide realistic oxygenation curves in five cases (correlation greater than 50%) out of seven. In case of nonresponsive tumors, the model failed in predicting the oxygenation trend while succeeded in reproducing the average oxygenation value according to the mean vascularization-flow index. Despite the need for deeper investigations, the outcomes of the present work support the hypothesis that a simple macroscale model of tumor response to radiotherapy is able to predict the tumor oxygenation. The possibility of an objective and quantitative validation on imaging data discloses the possibility to translate them as decision support tools in clinical practice and to move a step forward in the treatment personalization.

2D/3D reconstruction of the distal femur using statistical shape models addressing personalized surgical instruments in knee arthroplasty: A feasibility analysis

BACKGROUND

Personalized surgical instruments (PSI) have gained success in the domain of total knee replacement, demonstrating clinical outcomes similar or even superior to both traditional and navigated surgeries. The key requirement for prototyping PSI is the availability of the digital bony surface. In this paper, we aim at verifying whether the 2D/3D reconstruction of the distal femur, based on statistical shape models (SSM), grants sufficient accuracy, especially in the condylar regions, to support a PSI technique.

METHODS

Computed tomographic knee datasets acquired on 100 patients with severe cartilage damage were retrospectively considered in this work. All the patients underwent total knee replacement using the PSI-based surgical technique. Eighty out of 100 reconstructed distal femur surfaces were used to build the statistical model. The remaining 20 surfaces were used for testing. The 2D/3D reconstruction process was based on digital reconstructed radiographies (DRRs) obtained with a simulated X-ray projection process. An iterative optimization procedure, based on an evolutionary algorithm, systematically morphed the statistical model to decrease the difference between the DRR, obtained by the original CT dataset, and the DRR obtained from the morphed surface.

RESULTS

Over the 80 variations, the first ten modes were found sufficient to reconstruct the distal femur surface with accuracy. Using three DRR, the maximum Hausdorff and RMS distance errors were lower than 1.50 and 0.75 mm, respectively. As expected, the reconstruction quality improved by increasing the number of DRRs. Statistical difference (P < 0.001) was found in the 2 vs 3, 2 vs 4 and 2 vs 5 DRR, thus proving that adding just a single displaced projection to the two traditional sagittal and coronal X-ray images improved significantly the reconstruction quality. The effect of the PSI contact area errors on the distal cut direction featured a maximum median error lower than 2° and 0.5° on the sagittal and frontal plane, respectively. Statistical difference was found (P < 0.0001) in the reconstruction accuracy when comparing SSM built using pathologic with respect to non-pathologic shapes (cadavers), meaning that, to improve the patient-specific reconstruction, the morphologic anomalies, specific to the pathology, must be embedded into the SSM.

CONCLUSIONS

We showed that the X-ray based reconstruction of the distal femur is reasonable also in presence of pathologic bony conditions, featuring accuracy results similar to earlier reports in the literature that reconstructed normal femurs. This finding discloses the chance of applying the proposed methodology to the reconstruction of bony surfaces used in the PSI surgical approach.

The role of platelets in antiphospholipid syndrome

Despite extensive research, the pathogenesis of antiphospholipid syndrome (APS) remains obscure in many aspects. However, it is widely accepted that thrombosis is the result of a hypercoagulable state caused by antibodies directed against β2-glycoprotein I (β2-GPI), a protein whose physiological role is unknown. Although underestimated, platelets may be involved in APS and its thrombotic manifestations, especially arterial, in several ways. Thrombocytopenia is the most relevant non-criteria manifestation of APS, possibly caused by direct binding of anti-β2-GPI antibodies or anti-β2-GPI-β2-GPI complexes. On the other hand, platelets may have a key role in APS-related thrombosis due to the presence of multiple receptors that can interact with anti-β2-GPI antibodies (especially apolipoprotein E receptor 2' (apoER2') and glycoprotein Ibα (GPIbα)) with consequent release of different procoagulant mediators such as thromboxane B2, platelet factor 4 (PF4), and platelet factor 4 variant (CXCL4L1). The aim of this review is to put together evidence on the possible role of platelets in APS and to stimulate further research on the issue.

Patient-specific modeling of the trochlear morphologic anomalies by means of hyperbolic paraboloids

Diagnostic and therapeutic purposes are issuing pressing demands to improve the evaluation of the dysplasia condition of the femoral trochlea. The traditional clinical assessment of the dysplasia, based on Dejour classification, recognized 4 increasing (A, B, C, D) levels of severity. It has been extensively questioned in the literature that this classification methodology can be defective suggesting that quantitative measures can ensure more reliable criteria for the dysplasia severity assessment. This study reports on a novel technique to model the trochlear surface (TS), digitally reconstructed by 3D volumetric imaging, using three hyperbolic paraboloids (HP), one to describe the global trochlear aspect, two to represent the local aspects of the medial and lateral compartments, respectively. Results on a cohort of 43 patients, affected by aspecific anterior knee pain, demonstrate the consistency of the estimated model parameters with the morphologic aspect of the TS. The obtained small fitting error (on average lower than 0.80 mm) demonstrated that the ventral aspect of the trochlear morphology can be modeled with high accuracy by HPs. We also showed that HP modeling provides a continuous representation of morphologic variations in shape parameter space while we found that similar morphologic anomalies of the trochlear aspect are actually attributed to different severity grades in the Dejour classification. This finding is in agreement with recent works in the literature reporting that morphometric parameters can only optimistically be used to discriminate between the Grade A and the remaining three grades. In conclusion, we can assert that the proposed methodology is a further step toward modeling of anatomical surfaces that can be used to quantify deviations to normality on a patient-specific basis.

Mathematical modeling of tumor response to radiation: radio-sensitivity correlation with BOLD, TOLD, ΔR1 and ΔR2* investigated in large Dunning R3327-AT1 rat prostate tumors

Tumor response to radiation therapy can vary highly across patients. Several factors, both tumor- and environment-specific, can influence the radio-sensitivity, one of the most well-known being hypoxia. In this work, we investigated possible correlations between the radio-sensitivity parameters determined by means of a simple mathematical model of tumor volume evolution, and the MRI-based indicators of oxygenation in Dunning R3327-AT1 rats. Prior to irradiation the rats were subjected to an oxygen-breathing challenge, which was evaluated by MRI. The tumors were administered a single irradiation dose (30 Gy), while breathing air or oxygen. Despite a poor fitting performance, the model was able to identify two different tumor volume regression patterns. Moreover, the radio-sensitivity of the oxygen-breathing group was found to correlate with the variation of the transverse relaxation rate ΔR2* (-0.89). This suggests that MRI-based indices of tumor oxygenation may provide information about radio-sensitivity.

MRI quantification of pancreas motion as a function of patient setup for particle therapy -a preliminary study

Particle therapy (PT) has shown positive therapeutic results in local control of locally advanced pancreatic lesions. PT effectiveness is highly influenced by target localization accuracy both in space, since the pancreas is located in proximity to radiosensitive vital organs, and in time as it is subject to substantial breathing-related motion. The purpose of this preliminary study was to quantify pancreas range of motion under typical PT treatment conditions. Three common immobilization devices (vacuum cushion, thermoplastic mask, and compressor belt) were evaluated on five male patients in prone and supine positions. Retrospective four-dimensional magnetic resonance imaging data were reconstructed for each condition and the pancreas was manually segmented on each of six breathing phases. A k-means algorithm was then applied on the manually segmented map in order to obtain clusters representative of the three pancreas segments: head, body, and tail. Centers of mass (COM) for the pancreas and its segments were computed, as well as their displacements with respect to a reference breathing phase (beginning exhalation). The median three-dimensional COM displacements were in the range of 3 mm. Latero-lateral and superior-inferior directions had a higher range of motion than the anterior-posterior direction. Motion analysis of the pancreas segments showed slightly lower COM displacements for the head cluster compared to the tail cluster, especially in prone position. Statistically significant differences were found within patients among the investigated setups. Hence a patient-specific approach, rather than a general strategy, is suggested to define the optimal treatment setup in the frame of a millimeter positioning accuracy.

3D-printed applicators for high dose rate brachytherapy: Dosimetric assessment at different infill percentage

PURPOSE

Dosimetric assessment of high dose rate (HDR) brachytherapy applicators, printed in 3D with acrylonitrile butadiene styrene (ABS) at different infill percentage.

MATERIALS AND METHODS

A low-cost, desktop, 3D printer (Hamlet 3DX100, Hamlet, Dublin, IE) was used for manufacturing simple HDR applicators, reproducing typical geometries in brachytherapy: cylindrical (common in vaginal treatment) and flat configurations (generally used to treat superficial lesions). Printer accuracy was investigated through physical measurements. The dosimetric consequences of varying the applicator's density by tuning the printing infill percentage were analysed experimentally by measuring depth dose profiles and superficial dose distribution with Gafchromic EBT3 films (International Specialty Products, Wayne, NJ). Dose distributions were compared to those obtained with a commercial superficial applicator.

RESULTS

Measured printing accuracy was within 0.5mm. Dose attenuation was not sensitive to the density of the material. Surface dose distribution comparison of the 3D printed flat applicators with respect to the commercial superficial applicator showed an overall passing rate greater than 94% for gamma analysis with 3% dose difference criteria, 3mm distance-to-agreement criteria and 10% dose threshold.

CONCLUSION

Low-cost 3D printers are a promising solution for the customization of the HDR brachytherapy applicators. However, further assessment of 3D printing techniques and regulatory materials approval are required for clinical application.

Enhanced Surface Registration Techniques for Patient Positioning Control in Breast Cancer Radiotherapy

Experimental data, describing patient inter-fractional set-up errors in the clinical practice of breast radiotherapy, were exploited to simulate the performance in errors detection and correction of a constrained surface registration procedure, based on a hybrid configuration of control points (passive and laser). During 47 treatment sessions in three patients undergoing post-quadrantectomy radiotherapy, an opto-electronic localizer was used to acquire the three-dimensional coordinates of the hybrid control points, being two passive markers placed on selected skin landmarks on the sternum. Laser scanning technique was also applied for the acquisition of the 3-D surface model of the irradiated body area, which was used as reference for the automatic position correction procedure. A constrained surface registration algorithm was applied to estimate the rigid spatial transformation, describing the local errors affecting the control points. The improvement of the irradiation geometrical setup, by correcting the patient position according to the estimated spatial transformation parameters, was simulated. Results showed that the proposed surface registration method allowed us to detect and significantly (Wilcoxon signed rank analysis) reduce the initial misalignments, which exhibited overall median and 75th percentile values equal to 4.26 mm and 5.76 mm. Simulated residual errors dropped down to median and 75th percentile values measuring 2.95 mm and 3.87 mm, respectively. These results confirmed the high potentiality of surface registration techniques for the opto-electronic automatic patient positioning control in breast cancer radiotherapy.

Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

PURPOSE

Intensity modulated proton therapy (IMPT) of head and neck (H&N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation.

METHODS

Datasets of six H&N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT.

RESULTS

The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [-4, 2 Gy] for parotid glands (D(mean)), spinal cord (D(2%)), brainstem (D(2%)), and CTV (D(95%)). When using DIR generated contours for the vCT, those differences ranged within [-8, 11 Gy].

CONCLUSIONS

In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.