Modelling DTPA therapy following Am contamination in rats

A major challenge in modelling the decorporation of actinides (An), such as americium (Am), with DTPA (diethylenetriaminepentaacetic acid) is the fact that standard biokinetic models become inadequate for assessing radionuclide intake and estimating the resulting dose, as DTPA perturbs the regular biokinetics of the radionuclide. At present, most attempts existing in the literature are empirical and developed mainly for the interpretation of one or a limited number of specific incorporation cases. Recently, several approaches have been presented with the aim of developing a generic model, one of which reported the unperturbed biokinetics of plutonium (Pu), the chelation process and the behaviour of the chelated compound An-DTPA with a single model structure. The aim of the approach described in this present work is the development of a generic model that is able to describe the biokinetics of Am, DTPA and the chelate Am-DTPA simultaneously. Since accidental intakes in humans present many unknowns and large uncertainties, data from controlled studies in animals were used. In these studies, different amounts of DTPA were administered at different times after contamination with known quantities of Am. To account for the enhancement of faecal excretion and reduction in liver retention, DTPA is assumed to chelate Am not only in extracellular fluids, but also in hepatocytes. A good agreement was found between the predictions of the proposed model and the experimental results for urinary and faecal excretion and accumulation and retention in the liver. However, the decorporation from the skeletal compartment could not be reproduced satisfactorily under these simple assumptions.

Patient-specific biokinetics and hybrid 2D/3D approach integration in OEDIPE software: Application to radioiodine therapy

BACKGROUND

The progression of targeted radionuclide therapy requires the development of dosimetry software accounting for patient-specific biokinetics. New functionalities were thus developed in the OEDIPE software, to deal with multiple 3D images or multiple planar images and a SPECT image.

MATERIEL & METHOD

Methods were implemented to recover patient biokinetics in volumes of interest. If several 3D SPECT images are available, they are registered to a reference CT scan. When several planar images and a single SPECT are available, the planar images are registered to the SPECT and counts of the planar images converted to activity. To validate these developments, six SPECT/CT and planar images of a Jaszczak phantom containing I-131 were acquired at different dates. Cumulated activity was estimated in each sphere using the SPECT/CT images only or the planar series associated to one SPECT/CT. Biokinetics and doses in lesions and in the lungs of a patient treated with I-131 for differentiated thyroid cancer were then estimated using four planar images and a SPECT/CT scan. Whole-body retention data were used to compare the biokinetics obtained from the planar and SPECT data.

RESULTS

Activities and cumulated activities estimated using OEDIPE in the phantom spheres agreed well with the reference values for both approaches. Results obtained for the patient compared well with those derived from whole-body retention data.

CONCLUSION

The implemented features allow automatic evaluation of patient-specific biokinetics from different series of patient images, enabling patient-specific dosimetry without the need for external software to estimate the cumulated activities in different VOIs.

Breast cancer risk among thyroid cancer survivors and the role of I-131 treatment

BACKGROUND

Female thyroid cancer survivors are more likely to have a higher risk of breast cancer compared to the general population, and the underlying causes are yet to be understood. The potential role of I-131 treatment on this association remains controversial.

METHODS

We pooled individual data of women who were treated for differentiated thyroid cancer from 1934 to 2005 in France, Italy and Sweden. Standardized incidence ratios (SIRs) for breast cancer were estimated by comparison with age, sex and calendar-year expected values of the general population in each country. We estimated breast cancer risk in relation to I-131 treatment using time-dependent Poisson models.

RESULTS

Of 8475 women (mean age at diagnosis: 45 years, range 2-90 years), 335 were diagnosed with breast cancer [SIR = 1.52, 95% confidence interval (CI): 1.36-1.69] during a median follow-up time of 12.7 years since diagnosis. Overall, breast cancer risk did not differ between women treated or not with I-131 (relative risk=1.07, 95% CI 0.84-1.35). However, breast cancer risk increased with increasing cumulative I-131 activity, without significant departure from linearity (excess relative risk per 100 mCi=17%, 95% CI: 2% to 38%). The higher risk associated with a cumulative I-131 activity of ≥100 mCi and ≥400 mCi was translated into 4 (95% CI -4 to 13) and 42 (95% CI -8 to 93) excess breast cancer cases per 10,000 person-years, respectively.

CONCLUSIONS

An elevated risk was observed for the highest cumulative administered activity (>=400 mCi), and a significant dose-dependent association was observed among thyroid cancer survivors who were treated with I-131. However, overall, I-131 treatment might only explain partly the increase in breast cancer risk among female thyroid cancer survivors.

