Uncovering prostate cancer aggressiveness signal in T2-weighted MRI through a three-reference tissues normalization technique.

Quantitative T2-weighted MRI (T2W) interpretation is impeded by the variability of acquisition-related features, such as field strength, coil type, signal amplification, and pulse sequence parameters. The main purpose of this work is to develop an automated method for prostate T2W intensity normalization. The procedure includes the following: (i) a deep learning-based network utilizing MASK R-CNN for automatic segmentation of three reference tissues: gluteus maximus muscle, femur, and bladder; (ii) fitting a spline function between average intensities in these structures and reference values; and (iii) using the function to transform all T2W intensities. The T2W distributions in the prostate cancer regions of interest (ROIs) and normal appearing prostate tissue (NAT) were compared before and after normalization using Student's t-test. The ROIs' T2W associations with the Gleason Score (GS), Decipher genomic score, and a three-tier prostate cancer risk were evaluated with Spearman's correlation coefficient (rS ). T2W differences in indolent and aggressive prostate cancer lesions were also assessed. The MASK R-CNN was trained with manual contours from 32 patients. The normalization procedure was applied to an independent MRI dataset from 83 patients. T2W differences between ROIs and NAT significantly increased after normalization. T2W intensities in 231 biopsy ROIs were significantly negatively correlated with GS (rS = -0.21, p = 0.001), Decipher (rS = -0.193, p = 0.003), and three-tier risk (rS = -0.235, p < 0.001). The average T2W intensities in the aggressive ROIs were significantly lower than in the indolent ROIs after normalization. In conclusion, the automated triple-reference tissue normalization method significantly improved the discrimination between prostate cancer and normal prostate tissue. In addition, the normalized T2W intensities of cancer exhibited a significant association with tumor aggressiveness. By improving the quantitative utilization of the T2W in the assessment of prostate cancer on MRI, the new normalization method represents an important advance over clinical protocols that do not include sequences for the measurement of T2 relaxation times.

The role of radiomics in prostate cancer radiotherapy.

"Radiomics," as it refers to the extraction and analysis of a large number of advanced quantitative radiological features from medical images using high-throughput methods, is perfectly suited as an engine for effectively sifting through the multiple series of prostate images from before, during, and after radiotherapy (RT). Multiparametric (mp)MRI, planning CT, and cone beam CT (CBCT) routinely acquired throughout RT and the radiomics pipeline are developed for extraction of thousands of variables. Radiomics data are in a format that is appropriate for building descriptive and predictive models relating image features to diagnostic, prognostic, or predictive information. Prediction of Gleason score, the histopathologic cancer grade, has been the mainstay of the radiomic efforts in prostate cancer. While Gleason score (GS) is still the best predictor of treatment outcome, there are other novel applications of quantitative imaging that are tailored to RT. In this review, we summarize the radiomics efforts and discuss several promising concepts such as delta-radiomics and radiogenomics for utilizing image features for assessment of the aggressiveness of prostate cancer and its outcome. We also discuss opportunities for quantitative imaging with the advance of instrumentation in MRI-guided therapies.

Increasing the efficiency of cone-beam CT based delta-radiomics using automated contours to predict radiotherapy-related toxicities in prostate cancer.

Extracting longitudinal image quantitative data, known as delta-radiomics, has the potential to capture changes in a patient's anatomy throughout the course of radiation treatment for prostate cancer. Some of the major challenges of delta-radiomics studies are contouring the structures for individual fractions and accruing patients' data in an efficient manner. The manual contouring process is often time consuming and would limit the efficiency of accruing larger sample sizes for future studies. The problem is amplified because the contours are often made by highly trained radiation oncologists with limited time to dedicate to research studies of this nature. This work compares the use of automated prostate contours generated using a deformable image-based algorithm to make predictive models of genitourinary and changes in total international prostate symptom score in comparison to manually contours for a cohort of fifty patients. Area under the curve of manual and automated models were compared using the Delong test. This study demonstrated that the delta-radiomics models were similar for both automated and manual delta-radiomics models.

Dynamics of Daily Glioblastoma Evolution During Chemoradiation Therapy on the 0.35T Magnetic Resonance Imaging-Linear Accelerator.

