Cardiac Sparing with Volumetric Modulated Arc Therapy Enabled Total Body Irradiation (CS VMAT-TBI)

PURPOSE/OBJECTIVE(S)

Volumetric modulated arc therapy (VMAT) enabled total body irradiation (TBI) has replaced conventional TBI in our institution given the improved treatment accuracy, patient comfort, and dose modulation ability. The risk of cardiovascular disease is several folds higher among transplant patients who receive TBI, likely related to dose to the heart. We hypothesize that a cardiac-sparing (CS) VMAT-TBI technique is feasible and can meaningfully reduce dose to the heart while still adequately covering nearby lymphatic tissue.

MATERIALS/METHODS

VMAT-TBI is delivered via multi-isocentric external beams in a frame-based setup. Heart is contoured as per published guidelines. A lymph node contour, which includes tonsils, neck nodal stations, mediastinal, abdominal, retroperitoneal, and pelvic nodes is created. Coverage of the lymph node contour is prioritized over organ-sparing during inverse optimization; with a goal of V90% greater than 99.5% and mean dose less than 800 cGy for the lymph nodes and heart, respectively. An IRB-approved retrospective review was performed with mean heart dose collected for all patients treated with CS VMAT-TBI and compared to a representative cohort of five patients treated with VMAT-TBI without cardiac sparing.

RESULTS

Thirty-one patients were treated with CS VMAT-TBI between 2020-2022 with a median follow up time of 11.5 months. Mean heart dose was 796 ± 71 cGy in the CS VMAT-TBI compared to 1247 ± 29 cGy in the VMAT-TBI group without cardiac sparing (p < 0.001). Of those treated with CS VMAT-TBI, three patients relapsed; one relapse occurred in bone marrow only, one relapse occurred in bone marrow and cervical, thoracic, and intra-abdominal lymphoid tissues, and one patient was simulated but never received induction therapy due to overt progression. 100-day relapse-free survival and overall survival were 82.5% and 86.2%, respectively. Median survival time has not been met.

CONCLUSION

Cardiac sparing is feasible in VMAT-TBI and is associated with significant decrease in mean heart dose of ∼450 cGy. This is estimated to confer a 33.3% decreased absolute risk for lifetime major coronary events compared to patients treated with VMAT-TBI without cardiac sparing. Although limited by short follow-up time, there does not appear to be a significant risk for early relapse despite de-escalating cardiac tissue, likely due to prioritizing coverage of lymph nodes. Prospective clinical studies are needed to further validate cardiac and other organ at risk sparing VMAT-TBI techniques.

Towards Biology-Guided Radiotherapy Planning and Delivery on a Novel O-Ring PET-Linac Platform: Extended Beyond Bone and Lung Lesions

PURPOSE/OBJECTIVE(S)

Biology-guided radiotherapy (BgRT) with FDG signal collected via an on-board positron emission tomography (PET) system integrated in an O-ring gantry Linac was recently cleared by the FDA for lung and bone lesions. This study aims to determine if BgRT plans, guided via PET signal, are clinically acceptable for FDG-avid lesions in disease sites beyond bone and lung.

MATERIALS/METHODS

Ten patients previously treated for lesions in the liver, head and neck (HN), pancreas, renal and pelvic-abdominal lymph nodes were identified. Diagnostic PET/CT images of these treatment sites were first collected and processed/converted to mimic PET images that are acquired on PET-Linac and would be used to guide the delivery. For BgRT planning, the PTV was generated with 5 mm margin from GTV and a Biology Tracking Zone was generated including the anticipated full range of target motion. BgRT plans, guided by the emulated PET signal, were generated with 46Gy in 3 fractions for liver and 40Gy in 5 fractions for all other sites. BgRT plan deliverability was first assessed by evaluating the Activity Concentration (AC) and Normalized Target Signals (NTS) on converted PET images with the goal to meet NTS >2 (hard constraint) and AC >5kBq/ml (goal). BgRT plan quality was then evaluated with institutional guidelines on PTV coverage, OAR doses, conformity index (CI) and Heterogeneity index (HI).

RESULTS

BgRT plans were successfully generated for 11 target lesions among ten patients. The average diagnostic PET SUV, derived NTS and AC on converted PET images were 12.62, 9.33 and 12.10 kBq/ml, respectively. All images met the NTS constraints, and 8/11 plans met the AC goal for deliverability. All plans met the OAR hard constraints such as max dose on duodenum, small bowel, large bowel and spinal cord. Five of 11 plans had a limiting GI structure that resulted in an expected reduction in PTV coverage with an average PTV V100% = 77.9%, CI of 1.4, HI of 1.36 and max dose of 133.8%. The other 6 of 11 cases met the PTV V100% = 95%, had an average CI of 1.1, HI of 1.28 and Dmax of 127.67%. The estimated average time for BgRT delivery was 17 mins 25 secs. Although these plan parameters are deemed to be clinically acceptable, heterogeneity was detected inside the target region and suboptimal dose fall off was observed in some cases that may be caused by current implementation.

CONCLUSION

This preliminary study showed that BgRT plans were generated successfully with emulated PET images on 11 treatment sites covering HN, abdominal and pelvic regions. All plans met NTS constraints and 8 out of 11 met AC goals for deliverability. The plan quality of all BgRT plans were clinically acceptable based on institutional constraints. Further investigations are required to test more patients/sites for BgRT plan feasibility. Dosimetric benefit from margin reduction of BgRT target should also be investigated in future study.

Simulation-Omitted Replan with Cone Beam Computed Tomography based Adaptive Online Radiotherapy System - Transferring Adapted Plan to Non-Adaptive Ring Gantry Linear Accelerator for Image Guided Radiotherapy

PURPOSE/OBJECTIVE(S)

Artificial intelligence powered cone beam computed tomography (CBCT) based online adaptive radiotherapy (oART) system offers a streamlined and efficient process for daily ART as the default. In our prior work, we developed a workflow to utilize this oART system as simulation-omitted replan platform and treat the adapted plan on the oART system with image guided radiotherapy (IGRT) until next adaptation. However, the IGRT fractions will occupy the treatment slots of the machine. In this work, we aim to develop a semi-automatic workflow to allow the adapted plan to be treat on the non-adaptive ring gantry linear accelerator (non-ART Linac) and dedicate the oART system for adaptive treatments.

MATERIALS/METHODS

The oART system and the non-ART Linac were machine-matched to the same representative beam data. In the oART system, the initial plan is setup as 'adaptive' treatment and patients are only treated on the oART system for adaptive replan. The IGRT fractions are all treated on the non-ART Linac. An API script was developed to automatically (1) grab the adapted DICOM plan files from the secondary calculation system and write directly back to the database of the treatment management system (TMS), (2) change the DICOM tags to make the files compatible in the TMS system, (3) insert the kV-CBCT field to make the plan deliverable in the non-ART Linac. There are minimum remaining manual steps to setup the number of fractions to the intended number of IGRT fractions and to link plan to the prescription in TMS. We compare the required resources and the percentage of ART treatments on the oART system before and after the implementation of the proposed workflow to quantify the improvement of service.

RESULTS

The proposed workflow and automation eliminated the need to convert between IGRT/ART fractions in the Ethos system and reduced manual work by 25 minutes each adapted plan transfer. Table 1 summarizes the number of physics tasks and the percentage of ART fractions in oART system per month before and after the proposed workflow. This workflow reduced the physics IGRT/ART tasks from 107±31 to 65±21 tasks per month (p<0.05). Percentage of ART treatments on oART system increase from 30%±3% to 57%±13% (p<0.05). We also observed increased utilization of ART from 46% in the 1st month to 71% in the 6th month since it is easier to find a feasible time slot for the clinical team. The majority of the remaining IGRT on oART system are lung SBRT where the first fraction is not adapted due to being within a week of the simulation.

CONCLUSION

