Fixed-Field IMRT for Cervix Carcinoma Patients on an MR-LINAC Platform: Dosimetric Feasibility and Challenges

PURPOSE/OBJECTIVE(S)

To investigate the impact of MR-LINAC performance characteristics and inverse planning implementation on the feasibility of fixed-field IMRT for cervix carcinoma patients by benchmarking MR-LINAC plans against clinically used VMAT plans in a single institution study and multi-institutional treatment planning challenge.

MATERIALS/METHODS

For 10 cervix cancer patients who had previously received Linac-based VMAT, new treatment plans were optimized for MR-LINAC IMRT using 6X FFF fixed fields with maximum available field size of 27.4 x 24.1 cm2. Dose optimization was performed on the clinically used planning CT and structure set. Prescribed dose was 48.6 Gy in 27 fractions for all patients with 6 patients receiving an additional integrated boost for a total of 58.05 Gy to involved nodes. Constraints were based on our institutional protocol as per Table 1. IMRT delivery time was limited to 20 min. Original clinically used VMAT plans were generated on Eclipse (Varian Medical System) using 3 to 4 arcs. For the multi-institutional planning challenge, the data set from a single patient was anonymized and shared to participants. Participants used a single MR-based Linac planning platform to generate a plan based on our institutional constraints, with maximum treatment time limited to 20 min. For all analyses, a paired samples t-test was used to compare the significance defined at p < 0.05.

RESULTS

For MR-LINAC plans, the mean number of fields used was 23, mean number of segments 229, and the average estimated treatment delivery time was 17.3 minutes. MR-LINAC plans showed a significantly higher heterogeneity and dose to organs at risk compared to VMAT plans (Table 1). For the planning challenge, a total of 30 participants registered interest. Of this, seven plans were submitted to the challenge. On average, participants generated a plan that would be acceptable based on our institutional constraints (Table 1). However, the volumetric dose to bowel and pelvic bones were higher on MR-LINAC plans compared to the reference VMAT plan.

CONCLUSION

MR-LINAC fixed-field IMRT for cervix cancer patients is feasible but system constraints and optimization implementation result in greater dose heterogeneity and worse organ-at-risk sparing compared to Linac based VMAT. Further research is needed to determine if potential reduction of treatment margins, allowed by better MRI soft-tissue visualization, will result in MR-LINAC IMRT superior to Linac VMAT.

Analytical HDR prostate brachytherapy planning with automatic catheter and isotope selection

BACKGROUND

High dose rate (HDR) brachytherapy is commonly used to treat prostate cancer. Existing HDR planning systems solve the dwell time problem for predetermined catheters and a single energy source.

PURPOSE

Additional degrees of freedom can be obtained by relaxing the catheters' pre-designation and introducing more source types, and may have a dosimetric benefit, particularly in improving conformality to spare the urethra. This study presents a novel analytical approach to solving the corresponding HDR planning problem.

METHODS

The catheter and dual-energy source selection problem was formulated as a constrained optimization problem with a non-convex group sparsity regularization. The optimization problem was solved using the fast-iterative shrinkage-thresholding algorithm (FISTA). Two isotopes were considered. The dose rates for the HDR 4140 Ytterbium (Yb-169) source and the Elekta Iridium (Ir-192) HDR Flexisource were modeled according to the TG-43U1 formalism and benchmarked accordingly. Twenty-two retrospective HDR prostate brachytherapy patients treated with Ir-192 were considered. An Ir-192 only (IRO), Yb-169 only (YBO), and dual-source (DS) plan with optimized catheter location was created for each patient with N catheters, where N is the number of catheters used in the clinically delivered plans. The DS plans jointly optimized Yb-169 and Ir-192 dwell times. All plans and the clinical plans were normalized to deliver a 15 Gy prescription (Rx) dose to 95% of the clinical treatment volume (CTV) and evaluated for the CTV D90%, V150%, and V200%, urethra D0.1cc and D1cc, bladder V75%, and rectum V75%. Dose-volume histograms (DVHs) were generated for each structure.

RESULTS

The DS plans ubiquitously selected Ir-192 as the only treatment source. IRO outperformed YBO in organ at risk (OARs) OAR sparing, reducing the urethra D0.1cc and D1cc by 0.98% (p = 2.22 ∗ 10 - 9 $p\ = \ 2.22*{10^{ - 9}}$ ) and 1.09% (p = 1.22 ∗ 10 - 10 $p\ = \ 1.22*{10^{ - 10}}$ ) of the Rx dose, respectively, and reducing the bladder and rectum V75% by 0.09 (p = 0.0023 $p\ = \ 0.0023$ ) and 0.13 cubic centimeters (cc) (p = 0.033 $p\ = \ 0.033$ ), respectively. The YBO plans delivered a more homogenous dose to the CTV, with a smaller V150% and V200% by 3.20 (p = 4.67 ∗ 10 - 10 $p\ = \ 4.67*{10^{ - 10}}$ ) and 1.91 cc (p = 5.79 ∗ 10 - 10 $p\ = \ 5.79*{10^{ - 10}}$ ), respectively, and a lower CTV D90% by 0.49% (p = 0.0056 $p\ = \ 0.0056$ ) of the prescription dose. The IRO plans reduce the urethral D1cc by 2.82% (p = 1.38 ∗ 10 - 4 $p\ = \ 1.38*{10^{ - 4}}$ ) of the Rx dose compared to the clinical plans, at the cost of increased bladder and rectal V75% by 0.57 (p = 0.0022 $p\ = \ 0.0022$ ) and 0.21 cc (p = 0.019 $p\ = \ 0.019$ ), respectively, and increased CTV V150% by a mean of 1.46 cc (p = 0.010 $p\ = \ 0.010$ ) and CTV D90% by an average of 1.40% of the Rx dose (p = 8.80 ∗ 10 - 8 $p\ = \ 8.80*{10^{ - 8}}$ ). While these differences are statistically significant, the clinical differences between the plans are minimal.

CONCLUSIONS

The proposed analytical HDR planning algorithm integrates catheter and isotope selection with dwell time optimization for varying clinical goals, including urethra sparing. The planning method can guide HDR implants and identify promising isotopes for specific HDR clinical goals, such as target conformality or OAR sparing.