PURPOSE: Glioblastoma changes during chemoradiation therapy are inferred from magnetic resonance imaging (MRI) before and after treatment but are rarely investigated due to logistics of frequent MRI. Using a combination MRI-linear accelerator (MRI-linac), we evaluated changes during daily chemoradiation therapy. METHODS AND MATERIALS: Patients with glioblastoma were prospectively imaged daily during chemoradiation therapy on 0.35T MRI-linac and at 3 timepoints with and without contrast on standalone high-field MRI. Tumor or edema (lesion) and resection cavity dynamics throughout treatment were analyzed and compared with standalone T1 postcontrast (T1+C) and T2 volumes. RESULTS: Of 36 patients included in this analysis, 8 had cavity only, 12 had lesion only, and 16 had both cavity and lesion. Of these, 64% had lesion growth and 46% had cavity shrinkage during treatment on MRI-linac scans. The average MRI-linac migration distance was 1.3 cm (range, 0-4.1 cm) for lesion and 0.6 cm (range, 0.1-2.1 cm) for cavity. Standalone versus MRI-linac volumes correlated strongly with R2 values: 0.991 (T2 vs MRI-linac cavity), 0.972 (T1+C vs MRI-linac cavity), and 0.973 (T2 vs MRI-linac lesion). There was a moderate correlation between T1+C and MRI-linac lesion (R2 = 0.609), despite noncontrast MRI-linac inability to separate contrast enhancement from surrounding nonenhancing tumor and edema. From pretreatment to posttreatment in patients with all available scans (n = 35), T1+C and MRI-linac lesions changed together-shrank (n = 6), grew (n = 12), or unchanged (n = 8)-in 26 (74%) patients. Another 9 patients (26%) had growth on MRI-linac, although the T1+C component shrank. In no patient did T1+C lesion grow while MRI-linac lesion shrank. CONCLUSIONS: Anatomic changes are seen in patients with glioblastoma imaged daily on MRI-linac throughout the chemoradiation therapy course. As surgical resection cavities shrink, margins may be reduced to save normal brain. Patients with unresected or growing lesions may require margin expansions to cover changes. Limited volume glioblastoma boost trials could consider triggered gadolinium contrast administration for evaluation of adaptive radiation therapy when lesion growth is seen on noncontrast MRI-linac.

In vivo active-targeting fluorescence molecular imaging with adaptive background fluorescence subtraction.

Using active tumor-targeting nanoparticles, fluorescence imaging can provide highly sensitive and specific tumor detection, and precisely guide radiation in translational radiotherapy study. However, the inevitable presence of non-specific nanoparticle uptake throughout the body can result in high levels of heterogeneous background fluorescence, which limits the detection sensitivity of fluorescence imaging and further complicates the early detection of small cancers. In this study, background fluorescence emanating from the baseline fluorophores was estimated from the distribution of excitation light transmitting through tissues, by using linear mean square error estimation. An adaptive masked-based background subtraction strategy was then implemented to selectively refine the background fluorescence subtraction. First, an in vivo experiment was performed on a mouse intratumorally injected with passively targeted fluorescent nanoparticles, to validate the reliability and robustness of the proposed method in a stringent situation wherein the target fluorescence was overlapped with the strong background. Then, we conducted in vivo studies on 10 mice which were inoculated with orthotopic breast tumors and intravenously injected with actively targeted fluorescent nanoparticles. Results demonstrated that active targeting combined with the proposed background subtraction method synergistically increased the accuracy of fluorescence molecular imaging, affording sensitive tumor detection.

A Deep Learning Approach for Automatic Segmentation during Daily MRI-Linac Radiotherapy of Glioblastoma.