DTPA Treatment of Wound Contamination in Rats with Americium: Evaluation of Urinary Profiles Using STATBIODIS Shows Importance of Prompt Administration

ABSTRACT

In the nuclear industry, wound contamination with americium is expected to increase with decommissioning and waste management. Treatment of workers with diethylenetriaminepentaacetic acid (DTPA) requires optimization to reduce internal contamination and radiation exposure. This work aimed at evaluating and comparing different DTPA protocol efficacies after wound contamination of rats with americium. Wound contamination was simulated in rats by depositing americium nitrate in an incision in the hind limb. Different routes, times, and frequencies of DTPA administration were evaluated. Individual daily urinary americium excretion and tissue retention were analyzed using the statistical tool STATBIODIS. Urinary profiles, urinary enhancement factors, and inhibition percentages of tissue retention were calculated. A single DTPA administration the day of contamination induced a rapid increase in americium urinary excretion that decreased exponentially over 7 d, indicating that the first DTPA administration should be delivered as early as possible. DTPA treatment limited americium uptake in systemic tissues irrespective of the protocol. Liver and skeleton burdens were markedly reduced, which would drive reduction of radiation dose. Local or intravenous injections were equally effective. Inherent difficulties in wound site activity measurements did not allow identification of a significant decorporating effect at the wound site. Repeated intravenous injections of DTPA also increased americium urinary excretion, which supports the use of multiple DTPA administrations shortly after wound contamination. Results from these statistical analyses will contribute to a better understanding of americium behavior in the presence or absence of DTPA and may aid optimization of treatment for workers.

Evidence for a differential translocation of actinides across human lung epithelial cell monolayer in vitro according to their physicochemical properties and the presence of a chelating agent

The epithelial cell plays a key role in the transfer of radionuclides from lungs to blood following pulmonary exposure. The present study was designed to evaluate the transfer across human lung epithelial cells of various actinides (plutonium, americium and uranium), the influence of the physicochemical properties of plutonium compounds and of the chelating agent diethylene triamine pentaacetic acid (DTPA). To address this question, Calu-3 cells grown in a bicameral culture system were used. The integrity of the epithelial barrier was evaluated by measuring transepithelial electrical resistance (TEER) and the passage of a fluorescent marker, lucifer yellow. Activity measurement in basal compartment following periodic collection of culture medium was made from 2 h to seven days. To facilitate data handling and analysis, the statistical tool STATBIODIS was used. The results indicate differences in transfer for the different elements, and according to Pu physicochemical properties. Though to various extents, the chelating agent DTPA always increased the transfer of Pu and Am across the epithelial cells, without altering the integrity of the epithelial barrier. This in vitro cell culture model, by mimicking translocation of actinides from lungs to blood, can represent a valuable tool to further understand the underlying mechanisms and properties controlling this process.

Americium biodistribution in rats after wound contamination with different physicochemical forms in the presence or absence of plutonium: analyses using STATBIODIS

Americium (Am) biodistribution data obtained after wound contamination in rats were analysed to evaluate and quantify the influence of different physicochemical forms of Am in the presence or absence of plutonium (Pu). The biodistribution data were individual Am daily urinary excretion and tissue retention. The data were analysed with STATBIODIS, a statistical tool developed in the laboratory and based on the R language. Non-parametric methods were selected to comply with the data characteristics. Am systemic tissue retention and urinary excretion data were much greater for contamination with soluble physicochemical forms than insoluble forms. Meanwhile, Am relative biodistribution between the main retention tissues (skeleton, liver and kidney) remained the same. Hence, after absorption into blood the radionuclide behaviour was independent of the physicochemical form. The presence of Pu did not change the Am biodistribution. Comparisons of the biodistribution data from the laboratory with mean values published by other laboratories showed that soluble to moderately soluble forms of Am resulted in similar urine excretion after contamination, whether it was intravenous, intramuscular, subcutaneous injection or incision. Findings from this work will contribute to improve the understanding and interpretation of wound contamination cases with different physicochemical forms and mixtures of actinides including Am.

Actinide-contaminated Skin: Comparing Decontamination Efficacy of Water, Cleansing Gels, and DTPA Gels