[Analysis on missed diagnosis or misdiagnosis of anomalous origin of left coronary artery from pulmonary artery by echocardiography from one single medical center]

Objectives: To analyze the reasons of missed diagnosis or misdiagnosis on anomalous origin of left coronary artery from pulmonary artery (ALCAPA) by echocardiography. Methods: This is a retrospective study. Patients with ALCAPA who underwent surgical treatment in Union Hospital, Tongji Medical College, Huazhong University of Science and Technology from August 2008 to December 2021 were included. According to the results of preoperative echocardiography and surgical diagnosis, the patients were divided into confirmed group or missed diagnosis/misdiagnosis group. The results of preoperative echocardiography were collected, and the specific echocardiographic signs were analyzed. According to the experience of the doctors, the echocardiographic signs were divided into four types, namely clear displayed, vague/doubtful displayed, no display and no notice, and the display rate of each sign was calculated (display rate=number of clearly displayed cases/total number of cases×100%). By referring the surgical data, we analyzed and recorded the pathological anatomy and pathophysiological characteristics of the patients, and the rate of missed diagnosis/misdiagnosis of echocardiography in patients with different characteristics was compared. Results: A total of 21 patients were enrolled, including 11 males, aged 1.8 (0.8, 12.3) years (range 1 month to 47 years). Except for one patient with anomalous origin of left anterior descending artery, the others were all originated from the main left coronary artery (LCA). There were 13 cases of ALCAPA in infant and children, and 8 cases of adult ALCAPA. There were 15 cases in the confirmed group (diagnostic accuracy was 71.4% (15/21)), and 6 cases in the missed diagnosis/misdiagnosis group (three cases were misdiagnosed as primary endocardial fibroelastosis, two cases were misdiagnosed as coronary-pulmonary artery fistula; and one case was missed diagnosis). The working years of the physicians in the confirmed group were longer than those in the missed diagnosis/misdiagnosed group ((12.8±5.6) years vs. (8.3±4.7) years, P=0.045). In infants with ALCAPA, the detection rate of LCA-pulmonary shunt (8/10 vs. 0, P=0.035) and coronary collateral circulation (7/10 vs. 0, P=0.042) in confirmed group was higher than that in missed diagnosis/misdiagnosed group. In adult ALCAPA patients, the detection rate of LCA-pulmonary artery shunt was higher in confirmed group than that in missed diagnosis/misdiagnosed group (4/5 vs. 0, P=0.021). The missed diagnosis/misdiagnosis rate of adult type was higher than that of infant type (3/8 vs. 3/13, P=0.410). The rate of missed diagnosis/misdiagnosis was higher in patients with abnormal origin of branches than that of abnormal origin of main trunk (1/1 vs. 5/21, P=0.028). The rate of missed diagnosis/misdiagnosis in patients with LCA running between the main and pulmonary arteries was higher than that distant from the main pulmonary artery septum (4/7 vs. 2/14, P=0.064). The rate of missed diagnosis/misdiagnosis in patients with severe pulmonary hypertension was higher than that in patients without severe pulmonary hypertension (2/3 vs. 4/18, P=0.184). The reasons with an echocardiography missed diagnosis/misdiagnosis rate of≥50% included that (1) the proximal segment of LCA ran between the main and pulmonary arteries; (2) abnormal opening of LCA at the right posterior part of the pulmonary artery; (3) abnormal origin of LCA branches; (4) complicated with severe pulmonary hypertension. Conclusions: Echocardiography physicians' knowledge of ALCAPA and diagnostic vigilance are critical to the accuracy of diagnosis. Attention should be paid to the pediatric cases with no obvious precipitating factors of left ventricular enlargement, regardless of whether the left ventricular function is normal or not, the origin of coronary artery should be routinely explored.

Deep Transfer Learning Approach for Localization of Damage Area in Composite Laminates Using Acoustic Emission Signal

Intelligent composite structures with self-aware functions are preferable for future aircrafts. The real-time location of damaged areas of composites is a key step. In this study, deep transfer learning was used to achieve the real-time location of damaged areas. The sensor network obtained acoustic emission signals from different damaged areas of the aluminum alloy plate. The acoustic emission time-domain signal is transformed into the input image by continuous wavelet transform. The convolutional neural network-based model automatically localized the damaged area by extracting features from the input image. A small amount of composite acoustic emission data was used to fine-tune some network parameters of the basic model through transfer learning. This enabled the model to classify the damaged area of composites. The accuracy of the transfer learning model trained with 900 samples is 96.38%, which is comparable to the accuracy of the model trained directly with 1800 samples; the training time of the former is only 17.68% of that of the latter. The proposed method can be easily adapted to new composite structures using transfer learning and a small dataset, providing a new idea for structural health monitoring.

Tomographic detection of photon pairs produced from high-energy X-rays for the monitoring of radiotherapy dosing

Measuring the radiation dose reaching a patient's body is difficult. Here we report a technique for the tomographic reconstruction of the location of photon pairs originating from the annihilation of positron-electron pairs produced by high-energy X-rays travelling through tissue. We used Monte Carlo simulations on pre-recorded data from tissue-mimicking phantoms and from a patient with a brain tumour to show the feasibility of this imaging modality, which we named 'pair-production tomography', for the monitoring of radiotherapy dosing. We simulated three image-reconstruction methods, one applicable to a pencil X-ray beam scanning through a region of interest, and two applicable to the excitation of tissue volumes via broad beams (with temporal resolution sufficient to identify coincident photon pairs via filtered back projection, or with higher temporal resolution sufficient for the estimation of a photon's time-of-flight). In addition to the monitoring of radiotherapy dosing, we show that image contrast resulting from pair-production tomography is highly proportional to the material's atomic number. The technique may thus also allow for element mapping and for soft-tissue differentiation.

An efficient rectangular optimization method for sparse orthogonal collimator based small animal irradiation