Glioblastoma changes during chemoradiotherapy are inferred from high-field MRI before and after treatment but are rarely investigated during radiotherapy. The purpose of this study was to develop a deep learning network to automatically segment glioblastoma tumors on daily treatment set-up scans from the first glioblastoma patients treated on MRI-linac. Glioblastoma patients were prospectively imaged daily during chemoradiotherapy on 0.35T MRI-linac. Tumor and edema (tumor lesion) and resection cavity kinetics throughout the treatment were manually segmented on these daily MRI. Utilizing a convolutional neural network, an automatic segmentation deep learning network was built. A nine-fold cross-validation schema was used to train the network using 80:10:10 for training, validation, and testing. Thirty-six glioblastoma patients were imaged pre-treatment and 30 times during radiotherapy (n = 31 volumes, total of 930 MRIs). The average tumor lesion and resection cavity volumes were 94.56 ± 64.68 cc and 72.44 ± 35.08 cc, respectively. The average Dice similarity coefficient between manual and auto-segmentation for tumor lesion and resection cavity across all patients was 0.67 and 0.84, respectively. This is the first brain lesion segmentation network developed for MRI-linac. The network performed comparably to the only other published network for auto-segmentation of post-operative glioblastoma lesions. Segmented volumes can be utilized for adaptive radiotherapy and propagated across multiple MRI contrasts to create a prognostic model for glioblastoma based on multiparametric MRI.

Predictive Value of Delta-Radiomics Texture Features in 0.35 Tesla Magnetic Resonance Setup Images Acquired During Stereotactic Ablative Radiotherapy of Pancreatic Cancer.

Purpose: The purpose of this work is to explore delta-radiomics texture features for predicting response using setup images of pancreatic cancer patients treated with magnetic resonance image guided (MRI-guided) stereotactic ablative radiotherapy (SBRT). Methods: The total biological effective dose (BED) was calculated for 30 patients treated with MRI-guided SBRT that delivered physical doses of 30-60 Gy in three to five fractions. Texture features were then binned into groups based upon BED per fraction by dividing BED by the number of fractions. Delta-radiomics texture features were calculated after delivery of 20 Gy BED (BED20 features) and 40 Gy BED (BED40 features). A random forest (RF) model was constructed using BED20 and then BED40 features to predict binary outcome. During model training, the Gini Index, a measure of a variable's importance for accurate prediction, was calculated for all features, and the two features that ranked the highest were selected for internal validation. The two features selected from each bin were used in a bootstrapped logistic regression model to predict response and performance quantified using the area under the receiver operating characteristic curve (AUC). This process was an internal validation analysis. Results: After RF model training, the Gini Index was highest for gray-level co-occurrence matrix-based (GLCM) sum average, and neighborhood gray tone difference matrix-based (NGTDM) busyness for BED20 features and gray-level size zone matrix-based (GLSZM) large zones low gray-level emphasis and gray-level run length matrix-based (GLRLM) run percentage was selected from the BED40-based features. The mean AUC obtained using the two BED20 features was AUC = 0.845 with the 2.5 percentile and 97.5 percentile values ranging from 0.794 to 0.856. Internal validation of the BED40 delta-radiomics features resulted in a mean AUC = 0.567 with a 2.5 and 97.5 percentile range of 0.502-0.675. Conclusion: Early changes in treatment quantified with the BED20 delta-radiomics texture features in low field images acquired during MRI-guided SBRT demonstrated better performance in internal validation than features calculated later in treatment. Further analysis of delta-radiomics texture analysis in low field MRI is warranted.

MRI-based delta-radiomic features for prediction of local control in liver lesions treated with stereotactic body radiation therapy.

Real-time magnetic resonance image guided stereotactic ablative radiotherapy (MRgSBRT) is used to treat abdominal tumors. Longitudinal data is generated from daily setup images. Our study aimed to identify delta radiomic texture features extracted from these images to predict for local control in patients with liver tumors treated with MRgSBRT. Retrospective analysis of an IRB-approved database identified patients treated with MRgSBRT for primary liver and secondary metastasis histologies. Daily low field strength (0.35 T) images were retrieved, and the gross tumor volume was identified on each image. Next, images' gray levels were equalized, and 39 s-order texture features were extracted. Delta-radiomics were calculated as the difference between feature values on the initial scan and after delivered biological effective doses (BED, α/ß = 10) of 20 Gy and 40 Gy. Then, features were ranked by the Gini Index during training of a random forest model. Finally, the area under the receiver operating characteristic curve (AUC) was estimated using a bootstrapped logistic regression with the top two features. We identified 22 patients for analysis. The median dose delivered was 50 Gy in 5 fractions. The top two features identified after delivery of BED 20 Gy were gray level co-occurrence matrix features energy and gray level size zone matrix based large zone emphasis. The model generated an AUC = 0.9011 (0.752-1.0) during bootstrapped logistic regression. The same two features were selected after delivery of a BED 40 Gy, with an AUC = 0.716 (0.600-0.786). Delta-radiomic features after a single fraction of SBRT predicted local control in this exploratory cohort. If confirmed in larger studies, these features may identify patients with radioresistant disease and provide an opportunity for physicians to alter management much sooner than standard restaging after 3 months. Expansion of the patient database is warranted for further analysis of delta-radiomic features.