Skin contamination by alpha-emitting actinides is a risk to workers during nuclear fuel production and reactor decommissioning. Also, the list of items for potential use in radiological dispersal devices includes plutonium and americium. The actinide chemical form is important and solvents such as tributyl phosphate, used to extract plutonium, can influence plutonium behavior. This study investigated skin fixation and efficacy of decontamination products for these actinide forms using viable pig skin in the Franz cell diffusion system. Commonly used or recommended decontamination products such as water, cleansing gel, diethylenetriamine pentaacetic acid, or octadentate hydroxypyridinone compound 3,4,3-LI(1,2-HOPO), as well as diethylenetriamine pentaacetic acid hydrogel formulations, were tested after a 2-h contact time with the contaminant. Analysis of skin samples demonstrated that more plutonium nitrate is bound to skin as compared to plutonium-tributyl phosphate, and fixation of americium to skin was also significant. The data show that for plutonium-tributyl phosphate all the products are effective ranging from 80 to 90% removal of this contaminant. This may be associated with damage to the skin by this complex and suggests a mechanical/wash-out action rather than chelation. For removal of americium and plutonium, both Trait Rouge cleansing gel and diethylenetriamine pentaacetic acid are better than water, and diethylenetriamine pentaacetic acid hydrogel is better than Osmogel. The different treatments, however, did not significantly affect the activity in deeper skin layers, which suggests a need for further improvement of decontamination procedures. The new diethylenetriamine pentaacetic acid hydrogel preparation was effective in removing americium, plutonium, and plutonium-tributyl phosphate from skin; such a formulation offers advantages and thus merits further assessment.

Actinide bioimaging in tissues: Comparison of emulsion and solid track autoradiography techniques with the iQID camera

This work presents a comparison of three autoradiography techniques for imaging biological samples contaminated with actinides: emulsion-based, plastic-based autoradiography and a quantitative digital technique, the iQID camera, based on the numerical analysis of light from a scintillator screen. In radiation toxicology it has been important to develop means of imaging actinide distribution in tissues as these radionuclides may be heterogeneously distributed within and between tissues after internal contamination. Actinide distribution determines which cells are exposed to alpha radiation and is thus potentially critical for assessing absorbed dose. The comparison was carried out by generating autoradiographs of the same biological samples contaminated with actinides with the three autoradiography techniques. These samples were cell preparations or tissue sections collected from animals contaminated with different physico-chemical forms of actinides. The autoradiograph characteristics and the performances of the techniques were evaluated and discussed mainly in terms of acquisition process, activity distribution patterns, spatial resolution and feasibility of activity quantification. The obtained autoradiographs presented similar actinide distribution at low magnification. Out of the three techniques, emulsion autoradiography is the only one to provide a highly-resolved image of the actinide distribution inherently superimposed on the biological sample. Emulsion autoradiography is hence best interpreted at higher magnifications. However, this technique is destructive for the biological sample. Both emulsion- and plastic-based autoradiography record alpha tracks and thus enabled the differentiation between ionized forms of actinides and oxide particles. This feature can help in the evaluation of decorporation therapy efficacy. The most recent technique, the iQID camera, presents several additional features: real-time imaging, separate imaging of alpha particles and gamma rays, and alpha activity quantification. The comparison of these three autoradiography techniques showed that they are complementary and the choice of the technique depends on the purpose of the imaging experiment.

Analysis methodology and development of a statistical tool for biodistribution data from internal contamination with actinides

The aim of this work was to develop a computational tool that integrates several statistical analysis features for biodistribution data from internal contamination experiments. These data represent actinide levels in biological compartments as a function of time and are derived from activity measurements in tissues and excreta. These experiments aim at assessing the influence of different contamination conditions (e.g. intake route or radioelement) on the biological behavior of the contaminant. The ever increasing number of datasets and diversity of experimental conditions make the handling and analysis of biodistribution data difficult. This work sought to facilitate the statistical analysis of a large number of datasets and the comparison of results from diverse experimental conditions. Functional modules were developed using the open-source programming language R to facilitate specific operations: descriptive statistics, visual comparison, curve fitting, and implementation of biokinetic models. In addition, the structure of the datasets was harmonized using the same table format. Analysis outputs can be written in text files and updated data can be written in the consistent table format. Hence, a data repository is built progressively, which is essential for the optimal use of animal data. Graphical representations can be automatically generated and saved as image files. The resulting computational tool was applied using data derived from wound contamination experiments conducted under different conditions. In facilitating biodistribution data handling and statistical analyses, this computational tool ensures faster analyses and a better reproducibility compared with the use of multiple office software applications. Furthermore, re-analysis of archival data and comparison of data from different sources is made much easier. Hence this tool will help to understand better the influence of contamination characteristics on actinide biokinetics. Our approach can aid the optimization of treatment protocols and therefore contribute to the improvement of the medical response after internal contamination with actinides.