Leveraging CBCT based ART system as replan platform and non-ART Linac as IGRT platform is clinically feasible. This process significantly improved the turnaround time for replan, reduced the required resource and promotes the utilization of oART.

Intelligent Interactive Deformable Image Registration for Online Adaptive Radiotherapy

PURPOSE/OBJECTIVE(S)

The goal of this study is to streamline the time-consuming contouring process in online adaptive radiotherapy (ART) by utilizing a deep learning-based interactive deformable image registration (DIR) algorithm. The objective is to minimize manual review and editing of automatically generated initial contours of organs-at-risk (OARs) and targets, thereby improving the efficiency and effectiveness of the treatment process.

MATERIALS/METHODS

Our proposed method reforms the current DIR-based contour propagation method in clinical practice through the implementation of a deep learning-based interactive approach. The steps include: 1) generation of an initial deformable vector field (DVF) using a DL model, based on fixed and moving image pairs, resulting in the initial contours of OARs and targets; 2) clinician review/edit one the OAR/target contours as needed; 3) updated contour is sent to DL model to update the DVF and the remaining OARs/targets contours. Repeat this process until satisfactory contour qualities are achieved. We used the Open Access Series of Imaging Studies (OASIS) as the testbed, including 394 (train) and 20 (test) brain T1-weighted MRI scans, each containing 35 annotated organs. The U-Net architecture was employed to update the DVF from fixed/moving images, initial contours, and updated contours. We compared our approach to traditional manual editing without interaction and quantified the effort reduction using the added path length (APL) metric which is supposed to be proportional to the absolute time spent on the contour editing. We conducted paired t-test to show the significance. For comparison purpose, we assumed the clinicians edit the contours with the largest APL, i.e., the contours that require the most editing efforts.

RESULTS

The editing effort, as measured by APL, was reduced by 18.5% to 25.4% with a mean of 23.3%, median of 23.6%, and standard deviation of 1.9%. The significance of the results was confirmed with a p-value of 1.47e-24.

CONCLUSION

Our study demonstrates a significant reduction in editing effort, as measured by APL, compared to traditional manual contour editing. These results demonstrate the potential of our deep learning-based interactive approach to improve the efficiency and accuracy of the contouring process in clinical practice.

Deep Learning-Based Quality Assurance for Auto-Segmentation Masks in Radiotherapy

PURPOSE/OBJECTIVE(S)

Deep learning-based auto-contouring has shown great promise in several disease sites including GU and head and neck. However, quality assurance (QA) is key to identify poor auto-contours which is time consuming. We hypothesis that training a deep learning model to predict contour quality metrics, such as Dice coefficients (DSC) and associated uncertainties for QA.

MATERIALS/METHODS

We trained a 3D U-Net-based DL model for segmenting the target and three clinical-relevant OARs (bladder and rectum). To mimic the slice-by-slice review process in clinical practice, we then trained a 2D ResNet-based DL model to predict the 2D DSC for each 2D slice's contour, generated by the 3D segmentation model. Using the Monte Carlo dropout technique, we made 20 independent predictions per slice, with the final DSC calculated as their average and uncertainty estimated as 95% prediction intervals (PI). The study cohort consisted of 912 prostate cancer patients who received definitive radiotherapy. The 3D auto-segmentation model was trained on 129 patients and validated on 20, before being tested on 763 patients. The 2D DSC prediction model was trained on 293 patients with 11116 slices, validated on 73 patients with 2804 slices, and tested on 366 patients with 14117 slices. Rectum was chosen to test the 2D contour QA model as it is the most challenging OAR. We categorized 2D slices into three groups based on the lower and upper bounds of the prediction intervals. "no/minor edits" (lower bound > = 0.9), "major edits" (lower bound < 0.9 and upper bound > = 0.8), and "not acceptable" (upper bound < 0.8). The model performance was quantified by calculating correlation coefficients between predicted and ground truth DSC and the fraction of cases that were correctly identified in each category.

RESULTS

The results of the study showed that the overall correlation coefficient between predicted, and ground truth DSC was 0.842. The model was able to correctly identify 78.3%, 60.7%, and 53.4% of the "no/minor edits", "major edits", and "not acceptable" cases, respectively.

CONCLUSION

This study provides a valuable tool for clinicians in making quick decisions on the acceptance, rejection, or revision of auto-segmented masks during the radiation therapy planning process by providing quantitative results on predicted DSC and associated uncertainties.

Dosimetric Comparison of Adaptive Radiotherapy Modalities for Stereotactic Partial Breast Irradiation

PURPOSE/OBJECTIVE(S)

An increase in the availability of adaptive radiotherapy (ART) platforms have proven to be effective in the treatment of a variety of sites. In this study, we aim to evaluate the effectiveness of non-adaptive RT and 3 different ART platforms: (1) CBCT-based, (2) CT-based, and (3) MRI-based for stereotactic partial breast irradiation (SPBI).

MATERIALS/METHODS

Data were collected from 32 patients (16 left and 16 right breast) treated at a single institution. 16 patients (8 left and 8 right) treated using the non-ART platform were re-planned onto two different ART platforms, CBCT- and MRI-based. The remaining 16 patients treated using CT-based adaptive platform were not re-planned due to the prone patient treatment position (others systems supine). All cases were planned to 30 Gy in 5 fractions. Plan quality was evaluated based on pre-defined planning goals to the OARS: ipsilateral and contralateral lungs (Dmean, Dmax, V20 Gy, V9 Gy), ipsilateral (V15 Gy, V30 Gy) and contralateral breasts (Dmax), heart (Dmean, Dmax, V3 Gy, V1.5 Gy), skin (Dmax, V36.5 Gy), and rib (Dmax, V30 Gy). Target goals were defined by Dmax, Dmin, gradient index, and paddock conformality index. Re-planned cases were compared within the cohort using a paired t-test and a 2-sided t-test was used comparing to the CT-based platform.

RESULTS

Comparing the left and right breast cohort across all platforms, the CT-based ART system showed a signification dose reduction in Dmean (p<0.001 for all platforms), Dmax (p<0.001 for left breast, p<0.03 for right breast) and V9 Gy (p<0.004 for left breast, p<0.001 for right breast) to the ipsilateral lung, V15 Gy (p<0.004 for left breast cohort) to the ipsilateral breast, and Dmax to the contralateral breast (p<0.001) and ribs (p = 0.01, p<0.001, p = 0.01 for CBCT-ART, MRI-ART, and non-ART for left breast cohort only). On average, the MR-Linac platform showed the least degree of OAR sparing across nearly all dosimetric parameters evaluated when compared to all modalities, especially for contralateral lung Dmean and Dmax (p<0.05 for all dosimetric parameters for all platforms) and contralateral breast Dmax (p<0.003 for all platforms). The CBCT-based platform showed superior dose reduction in contralateral lung mean (p<0.03 for all platforms) and heart Dmean (p = 0.065, p<0.001, p = 0.045 for non-adaptive, MRI-ART, and CT-ART for left breast and p<0.008 for right breast). PTV coverage was comparable across all platforms, averaging at approximately 95%. The CT-based ART platform showed a significantly reduced gradient index relative to the CBCT- and MRI-based platforms (p<0.001).

CONCLUSION

For SPBI treatments, the CT-based ART platforms displayed a higher degree of OAR sparing for many of the dosimetric parameters recorded relative to the other ART and non-ART platforms presented. The MRI-based system typically showed less reduced OAR sparing; however, the advantage of the system is shown if soft tissue contrast is needed. PTV coverage remained comparable across all platforms.

Real-Time 3D Liver Tumor Localization via Combined Optical Surface Imaging and an X-Ray Projection from Arbitrary Imaging Angles

PURPOSE/OBJECTIVE(S)

A deep learning (DL)-based, deformable registration-driven liver tumor localization technique was developed for onboard deformable motion tracking. The technique (Surf-X360-Bio) uses in-room optical surface imaging and an x-ray projection at an arbitrary scan angle to solve volumetric liver and liver tumor motion in real-time.