Objective.Intensity-modulated radiotherapy (IMRT) is widely used in clinical radiotherapy, treating varying malignancies with conformal doses. As the test field for clinical translation, preclinical small animal experiments need to mimic the human radiotherapy condition, including IMRT. However, small animal IMRT is a systematic challenge due to the lack of corresponding hardware and software for miniaturized targets.Approach.The sparse orthogonal collimators (SOC) based on the direct rectangular aperture optimization (RAO) substantially simplified the hardware for miniaturization. This study investigates and evaluates a significantly improved RAO algorithm for complex mouse irradiation using SOC. Because the Kronecker product representation of the rectangular aperture is the main limitation of the computational performance, we reformulated matrix multiplication in the data fidelity term using multiplication with small matrices instead of the Kronecker product of the dose loading matrices. Solving the optimization problem was further accelerated using the Fast Iterative Shrinkage-Thresholding Algorithm (FISTA).Main results.Four mouse cases, including a liver, a brain tumor, a concave U-target, and a complex total marrow irradiation (TMI) case, were included in this study with manually delineated targets and OARs. Seven coplanar-field SOC IMRT (sIMRT) plans were compared with idealistic fluence map based IMRT (iIMRT) plans. For the first three cases with simpler and smaller targets, the differences between sIMRT plans and iIMRT plans in the planning target volumes (PTV) statistics are within 1%. For the TMI case, the sIMRT plans are superior in reducing hot spots (also termedDmax) of PTV, kidneys, lungs, heart, and bowel by 20.5%, 31.5%, 24.67%, 20.13%, and 17.78%, respectively. On average, in four cases in this study, the sIMRT plan conformity is comparable to that of the iIMRT's with lightly increased R50 and Integral Dose by 2.23% and 2.78%.Significance.The significantly improved sIMRT optimization method allows fast plan creation in under 1 min for smaller targets and makes complex TMI planning feasible while achieving comparable dosimetry to idealistic IMRT with fluence map optimization.

[Effect analysis of treating intracranial wide-neck bifurcation aneurysms through Woven EndoBridge]

Objective: To explore the clinical effect of Woven EndoBridge (WEB) in the treatment of wide-neck bifurcation aneurysms. Methods: The clinical and imaging data of 11 patients with intracranial wide-neck bifurcation aneurysms treated by WEB alone at Department of Neurosurgery of the Northern Theater General Hospital from September 2017 to May 2018, were retrospectively analyzed. The patients were 7 males and 4 females, aged (54±11) years (ranged from 31 to 66 years). The aneurysms of 5 patients were located in the anterior communicating artery, 3 in the top of the basilar artery, and 3 in the bifurcation of the middle cerebral artery. The intraoperative and postoperative conditions of the patients were recorded, and the degree of aneurysm embolization was evaluated by WEB embolization aneurysm occlusion scale (WOS). Results: The intraoperative WEB release of all the 11 patients was good, with 3 cases of WOS grade A, 1 of grade B and 7 of grade C, with no intraoperative acute complications occurring. The imaging follow-up was not carried out in 1 patient due to economic reason, and the clinical follow-up was good until 3 years after the operation; 10 patients were followed up by imaging for 6 months to 3 years, and no postoperative complications occurred in the target treatment area. Among the 2 patients with WOS grade A and 1 patient with grade B during operation, according to the postoperative follow-up, all were WOS grade A; among the 7 patients with WOS grade C during operation, 4 were still of grade C and 3 were of grade D according to the follow-up. Among the 3 patients with WOS grade D, 1 patient received secondary embolization due to poor recurrence morphology, unstable hemodynamics and high possibility of rupture of aneurysm, stent assisted coil embolization was adopted, with good immediate effect; the other 2 cases had recurrent aneurysms, but the aneurysms had good morphology and stable hemodynamics, therefore, clinical follow-up was continued and no secondary surgery was performed. No complications occurred in all these 11 patients. Conclusions: The operation of treating unruptured intracranial wide-neck bifurcation aneurysms with WEB device alone is simple, and there is no need for anticoagulation and antiplatelet treatment before and after the operation, the clinical effect is being good. WEB device provides a new treatment option for intracranial wide-neck bifurcation aneurysms.

Reformulated McNamara RBE-weighted beam orientation optimization for intensity modulated proton therapy

PURPOSE

Empirical relative biological effectiveness (RBE) models have been used to estimate the biological dose in proton therapy but do not adequately capture the factors influencing RBE values for treatment planning. We reformulate the McNamara RBE model such that it can be added as a linear biological dose fidelity term within our previously developed sensitivity-regularized and heterogeneity-weighted beam orientation optimization (SHBOO) framework.

METHODS

Based on our SHBOO framework, we formulated the biological optimization problem to minimize total McNamara RBE dose to OARs. We solve this problem using two optimization algorithms: FISTA (McNam-FISTA) and Chambolle-Pock (McNam-CP). We compare their performances with a physical dose optimizer assuming RBE = 1.1 in all structures (PHYS-FISTA) and an LET-weighted dose model (LET-FISTA). Three head and neck patients were planned with the four techniques and compared on dosimetry and robustness.

RESULTS

Compared to Phys-FISTA, McNam-CP was able to match CTV [HI, Dmax, D95%, D98%] by [0.00, 0.05%, 1.4%, 0.8%]. McNam-FISTA and McNam-CP were able to significantly improve overall OAR [Dmean, Dmax] by an average of [36.1%,26.4%] and [29.6%, 20.3%], respectively. Regarding CTV robustness, worst [Dmax, V95%, D95%, D98%] improvement of [-6.6%, 6.2%, 6.0%, 4.8%] was reported for McNam-FISTA and [2.7%, 2.7%, 5.3%, -4.3%] for McNam-CP under combinations of range and setup uncertainties. For OARs, worst [Dmax, Dmean] were improved by McNam-FISTA and McNam-CP by an average of [25.0%, 19.2%] and [29.5%, 36.5%], respectively. McNam-FISTA considerably improved dosimetry and CTV robustness compared to LET-FISTA, which achieved better worst-case OAR doses.

CONCLUSION

The four optimization techniques deliver comparable biological doses for the head and neck cases. Besides modest CTV coverage and robustness improvement, OAR biological dose and robustness were substantially improved with both McNam-FISTA and McNam-CP, showing potential benefit for directly incorporating McNamara RBE in proton treatment planning.

[Comparison of single incision robot-assisted laparoscopic radical prostatectomy with and without extraperitoneal special channel device]