BODY SIZE-SPECIFIC EFFECTIVE DOSE CONVERSION COEFFICIENTS FOR CT SCANS

Effective dose from computed tomography (CT) examinations is usually estimated using the scanner-provided dose-length product and using conversion factors, also known as k-factors, which correspond to scan regions and differ by age according to five categories: 0, 1, 5, 10 y and adult. However, patients often deviate from the standard body size on which the conversion factor is based. In this study, a method for deriving body size-specific k-factors is presented, which can be determined from a simple regression curve based on patient diameter at the centre of the scan range. Using the International Commission on Radiological Protection reference paediatric and adult computational phantoms paired with Monte Carlo simulation of CT X-ray beams, the authors derived a regression-based k-factor model for the following CT scan types: head-neck, head, neck, chest, abdomen, pelvis, abdomen-pelvis (AP) and chest-abdomen-pelvis (CAP). The resulting regression functions were applied to a total of 105 paediatric and 279 adult CT scans randomly sampled from patients who underwent chest, AP and CAP scans at the National Institutes of Health Clinical Center. The authors have calculated and compared the effective doses derived from the conventional age-specific k-factors with the values computed using their body size-specific k-factor. They found that by using the age-specific k-factor, paediatric patients tend to have underestimates (up to 3-fold) of effective dose, while underweight and overweight adult patients tend to have underestimates (up to 2.6-fold) and overestimates (up to 4.6-fold) of effective dose, respectively, compared with the effective dose determined from their body size-dependent factors. The authors present these size-specific k-factors as an alternative to the existing age-specific factors. The body size-specific k-factor will assess effective dose more precisely and on a more individual level than the conventional age-specific k-factors and, hence, improve awareness of the true exposure, which is important for the clinical community to understand.

S values for 131I based on the ICRP adult voxel phantoms

To improve the estimates of organ doses from nuclear medicine procedures using (131)I, the authors calculated a comprehensive set of (131)I S values, defined as absorbed doses in target tissues per unit of nuclear transition in source regions, for different source and target combinations. The authors used the latest reference adult male and female voxel phantoms published by the International Commission on Radiological Protection (ICRP Publication 110) and the (131)I photon and electron spectra from the ICRP Publication 107 to perform Monte Carlo radiation transport calculations using MCNPX2.7 to compute the S values. For each phantom, the authors simulated 55 source regions with an assumed uniform distribution of (131)I. They computed the S values for 42 target tissues directly, without calculating specific absorbed fractions. From these calculations, the authors derived a comprehensive set of S values for (131)I for 55 source regions and 42 target tissues in the ICRP male and female voxel phantoms. Compared with the stylised phantoms from Oak Ridge National Laboratory (ORNL) that consist of 22 source regions and 24 target regions, the new data set includes 1662 additional S values corresponding to additional combinations of source-target tissues that are not available in the stylised phantoms. In a comparison of S values derived from the ICRP and ORNL phantoms, the authors found that the S values to the radiosensitive tissues in the ICRP phantoms were 1.1 (median, female) and 1.3 (median, male) times greater than the values based on the ORNL phantoms. However, for several source-target pairs, the difference was up to 10-fold. The new set of S values can be applied prospectively or retrospectively to the calculation of radiation doses in adults internally exposed to (131)I, including nuclear medicine patients treated for thyroid cancer or hyperthyroidism.

Organ Dose Estimates for Hyperthyroid Patients Treated with (131)I: An Update of the Thyrotoxicosis Follow-Up Study