MATERIALS/METHODS

Surf-X360-Bio solves the volumetric motion of the liver and localizes the liver tumor, through deforming liver and liver tumor meshes segmented on prior 4D-CTs/CBCTs. It uses real-time onboard information from an optical surface image and a simultaneously-acquired x-ray projection (from an arbitrary scan angle). Surf-X360-Bio localizes tumors via two steps: liver boundary motion estimation and intra-liver motion derivation. Surf-X360-Bio first estimates liver boundary motion by a patient-specific surface imaging model (Surf), utilizing the correlation between the external body surface motion and the internal liver boundary motion. As the correlation can be imperfect, the residual motion estimation errors were corrected by a patient-specific, angle-agnostic x-ray imaging model (X360). X360 deformed the liver boundary to match motion-related imaging features encoded in an arbitrarily-angled x-ray projection, using the Surf output for initialization. X360 adopted a geometry-aware learning module to extract and adapt to angle-varying features of arbitrarily-angled x-ray images, by calculating the projection system matrix of each x-ray image on-the-fly during model training and inference. After the liver boundary motion estimation by Surf and X360, intra-liver deformation was solved by a biomechanical model (Bio) to propagate the liver boundary motion inside to localize the tumors. The DL-based Bio model used domain knowledge of tissue biomechanics and finite element analysis (FEA) to solve intra-liver motion, with a much faster speed than conventional FEA methods to meet the real-time constraint. Surf-X360-Bio was evaluated using a dataset of 34 liver patients. Liver tumor localization accuracy was evaluated by the Dice similarity coefficient (DSC), the 95-percentile Hausdorff distance (HD95), and the center-of-mass error (COME).

RESULTS

Using 3,306 motion scenarios spanning the 360 degree x-ray scan angles for each testing patient, Surf-X360-Bio localized the liver tumors to 0.81 (mean) ± 0.16 (s.d.) in DSC, 2.5 ± 1.7 mm in HD95, and 2.1 ± 1.8 mm in COME. In comparison, the prior reference image without deformable registration-driven localization yielded 0.42 ± 0.29 in DSC, 8.1 ± 5.2 mm in HD95, and 8.5 ± 5.2 mm in COME. Via Surf-X360-Bio, the overall inference time was below 230 ms for each case.

CONCLUSION

Combining optical surface imaging and x-ray imaging, Surf-X360-Bio achieved fast (<250 ms inference time), accurate (mean error < 2.1 mm), and robust liver tumor localizations at arbitrary x-ray scan angles for real-time, marker-less 3D deformable motion tracking.

[The efficacy and safety of ultrafiltration for patients with heart failure: results from a single-center randomized controlled study]

Objective: To evaluate the efficacy and safety of ultrafiltration in patients with heart failure. Methods: One hundred and thirty four cases of patients with heart failure, who hospitalized in our hospital from June 2010 to June 2016 were enrolled in this study. Random serial number was generated using SPSS 22.0 software, patients were then randomly divided into control group and ultrafiltration group with the proportion of 1∶1 (67 cases in each group). Patients in the control group received standard therapy. Patients in the ultrafiltration group received ultrafiltration therapy for 8 hours. Curative effect was evaluated after 8 hours treatment in the control group and after 12 hours in the ultrafiltration group. Following parameters were compared between the two groups: body weight, dyspnea score and 6 minutes walking distance as well as blood pressure, heart rate, Na(+) , K(+) , Cl(-), pH, HCO(3)(-), Hb, PLT, Cr, BUN levels. Results: (1)Two patients died during run-in process and eventually 132 cases were chosen for final analysis (65 cases in control group and 67 cases in the ultrafiltration group). Gender, age, type of heart failure, dyspnea score, body weight at baseline were similar between the two groups. (2)Post therapy, patients' body weight decreased obviously, while dyspnea score and 6 minutes walking distance increased significantly in the ultrafiltration group compared to baseline(all P<0.05), and the improvement was significantly greater compared to control group(all P<0.05). (3)The safety index comparison of two groups: blood pressure, heart rate, Na(+) , K(+) , Cl(-), pH, HCO(3)(-), Hb, PLT, Cr, and BUN were similar between the two groups at baseline and post therapy. Conclusion: Ultrafiltration therapy is safe and effective to treat patients with heart failure.

Monte Carlo dose calculations for high-dose-rate brachytherapy using GPU-accelerated processing

PURPOSE

Current clinical brachytherapy dose calculations are typically based on the Association of American Physicists in Medicine Task Group report 43 (TG-43) guidelines, which approximate patient geometry as an infinitely large water phantom. This ignores patient and applicator geometries and heterogeneities, causing dosimetric errors. Although Monte Carlo (MC) dose calculation is commonly recognized as the most accurate method, its associated long computational time is a major bottleneck for routine clinical applications. This article presents our recent developments of a fast MC dose calculation package for high-dose-rate (HDR) brachytherapy, gBMC, built on a graphics processing unit (GPU) platform.

METHODS AND MATERIALS

gBMC-simulated photon transport in voxelized geometry with physics in (192)Ir HDR brachytherapy energy range considered. A phase-space file was used as a source model. GPU-based parallel computation was used to simultaneously transport multiple photons, one on a GPU thread. We validated gBMC by comparing the dose calculation results in water with that computed TG-43. We also studied heterogeneous phantom cases and a patient case and compared gBMC results with Acuros BV results.

RESULTS

Radial dose function in water calculated by gBMC showed <0.6% relative difference from that of the TG-43 data. Difference in anisotropy function was <1%. In two heterogeneous slab phantoms and one shielded cylinder applicator case, average dose discrepancy between gBMC and Acuros BV was <0.87%. For a tandem and ovoid patient case, good agreement between gBMC and Acruos BV results was observed in both isodose lines and dose-volume histograms. In terms of the efficiency, it took ∼47.5 seconds for gBMC to reach 0.15% statistical uncertainty within the 5% isodose line for the patient case.

CONCLUSIONS

The accuracy and efficiency of a new GPU-based MC dose calculation package, gBMC, for HDR brachytherapy make it attractive for clinical applications.

[Efficacy and safety of a novel ultrafiltration device for treating patients with refractory heart failure]

OBJECTIVE

To evaluate the efficacy and safety of a new ultrafiltration device for treating refractory heart failure patients.

METHODS

A total of 52 patients (37 male, age 29-85(33±44)years) with refractory heart failure were treated using a new ultrafiltration device (FQ-16). Body weight, dyspnea score, oxygen saturation (SatO2), left ventricular ejection fraction(LVEF), BUN, creatinine, electrolytes and blood gas analysis were assessed before and after the treatment. Hypotension event and other main adverse events were recorded.

RESULTS

Ultrafiltration duration ranged between 8-22 hours. Total ultrafiltration volume was (4 489±1 548) ml. Compared with baseline, patients' body weight decreased from (75.3±8.74) kg to (69.8±8.39) kg (P<0.01), dyspnea score improved from 2.47±1.55 to 12.87±3.61 (P<0.01) and SatO2 increased from 91.0±6.01 to 96.4±2.52 (P<0.01) and LVEF increased from (30.0±4.1)% to (36.0±4.3)% (P<0.01) after ultrafiltration. Blood creatinine, BUN, electrolytes and blood gas analysis values were similar at baseline and post ultrafiltration. No hypotension event and other main adverse events occurred during the ultrafiltration treatment.

CONCLUSIONS

The novel ultrafiltration device adequately relieved hypervolemia and dyspnea in patients with refractory heart failure and the treatment process is safe in this patient cohort.

Racial disparities in the association between diabetes mellitus-associated polymorphic locus rs4430796 of the HNF1β gene and prostate cancer: a systematic review and meta-analysis