Quantifying Radiosensitization of PSMA-Targeted Gold Nanoparticles on Prostate Cancer Cells at Megavoltage Radiation Energies by Monte Carlo Simulation and Local Effect Model.

Active targeting gold nanoparticles (AuNPs) are a very promising avenue for cancer treatment with many publications on AuNP mediated radiosensitization at kilovoltage (kV) photon energies. However, uncertainty on the effectiveness of AuNPs under clinically relevant megavoltage (MV) radiation energies hinders the clinical translation of AuNP-assisted radiation therapy (RT) paradigm. The aim of this study was to investigate radiosensitization mediated by PSMA-targeted AuNPs irradiated by a 6 MV radiation beam at different depths to explore feasibility of AuNP-assisted prostate cancer RT under clinically relevant conditions. PSMA-targeted AuNPs (PSMA-AuNPs) were synthesized by conjugating PSMA antibodies onto PEGylated AuNPs through EDC/NHS chemistry. Confocal fluorescence microscopy was used to verify the active targeting of the developed PSMA-AuNPs. Transmission electron microscopy (TEM) was used to demonstrate the intracellular biodistribution of PSMA-AuNPs. LNCaP prostate cancer cells treated with PSMA-AuNPs were irradiated on a Varian 6 MV LINAC under varying depths (2.5 cm, 10 cm, 20 cm, 30 cm) of solid water. Clonogenic assays were carried out to determine the in vitro cell survival fractions. A Monte Carlo (MC) model developed on TOPAS platform was then employed to determine the nano-scale radial dose distribution around AuNPs, which was subsequently used to predict the radiation dose response of LNCaP cells treated with AuNPs. Two different cell models, with AuNPs located within the whole cell or only in the cytoplasm, were used to assess how the intracellular PSMA-AuNP biodistribution impacts the prostate cancer radiosensitization. Then, MC-based microdosimetry was combined with the local effect model (LEM) to calculate cell survival fraction, which was benchmarked against the in vitro clonogenic assays at different depths. In vitro clonogenic assay of LNCaP cells demonstrated the depth dependence of AuNP radiosensitization under clinical megavoltage beams, with sensitization enhancement ratio (SER) of 1.14 ± 0.03 and 1.55 ± 0.05 at 2.5 cm depth and 30 cm depth, respectively. The MC microdosimetry model showed the elevated percent of low-energy photons in the MV beams at greater depth, consequently resulting in increased dose enhancement ratio (DER) of AuNPs with depth. The AuNP-induced DER reached ~5.7 and ~8.1 at depths of 2.5 cm and 30 cm, respectively. Microdosimetry based LEM accurately predicted the cell survival under 6 MV beams at different depths, for the cell model with AuNPs placed only in the cell cytoplasm. TEM results demonstrated the distribution of PSMA-AuNPs in the cytoplasm, confirming the accuracy of MC microdosimetry based LEM with modelled AuNPs distributed within the cytoplasm. We conclude that AuNP radiosensitization can be achieved under megavoltage clinical radiotherapy energies with a dependence on tumor depth. Furthermore, the combination of Monte Carlo microdosimetry and LEM will be a valuable tool to assist with developing AuNP-aided radiotherapy paradigm and drive clinical translation.