The Thyrotoxicosis Therapy Follow-up Study (TTFUS) is comprised of 35,593 hyperthyroid patients treated from the mid-1940s through the mid-1960s. One objective of the TTFUS was to evaluate the long-term effects of high-dose iodine-131 ((131)I) treatment (1-4). In the TTFUS cohort, 23,020 patients were treated with (131)I, including 21,536 patients with Graves disease (GD), 1,203 patients with toxic nodular goiter (TNG) and 281 patients with unknown disease. The study population constituted the largest group of hyperthyroid patients ever examined in a single health risk study. The average number (± 1 standard deviation) of (131)I treatments per patient was 1.7 ± 1.4 for the GD patients and 2.1 ± 2.1 for the TNG patients. The average total (131)I administered activity was 380 ± 360 MBq for GD patients and 640 ± 550 MBq for TNG patients. In this work, a biokinetic model for iodine was developed to derive organ residence times and to reconstruct the radiation-absorbed doses to the thyroid gland and to other organs resulting from administration of (131)I to hyperthyroid patients. Based on (131)I data for a small, kinetically well-characterized sub-cohort of patients, multivariate regression equations were developed to relate the numbers of disintegrations of (131)I in more than 50 organs and tissues to anatomical (thyroid mass) and clinical (percentage thyroid uptake and pulse rate) parameters. These equations were then applied to estimate the numbers of (131)I disintegrations in the organs and tissues of all other hyperthyroid patients in the TTFUS who were treated with (131)I. The reference voxel phantoms adopted by the International Commission on Radiological Protection (ICRP) were then used to calculate the absorbed doses in more than 20 organs and tissues of the body. As expected, the absorbed doses were found to be highest in the thyroid (arithmetic means of 120 and 140 Gy for GD and TNG patients, respectively). Absorbed doses in organs other than the thyroid were much smaller, with arithmetic means of 1.6 Gy, 1.5 Gy and 0.65 Gy for esophagus, thymus and salivary glands, respectively. The arithmetic mean doses to all other organs and tissues were more than 100 times less than those to the thyroid gland. To our knowledge, this work represents the most comprehensive study to date of the doses received by persons treated with (131)I for hyperthyroidism.

Reconstruction of organ dose for external radiotherapy patients in retrospective epidemiologic studies

Organ dose estimation for retrospective epidemiological studies of late effects in radiotherapy patients involves two challenges: radiological images to represent patient anatomy are not usually available for patient cohorts who were treated years ago, and efficient dose reconstruction methods for large-scale patient cohorts are not well established. In the current study, we developed methods to reconstruct organ doses for radiotherapy patients by using a series of computational human phantoms coupled with a commercial treatment planning system (TPS) and a radiotherapy-dedicated Monte Carlo transport code, and performed illustrative dose calculations. First, we developed methods to convert the anatomy and organ contours of the pediatric and adult hybrid computational phantom series to Digital Imaging and Communications in Medicine (DICOM)-image and DICOM-structure files, respectively. The resulting DICOM files were imported to a commercial TPS for simulating radiotherapy and dose calculation for in-field organs. The conversion process was validated by comparing electron densities relative to water and organ volumes between the hybrid phantoms and the DICOM files imported in TPS, which showed agreements within 0.1 and 2%, respectively. Second, we developed a procedure to transfer DICOM-RT files generated from the TPS directly to a Monte Carlo transport code, x-ray Voxel Monte Carlo (XVMC) for more accurate dose calculations. Third, to illustrate the performance of the established methods, we simulated a whole brain treatment for the 10 year-old male phantom and a prostate treatment for the adult male phantom. Radiation doses to selected organs were calculated using the TPS and XVMC, and compared to each other. Organ average doses from the two methods matched within 7%, whereas maximum and minimum point doses differed up to 45%. The dosimetry methods and procedures established in this study will be useful for the reconstruction of organ dose to support retrospective epidemiological studies of late effects in radiotherapy patients.

Database of normalised computed tomography dose index for retrospective CT dosimetry

Volumetric computed tomography dose index (CTDIvol) is an important dose descriptor to reconstruct organ doses for patients combined with the organ dose calculated from computational human phantoms coupled with Monte Carlo transport techniques. CTDIvol can be derived from weighted CTDI (CTDIw) normalised to the tube current-time product (mGy/100 mAs), using knowledge of tube current-time product (mAs), tube potential (kVp), type of CTDI phantoms (head or body), and pitch. The normalised CTDIw is one of the characteristics of a CT scanner but not readily available from the literature. In the current study, we established a comprehensive database of normalised CTDIw values based on multiple data sources: the ImPACT dose survey from the United Kingdom, the CT-Expo dose calculation program, and surveys performed by the US Food and Drug Administration (FDA) and the National Lung Screening Trial (NLST). From the sources, the CTDIw values for a total of 68, 138, 30, and 13 scanner model groups were collected, respectively. The different scanner groups from the four data sources were sorted and merged into 162 scanner groups for eight manufacturers including General Electric (GE), Siemens, Philips, Toshiba, Elscint, Picker, Shimadzu, and Hitachi. To fill in missing CTDI values, a method based on exponential regression analysis was developed based on the existing data. Once the database was completed, two different analyses of data variability were performed. First, we averaged CTDI values for each scanner in the different data sources and analysed the variability of the average CTDI values across the different scanner models within a given manufacturer. Among the four major manufacturers, Toshiba and Philips showed the greatest coefficient of variation (COV) (=standard deviation/mean) for the head and body normalised CTDIw values, 39% and 54%, respectively. Second, the variation across the different data sources was analysed for CT scanners where more than two data sources were involved. The CTDI values for the scanners from Siemens showed the greatest variation across the data sources, being about four times greater than the variation of Toshiba scanners. The established CTDI database will be used for the reconstruction of CTDIvol and then the estimation of individualised organ doses for retrospective patient cohorts in epidemiologic studies.