Objective: To compare the clinical effects of single-incision robot-assisted laparoscopic radical prostatectomy (RARP) with and without extraperitoneal special channel device. Methods: The clinical data of 70 patients who had undergone RARP in the Robotic Minimally Invasive Surgery Center of Sichuan Provincial People's Hospital from September 2020 to February 2021 were analyzed retrospectively, including 29 cases who were operated on without special channel device (group A) and 41 cases with special channel device (group B). All operations were performed by robot-assisted single-incision retrograde bladder neck exfoliation via extraperitoneal approach in patients by the same operator. The operation time, intraoperative blood loss, the bladder neck urethral anastomosis time, postoperative hospital stay, postoperative exhaust time, positive rate of incisal margin, indwelling time of urinary catheter, retention rate of postoperative erectile function, satisfaction rate of immediate postoperative urine control, positive rate of postoperative lymph node pathology, incision length, treatment cost and the rate of prostate specific antigen (PSA)lower than 0.2 μg/L at 6 weeks after operation were compared between the two groups. Results: All 70 cases were operated successfully. The difference of age[ (68.9±3.9) vs (69.4±5.4) years], preoperative PSA level[14.1(6.3, 19.8)vs13.7(5.8, 18.1)μg/L], prostate volume[44.8(30.7,172.6)vs 56.3(40.9,163.4)ml ] of the two groups was not statistically significant(all P>0.05). The difference of operation time [ (59.1±18.5) vs (59.6±18.0) min ], intraoperative blood loss [93(66,198)vs 95(68,203) ml ], bladder neck urethral anastomosis time [ (12.6±1.3) vs (13.7±2.8) min ], postoperative hospital stay [ (8.1±2.3) vs (9.1±1.3) d], postoperative exhaust time [ (1.4±0.6) vs (1.3±0.6) d], positive rate of incisal margin (20.7% vs 19.5%), indwelling time of the urinary catheter after operation [ (6.8±1.5) vs (7.1±2.0) d ], the retention rate of postoperative erectile function (31.0% vs 27.0%), the satisfaction rate of immediate postoperative urine control (79.3% vs 75.6%), the positive rate of postoperative lymph node pathology (17.2% vs 14.6%), the length of incision [ (5.1±0.5) vs (6.1±0.4) cm ], the rate of PSA lower than 0.2 μg/L at 6 weeks after operation (86.2% vs 83.0%) of the two groups was not statistically significant(all P>0.05). The operation cost of group A[(62 000±4 000) yuan]was lower than group B[(68 000±4 000) yuan] (P<0.05). Conclusion: Extraperitoneal non-special channel device single-incision RARP is safe and feasible.

ROAD: ROtational direct Aperture optimization with a Decoupled ring-collimator for FLASH radiotherapy

Ultra-high dose rate in radiotherapy (FLASH) has been shown to increase the therapeutic index with markedly reduced normal tissue toxicity and the same or better tumor cell killing. The challenge to achieve FLASH using x-rays, besides developing a high output linac, is to intensity-modulate the high-dose-rate x-rays so that the biological gain is not offset by the lack of physical dose conformity. In this study, we develop the ROtational direct Aperture optimization with a Decoupled ring-collimator (ROAD) to achieve simultaneous ultrafast delivery and complex dose modulation. The ROAD design includes a fast-rotating slip-ring linac and a decoupled collimator-ring with 75 pre-shaped multi-leaf-collimator (MLC) modules. The ring-source rotates at 1 rotation per second (rps) clockwise while the ring-collimator is either static or rotating at 1 rps counterclockwise, achieving 75 (ROAD-75) or 150 (ROAD-150) equal-angular beams for one full arc. The Direct Aperture Optimization (DAO) for ROAD was formulated to include a least-square dose fidelity, an anisotropic total variation term, and a single segment term. The FLASH dose (FD) and FLASH biological equivalent dose (FBED) were computed voxelwise, with the latter using a spatiotemporal model accounting for radiolytic oxygen depletion. ROAD was compared with clinical volumetric modulated arc therapy (VMAT) on a brain, a lung, a prostate, and a head and neck cancer patient. The mean dose rate of ROAD-75 and ROAD-150 are 76.2 Gy s^-1 and 112 Gy s^-1 respectively to deliver 25 Gy single-fraction dose in 1 s. With improved PTV homogeneity, ROAD-150 reduced (max, mean) OAR physical dose by (4.8 Gy, 6.3 Gy). The average R50 and integral dose of (VMAT, ROAD-75, ROAD-150) are (4.8, 3.2, 3.2) and (89, 57, 56) Gy×Liter, respectively. The FD and FBED showed model dependent FLASH effects. The novel ROAD design achieves ultrafast dose delivery and improves physical dosimetry compared with clinical VMAT, providing a potentially viable engineering solution for x-ray FLASH radiotherapy.

DeepMC: a deep learning method for efficient Monte Carlo beamlet dose calculation by predictive denoising in magnetic resonance-guided radiotherapy

Emerging magnetic resonance (MR) guided radiotherapy affords significantly improved anatomy visualization and, subsequently, more effective personalized treatment. The new therapy paradigm imposes significant demands on radiation dose calculation quality and speed, creating an unmet need for the acceleration of Monte Carlo (MC) dose calculation. Existing deep learning approaches to denoise the final plan MC dose fail to achieve the accuracy and speed requirements of large-scale beamlet dose calculation in the presence of a strong magnetic field for online adaptive radiotherapy planning. Our deep learning dose calculation method, DeepMC, addresses these needs by predicting low-noise dose from extremely noisy (but fast) MC-simulated dose and anatomical inputs, thus enabling significant acceleration. DeepMC simultaneously reduces MC sampling noise and predicts corrupted dose buildup at tissue-air material interfaces resulting from MR-field induced electron return effects. Here we demonstrate our model's ability to accelerate dose calculation for daily treatment planning by a factor of 38 over traditional low-noise MC simulation with clinically meaningful accuracy in deliverable dose and treatment delivery parameters. As a post-processing approach, DeepMC provides compounded acceleration of large-scale dose calculation when used alongside established MC acceleration techniques in variance reduction and graphics processing unit-based MC simulation.

A novel energy layer optimization framework for spot-scanning proton arc therapy

PURPOSE

