Research Article | Open Access
Shear Wave Elastography: A New Noninvasive Tool to Assess the Intensity of Fibrosis of Irradiated Salivary Glands in Head and Neck Cancer Patients
The aim of the study was to assess salivary gland parenchyma by means of sonoelastography in patients irradiated for head and neck squamous cell carcinoma (HNSCC). The studied group consisted of 52 patients after radiotherapy (RT) and 54 healthy volunteers. All of the former were treated for advanced larynx (40), oropharynx (9), or maxilla (3) squamous cancers and suffered from chronic dryness. Ultrasonography (US) and elastography (ES) were performed, as well as an assessment of the amount of saliva and Common Terminology Criteria for Adverse Events (CTCAE) scale. There was a statistical difference between ES values in the RT group and in the controls for parotid glands (41.7 kPa versus 26.03 kPa, ) and for submandibular glands (37.6 kPa versus 22.4 kPa; ). There was a significant correlation between the CTCAE scores and objective saliva amount (), and the median amount of saliva in the examined group was lower than in the reference group (1.86 g versus 2.75 g, ). In conclusion sonoelastography adds a new parameter to ultrasonography in “one touch examination” and may be a useful tool for major salivary gland evaluation during the radiotherapy course and follow-up period.
Radiotherapy (RT) is the mainstay in head and neck squamous cell cancer (HNSCC) treatment. Radiation-induced salivary gland injuries are very common and occur in up to 80% of cases [1, 2]. Two aspects of this problem can be listed: profound salivary gland hypofunction (diminished salivary flow) and xerostomia (subjective sensation of a dry mouth). Since the National Institutes of Health Consensus Development Conference on the Oral Complications of Cancer Therapies (1989) salivary gland hypofunction and xerostomia have been considered to be basic sequelae of radiotherapy. What is important is that reduced gland function and parallel oral health-related quality of life decrease and are typical for both acute and chronic phases of sialoadenitis [3, 4]. One of the most important tools in xerostomia assessment is the Common Terminology Criteria for Adverse Events (CTCAE v. 4.0) scale, which indicates the intensity of a patient’s subjective complaints.
The apparent radiosensitivity of the salivary glands is still a “radiobiology enigma,” because their cells are highly differentiated, have a slow turnover, and should therefore not be particularly radiosensitive . The histological changes in the gland tissue during irradiation are not accessible in vivo . The ongoing trials conducted to develop radioprotective agents for salivary glands will probably improve gland parenchyma assessment . The diagnostic techniques for salivary gland evaluation are well established; however efforts to introduce new methods are still being made.
Elastography is a term referring to the imaging technique that assesses tissue elasticity in vivo. Qualitative real time ultrasound elastography (ES) is a newly available modality, overlapping elasticity to the grey-scale images. Sonoelastography was first used for predicting malignancy in thyroid, prostate, pancreatic, and breast nodules [8–10]. The usefulness of this novel technique was demonstrated in the evaluation of the status of thyroid parenchyma in de Quervain thyroiditis and liver fibrosis in chronic hepatitis and the assessment of obstructive diseases of the salivary glands [11–14]. The results of elastographic evaluation of parotid and submandibular focal gland masses were then published [15–18].
The question arises of whether elastography can deliver additional information about the irradiated salivary gland and whether this method could therefore help in better qualification of patients for subsequent intensive postradiation treatment. The subject of our research implements the NIH Development Consensus Conference recommendation to develop future research applicable to the management strategies of salivary gland hypofunction and xerostomia . The authors hypothesize that sonoelastography introduced to salivary gland routine repeatable imaging may constitute an objective measure that would correlate with dryness of the mouth. The aim of the study was to assess the features of salivary gland parenchyma by means of high resolution sonoelastography in patients irradiated for HNSCC and to correlate the findings with objective saliva amount and subjective complaints caused by irradiation-induced sialadenitis.
2. Material and Methods
Informed patient consent and Bioethical Committee of Poznan University of Medical Sciences approval (nr 494/12) were obtained. To address the research purpose the investigators designed and implemented the prospective cross-sectioned study of patients with HNSCC aged 48–77, mean 63, median 62 years. Fifty-four healthy volunteers aged 43–78 years, mean 60, constituted the reference group for the elastography values. The patients were matched for age and gender with the control group. To be included in the study sample patients needed to have HNSCC treated with a full course of radiotherapy. All 52 patients were treated for advanced larynx (40), oropharynx (9), and maxilla (3) squamous cancers and suffered from chronic dryness. Forty-three out of the 52 patients were operated on, and adjuvant RT was administered (dose 60–64 Gy); 9/52 were treated by primary radiochemotherapy (RT/CT) (dose 62–64 Gy + 3 Cisplatin courses). Two RT modalities were applied: 33/52 received conventional external beam radiotherapy, and 19/52 received IMRT.
The intensity of xerostomia was assessed immediately after RT and during follow-up after 6, 12, and 24 months. Subjective complaints after RT concerning dry mouth were assessed twice (after RT and during follow-up) in the CTCAE scale. According to this scale, dry mouth sensation is divided into 3 categories: (1) symptomatic (e.g., dry or thick saliva) without significant dietary alteration: unstimulated saliva flow >0.2 mL/min; (2) moderate symptoms: oral intake alterations (e.g., copious water, other lubricants, diet limited to purees and/or soft, moist foods), unstimulated saliva 0.1 to 0.2 mL/min; and (3) inability to adequately aliment orally: tube feeding or total parenteral nutrition (TPN) indicated, unstimulated saliva <0.1 mL/min. Objective salivary volume measurement at follow-up visit encompassed total parotid and submandibular gland secretion. Saliva was collected after 2 minutes of chewing, and the amount was given in grams (g).
Ultrasonography and sonoelastography as a real time “one touch examination” using Aixplorer (SuperSonic Imagine) were performed. Two operators with four years’ experience in the technique conducted the examination. The patients were examined in the supine position. Transverse scans of salivary glands were acquired. A high resolution grey scale US 5–15 MHz linear array transducer was placed both in the middle of the unaffected gland and in heterogenic areas, in the case of alterations in the US picture. Elasticity values in kiloPascal’s (kPa) were collected; the minimum, maximum, and mean values at 3 measurement points in different parts of the glands were taken (Figure 1). The US images were categorized as isoechogenic (50 patients) and with marked fibrosis (2 patients).
The following variables were analyzed: localization of HNSCC primary site, time period from RT. Objective amount of saliva and subjective complaints in the CTCAE scale were taken into consideration. The primary outcome variable for this research was the stiffness of salivary gland parenchyma, measured in kPa. Statistical analysis was performed by means of Spearman’s rank correlation coefficient, Levene’s, Newman-Keuls, Kruskal-Wallis, and Mann-Whitney tests.
In conventional ultrasonography, in all 52 patients except two the submandibular and parotid glands were isoechogenic, without any alterations. In two patients moderate fibrosis was found.
Directly after radiotherapy 16, 29, and 2 patients gained CTCAE scores of 1, 2, and 3, respectively; 5 had no complaints (score “0”). At the follow-up visit 32 and 7 patients had CTCAE scores of 1 and 2; 13 had no dryness (score “0”). The variability of complaints dependent on the time that elapsed from RT was not statistically significant. Taking into consideration the type of RT, the difference between the CTCAE scores in 9 primary (all staged as 2) and 43 adjuvant RT (stages 0, 1, 2, and 3 in 5, 16, 20, and 2, resp.) was significant immediately after treatment (). The difference after one year of follow-up was not noticeable. There was no statistical difference between CTCAE scores in 19 cases after IMRT (1, 4, 14, and 0 patients gained scores of 0, 1, 2, and 3, resp.) and 33 subjects after conventional RT (4, 12, 15, and 3 patients, resp.). The localization of the tumor had no impact on CTCAE scale result. There was no difference between the 37 patients who had undergone RT for the lateral (lymph nodes regions) and the 15 irradiated exclusively for central neck fields (the primary tumor area).
Objective saliva amount collected from the examined group ranged from 0.4 g to 4.75 g (mean 2.15 g, median 1.85 g) and from the control group 2.75 g (). Time from RT had no impact on the mean amount of saliva, which was equal to 1.6 g, 2.3 g, and 2.2 g (median 1.4 g, 1.7 g, 2.1 g) after 6, 12, and 24 months, respectively. There was no statistical difference between 9 primary (mean 1.6 g, median 1.6 g) and 43 adjuvant RT (mean 2.27 g, median 1.9 g) patients. Regarding particular primaries, saliva amount was statistically lower for patients irradiated for larynx (mean 2.13 g, median 1.87 g) than maxilla (mean 4.4 g, median 5.5) () and even lower for the oral cavity (mean 1.5 g, median 1.6 g) (); the patients irradiated for tonsil cancer turned out to be the most prone to dryness. There was no difference between patients who had undergone RT for both neck fields (2.18 g), one side of the neck (2.02 g), central compartment (2.2 g), and exclusively for the primary focus outside the neck (2.0 g).
There was a significant correlation between the obtained subjective CTCAE scale scores and objective saliva amount directly after radiotherapy () and in follow-up (). The mean amount of saliva directly after radiotherapy was 3.5 g, 2.5 g, 1.8 g, and 0.66 g (median 2.93 g, 2.07 g, 1.60 g, and 0.66 g) for the patients with scores 0, 1, 2, and 3; during the follow-up the saliva amount was 3.4 g, 2.0 g, and 0.55 g (median 2.9 g, 1.80 g, 0.53 g) for scores 0, 1, and 2.
There was a statistical difference between elastography values in the irradiated group and in controls in terms of parotid glands (41.7 kPa versus 26.0 kPa; ) and submandibular glands (37.6 kPa versus 22.4 kPa; ). Sonoelastography results of 52 HNSCC patients after RT are shown in Table 1. The values in kPa were taken from the 3 different areas of parenchyma. The values ranged from 4 kPa to 160 kPa. In fact they differed dramatically, although the parenchyma was sonographically isoechogenic. The mean and the median values of the changed parenchyma for the whole group (representing the average highest stiffness value) were 39.6 kPa and 34.9 kPa, respectively; these two variables were taken into consideration in further analysis. The mean and the median value for both parotid glands in the examined group were 41.7 kPa and 40 kPa (Figure 2), while for both submandibular glands they were 37.6 kPa and 34.5 kPa, respectively.
Stiffness of parotid and submandibular glands was more pronounced in patients with primaries localized in the oral cavity (mean 50.7 kPa, median 47.5 kPa) than the larynx (mean 39.73 kPa, median 34.6 kPa) or maxilla (mean 41.2 kPa, median 39.5 kPa) but did not achieve statistical significance. The impact of time since RT for elasticity values is unclear. Directly after RT, in 6, 12, and 24 months the mean elastography values were 34.2 kPa, 36.6 kPa, 46.0 kPa, and 39.7 kPa, respectively, with differences that were not significant.
The mean elasticity values did not significantly differ between the IMRT technique and conventional RT for parotid glands (44.4 kPa versus 40 kPa) and submandibular glands (34.5 kPa versus 39.4 kPa). Noticeable disparities between the primary treatment and adjuvant RT for parotid glands (46.2 kPa versus 33.6 kPa) and submandibular glands (39.5 kPa versus 30.3 kPa) were noted; nevertheless the scores did not achieve statistical significance. There was no correlation between the CTCAE score (subjective complaints) and elastography values. Patients with scores 0, 1, 2, and 3 received mean elastographic values 37.3 kPa, 45.4 kPa, 37.1 kPa, and 35.5 kPa, respectively, from all major glands. Finally, there was no correlation between the amount of saliva and stiffness.
The clinical problem, radiation-induced sialadenitis, causing the disabled function of salivary glands and dryness syndrome, is the most common cause of the decrease in life quality after external radiotherapy in the majority of HNSCC patients , and thus the significance of appropriate diagnostics is crucial. Reliable radiological methods to evaluate xerostomia are also necessary; ultrasonography, computed tomography, and magnetic resonance are the methods of choice in imaging the status of parenchyma . Sialography, sialoendoscopy, and MR sialography are useful in Stensen’s duct visualization. The measurements of saliva amount and examination of sediment indicate the secretory function of the glands . Nonetheless, binding the parenchyma features with secretory abilities is difficult. The objective methods, scintiscan, single-photon emission computed tomography (SPECT) , and dynamic MR sialography [23, 24], are more difficult to perform, especially in acute phase of radiation sialadenitis. One currently available technique that can visualize the major salivary glands and does not require any invasive procedures, like cannulation and contrast media infusion, is MR sialography . In addition, noninvasiveness of this method is suitable for chronologic observation of gland hypofunction. However, for the prediction of radiation-induced xerostomia these techniques have disadvantages, including unavailability, high cost, and long duration of the test . Katz et al. have shown that ultrasonography with high resolution transducers, routinely used in neck node monitoring in head and neck follow-up, is a highly sensitive method of evaluating the major salivary glands . Nevertheless, in chronic inflammation the glandular modifications seen on US are often less prominent than in acute diseases; unfortunately in this technique no functional data are available. In the early, acute phase of radiation sialadenitis the gland is anechoic . Furthermore, if later glandular atrophy occurs, marked reduction of volume with higher echogenicity appears, but still the sonographic picture does not represent the secreting ability of the glands, and thus its correlation with subjective patient complaints is not possible. According to Zenk , in the first, acute phase the gland is swollen and anechoic. If later glandular atrophy occurs, a marked reduction of volume with higher echogenicity may be found. Our results confirmed that in fact the sonographic picture does not reflect the secreting ability of the glands, and discrete sonographic changes rarely correspond with the subjective patient complaints. In all our patients except two, normoechogenic tissue was visualized. To conclude, the imaging possibilities of ultrasonography in radiation-induced sialadenitis are not in our opinion sufficient, and the usefulness of this method in predicting the xerostomia is very limited.
Elastography, a new development in ultrasound imaging, may deliver additional information about parenchyma status . Elastography is a term referring to an imaging technique that assesses tissue elasticity in vivo. It provides a new insight into the evaluation of gland parenchyma that may be translated into clinical applications in the field of irradiation consequences in HNSCC treatment, but there are no literature data concerning this issue. The application of sonoelastography in our clinical practice tries to answer the need for the development of an accurate, quantifiable, and reproducible assessment of salivary glands in HNSCC patients. Tissue elasticity in salivary glands probably reflects a combination of factors, including the quantity and quality of salivary glandular tissue, and other components, such as fatty tissue or fibrosis, which may be increased in irradiated glands.
The presented study confirmed a correlation between tissue elasticity and postradiative salivary gland hypofunction. The most important finding in this study was the statistical difference in parenchymal elasticity between the RT and the reference group. The mean values from healthy volunteers (22 kPa for submandibular and 26 kPa for parotid) were statistically lower than in irradiated glands (38 kPa and 42 kPa). Unfortunately, we were unable to find correlations between the specified classes inside the RT group. Elasticity values did not correlate with CTCAE score, which can probably be explained by the small number of patients in particular the subgroups. As the time in radiotherapy plays a great role in radiobiological processes and patients’ wellbeing, we pointed this out separately. The biggest difference in the mean elastography values was directly after RT (34.2 kPa) and after 12 months (46 kPa); in the further follow-up (24 months) the value became lower (40 kPa), but this trend needs further examination. The second issue is the localization of the irradiated primary. The amount of saliva indicated that patients irradiated for tonsil cancer were most liable to experience dryness. However, this fact was not reflected in CTCAE scores. Stiffness in the parotid and submandibular glands was higher in patients with primaries localized in the oral cavity (50.7 kPa and 38.55 kPa) as compared to the larynx (39.73 kPa and 34.6 kPa) and maxilla (41.2 kPa, 30 kPa), but the differences did not reach statistical significance. Similarly, the mean elastography value for the glands on the irradiated side of the neck was 40.64 kPa, and for the nonirradiated side it was 35.19 kPa (parotid 43.2 kPa versus 36.8 kPa, submandibular 38.4 kPa versus 33.5 kPa, resp.), and the trend is noticeable.
Unfortunately, our study and the method itself have some limitations. Xerostomia refers to salivary secretory quantity. It does not follow the elasticity of the glandular parenchyma. Based on these facts we can assume that the elasticity of salivary tissue depends on many antagonistic factors (e.g., fatty replacement and fibrosis) and does not correspond directly to salivary function. Increased fibrosis, which is difficult to predict in irradiated glands, may or may not influence the secretory properties of the gland. There is also a wide overlap in the values of elasticity obtained from healthy and irradiated parenchyma; therefore, it is difficult to indicate cut-off values differentiating between normal and affected glands. The number and size of regions of interests need also to be standardized. The limitation of our group is quite heterogeneous in the character of the studied patients and in the prospective assessment of all studied subjects in every point of the follow-up that were not always available. Larger prospective studies are needed to fully assess the changes of elasticity during the follow-up after radiation therapy.
To conclude, this study suggests that sonoelastography adds a new parameter to the grayscale ultrasound images in “one touch examination” and may be a useful tool for salivary gland evaluation during the radiotherapy course and the follow-up. In spite of the fact that this method has not provided detailed information about the internal architecture and morphological changes of the gland, increased elasticity corresponds to the radiobiological effect and may potentially identify individuals more liable for sialadenitis during RT course.
Conflict of Interests
There is no potential conflict of interests regarding this paper.
The study was approved by the Bioethical Committee of Poznan University of Medical Sciences (nr 494/12).
- A. Wada, N. Uchida, M. Yokokawa, T. Yoshizako, and H. Kitagaki, “Radiation-induced xerostomia: objective evaluation of salivary gland injury using MR sialography,” American Journal of Neuroradiology, vol. 30, no. 1, pp. 53–58, 2009.
- S. B. Jensen, A. M. L. Pedersen, A. Vissink et al., “A systematic review of salivary gland hypofunction and xerostomia induced by cancer therapies: prevalence, severity and impact on quality of life,” Supportive Care in Cancer, vol. 18, no. 8, pp. 1039–1060, 2010.
- S. B. Jensen, A. M. L. Pedersen, A. Vissink et al., “A systematic review of salivary gland hypofunction and xerostomia induced by cancer therapies: management strategies and economic impact,” Supportive Care in Cancer, vol. 18, no. 8, pp. 1061–1079, 2010.
- M. S. Chambers, A. S. Garden, M. S. Kies, and J. W. Martin, “Radiation-induced xerostomia in patients with head and neck cancer: pathogenesis, impact on quality of life, and management,” Head and Neck, vol. 26, no. 9, pp. 796–807, 2004.
- R. M. Nagler, “The enigmatic mechanism of irradiation-induced damage to the major salivary glands,” Oral Diseases, vol. 8, no. 3, pp. 141–146, 2002.
- M. J. van Holten, J. M. Roesink, C. H. J. Terhaard, and P. M. Braam, “New insights in the vascular supply of the human parotid gland—consequences for parotid gland-sparing irradiation,” Head and Neck, vol. 32, no. 7, pp. 837–843, 2010.
- A. Teymoortash, F. Müller, J. Juricko et al., “Botulinum toxin prevents radiotherapy-induced salivary gland damage,” Oral Oncology, vol. 45, no. 8, pp. 737–739, 2009.
- K. S. S. Bhatia, D. P. Rasalkar, Y. P. Lee et al., “Cystic change in thyroid nodules: a confounding factor for real-time qualitative thyroid ultrasound elastography,” Clinical Radiology, vol. 66, no. 9, pp. 799–807, 2011.
- N. Kawada, S. Tanaka, H. Uehara et al., “Feasibility of second-generation transabdominal ultrasound-elastography to evaluate solid pancreatic tumors: preliminary report of 36 cases,” Pancreas, vol. 41, no. 6, pp. 978–980, 2012.
- A. Evans, P. Whelehan, K. Thomson et al., “Differentiating benign from malignant solid breast masses: value of shear wave elastography according to lesion stiffness combined with greyscale ultrasound according to BI-RADS classification,” British Journal of Cancer, vol. 107, no. 2, pp. 224–229, 2012.
- M. Ruchala, E. Szczepanek, and J. Sowinski, “Sonoelastography in de Quervain thyroiditis,” Journal of Clinical Endocrinology and Metabolism, vol. 96, no. 2, pp. 289–290, 2011.
- X. P. Ye, H. T. Ran, J. Cheng et al., “Liver and spleen stiffness measured by acoustic radiation force impulse elastography for noninvasive assessment of liver fibrosis and esophageal varices in patients with chronic hepatitis B,” Journal of Ultrasound in Medicine, vol. 31, no. 8, pp. 1245–1253, 2012.
- P. Zengel, F. Schrötzlmair, F. Schwarz et al., “Elastography: a new diagnostic tool for evaluation of obstructive diseases of the salivary glands; primary results,” Clinical Hemorheology and Microcirculation, vol. 50, no. 1-2, pp. 91–99, 2012.
- G. Ferraioli, C. Tinelli, B. Dal Bello, M. Zicchetti, G. Filice, and C. Filice, “Accuracy of real-time shear wave elastography for assessing liver fibrosis in chronic hepatitis C: a pilot study,” Hepatology, vol. 56, no. 6, pp. 2125–2133, 2012.
- K. S. S. Bhatia, C. C. M. Cho, C. S. L. Tong, Y. Y. P. Lee, E. H. Y. Yuen, and A. T. Ahuja, “Shear wave elastography of focal salivary gland lesions: preliminary experience in a routine head and neck US clinic,” European Radiology, vol. 22, no. 5, pp. 957–965, 2012.
- D. Dumitriu, S. Dudea, C. Botar-Jid, M. Bǎciuţ, and G. Bǎciuť, “Real-time sonoelastography of major salivary gland tumors,” American Journal of Roentgenology, vol. 197, no. 5, pp. W924–W930, 2011.
- N. Klintworth, K. Mantsopoulos, J. Zenk, G. Psychogios, H. Iro, and A. Bozzato, “Sonoelastography of parotid gland tumours: initial experience and identification of characteristic patterns,” European Radiology, vol. 22, no. 5, pp. 947–956, 2012.
- O. Westerland and D. Howlett, “Sonoelastography techniques in the evaluation and diagnosis of parotid neoplasms,” European Radiology, vol. 22, no. 5, pp. 966–969, 2012.
- P. Dirix, S. Nuyts, V. V. Poorten, P. Delaere, and W. van den Bogaert, “The influence of xerostomia after radiotherapy on quality of life: results of a questionnaire in head and neck cancer,” Supportive Care in Cancer, vol. 16, no. 2, pp. 171–179, 2008.
- T. Kan, K. Kodani, K. Michimoto, S. Fujii, and T. Ogawa, “Radiation-induced damage to microstructure of parotid gland: evaluation using high-resolution magnetic resonance imaging,” International Journal of Radiation Oncology, Biology, Physics, vol. 77, no. 4, pp. 1030–1038, 2010.
- M. S. Chambers, K. L. Tomsett, I. I. Artopoulou et al., “Salivary flow rates measured during radiation therapy in head and neck cancer patients: a pilot study assessing salivary sediment formation,” Journal of Prosthetic Dentistry, vol. 100, no. 2, pp. 142–146, 2008.
- A. Nömayr, W. Römer, D. Strobel, W. Bautz, and T. Kuwert, “Anatomical accuracy of hybrid SPECT/spiral CT in the lower spine,” Nuclear Medicine Communications, vol. 27, no. 6, pp. 521–528, 2006.
- M. Habu, T. Tanaka, T. Tomoyose et al., “Significance of dynamic magnetic resonance sialography in prognostic evaluation of saline solution irrigation of the parotid gland for the treatment of xerostomia,” Journal of Oral and Maxillofacial Surgery, vol. 68, no. 4, pp. 768–776, 2010.
- E. Astreinidou, J. M. Roesink, C. P. J. Raaijmakers et al., “3D MR sialography as a tool to investigate radiation-induced xerostomia: feasibility study,” International Journal of Radiation Oncology Biology Physics, vol. 68, no. 5, pp. 1310–1319, 2007.
- M. Kalinowski, J. T. Heverhagen, E. Rehberg, K. J. Klose, and H. Wagner, “Comparative study of MR sialography and digital subtraction sialography for benign salivary gland disorders,” American Journal of Neuroradiology, vol. 23, no. 9, pp. 1485–1492, 2002.
- P. Katz, D. M. Hartl, and A. Guerre, “Clinical ultrasound of the salivary glands,” Otolaryngologic Clinics of North America, vol. 42, no. 6, pp. 973–1000, 2009.
- J. L. Zenk, “New therapy for the prevention and prophylactic treatment of acute radiation syndrome,” Expert Opinion on Investigational Drugs, vol. 16, no. 6, pp. 767–770, 2007.
- M. Ying, K. S. Bhatia, Y. P. Lee et al., “Review of ultrasonography of malignant neck nodes: greyscale, Doppler, contrast enhancement and elastography,” Cancer Imaging, vol. 13, no. 4, pp. 658–669, 2013.
Copyright © 2014 Jarosław Kałużny et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.