Prediction of the location and size of the stomach using patient characteristics for retrospective radiation dose estimation following radiotherapy

Following cancer radiotherapy, reconstruction of doses to organs, other than the target organ, is of interest for retrospective health risk studies. Reliable estimation of doses to organs that may be partially within or fully outside the treatment field requires reliable knowledge of the location and size of the organs, e.g., the stomach, which is at risk from abdominal irradiation. The stomach location and size are known to be highly variable between individuals, but have been little studied. Moreover, for treatments conducted years ago, medical images of patients are usually not available in medical records to locate the stomach. In light of the poor information available to locate the stomach in historical dose reconstructions, the purpose of this work was to investigate the variability of stomach location and size among adult male patients and to develop prediction models for the stomach location and size using predictor variables generally available in medical records of radiotherapy patients treated in the past. To collect data on stomach size and position, we segmented the contours of the stomach and of the skeleton on contemporary computed tomography (CT) images for 30 male patients in supine position. The location and size of the stomach was found to depend on body mass index (BMI), ponderal index (PI), and age. For example, the anteroposterior dimension of the stomach was found to increase with increasing BMI (≈0.25 cm kg(-1) m(2)) whereas its craniocaudal dimension decreased with increasing PI (≈-3.3 cm kg(-1) m(3)) and its transverse dimension increased with increasing PI (≈2.5 cm kg(-1) m(3)). Using the prediction models, we generated three-dimensional computational stomach models from a deformable hybrid phantom for three patients of different BMI. Based on a typical radiotherapy treatment, we simulated radiotherapy treatments on the predicted stomach models and on the CT images of the corresponding patients. Those dose calculations demonstrated good agreement between predicted and actual stomachs compared with doses derived from a reference model of the body that might be used in the absence of individual CT scan data.

Evaluation of the use of surrogate tissues for calculating radiation dose to lymphatic nodes from external photon beams

Lymphatic node chains of the human body are particularly difficult to realistically model in computational human phantoms. In the absence of a lymphatic node model, researchers have used the following surrogate tissues to calculate the radiation dose to the lymphatic nodes: blood vessels, muscle and the combination of the muscle and adipose tissues. In the present work, the authors investigated whether and in which extent the use of different surrogate tissues is appropriate to assess the lymph node dose, using a realistic model of lymphatic nodes that the authors recently reported. Using a Monte Carlo radiation transport method coupled with the adult male hybrid phantom that included the lymph node model, the air kerma-to-absorbed dose conversion coefficients (Gy Gy(-1)) to the lymph nodes and other tissues used as surrogates for external photon beams of 15 discrete energies (0.015-10 MeV) were computed using the following six idealised geometries: anterior-posterior (AP), posterior-anterior (PA), right lateral, left lateral, rotational and isotropic. To validate the results of this study, the lymph node dose calculated here was compared with the dose published by the International Commission on Radiological Protection for the adult male reference phantom. The lymph node dose conversion coefficients with the values calculated for the blood vessels, muscle, adipose tissue and the combination of muscle and adipose tissues were then compared. It was found that muscle was the best estimator for the lymph nodes, with a dose difference averaged across energies >0.08 MeV of <8 % in all irradiation geometries excluding the AP and PA geometries for which the blood vessels were found to be the best estimator. In conclusion, muscle and blood vessels may preferably be used as surrogate tissues in the absence of lymphatic nodes in a given voxel phantom. For energies <0.08 MeV, for which the authors observed a difference of up to 30-fold, an explicit lymph node model may be required to prevent increasing differences with the lymph node dose as the photon energy decreases, though the absolute values of the dose conversion coefficients are smaller than at higher energy.

Computational lymphatic node models in pediatric and adult hybrid phantoms for radiation dosimetry