Spot-scanning proton arc therapy (SPAT) is an emerging modality to improve plan conformality and delivery efficiency. A greedy and heuristic method is proposed in the existing SPAT algorithm to select energy layers and sequence energy switching with gantry rotation, which does not promise optimality in either dosimetry or efficiency. We aim to develop a method to solve the energy layer switching and dosimetry optimization problems in an integrated framework for SPAT.

METHODS

In an integrated approach, energy layer optimization for spot-scanning proton arc therapy (ELO-SPAT) is formulated with a dose fidelity term, a group sparsity regularization, a log barrier regularization, and an energy sequencing (ES) penalty. The combination of L2,1/2-norm group sparsity regularization and log barrier function allows one energy layer being selected per control point. The ES regularization term sorts the delivery sequence from high energy to low energy to reduce the total energy layer switching time (ELST) and subsequently the total delivery time. Within the ES penalty, the gradient of layer weights between adjacent beams is first calculated along beam direction and then along energy direction. The gradients indicate energy switch patterns between two adjacent beams. The time-wise costly energy switch-up is more heavily penalized in the ES term. This ELO-SPAT method was tested on one frontal base-of-skull (BOS) patient, one chordoma (CHDM) patient with a simultaneous integrated boost, one bilateral head-and-neck (H&N) patient, and one lung (LNG) patient. We compared ELO-SPAT with intensity-modulated proton therapy (IMPT) using discrete beams and SPArc by Ding et al. For the two arc algorithms, both the plans with and without energy sequencing were created and compared.

RESULTS

Energy layer optimization for spot-scanning proton arc therapy reduced the runtime of optimization by 84% on average compared with the greedy SPArc method. In both the ELO-SPAT plans with and without ES, one energy layer per control point was selected. Without ES regularization, the energy sequence was arbitrary, with around 40-60 switch-up for the tested cases. After adding ES regularization, the number of energy switch-up was reduced to less than 20. Compared with the energy sequenced SPArc plans, the ELO-SPAT plans with ES led to 24% less total ELST for synchrotron plans and 14% less for cyclotron plans. Both the ELO-SPAT and SPArc plans achieved better sparing compared with the IMPT plans for most Organs-at-risks (OARs), with or without ES. Without ES, the ELO-SPAT plans achieved further improvement of the OARs compared with the SPArc plans, with an averaged reduction of OAR [Dmean, Dmax] by [1.57, 3.34] GyRBE. Adding the ES regularization degraded the plan quality, but the ELO-SPAT plans still had comparable or slightly better sparing than the SPArc plans with ES, with an averaged reduction of OAR [Dmean, Dmax] by [1.42, 2.34] GyRBE.

CONCLUSION

We developed a computationally efficient spot-scanning proton arc optimization method, which solved energy layer selection and sequencing in an integrated framework, generating plans with good dosimetry and high delivery efficiency.

Many-isocenter optimization for robotic radiotherapy

Despite significant dosimetric gains, clinical implementation of the 4π non-coplanar radiotherapy on the widely available C-arm gantry system is hindered by limited clearance, and the need to perform complex coordinated gantry and couch motion. A robotic radiotherapy platform would be conducive to such treatment but a new conflict between field size and MLC modulation resolution needs to be managed for versatile applications. This study investigates the dosimetry and delivery efficiency of purposefully creating many isocenters to achieve simultaneously high MLC modulation resolution and large tumor coverage. An integrated optimization framework was proposed for simultaneous beam orientation optimization (BOO), isocenter selection, and fluence map optimization (FMO). The framework includes a least-square dose fidelity objective, a total variation term for regularizing the fluence smoothness, and a group sparsity term for beam selection. A minimal number of isocenters were identified for efficient target coverage. Colliding beams excluded, high-resolution small-field 4π intensity-modulated radiotherapy (IMRT) treatment plans with 50 cm source-to-isocenter distance (SID-50) on 10 Head and Neck (H&N) cancer patients were compared with low-resolution large-field plans with 100 cm SID (SID-100). With the same or better target coverage, the average reduction of [Dmean, Dmax] of 20-beam SID-50 plans from 20-beam SID-100 plans were [2.09 Gy, 1.19 Gy] for organs at risk (OARs) overall, [3.05 Gy, 0.04 Gy] for parotid gland, [3.62 Gy, 5.19 Gy] for larynx, and [3.27 Gy, 1.10 Gy] for mandible. R50 and integral dose were reduced by 5.3% and 9.6%, respectively. Wilcoxon signed-rank test showed significant difference (p   <  0.05) in planning target volume (PTV) homogeneity, PTV Dmax, R50, Integral dose, and OAR Dmean and Dmax. The estimated delivery time of 20-beam [SID-50, SID-100] plans were [19, 18] min and [14, 9] min, assuming 5 fractions and 30 fractions, respectively. With clinically acceptable delivery efficiency, many-isocenter optimization is dosimetrically desirable for treating large targets with high modulation resolution on the robotic platform.

Image-domain multimaterial decomposition for dual-energy computed tomography with nonconvex sparsity regularization

Dual-energy computed tomography (CT) has the potential to decompose tissues into different materials. However, the classic direct inversion (DI) method for multimaterial decomposition (MMD) cannot accurately separate more than two basis materials due to the ill-posed problem and amplified image noise. We propose an integrated MMD method that addresses the piecewise smoothness and intrinsic sparsity property of the decomposition image. The proposed MMD was formulated as an optimization problem including a quadratic data fidelity term, an isotropic total variation term that encourages image smoothness, and a nonconvex penalty function that promotes decomposition image sparseness. The mass and volume conservation rule was formulated as the probability simplex constraint. An accelerated primal-dual splitting approach with line search was applied to solve the optimization problem. The proposed method with different penalty functions was compared against DI on a digital phantom, a Catphan^® 600 phantom, a quantitative imaging phantom, and a pelvis patient. The proposed framework distinctly separated the CT image up to 12 basis materials plus air with high decomposition accuracy. The cross talks between two different materials are substantially reduced, as shown by the decreased nondiagonal elements of the normalized cross correlation (NCC) matrix. The mean square error of the measured electron densities was reduced by 72.6%. Across all datasets, the proposed method improved the average volume fraction accuracy from 61.2% to 99.9% and increased the diagonality of the NCC matrix from 0.73 to 0.96. Compared with DI, the proposed MMD framework improved decomposition accuracy and material separation.