Polymorphism 17q12 rs4430796 within HNF1β is a genetic variant associated with both diabetes mellitus and prostate cancer, but findings on the correlations of rs4430796 with prostate cancer risk specifically are not in agreement, especially among diverse populations. To shed some light on the contradictory findings, therefore, we carried out a meta-analysis by pooling the odds ratios (ORs) with corresponding 95% confidence intervals (CIs) of all currently available case-control studies located within PubMed and Embase databases up to December 2012. A total of 16 studies comprising 30 datasets that collectively involved 25,535 prostate cancer patients and 25,726 controls were ultimately included in this analysis. The meta-analysis of all the studies revealed that the rs4430796 polymorphism was significantly associated with an increased risk of prostate cancer in all contrast models (ORA vs G = 1.25, 95%CI = 1.21-1.30, POR < 0.001; ORAA vs GG = 1.53, 95%CI = 1.45-1.62, POR < 0.001; ORAG vs GG = 1.24, 95%CI = 1.16-1.34, POR < 0.001; ORAA vs AG+GG = 1.36, 95%CI = 1.30-1.42, POR < 0.001; ORAA+AG vs GG = 1.37, 95%CI = 1.30-1.44, POR < 0.001). After subgroup analyses stratified by ethnicity, however, the rs4430796 polymorphism was significantly associated with prostate cancer in both Caucasians and Asians but not in African-Americans. In conclusion, our meta-analysis identified a significant association between the 17q12 rs4430796 polymorphism and prostate cancer risk, although the degree of this association and frequency of the causative allele varied among men of different races.

Genetic diversity analysis of okra (Abelmoschus esculentus L.) by inter-simple sequence repeat (ISSR) markers

Okra (Abelmoschus esculentus L.) is not only a nutrient-rich vegetable but also an important medicinal herb. Inter-simple sequence repeat (ISSR) markers were employed to investigate the genetic diversity and differentiation of 24 okra genotypes. In this study, the PCR products were separated by electrophoresis on 8% nondenaturing polyacrylamide gel and visualized by silver staining. The 22 ISSR primers produced 289 amplified DNA fragments, and 145 (50%) fragments were polymorphic. The 289 markers were used to construct the dendrogram based on the unweighted pair-group method with arithmetic average (UPGMA) cluster analysis. The dendrogram indicated that 24 okras were clustered into 4 geographically distinct groups. The average polymorphism information content (PIC) was 0.531929, which showed that the majority of primers were informative. The high values of allele frequency, genetic diversity, and heterozygosity showed that primer-sample combinations produced measurable fragments. The mean distances ranged from 0.045455 to 0.454545. The dendrogram indicated that the ISSR markers succeeded in distinguishing most of the 24 varieties in relation to their genetic backgrounds and geographical origins.

Target Localization and Toxicity in Dose-Escalated Prostate Radiotherapy with Image-Guided Approach using Daily Planar Kilovoltage Imaging

Dose escalation with intensity-modulated radiation therapy (IMRT) for carcinoma of the prostate has augmented the need for accurate prostate localization prior to dose delivery. Daily planar kilovoltage (kV) imaging is a low-dose image-guidance technique that is prevalent among radiation oncologists. However, clinical outcomes evaluating the benefit of daily kV imaging are lacking. The purpose of this study was to report our clinical experience, including prostate motion and gastrointestinal (GI) and genitourinary (GU) toxicities, using this modality. A retrospective analysis of 100 patients treated consecutively between December 2005 and March 2008 with definitive external beam IMRT for T1c-T4 disease were included in this analysis. Prescription doses ranged from 74–78 Gy (median, 76) in 2 Gy fractions and were delivered following daily prostate localization using on-board kV imaging (OBI) to localize gold seed fiducial markers within the prostate. Acute and late toxicities were graded as per the NCI CTCAEv3.0. The median follow-up was 22 months. The magnitude and direction of prostate displacement and daily shifts in three axes are reported. Of note, 9.1% and 12.9% of prostate displacements were ≥ 5 mm in the anterior-posterior and superior-inferior directions, respectively. Acute grade 2 GI and GU events occurred in 11% and 39% of patients, respectively, however no grade 3 or higher acute GI or GU events were observed. Regarding late toxicity, 2% and 17% of patients developed grade 2 toxicities, and similarly no grade 3 or higher events had occurred by last follow-up. Thus, kV imaging detected a substantial amount of inter-fractional displacement and may help reduce toxicity profiles, especially high grade events, by improving the accuracy of dose delivery.

Derivation of the tumor position from external respiratory surrogates with periodical updating of the internal/external correlation

In this work we develop techniques that can derive the tumor position from external respiratory surrogates (abdominal surface motion) through periodically updated internal/external correlation. A simple linear function is used to express the correlation between the tumor and surrogate motion. The function parameters are established during a patient setup session with the tumor and surrogate positions simultaneously measured at a 30 Hz rate. During treatment, the surrogate position, constantly acquired at 30 Hz, is used to derive the tumor position. Occasionally, a pair of radiographic images is acquired to enable the updating of the linear correlation function. Four update methods, two aggressive and two conservative, are investigated: (A1) shift line through the update point; (A2) re-fit line through the update point; (C1) re-fit line with extra weight to the update point; (C2) minimize the distances to the update point and previous line fit point. In the present study of eight lung cancer patients, tumor and external surrogate motion demonstrate a high degree of linear correlation which changes dynamically over time. It was found that occasionally updating the correlation function leads to more accurate predictions than using external surrogates alone. In the case of high imaging rates during treatment (greater than 2 Hz) the aggressive update methods (A1 and A2) are more accurate than the conservative ones (C1 and C2). The opposite is observed in the case of low imaging rates.

Effect of statistical uncertainties on Monte Carlo treatment planning

This paper reviews the effect of statistical uncertainties on radiotherapy treatment planning using Monte Carlo simulations. We discuss issues related to the statistical analysis of Monte Carlo dose calculations for realistic clinical beams using various variance reduction or time saving techniques. We discuss the effect of statistical uncertainties on dose prescription and monitor unit calculation for conventional treatment and intensity-modulated radiotherapy (IMRT) based on Monte Carlo simulations. We show the effect of statistical uncertainties on beamlet dose calculation and plan optimization for IMRT and other advanced treatment techniques such as modulated electron radiotherapy (MERT). We provide practical guidelines for the clinical implementation of Monte Carlo treatment planning and show realistic examples of Monte Carlo based IMRT and MERT plans.

A quality assurance phantom for IMRT dose verification

This paper investigates a quality assurance (QA) phantom specially designed to verify the accuracy of dose distributions and monitor units (MU) calculated by clinical treatment planning optimization systems and by the Monte Carlo method for intensity-modulated radiotherapy (IMRT). The QA phantom is a PMMA cylinder of 30 cm diameter and 40 cm length with various bone and lung inserts. A procedure (and formalism) has been developed to measure the absolute dose to water in the PMMA phantom. Another cylindrical phantom of the same dimensions, but made of water, was used to confirm the results obtained with the PMMA phantom. The PMMA phantom was irradiated by 4, 6 and 15 MV photon beams and the dose was measured using an ionization chamber and compared to the results calculated by a commercial inverse planning system (CORVUS, NOMOS, Sewickley, PA) and by the Monte Carlo method. The results show that the dose distributions calculated by both CORVUS and Monte Carlo agreed to within 2% of dose maximum with measured results in the uniform PMMA phantom for both open and intensity-modulated fields. Similar agreement was obtained between Monte Carlo calculations and measured results with the bone and lung heterogeneity inside the PMMA phantom while the CORVUS results were 4% different. The QA phantom has been integrated as a routine QA procedure for the patient's IMRT dose verification at Stanford since 1999.

A Monte Carlo dose calculation tool for radiotherapy treatment planning