Cone-beam CT delta-radiomics to predict genitourinary toxicities and international prostate symptom of prostate cancer patients: a pilot study.

For prostate cancer (PCa) patients treated with definitive radiotherapy (RT), acute and late RT-related genitourinary (GU) toxicities adversely impact disease-specific quality of life. Early warning of potential RT toxicities can prompt interventions that may prevent or mitigate future adverse events. During intensity modulated RT (IMRT) of PCa, daily cone-beam computed tomography (CBCT) images are used to improve treatment accuracy through image guidance. This work investigated the performance of CBCT-based delta-radiomic features (DRF) models to predict acute and sub-acute International Prostate Symptom Scores (IPSS) and Common Terminology Criteria for Adverse Events (CTCAE) version 5 GU toxicity grades for 50 PCa patients treated with definitive RT. Delta-radiomics models were built using logistic regression, random forest for feature selection, and a 1000 iteration bootstrapping leave one analysis for cross validation. To our knowledge, no prior studies of PCa have used DRF models based on daily CBCT images. AUC of 0.83 for IPSS and greater than 0.7 for CTCAE grades were achieved as early as week 1 of treatment. DRF extracted from CBCT images showed promise for the development of models predictive of RT outcomes. Future studies will include using artificial intelligence and machine learning to expand CBCT sample sizes available for radiomics analysis.

Deep Learning for Per-Fraction Automatic Segmentation of Gross Tumor Volume (GTV) and Organs at Risk (OARs) in Adaptive Radiotherapy of Cervical Cancer.

Background/Hypothesis: MRI-guided online adaptive radiotherapy (MRI-g-OART) improves target coverage and organs-at-risk (OARs) sparing in radiation therapy (RT). For patients with locally advanced cervical cancer (LACC) undergoing RT, changes in bladder and rectal filling contribute to large inter-fraction target volume motion. We hypothesized that deep learning (DL) convolutional neural networks (CNN) can be trained to accurately segment gross tumor volume (GTV) and OARs both in planning and daily fractions' MRI scans. Materials/Methods: We utilized planning and daily treatment fraction setup (RT-Fr) MRIs from LACC patients, treated with stereotactic body RT to a dose of 45-54 Gy in 25 fractions. Nine structures were manually contoured. MASK R-CNN network was trained and tested under three scenarios: (i) Leave-one-out (LOO), using the planning images of N- 1 patients for training; (ii) the same network, tested on the RT-Fr MRIs of the "left-out" patient, (iii) including the planning MRI of the "left-out" patient as an additional training sample, and tested on RT-Fr MRIs. The network performance was evaluated using the Dice Similarity Coefficient (DSC) and Hausdorff distances. The association between the structures' volume and corresponding DSCs was investigated using Pearson's Correlation Coefficient, r. Results: MRIs from fifteen LACC patients were analyzed. In the LOO scenario the DSC for Rectum, Femur, and Bladder was >0.8, followed by the GTV, Uterus, Mesorectum and Parametrium (0.6-0.7). The results for Vagina and Sigmoid were suboptimal. The performance of the network was similar for most organs when tested on RT-Fr MRI. Including the planning MRI in the training did not improve the segmentation of the RT-Fr MRI. There was a significant correlation between the average organ volume and the corresponding DSC (r = 0.759, p = 0.018). Conclusion: We have established a robust workflow for training MASK R-CNN to automatically segment GTV and OARs in MRI-g-OART of LACC. Albeit the small number of patients in this pilot project, the network was trained to successfully identify several structures while challenges remain, especially in relatively small organs. With the increase of the LACC cases, the performance of the network will improve. A robust auto-contouring tool would improve workflow efficiency and patient tolerance of the OART process.

Bow-tie wobble artifact: effect of source assembly motion on cone-beam CT.