[Streptococcal toxic shock syndrome caused by Streptococcus pyogenes: a retrospective study of 15 pediatric cases]

Objective: To improve the understanding of clinical characteristics of streptococcal toxic shock syndrome (STSS) caused by Streptococcus pyogenes (S. pyogenes) in children. Methods: A retrospective study was conducted to analyze the clinical data of STSS caused by S. pyogenes (culture-confirmed) in 7 tertiary hospitals during 2010-2017 in China. Clinical and laboratory data were collected by reviewing the medical records. Results: Fifteen cases of STSS, including 9 males, were confirmed and the ages of the patients ranged from 6 months to 15 years, with median age of 3 years. All cases had the positive blood culture for S. pyogenes and only 3 cases had short course of β-lactam treatment before blood culture. Medical evaluation was initiated within (5.1±4.6) days after symptom onset. All patients had fever, and 13 patients had multiple organ dysfunction and 10 patients had disseminated intravascular coagulationl (DIC). Twelve cases had severe pneumonia with or without skin and (or) soft tissue infections. Underlying conditions included giant hemangioma of the skin in 2 patients and varicella in 1 patient. All isolated strains in 14 cases were sensitive to penicillin G, ceftriaxone/cefotaxime, vancomycin, but 12 and 13 isolates were resistant to clindamycin and erythromycin, respectively. Eight patients died, and 5 of them died within 24 hours after admission. One patient was lost to follow-up after intended discharge against medical advice. Conclusion: STSS caused by S. pyogenes in children is a severe syndrome with rapid clinical progression and high mortality rate, and thus the pediatricians should be aware of STSS and immediately initiate aggressive treatment for the suspected cases.

[Association of CYP19A1 gene rs7176005 single nucleotide polymorphism with breast cancer risk and clinicopathologic features of tumor]

Objective: The aim of this study was to investigate the association of the CYP19A1 rs7176005 single nucleotide polymorphism (SNP) with breast cancer risk and with clinicopathologic features of tumors. Methods: This study was conducted by including 138 patients with breast cancer (cancer group), those who diagnosed as primary breast cancer after operation by pathology. There were 293 cases in the group of benign breast disease which was presented as a solid mass by the color ultrasound and pathologically diagnosed as "fibroadenoma or adenosis" (benign breast disease group), the cases were paired with breast cancer patients by age±5 in the same period, and there were 259 cases in the group of healthy control who received routine physical examination during the same period and were paired with breast cancer patients by age±5 without any detection of breast related diseases (healthy control group) at West China hospital between September 2012 and November 2016. The CYP19A1 rs7176005 SNP was detected by a direct sequencing method. Hardy-Weinberg test was used to analyze the genetic balance of the 3 groups. Chi square test was used to compare the distribution of rs7176005 genotypes between the 3 groups, and the differences of clinicopathological features in breast cancer patients carrying different genotypes. Results: The ages of the breast cancer cases, the benign breast disease group and the healthy control group were (44.69±8.09), (42.33±11.44) and (41.92±9.61) years old, respectively. Hardy-Weinberg equilibrium test identified that the composition ratios of alleles C and T in breast cancer group, benign breast disease group and healthy group were not statistically significant (χ(2) values were 0.83, 0.34 and 0.04, respectively, P values were 0.363, 0.561, and 0.852, respectively). All the three groups met the genetic balance, had consistency and could represent the population. Among the 138 cases of breast cancer, the CYP19A1 rs7176005 SNP was significantly associated with the diameter of the tumor (P=0.031). The majority of tumor size was <2 cm in patients who carrying TT and CT genotypes, and the proportion was 75% (12/16) and 58% (40/69), respectively. While those patients with TT genotype were mainly >2 cm and ≤5 cm, and the proportion was 51% (27/53). The distribution of TNM stage among patients with different genotypes was also statistically significant (χ(2)=11.19, P=0.025). The most common stage was Ⅱ in Patients who carrying CC and CT genotypes, and the proportion was 45.3% (24/53) and 52.2% (36/69), respectively. While those patients with TT genotype was mainly in stage Ⅰ and the proportion was 56.3% (9/16). Conclusion: Though the CYP19A1 rs7176005 SNP is not associated with breast cancer development, breast cancer patients with the C allele exhibit a high tumor growth rate and large diameters.

A novel optimization framework for VMAT with dynamic gantry couch rotation

Existing volumetric modulated arc therapy (VMAT) optimization using coplanar arcs is highly efficient but usually dosimetrically inferior to intensity modulated radiation therapy (IMRT) with optimized non-coplanar beams. To achieve both dosimetric quality and delivery efficiency, we proposed in this study, a novel integrated optimization method for non-coplanar VMAT (4πVMAT). 4πVMAT with direct aperture optimization (DAO) was achieved by utilizing a least square dose fidelity objective, along with an anisotropic total variation term for regularizing the fluence smoothness, a single segment term for imposing simple apertures, and a group sparsity term for selecting beam angles. Continuous gantry/couch angle trajectories were selected using the Dijkstra's algorithm, where the edge and node costs were determined based on the maximal gantry rotation speed and the estimated fluence map at the current iteration, respectively. The couch-gantry-patient collision space was calculated based on actual machine geometry and a human subject 3D surface. Beams leading to collision are excluded from the DAO and beam trajectory selection (BTS). An alternating optimization strategy was implemented to solve the integrated DAO and BTS problem. The feasibility of 4πVMAT using one full-arc or two full-arcs was tested on nine patients with brain, lung, or prostate cancer. The plan was compared against a coplanar VMAT (2πVMAT) plan using one additional arc and collimator rotation. Compared to 2πVMAT, 4πVMAT reduced the average maximum and mean organs-at-risk dose by 9.63% and 3.08% of the prescription dose with the same target coverage. R50 was reduced by 23.0%. Maximum doses to the dose limiting organs, such as the brainstem, the major vessels, and the proximal bronchus, were reduced by 8.1 Gy (64.8%), 16.3 Gy (41.5%), and 19.83 Gy (55.5%), respectively. The novel 4πVMAT approach affords efficient delivery of non-coplanar arc trajectories that lead to dosimetric improvements compared with coplanar VMAT using more arcs.

