The influence of MRI scan position on patients with oropharyngeal cancer undergoing radical radiotherapy
© Hanvey et al.; licensee BioMed Central Ltd. 2013
Received: 19 December 2012
Accepted: 24 May 2013
Published: 28 May 2013
The purpose of this study was to demonstrate how magnetic resonance imaging (MRI) patient position protocols influence registration quality in patients with oropharyngeal cancer undergoing radical radiotherapy and the consequences for gross tumour volume (GTV) definition and radiotherapy planning.
Methods and materials
Twenty-two oropharyngeal patients underwent a computed tomography (CT), a diagnostic MRI (MRID) and an MRI in the radiotherapy position within an immobilization mask (MRIRT). Clinicians delineated the GTV on the CT viewing the MRID separately (GTVC); on the CT registered to MRID (GTVD) and on the CT registered to MRIRT (GTVRT). Planning target volumes (PTVs) were denoted similarly. Registration quality was assessed by measuring disparity between structures in the three set-ups. Volumetric modulated arc therapy (VMAT) radiotherapy planning was performed for PTVC, PTVD and PTVRT. To determine the dose received by the reference PTVRT, we optimized for PTVC and PTVD while calculating the dose to PTVRT. Statistical significance was determined using the two-tailed Mann–Whitney or two-tailed paired student t-tests.
A significant improvement in registration accuracy was found between CT and MRIRT versus the MRID measuring distances from the centre of structures (geometric mean error of 2.2 mm versus 6.6 mm). The mean GTVC (44.1 cm3) was significantly larger than GTVD (33.7 cm3, p value = 0.027) or GTVRT (30.5 cm3, p value = 0.014). When optimizing the VMAT plans for PTVC and investigating the mean dose to PTVRT neither the dose to 99% (58.8%) nor 95% of the PTV (84.7%) were found to meet the required clinical dose constraints of 90% and 95% respectively. Similarly, when optimizing for PTVD the mean dose to PTVRT did not meet clinical dose constraints for 99% (14.9%) nor 95% of the PTV (66.2%). Only by optimizing for PTVRT were all clinical dose constraints achieved.
When oropharyngeal patients MRI scans are performed in the radiotherapy position there are significant improvements in CT-MR image registration, target definition and PTV dose coverage.
KeywordsMRI CT Oropharyngeal Radiotherapy Image registration Patient setup
With intensity modulated radiotherapy and volumetric modulated arc therapy (VMAT), it is possible to deliver high doses of radiation to irregular target volumes whilst sparing normal tissue, which can result in reduced severity of radiation toxicities . Increased dose delivery and dose conformity has led to a greater significance on the accurate localization of the gross tumour volume (GTV) and neighbouring structures.
In radiotherapy (RT) planning, computed tomography (CT) remains the first choice since it provides accurate dosimetric information. Imaging modalities such as magnetic resonance imaging (MRI) and positron emission tomography (PET) may present advantages compared to CT in terms of target definition. In particular, MRI offers improved soft tissue contrast [2, 3] and reduced artefacts from dental amalgam; hence, MRI is the imaging modality of choice for oropharyngeal cancers . MRI improves target definition for patients with head and neck [5–7], prostate [8, 9] and brain [10, 11] cancer.
Challenges with MRI in clinical practice include: geometric distortion [12, 13]; motion related artefacts  and magnetic susceptibility [15, 16]. Patient positioning must also be considered, since differing patient set-ups between modalities can result in uncertainties in the location and magnitude of the GTV and thus the dose received by the patient .
Patients receiving RT to the brain can undergo an MRI in the RT position within an immobilization mask, without loss of image quality over standard imaging methods . An approximate image registration uncertainty of 2 mm occurs when registering CT image sets in the treatment position to a diagnostic MRI for patients with brain cancer  and therefore many RT centres routinely register these scans for RT planning. While it could be anticipated that imaging patients with oropharyngeal cancer in the treatment position would result in improved registration and target volume definition, to our knowledge there is no published evidence to suggest that using the diagnostic MRI provides results to the contrary. In addition, no publications have analyzed the dosimetric consequences of patient positioning during MRI scan acquisition.
The purpose of this study was to compare an MRI acquired with a diagnostic patient position on a standard MRI table to an MRI acquired on a flat table with custom immobilization and to determine how this affects CT-MRI registration, GTV definition, the resulting VMAT RT plans and the magnitude of the associated geometric and dosimetric errors. The results obtained would then provide an answer to whether it was necessary to obtain a planning MRI scan for patients with oropharyngeal cancer in the RT treatment position using an immobilization mask or whether a diagnostic patient position protocol would suffice.
Patient group and study overview
Twenty two patients with oropharyngeal cancer (age 37–72), being worked up for radical RT (patients being treated with curative intent), were identified for the study, regardless of tumour or nodal stage. The study protocol was approved by Local Ethics Committee (West of Scotland Research Ethics Service, Western Infirmary, Glasgow G11 6NT, Scotland, UK) and informed written consent was obtained from all patients. Registration quality assessment was conducted on all patients, however, in three patients a GTV evaluation was not possible, since two patients had undergone primary surgical resection and another had received induction chemotherapy, resulting in complete response.
CT and MRI scanning protocol
Patients were scanned on a GE Light-speed RT 16 slice CT scanner (GE Healthcare, WI, USA), using the current clinical scanning protocol, and immobilized within a full face and neck five point fixation thermoplastic beam directional shell (BDS) with appropriate head rest (CIVCO Medical Solutions, IA, USA). The scan extent was from superior orbital ridge to carina. A helical scan was acquired with a detector configuration of 16 × 1.25, pitch 0.938, matrix 512 × 512 and speed 18.75 mm/rot with a slice thickness of 2.5 mm.
MRI imaging parameters
Slice thick (mm)
Bandwidth (± kHz)
2D Driven-equilibrium FSE (T2-weighted)
512 × 256
3D Spoiled gradient echo (T1-weighted post-contrast)
256 × 256
The scan was repeated for the MRID on a standard curved diagnostic table without BDS using a 16-channel head, neck and spine coil, with the CT scan extent. The 16-channel head, neck and spine coil was not used for the MRIRT scan since it is not compatible with the BDS. A previous investigation using test objects has shown the image quality obtained using the 4-channel cardiac coil is of diagnostic quality .
The treatment planning system (Eclipse, Varian Medical Systems, Inc. CA, USA, version 10.0.39) was used to register the MRIRT and MRID T1 and T2-weighted scans to the CT data. This software enabled fully automated, mutual-information based rigid registration to be performed. While deformable registration is available for Eclipse, it does not allow inter-modality registration. Since commercial inter-modality deformable image registration is generally not available, nor employed clinically in our centre at present, this was not subject to further investigation. The volume over which the registration was performed was centred over the oropharyngeal region and chosen to include as much information for the registration algorithm as possible . Fully automated registration was used without manual adjustment since, for this site, satisfactory registration for both target and nodes is not always achievable and would introduce subjective errors.
To quantitatively assess the quality of registration the orbits and the odontoid process were outlined on the CT and T1-weighted MR image sets. These three structures were delineated on the CT, for each patient, and then repeated on the MR datasets registered to the CT.
where CT ∩ MRI is the volume of intersection between the CT and MRI structures. The value of the spatial overlap can range from zero, which indicates no spatial overlap between the CT and MRI volumes, to one, which indicates complete overlap . Since structures outlined on CT and MRI may differ, even with perfect registration, a spatial overlap of one may not be achieved in practice, but will still be dependent on different patient set-up and registration quality.
Gross tumour volume, lower risk clinical target volume and organ at risk delineation
Three Oncologists were assigned five patients each and a fourth was assigned four patients. The Oncologists delineated the GTV on their patients using the three set-ups PPC, PPD and PPRT, which in this study are referred to as GTVC, GTVD and GTVRT respectively. Anonymized information sheets, containing the patient’s clinical history and radiology report, were available to the Oncologists. The clinicians generally utilized the T2-weighted MR datasets while contouring the GTV, although T1-weighted images were also referenced. A period of at least a week was given between delineations of the GTV for the same patient using a different imaging protocol and the Oncologists were blinded to previous delineations. Changes in the magnitude of the GTV were assessed. Contouring was also performed on the lower risk clinical target volume (CTV LR) of the nodal areas at risk of microscopic involvement for a randomly selected cohort of ten patients, in each of the set-ups. Nodal delineation was performed according to international consensus guidelines . For these patients, the organs at risk (OARs), which included the left and right parotids, larynx, spinal cord and the brainstem, were also contoured.
The GTV and CTV LR were expanded to obtain planning target volumes, (PTV) and lower risk planning target volumes (PTV LR) respectively. To create the clinical target volume (CTV) the GTV was expanded by 1 cm isotropically, removing any overlap with bone or air cavities. The CTV was then enlarged by 3 mm isotropically and cropped from the external outline of the body to create the PTV. It is necessary to crop the PTV from the body outline to assist in the VMAT optimization process. The PTVs were generated from the GTVs using the three set-ups PPC, PPD and PPRT, which are denoted as PTVC, PTVD and PTVRT respectively. Similarly, to generate the PTV LR the CTV LR was expanded by 3 mm isotropically and cropped by an appropriate margin from the body outline. Target volumes and OARs were generated in accordance with local protocol.
VMAT plans were calculated for the ten patients for whom the GTV, CTV LR and OARs were contoured, to determine the impact that changes in target volume definition have on RT planning. The VMAT plans were calculated in accordance with our centre’s clinical dose constraints using the Anisotropic Analytical Algorithm and Progressive Resolution Optimizer VMAT algorithm in Eclipse (Varian Medical Systems, Inc, CA, USA) version 10.0.28. A VMAT plan was calculated for PTVC, PTVD and PTVRT. Dose volume histograms (DVHs) were generated for PTVC, PTVD and PTVRT. A mean DVH for all ten patients was then calculated for PPC, PPD and PPRT. These results are presented using our centre’s dose constraint protocol i.e. D99 > 90%, D95 > 95%, D5 < 105% and D2 < 107%, where D99 > 90% means 99% of the total PTV volume should receive a dose > 90% of the prescribed dose. The other dose constraints are defined similarly.
Of particular interest in this study was to establish the quality of the RT plan for PPC and PPD with reference to PPRT, since our working hypothesis postulates that the optimum target volume definition is achieved with PPRT. To achieve this we optimized for PTVC and PTVD but investigated the dose coverage of PTVRT at D99 > 90%, D95 > 95%, D5 < 105% and D2 < 107%.
A quantitative comparison of the homogeneity of the dose to the PTV was completed using the sigma index. The sigma index was compared individually for PTVC, PTVD and PTVRT as well as for PTVC, PTVD with reference to PTVRT. The sigma index is equal to the standard deviation of the dose throughout the PTV, calculated on a voxel by voxel basis , thus the higher the sigma index, the greater the dose inhomogeneity.
Finally, the dose to the OARs was assessed by evaluating the mean and maximum dose received using PPC, PPD and PPRT.
Two-tailed paired student t-tests were performed to examine the statistical differences of the registration quality and dosimetric indices, except for the geometric mean distance from the centre of structures in CT and MRI where two-tailed Mann–Whitney tests were performed. A Mann–Whitney test was used to assess the geometric mean distance from the centre of the structures because the data was not normally distributed. The null hypothesis was rejected when the p value was less than 0.05.
Mean spatial overlap of the anatomical landmarks for the two registration set-ups
CT-MRID* mean spatial overlap
CT-MRIRT† mean spatial overlap
Mean GTV* (cm 3 ) delineated with the different patient position protocols
Mean GTVD†† (p value)
Mean GTVRT#(p value)
The justification for choosing only ten patients to perform the dose analysis was that none of the ten patients met the 90% dose constraint for PTVC: PTVRT and PTVD: PTVRT and only one of the ten patients met the 95% dose constraint for both PTVC: PTVRT and PTVD: PTVRT. This relates to a 95% confidence interval of 0.0% to 30.9% and 0.3% to 44.5% respectively. These confidence intervals demonstrate that no further cases are statistically necessary in the analysis.
Mean sigma indices for the PTVs
Mean sigma index (%)
Mean dose for the (a) parotid and (b) larynx
Parotid mean dose (cGy)
Mean dose constraint 2400 cGy
Number of patients which exceed the dose constraint
Larynx mean dose (cGy)
Mean dose constraint 4000 cGy
PP C : PP C
PP D : PP D
PP RT : PP RT
PP C : PP RT
PP D : PP RT
Number of patients which exceed the dose constraint
When the mean dose to the parotid and larynx for all 10 patients was calculated, an incrementally smaller value was seen for PTVC, PTVD and PTVRT (first 3 columns of Table 5), which can be explained by the decrease in the magnitude of the GTV. Dose to the spinal cord and brainstem were also calculated but were not found to exceed clinically relevant tolerances for any of the patient position protocols.
Advanced imaging techniques have been shown to improve tumour and nodal staging  and the benefits of integrating MRI are well known [5–7, 23, 24]. While the advantages in positioning patients in a similar way to their CT planning scan when acquiring MR images are also known [6, 17, 23], this is the first study to compare the CT-MRI registration accuracy and dosimetric effects of a diagnostic versus RT positioned MRI scan in patients with oropharyngeal cancer.
During the process of registering CT to MR images, there was improved registration with fewer discrepancies using MRIRT compared to MRID, as demonstrated in Figure 1. In this example, there are registration discrepancies at the body outline and spinal cord for CT-MRID, as indicated by the arrows (Figure 1(a)), but not with MRIRT (Figure 1(b)).
A significant improvement in the registration quality of CT to MRIRT versus the MRID was demonstrated by a reduction in the geometric mean distance from the centre of the orbits and odontoid process delineated on CT and the MRIRT and an increase in the spatial overlap of these structures. These results show that patient setup significantly influences CT-MRI registration accuracy. With increased interest in the use of dose escalation and dose painting techniques within RT planning the importance of improved image registration becomes ever more relevant.
The significant difference between the magnitude of the GTVC and GTVD and between the GTVC and GTVRT implies that the GTV is significantly smaller when using registered rather than unregistered CT-MR images. This underlines the importance of registering CT to MRI for patients with oropharyngeal cancer, rather than viewing them separately. It also highlights the difficulty in delineating oropharyngeal cancers with CT due to the similarity in Hounsfield Units of tumour and surrounding tissue as well as artefacts caused by dental amalgam. While there was no significant difference between the mean GTVD and GTVRT, there were important differences in the VMAT plans, as discussed below. To achieve the clinical goal of reduced late toxicities and improved tumour control using dose escalation with tighter PTV margins, uncertainties in GTV delineation need to be minimized and our data suggest PPRT offers the optimal of the three set-ups.
The results reveal that there are potentially clinically relevant improvements to the quality of the VMAT plans when using PPRT rather than PPC or PPD. This is demonstrated by the PTV dose coverage, PTV dose homogeneity and instances of unmet dose constraints by the OARs. Due to improved registration accuracy and MRI being the recommended imaging modality for soft tissue oropharyngeal cancers , it may be assumed that PTVRT would be the reference PTV. To determine the dose received by the reference, PTVRT, we optimized for PTVC and PTVD while calculating the dose to PTVRT. When optimizing for PTVRT it was shown that the mean DVH for PTVD had poorer dose coverage than PTVC (Figure 4). Despite the magnitude of the mean GTVRT and GTVD being similar there may be differences in the shape and location of the GTV using these patient set-ups which would result in the DVHs of Figure 4. Neither the PPC nor the PPD VMAT plans were able to meet the clinical dose constraints of D99 and D95 for the reference, PTVRT. Furthermore, it has been argued that tumour control probability can be considerably compromised by an inhomogeneous dose to the PTV . It is therefore suboptimal to use PPC or PPD rather than PPRT for RT planning of patients with oropharyngeal cancer.
Obtaining a further MRI in the RT position, rather than using the original diagnostic MRI, may place greater demands on increasingly stretched healthcare resources. However, this must be weighed against the potential advantages of improved image registration and, by consequence, superior target volume definition and dose coverage of the PTV, as these results have demonstrated. Our study suggests further research, particularly in correlating dosimetric investigations with clinical outcome data, would be warranted.
When MRI scans are performed in the RT position, as opposed to using diagnostic MR images not obtained in RT position, there are significant improvements in the quality of CT-MR registration. This study has also shown that RT positioned MRI scans offer improvements in target definition, dose coverage and dose homogeneity, which could have significant implications for tumour control rates. To our knowledge, this is the first study in the literature to confirm these advantages.
Beam directional shell
Clinical target volume
- CTV LR:
Lower risk clinical target volume
Dose volume histograms
Gross tumour volume
Gross tumour volume delineated using PPC
Gross tumour volume delineated using PPD
Gross tumour volume delineated using PPRT
Magnetic resonance imaging
MRI with the patient in the standard diagnostic position
MRI with the patient in the radiotherapy position in an immobilization mask
Organ at risk
Positron emission tomography
Patient position protocol with the CT and MRID on a separate console
Patient position protocol with the CT and MRID registered
Patient position protocol with the CT and MRIRT registered
Planning target volume
LR Lower risk planning target volume
Planning target volume generated from GTVC
Gross tumour volume generated from GTVD
Gross tumour volume generated from GTVRT
Volumetric modulated arc therapy.
The authors would like to express their gratitude to the Beatson radiotherapy and diagnostic radiography staff and to Mr James Paul, Head of Biostatistics at the Clinical Trials Unit, Beatson West of Scotland Cancer Centre. We would also like to acknowledge financial support from the NHS Learning and Education and The Beatson Oncology Centre Fund.
- Staffurth JN: A review of clinical evidence for intensity-modulated radiotherapy. Clin Oncol 2010, 22: 643-657. 10.1016/j.clon.2010.06.013View Article
- Lam P, Au-Yeung KM, Cheng PW, Wei WI, Yuen AP, Trendell-Smith N, Li JH, Li R: Correlating MRI and histologic tumor thickness in the assessment of oral tongue cancer. Am J Roentgenol 2004, 182: 803-808. 10.2214/ajr.182.3.1820803View Article
- Tetsumura A, Yoshino N, Amagasa T, Nagumo K, Okada N, Sasaki T: High-resolution magnetic resonance imaging of squamous cell carcinoma of the tongue: an in vitro study. Dentomaxillofac Radiol 2001, 30: 14-21. 10.1038/sj.dmfr.4600565View ArticlePubMed
- Roland NJ, Paleri V: Head and Neck Cancer: Multidisciplinary Management Guidelines. 4th edition. London: ENT UK; 2011.
- Chung NN, Ting LL, Hsu WC, Lui LT, Wang PM: Impact of magnetic resonance imaging versus CT on nasopharyngeal carcinoma: primary tumor target delineation for radiotherapy. Head Neck 2004, 26: 241-246. 10.1002/hed.10378View ArticlePubMed
- Ahmed M, Schmidt M, Sohaib A, Kong C, Burke K, Richardson C, Usher M, Brennan S, Riddell A, Davies M, Newbold K, Harrington KJ, Nutting CM: The value of magnetic resonance imaging in target volume delineation of base of tongue tumours – a study using flexible surface coils. Radiother Oncol 2010, 94: 161-167. 10.1016/j.radonc.2009.12.021View ArticlePubMed
- Newbold K, Partidge M, Cook G, Sohaib SA, Charles-Edwards E, Rhys-Evans P, Harrington K, Nutting C: Advanced imaging applied to radiotherapy planning in head and neck cancer: a clinical review. Br J Radiol 2006, 79: 554-561. 10.1259/bjr/48822193View ArticlePubMed
- Sefrova J, Odrazka K, Paluska P, Belobradek Z, Brodak M, Dolezel M, Prosvic P, Macingova Z, Vosmik M, Hoffmann P, Louda M, Nejedla A: Magnetic resonance imaging in postprostatectomy radiotherapy planning. Int J Radiat Oncol Biol Phys 2012, 82: 911-918. 10.1016/j.ijrobp.2010.11.004View ArticlePubMed
- Smith WL, Lewis C, Bauman G, Rodrigues G, D’Souza D, Ash R, Ho D, Venkatesan V, Downey D, Fenster A: Prostate volume contouring: a 3D analysis of segmentation using 3DTRUS, CT, and MR. Int J Radiat Oncol Biol Phys 2007, 67: 1238-1247. 10.1016/j.ijrobp.2006.11.027View ArticlePubMed
- Khoo VS, Adams EJ, Saran F, Bedford JL, Perks JR, Warrington AP, Brada M: A comparison of clinical target volumes determined by CT and MRI for the radiotherapy planning of base of skull meningiomas. Int J Radiat Oncol Biol Phys 2000, 46: 1309-1317. 10.1016/S0360-3016(99)00541-6View ArticlePubMed
- Prabhakara R, Haresh KP, Ganesh T, Joshi RC, Julka PK, Rath GK: Comparison of computed tomography and magnetic resonance based target volume in brain tumors. J Canc Res Ther 2007, 3: 121-123. 10.4103/0973-1482.34694View Article
- Doran SJ, Charles-Edwards L, Reinsberg SA, Leach MO: A complete distortion correction for MR images: I. Gradient warp correction. Phys Med Biol 2005, 50: 1343-1361. 10.1088/0031-9155/50/7/001View ArticlePubMed
- Khoo VS, Joon DL: New developments in MRI for target volume delineation in radiotherapy. Br J Radiol 2006, 79: S2-S15. 10.1259/bjr/41321492View ArticlePubMed
- Tamhane1 AA, Arfanakis K: Motion correction in PROPELLER and Turboprop-MRI. Magn Reson Med 2009, 62: 174-182. 10.1002/mrm.22004View Article
- Fransson A, Andreo P, Pötter R: Aspects of MR image distortions in radiotherapy treatment planning. Strahlenther Onkol 2001, 177: 59-73. 10.1007/PL00002385View ArticlePubMed
- Schenck JF: The role of magnetic susceptibility in magnetic resonance imaging: MRI magnetic compatibility of the first and second kinds. Med Phys 1996, 23: 815-850. 10.1118/1.597854View ArticlePubMed
- Hanvey S, Sadozye AH, McJury M, Glegg M, Foster J: The influence of MRI scan position on image registration accuracy, target delineation and calculated dose in prostatic radiotherapy. Br J Radiol 2012, 85: e1256-e1262. 10.1259/bjr/26802977View ArticlePubMedPubMed Central
- Hanvey S, Glegg M, Foster J: Magnetic resonance imaging for radiotherapy planning of brain cancer patients using immobilization and surface coils. Phys Med Biol 2009, 54: 5381-5394. 10.1088/0031-9155/54/18/002View ArticlePubMed
- Kenneth U, Urie MM, Cherlow JM: Results of a multi-institutional benchmark test for cranial CT/MR image registration. Int J Radiat Oncol Biol Phys 2010, 77: 1584-1589. 10.1016/j.ijrobp.2009.10.017View Article
- Zou KH, Warfield SK, Bharatha A, Tempany CM, Kaus MR, Haker SJ, Wells WM 3rd, Jolesz FA, Kinis R: Statistical validation of image segmentation quality based on a spatial overlap index. Acad Radiol 2004, 11: 178-178. 10.1016/S1076-6332(03)00671-8View ArticlePubMedPubMed Central
- Grégoire V, Levendag P, Ang KK, Bernier J, Braaksma M, Budach V, Chao C, Coche E, Cooper JS, Cosnard G, Eisbruch A, El-Sayed S, Emami B, Grau C, Hamoir M, Lee N, Maingon P, Muller K, Reychler H: CT-based delineation of lymph node levels and related CTVs in the node-negative neck: DAHANCA, EORTC, GORTEC, NCIC, RTOG consensus guidelines. Radiother Oncol 2003, 69: 227-236. 10.1016/j.radonc.2003.09.011View ArticlePubMed
- Yoon M, Park SY, Shin D, Lee SB, Pyo HR, Kim DY, Cho KH: A new homogeneity index based on statistical analysis of the dose–volume histogram. J Appl Clin Med Phys 2007, 8: 9-17.PubMed
- Prestwich RJD, Sykes J, Carey B, Sen M, Dyker KE, Scarsbrook AF: Improving target definition for head and neck radiotherapy: a place for magnetic resonance imaging and 18-fluoride fluorodeoxyglucose positron emission tomography? Clin Oncol 2012, 24: 577-589. 10.1016/j.clon.2012.04.002View Article
- Rosenman JG, Miller EP, Tracton G, Cullip TJ: Image registration: an essential part of radiation therapy planning. Int J Radiat Oncol Biol Phys 1998, 40: 197-205. 10.1016/S0360-3016(97)00546-4View ArticlePubMed
- Goitein M, Niemierko A: Intensity-modulated therapy and inhomogeneous dose to the tumor: a note of caution. Int J Radiat Oncol Biol Phys 1996, 36: 519-522. 10.1016/S0360-3016(96)00348-3View ArticlePubMed
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.