PURPOSE: To investigate the cause of a bow-tie wobble artifact (BWA) discovered on Varian OBI CBCT images and to develop practical correction strategies. METHOD AND MATERIALS: The dependence of the BWA on phantom geometry, phantom position, specific system, and reconstruction algorithm was investigated. Simulations were conducted to study the dependence of the BWA on scatter and beam hardening corrections. Geometric calibration was performed to rule out other gantry-angle dependent mechanical non-idealities as BWA causes. Air scans were acquired with ball-bearing markers to study the motions of the x-ray head assembly as functions of gantry angle. Based on measurements, we developed hypothesis regarding the BWA cause. Simulations were performed to validate our hypothesis. Two correction strategies were implemented: a measurement-based method, which acquires gantry-dependent normalization projections (NPs); and a model-based method that involves numerically shifting the single-angle NP to compensate for the previously-measured bow-tie-filter (BTF) motion. RESULTS: The BWA has a diameter of approximately 15 cm, is centered at the isocenter, and is reproducible independent of phantom, position, system, reconstruction, and standard corrections, but only when the BTF is used. Measurements identified a 2D sinusoidal gantry-angle-dependent motion of the x-ray head assembly, and it was the BTF motion (>3 mm amplitude projected onto the detector) resulting an intensity mismatch between the all-angle CBCT projections and a single-angle NP that caused the BWA. Both correction strategies were demonstrated effective. CONCLUSIONS: A geometric mismatch between the BTF modulation patterns on CBCT projections and on the NP causes the BWA. The BTF wobble requires additional degrees of freedom in CBCT geometric calibration to characterize.

MR-Guided Radiotherapy for Brain and Spine Tumors.

MRI is the standard modality to assess anatomy and response to treatment in brain and spine tumors given its superb anatomic soft tissue contrast (e.g., T1 and T2) and numerous additional intrinsic contrast mechanisms that can be used to investigate physiology (e.g., diffusion, perfusion, spectroscopy). As such, hybrid MRI and radiotherapy (RT) devices hold unique promise for Magnetic Resonance guided Radiation Therapy (MRgRT). In the brain, MRgRT provides daily visualizations of evolving tumors that are not seen with cone beam CT guidance and cannot be fully characterized with occasional standalone MRI scans. Significant evolving anatomic changes during radiotherapy can be observed in patients with glioblastoma during the 6-week fractionated MRIgRT course. In this review, a case of rapidly changing symptomatic tumor is demonstrated for possible therapy adaptation. For stereotactic body RT of the spine, MRgRT acquires clear isotropic images of tumor in relation to spinal cord, cerebral spinal fluid, and nearby moving organs at risk such as bowel. This visualization allows for setup reassurance and the possibility of adaptive radiotherapy based on anatomy in difficult cases. A review of the literature for MR relaxometry, diffusion, perfusion, and spectroscopy during RT is also presented. These techniques are known to correlate with physiologic changes in the tumor such as cellularity, necrosis, and metabolism, and serve as early biomarkers of chemotherapy and RT response correlating with patient survival. While physiologic tumor investigations during RT have been limited by the feasibility and cost of obtaining frequent standalone MRIs, MRIgRT systems have enabled daily and widespread physiologic measurements. We demonstrate an example case of a poorly responding tumor on the 0.35 T MRIgRT system with relaxometry and diffusion measured several times per week. Future studies must elucidate which changes in MR-based physiologic metrics and at which timepoints best predict patient outcomes. This will lead to early treatment intensification for tumors identified to have the worst physiologic responses during RT in efforts to improve glioblastoma survival.

Assessment of CT to CBCT contour mapping for radiomic feature analysis in prostate cancer.

This study provides a quantitative assessment of the accuracy of a commercially available deformable image registration (DIR) algorithm to automatically generate prostate contours and additionally investigates the robustness of radiomic features to differing contours. Twenty-eight prostate cancer patients enrolled on an institutional review board (IRB) approved protocol were selected. Planning CTs (pCTs) were deformably registered to daily cone-beam CTs (CBCTs) to generate prostate contours (auto contours). The prostate contours were also manually drawn by a physician. Quantitative assessment of deformed versus manually drawn prostate contours on daily CBCT images was performed using Dice similarity coefficient (DSC), mean distance-to-agreement (MDA), difference in center-of-mass position (ΔCM) and difference in volume (ΔVol). Radiomic features from 6 classes were extracted from each contour. Lin's concordance correlation coefficient (CCC) and mean absolute percent difference in radiomic feature-derived data (mean |%Δ|RF) between auto and manual contours were calculated. The mean (± SD) DSC, MDA, ΔCM and ΔVol between the auto and manual prostate contours were 0.90 ± 0.04, 1.81 ± 0.47 mm, 2.17 ± 1.26 mm and 5.1 ± 4.1% respectively. Of the 1,010 fractions under consideration, 94.8% of DIRs were within TG-132 recommended tolerance. 30 radiomic features had a CCC > 0.90 and 21 had a mean |%∆|RF < 5%. Auto-propagation of prostate contours resulted in nearly 95% of DIRs within tolerance recommendations of TG-132, leading to the majority of features being regarded as acceptably robust. The use of auto contours for radiomic feature analysis is promising but must be done with caution.