VMAT optimization with dynamic collimator rotation

PURPOSE

Although collimator rotation is an optimization variable that can be exploited for dosimetric advantages, existing Volumetric Modulated Arc Therapy (VMAT) optimization uses a fixed collimator angle in each arc and only rotates the collimator between arcs. In this study, we develop a novel integrated optimization method for VMAT, accounting for dynamic collimator angles during the arc motion.

METHODS

Direct Aperture Optimization (DAO) for Dynamic Collimator in VMAT (DC-VMAT) was achieved by adding to the existing dose fidelity objective an anisotropic total variation term for regulating the fluence smoothness, a binary variable for forming simple apertures, and a group sparsity term for controlling collimator rotation. The optimal collimator angle for each beam angle was selected using the Dijkstra's algorithm, where the node costs depend on the estimated fluence map at the current iteration and the edge costs account for the mechanical constraints of multi-leaf collimator (MLC). An alternating optimization strategy was implemented to solve the DAO and collimator angle selection (CAS). Feasibility of DC-VMAT using one full-arc with dynamic collimator rotation was tested on a phantom with two small spherical targets, a brain, a lung and a prostate cancer patient. The plan was compared against a static collimator VMAT (SC-VMAT) plan using three full arcs with 60 degrees of collimator angle separation in patient studies.

RESULTS

With the same target coverage, DC-VMAT achieved 20.3% reduction of R50 in the phantom study, and reduced the average max and mean OAR dose by 4.49% and 2.53% of the prescription dose in patient studies, as compared with SC-VMAT. The collimator rotation co-ordinated with the gantry rotation in DC-VMAT plans for deliverability. There were 13 beam angles in the single-arc DC-VMAT plan in patient studies that requires slower gantry rotation to accommodate multiple collimator angles.

CONCLUSIONS

The novel DC-VMAT approach utilizes the dynamic collimator rotation during arc delivery. In doing so, DC-VMAT affords more sophisticated intensity modulation, alleviating the limitation previously imposed by the square beamlet from the MLC leaf thickness and achieves higher effective modulation resolution. Consequently, DC-VMAT with a single arc manages to achieve superior dosimetry than SC-VMAT with three full arcs.

Technical Note: Iterative megavoltage CT (MVCT) reconstruction using block-matching 3D-transform (BM3D) regularization

PURPOSE

Megavoltage CT (MVCT) images are noisier than kilovoltage CT (KVCT) due to low detector efficiency to high-energy x rays. Conventional denoising methods compromise edge resolution and low-contrast object visibility. In this work, we incorporated block-matching 3D-transform shrinkage (BM3D) transformation into MVCT iterative reconstruction as nonlocal patch-wise regularization.

METHODS

The iterative reconstruction was achieved by adding to the existing least square data fidelity objective a regularization term, formulated as the L1 norm of the BM3D transformed image. A Fast Iterative Shrinkage-Thresholding Algorithm (FISTA) was adopted to accelerate CT reconstruction. The proposed method was compared against total variation (TV) regularization, BM3D postprocess method, and filtered back projection (FBP).

RESULTS

In the Catphan phantom study, BM3D regularization better enhances low-contrast objects compared with TV regularization and BM3D postprocess method at the same noise level. The spatial resolution using BM3D regularization is 2.79 and 2.55 times higher than that using the TV regularization at 50% of the modulation transfer function (MTF) magnitude, for the fully sampled reconstruction and down-sampled reconstruction, respectively. The BM3D regularization images show better bony details and low-contrast soft tissues, on the head and neck (H&N) and prostate patient images.

CONCLUSIONS

The proposed iterative BM3D regularization CT reconstruction method takes advantage of both the BM3D denoising capability and iterative reconstruction data fidelity consistency. This novel approach is superior to TV regularized iterative reconstruction or BM3D postprocess for improving noisy MVCT image quality.

Assessment of dynamic material properties of intact rocks using seismic wave attenuation: an experimental study

The mechanical properties of any substance are essential facts to understand its behaviour and make the maximum use of the particular substance. Rocks are indeed an important substance, as they are of significant use in the energy industry, specifically for fossil fuels and geothermal energy. Attenuation of seismic waves is a non-destructive technique to investigate mechanical properties of reservoir rocks under different conditions. The attenuation characteristics of five different rock types, siltstone, shale, Australian sandstone, Indian sandstone and granite, were investigated in the laboratory using ultrasonic and acoustic emission instruments in a frequency range of 0.1-1 MHz. The pulse transmission technique and spectral ratios were used to calculate the attenuation coefficient (α) and quality factor (Q) values for the five selected rock types for both primary (P) and secondary (S) waves, relative to the reference steel sample. For all the rock types, the attenuation coefficient was linearly proportional to the frequency of both the P and S waves. Interestingly, the attenuation coefficient of granite is more than 22% higher than that of siltstone, sandstone and shale for both P and S waves. The P and S wave velocities were calculated based on their recorded travel time, and these velocities were then used to calculate the dynamic mechanical properties including elastic modulus (E), bulk modulus (K), shear modulus (µ) and Poisson's ratio (ν). The P and S wave velocities for the selected rock types varied in the ranges of 2.43-4.61 km s^-1 and 1.43-2.41 km h^-1, respectively. Furthermore, it was observed that the P wave velocity was always greater than the S wave velocity, and this confirmed the first arrival of P waves to the sensor. According to the experimental results, the dynamic E value is generally higher than the static E value obtained by unconfined compressive strength tests.

A comprehensive formulation for volumetric modulated arc therapy planning

PURPOSE

Volumetric modulated arc therapy (VMAT) is a widely employed radiation therapy technique, showing comparable dosimetry to static beam intensity modulated radiation therapy (IMRT) with reduced monitor units and treatment time. However, the current VMAT optimization has various greedy heuristics employed for an empirical solution, which jeopardizes plan consistency and quality. The authors introduce a novel direct aperture optimization method for VMAT to overcome these limitations.

METHODS

