Dose-volume histogram prediction in volumetric modulated arc therapy for nasopharyngeal carcinomas based on uniform-intensity radiation with equal angle intervals

Application of Artificial Intelligence to the Diagnosis and Therapy of Nasopharyngeal Carcinoma

Artificial intelligence (AI) is an interdisciplinary field that encompasses a wide range of computer science disciplines, including image recognition, machine learning, human-computer interaction, robotics and so on. Recently, AI, especially deep learning algorithms, has shown excellent performance in the field of image recognition, being able to automatically perform quantitative evaluation of complex medical image features to improve diagnostic accuracy and efficiency. AI has a wider and deeper application in the medical field of diagnosis, treatment and prognosis. Nasopharyngeal carcinoma (NPC) occurs frequently in southern China and Southeast Asian countries and is the most common head and neck cancer in the region. Detecting and treating NPC early is crucial for a good prognosis. This paper describes the basic concepts of AI, including traditional machine learning and deep learning algorithms, and their clinical applications of detecting and assessing NPC lesions, facilitating treatment and predicting prognosis. The main limitations of current AI technologies are briefly described, including interpretability issues, privacy and security and the need for large amounts of annotated data. Finally, we discuss the remaining challenges and the promising future of using AI to diagnose and treat NPC.

Deep learning architecture with transformer and semantic field alignment for voxel-level dose prediction on brain tumors

PURPOSE

The use of convolution neural networks (CNN) to accurately predict dose distributions can accelerate intensity-modulated radiation therapy (IMRT) planning. The purpose of our study is to develop a novel deep learning architecture for precise voxel-level dose prediction on brain tumors.

METHODS

A dataset of 120 patients with brain tumors is built for the retrospective study. The dose distributions are predicted by a designed end-to-end model called TS-Net, in which the transformer encoder module is utilized to obtain abundant global features by learning long-range correlations of the input sequence. In addition, semantic field alignment (SFA) block is proposed in decoding path to ensure effective propagation of strong semantic information from deep to shallow. Five images from different channels are fed into the architecture, including a computed tomography (CT) image, a planning target volumes (PTV) image, an organs-at-risk (OARs) image, a beam configuration image, and a distance image, and the predicted dose distributions are taken as outputs. We use different evaluation metrics to evaluate the performance of the model and discuss the role of the auxiliary beam configuration information provided by non-modulated dose distributions.

RESULTS

The TS-Net prediction accuracies in terms of mean absolute error (MAE) are 2.98% for PTV, 7.19% for brainstem, 1.88% for left len, 2.48% for right len, 9.61% for left optic nerve, 9.10% for right optic nerve, 8.99% for optic chiasma, and 8.28% for pituitary. There is no statistically significant difference between the predicted results and clinical dose distributions for clinical indexes including homogeneity index (HI), D50, and D95 for PTV; V40, mean dose, and max dose for OARs; except for conformation index (CI) and D2 for PTV. The model has dice similarity coefficient (DSC) values of above 0.91 for most isodose volumes, clearly outperforming HD U-Net, and being slightly better than U-Net and DCNN.

CONCLUSION

The proposed TS-Net with beam configuration input can achieve accurate voxel-level dose prediction for brain tumors, and is a usable tool for improving the efficiency and quality of radiotherapy.

Knowledge-based planning using both the predicted DVH of organ-at risk and planning target volume

PURPOSE

The purpose of this study was to evaluate the performance of a knowledge-based planning (KBP) method in nasopharyngeal cancer radiotherapy using the predicted dose-volume histogram (DVH) of organ-at risk (OAR) and planning target volume (PTV).

METHODS AND MATERIALS

A total of 85 patients previously treated for nasopharyngeal cancer using 9-field 6-MV intensity-modulated radiation therapy (IMRT) were identified for training and 30 similar patients were identified for testing. The dosimetric deposition information, individual dose-volume histograms (IDVHs) induced by a series of fields with uniform-intensity irradiation, was used to predict both OAR and PTV DVH. Two KBP methods (KBP_OAR and KBP_OAR+PTV) were established for plan generation based on the DVH prediction. The KBP_OAR method utilized the dose constraints based on the predicted OAR DVH and the PTV dose constraints obtained according to the planning experience, while the KBP_OAR+PTV method applied the dose constraints based on the predicted OAR and PTV DVH. For the plan evaluation, the PTV dose coverage was used D₉₈ and D₂, and the maximum dose, mean dose or dose-volume parameters were used for the OARs. Statistical differences of the two KBP methods were tested with the Wilcoxon signed rank test.

RESULTS

For patients with T3 tumors, there was no significant difference between the KBP_OAR and KBP_OAR+PTV methods in dosimetric results at most OARs and PTVs. Both KBP methods achieved a similar number of plans meeting the dose requirements. For patients with T4 tumors, KBP_OAR+PTV reduced the maximum dose by more than 1 Gy in the body, spinal cord, optic nerve, eye and temporal lobes and reduced the V₅₀ value by more than 3.9% in the larynx and tongue without reducing the PTV dose compared with KBP_OAR. The KBP_OAR+PTV method increased the plans by more than 14.2% in meeting the maximum dose requirements at the body, optic nerve, mandible and eye and increased the plans by more than 21.4% in meeting the V₅₀ of the larynx and V₅₀ of the tongue when compared with the KBP_OAR method.