Repeatability of CBCT radiomic features and their correlation with CT radiomic features for prostate cancer.

PURPOSE: Radiomic features of cone-beam CT (CBCT) images have potential as biomarkers to predict treatment response and prognosis for patients of prostate cancer. Previous studies of radiomic feature analysis for prostate cancer were assessed in a variety of imaging modalities, including MRI, PET, and CT, but usually limited to a pretreatment setting. However, CBCT images may provide an opportunity to capture early morphological changes to the tumor during treatment that could lead to timely treatment adaptation. This work investigated the quality of CBCT-based radiomic features and their relationship with reconstruction methods applied to the CBCT projections and the preprocessing methods used in feature extraction. Moreover, CBCT features were correlated with planning CT (pCT) features to further assess the viability of CBCT radiomic features. METHODS: The quality of 42 CBCT-based radiomic features was assessed according to their repeatability and reproducibility. Repeatability was quantified by correlating radiomic features between 20 CBCT scans that also had repeated scans within 15 minutes. Reproducibility was quantified by correlating radiomic features between the planning CT (pCT) and the first fraction CBCT for 20 patients. Concordance correlation coefficients (CCC) of radiomic features were used to estimate the repeatability and reproducibility of radiomic features. The same patient dataset was assessed using different reconstruction methods applied to the CBCT projections. CBCT images were generated using 18 reconstruction methods using iterative (iCBCT) and standard (sCBCT) reconstructions, three convolution filters, and five noise suppression filters. Eighteen preprocessing settings were also considered. RESULTS: Overall, CBCT radiomic features were more repeatable than reproducible. Five radiomic features are repeatable in > 97% of the reconstruction and preprocessing methods, and come from the gray-level size zone matrix (GLSZM), neighborhood gray-tone difference matrix (NGTDM), and gray-level-run length matrix (GLRLM) radiomic feature classes. These radiomic features were reproducible in > 9.8% of the reconstruction and preprocessing methods. Noise suppression and convolution filter smoothing increased radiomic features repeatability, but decreased reproducibility. The top-repeatable iCBCT method (iCBCT-Sharp-VeryHigh) is more repeatable than the top-repeatable sCBCT method (sCBCT-Smooth) in 64% of the radiomic features. CONCLUSION: Methods for reconstruction and preprocessing that improve CBCT radiomic feature repeatability often decrease reproducibility. The best approach may be to use methods that strike a balance repeatability and reproducibility such as iCBCT-Sharp-VeryLow-1-Lloyd-256 that has 17 repeatable and eight reproducible radiomic features. Previous radiomic studies that only used pCT radiomic features have generated prognostic models of prostate cancer outcome. Since our study indicates that CBCT radiomic features correlated well with a subset of pCT radiomic features, one may expect CBCT radiomics to also generate prognostic models for prostate cancer.

Impact of quantization algorithm and number of gray level intensities on variability and repeatability of low field strength magnetic resonance image-based radiomics texture features.

