Real-imaging cDNA-AFLP transcript profiling of pancreatic cancer patients: Egr-1 as a potential key regulator of muscle cachexia
- Alexander Skorokhod†1, 4Email author,
- Jeannine Bachmann†2,
- Nathalia A Giese3,
- Marc E Martignoni2 and
- Holger Krakowski-Roosen1
© Skorokhod et al.; licensee BioMed Central Ltd. 2012
Received: 19 December 2011
Accepted: 24 May 2012
Published: 21 June 2012
Cancer cachexia is a progressive wasting syndrome and the most prevalent characteristic of cancer in patients with advanced pancreatic adenocarcinoma. We hypothesize that genes expressed in wasted skeletal muscle of pancreatic cancer patients may determine the initiation and severity of cachexia syndrome.
We studied gene expression in skeletal muscle biopsies from pancreatic cancer patients with and without cachexia utilizing Real-Imaging cDNA-AFLP-based transcript profiling for genome-wide expression analysis.
Our approach yielded 183 cachexia-associated genes. Ontology analysis revealed characteristic changes for a number of genes involved in muscle contraction, actin cytoskeleton rearrangement, protein degradation, tissue hypoxia, immediate early response and acute-phase response.
We demonstrate that Real-Imaging cDNA-AFLP analysis is a robust method for high-throughput gene expression studies of cancer cachexia syndrome in patients with pancreatic cancer. According to quantitative RT-PCR validation, the expression levels of genes encoding the acute-phase proteins α-antitrypsin and fibrinogen α and the immediate early response genes Egr-1 and IER-5 were significantly elevated in the skeletal muscle of wasted patients. By immunohistochemical and Western immunoblotting analysis it was shown, that Egr-1 expression is significantly increased in patients with cachexia and cancer. This provides new evidence that chronic activation of systemic inflammatory response might be a common and unifying factor of muscle cachexia.
Pancreatic cancer remains one of the most deadly tumor types, with the highest death-to-incidence ratio of all cancers. With a 5-year survival rate of less than 5% and a death-to-incidence ratio of 0.99, pancreatic ductal adenocarcinoma is currently one of the most aggressive gastrointestinal carcinomas . The development of cachexia occurs in most patients with resectable and with unresectable pancreatic cancer, whose resting energy expenditure is generally increased. The term “cachexia” is used in the literature with different definitions and descriptions of clinical symptoms of a patient. It has been shown that a loss of weight exceeding 10% of the stable pre illness weight is correlated to a worse prognosis .
In the case of muscle wasting associated with cancer, the balance between protein synthesis and degradation is shifted markedly toward degradation, and there could be selective targeting of skeletal muscle gene products by different cachectic factors such as cytokines and tumor-derived compounds [3, 4]. Weight loss of pancreatic cancer patients has been also related to the induction of acute-phase response, resulting in the loss of both fat and lean tissue. It has been demonstrated that a complex cascade of inflammatory responses may mediate alteration in metabolic processes, contributing to the development of cancer cachexia .
Unfortunately, patients with pancreatic cancer generally present in late stage of the disease, and frequently develop cachexia syndrome. This emphasizes the importance of early diagnosis, especially with the hope that intervention may be more effective in the early stages of disease. This also underscores the need for biomarkers, which may also contribute to the understanding of cancer cachexia syndrome . Large-scale analysis of gene expression has been widely proposed as a powerful method for malignancy diagnosis, predicting invasion and metastasis through the identification of biomarkers. To our knowledge, however, transcriptome-wide analysis of muscle cachexia in patients has not yet been performed.
Here we present our study that focuses on the identification of common differentially expressed genes involved in muscle wasting progression. We hypothesized that gene expression patterns in the skeletal muscle of cancer patients will change under muscle wasting conditions, and that these changes may determine the initiation and severity of cachexia syndrome. We suggest that characterizing the skeletal muscle in response to cancer cachexia will provide new insights into gene regulation.
Patients with pancreatic cancer were selected for this study, as these patients usually experience severe weight loss (weight loss exceeding 10% of the stable body weight). All included patients gave written informed consent for tissue and data collection under the local ethics committee’s permission of the University of Heidelberg.
On admission to hospital, each patient was asked for his height and current weight. When a patient gave a history of weight loss, the corresponding time period was also documented. A patient was classified as cachectic when weight loss exceeded 10% of the stable pre-illness weight. It has been shown, that in patients with a history of weight loss exceeding 10% of the stable weight, the survival is significantly reduced [2, 7].
Histological diagnosis: Two independent pathologists from the University of Heidelberg’s Department of Pathology confirmed the diagnosis according to the UICC guidelines .
Biopsy specimens were obtained from the rectus abdominis muscle of each patient during the initial phase of the operation. After skin incision and dissection through the subcutaneous fat, the anterior sheet of the rectus abdominis muscle was opened with scissors and a muscle biopsy specimen weighing about 50 mg was obtained and immediately frozen in liquid nitrogen.
Real-Imaging cDNA-AFLP transcript profiling
mRNA was extracted from 20 mg rectus abdominis muscle using MagNA pure LC isolation system according to manufacturer’s instructions (Roche). After extraction, samples were immediately stored at −80°C until usage. mRNA quality was assessed by reverse transcription PCR (RT-PCR) with two housekeeping genes encoding ribosomal protein L13a (RPL13A) and β-actin. Negative controls for genomic contamination and reverse transcription inhibition were included.
The cDNA-AFLP procedure was performed as described previously with some modifications . cDNA was synthesized from 100 ng mRNA according to the manufacturer’s instructions (cDNA Synthesis Kit, Roche). After restriction enzyme digestion, the cDNA fragments were ligated to oligomers introducing the priming sites for TaqI and MseI. The cDNA pools were amplified with primers TaqI-Ø (5’-CTCGTAGACTGCGTACACGA-3’) and MseI-Ø (5’-GACGATGAGTCCTGAGTAA-3’).
cDNA-AFLP amplification was performed with 16 primer pairs MseI+2 and TaqI+2IRD700 (+2 states for additional two-nucleotide combination at 3’-end of primer) generating approximately 30,000 amplified fragments possessing 256 probe sets. IRD-labeled PCR products were resolved in 40 cm long 6% SequaGel XR gel (Biozym) using automated DNA sequencer DNAIR 4200 (Li-COR). After 1 h electrophoresis the gel was transferred onto the green grid paper and vacuum dried on 1.2-mm-thick Whatman paper. The dried gel was scanned at 700/800 nm (intensity 4.0, resolution 84 μm), obtaining binary data for differentially expressed tags (Odyssey Imaging System, Li-COR).
Isolation and analysis of cDNA-AFLP tag fragments
The cDNA fragments of interest were excised from dried gel, followed by elution in 100 μl TE for 10 min at 95°C. 5 μl were re-amplified with TaqI-Ø and MseI-Ø primers. Specific fragments were verified by visualizing 5 μl PCR product on a 1.5% agarose gel. The amplified fragments were directly sequenced with the primers TaqI-ØIRD700 and MseI-ØIRD800 using the SequiTherm EXCEL II DNA sequencing kit (Epicentre) on the automated infrared sequencer DNAIR 4200 (Li-COR).
Sequencing data were aligned against BlastX and BlastN in the publicly available databases and clustered according to gene ontology using the nodes of the DNASpace software package (DNASIS MAX 2.6, Miraibio).
Single-stranded cDNA was synthesized from 100 ng of mRNA according to the manufacturer’s instructions (RevertAid, Fermentas). Primers complementary to target genes were designed using Primer Express Software (Life Technologies) to amplify fragments in the range from 75 bp to 200 bp, close to the 3’-end of transcripts. A quantitative RT-PCR (qPCR) was performed with SYBR Green PCR Master Mix using GeneAmp 7300 (Life Technologies). qPCR reactions were performed two times in doublets, and the threshold cycle numbers were averaged. Levels of expression were normalized against RPL13A and β-actin.
For immunohistochemistry, 5-μm sections of FFPE tissue were incubated in serum-free protein blocking buffer for 30 min, followed by overnight incubation at 4°C with 1:200 dilution of Egr-1 antibodies (588, Santa Cruz). A negative control for non-specific staining was performed with rabbit IgG. DAKO anti-rabbit polymer-HRPO was applied for signal visualization according to the manufacturer’s instructions (DAKO). Finally, sections were analyzed with an Olympus BX50 microscope using 40x-lens and software program “Cell^F 2.5” (Olympus).
Protein extraction and immunoblotting
30 μg total protein lysate from the biopsy homogenate were electrophoresed on gradient 4-15% SDS-PAA gels (Bio-RAD). After undergoing electrophoresis, proteins were transferred semi-dried onto PVDF membrane using the Semi-Dry Blotter (Bio-RAD). After 1 h blocking, the membrane was exposed overnight at 4°C to primary anti-human rabbit Egr-1 antibody (588, Santa Cruz) and anti-human mouse β-tubulin antibody (3 F3-G2, Santa Cruz) at 1:200.
After washing the membrane was exposed to appropriate IRDye-700/800 labeled secondary antibodies: Alexa Fluor 680 goat anti-mouse for β-tubulin (Invitrogen) and IRDye-800 goat anti-rabbit for Egr-1 (Rockland) at 1:10,000. The amounts of Egr-1 and β-tubulin were determined and quantified by densitometric scanning of PVDF membranes at 700/800 nm (intensity 5.0, resolution 84 μm) using the ODYSSEY Infrared Imaging System (Li-COR). Individual band intensities were normalized against β-tubulin.
All data are presented as means ± standard deviation. The significance of differences between groups (cachexia vs. no cachexia) was tested with Student’s t-test and Mann–Whitney U test. The statistical package SPSS for Windows (SPSS) was used for all statistical analyses. Differences were considered significant at the 0.05 level of confidence.
Clinical characteristics of patients with pancreatic cancer with (N = 10) and without cachexia (N = 13): gender and performance status
No cachexia (N = 13) N [%]
Cachexia (N = 10) N [%]
Characteristics of patients with pancreatic cancer with (N = 10) and without cachexia (N = 13)
No cachexia (N = 13)
Cachexia (N = 10)
Real-Imaging cDNA-AFLP screening
Analysis of differentially expressed genes
183 cDNA sequences were identified after a redundancy check of about 400 isolated tag fragments. They were clustered into groups according to their annotation fields. 127 cDNA fragments were referred to the gene transcripts that predominantly activated in cachectic patients, as well as 56 individual genes whose transcripts accumulated predominantly in the non-cachectic control group (Additional file 2: Table S1).
According to group classification, the identified genes provide molecular counterparts of muscle cachexia responsible for actin cytoskeleton rearrangement, actin filament disruption, extracellular matrix dysregulation, protein degradation, energy expenditure, tissue hypoxia, immediate early response, and acute-phase response. We identified a number of differentially expressed genes, which have already been characterized in relation to cancer cachexia or muscle atrophy. However, our results also indicate that there are other differentially expressed genes, which play unknown roles in muscle cachexia.
In wasted patients we observed the upregulation in gene expression of the contractive muscle proteins (Additional file 2: Table S1), correlating directly with other findings . Individual differences were found in expression of the genes belonging to different proteolytic pathways (Additional file 2: Table S1, “Protein Degradation”). A number of genes encoding both transcription and translation factors as well as different ribosomal proteins were induced in response to muscle cachexia. In excellent agreement with previous findings, we observed the upregulation of class O type forkhead transcription factor 1 (FOXO1), which has been implicated as a key regulator of gene expression during skeletal muscle atrophy . We have also identified alteration of the genes that were previously described in association with hypoxia (Additional file 2: Table S1, “Hypoxia”). Results of our cDNA-AFLP screening revealed that the genes encoding acute-phase response (APR) proteins fibrinogen α (FGA), fibrinogen γ (FGG), α1-antitrypsin, and α2-macroglobulin are induced in the skeletal muscle of pancreatic cancer patients with cachexia.
Validation of cachexia-associated genes with qPCR
mRNA expression level of cDNA-AFLP candidate genes – cachexia vs. no cachexia
forkhead transcription factor FOXO1A
immediate early response 5 (IER-5)
karyopherin alpha (importin 5)
breast cancer-associated protein
mt translational release factor 1 (MTRF 1)
myocyte enhancer factor 2 (MEF2)
cardiomyopathy associated protein 5
myosin regulatory light chain (MYLC)
myosin heavy chain 7 (MYH7)
cyclin-dependent kinase inhibitor 1A (p21)
NK transcript (IL-32)
bHLH transcription factor DEC1
bHLH transcription factor DEC2
bHLH transcription factor HES1
early growth response 1 (Egr-1)
Egr-1 protein level is highly elevated in wasted skeletal muscle
Finally, immunohistochemistry and Western immunoblotting analysis confirmed the strong correlation of Egr-1 within the skeletal muscle of pancreatic cancer patients in response to cachexia.
Cachexia is present in around 80% of patients with cancer in advanced stages, inducing marked depletion of lean and fat mass, and in more than 20% of patients cachexia is responsible for death . Whereas starvation and chronic malnutrition are associated with adaptive decreases in metabolic rate, cancer cachexia is characterized by altered metabolic processes and activation of catabolic pathways .
We hypothesized that alteration in gene expression within the skeletal muscle in response to advanced pancreatic cancer plays a crucial role in the development of the muscle cachexia syndrome. To test this hypothesis, we carried out the transcriptome-wide characterization of skeletal muscle biopsies utilizing the Real-Imaging cDNA-AFLP approach to identify those genes that are commonly induced in the skeletal muscle of pancreatic cancer patients either with or without muscle cachexia.
In contrast to other studies of skeletal muscle cancer cachexia based on cDNA microarray technique[14, 15], our Real-Imaging cDNA-AFLP approach gives the ability to identify and annotate expressed genes without prior sequence knowledge making possible to distinguish among different transcripts from genes belonging to the same gene family or splice variants of the same gene. Importantly, our results indicate that molecular signatures of gene expression appear to be predictive of muscle cachexia, before we can measure this in the preoperative collected blood samples as there is no significant difference in protein levels or CrP-levels to detect. We suggest that in characterizing the skeletal muscle transcriptome, one might expect that different signal systems may work in combination to shift the disease stage of pancreatic cancer patients towards cachectic conditions.
We observed that transcript levels of contractile proteins in skeletal muscle of patients with cachexia change in the same direction. We observed more than 6-fold upregulation of cyclin-dependent kinase inhibitor p21 (Table 3). The mRNA levels of increased p21 mRNA expression in cachectic patients may reflect increased activity of myostatin, a negative regulator of muscle growth and development which is thought to induce cachexia . The skeletal muscle mass appears to be significantly involved in the process of cancer cachexia in patients with pancreatic cancer, being the potential target of various tumor products which can activate protein degradation . All these changes may indicate the abnormal metabolic processes shifting towards hypermetabolism and negative energy expenditure in patients with cachexia.
We have also observed the induction of genes encoding basic helix-loop-helix (bHLH) transcription factors HES1, DEC1, and DEC2, which are primarily associated with hypoxia (Table 3). As it was reported previously, hypoxia enhances Notch signaling and activates Notch-responsive promoters and increases expression of bHLH transcription factors . HES1 and DEC1 (Table 3) are immediate downstream targets of Notch signaling inhibiting the transcriptional expression and/or function of lineage-specifying genes involved in myogenesis .
Our study demonstrates for the first time that the genes encoding acute-phase response (APR) proteins fibrinogen α, α1-antitrypsin and α2-macroglobulin are markedly upregulated not only in the liver but also in the skeletal muscle of wasted patients with pancreatic cancer (Table 3), which are associated with altered energy expenditure and increased protein metabolism of pancreatic cancer patients . It has been shown that approximately 40 percent of pancreatic cancer patients show an acute-phase response at the time of diagnosis, and this amount increases to around 80 percent at the time of death .
While acute-phase response indicates systemic inflammation, the pro-inflammatory cytokines appear to mediate this response. Furthermore, it has been reported that non-immune tissues such as cardiac and skeletal muscle appears to express pro-inflammatory cytokines . In excellent agreement with these facts, we observed 6-fold upregulation of interleukin IL-32 (Table 3), play an important role in inflammatory diseases as direct inducer of TNF-α .
Notably, in contrast to the previous study on human skeletal muscle cancer cachexia , we observed a more than 5-fold upregulation of FOXO1, a key regulator of gene expression during skeletal muscle atrophy  (p<0.05, Table 3). We suppose that it may be explained by more stringent selection of cachectic patients in our study with weigh loss criteria from 10 to 19% within 6 months (Table 1).
Importantly, together with an ongoing inflammatory response we observed the significant induction of early response genes Egr-1 and IER-5 in wasted pancreatic cancer patients (p<0.05, Table 3). Conceptually, Egr-1, the master switch regulating inflammatory parameters, can be viewed as a main initiator of inflammatory response, as it is directly involved in “early signal events”[24, 25]. Egr-1 is also recognized as playing a critical role in the regulation of over 40 target genes, including tumor necrosis factor α (TNF-α), interleukin (IL)-1β, and IL-6, all of which can induce various features of cachexia syndrome in the skeletal muscle in animal models . Together with elevated acute phase response gene expression, the Egr-1 may serve as a biomarker of muscle cachexia syndrome.
This report provides the first evidence that Egr-1 is aberrantly upregulated under conditions of muscle wasting, underlining the general importance of pro-inflammatory response in cancer cachexia (Figures 4 and 5). Together with altered IER-5 and ongoing acute-phase response, it adds to the body of evidence that systemic inflammatory response contributes mainly to cancer cachexia progression mediating altered metabolic processes, considering this condition as a chronic systemic inflammatory disease.
In summary, taken together with previous works in the field, the current study supports a model of muscle wasting in which the major influence of cancer cachexia is based on the chronic inflammatory effects resulting in onset of acute-phase response. Furthermore, this report provides new insights into the inflammatory response involved in initiation and progression of cachexia by the activation of pro-inflammatory transcription factor Egr-1 and acute-phase response. How they are transduced to initiate the muscle cachexia remains to be investigated.
cDNA amplified fragment-length polymorphism
early growth response
immediate early response 5
We thank Brunhilde Bentzinger and Inna Babushkina for their technical assistance in setting up the immunohistochemical experiments.
This work was supported by German Cancer Research Center.
Grant funding for research
A. Skorokhod and H. Krakowski-Roosen declare:
These two authors had financial support from the German Cancer Research Centre for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous 3 years; no other relationships or activities that could appear to have influenced the submitted work.
- Jemal A, Siegel R, Ward E, Hao Y, Xu J, Thun MJ: Cancer statistics, 2009. CA Cancer J Clin. 2009, 59: 225-249. 10.3322/caac.20006.View ArticlePubMed
- Bachmann J, Heiligensetzer M, Krakowski-Roosen H, Buchler MW, Friess H, Martignoni ME: Cachexia worsens prognosis in patients with resectable pancreatic cancer. J GastrointestSurg. 2008, 12: 1193-1201.View Article
- Sacheck JM, Hyatt JP, Raffaello A, et al: Rapid disuse and denervation atrophy involve transcriptional changes similar to those of muscle wasting during systemic diseases. FASEB J. 2007, 21: 140-155.View ArticlePubMed
- Argiles JM, Busquets S, Felipe A, Lopez-Soriano FJ: Molecular mechanisms involved in muscle wasting in cancer and ageing: cachexia versus sarcopenia. Int J Biochem Cell Biol. 2005, 37: 1084-1104. 10.1016/j.biocel.2004.10.003.View ArticlePubMed
- Deans C, Wigmore SJ: Systemic inflammation, cachexia and prognosis in patients with cancer. CurrOpinClinNutrMetab Care. 2005, 8: 265-269.
- Fortunati N, Manti R, Birocco N, et al: Pro-inflammatory cytokines and oxidative stress/antioxidant parameters characterize the bio-humoral profile of early cachexia in lung cancer patients. Oncol Rep. 2007, 18: 1521-1527.PubMed
- Bachmann J, Ketterer K, Marsch C, et al: Pancreatic cancer related cachexia: influence on metabolism and correlation to weight loss and pulmonary function. BMC Cancer. 2009, 9: 255-10.1186/1471-2407-9-255.PubMed CentralView ArticlePubMed
- Wittekind C: 2010 TNM system: on the 7th edition of TNM classification of malignant tumors. Pathologe. 2010, 31: 331-332.View ArticlePubMed
- Dellagi A, Birch PR, Heilbronn J, Lyon GD, Toth IK: cDNA-AFLP analysis of differential gene expression in the prokaryotic plant pathogen Erwiniacarotovora. Microbiology. 2000, 146 (Pt 1): 165-171.View ArticlePubMed
- Ferrando AA, Wolfe RR: Restoration of hormonal action and muscle protein. Crit Care Med. 2007, 35: S630-S634. 10.1097/01.CCM.0000278529.44899.57.View ArticlePubMed
- Sandri M, Sandri C, Gilbert A, et al: Foxo transcription factors induce the atrophy-related ubiquitin ligase atrogin-1 and cause skeletal muscle atrophy. Cell. 2004, 117: 399-412. 10.1016/S0092-8674(04)00400-3.PubMed CentralView ArticlePubMed
- Inui A, Meguid MM: Cachexia and obesity: two sides of one coin?. CurrOpinClinNutrMetab Care. 2003, 6: 395-399.
- Baracos VE: Hypercatabolism and hypermetabolism in wasting states. CurrOpinClinNutrMetab Care. 2002, 5: 237-239.
- Stephens NA, Gallagher IJ, Rooyackers O, et al: Using transcriptomics to identify and validate novel biomarkers of human skeletal muscle cancer cachexia. Genome Med. 2010, 2: 1-10.1186/gm122.PubMed CentralView ArticlePubMed
- Stevenson EJ, Koncarevic A, Giresi PG, Jackman RW, Kandarian SC: Transcriptional profile of a myotube starvation model of atrophy. J ApplPhysiol. 2005, 98: 1396-1406.
- Jespersen J, Kjaer M, Schjerling P: The possible role of myostatin in skeletal muscle atrophy and cachexia. Scand J Med Sci Sports. 2006, 16: 74-82. 10.1111/j.1600-0838.2005.00498.x.View ArticlePubMed
- Martignoni ME, Kunze P, Hildebrandt W, et al: Role of mononuclear cells and inflammatory cytokines in pancreatic cancer-related cachexia. Clin Cancer Res. 2005, 11: 5802-5808. 10.1158/1078-0432.CCR-05-0185.View ArticlePubMed
- Gustafsson MV, Zheng X, Pereira T, et al: Hypoxia requires notch signaling to maintain the undifferentiated cell state. Dev Cell. 2005, 9: 617-628. 10.1016/j.devcel.2005.09.010.View ArticlePubMed
- Carey KA, Farnfield MM, Tarquinio SD, Cameron-Smith D: Impaired expression of Notch signaling genes in aged human skeletal muscle. J Gerontol A BiolSci Med Sci. 2007, 62: 9-17. 10.1093/gerona/62.1.9.View Article
- Moses AG, Maingay J, Sangster K, Fearon KC, Ross JA: Pro-inflammatory cytokine release by peripheral blood mononuclear cells from patients with advanced pancreatic cancer: Relationship to acute phase response and survival. Oncol Rep. 2009, 21: 1091-1095.PubMed
- Fearon KC, Voss AC, Hustead DS: Definition of cancer cachexia: effect of weight loss, reduced food intake, and systemic inflammation on functional status and prognosis. Am J Clin Nutr. 2006, 83: 1345-1350.PubMed
- Carballo-Jane E, Pandit S, Santoro JC, et al: Skeletal muscle: a dual system to measure glucocorticoid-dependent transactivation and transrepression of gene regulation. J Steroid BiochemMolBiol. 2004, 88: 191-201. 10.1016/j.jsbmb.2003.11.012.View Article
- Dinarello CA, Kim SH: IL-32, a novel cytokine with a possible role in disease. Ann Rheum Dis. 2006, 65 (Suppl 3): iii61-iii64.PubMed CentralPubMed
- Schmidt J, Stoffels B, Moore BA, et al: Proinflammatory role of leukocyte-derived Egr-1 in the development of murine postoperative ileus. Gastroenterology. 2008, 135: 926-936. 10.1053/j.gastro.2008.05.079. 936PubMed CentralView ArticlePubMed
- Yan SF, Fujita T, Lu J, et al: Egr-1, a master switch coordinating upregulation of divergent gene families underlying ischemic stress. Nat Med. 2000, 6: 1355-1361. 10.1038/82168.View ArticlePubMed
- McNamara MJ, Alexander HR, Norton JA: Cytokines and their role in the pathophysiology of cancer cachexia. JPEN J Parenter Enteral Nutr. 1992, 16: 50S-55S. 10.1177/014860719201600603.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/265/prepub
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.