CONCLUSIONS

For patients with T3 tumors, no significant difference was found between the KBP_OAR and KBP_OAR+PTV methods in dosimetric results at most OARs and PTVs. For patients with T4 tumors, the KBP_OAR+PTV method performs better than the KBP_OAR method in improving the quality of the plans. Compared with the KBP_OAR method, dose sparing of some OARs was achieved without reducing PTV dose coverage and helped to increase the number of plans meeting the dose requirements when the KBP_OAR+PTV method was utilized.

A recurrent neural network for rapid detection of delivery errors during real-time portal dosimetry

Background and purpose

Real-time portal dosimetry compares measured images with predicted images to detect delivery errors as the radiotherapy treatment proceeds. This work aimed to investigate the performance of a recurrent neural network for processing image metrics so as to detect delivery errors as early as possible in the treatment.

Materials and methods

Volumetric modulated arc therapy (VMAT) plans of six prostate patients were used to generate sequences of predicted portal images. Errors were introduced into the treatment plans and the modified plans were delivered to a water-equivalent phantom. Four different metrics were used to detect errors. These metrics were applied to a threshold-based method to detect the errors as soon as possible during the delivery, and also to a recurrent neural network consisting of four layers. A leave-two-out approach was used to set thresholds and train the neural network then test the resulting systems.

Results

When using a combination of metrics in conjunction with optimal thresholds, the median segment index at which the errors were detected was 107 out of 180. When using the neural network, the median segment index for error detection was 66 out of 180, with no false positives. The neural network reduced the rate of false negative results from 0.36 to 0.24.

Conclusions

The recurrent neural network allowed the detection of errors around 30% earlier than when using conventional threshold techniques. By appropriate training of the network, false positive alerts could be prevented, thereby avoiding unnecessary disruption to the patient workflow.

The status of medical physics in radiotherapy in China

PURPOSE

To present an overview of the status of medical physics in radiotherapy in China, including facilities and devices, occupation, education, research, etc. MATERIALS AND METHODS: The information about medical physics in clinics was obtained from the 9-th nationwide survey conducted by the China Society for Radiation Oncology in 2019. The data of medical physics in education and research was collected from the publications of the official and professional organizations.

RESULTS

By 2019, there were 1463 hospitals or institutes registered to practice radiotherapy and the number of accelerators per million population was 1.5. There were 4172 medical physicists working in clinics of radiation oncology. The ratio between the numbers of radiation oncologists and medical physicists is 3.51. Approximately, 95% of medical physicists have an undergraduate or graduate degrees in nuclear physics and biomedical engineering. 86% of medical physicists have certificates issued by the Chinese Society of Medical Physics. There has been a fast growth of publications by authors from mainland of China in the top international medical physics and radiotherapy journals since 2018.

CONCLUSIONS

Demand for medical physicists in radiotherapy increased quickly in the past decade. The distribution of radiotherapy facilities in China became more balanced. High quality continuing education and training programs for medical physicists are deficient in most areas. The role of medical physicists in the clinic has not been clearly defined and their contributions have not been fully recognized by the community.

Evaluation of dose-volume histogram prediction for organ-at risk and planning target volume based on machine learning

The purpose of this work is to evaluate the performance of applying patient dosimetric information induced by individual uniform-intensity radiation fields in organ-at risk (OAR) dose-volume histogram (DVH) prediction, and extend to DVH prediction of planning target volume (PTV). Ninety nasopharyngeal cancer intensity-modulated radiation therapy (IMRT) plans and 60 rectal cancer volumetric modulated arc therapy (VMAT) plans were employed in this study. Of these, 20 nasopharyngeal cancer cases and 15 rectal cancer cases were randomly selected as the testing data. The DVH prediction was performed using two methods. One method applied the individual dose-volume histograms (IDVHs) induced by a series of fields with uniform-intensity irradiation and the other method applied the distance-to-target histogram and the conformal-plan-dose-volume histogram (DTH + CPDVH). The determination coefficient R² and mean absolute error (MAE) were used to evaluate DVH prediction accuracy. The PTV DVH prediction was performed using the IDVHs. The PTV dose coverage was evaluated using D₉₈, D₉₅, D₁ and uniformity index (UI). The OAR dose was compared using the maximum dose, V₃₀ and V₄₀. The significance of the results was examined with the Wilcoxon signed rank test. For PTV DVH prediction using IDVHs, the clinical plan and IDVHs prediction method achieved mean UI values of 1.07 and 1.06 for nasopharyngeal cancer, and 1.04 and 1.05 for rectal cancer, respectively. No significant difference was found between the clinical plan results and predicted results using the IDVHs method in achieving PTV dose coverage (D_98,D_95,D₁ and UI) for both nasopharyngeal cancer and rectal cancer (p-values ≥ 0.052). For OAR DVH prediction, no significant difference was found between the IDVHs and DTH + CPDVH methods for the R², MAE, the maximum dose, V₃₀ and V₄₀ (p-values ≥ 0.087 for all OARs). This work evaluates the performance of dosimetric information of several individual fields with uniform-intensity radiation for DVH prediction, and extends its application to PTV DVH prediction. The results indicated that the IDVHs method is comparable to the DTH + CPDVH method in accurately predicting the OAR DVH. The IDVHs method quantified the input features of the PTV and showed reliable PTV DVH prediction, which is helpful for plan quality evaluation and plan generation.