The comprehensive VMAT (comVMAT) planning was formulated as an optimization problem with an L2-norm fidelity term to penalize the difference between the optimized dose and the prescribed dose, as well as an anisotropic total variation term to promote piecewise continuity in the fluence maps, preparing it for direct aperture optimization. A level set function was used to describe the aperture shapes and the difference between aperture shapes at adjacent angles was penalized to control MLC motion range. A proximal-class optimization solver was adopted to solve the large scale optimization problem, and an alternating optimization strategy was implemented to solve the fluence intensity and aperture shapes simultaneously. Single arc comVMAT plans, utilizing 180 beams with 2° angular resolution, were generated for a glioblastoma multiforme case, a lung (LNG) case, and two head and neck cases-one with three PTVs (H&N3PTV) and one with foue PTVs (H&N4PTV)-to test the efficacy. The plans were optimized using an alternating optimization strategy. The plans were compared against the clinical VMAT (clnVMAT) plans utilizing two overlapping coplanar arcs for treatment.

RESULTS

The optimization of the comVMAT plans had converged within 600 iterations of the block minimization algorithm. comVMAT plans were able to consistently reduce the dose to all organs-at-risk (OARs) as compared to the clnVMAT plans. On average, comVMAT plans reduced the max and mean OAR dose by 6.59% and 7.45%, respectively, of the prescription dose. Reductions in max dose and mean dose were as high as 14.5 Gy in the LNG case and 15.3 Gy in the H&N3PTV case. PTV coverages measured by D95, D98, and D99 were within 0.25% of the prescription dose. By comprehensively optimizing all beams, the comVMAT optimizer gained the freedom to allow some selected beams to deliver higher intensities, yielding a dose distribution that resembles a static beam IMRT plan with beam orientation optimization.

CONCLUSIONS

The novel nongreedy VMAT approach simultaneously optimizes all beams in an arc and then directly generates deliverable apertures. The single arc VMAT approach thus fully utilizes the digital Linac's capability in dose rate and gantry rotation speed modulation. In practice, the new single VMAT algorithm generates plans superior to existing VMAT algorithms utilizing two arcs.

[Inhibitory effect of migration-inducing gene-7-shRNA recombinant retrovirus combined with endostatin on growth and metastasis of hepatoma xenograft]

Objective: To investigate the inhibitory effect of migration-inducing gene-7(Mig-7)interfered with retrovirus-mediated RNA(shRNA)combined with recombinant human endostatin(ES)on the growth and metastasis of subcutaneous xenograft of human hepatoma cells in nude mice. Methods: Two Mig-7-mRNA oligonucleotide sequences(Mig-7-shRNA-1 and Mig-7-shRNA-2)and one sequence as a negative control(Mig-7-shRNA-N)were designed. The specific Mig-7-shRNA recombinant retrovirus expression vector plasmid was constructed and used for the transfection of human hepatoma MHCC-97H cells with high expression of Mig-7. The subcutaneous xenograft tumor model of human hepatocellular carcinoma(HCC)in nude mice was established, and according to the condition of transfection and administration, the nude mice were divided into pSIREN-M1 group, pSIREN-MN group, ES group, and pSIREN-M1+ES group. The xenograft tumor volume, mass, and metastasis were compared between groups. Immunohistochemistry was used to observe the formation of vasculogenic mimicry(VM)in xenograft tumor and the difference in tumor microvascular density(MVD), and Western blot was used to measure the expression of Mig-7 and vascular endothelial growth factor(VEGF)in each group. A one-way analysis of variance was used for comparison between groups, and the Fisher's exact test was used for comparison of continuous data between groups. Results: Compared with the pSIREN-MN group, the pSIREN-M1 group had significantly lower xenograft tumor volume, mass, and metastasis rate, Mig-7 expression, and formation of VM(P < 0.05), as well as significantly higher VEGF expression and MVD(P < 0.05). Compared with the pSIREN-MN group, the ES group had significantly lower xenograft tumor volume, mass, and metastasis rate, VEGF expression, and MVD(P < 0.05), as well as significantly higher Mig-7 expression and formation of VM(P < 0.05). Compared with the pSIREN-M1 group and the ES group, the pSIREN-M1+ES group had significantly lower xenograft tumor volume, mass, and metastasis rate, Mig-7 expression, formation of VM, VEGF expression, and MVD(P < 0.05). Conclusion: Mig-7-shRNA recombinant retrovirus combined with ES has a better inhibitory effect on the growth and metastasis of HCC xenograft tumor than Mig-7-shRNA recombinant retrovirus or ES alone. The anti-tumor angiogenesis therapy alone, which targets vascular endothelial cells in vivo, has a limited effect, since it may promote the formation of VM.

[Echocardiographic features of sinus of Valsalva aneurysm extending into left ventricle]

OBJECTIVE

To observe image features of sinus of Valsalva aneurysm (SVA) extending into left ventricle by echocardiography.

METHODS

Echocardiographic features of 5 cases of SVA extending into left ventricle and underwent surgery at Union Hospital from July 1995 to September 2015 were reviewed retrospectively and image and surgical findings were compared. A total of 4 patients were diagnosed by conventional and 3D echocardiography before surgery, and 1 patient was diagnosed by conventional echocardiography.

RESULTS

The origin, extending position, rupture status, complications of the SVA and associated cardiovascular lesions determined by echocardiography were entirely consistent with surgical findings in all cases, with the exception of one failed diagnosis of hypoplasty of an adjacent aortic cusp. Besides the common features of SVA, echocardiographic features of SVA extending into left ventricle included a thin-walled saccular lesion arising from the aortic root in continuation with the aortic annulus, with significant morphological changes and movement. Moreover, following features were observed: very low origin of the saccular lesion arising between the sinus base and the aortic annulus; the aneurysm going back and forth between the aortic root and the left ventricular outflow tract in 4 cases with an intact interventricular septum, and between the left ventricle and the right ventricle through the septal defect in another case complicated by a huge ventricular septal defect; diastolic shunt into the left ventricle when ruptured; displacement of the adjacent aortic annulus due to compression of aneurismal origin and prolapse of both aortic valve and annulus observed in all cases, resulting in a severe aortic regurgitation; obstructions of the left ventricular outflow tract due to the space-occupying effect were also found in 2 cases.

CONCLUSION

The SVA extending into left ventricle has distinguished echocardiographic characteristics, which could be accurately diagnosed either by conventional or real-time 3D echocardiography.