A Monte Carlo user code, MCDOSE, has been developed for radiotherapy treatment planning (RTP) dose calculations. MCDOSE is designed as a dose calculation module suitable for adaptation to host RTP systems. MCDOSE can be used for both conventional photon/electron beam calculation and intensity modulated radiotherapy (IMRT) treatment planning. MCDOSE uses a multiple-source model to reconstruct the treatment beam phase space. Based on Monte Carlo simulated or measured beam data acquired during commissioning, source-model parameters are adjusted through an automated procedure. Beam modifiers such as jaws, physical and dynamic wedges, compensators, blocks, electron cut-outs and bolus are simulated by MCDOSE together with a 3D rectilinear patient geometry model built from CT data. Dose distributions calculated using MCDOSE agreed well with those calculated by the EGS4/DOSXYZ code using different beam set-ups and beam modifiers. Heterogeneity correction factors for layered-lung or layered-bone phantoms as calculated by both codes were consistent with measured data to within 1%. The effect of energy cut-offs for particle transport was investigated. Variance reduction techniques were implemented in MCDOSE to achieve a speedup factor of 10-30 compared to DOSXYZ.

Monte Carlo based treatment planning for modulated electron beam radiation therapy

A Monte Carlo based treatment planning system for modulated electron radiation therapy (MERT) is presented. This new variation of intensity modulated radiation therapy (IMRT) utilizes an electron multileaf collimator (eMLC) to deliver non-uniform intensity maps at several electron energies. In this way, conformal dose distributions are delivered to irregular targets located a few centimetres below the surface while sparing deeper-lying normal anatomy. Planning for MERT begins with Monte Carlo generation of electron beamlets. Electrons are transported with proper in-air scattering and the dose is tallied in the phantom for each beamlet. An optimized beamlet plan may be calculated using inverse-planning methods. Step-and-shoot leaf sequences are generated for the intensity maps and dose distributions recalculated using Monte Carlo simulations. Here, scatter and leakage from the leaves are properly accounted for by transporting electrons through the eMLC geometry. The weights for the segments of the plan are re-optimized with the leaf positions fixed and bremsstrahlung leakage and electron scatter doses included. This optimization gives the final optimized plan. It is shown that a significant portion of the calculation time is spent transporting particles in the leaves. However, this is necessary since optimizing segment weights based on a model in which leaf transport is ignored results in an improperly optimized plan with overdosing of target and critical structures. A method of rapidly calculating the bremsstrahlung contribution is presented and shown to be an efficient solution to this problem. A homogeneous model target and a 2D breast plan are presented. The potential use of this tool in clinical planning is discussed.

Derivation of electron and photon energy spectra from electron beam central axis depth dose curves

A method for deriving the electron and photon energy spectra from electron beam central axis percentage depth dose (PDD) curves has been investigated. The PDD curves of 6, 12 and 20 MeV electron beams obtained from the Monte Carlo full phase space simulations of the Varian linear accelerator treatment head have been used to test the method. We have employed a 'random creep' algorithm to determine the energy spectra of electrons and photons in a clinical electron beam. The fitted electron and photon energy spectra have been compared with the corresponding spectra obtained from the Monte Carlo full phase space simulations. Our fitted energy spectra are in good agreement with the Monte Carlo simulated spectra in terms of peak location, peak width, amplitude and smoothness of the spectrum. In addition, the derived depth dose curves of head-generated photons agree well in both shape and amplitude with those calculated using the full phase space data. The central axis depth dose curves and dose profiles at various depths have been compared using an automated electron beam commissioning procedure. The comparison has demonstrated that our method is capable of deriving the energy spectra for the Varian accelerator electron beams investigated. We have implemented this method in the electron beam commissioning procedure for Monte Carlo electron beam dose calculations.

The MLC tongue-and-groove effect on IMRT dose distributions

We have investigated the tongue-and-groove effect on the IMRT dose distributions for a Varian MLC. We have compared the dose distributions calculated using the intensity maps with and without the tongue-and-groove effect. Our results showed that, for one intensity-modulated treatment field, the maximum tongue-and-groove effect could be up to 10% of the maximum dose in the dose distributions. For an IMRT treatment with multiple gantry angles (> or = 5), the difference between the dose distributions with and without the tongue-and-groove effect was hardly visible, less than 1.6% for the two typical clinical cases studied. After considering the patient setup errors, the dose distributions were smoothed with reduced and insignificant differences between plans with and without the tongue-and-groove effect. Therefore, for a multiple-field IMRT plan (> or = 5), the tongue-and-groove effect on the IMRT dose distributions will be generally clinically insignificant due to the smearing effect of individual fields. The tongue-and-groove effect on an IMRT plan with small number of fields (< 5) will vary depending on the number of fields in a plan (coplanar or non-coplanar), the MLC leaf sequences and the patient setup uncertainty, and may be significant (> 5% of maximum dose) in some cases, especially when the patient setup uncertainty is small (< or = 2 mm).

Monte Carlo simulation of electron beams from an accelerator head using PENELOPE

The Monte Carlo code PENELOPE has been used to simulate electron beams from a Siemens Mevatron KDS linac with nominal energies of 6, 12 and 18 MeV. Owing to its accuracy, which stems from that of the underlying physical interaction models, PENELOPE is suitable for simulating problems of interest to the medical physics community. It includes a geometry package that allows the definition of complex quadric geometries, such as those of irradiation instruments, in a straightforward manner. Dose distributions in water simulated with PENELOPE agree well with experimental measurements using a silicon detector and a monitoring ionization chamber. Insertion of a lead slab in the incident beam at the surface of the water phantom produces sharp variations in the dose distributions, which are correctly reproduced by the simulation code. Results from PENELOPE are also compared with those of equivalent simulations with the EGS4-based user codes BEAM and DOSXYZ. Angular and energy distributions of electrons and photons in the phase-space plane (at the downstream end of the applicator) obtained from both simulation codes are similar, although significant differences do appear in some cases. These differences, however, are shown to have a negligible effect on the calculated dose distributions. Various practical aspects of the simulations, such as the calculation of statistical uncertainties and the effect of the 'latent' variance in the phase-space file, are discussed in detail.

Modeling the extrafocal radiation and monitor chamber backscatter for photon beam dose calculation

A simple analytical approach has been developed to model extrafocal radiation and monitor chamber backscatter for clinical photon beams. Model parameters for both the extrafocal source and monitor chamber backscatter are determined simultaneously using conventional measured data, i.e., in-air output factors for square and rectangular fields defined by the photon jaws. The model has been applied to 6 MV and 15 MV photon beams produced by a Varian Clinac 2300C/D accelerator. Contributions to the in-air output factor from the extrafocal radiation and monitor chamber backscatter, as predicted by the model, are in good agreement with the measurements. The model can be used to calculate the in-air output factors analytically, with an accuracy of 0.2% for symmetric or asymmetric rectangular fields defined by jaws when the calculation point is at the isocenter and 0.5% when the calculation point is at an extended SSD. For MLC-defined fields, with the jaws at the recommended positions, calculated in-air output factors agree with the measured data to within 0.3% at the isocenter and 0.7% at off-axis positions. The model has been incorporated into a Monte Carlo dose algorithm to calculate the absolute dose distributions in patients or phantoms. For three MLC-defined irregular fields (triangle shape, C-shape, and L-shape), the calculations agree with the measurements to about 1% even for points at off-axis positions. The model will be particularly useful for IMRT dose calculations because it accurately predicts beam output and penumbra dose.

Monte Carlo and experimental investigations of multileaf collimated electron beams for modulated electron radiation therapy