PURPOSE: The purpose of this work was to investigate the impact of quantization preprocessing parameter selection on variability and repeatability of texture features derived from low field strength magnetic resonance (MR) images. METHODS: Texture features were extracted from low field strength images of a daily image QA phantom with four texture inserts. Feature variability over time was quantified using all combinations of three quantization algorithms and four different numbers of gray level intensities. In addition, texture features were extracted using the same combinations from the low field strength MR images of the gross tumor volume (GTV) and left kidney of patients with repeated set up scans. The impact of region of interest (ROI) preprocessing on repeatability was investigated with a test-retest study design. RESULTS: The phantom ROIs quantized to 64 Gy level intensities using the histogram equalization method resulted in the greatest number of features with the least variability. There was no clear method that resulted in the highest repeatability in the GTV or left kidney. However, eight texture features extracted from the GTV were repeatable regardless of ROI processing combination. CONCLUSION: Low field strength MR images can provide a stable basis for texture analysis with ROIs quantized to 64 Gy levels using histogram equalization, but there is no clear optimal combination for repeatability.

Predictive value of 0.35 T magnetic resonance imaging radiomic features in stereotactic ablative body radiotherapy of pancreatic cancer: A pilot study.

PURPOSE: The aim of this study was to evaluate the potential and feasibility of radiomic features extracted from low field strength (0.35 T) magnetic resonance images (MRIs) in predicting treatment response for patients with pancreatic cancer undergoing stereotactic body radiotherapy (SBRT). METHODS: Twenty patients with unresected, non-metastatic pancreatic ductal adenocarcinoma (PDAC) were enrolled, all of whom received neoadjuvant chemotherapy followed by five-fraction MR-guided SBRT with a radiation dose range of 33-50 Gy. For each patient, five daily setup scans were acquired from a hybrid 0.35 T MRI/radiotherapy unit. Tumor heterogeneity quantified with radiomic features extracted from the gross tumor volume (GTV) was averaged over the course of treatment. Random forest (RF) and adaptive least absolute shrinkage and selection operator (LASSO) classification models were constructed to identify radiomics features predictive of treatment response. Predictive capability of the top-performing features was then evaluated using the receiver operating characteristic area under curve (AUC) obtained using leave-one-out cross-validation. RESULTS: Half of the 20 patients showed response to treatment, defined by tumor regression on histopathology or tumor response on follow-up dynamic contrast-enhanced computed tomography (CT). The most predictive features selected by the RF method were GLCM energy and GLSZM gray-level variance. The RF-based model achieved an AUC = 0.81 with a 95% confidence interval of [0.594 to 1] The LASSO algorithm selected GLCM energy as the only predictive feature, achieving an AUC = 0.81 with 95% confidence interval of [0.596 to 1]. CONCLUSION: The findings of this study suggest that radiomic features extracted during MR-guided SBRT may contain predictive information about response of PDAC patients to treatment. Using the images acquired during treatment of PDAC patients supports continued expansion of radiomic analysis based on low field strength MR images and may hold the potential for providing timely indications of response to treatment.

MRI-guided stereotactic ablative radiation therapy of spinal bone metastases: a preliminary experience.

OBJECTIVE: MRI provides clear visualization of spinal cord, tumor, and bone for patient positioning and verification during MRI-guided radiotherapy (MRI-RT). Therefore, we wished to evaluate spine stereotactic ablative radiotherapy (SABR) feasibility with MRI-RT. Given dosimetric limitations of first generation Co-60 MRI-RT, we then evaluated improvements by newer linear accelerator (linac) MRI-RT. METHODS: Nine spinal metastases were treated with Co-60 MRI-RT. Seven received a single 16 Gy fraction, and two received three fractions totaling 24 or 30 Gy. After replanning with linac MRI-RT software, comparisons of organ at risk and dose spillage objectives between Co-60 and linac plans were performed. RESULTS: Spinal cord and cauda equina dose constraints were met in all Co-60 cases. Treatments were delivered successfully with real-time imaging during treatment and no treatment-related toxicities. While limits for dose spillage into surrounding soft tissues were not achieved due to the limitations of the Co-60 system, this could be corrected with linac MRI-RT delivery. CONCLUSIONS: MRI-RT SABR of spinal metastases is feasible with Co-60 MRI-RT. Dose delivery is improved by linac MRI-RT. ADVANCES IN KNOWLEDGE: This is the first report of MRI-RT for SABR of spinal metastases. The enhanced visualization of anatomy by MRI may facilitate RT dose escalation for spine SABR.