We developed models of lymphatic nodes for six pediatric and two adult hybrid computational phantoms to calculate the lymphatic node dose estimates from external and internal radiation exposures. We derived the number of lymphatic nodes from the recommendations in International Commission on Radiological Protection (ICRP) Publications 23 and 89 at 16 cluster locations for the lymphatic nodes: extrathoracic, cervical, thoracic (upper and lower), breast (left and right), mesentery (left and right), axillary (left and right), cubital (left and right), inguinal (left and right) and popliteal (left and right), for different ages (newborn, 1-, 5-, 10-, 15-year-old and adult). We modeled each lymphatic node within the voxel format of the hybrid phantoms by assuming that all nodes have identical size derived from published data except narrow cluster sites. The lymph nodes were generated by the following algorithm: (1) selection of the lymph node site among the 16 cluster sites; (2) random sampling of the location of the lymph node within a spherical space centered at the chosen cluster site; (3) creation of the sphere or ovoid of tissue representing the node based on lymphatic node characteristics defined in ICRP Publications 23 and 89. We created lymph nodes until the pre-defined number of lymphatic nodes at the selected cluster site was reached. This algorithm was applied to pediatric (newborn, 1-, 5-and 10-year-old male, and 15-year-old males) and adult male and female ICRP-compliant hybrid phantoms after voxelization. To assess the performance of our models for internal dosimetry, we calculated dose conversion coefficients, called S values, for selected organs and tissues with Iodine-131 distributed in six lymphatic node cluster sites using MCNPX2.6, a well validated Monte Carlo radiation transport code. Our analysis of the calculations indicates that the S values were significantly affected by the location of the lymph node clusters and that the values increased for smaller phantoms due to the shorter inter-organ distances compared to the bigger phantoms. By testing sensitivity of S values to random sampling and voxel resolution, we confirmed that the lymph node model is reasonably stable and consistent for different random samplings and voxel resolutions.

Risk of treatment-related esophageal cancer among breast cancer survivors

BACKGROUND

Radiotherapy for breast cancer may expose the esophagus to ionizing radiation, but no study has evaluated esophageal cancer risk after breast cancer associated with radiation dose or systemic therapy use.

DESIGN

Nested case-control study of esophageal cancer among 289 748 ≥5-year survivors of female breast cancer from five population-based cancer registries (252 cases, 488 individually matched controls), with individualized radiation dosimetry and information abstracted from medical records.

RESULTS

The largest contributors to esophageal radiation exposure were supraclavicular and internal mammary chain treatments. Esophageal cancer risk increased with increasing radiation dose to the esophageal tumor location (P(trend )< 0.001), with doses of ≥35 Gy associated with an odds ratio (OR) of 8.3 [95% confidence interval (CI) 2.7-28]. Patients with hormonal therapy ≤5 years preceding esophageal cancer diagnosis had lower risk (OR = 0.4, 95% CI 0.2-0.8). Based on few cases, alkylating agent chemotherapy did not appear to affect risk. Our data were consistent with a multiplicative effect of radiation and other esophageal cancer risk factors (e.g. smoking).

CONCLUSIONS

Esophageal cancer is a radiation dose-related complication of radiotherapy for breast cancer, but absolute risk is low. At higher esophageal doses, the risk warrants consideration in radiotherapy risk assessment and long-term follow-up.

Comparison of internal dosimetry factors for three classes of adult computational phantoms with emphasis on I-131 in the thyroid

The S values for 11 major target organs for I-131 in the thyroid were compared for three classes of adult computational human phantoms: stylized, voxel and hybrid phantoms. In addition, we compared specific absorbed fractions (SAFs) with the thyroid as a source region over a broader photon energy range than the x- and gamma-rays of I-131. The S and SAF values were calculated for the International Commission on Radiological Protection (ICRP) reference voxel phantoms and the University of Florida (UF) hybrid phantoms by using the Monte Carlo transport method, while the S and SAF values for the Oak Ridge National Laboratory (ORNL) stylized phantoms were obtained from earlier publications. Phantoms in our calculations were for adults of both genders. The 11 target organs and tissues that were selected for the comparison of S values are brain, breast, stomach wall, small intestine wall, colon wall, heart wall, pancreas, salivary glands, thyroid, lungs and active marrow for I-131 and thyroid as a source region. The comparisons showed, in general, an underestimation of S values reported for the stylized phantoms compared to the values based on the ICRP voxel and UF hybrid phantoms and relatively good agreement between the S values obtained for the ICRP and UF phantoms. Substantial differences were observed for some organs between the three types of phantoms. For example, the small intestine wall of ICRP male phantom and heart wall of ICRP female phantom showed up to eightfold and fourfold greater S values, respectively, compared to the reported values for the ORNL phantoms. UF male and female phantoms also showed significant differences compared to the ORNL phantom, 4.0-fold greater for the small intestine wall and 3.3-fold greater for the heart wall. In our method, we directly calculated the S values without using the SAFs as commonly done. Hence, we sought to confirm the differences observed in our S values by comparing the SAFs among the phantoms with the thyroid as a source region for selected target organs--small intestine wall, lungs, pancreas and breast--as well as illustrate differences in energy deposition across the energy range (12 photon energies from 0.01 to 4 MeV). Differences were found in the SAFs between phantoms in a similar manner as the differences observed in S values but with larger differences at lower photon energies. To investigate the differences observed in the S and SAF values, the chord length distributions (CLDs) were computed for the selected source--target pairs and compared across the phantoms. As demonstrated by the CLDs, we found that the differences between phantoms in those factors used in internal dosimetry were governed to a significant degree by inter-organ distances which are a function of organ shape as well as organ location.