Modulated electron radiation therapy (MERT) has been proposed as a means of delivering conformal dose to shallow tumors while sparing distal structures and surrounding tissues. Conventional systems for electron beam collimation are labor and time intensive in their construction and are therefore inadequate for use in the sequential delivery of multiple complex fields required by MERT. This study investigates two proposed methods of electron beam collimation: the use of existing photon multileaf collimators (MLC) in a helium atmosphere to reduce in-air electron scatter, and a MLC specifically designed for electron beam collimation. Monte Carlo simulations of a Varian Clinac 2100C were performed using the EGS4/BEAM system and dose calculations performed with the MCDOSE code. Dose penumbras from fields collimated by photon MLCs both with air and with helium at 6, 12, and 20 MeV at a range of SSDs from 70 to 90 cm were examined. Significant improvements were observed for the helium based system. Simulations were also performed on an electron specific MLC located at the level of the last scraper of a 25x25 cm2 applicator. A number of leaf materials, thicknesses, end shapes, and widths were simulated to determine optimal construction parameters. The results demonstrated that tungsten leaves 15 mm thick and 5 mm wide with unfocused ends would provide sufficient collimation for MERT fields. A prototype electron MLC was constructed and comparisons between film measurements and simulation demonstrate the validity of the Monte Carlo model. Further simulations of dose penumbras demonstrate that such an electron MLC would provide improvements over the helium filled photon MLC at all energies, and improvements in the 90-10 penumbra of 12% to 45% at 20 MeV and 6 MeV, respectively. These improvements were also seen in isodose curves when a complex field shape was simulated. It is thus concluded that an MLC specific for electron beam collimation is required for MERT.

Validation of a Monte Carlo dose calculation tool for radiotherapy treatment planning

A new EGS4/PRESTA Monte Carlo user code, MCDOSE, has been developed as a routine dose calculation tool for radiotherapy treatment planning. It is suitable for both conventional and intensity modulated radiation therapy. Two important features of MCDOSE are the inclusion of beam modifiers in the patient simulation and the implementation of several variance reduction techniques. Before this tool can be used reliably for clinical dose calculation, it must be properly validated. The validation for beam modifiers has been performed by comparing the dose distributions calculated by MCDOSE and the well-benchmarked EGS4 user codes BEAM and DOSXYZ. Various beam modifiers were simulated. Good agreement in the dose distributions was observed. The differences in electron cutout factors between the results of MCDOSE and measurements were within 2%. The accuracy of MCDOSE with various variance reduction techniques was tested by comparing the dose distributions in different inhomogeneous phantoms with those calculated by DOSXYZ without variance reduction. The agreement was within 1.0%. Our results demonstrate that MCDOSE is accurate and efficient for routine dose calculation in radiotherapy treatment planning, with or without beam modifiers.

Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system

The purpose of this work was to use Monte Carlo simulations to verify the accuracy of the dose distributions from a commercial treatment planning optimization system (Corvus, Nomos Corp., Sewickley, PA) for intensity-modulated radiotherapy (IMRT). A Monte Carlo treatment planning system has been implemented clinically to improve and verify the accuracy of radiotherapy dose calculations. Further modifications to the system were made to compute the dose in a patient for multiple fixed-gantry IMRT fields. The dose distributions in the experimental phantoms and in the patients were calculated and used to verify the optimized treatment plans generated by the Corvus system. The Monte Carlo calculated IMRT dose distributions agreed with the measurements to within 2% of the maximum dose for all the beam energies and field sizes for both the homogeneous and heterogeneous phantoms. The dose distributions predicted by the Corvus system, which employs a finite-size pencil beam (FSPB) algorithm, agreed with the Monte Carlo simulations and measurements to within 4% in a cylindrical water phantom with various hypothetical target shapes. Discrepancies of more than 5% (relative to the prescribed target dose) in the target region and over 20% in the critical structures were found in some IMRT patient calculations. The FSPB algorithm as implemented in the Corvus system is adequate for homogeneous phantoms (such as prostate) but may result in significant under or over-estimation of the dose in some cases involving heterogeneities such as the air-tissue, lung-tissue and tissue-bone interfaces.

Removing the effect of statistical uncertainty on dose-volume histograms from Monte Carlo dose calculations

Dose-volume histograms (DVHs) of the dose distributions calculated by the Monte Carlo method contain statistical uncertainties. The Monte Carlo DVH can be considered as blurred from the noiseless DVH by the statistical uncertainty. The focus of the present work is on the removal of the statistical uncertainty effect on the Monte Carlo DVHs and the reconstruction of the noiseless DVHs. We first study the effect of statistical uncertainty. It is found that the steeper the DVH, the more significant the effect. For typical critical structure DVHs the effect is usually negligible. For the target DVHs the effect could be clinically significant, depending on the value of uncertainty and the slope of the DVH. We then propose an iterative reconstruction algorithm. Using the DVHs and statistical uncertainties from the Monte Carlo simulations, we are able to reconstruct the noiseless DVHs. A hypothetical example and a number of clinical cases have been used to test the proposed algorithm. For each clinical case, two Monte Carlo simulations (denoted A and B) were performed. Simulation A has very large statistical uncertainties (about 10% of dose in the target volume) while simulation B has very small uncertainties (about 1%). DVHs from simulation B were used to approximate the noiseless DVHs. Using the proposed algorithm, the effect of statistical uncertainty can be removed from the DVHs of simulation A. The reconstructed DVHs were in good agreement with the DVHs from simulation B. The proposed approach is expected to be useful in removing the blurring effect on a quickly calculated Monte Carlo DVH when performing the iterative forward treatment planning.

Energy- and intensity-modulated electron beams for radiotherapy

This work investigates the feasibility of optimizing energy- and intensity-modulated electron beams for radiation therapy. A multileaf collimator (MLC) specially designed for modulated electron radiotherapy (MERT) was investigated both experimentally and by Monte Carlo simulations. An inverse-planning system based on Monte Carlo dose calculations was developed to optimize electron beam energy and intensity to achieve dose conformity for target volumes near the surface. The results showed that an MLC with 5 mm leaf widths could produce complex field shapes for MERT. Electron intra- and inter-leaf leakage had negligible effects on the dose distributions delivered with the MLC, even at shallow depths. Focused leaf ends reduced the electron scattering contributions to the dose compared with straight leaf ends. As anticipated, moving the MLC position toward the patient surface reduced the penumbra significantly. There were significant differences in the beamlet distributions calculated by an analytic 3-D pencil beam algorithm and the Monte Carlo method. The Monte Carlo calculated beamlet distributions were essential to the accuracy of the MERT dose distribution in cases involving large air gaps, oblique incidence and heterogeneous treatment targets (at the tissue-bone and bone-lung interfaces). To demonstrate the potential of MERT for target dose coverage and normal tissue sparing for treatment of superficial targets, treatment plans for a hypothetical treatment were compared using photon beams and MERT.

Er:YAG laser for the treatment of actinic keratoses

BACKGROUND

There is no single optimal treatment for multiple facial actinic keratoses. The existing therapies such as topical 5-fluorouracil, chemical peels, cryotherapy, dermabrasion, and CO2 laser resurfacing can produce prolonged recovery time or are often operator dependent.

OBJECTIVE

The purpose of this study was to investigate another therapeutic modality which provides a shorter recovery time with uniform results. We performed a prospective pilot study investigating the use of the Er:YAG laser for the treatment of multiple facial actinic keratoses.

METHODS

Five patients with multiple facial actinic keratoses were treated with two to three passes of Er:YAG laser. Anesthesia was achieved in all cases by topical application and local infiltration when indicated. All patients were treated with 2.0 J, 5 mm spot size, and a fluence of 10 J/cm2. Clinical and histologic evaluations were performed both pre- and postoperatively.

RESULTS

All patients showed a decrease in the total number of clinical actinic keratoses on the face ranging from 86 to 96%. In addition to the reversal of actinic damage in the epidermis, histologic evidence revealed increased fibroplasia and decreased superficial solar elastosis 3 months after the laser resurfacing. Reepithelialization occurred in 5-8 days, and erythema lasted for about 3-6 weeks after the procedure. There was no evidence of scarring or pigmentary changes in any of the patients during the follow-up period.

CONCLUSION

