Hypoxia increases membrane metallo-endopeptidase expression in a novel lung cancer ex vivo model – role of tumor stroma cells
- Katharina Leithner1,
- Christoph Wohlkoenig1,
- Elvira Stacher2,
- Jörg Lindenmann3,
- Nicole A Hofmann4, 5,
- Birgit Gallé6,
- Christian Guelly6,
- Franz Quehenberger7,
- Philipp Stiegler8,
- Freyja-Maria Smolle-Jüttner3,
- Sjaak Philipsen9,
- Helmut H Popper2,
- Andelko Hrzenjak1, 11,
- Andrea Olschewski10, 11 and
- Horst Olschewski1Email author
© Leithner et al.; licensee BioMed Central Ltd. 2014
Received: 1 July 2013
Accepted: 23 January 2014
Published: 25 January 2014
Hypoxia-induced genes are potential targets in cancer therapy. Responses to hypoxia have been extensively studied in vitro, however, they may differ in vivo due to the specific tumor microenvironment. In this study gene expression profiles were obtained from fresh human lung cancer tissue fragments cultured ex vivo under different oxygen concentrations in order to study responses to hypoxia in a model that mimics human lung cancer in vivo.
Non-small cell lung cancer (NSCLC) fragments from altogether 70 patients were maintained ex vivo in normoxia or hypoxia in short-term culture. Viability, apoptosis rates and tissue hypoxia were assessed. Gene expression profiles were studied using Affymetrix GeneChip 1.0 ST microarrays.
Apoptosis rates were comparable in normoxia and hypoxia despite different oxygenation levels, suggesting adaptation of tumor cells to hypoxia. Gene expression profiles in hypoxic compared to normoxic fragments largely overlapped with published hypoxia-signatures. While most of these genes were up-regulated by hypoxia also in NSCLC cell lines, membrane metallo-endopeptidase (MME, neprilysin, CD10) expression was not increased in hypoxia in NSCLC cell lines, but in carcinoma-associated fibroblasts isolated from non-small cell lung cancers. High MME expression was significantly associated with poor overall survival in 342 NSCLC patients in a meta-analysis of published microarray datasets.
The novel ex vivo model allowed for the first time to analyze hypoxia-regulated gene expression in preserved human lung cancer tissue. Gene expression profiles in human hypoxic lung cancer tissue overlapped with hypoxia-signatures from cancer cell lines, however, the elastase MME was identified as a novel hypoxia-induced gene in lung cancer. Due to the lack of hypoxia effects on MME expression in NSCLC cell lines in contrast to carcinoma-associated fibroblasts, a direct up-regulation of stroma fibroblast MME expression under hypoxia might contribute to enhanced aggressiveness of hypoxic cancers.
KeywordsHypoxia Tumor Expression array Prognosis
Survival following diagnosis of non-small cell lung cancer (NSCLC) is poor despite therapy . Hypoxia is typically present in solid tumors like lung cancer and is known to enhance tumor progression and therapy resistance . The effects of hypoxia are largely mediated by the hypoxia-inducible factors (HIFs) HIF-1α [3, 4] and HIF-2α . HIFs induce the expression of many different proteins that are involved in key functions of cancer cells, including cell survival, metabolic reprogramming, angiogenesis, invasion, and metastasis. Under normoxic conditions, HIFs are rapidly degraded, while under hypoxia they are stabilized [3, 4]. In addition to oxygen-dependent regulation, HIFs can be up-regulated by other mechanisms, e.g. growth factor induced pathways [3, 4]. The biological response of tumors to hypoxia is influenced by the interplay of neoplastic cancer cells and the surrounding stroma cells, e.g. cancer-associated fibroblasts (CAFs) . Ex vivo human cancer models based on the short-term culture of small tumor fragments or slices are suitable to study tumor responses within the natural in situ microenvironment, comprising a close contact between tumor cells and the accompanying stroma cells. Such models have been used e.g. for the study of drug effects in lung cancer  and other cancers [8, 9]. Here we used a human ex vivo lung cancer model involving culture of fresh tumor fragments in a hypoxic atmosphere to mimic in vivo tumor hypoxia and performed a comparative expression profiling study. We found that hypoxia led to overexpression of a stem-cell marker with elastase activity, membrane metallo-endopeptidase (MME), in tumor fragments, which was attributable to carcinoma-associated fibroblasts, not the neoplastic cancer cells.
Lung cancer fragments
Tumor tissue samples from 70 consecutive patients with NSCLC who were referred for surgical resection to the Division of Thoracic and Hyperbaric Surgery, Medical University of Graz, from May 2007 to May 2013, were included in the study. Patients with pre-operative chemotherapy were excluded from the study. Surgical specimens were dissected into small fragments using a razor blade and fragments were incubated in 35 mm Petri dishes (up to ten fragments per well) in 2 ml of DMEM/F-12 growth medium (Gibco, Carlsbad, CA) containing 10% fetal calf serum (Biowest Ltd, Ringmer, UK), 2 mM L-glutamine (Gibco), 100 U/ml penicillin, and 100 μg/ml streptomycin (Gibco). The study protocol was approved by the ethics review board of the Medical University of Graz. Signed informed consent was obtained from all patients prior to surgery.
The human NSCLC cell lines A549 and A427 were purchased from Cell Lines Service (Eppelheim, Germany) and cultured in DMEM/F-12 medium containing the supplements described above. The human NSCLC cell lines NCI-H23, NCI-H358, NCI-H1299, and NCI-H441 were purchased from American Type Culture Collection (ATCC, Manassas, VA) and cultured in RPMI (Gibco), supplemented with 10% fetal calf serum (Biowest) and antibiotics.
Carcinoma-associated fibroblasts (CAFs) were isolated from three fresh NSCLC samples as described  and cultured in DMEM supplemented with 10% fetal calf serum (Biowest) and antibiotics. CAFs were identified to be positive for vimentin and negative for cytokeratin using immunofluorescence. The purity of the cells was 97-99%. Human lung fibroblasts were cultured from donor lungs that could not be used for transplantation as previously described .
Fragments were cultured for three days at 37°C in ambient (21%) oxygen or 1% oxygen in the automated Xvivo System G300CL (BioSpherix, Lacona, NY). NSCLC cells or fibroblasts were plated into cell culture flasks at 13,000/cm2 and let attach, thereafter cells were cultured for three days in ambient oxygen or 1% oxygen as described above. Exposure to oxygen was controlled throughout the experiments in the hypoxic workstation.
The MTT assay (Chemicon, Billerica, MA) was performed on cultured fragments according to the manufacturer’s instructions. Briefly fragments were incubated in the MTT substrate solution for one hour and formazan was dissolved in isopropanol. After dissolving the formazan 100 μL of sample was analyzed on a colorimetric microplate reader at 570 nm. A549 cells were used as a positive control.
The assay (Hypoxyprobe™, HPI, Burlington, MA) was performed essentially according to the manufacturer’s instructions. Fragments were incubated for one or three days in hypoxia or normoxia. Thereafter fragments were treated with 100 μM pimonidazole HCl (HPI) in hypoxia in the closed Xvivo hypoxic working chamber (BioSpherix) or in normoxia and incubated for one hour, fixed and paraffin embedded. Bound pimonidazole was visualized using mouse monoclonal pimonidazole antibody (1:50 dilution, HPI).
RNA extraction and cDNA synthesis
Total RNA was extracted using the Qiagen RNeasy Mini kit (Qiagen, Hilden, Germany) and DNase digestion (Qiagen) according to the manufacturer’s instructions. RNA integrity was assessed using the Agilent 2100 Bioanalyzer and the Agilent RNA 6000 Nano Kit (Agilent, Palo Alto, CA). All samples exhibited a RIN (RNA Integrity Number) >5. Samples with RIN > 8 were eligible for microarray analysis. Total RNA (1 μg) was reverse transcribed using the RevertAid H Minus First Strand cDNA synthesis kit (Fermentas, Burlington, Canada).
Quantitative real-time PCR
For single gene quantitative polymerase chain reactions (PCR) the 7900 Real-Time PCR System (Applied Biosystems, Foster City, CA) was used. Gene expression assays (TaqMan® Gene Expression Assays, Applied Biosystems) suitable for this system were used for the detection of carbonic anhydrase IX, PPP1R3C, MME, KCTD11, FAM115C, and hexokinase 2. ACTB (ß-actin) was used as a reference gene. Primer data are indicated in Additional file 1: Table S2. The PCR was performed in 10 μl reactions containing cDNA (equal to 2.5 ng or 12.5 ng total RNA), 1× TaqMan® Gene Expression Mastermix (Applied Biosystems) and 1× TaqMan® Gene Expression Assay (Applied Biosystems). The mean threshold cycle (Ct) number of triplicate runs was used for data analysis. ΔCt was calculated by subtracting the Ct number of the gene of interest from that of the reference gene β-actin (ACTB). For calculation of differences between two groups, ∆Ct-values of the control group (normoxia) were substracted from ∆Ct-values of the treated group (hypoxia).
The microarray analysis was performed using GeneChip Human Gene 1.0 ST Arrays (Affymetrix, Santa Clara, CA). Manufacturer’s instructions were followed for the hybridization, washing, and scanning steps. Pre-labelled spike-in controls, unlabelled spike controls, and background probes were included in the analysis. All the microarray data are available at Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/; accession number GSE30979).
Processing of microarray data
Statistical analysis of the microarray data was performed using Partek Genomic Suite Software (Partek, St. Louis, MO). RMA (Robust Multi Chip Analysis) background correction of raw microarray data and normalization of expression values were performed using Partek Genomic Suite Software (Partek). Fold-changes of expression values were calculated as the ratio of the mean RMA corrected expression value in the hypoxic group to the normoxic group. Fold-change values <1 were converted to the negative of the inverse ratio. Hypoxic and normoxic samples were compared using the paired Student’s t-test. The false discovery rate (FDR) was set to 5% to correct for multiple testing. In the case of subgroup analyses, the threshold was set to P<0.005. A gene was considered modulated when at least one of the corresponding probe sets showed significantly different expression levels after correction for multiple testing with a minimal two-fold change.
Meta-analysis of lung cancer transcriptome studies
Expression values for the genes of interest were obtained from four eligible lung cancer datasets published at Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/). Details on data processing and patient characteristics are reported at GEO and in the cited literature. Details on data retrieval are indicated in Additional file 1.
Meta-analysis of the effect of MME on patient survival after surgery was performed with a proportional hazards model with Gaussian random effects [12, 13] using the package coxme 2.1-3 of R 2.13.2 statistical software (http://www.r-project.org). For details see Additional file 1. All other data were compiled and analyzed with the SPSS software package, version 18.0 (Chicago, IL). Group differences were calculated with the paired Student’s t-test, one-sample Student’s t-test, Mann–Whitney-U test, or Wilcoxon signed rank test as applicable. P-values smaller than 0.05 were considered significant.
Apoptosis and hypoxia markers in NSCLC fragments
HIF-1α and HIF-2α immunohistochemistry was performed in NSCLC fragments cultured for three days under normoxia or hypoxia (Figure 1B). HIF-1α was localized predominantly in the nucleus, while HIF-2α was found in the cytoplasm. Both, cytoplasmic and nuclear localization of HIF-1α and HIF-2α, have been reported [14, 15]. Hypoxic fragments displayed more pronounced staining for HIF-1α than normoxic fragments, though the difference was significant only in stroma cells, not in tumor cells (Figure 1B). For HIF-2α no difference between fragments cultured in hypoxia or normoxia was found, neither in tumor cells, nor in stroma cells (Figure 1B). Next we assessed the presence of hypoxia in cultured fragments using pimonidazole. Figure 1C shows examples of NSCLC fragments cultured in normoxia or hypoxia for one and three days. Pimonidazole was bound almost to the entire hypoxic fragments, while only focal pimonidazole binding occurred in normoxic fragments, obviously due to diminished oxygen concentrations in central fragment areas. In several hypoxic fragments some cells showed higher pimonidazole binding than others (Figure 1C), which might be caused by a different content of redox enzymes (J. Raleigh, Hypoxyprobe Inc. and UNC School of Medicine, Chapel Hill, North Carolina, USA, personal communication) or due to other cell-related causes, such as differences in pimonidazole uptake or pH . Expression of the HIF-1α target carbonic anhydrase IX (CA IX), which was shown to be linked to hypoxia in NSCLCs in vivo, was analyzed by quantitative PCR. CA IX mRNA levels were significantly higher in hypoxic fragments compared to normoxic fragments (Figure 1D). Taken together, NSCLC fragments remained viable for the duration of the experiments and hypoxia markers were increased under hypoxic treatment.
Gene regulation by hypoxia in NSCLC fragments
Genes up-regulated by hypoxia
Protein phosphatase 1, regulatory (inhibitor) subunit 3C
Glucan (1,4-alpha-), branching enzyme 1
Lysine (K)-specific demethylase 3A
Potassium channel tetramerisation domain containing 11
Family with sequence similarity 115, member C
BCL2/adenovirus E1B 19 kDa interacting protein 3
DnaJ (Hsp40) homolog, subfamily B, member 9
Chromosome 7 open reading frame 68
Heat shock protein 70 kDa family, member 13
Leucine rich repeat containing 49
Triggering receptor expressed on myeloid cells 1
Cerebellar degeneration-related protein 1
DNA-damage-inducible transcript 3
Chromosome 1 open reading frame 129
FIC domain containing
Fatty acid binding protein 3
ABI family, member 3 (NESH) binding protein
SEC11 homolog C
Interestingly, the overall impact of hypoxia on gene expression was lower than the impact of histology or inter-patient variability (Additional file 3: Figure S1). Normoxic and hypoxic fragments derived from each patient clustered together significantly in 9 of 10 patients in pvclust analysis (Additional file 3: Figure S1). Both clusters on the top of the hierarchy were significant in pvclust analysis. One cluster contained four squamous cell carcinomas, the other cluster contained all adenocarcinomas and one squamous cell carcinoma (Additional file 3: Figure S1).
Expression of hypoxia-regulated genes in NSCLC cells and carcinoma-associated fibroblasts (CAFs)
We further analyzed the expression of MME under hypoxia in both, NSCLC cell lines and fibroblasts, which are the predominant cell type in lung cancer stroma, using quantitative PCR. PPP1R3C, KCTD11, FAM115C, and HK2, a well-known hypoxia-regulated gene, were up-regulated by hypoxia in a panel of NSCLC cell lines to variable degrees, while MME mRNA showed no increase in expression under hypoxia in any of the cell lines (Figure 3C). On the contrary in carcinoma-associated fibroblasts (CAFs) from NSCLC and, to a lesser extent, in primary lung fibroblasts MME mRNA was significantly up-regulated by hypoxia (Figure 3D).
MME expression is an adverse prognostic factor in lung adenocarcinoma patients
Next, we examined whether expression of the four hypoxia genes was associated with survival in patients with NSCLC. Due to the relatively short observation period in our patient cohort, we used large published microarray datasets containing gene expression data linked to clinical and prognostic information in NSCLC patients. The Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) is one of the largest microarray databases. A search for GEO datasets/series using the search criteria „lung cancer 50:500[Number of Samples]” yielded 84 results (status June 2011). Of these 84 datasets/series, 68 contained expression profiling data. Four of these series included expression data of a minimal number of 50 NSCLC patients treated by surgery with linked information on survival, GSE11969 , GSE13213 , GSE14814 , and GSE19188 . Altogether 342 patients were included in the meta-analysis.
Identifying hypoxia-regulated genes may promote understanding of the molecular response to hypoxic stress in cancers. Changes in gene expression in hypoxic cancer cells have been studied extensively in vitro. However, hypoxia-responses in vivo may differ from the in vitro situation due to the complex tumor microenvironment. In fact, hypoxia activates tumor promoting stroma cells and HIF-1α has been identified as the major driver of tumor-stroma “co-evolution” . Here we studied hypoxia-induced gene expression experimentally in human cancer tissue in its preserved 3D-structure. In this fragment model the tissue contains both tumor and stroma cells and mimics the in vivo situation.
The model has several advantages compared to in vitro cancer cell lines. The tumor cells remain in contact with their original tumor microenvironment (stroma cells, extracellular matrix), the 3D-morphology is preserved, and inter-patient variability is taken into account by using material derived from different patients. The major limitation of our study is that the exact oxygen concentration could only be controlled on the surface of the fragments. Inside the tumor fragments there are supposed to be oxygen gradients, depending on the size and composition of the tissue fragment. The size of normoxic and hypoxic fragments did not differ in our study. In fact, using pimonidazole staining, fragments cultured in hypoxia were found to be entirely hypoxic, while only a core of hypoxia was found in fragments cultured in normoxia. In addition to pimonidazole, also other major hypoxia-markers were significantly increased in the hypoxic fragments, such as HIF-1α and CA IX. However, HIF-2α, which is known to be stabilized by hypoxia similarly to HIF-1α, was expressed only at low levels, both in normoxia and hypoxia, and was not elevated in hypoxic fragments. Different co-activators and different kinetics of activation under hypoxia  might play a role. This indicates that the difference in oxygen concentration was preserved despite the expected oxygen gradients inside the fragments. Furthermore the oxygen decline is supposed to occur in both, normoxic and hypoxic fragments. Thus our approach is feasible to study differential gene expression under high and low oxygen concentrations.
Apoptosis rates were comparable in NSCLC fragments cultured in 1% O2 or normoxia for three days. This agrees with our previous study where we showed that hypoxia-induced adaptation and cisplatin-resistance are reversible in lung cancer cells and occur without hypoxia-induced cell death and selection . In an attempt to identify common hypoxia-regulated genes, Ortiz-Barahona et al.  identified 17 genes consistently up-regulated by hypoxia, hypoxia-mimetics, or HIF-1α using a meta-analysis of expression data from 16 GEO datasets. Of these 17, mostly well-known hypoxia-regulated genes, 65% appear among the significantly regulated genes in our study (after correction for multiple testing). When we compared a hypoxia signature found to be prognostically relevant in many cancers (the “hypoxia metagene” ) with our hypoxia profile, we also found a considerable overlap. Approximately half of the top-ranked hypoxia-induced genes with prognostic relevance identified by Buffa et al.  were significantly up-regulated by hypoxia in our study.
Four genes were significantly up-regulated by hypoxia in both adenocarcinoma and squamous cell carcinoma fragments in our setting. We confirmed the differential expression of the four overlapping hypoxia genes under hypoxia in an independent validation set using quantitative PCR (qPCR). Also the well-established hypoxia-responsive gene HK2, which phosphorylates glucose and thus contributes to the glycolytic flux in cancer cells, was significantly up-regulated by hypoxia in the fragments, both in the microarray analysis and by qPCR.
The four hypoxia-genes identified in our study have been found to be up-regulated by hypoxia in several microarray studies, however these findings were not validated e.g. by qPCR [27–32]. To the best of our knowledge, validated data on hypoxia-regulation of the four hypoxia-regulated genes exist for PPP1R3C (up-regulated in MCF-7 breast cancer cells)  and on MME, which was shown to be up-regulated in primary rat astrocytes  and down-regulated in pulmonary artery smooth muscle cells , human neuroblastoma cells , rat neurons , and mouse neurons . Cobalt chloride, a hypoxia mimetic, was shown to reduce MME expression in prostate cancer cell lines , and human umbilical vein endothelial cells . In addition, exposure of rats and mice to a hypoxic atmosphere led to down-regulation of MME expression [38, 39].
In our study we found MME localized to neoplastic tumor cells, but also to stroma cells in fresh NSCLC tissue, which is in line with published data [40, 41]. The observed up-regulation of MME under hypoxia in NSCLC fragments might thus be attributable to tumor cells or stroma cells, or both. While the hypoxic regulation of KCTD11, FAM115C, PPP1R3C and HK2 was also observed to a variable degree in a panel of NSCLC cell lines cultured as a monolayer, MME was not regulated by hypoxia in the cell lines in our study. Fibroblasts are the predominant cell type in lung cancer stroma . When we studied MME mRNA in CAFs we found a significant induction by hypoxia. A similar effect was found in normal lung fibroblasts, however to a lesser extent. The exact mechanism of MME regulation by hypoxia in fibroblasts remains to be elucidated. The proximal promoter regions of the different MME splice variants have been shown to harbour binding sites for the transcription factors Sp1, PEA3 and PU.1 . PEA3 (also known as E1AF and ETV4) is a member of the Ets-family of transcription factors. PEA3 was shown enhance cancer metastasis . Recently, PEA3 has been shown to interact with HIF-1α . This might at least partially be responsible for the observed effect of hypoxia on MME expression.
MME, which is identical to common acute leukemia antigen (CALLA), is a 90–110 kDa zinc binding cell surface peptidase, which cleaves small peptides, such as atrial natriuretic peptide, substance P, endothelin-1, and bombesin (for review see [46, 47]). It also possesses elastase activity . MME is a membrane-bound protein, however, as was recently shown, MME can be released to the microenvironment of cells in exosomes . MME is expressed in a variety of non-malignant and malignant tissues (for review see [46, 47]) including lung cancer [40, 50–52]. In small-cell lung carcinoma (SCLC) cells, bombesin-like peptides, substrates for MME, are autocrine growth factors. Cleaving these peptides by recombinant MME has been shown to inhibit SCLC cell proliferation [53, 54]. In NSCLC cells, recombinant MME inhibited tumor cell proliferation in vitro, but only at very high concentrations and after long exposure . On the contrary, MME inhibitors have been found to decrease cell proliferation in the airway wall in response to cigarette smoke in rats . While the role of MME in neoplastic tumor cells is still unclear, several reports suggest that stroma cell MME expression plays a role in tumor progression. MME-positive stroma cells, including mesenchymal stem cells and fibroblasts, have been shown to promote tumor aggressiveness and metastasis [56, 57]. Elastin is degraded by MME , which might facilitate tumor and/or stroma cell invasion.
In order to analyze, whether levels of the common hypoxia-genes identified in our study are associated with overall survival in NSCLC patients we used all eligible studies deposited in one of the largest microarray depositories, the GEO database. We were able to show that MME expression is a highly significant, independent adverse prognostic factor in surgically treated lung adenocarcinoma patients in multivariate analysis involving tumor stage and MME status. No association was found in the subgroup of non-adenocarcinoma patients. The reason for the different results in the histological subgroups is unknown, however, lung adenocarcinomas have been shown to possess more elastin than squamous cell carcinomas . Since the largest study with 116 adenocarcinoma patients (GSE13213) contained only adenocarcinomas, a study-bias cannot be excluded.
To the best of our knowledge, three other studies examined the association of MME expression and survival in lung cancer [40, 41, 51]. All studies are immunohistochemical studies. In a study by Kristiansen et al.  in 114 NSCLC patients no association of MME immunostaining and survival was found. Only neoplastic cancer cells were evaluated in that study (G. Kristiansen, personal communication). In a recent study by Ono et al.  on 142 stage I squamous cell lung carcinoma patients MME expression was examined in tumor cells and stroma cells separately. Patients with low MME expression in stroma or in tumor cells survived slightly longer, but the differences were not significant. In a study by Gurel et al.  MME expression was studied in tumor cells and stroma cells in 66 patients with NSCLC using immunohistochemistry. In the squamous cell carcinoma subgroup high tumor cell and stroma cell MME expression were both associated with poor overall survival. In non-squamous cell NSCLC (35 patients, adenocarcinoma, large cell carcinoma, sarcomatoid carcinoma, and mixed types) the opposite association was found. No stroma cell MME expression was found in that subgroup . The low number of patients may make the interpretation of these results difficult.
All of these studies were immunohistochemical studies, in contrast to our MME mRNA based survival analysis. Since MME may be excreted in exosomes , which has been shown e.g. for mesenchymal stem cells , the question arises, whether MME was excreted and then lost during conventional tissue fixation and immunohistochemistry. With immunohistochemistry thus the extent of MME expression in cancer tissue may be underestimated. This is supported by the fact that CAFs isolated from all three NSCLCs expressed MME mRNA in our study, while in the studies mentioned above high MME staining in stroma cells was found only in 11% to 19% of cases. This underestimation may partly explain the lack of association between MME expression and worse prognosis in the mentioned studies, as opposed to our mRNA based study.
Additional studies examined the expression of MME in combination with other factors and survival. In the study by Tokuhara et al.  132 NSCLC patients were grouped according to their tumor MME mRNA and aminopeptidase N mRNA expression. Patients assigned to the group with high MME and low aminopeptidase N mRNA showed significantly improved survival. No analysis on MME expression alone was performed. Tumor tissue samples were selected to contain primarily cancer cells in that study. In a study by Navab et al.  MME was among a subset of eleven genes identified to be up-regulated in cancer associated fibroblasts, forming a prognostic gene-expression signature in NSCLC.
The novel ex vivo model allowed for the first time to analyze hypoxia-regulated gene expression in preserved human lung cancer tissue. The study shows that gene expression profiles in human hypoxic lung cancer tissue overlap with hypoxia-signatures from cancer cell lines, however, MME was identified as a novel hypoxia-induced gene in lung cancer. Despite the advantages of ex vivo tissue culture, cell monolayers still appear to be the method of choice to study mechanisms of adaptation of individual cell types to hypoxia, since the oxygen concentration can be controlled only on the surface of such three-dimensional structures. Thus we analyzed expression of the hypoxia-regulated genes identified in the NSCLC fragments in different NSCLC cell lines and primary CAFs isolated from NSCLC tissue. We show that MME expression is up-regulated by hypoxia in CAFs, not in NSCLC cells. High global levels of MME mRNA in NSCLC tissue were shown in our study to predict poor survival. A direct effect of hypoxia on stromal fibroblast MME expression might thus contribute to enhanced aggressiveness of hypoxic cancers.
Non-small cell lung cancer
Common acute leukemia antigen
- CA IX:
Carbonic anhydrase IX
False discovery rate
Protein phosphatase 1 regulatory subunit 3C
Potassium channel tetramerisation domain containing 11
Family with sequence similarity 115 member C
We are grateful to Alexandra Bertsch and Elisabeth Pöllitzer, Division of Pulmonology, and Karin Wagner, Core Facility Molecular Biology, Medical University of Graz, Austria, for excellent technical assistance. We thank Dr. Armin Frille, University Hospital Leipzig, Leipzig, Germany, for helpful discussions and for his help with experimental work. We thank Dr. Grazyna Kwapiszewska, Dr. Slaven Crnkovic, and all other members of the Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria, for providing lung fibroblasts and for their support. We would like to thank Eugenia Lamont BA, Department of Surgery, Medical University of Graz, Austria, for critical reading of the manuscript. We are grateful to Dr. Franz Gollowitsch, Institute of Pathology, Medical University of Graz, Austria, for his help. We would like to thank Dr. Chang Qi Zhu, Princess Margaret Hospital/Ontario Cancer Institute, Toronto, Canada, for providing additional clinical information for the meta-analysis and Prof. Dr. Glen Kristiansen, University Hospital of Bonn, Bonn, Germany, for his comment on MME immunohistochemistry in NSCLC tissues. The valuable advice by Prof. Dr. Adrian L. Harris, Weatherall Institute of Molecular Medicine, Oxford, UK, is highly appreciated.
The study was supported by funds of the Oesterreichische Nationalbank (Anniversary Fund, project number 12713 to HO) and the Start Funding Grant of the Medical University of Graz (ASO212009100 to KL). CW was supported by the Medical University of Graz (PhD Program Molecular Medicine). SP was supported by the Netherlands Genomics Initiative (NGI) through the Cancer Genomics Centre and the Netherlands Consortium for Systems Biology. The content is solely in the responsibility of the authors.
- Manegold C, Thatcher N: Survival improvement in thoracic cancer: progress from the last decade and beyond. Lung Cancer. 2007, 57 (Suppl 2): S3-S5.View ArticlePubMedGoogle Scholar
- Höckel M, Vaupel P: Tumor hypoxia: definitions and current clinical, biologic, and molecular aspects. J Natl Cancer Inst. 2001, 93 (4): 266-276. 10.1093/jnci/93.4.266.View ArticlePubMedGoogle Scholar
- Semenza GL: Defining the role of hypoxia-inducible factor 1 in cancer biology and therapeutics. Oncogene. 2010, 29 (5): 625-634. 10.1038/onc.2009.441.View ArticlePubMedGoogle Scholar
- Harris AL: Hypoxia - a key regulatory factor in tumour growth. Nat Rev Cancer. 2002, 2 (1): 38-47. 10.1038/nrc704.View ArticlePubMedGoogle Scholar
- Shimoda LA, Semenza GL: HIF and the lung: role of hypoxia-inducible factors in pulmonary development and disease. Am J Respir Crit Care Med. 2011, 183 (2): 152-156. 10.1164/rccm.201009-1393PP.View ArticlePubMedPubMed CentralGoogle Scholar
- Giaccia AJ, Schipani E: Role of carcinoma-associated fibroblasts and hypoxia in tumor progression. Curr Top Microbiol Immunol. 2010, 345: 31-45.PubMedGoogle Scholar
- Schmid JO, Dong M, Haubeiss S, Friedel G, Bode S, Grabner A, Ott G, Murdter TE, Oren M, Aulitzky WE, van der Kuip H: Cancer cells Cue the p53 response of cancer-associated fibroblasts to cisplatin. Cancer Res. 2012, 72 (22): 5824-5832. 10.1158/0008-5472.CAN-12-1201.View ArticlePubMedGoogle Scholar
- Hougardy BM, Reesink-Peters N, van den Heuvel FA, ten Hoor KA, Hollema H, de Vries EG, de Jong S, van der Zee AG: A robust ex vivo model for evaluation of induction of apoptosis by rhTRAIL in combination with proteasome inhibitor MG132 in human premalignant cervical explants. Int J Cancer. 2008, 123 (6): 1457-1465. 10.1002/ijc.23684.View ArticlePubMedGoogle Scholar
- Kirby TO, Rivera A, Rein D, Wang M, Ulasov I, Breidenbach M, Kataram M, Contreras JL, Krumdieck C, Yamamoto M, Rots MG, Haisma HJ, Alvarez RD, Mahasreshti PJ, Curiel DT: A novel ex vivo model system for evaluation of conditionally replicative adenoviruses therapeutic efficacy and toxicity. Clin Cancer Res. 2004, 10 (24): 8697-8703. 10.1158/1078-0432.CCR-04-1166.View ArticlePubMedGoogle Scholar
- Navab R, Strumpf D, Bandarchi B, Zhu CQ, Pintilie M, Ramnarine VR, Ibrahimov E, Radulovich N, Leung L, Barczyk M, Panchal D, To C, Yun JJ, Der S, Shepherd FA, Jurisica I, Tsao MS: Prognostic gene-expression signature of carcinoma-associated fibroblasts in non-small cell lung cancer. Proc Natl Acad Sci USA. 2011, 108 (17): 7160-7165. 10.1073/pnas.1014506108.View ArticlePubMedPubMed CentralGoogle Scholar
- Avcuoglu S, Wygrecka M, Marsh LM, Gunther A, Seeger W, Weissmann N, Fink L, Morty RE, Kwapiszewska G: Neurotrophic tyrosine kinase receptor B/neurotrophin 4 signaling axis is perturbed in clinical and experimental pulmonary fibrosis. Am J Respir Cell Mol Biol. 2011, 45 (4): 768-780. 10.1165/rcmb.2010-0195OC.View ArticlePubMedGoogle Scholar
- Michiels S, Baujat B, Mahe C, Sargent DJ, Pignon JP: Random effects survival models gave a better understanding of heterogeneity in individual patient data meta-analyses. J Clin Epidemiol. 2005, 58 (3): 238-245. 10.1016/j.jclinepi.2004.08.013.View ArticlePubMedGoogle Scholar
- Ripatti S, Palmgren J: Estimation of multivariate frailty models using penalized partial likelihood. Biometrics. 2000, 56 (4): 1016-1022. 10.1111/j.0006-341X.2000.01016.x.View ArticlePubMedGoogle Scholar
- Talks KL, Turley H, Gatter KC, Maxwell PH, Pugh CW, Ratcliffe PJ, Harris AL: The expression and distribution of the hypoxia-inducible factors HIF-1alpha and HIF-2alpha in normal human tissues, cancers, and tumor-associated macrophages. Am J Pathol. 2000, 157 (2): 411-421. 10.1016/S0002-9440(10)64554-3.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersen S, Eilertsen M, Donnem T, Al-Shibli K, Al-Saad S, Busund LT, Bremnes RM: Diverging prognostic impacts of hypoxic markers according to NSCLC histology. Lung Cancer. 2011, 72 (3): 294-302. 10.1016/j.lungcan.2010.10.006.View ArticlePubMedGoogle Scholar
- Kleiter MM, Thrall DE, Malarkey DE, Ji X, Lee DY, Chou SC, Raleigh JA: A comparison of oral and intravenous pimonidazole in canine tumors using intravenous CCI-103 F as a control hypoxia marker. Int J Radiat Oncol Biol Phys. 2006, 64 (2): 592-602. 10.1016/j.ijrobp.2005.09.010.View ArticlePubMedGoogle Scholar
- Le QT, Chen E, Salim A, Cao H, Kong CS, Whyte R, Donington J, Cannon W, Wakelee H, Tibshirani R, Mitchell JD, Richardson D, O’Byrne KJ, Koong AC, Giaccia AJ: An evaluation of tumor oxygenation and gene expression in patients with early stage non-small cell lung cancers. Clin Cancer Res. 2006, 12 (5): 1507-1514. 10.1158/1078-0432.CCR-05-2049.View ArticlePubMedGoogle Scholar
- Takeuchi T, Tomida S, Yatabe Y, Kosaka T, Osada H, Yanagisawa K, Mitsudomi T, Takahashi T: Expression profile-defined classification of lung adenocarcinoma shows close relationship with underlying major genetic changes and clinicopathologic behaviors. J Clin Oncol. 2006, 24 (11): 1679-1688. 10.1200/JCO.2005.03.8224.View ArticlePubMedGoogle Scholar
- Tomida S, Takeuchi T, Shimada Y, Arima C, Matsuo K, Mitsudomi T, Yatabe Y, Takahashi T: Relapse-related molecular signature in lung adenocarcinomas identifies patients with dismal prognosis. J Clin Oncol. 2009, 27 (17): 2793-2799. 10.1200/JCO.2008.19.7053.View ArticlePubMedGoogle Scholar
- Zhu CQ, Ding K, Strumpf D, Weir BA, Meyerson M, Pennell N, Thomas RK, Naoki K, Ladd-Acosta C, Liu N, Pintilie M, Der S, Seymour L, Jurisica I, Shepherd FA, Tsao MS: Prognostic and predictive gene signature for adjuvant chemotherapy in resected non-small-cell lung cancer. J Clin Oncol. 2010, 28 (29): 4417-4424. 10.1200/JCO.2009.26.4325.View ArticlePubMedPubMed CentralGoogle Scholar
- Hou J, Aerts J, den Hamer B, van Ijcken W, den Bakker M, Riegman P, van der Leest C, van der Spek P, Foekens JA, Hoogsteden HC, Grosveld F, Philipsen S: Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS One. 2010, 5 (4): e10312-10.1371/journal.pone.0010312.View ArticlePubMedPubMed CentralGoogle Scholar
- Lisanti MP, Martinez-Outschoorn UE, Chiavarina B, Pavlides S, Whitaker-Menezes D, Tsirigos A, Witkiewicz A, Lin Z, Balliet R, Howell A, Sotgia F: Understanding the “lethal” drivers of tumor-stroma co-evolution: emerging role(s) for hypoxia, oxidative stress and autophagy/mitophagy in the tumor micro-environment. Cancer Biol Ther. 2010, 10 (6): 537-542. 10.4161/cbt.10.6.13370.View ArticlePubMedPubMed CentralGoogle Scholar
- Loboda A, Jozkowicz A, Dulak J: HIF-1 versus HIF-2–is one more important than the other?. Vascul Pharmacol. 2012, 56 (5–6): 245-251.View ArticlePubMedGoogle Scholar
- Wohlkoenig C, Leithner K, Deutsch A, Hrzenjak A, Olschewski A, Olschewski H: Hypoxia-induced cisplatin resistance is reversible and growth rate independent in lung cancer cells. Cancer Lett. 2011, 308 (2): 134-143. 10.1016/j.canlet.2011.03.014.View ArticlePubMedGoogle Scholar
- Ortiz-Barahona A, Villar D, Pescador N, Amigo J, del Peso L: Genome-wide identification of hypoxia-inducible factor binding sites and target genes by a probabilistic model integrating transcription-profiling data and in silico binding site prediction. Nucleic Acids Res. 2010, 38 (7): 2332-2345. 10.1093/nar/gkp1205.View ArticlePubMedPubMed CentralGoogle Scholar
- Buffa FM, Harris AL, West CM, Miller CJ: Large meta-analysis of multiple cancers reveals a common, compact and highly prognostic hypoxia metagene. Br J Cancer. 2010, 102 (2): 428-435. 10.1038/sj.bjc.6605450.View ArticlePubMedPubMed CentralGoogle Scholar
- Koritzinsky M, Seigneuric R, Magagnin MG, van den Beucken T, Lambin P, Wouters BG: The hypoxic proteome is influenced by gene-specific changes in mRNA translation. Radiother Oncol. 2005, 76 (2): 177-186. 10.1016/j.radonc.2005.06.036.View ArticlePubMedGoogle Scholar
- Mole DR, Blancher C, Copley RR, Pollard PJ, Gleadle JM, Ragoussis J, Ratcliffe PJ: Genome-wide association of hypoxia-inducible factor (HIF)-1alpha and HIF-2alpha DNA binding with expression profiling of hypoxia-inducible transcripts. J Biol Chem. 2009, 284 (25): 16767-16775. 10.1074/jbc.M901790200.View ArticlePubMedPubMed CentralGoogle Scholar
- Winter SC, Buffa FM, Silva P, Miller C, Valentine HR, Turley H, Shah KA, Cox GJ, Corbridge RJ, Homer JJ, Musgrove B, Slevin N, Sloan P, Price P, West CM, Harris AL: Relation of a hypoxia metagene derived from head and neck cancer to prognosis of multiple cancers. Cancer Res. 2007, 67 (7): 3441-3449. 10.1158/0008-5472.CAN-06-3322.View ArticlePubMedGoogle Scholar
- Loewen N, Chen J, Dudley VJ, Sarthy VP, Mathura JR: Genomic response of hypoxic Müller cells involves the very low density lipoprotein receptor as part of an angiogenic network. Exp Eye Res. 2009, 88 (5): 928-937. 10.1016/j.exer.2008.11.037.View ArticlePubMedGoogle Scholar
- Costello CM, Howell K, Cahill E, McBryan J, Konigshoff M, Eickelberg O, Gaine S, Martin F, McLoughlin P: Lung-selective gene responses to alveolar hypoxia: potential role for the bone morphogenetic antagonist gremlin in pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol. 2008, 295 (2): L272-L284. 10.1152/ajplung.00358.2007.View ArticlePubMedGoogle Scholar
- Koklanaris N, Nwachukwu JC, Huang SJ, Guller S, Karpisheva K, Garabedian M, Lee MJ: First-trimester trophoblast cell model gene response to hypoxia. Am J Obstet Gynecol. 2006, 194 (3): 687-693. 10.1016/j.ajog.2006.01.067.View ArticlePubMedGoogle Scholar
- Shen GM, Zhang FL, Liu XL, Zhang JW: Hypoxia-inducible factor 1-mediated regulation of PPP1R3C promotes glycogen accumulation in human MCF-7 cells under hypoxia. FEBS Lett. 2010, 584 (20): 4366-4372. 10.1016/j.febslet.2010.09.040.View ArticlePubMedGoogle Scholar
- Fisk L, Nalivaeva NN, Boyle JP, Peers CS, Turner AJ: Effects of hypoxia and oxidative stress on expression of neprilysin in human neuroblastoma cells and rat cortical neurones and astrocytes. Neurochem Res. 2007, 32 (10): 1741-1748. 10.1007/s11064-007-9349-2.View ArticlePubMedGoogle Scholar
- Wick MJ, Buesing EJ, Wehling CA, Loomis ZL, Cool CD, Zamora MR, Miller YE, Colgan SP, Hersh LB, Voelkel NF, Dempsey EC: Decreased neprilysin and pulmonary vascular remodeling in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2011, 183 (3): 330-340. 10.1164/rccm.201002-0154OC.View ArticlePubMedGoogle Scholar
- Wang Z, Yang D, Zhang X, Li T, Li J, Tang Y, Le W: Hypoxia-induced down-regulation of neprilysin by histone modification in mouse primary cortical and hippocampal neurons. PLoS One. 2011, 6 (4): e19229-10.1371/journal.pone.0019229.View ArticlePubMedPubMed CentralGoogle Scholar
- Mitra R, Chao OS, Nanus DM, Goodman OB: Negative regulation of NEP expression by hypoxia. Prostate. 2013, 73 (7): 706-714. 10.1002/pros.22613.View ArticlePubMedGoogle Scholar
- Dempsey EC, Wick MJ, Karoor V, Barr EJ, Tallman DW, Wehling CA, Walchak SJ, Laudi S, Le M, Oka M, Majka S, Cool CD, Fagan KA, Klemm DJ, Hersh LB, Gerard NP, Gerard C, Miller YE: Neprilysin null mice develop exaggerated pulmonary vascular remodeling in response to chronic hypoxia. Am J Pathol. 2009, 174 (3): 782-796. 10.2353/ajpath.2009.080345.View ArticlePubMedPubMed CentralGoogle Scholar
- Carpenter TC, Stenmark KR: Hypoxia decreases lung neprilysin expression and increases pulmonary vascular leak. Am J Physiol Lung Cell Mol Physiol. 2001, 281 (4): L941-L948.PubMedGoogle Scholar
- Gurel D, Kargi A, Karaman I, Onen A, Unlu M: CD10 expression in epithelial and stromal cells of non-small cell lung carcinoma (NSCLC): a clinic and pathologic correlation. Pathol Oncol Res. 2012, 18 (2): 153-160. 10.1007/s12253-011-9421-8.View ArticlePubMedGoogle Scholar
- Ono S, Ishii G, Nagai K, Takuwa T, Yoshida J, Nishimura M, Hishida T, Aokage K, Fujii S, Ikeda N, Ochiai A: Podoplanin-positive cancer-associated fibroblasts could have prognostic value independent of cancer cell phenotype in stage I lung squamous cell carcinoma: usefulness of combining analysis of both cancer cell phenotype and cancer-associated fibroblast phenotype. Chest. 2013, 143 (4): 963-970. 10.1378/chest.12-0913.View ArticlePubMedGoogle Scholar
- Bremnes RM, Donnem T, Al-Saad S, Al-Shibli K, Andersen S, Sirera R, Camps C, Marinez I, Busund LT: The role of tumor stroma in cancer progression and prognosis: emphasis on carcinoma-associated fibroblasts and non-small cell lung cancer. J Thorac Oncol. 2011, 6 (1): 209-217. 10.1097/JTO.0b013e3181f8a1bd.View ArticlePubMedGoogle Scholar
- Ishimaru F, Shipp MA: Analysis of the human CD10/neutral endopeptidase 24.11 promoter region: two separate regulatory elements. Blood. 1995, 85 (11): 3199-3207.PubMedGoogle Scholar
- de Launoit Y, Baert JL, Chotteau-Lelievre A, Monte D, Coutte L, Mauen S, Firlej V, Degerny C, Verreman K: The Ets transcription factors of the PEA3 group: transcriptional regulators in metastasis. Biochim Biophys Acta. 2006, 1766 (1): 79-87.PubMedGoogle Scholar
- Wollenick K, Hu J, Kristiansen G, Schraml P, Rehrauer H, Berchner-Pfannschmidt U, Fandrey J, Wenger RH, Stiehl DP: Synthetic transactivation screening reveals ETV4 as broad coactivator of hypoxia-inducible factor signaling. Nucleic Acids Res. 2012, 40 (5): 1928-1943. 10.1093/nar/gkr978.View ArticlePubMedGoogle Scholar
- Sumitomo M, Shen R, Nanus DM: Involvement of neutral endopeptidase in neoplastic progression. Biochim Biophys Acta. 2005, 1751 (1): 52-59. 10.1016/j.bbapap.2004.11.001.View ArticlePubMedGoogle Scholar
- Maguer-Satta V, Besancon R, Bachelard-Cascales E: Concise review: neutral endopeptidase (CD10): a multifaceted environment actor in stem cells, physiological mechanisms, and cancer. Stem Cells. 2011, 29 (3): 389-396. 10.1002/stem.592.View ArticlePubMedGoogle Scholar
- Morisaki N, Moriwaki S, Sugiyama-Nakagiri Y, Haketa K, Takema Y, Imokawa G: Neprilysin is identical to skin fibroblast elastase: its role in skin aging and UV responses. J Biol Chem. 2010, 285 (51): 39819-39827. 10.1074/jbc.M110.161547.View ArticlePubMedPubMed CentralGoogle Scholar
- Mazurov D, Barbashova L, Filatov A: Tetraspanin protein CD9 interacts with metalloprotease CD10 and enhances its release via exosomes. FEBS J. 2013, 280 (5): 1200-1213. 10.1111/febs.12110.View ArticlePubMedGoogle Scholar
- Tokuhara T, Adachi M, Hashida H, Ishida H, Taki T, Higashiyama M, Kodama K, Tachibana S, Sasaki S, Miyake M: Neutral endopeptidase/CD10 and aminopeptidase N/CD13 gene expression as a prognostic factor in non-small cell lung cancer. Jpn J Thorac Cardiovasc Surg. 2001, 49 (8): 489-496. 10.1007/BF02919543.View ArticlePubMedGoogle Scholar
- Kristiansen G, Schlüns K, Yongwei Y, Dietel M, Petersen I: CD10 expression in non-small cell lung cancer. Anal Cell Pathol. 2002, 24 (1): 41-46.View ArticlePubMedPubMed CentralGoogle Scholar
- Cohen AJ, Bunn PA, Franklin W, Magill-Solc C, Hartmann C, Helfrich B, Gilman L, Folkvord J, Helm K, Miller YE: Neutral endopeptidase: variable expression in human lung, inactivation in lung cancer, and modulation of peptide-induced calcium flux. Cancer Res. 1996, 56 (4): 831-839.PubMedGoogle Scholar
- Shipp MA, Tarr GE, Chen CY, Switzer SN, Hersh LB, Stein H, Sunday ME, Reinherz EL: CD10/neutral endopeptidase 24.11 hydrolyzes bombesin-like peptides and regulates the growth of small cell carcinomas of the lung. Proc Natl Acad Sci USA. 1991, 88 (23): 10662-10666. 10.1073/pnas.88.23.10662.View ArticlePubMedPubMed CentralGoogle Scholar
- Bunn PA, Helfrich BA, Brenner DG, Chan DC, Dykes DJ, Cohen AJ, Miller YE: Effects of recombinant neutral endopeptidase (EC 126.96.36.199) on the growth of lung cancer cell lines in vitro and in vivo. Clin Cancer Res. 1998, 4 (11): 2849-2858.PubMedGoogle Scholar
- Wright JL, Jeng AY, Battistini B: Effect of ECE and NEP inhibition on cigarette smoke-induced cell proliferation in the rat lung. Inhal Toxicol. 2001, 13 (6): 497-511.View ArticlePubMedGoogle Scholar
- Karnoub AE, Dash AB, Vo AP, Sullivan A, Brooks MW, Bell GW, Richardson AL, Polyak K, Tubo R, Weinberg RA: Mesenchymal stem cells within tumour stroma promote breast cancer metastasis. Nature. 2007, 449 (7162): 557-563. 10.1038/nature06188.View ArticlePubMedGoogle Scholar
- Cui L, Ohuchida K, Mizumoto K, Moriyama T, Onimaru M, Nakata K, Nabae T, Ueki T, Sato N, Tominaga Y, Tanaka M: Prospectively isolated cancer-associated CD10(+) fibroblasts have stronger interactions with CD133(+) colon cancer cells than with CD133(-) cancer cells. PLoS One. 2010, 5 (8): e12121-10.1371/journal.pone.0012121.View ArticlePubMedPubMed CentralGoogle Scholar
- Soltermann A, Tischler V, Arbogast S, Braun J, Probst-Hensch N, Weder W, Moch H, Kristiansen G: Prognostic significance of epithelial-mesenchymal and mesenchymal-epithelial transition protein expression in non-small cell lung cancer. Clin Cancer Res. 2008, 14 (22): 7430-7437. 10.1158/1078-0432.CCR-08-0935.View ArticlePubMedGoogle Scholar
- Katsuda T, Tsuchiya R, Kosaka N, Yoshioka Y, Takagaki K, Oki K, Takeshita F, Sakai Y, Kuroda M, Ochiya T: Human adipose tissue-derived mesenchymal stem cells secrete functional neprilysin-bound exosomes. Sci Rep. 2013, 3: 1197-View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/40/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.