Study of the influence of radionuclide biokinetics on the efficiency of in vivo counting using Monte Carlo simulation

To improve calibration methods of in vivo counting, our laboratory has developed a computer tool to model internal contamination and assess in vivo activity and corresponding organ absorbed doses. The aim of the recent work was to define a more realistic source based on biokinetic models. The influence of the biokinetic parameters on the in vivo counting was studied through the simulation of an acute inhalation intake of (241)Am. The tissue distribution of activity predicted by the biokinetic model was visualized. Two equivalent methods for determination of the efficiency related to the total activity distributed in the body were used. The comparison between the efficiency taking the biokinetics into account and the classically estimated efficiency quantifies the influence of the activity distribution in the body and provides conversion factors for correcting the classical efficiency to account for biokinetics.

Analytical and Monte Carlo assessment of activity and local dose after a wound contamination by activation products

The activity and local dose following a right index finger wound contamination by activation products are assessed. Measurements with a high purity germanium detector and a four positions measurement protocol enabled a better localization of the contaminant source. From the source location and detector calibration, the remaining wound activity and local absorbed dose were deduced. An analytical model, based on a two dimensional simplification of the problem, is presented. It is shown to provide a fast and quite accurate activity assessment when the contaminants are described as a point source. The contaminants' location and activity were then more accurately assessed using Monte Carlo calculations based on the OEDIPE software and a voxelized phantom of the index finger. Describing the contaminant mixture as a point source resulted in an agreement of experimental and computed data around 6% for most of the radionuclides. The total activity, due to 11 radionuclides, was estimated to be (9.5 +/- 0.4) kBq at measurement day. Since the point source is found to be less than 1 mm under the skin, the equivalent skin dose is calculated and found to be around 680 mSv in the first year after the contamination, and this value decreases to 250 mSv in the second year. The relevance of equivalent skin dose as an estimate of the sanitary risk is discussed, and it is concluded that for this case it gives the upper end estimate of the risk.

OEDIPE: a new graphical user interface for fast construction of numerical phantoms and MCNP calculations

Although great efforts have been made to improve the physical phantoms used to calibrate in vivo measurement systems, these phantoms represent a single average counting geometry and usually contain a uniform distribution of the radionuclide over the tissue substitute. As a matter of fact, significant corrections must be made to phantom-based calibration factors in order to obtain absolute calibration efficiencies applicable to a given individual. The importance of these corrections is particularly crucial when considering in vivo measurements of low energy photons emitted by radionuclides deposited in the lung such as actinides. Thus, it was desirable to develop a method for calibrating in vivo measurement systems that is more sensitive to these types of variability. Previous works have demonstrated the possibility of such a calibration using the Monte Carlo technique. Our research programme extended such investigations to the reconstruction of numerical anthropomorphic phantoms based on personal physiological data obtained by computed tomography. New procedures based on a new graphical user interface (GUI) for development of computational phantoms for Monte Carlo calculations and data analysis are being developed to take advantage of recent progress in image-processing codes. This paper presents the principal features of this new GUI. Results of calculations and comparison with experimental data are also presented and discussed in this work.

Automatic application of ICRP biokinetic models in voxel phantoms for in vivo counting and internal dose assessment

As part of the improvement of calibration techniques of in vivo counting, the Laboratory of Internal Dose Assessment of the Institute of Radiological Protection and Nuclear Safety has developed a computer tool, 'OEDIPE', to model internal contamination, to simulate in vivo counting and to calculate internal dose. The first version of this software could model sources located in a single organ. As the distribution of the contamination evolves from the time of intake according to the biokinetics of the radionuclide, a new facility has been added to the software first to allow complex heterogeneous source modelling and then to automatically integrate the distribution of the contamination in the different tissues estimated by biokinetic calculation at any time since the intake. These new developments give the opportunity to study the influence of the biokinetics on the in vivo counting, leading to a better assessment of the calibration factors and the corresponding uncertainties.