Er:YAG laser skin resurfacing is a safe and effective treatment for multiple facial actinic keratoses. Histologic data suggest a new zone of collagen deposition occurs in the superficial papillary dermis. Under our current parameters, Er:YAG laser skin resurfacing has a relatively short recovery period and a low risk of scarring. Unlike the CO2 laser, Er:YAG laser skin resurfacing can be performed with topical anesthesia alone.

Photon beam characterization and modelling for Monte Carlo treatment planning

Photon beams of 4, 6 and 15 MV from Varian Clinac 2100C and 2300C/D accelerators were simulated using the EGS4/BEAM code system. The accelerators were modelled as a combination of component modules (CMs) consisting of a target, primary collimator, exit window, flattening filter, monitor chamber, secondary collimator, ring collimator, photon jaws and protection window. A full phase space file was scored directly above the upper photon jaws and analysed using beam data processing software, BEAMDP, to derive the beam characteristics, such as planar fluence, angular distribution, energy spectrum and the fractional contributions of each individual CM. A multiple-source model has been further developed to reconstruct the original phase space. Separate sources were created with accurate source intensity, energy, fluence and angular distributions for the target, primary collimator and flattening filter. Good agreement (within 2%) between the Monte Carlo calculations with the source model and those with the original phase space was achieved in the dose distributions for field sizes of 4 cm x 4 cm to 40 cm x 40 cm at source surface distances (SSDs) of 80-120 cm. The dose distributions in lung and bone heterogeneous phantoms have also been found to be in good agreement (within 2%) for 4, 6 and 15 MV photon beams for various field sizes between the Monte Carlo calculations with the source model and those with the original phase space.

Electron beam modeling and commissioning for Monte Carlo treatment planning

A hybrid approach for commissioning electron beam Monte Carlo treatment planning systems has been studied. The approach is based on the assumption that accelerators of the same type have very similar electron beam characteristics and the major difference comes from the on-site tuning of the electron incident energy at the exit window. For one type of accelerator, a reference machine can be selected and simulated with the Monte Carlo method. A multiple source model can be built on the full Monte Carlo simulation of the reference beam. When commissioning electron beams from other accelerators of the same type, the energy spectra in the source model are tuned to match the measured dose distributions. A Varian Clinac 2100C accelerator was chosen as the reference machine and a four-source beam model was established based on the Monte Carlo simulations. This simplified beam model can be used to generate Monte Carlo dose distributions accurately (within 2%/2 mm compared to those calculated with full phase space data) for electron beams from the reference machine with various nominal energies, applicator sizes, and SSDs. Three electron beams were commissioned by adjusting the energy spectra in the source model. The dose distributions calculated with the adjusted source model were compared with the dose distributions calculated using the phase space data for these beams. The agreement is within 1% in most of cases and 2% in all situations. This preliminary study has shown the capability of the commissioning approach for handling large variation in the electron incident energy. The possibility of making the approach more versatile is also discussed.

Monte Carlo modelling of electron beams from medical accelerators

Monte Carlo simulation of radiation transport is considered to be one of the most accurate methods of radiation therapy dose calculation. With the rapid development of computer technology, Monte Carlo based treatment planning for radiation therapy is becoming practical. A basic requirement for Monte Carlo treatment planning is a detailed knowledge of the radiation beams from medical accelerators. A practical approach to obtain the above is to perform Monte Carlo simulation of radiation transport in the medical accelerator. Additionally, Monte Carlo modelling of the treatment machine head can also improve our understanding of clinical beam characteristics, help accelerator design and improve the accuracy of clinical dosimetry by providing more realistic beam data. This paper summarizes work over the past two decades on Monte Carlo simulation of clinical electron beams from medical accelerators.

Extracorporeal photochemotherapy for cutaneous T-cell lymphoma: a 9.7-year experience

Cutaneous T-cell lymphoma (CTCL) is an indolent lymphoma usually of CD4+ T lymphocytes in which the aggressive treatment for the advanced stages does not increase survival. Photopheresis has been established as an alternative modality for the therapy of erythrodermic CTCL and reportedly improves survival in patients with advanced stages of the disease. The objective of this study is to review the experience of treating patients with erythrodermic CTCL with extracorporeal photochemotherapy (ECP) at the New York Veteran Affairs Medical Center/NYU Medical Center between September 1987 and April 1997. Forty-one patients with erythrodermic CTCL (stages III and IV) received photopheresis; 25 of them fulfilled the inclusion criterion, i.e., the completion of greater than or equal to 6 cycles of photopheresis. Skin score was defined as a product of severity and percentage of involved surface area. Complete clinical response was defined as disappearance of measurable disease for at least one month, and partial response was defined as greater than or equal to 50% clearance of measurable disease for at least one month. The profile of the patients was: 20 men, 5 women; average age: 64.2 years; 17 patients had stage III disease, and 8 had stage IV disease. Five of the 25 patients (20%) achieved complete clinical response, another 15 (60%) had partial response, and 5 (20%) had no response. The mean time (+/- SD) to achieve complete clinical clearance was 12.6 +/- 10 months (range: 4-30 months) and the mean time (+/- SD) to obtain partial clinical response, including complete response, was 9.7 +/- 5.3 months (range: 4-17 months). Remission duration ranged from 9 to 67 months. The median survival time from the time of initiation of photopheresis is estimated at 70 months. The complete responder group had a lower median CD4/CD8 ratio compared to the non-responders at baseline (3.8 vs 7.2, respectively), although the difference was not statistically significant (P = 0.40). At the time of maximal response, the CD4/CD8 ratio of the complete responder group decreased towards normal values (median = 1.2), whereas this ratio increased among the non-responders (median = 11.0; P = 0.04). Side effects were minimal. Extracorporeal photochemotherapy is an effective and safe treatment for erythrodermic CTCL. In some of these patients, it can induce a long-term and complete clinical remission.

Microcystic adnexal carcinoma: collaborative series review and update

BACKGROUND

Microcystic adnexal carcinoma (MAC) is a malignant appendageal tumor first described in 1982. It can be clinically and histologically confused with other malignant and benign cutaneous neoplasms, leading to inadequate initial treatment. This neoplasm is locally aggressive and deeply infiltrating, characterized by high morbidity and frequent recurrence. Mohs micrographic surgery has been used to conserve tissue and improve the likelihood for cure.

OBJECTIVE

We report our experience using Mohs micrographic surgery for the treatment of MAC and compare with earlier reports in the literature. In addition, we review the epidemiology, clinical and histologic characteristics, and optimal treatment of this rare neoplasm. We also describe a 15-year-old white male patient with MAC on the scalp occurring only 7 years after radiation exposure.

METHODS

The medical records of 11 patients with MAC who were treated by Mohs micrographic surgery were reviewed at both departments, and follow-up data were obtained.

RESULTS

In all patients treated with Mohs micrographic surgery, there were no recurrences after a mean follow-up of 5 years.

CONCLUSION

Mohs technique enables the detection of clinically unrecognizable tumor spread and perineural invasion often encountered with MAC. Aggressive initial treatment by microscopically controlled excision appears to offer the greatest likelihood of cure for this neoplasm, while providing conservation of normal tissue. In addition, we describe the second youngest patient with MAC and readdress the issue of previous radiotherapy as an important predisposing factor.

Atypical erythema toxicum neonatorum of delayed onset in a term infant

We describe a 10-day-old term infant who presented to the emergency room with an acute pustular eruption. Laboratory tests and clinical outcome confirmed the diagnosis of erythema toxicum neonatorum. A full septic workup was performed and all cultures were negative. Wright-stained smear of pustular contents showed a predominance of neutrophils with 10% eosinophils. The white blood cell count was 19,000/mm3 with 10% eosinophils. The eruption resolved spontaneously at 15 days of age leaving no sequelae. This is the first fully documented case of erythema toxicum in a term infant occurring as late as 10 days of age. When erythema toxicum presents in an atypical fashion, diagnostic tests are important to exclude other causes of pustular dermatoses of the neonate.

Dose perturbation caused by high-density inhomogeneities in small beams in stereotactic radiosurgery