DVH Prediction for VMAT in NPC with GRU-RNN: An Improved Method by Considering Biological Effects

Purpose

A recurrent neural network (RNN) and its variants such as gated recurrent unit-based RNN (GRU-RNN) were found to be very suitable for dose-volume histogram (DVH) prediction in our previously published work. Using the dosimetric information generated by nonmodulated beams of different orientations, the GRU-RNN model was capable of accurate DVH prediction for nasopharyngeal carcinoma (NPC) treatment planning. On the basis of our previous work, we proposed an improved approach and aimed to further improve the DVH prediction accuracy as well as study the feasibility of applying the proposed method to relatively small-size patient data.

Methods

Eighty NPC volumetric modulated arc therapy (VMAT) plans with local IRB's approval in recent two years were retrospectively and randomly selected in this study. All these original plans were created using the Eclipse treatment planning system (V13.5, Varian Medical Systems, USA) with ≥95% of PGTVnx receiving the prescribed doses of 70 Gy, ≥95% of PGTVnd receiving 66 Gy, and ≥95% of PTV receiving 60 Gy. Among them, fifty plans were used to train the DVH prediction model, and the remaining were used for testing. On the basis of our previously published work, we simplified the 3-layer GRU-RNN model to a single-layer model and further trained every organ at risk (OAR) separately with an OAR-specific equivalent uniform dose- (EUD-) based loss function.

Results

The results of linear least squares regression obtained by the new proposed method showed the excellent agreements between the predictions and the original plans with the correlation coefficient r = 0.976 and 0.968 for EUD results and maximum dose results, respectively, and the coefficient r of our previously published method was 0.957 and 0.946, respectively. The Wilcoxon signed-rank test results between the proposed and the previous work showed that the proposed method could significantly improve the EUD prediction accuracy for the brainstem, spinal cord, and temporal lobes with a p value < 0.01.

Conclusions

The accuracy of DVH prediction achieved in different OARs showed the great improvements compared to the previous works, and more importantly, the effectiveness and robustness showed by the simplified GRU-RNN trained from relatively small-size DVH samples, fully demonstrated the feasibility of applying the proposed method to small-size patient data. Excellent agreements in both EUD results and maximum dose results between the predictions and original plans indicated the application prospect in a physically and biologically related (or a mixture of both) model for treatment planning.

Application of dose-volume histogram prediction in biologically related models for nasopharyngeal carcinomas treatment planning

Purpose

In this study, we employed a gated recurrent unit (GRU)-based recurrent neural network (RNN) using dosimetric information induced by individual beam to predict the dose-volume histogram (DVH) and investigated the feasibility and usefulness of this method in biologically related models for nasopharyngeal carcinomas (NPC) treatment planning.

Methods and materials

One hundred patients with NPC undergoing volumetric modulated arc therapy (VMAT) between 2018 and 2019 were randomly selected for this study. All the VMAT plans were created using the Monaco treatment planning system (Elekta, Sweden) and clinically approved: > 98% of PGTVnx received the prescribed doses of 70 Gy, > 98% of PGTVnd received the prescribed doses of 66 Gy and > 98% of PCTV received 60 Gy. Of these, the data from 80 patients were used to train the GRU-RNN, and the data from the other 20 patients were used for testing. For each NPC patient, the DVHs of different organs at risk were predicted by a trained GRU-based RNN using the information given by individual conformal beams. Based on the predicted DVHs, the equivalent uniform doses (EUD) were calculated and applied as dose constraints during treatment planning optimization. The regenerated VMAT experimental plans (EPs) were evaluated by comparing them with the clinical plans (CPs).

Results

For the 20 test patients, the regenerated EPs guided by the GRU-RNN predictive model achieved good consistency relative to the CPs. The EPs showed better consistency in PTV dose distribution and better dose sparing for many organs at risk, and significant differences were found in the maximum/mean doses to the brainstem, brainstem PRV, spinal cord, lenses, temporal lobes, parotid glands and larynx with P-values < 0.05. On average, compared with the CPs, the maximum/mean doses to these OARs were altered by − 3.44 Gy, − 1.94 Gy, − 1.88 Gy, 0.44 Gy, 1.98 Gy, − 1.82 Gy and 2.27 Gy, respectively. In addition, significant differences were also found in brainstem and spinal cord for the dose received by 1 cc volume with 4.11 and 1.67 Gy dose reduction in EPs on average.

Conclusion

The GRU-RNN-based DVH prediction method was capable of accurate DVH prediction. The regenerated plans guided by the predicted EUDs were not inferior to the manual plans, had better consistency in PTVs and better dose sparing in critical OARs, indicating the usefulness and effectiveness of biologically related model in knowledge-based planning.