- Open Access
Doses to internal organs for various breast radiation techniques - implications on the risk of secondary cancers and cardiomyopathy
© Pignol et al; licensee BioMed Central Ltd. 2011
Received: 4 November 2010
Accepted: 14 January 2011
Published: 14 January 2011
Breast cancers are more frequently diagnosed at an early stage and currently have improved long term outcomes. Late normal tissue complications induced by adjuvant radiotherapy like secondary cancers or cardiomyopathy must now be avoided at all cost. Several new breast radiotherapy techniques have been developed and this work aims at comparing the scatter doses of internal organs for those techniques.
A CT-scan of a typical early stage left breast cancer patient was used to describe a realistic anthropomorphic phantom in the MCNP Monte Carlo code. Dose tally detectors were placed in breasts, the heart, the ipsilateral lung, and the spleen. Five irradiation techniques were simulated: whole breast radiotherapy 50 Gy in 25 fractions using physical wedge or breast IMRT, 3D-CRT partial breast radiotherapy 38.5 Gy in 10 fractions, HDR brachytherapy delivering 34 Gy in 10 treatments, or Permanent Breast 103Pd Seed Implant delivering 90 Gy.
For external beam radiotherapy the wedge compensation technique yielded the largest doses to internal organs like the spleen or the heart, respectively 2,300 mSv and 2.7 Gy. Smaller scatter dose are induced using breast IMRT, respectively 810 mSv and 1.1 Gy, or 3D-CRT partial breast irradiation, respectively 130 mSv and 0.7 Gy. Dose to the lung is also smaller for IMRT and 3D-CRT compared to the wedge technique. For multicatheter HDR brachytherapy a large dose is delivered to the heart, 3.6 Gy, the spleen receives 1,171 mSv and the lung receives 2,471 mSv. These values are 44% higher in case of a balloon catheter. In contrast, breast seeds implant is associated with low dose to most internal organs.
The present data support the use of breast IMRT or virtual wedge technique instead of physical wedges for whole breast radiotherapy. Regarding partial breast irradiation techniques, low energy source brachytherapy and external beam 3D-CRT appear safer than 192Ir HDR techniques.
Breast is the most common site of cancer in women and with the wide-spread use of mammography more than two-thirds of breast cancers are diagnosed at an early stage [1, 2]. Early stage breast cancer carries a better prognosis, with outcomes having improved dramatically over the last two decades with a 25% reduction of breast cancer mortality . As patients diagnosed with breast cancer are more likely to survive longer, it is essential to prevent treatment induced fatalities. The main types of radiation therapy induced fatalities that have been widely reported are cardiomyopathy and secondary cancers . Though their occurrence is also influenced by lifestyle and/or a predisposing genetic condition [5, 6], it is primarily related to the amount of dose deposited in specific organs [6, 7]. So the most efficient way to prevent these sequelae is to reduce the amount of dose scattered to internal organs; for example, choosing a radiation technique that minimizes the exposure of internal organs [5–7]. In regards to secondary cancers, a recent review from Xu et al. showed that secondary tumors occur more frequently in organs that are close to radiation fields, in the high/intermediate dose zones , and that it is important to assess the scattered dose to those internal organs along with their secondary cancer susceptibility in selecting a radiation technique. In regards to the cardiomyopathy risk, a critical review published by Schultz-Hector stresses the risk of acute dose as low as 1 ~ 2 Gy and a dose-dependent cardiac mortality below 10 Gy .
On the other hand, for adjuvant breast radiotherapy, several innovations and new paradigm have been introduced over the last decade. Physical wedge were replaced by virtual wedges, and eventually the dose distribution homogeneity was improved using breast Intensity Modulated Radiation Therapy (IMRT) . Multiple techniques have been proposed for Accelerated Partial Breast Irradiation (APBI) of early stage breast cancer [10, 11], which include high dose rate multi-catheter brachytherapy and permanent breast seed implant (PBSI), intra-operative radiotherapy using kilovoltage generator or direct electron beam, and 3D-conformal radiotherapy [12–17].
All of these techniques deliver different levels of scatter doses to internal organs and hence may induce different risks of secondary cancers or cardiomyopathy. The purpose of this paper is to evaluate the amount of scattered dose to internal organs situated in the intermediate/high dose region including the heart, the lung, the contralateral breast and the spleen for different techniques of adjuvant radiotherapy for a typical left sided breast cancer. To avoid confounding factors linked to patient's anatomical characteristics and assess internal organ dose deposition accurately, we used Monte Carlo simulation in an anthropomorphic phantom based on a realistic patient anatomy.
1 Radiotherapy protocols
Five different breast irradiation protocols were selected: a standard whole breast radiotherapy delivering 50 Gy in 25 treatments to the breast alone, using either physical wedge or virtual wedge/breast IMRT for missing tissue compensation [9, 18], partial breast 3D-conformal radiotherapy (3D-CRT) delivering 38.5 Gy in 10 treatments , multi-catheter High Dose Rate (HDR) brachytherapy delivering 34 Gy in 10 treatments to the 85% isodose [11, 12], and permanent breast seed implants with 103Pd seeds delivering a dose of 90 Gy on the Planning Target Volume (PTV) .
2 Realistic anthropomorphic phantom
A realistic anthropomorphic phantom of a female chest was described in the data entry card of the MCNP Monte Carlo code . This phantom mimicked the planning CT of a small breasted patient randomly selected from the treatment planning database. The geometry modeled was of a typical early stage cancer in the left breast. Complex volumes were build using elementary surfaces combination to create breasts, lungs, heart, chest walls, spleen and other body volumes. Small spherical tally volumes (0.5 to 0.8 cc) were placed in the left and right breasts, on the anterior part of the heart corresponding to the left anterior descending coronary artery , and in the posterior part of the ipsilateral lung. A larger spherical tally volume (150 cc) was placed at the position of the spleen, about 5 cm inferiorly to the breast field edge. The MCNP *F8 pulsed height tally function corrected for energy deposition was used to calculate the amount of energy absorbed in each tally volume. This function calculates for each tally the amount of energy deposited minus the energy leaving the volume. Previous work done by our group demonstrated the accuracy of this method in estimating the absorbed dose . These values were converted into dose, accounting for the energy absorbed in the treated breast and the treatment protocol. To facilitate comparison with previously published data, the doses were expressed in Gy (J kg-1) when discussing the risk of cardiomyopathy, and in mSv when discussing the risk of secondary cancers.
3 External beam radiotherapy
3.1 Hybrid method
Head leakage and room scatter contributions are challenging to assess using Monte Carlo simulation because of the very low probability for a photon to reach a detector inside the phantom. So a hybrid method was used to calculate the scatter dose for external beam radiotherapy techniques. This method adds the dose corresponding to head leakage and room back-scatter measured in a water phantom to the scatter dose produced in beam modifiers and internal phantom scatter calculated using Monte Carlo simulation.
3.2 Head leakage and room back-scatter
The head leakage and room back-scatter contributions were measured in a solid water phantom (Gammax RMI, Middleton, WI) using a Farmer ionization chamber (model 2571). The phantom was placed at a source-axis distance of 100 cm, laterally abutting the central axis of half beam irradiation fields of various sizes: 16 × 20 and 8 × 20 cm2. Doses were measured at 5 cm depth in the phantom and at 2.5, 7, 10, 19 and 28 cm away from the beam axis. These scatter doses were interpolated for each field size using a power law.
3.3 Scatter contribution
Dose contributions due to photons scattering from beam modifiers and/or inside the phantom were simulated using the MCNP Monte Carlo code . The photon energy phase-space from a Siemens Primus (Walnut Creek, CA) 6 MV accelerator was pre-calculated . Two opposed parallel beams described as being tangential to the chest wall with a 1 cm lung margin. Field sizes were 16 × 20 cm2 for whole breast irradiation, and 8 × 20 cm2 for the 3D-CRT partial breast irradiation technique. Missing tissue compensation technique used either 30°steel wedges (ρ = 7.81 g.cm-3), or field in field segments for about 20% of the dose, the remaining 80% was delivered using open beams. This was done to simulate a virtual wedge/breast IMRT technique.
4.1 Catheter 192Ir HDR brachytherapy
A photon energy spectrum with discrete energy probabilities corresponding to 192Ir decay was described in the source card. Photons were emitted in 4 π starting randomly from the source placed in the middle of the left breast. The number of photons that were generated was calculated based on the dwell time needed to treat a target volume of 3 cm radius (113 cc), corresponding to a volume of 113 cc, using a 10 Ci source. The Nucletron Plato treatment planning system (Veenendaal, Netherland) was used to calculate the total dwell time, placing catheter evenly spaced every cm across the target volume. In this later case the IPSA dose optimization algorithm was used to generate the dwell positions, to deliver the prescribed dose to the target volume and calculate the total treatment time .
4.1 Permanent breast seed implants (PBSI)
The same target volume geometry was used to simulate the PBSI case. A target volume of 113 cc requires a hundred 103Pd seeds of 2.7 U, corresponding to a total activity of 0.2088 Ci to deliver a dose of 90 Gy on the minimal peripheral dose .
5 Risk of secondary cancers estimation
The lifetime probabilities of developing fatal secondary malignancies were calculated per Sv absorbed in breast and lung using the National Council on Radiation Protection and Measurements (NCRP) report 116 Table Seven Part Two page 32 .
6 Estimation of statistical errors
A typical Monte Carlo result represents the average of the contributions from many particles histories. To calculate this average and the standard deviation the initial problem is divided in several smaller batches. A standard error, R, is then calculated as being the ratio between the standard deviation and the average: . A standard error below 5% is generally considered reliable for most calculation. For the current study, the transport of 109 photons sources was simulated for each opposed beam in order to get reliable estimation of the scattered dose, i.e. with standard error below 1%.
Relative contribution from head compensator, leakage and internal scatter to the dose to various organs.
Heart (anterior 1/3)
Dose to various organs for various breast radiotherapy techniques.
For external beam radiotherapy the physical wedge compensation technique yields the largest dose to neighboring solid organs like the spleen or the heart giving respectively 2,356 mSv and 3.0 Gy respectively. Breast IMRT reduces the dose these neighboring organs to 866 mSv and 1.4 Gy respectively, and partial breast irradiation using 3D-CRT is the safest technique with doses of 130 mSv and 0.7 Gy respectively. The dose scattered in the lung is small for IMRT and 3D-CRT, but higher for the wedge technique.
For partial left breast irradiation using 192Ir HDR brachytherapy large doses are scattered to the heart (3.6 Gy), the spleen (1,171 mSv), and the lung (2,471 mSv). Using a balloon catheter these doses are increased by 44% reaching 5.2 Gy to the heart, 1,686 mSv to the spleen and 3,558 mSv to the posterior part of the ipsilateral lung. In contrast, permanent breast seeds implant brachytherapy using low energy source is associated with low doses to most organs despite a higher physical dose is delivered to the target volume. The brachytherapy techniques tend to deliver higher dose to the lung compared to external beam techniques where shielding is used.
This report shows that depending on the radiotherapy techniques large variations, e.g. up to 20 fold for the ipsilateral lung and 800 fold for the contralateral breast, are found in the amount of scattered dose to the organs depending on the adjuvant breast radiation technique. The objective of this work was not to describe the range of scatter doses received by adjuvant breast radiotherapy, since this amount is also highly dependant on other factors including the breast size and side, the location of the surgical cavity for brachytherapy techniques, and the patient body shape and size [5, 18]. For example we previously reported up to a 10 fold variation in the dose scattered in the contralateral breast in a prospective study measuring the scatter dose to various body locations in patients receiving standard external beam radiotherapy . To evaluate the long term risks of breast radiotherapy, we compared the scattered dose produced by various radiotherapy techniques while keeping the patient geometry constant. We purposely selected a small left breasted patient to compare the amount of scattered dose for partial breast radiotherapy techniques versus standard whole breast radiotherapy in a worse case scenario.
Lifetime risk of secondary cancers for various breast radiotherapy techniques using the likelihoods from the National Council on Radiation Protection and Measurements (NCRP) Report 116 Table 7.2, page 32.
For the clinical case used in this study the incremental risk of secondary cancer breast cancer is calculated 0.34% for a whole breast technique and wedge compensators. This is likely undetectable compared to the observed frequency of contralateral breast cancer which is about 7% at 10 years and 10% at 15 years [25, 26]. For example Obedian did not find significant difference in the occurrence of contralateral breast cancer at 15 years in a retrospective series of 2,416 patients treated with breast conserving surgery and adjuvant radiotherapy or mastectomy without radiotherapy . Though this risk might be higher for younger women or patients with predisposing genetic risks [6, 25, 27], it remains difficult to detect. Moreover, compared to physical wedge compensation radiotherapy the other techniques, especially the ones delivering partial breast irradiation, yield at least 7 times less scatter dose. So the risk of developing a contralateral breast cancer should be truly undetectable.
For a whole breast technique using physical wedge compensation the lifetime incremental risk of lung cancer is calculated at 0.49%. This value is little higher but of the same order of magnitude than the 0.30% increased risk for adjuvant radiotherapy found by Zablotska on a cohort of 260,000 patients included in the Surveillance Epidemiology and End Results (SEER) database . This difference could be due to the high dose gradient in the lung, the choice of a small breasted women and the position of the detector in the ipsilateral lung that all could increase the amount of scatter dose detected. Nevertheless, from a clinical perspective those rates are small and the risk remains acceptable. Since most radiotherapy techniques except the HDR brachytherapy are yielding similar or lower amount of radiation scatter to the lung they should also be deemed acceptable. The only scenario where a large scatter dose is found in the lung is 192Ir HDR brachytherapy. This is likely due to the limited absorption in the lung tissue of the high energy photons (average energy 367 keV) that are emitted in all directions and without shielding from the 192Ir source. The risk of secondary lung cancer calculated in this case is increased by a factor 4, with 2 in 100 women at risk of developing lung cancer.
Regarding cardiac risk, a recent critical review published by Schultz-Hector suggest that acute single dose of 1~2 Gy to the heart increased the risk of developing ischemic heart disease significantly . And the excess relative risk could be linearly fitted with a slope of 17% per Gy. Bearing in mind those values, external beam radiotherapy with physical wedge compensation and HDR breast brachytherapy which yield excess dose to the heart are deemed inappropriate breast adjuvant radiotherapy techniques. Since the use of 103Pd has a strong protecting effect on the heart dose, the low energy photons being absorbed rapidly in the tissue, alternative sources like low energy electronic or 169Yb sources should be considered for HDR applications [29, 30].
Since the majority of women eligible for breast conserving therapy have improved outcomes, they are likely to live long enough to develop secondary cancers or cardiac failures and it is important to prevent those morbidities when considering a new technique. Whole breast radiotherapy, breast IMRT and virtual wedges appears safer than physical wedge compensation, and for partial breast irradiation techniques, external beam 3D-CRT and low energy source brachytherapy appear safer than 192Ir HDR techniques.
This project was made possible with the generous support from the Canadian Breast Cancer Foundation - Ontario Chapter.
- Elkin EB, Hudis C, Begg CB, Schrag D: The effect of changes in tumor size on breast carcinoma survival in the U.S.: 1975-1999. Cancer 2005, 104: 1149-1157. 10.1002/cncr.21285View ArticlePubMedGoogle Scholar
- Anderson WF, Jatoi I, Devesa SS: Assessing the impact of screening mammography: breast cancer incidence and mortality rates in Connecticut (1943-2002). Breast Cancer Res Treat 2006, 99: 333-340. 10.1007/s10549-006-9214-zView ArticlePubMedGoogle Scholar
- Peto R, Boreham J, Clarke M, Davies C, Beral V: UK and USA breast cancer deaths down 25% in year 2000 at ages 20-69 years. Lancet 2000, 355: 1822. 10.1016/S0140-6736(00)02277-7View ArticlePubMedGoogle Scholar
- Darby SC, McGale P, Taylor CW, Peto R: Long-term mortality from heart disease and lung cancer after radiotherapy for early breast cancer: prospective cohort study of about 300,000 women in US SEER cancer registries. Lancet Oncol 2005, 6: 557-65. 10.1016/S1470-2045(05)70251-5View ArticlePubMedGoogle Scholar
- Hall EJ: Intensity-modulated radiation therapy, protons, and the risk of second cancers. Int J Radiat Oncol Biol Phys 2006, 65: 1-7. 10.1016/j.ijrobp.2006.01.027View ArticlePubMedGoogle Scholar
- Tubiana M: Can we reduce the incidence of second primary malignancies occurring after radiotherapy? A critical review. Radiother Oncol 2009, 91: 4-15. 10.1016/j.radonc.2008.12.016View ArticlePubMedGoogle Scholar
- Xu XG, Bednarz B, Paganetti H: A review of dosimetry studies on external-beam radiation treatment with respect to second cancer induction. Phys Med Biol 2008, 53: R193-241. 10.1088/0031-9155/53/13/R01PubMed CentralView ArticlePubMedGoogle Scholar
- Schultz-Hector S, Trott K: Radiation-induced cardiovascular diseases: Is the epidemiologic evidence compatible with the radiobiologic data? Int J Radiat Oncol Biol Phys 2007, 67: 10-8. 10.1016/j.ijrobp.2006.08.071View ArticlePubMedGoogle Scholar
- Pignol JP, Olivoto I, Rakovitch E, Gardner S, Sixel K, Beckham W, Vu TT, Truong P, Ackerman I, Paszat L: A Phase III multicentre clinical trial of Breast Intensity Modulated Radiation Therapy (IMRT) to reduce skin radiation side effects. J Clin Oncol 2008, 26: 2085-2092. 10.1200/JCO.2007.15.2488View ArticlePubMedGoogle Scholar
- Bethune WA: Partial breast irradiation for early stage breast cancer. J Natl Med Assoc 1991, 83: 768.PubMed CentralPubMedGoogle Scholar
- Vicini FA, Kestin L, Chen P, Benitez P, Goldstein NS, Martinez : Limited-field radiation therapy in the management of early stage breast cancer. J Natl Canc Inst 2003, 95: 1205-1211. 10.1093/jnci/djg023View ArticleGoogle Scholar
- Wazer DE, Berle L, Graham R, Chung M, Rothschild J, Graves T, Cady B, Ulin K, Ruthazer R, DiPetrillo TA: Preliminary results of a Phase I/II study of HDR sbrachytherapy alone for T1/T2 breast cancer. Int J Radiat Oncol Biol Phys 2002, 53: 889-897. 10.1016/S0360-3016(02)02824-9View ArticlePubMedGoogle Scholar
- Keisch M, Vicini F, Kuske RR, Hebert M, White J, Quiet C, Arthur D, Scroggins T, Streeter O: Initial clinical experience with the MammoSite breast brachytherapy applicator in women with early-stage breast cancer treated with breast-conserving therapy. Int J Radiat Oncol Biol Phys 2003, 55: 289-93. 10.1016/S0360-3016(02)04277-3View ArticlePubMedGoogle Scholar
- Pignol JP, Rakovitch E, Keller B, Sankreacha R, Chartier C: Tolerance and acceptance results of a palladium-103 permanent breast seed implant Phase I/II study. Int J Radiat Oncol Biol Phys 2009, 73: 1482-1488. 10.1016/j.ijrobp.2008.06.1945View ArticlePubMedGoogle Scholar
- Vaidya JS, Joseph DJ, Tobias JS, Bulsara M, Wenz F, Saunders C, Alvarado M, Flyger HL, Massarut S, Eiermann W, Keshtgar M, Dewar J, Kraus-Tiefenbacher U, Sütterlin M, Esserman L, Holtveg HM, Roncadin M, Pigorsch S, Metaxas M, Falzon M, Matthews A, Corica T, Williams NR, Baum M: Targeted intraoperative radiotherapy versus whole breast radiotherapy for breast cancer (TARGIT-A trial): an international, prospective, randomised, non-inferiority phase 3 trial. Lancet 2010, 376: 91-102. 10.1016/S0140-6736(10)60837-9View ArticlePubMedGoogle Scholar
- Veronesi U, Orecchia R, Luini A, Galimberti V, Gatti G, Intra M, Veronesi P, Leonardi MC, Ciocca M, Lazzari R, Caldarella P, Simsek S, Silva LS, Sances D: Full-dose intraoperative radiotherapy with electrons during breast-conserving surgery: experience with 590 cases. Ann Surg 2005, 242: 101-6. 10.1097/01.sla.0000167927.82353.bcPubMed CentralView ArticlePubMedGoogle Scholar
- Vicini FA, Remouchamps V, Wallace M, Sharpe M, Fayad J, Tyburski L, Letts N, Kestin L, Edmundson G, Pettinga J, Goldstein NS, Wong J: Ongoing clinical experience utilizing 3D conformal external beam radiotherapy to deliver partial-breast irradiation in patients with early stage breast cancer treated with breast conserving surgery. Int J Radiat Oncol Biol Phys 2003, 57: 1247-1253. 10.1016/S0360-3016(03)01573-6View ArticlePubMedGoogle Scholar
- Woo TC, Pignol JP, Rakovitch E, Vu T, Hicks D, O'Brien P, Pritchard K: A Prospective Study of Scattered Radiation During Breast Radiotherapy. Int J Radiat Oncol Biol Phys 2006, 65: 52-58. 10.1016/j.ijrobp.2005.11.023View ArticlePubMedGoogle Scholar
- Briesmeister JF: MCNP - A General Monte Carlo N-Particle Transport Code, Version 5. Los Alamos National Laboratory: Los Alamos; 2003.Google Scholar
- Fuller SA, Haybittle JL, Smith RE, Dobbs HJ: Cardiac doses in post-operative breast irradiation. Radiother Oncol 1992, 25: 19-24. 10.1016/0167-8140(92)90190-6View ArticlePubMedGoogle Scholar
- Keller B, Beachey D, Pignol JP: Experimental Measurement of Radiological Penumbra Associated with Intermediate Energy X-Rays (1 MV): Implications for Small Field Stereotactic Radiosurgery. Med Phys 2007, 34: 3996-4002. 10.1118/1.2775666View ArticlePubMedGoogle Scholar
- Pignol JP, Keller B: Electron and photon spread contributions to the radiological penumbra for small monoenergetic x-ray beam (2 MeV). J Appl Phys 2009, 105: 111020. 10.1063/1.3115667View ArticleGoogle Scholar
- Lessard E, Pouliot J: Inverse planning anatomy-based dose optimization for HDR-brachytherapy of the prostate using fast simulated annealing algorithm and dedicated objective function. Med Phys 2001, 28: 773-9. 10.1118/1.1368127View ArticlePubMedGoogle Scholar
- National Council on Radiation Protection and Measurements (NCRP) report 116[http://www.ncrponline.org/Publications/116press.html]
- Fowble B, Hanlon A, Freedman G, Nicolaou N, Anderson P: Second cancers after conservative surgery and radiation for stages I-II breast cancer: identifying a subset of women at increased risk. Int J Radiat Oncol Biol Phys 2001, 51: 679-90. 10.1016/S0360-3016(01)01665-0View ArticlePubMedGoogle Scholar
- Obedian E, Fischer DB, Haffty BG: Second malignancies after treatment of early-stage breast cancer: lumpectomy and radiation therapy versus mastectomy. J Clin Oncol 2000, 18: 2406-12.PubMedGoogle Scholar
- Stovall M, Smith SA, Langholz BM, Boice JD Jr, Shore RE, Andersson M, Buchholz TA, Capanu M, Bernstein L, Lynch CF, Malone KE, Anton-Culver H, Haile RW, Rosenstein BS, Reiner AS, Thomas DC, Bernstein JL: Dose to the contralateral breast from radiotherapy and risk of second primary breast cancer in the WECARE study. Int J Radiat Oncol Biol Phys 2008, 72: 1021-30. 10.1016/j.ijrobp.2008.02.040PubMed CentralView ArticlePubMedGoogle Scholar
- Zablotska LB, Neugut AI: Lung carcinoma after radiation therapy in women treated with lumpectomy or mastectomy for primary breast carcinoma. Cancer 2003, 97: 1404-11. 10.1002/cncr.11214View ArticlePubMedGoogle Scholar
- Rivard MJ, Davis SD, DeWerd LA, Rusch TW, Axelrod S: Calculated and measured brachytherapy dosimetry parameters in water for the Xoft Axxent X-Ray Source: an electronic brachytherapy source. Med Phys 2006, 33: 4020-32. 10.1118/1.2357021View ArticlePubMedGoogle Scholar
- Lymperopoulou G, Papagiannis P, Angelopoulos A, Karaiskos P, Georgiou E, Baltas D: A dosimetric comparison of 169 Yb and 192 Ir for HDR brachytherapy of the breast, accounting for the effect of finite patient dimensions and tissue inhomogeneities. Med Phys 2006, 33: 4583-9. 10.1118/1.2392408View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.