The influence of high-density tissue heterogeneities in small-diameter beams used in stereotactic radiosurgery has been investigated. Dose perturbation immediately behind aluminium sheets, used to simulate a high-density tissue inhomogeneity such as bone, was studied in a solid water phantom. Dose reduction factors (DRFs), which are the ratios of the dose in the presence of the inhomogeneity to dose in a uniform density solid water phantom, were measured with a diamond detector for three thicknesses of aluminium. DRFs exhibit dependence on both the inhomogeneity thickness and the beam diameter. The DRF decreases with inhomogeneity thickness. The DRF initially decreases with increase in the beam diameter from 12.5 to 25 mm. For fields greater than 25 mm, the DRFs are nearly constant. The commonly used algorithms such as the TAR ratio method underestimate the magnitude of the measured effect. A good agreement between these measurements and Monte Carlo calculations is obtained. The influence of the high-density inhomogeneity on the tissue maximum ratio (TMR) was also measured with the inhomogeneity at a fixed depth dmax from the entrance surface. The TMR is reduced for all detector-inhomogeneity distances investigated. The dose build-up phenomenon observed in the presence of low-density air inhomogeneity is absent in the presence of a high-density inhomogeneity. The beam width (defined by 50% dose points) immediately beyond the inhomogeneity is unaffected by the high-density inhomogeneity. However, the 90%-10% and 80%-20% dose penumbra widths and the dose outside the beam edge (beyond the 50% dose point) are reduced. This reduction in dose outside the beam edge is caused by the reduced range of the secondary radiation (photons and electrons) in the high-density medium.

Do we need Monte Carlo treatment planning for linac based radiosurgery? A case study

The accuracy of conventional empirical and semi-empirical dose calculation algorithms for radiation therapy treatment planning is limited. The main problem is that these algorithms fail to adequately consider the lateral transport of radiation. Most conventional algorithms use measured dose distribution data as input. These data induce an added inaccuracy to stereotactic radiosurgery dose calculations due to the difficulty of acquiring accurate dosimetric data for very small beams; however, since multiple arcs of large solid angles are usually used in stereotactic radiosurgery, the errors introduced by conventional dose algorithms are quite likely to be diluted. The use of Monte Carlo treatment planning for stereotactic radiosurgery has been investigated and described in the present paper. The OMEGA Monte Carlo code system is used as the dose engine in an in-house developed radiosurgery treatment planning system. The Monte Carlo treatment plans are done for two typical clinical cases. In one case, the collimator of 20 mm diameter is used and the lesion is located in the peripheral part of the brain. In the other case, the collimator diameter is 30 mm and the lesion is in the central part of the brain. The resultant dose distributions are compared with those calculated with a conventional dose algorithm which is based on the standard Tissue Maximum Ratio (TMR)/Off Axis Ratio (OAR) formalism. Without the inhomogeneity correction, the conventional algorithm yields accurate relative dose distributions for both cases compared with the Monte Carlo calculations. The absolute dose at the isocenter may be overestimated by the conventional algorithm by 1.5% for the first case and 2.6% for the second case; however, using the method of ratio of TMRs for inhomogeneity correction, the overestimation can be greatly reduced for both cases. The inclusion of the inhomogeneity correction into the conventional dose algorithm does not alter the relative dose distributions. Based on the clinical cases studied, it may be concluded that the conventional dose algorithm is sufficient for radiosurgery treatment planning and the Monte Carlo based radiosurgery treatment planning is unwarranted.

On compensator design for photon beam intensity-modulated conformal therapy

Recently the compensator has been shown to be an in expensive and reliable dose delivery device for photon beam intensity-modulated radiation therapy (IMRT). The goal of IMRT compensator design is to produce an optimized primary fluence profile at the patient's surface obtained from the optimization procedure. In this paper some of the problems associated with IMRT compensator design, specifically the beam perturbations caused by the compensator, are discussed. A simple formula is derived to calculate the optimal compensator thickness profile from an optimized primary fluence profile. The change of characteristics of a 6 MV beam caused by the introduction of cerrobend compensators in the beam is investigated using OMEGA Monte Carlo codes. It is found that the compensator significantly changes the energy spectrum and the mean energy of the primary photons at the patient's surface. However, beam hardening does not have as significant an effect on the percent depth dose as it does on the energy spectrum. We conclude that in most situations the beam hardening effect can be ignored during compensator design and dose calculation. The influence of the compensator on the contaminant electron buildup dose is found to be small and independent of the compensator thickness of interest. Therefore, it can be ignored in the compensator design and included as a correction into the final dose distribution. The scattered photons from the compensator are found to have no effect on the surface dose. These photons produce a uniform low fluence distribution at the patient's surface, which is independent of compensator shape. This is also true for very large fields and extremely asymmetric and nonuniform compensator thickness profiles. Compared to the primary photons, the scattered photons have much larger angular spread and similar energy spectrum at the patient's surface. These characteristics allow the compensator thickness profile and the dose distribution to be calculated from the optimized fluence profile of primary photons, without considering the scattered photons.

A fast numerical algorithm for electron mean energy calculation in radiation therapy

A numerical algorithm for calculating the mean energy of radiotherapy electron beams has been developed. This algorithm is fast and accurate which makes it suitable for routine clinical use. First, a Gaussian distribution of the electron energy spectrum is derived from the linear Boltzmann equation. Based on this Gaussian spectrum, and after introducing a correction to the CSDA mean energy, a recursive-iterative algorithm for the mean energy calculation is developed. The multiple-scattering correction is taken into account using Yang's pathlength distribution theory. Numerical results of the present algorithm are compared with the results obtained through Monte Carlo simulation as well as Harder's formula. Good agreement with Monte Carlo simulation is achieved. Also the new algorithm is much more accurate than the commonly-used empirical formula of Harder.

Efficacy of ultraviolet B phototherapy for psoriasis in patients infected with human immunodeficiency virus

To evaluate the efficacy of ultraviolet B (UVB) phototherapy for the treatment of psoriasis in patients infected with human immunodeficiency virus (HIV), the response of 14 patients was compared to that of matched seronegative control individuals. All patients were evaluated prior to treatment (baseline) and after 21 treatments for the extent of total body surface area (TBSA) involvement and the quantification of scale, erythema, and thickness of plaques using a scale of 0 (absent) to 4 (severe). The only concomitant medication allowed was salicylic acid in petrolatum. The cumulative score for scale, erythema, and thickness improved 1.9 +/- 0.5 [mean +/- standard error of mean (SEM)] in the HIV group and 2.4 +/- 0.3 in controls. There was 40.9 +/- 7.3% reduction of TBSA involvement in the former and 38.4 +/- 7.6% reduction in the latter group. None of the differences was statistically significant. There was no statistically significant difference in the response to therapy among various stages of immunosuppression in the HIV group. There was also no deterioration of immune status in this group. UVB phototherapy is an effective treatment for psoriasis in patients infected with HIV. The response is identical to that of matched control individuals.

Effects of tubeimoside-1 on HIV core protein p24 and cytopathogenesis in vitro

Tubeimoside-1 (Tub) is a triterpenoid saponin first isolated in China from Bolbostemma paniculatum (Maxim) Franquet, Cucurbitaceae. To find out whether Tub has any anti-infective activity on human immunodeficiency virus (HIV), we tested its effects on HIV core protein p24 (with an ELISA) and HIV-mediated cytopathogenesis (using a colorimetric assay). The results showed that Tub inhibited both p24 production and cytopathogenesis mediated by HTLV-IIIB, and the median effective concentrations (EC50) were 24.1 and 22.9 micrograms.ml-1, respectively. Tub also effectively neutralized the infection of 2 other isolates, HTLV-IIIRF and HTLV-IIIMN. It is concluded that Tub had an inhibitory action on the infection of HIV-1 isolates and would be a promising candidate for treatment of AIDS.