Open Access

Changes in circulating microRNAs after radiochemotherapy in head and neck cancer patients

  • Isolde Summerer1,
  • Maximilian Niyazi2,
  • Kristian Unger1, 4,
  • Adriana Pitea1,
  • Verena Zangen1, 4,
  • Julia Hess1, 4,
  • Michael J Atkinson3,
  • Claus Belka2, 4,
  • Simone Moertl3 and
  • Horst Zitzelsberger1, 4Email author
Radiation Oncology20138:296

https://doi.org/10.1186/1748-717X-8-296

Received: 13 November 2013

Accepted: 23 December 2013

Published: 28 December 2013

The Erratum to this article has been published in Radiation Oncology 2015 10:102

Abstract

Introduction

Circulating microRNAs (miRNAs) are easily accessible and have already proven to be useful as prognostic markers in cancer patients. However, their origin and function in the circulation is still under discussion. In the present study we analyzed changes in the miRNAs in blood plasma of head and neck squamous cell carcinoma (HNSCC) patients in response to radiochemotherapy and compared them to the changes in a cell culture model of primary HNSCC cells undergoing simulated anti-cancer therapy.

Materials and methods

MiRNA-profiles were analyzed by qRT-PCR arrays in paired blood plasma samples of HNSCC patients before therapy and after two days of treatment. Candidate miRNAs were validated by single qRT-PCR assays. An in vitro radiochemotherapy model using primary HNSCC cell cultures was established to test the possible tumor origin of the circulating miRNAs. Microarray analysis was performed on primary HNSCC cell cultures followed by validation of deregulated miRNAs via qRT-PCR.

Results

Unsupervised clustering of the expression profiles using the six most regulated miRNAs (miR-425-5p, miR-21-5p, miR-106b-5p, miR-590-5p, miR-574-3p, miR-885-3p) significantly (p = 0.012) separated plasma samples collected prior to treatment from plasma samples collected after two days of radiochemotherapy. MiRNA profiling of primary HNSCC cell cultures treated in vitro with radiochemotherapy revealed differentially expressed miRNAs that were also observed to be therapy-responsive in blood plasma of the patients (miR-425-5p, miR-21-5p, miR-106b-5p, miR-93-5p) and are therefore likely to stem from the tumor. Of these candidate marker miRNAs we were able to validate by qRT-PCR a deregulation of eight plasma miRNAs as well as miR-425-5p and miR-93-5p in primary HNSCC cultures after radiochemotherapy.

Conclusion

Changes in the abundance of circulating miRNAs during radiochemotherapy reflect the therapy response of primary HNSCC cells after an in vitro treatment. Therefore, the responsive miRNAs (miR-425-5p, miR-93-5p) may represent novel biomarkers for therapy monitoring. The prognostic value of this exciting observation requires confirmation using an independent patient cohort that includes clinical follow-up data.

Keywords

Head and neck cancer Circulating non-coding RNA Biomarker Radiotherapy outcome HNSCC cell culture model

Introduction

Surgical treatment of head and neck squamous cell carcinoma (HNSCC) is limited by the complex anatomy of the tumors and the associated risk of morbidity. Hence, for many patients an alternative treatment strategy of definitive radiotherapy alone or in conjunction with chemotherapy or immunotherapy is required[1]. Despite considerable progress in treatment options, disease recurrence and metastasis with a very strong impact on long-term survival are still the major challenges[2].

MicroRNAs (miRNAs) are evolutionarily conserved small RNAs, representing a class of regulators of posttranscriptional gene expression. There are more than 2000 mature miRNAs currently annotated in the human genome miRBase (release 19.0)[3] potentially targeting over 60% of all proteins[4]. Fundamental cellular processes including development, apoptosis[5], cell cycle control, proliferation[6] and DNA-damage repair[7] are influenced by miRNAs. During the last decade alterations in miRNA expression have been associated with a number of human diseases, including cancer (reviewed in[8]). The recent discovery of miRNAs in body fluids such as cerebrospinal fluid[9] and blood plasma[10] opens up the possibility of using miRNAs as minimally invasive biomarkers for the prediction of clinical endpoints such as overall survival[11, 12]. The well described ‘onco-miR’ miR-21 has already been identified as a valuable plasma biomarker with high prognostic power in esophageal squamous cell carcinoma[13] and gastric cancer[14].

Since the expression of miRNAs is known to be altered by ionizing radiation at both the cellular level[15] and plasma levels of mice[16], they might serve as easily accessible predictors of the individual response to radiation therapy. Additionally, miRNAs regulate drug sensitivity[17] and influence radioresistance[1820]. This offers the possibility of using specific plasma miRNAs as biomarkers for optimized treatment decisions.

The present study aimed to identify blood plasma miRNAs showing a response to radiochemotherapy and further to clarify the origin of these miRNAs in order to use them as minimally invasive tools for therapy monitoring.

For this purpose we compared miRNA levels in samples of blood plasma from HNSCC patients prior to treatment and after the completion of the first two fractions of therapeutic irradiation. MiRNAs displaying altered concentration levels after treatment were further analyzed in peripheral blood mononuclear cells (PBMC) of the same patients to test the hypothesis of PBMC-derived alterations of the plasma miRNAs. Moreover, we established an in vitro radiochemotherapy model using primary HNSCC cells in order to investigate a potential tumor origin of the therapy-responsive miRNAs in blood plasma.

Materials and methods

Patient samples

Plasma miRNA analysis was performed on 18 patients (17 HNSCC, 1 esophageal adenocarcinoma) treated with local X-ray-irradiation using a linear accelerator (6 MV, Siemens Mevatron M or ELEKTA Synergy®). After a planning (PET-) CT scan 70 Gy were applied to the macroscopic tumor and involved lymph nodes in daily dose fractions of 2 Gy five days per week. The adjuvant lymphatics were irradiated with up to 50 Gy and the high-risk lymphatics (adjacent to the involved lymph node levels) with up to 60 Gy.

16 out of 18 patients received concurrent chemotherapy (12 patients received 5-fluorouracil (5-FU) plus mitomycin C (MMC) [21], 3 patients MMC and 1 patient cisplatin weekly). None of the patients underwent surgical treatment. 5-FU treatment was usually applied on each of the first 5 days of therapeutic irradiation whereas MMC was applied only on day 5 and day 36 during radiotherapy. Patient characteristics are listed in Table  1.
Table 1

Patient characteristics

Characteristic

Number of patients

Gender

 
 

male

14

 

female

4

Median age, years

57.9

Age range, years

45.1–80.6

Tumor site

 
 

Larynx

5

 

Oropharynx

3

 

Mouth floor

2

 

Tongue

2

 

Esophagus

1

 

Hypopharynx

1

 

Maxilla

1

 

Nasopharynx

1

 

Sinuses

1

 

Soft palate

1

T-Stage

 
 

I

4

 

II

2

 

III

6

 

IV

6

N-Stage

 
 

N0

4

 

N1

4

 

N2

10

M-Stage

 
 

M0

16

 

M1

2

Concomitant therapy (in addition to radiotherapy)

 
 

  5-FU + MMC

12 (patient 3, 5–7, 9–13, 15–17)

 

  MMC

3 (patient 1, 14, 18)

 

  Cisplatin

1 (patient 8)

 

  Cetuximab

1 (patient 4)

 

  none

1 (patient 2)

Acute Toxicity

 
 

severe

11

 

moderate

5

 

n.a.

2

5-FU = 5-fluorouracil; MMC = mitomycin C;

n.a. = not available.

After obtaining ethical approval and informed consent, 15 ml of EDTA-peripheral blood were collected from each patient prior to the first fraction of therapy and within one hour after the second fraction of therapeutic irradiation. EDTA-blood samples were centrifuged at 350 × g for 10 min within two hours after collection to obtain plasma. To avoid cellular contamination the plasma samples were re-centrifuged at 1,200 × g for 3 min and subsequently at 14,000 × g for 10 min to remove cell debris.

PBMC were isolated using Ficoll gradient centrifugation.

Samples were stored at -20°C until further analysis or at -80°C for long term storage.

The study was approved by the ethics committee of the University of Munich (Germany).

Primary HNSCC cells

Primary tumor cells were obtained from the fresh tumor biopsies of two HNSCC cases, HN1957 (left maxilla/left nasal floor) and HN2092 (right floor of mouth) received from the Wales Cancer Bank, UK[22]. The tumor tissues were rinsed with PBS, minced and cultivated in keratinocyte media supplemented with L-glutamine, penicillin/streptomycin, EGF, BPE, CaCl2, F-12, DMEM (all components: Life Technologies) and 10% FBS at 37°C and 5% CO2. To ensure the exclusive cultivation of epithelial cells from tissue particles, fibroblast cells were identified via frequent microscopy and removed by scraping. The FBS-content of the media was continuously reduced down to serum-free cultivation. The epithelial origin of both cultures was assessed by positive cytokeratin staining via immunohistochemistry (Additional file1). Characteristics of the two primary HNSCC cultures are listed in Table 2.
Table 2

Characteristics of primary HNSCC cell cultures

Characteristic

HN1957

HN2092

Gender of patient

female

male

Age at diagnosis, years

85

73

Tumor site

left maxilla/left nasal floor

right floor of mouth

TNM

n.a.

pT4pN0

Cell type

epithelial

epithelial

n.a. = not available.

Treatment of HNSCC cells

To model as closely as possible the most frequent chemotherapy treatment we tested different 5-FU treatment options using the primary HNSCC cultures. Cell viability was measured using the XTT cell proliferation kit II (Roche) for both cell cultures after applying a wide range of 5-FU doses (Figure 1a, b). The 5-FU concentrations were selected to be relatively high, however lying within the linear region of the dose–response curve. Thus, HN1957 cells were treated with 50 μM 5-FU and HN2092 cells, as they showed less sensitivity towards the agent, were treated with 100 μM 5-FU.
Figure 1

Sensitivity of primary HNSCC cell cultures to 5-FU treatment and ionizing radiation. Cell viability (XTT-assay) of HN1957 (a) and HN2092 (b) 48 h after 5-FU treatment relative to DMSO-controls. Each data point represents the mean of three biological replicates. Dashed lines mark the 5-FU concentrations selected for in vitro radiochemotherapy modeling. Survival curves were generated by colony forming assay for HN1957 (c) 6 days after irradiation and HN2092 (d) 12 days after irradiation. Error bars represent the standard deviation of three biological replicates.

Cells were seeded in 60 mm-dishes and 24-well plates for RNA assays and cell viability assays, respectively. On the following day, cells were irradiated with 2 Gy using a 137Cs source and treated with 5-FU (Sigma; solved in DMSO). Controls were treated with the corresponding volumes of DMSO and sham-irradiated. 24 h after the first irradiation a second fraction of 2 Gy was applied to the 5-FU-treated cells followed by incubation for 1 h at 37°C. Cells were harvested by trypsinization and stored at -20°C until further processing.

Cell viability assay (XTT)

In order to determine cell viability of cells treated with in vitro radiochemotherapy or 5-FU treatment alone, 24-well plates were used for three biological and two technical replicates of each condition. XTT-assay was conducted using 350 μl of Medium and 150 μl of XTT-labeling mixture. Measurement was performed 24 h after 5-FU treatment.

Colony forming assay

In order to determine the radiation sensitivity of the primary HNSCC cells, the colony forming assay was performed in 6-well plates (Figure 1c, d). To ensure stable growth of the primary cells, a feeder layer consisting of the same HNSCC cells was used. Feeder cells were irradiated with 50 Gy in a 60Co source before seeding in the 6-well plates (10,000 cells/well). The following day, fresh cells for colony formation were plated on the feeder layer. On the following day the cells were irradiated using a 137Cs source with 0, 0.5, 1, 2, 3, 5 or 7.5 Gy. HN1957 colonies were stained with Giemsa six days and HN2092 colonies 12 days after irradiation. Colonies consisting of at least 50 cells were scored. Each colony formation assay was carried out in triplicate and repeated three times.

RNA extraction

Total RNA was extracted from 300 μl plasma using the mirVana miRNA Isolation Kit (Ambion) according to the manufacturer’s protocol with the following modification: 1,000 μl of lysis buffer was added to 300 μl of plasma. 1.04 fmol of synthetic cel-miR-39 (C. elegans) was added to each plasma sample (300 μl) after the protein denaturation step to normalize sample-to-sample variation in RNA recovery. RNA was finally eluted into 50 μl of pre-heated (95°C) nuclease-free water.

RNA from PBMC was extracted with the mirVana miRNA Isolation Kit (Ambion) according to the manufacturer’s protocol. RNA purity was measured by spectrophotometry (OD 260/280 ratio) using a Nanodrop ND-1000 (Thermo Scientific). Ratios were in the range of 1.89 to 2.09.

Extraction of total RNA from primary HNSCC cells was performed using the miRNeasy mini kit (Qiagen) according to the manufacturer’s protocol without DNase digest or small RNA enrichment. OD 260/280 ratios, measured on a Nanodrop ND-1000 (Thermo Scientific), were in the range of 1.92 to 2.04. Additionally, RNA quality was assessed prior to the Agilent microarray experiments using an Agilent 2100 Bioanalyzer (Agilent Technologies). The computed RNA integrity numbers (RINs) ranged from 9.3 to 10.0.

MiRNA profiling in patient plasma

TaqMan Array Human MicroRNA A Cards v2.0 (Applied Biosystems) representing 377 mature human miRNAs were used to obtain miRNA profiles of the blood plasma samples. Reverse transcription was performed using the TaqMan miRNA reverse transcription kit (Applied Biosystems) in combination with the stem-loop megaplex primer pool set A v2.1. Because of the low abundance of circulating miRNA in the starting material we used a fourfold volume of the reaction mix, prepared according to the manufacturer’s protocol adding 12 μl of total plasma RNA.

For the subsequent quantitative real-time PCR 30 μl of cDNA was used for each array and PCR was carried out on an Applied Biosystems 7900HT. The reaction mixtures were incubated at 50°C for 2 min and 95°C for 10 min, followed by 45 cycles of 95°C for 30 s and 60°C for 1 min. All Ct values were normalized using the median Ct value of all detectable miRNAs on the array.

MicroRNA profiling in HNSCC cells

To analyze the effects of in vitro radiochemotherapy on the cellular miRNA expression levels Sure Print G3 human 8x60k miRNA microarrays (Agilent Technologies) were used covering 1205 human miRNAs (Sanger miRBase release 16). 100 ng of total RNA was dephosphorylated and labeled with cyanine 3-cytidine biphosphate including a labeling spike-in solution (Agilent Technologies) to assess the labeling efficiency. After purification of the labeled RNA, the samples were hybridized on the arrays including a hybridization spike-in solution (Agilent Technologies) to monitor hybridization efficiency. Arrays were scanned with a G2505C Sure Scan Microarray Scanner (Agilent Technologies) using Scan Control software. For data processing the following software was used: Feature extraction 10.7 (Agilent Technologies) and GeneSpring 11.5. MiRNA analysis was conducted with three biological and two technical replicates for each data point. All steps were performed according to the manufacturer’s protocol. The raw data were imported into the R statistical platform[23]. The total microRNA gene signal (TGS) was computed by the Agilent Extraction feature. The resulting signal was normalized to the mean signal of all arrays for compensation of systematic technical slide-to-slide variations. A differential expression analysis was performed using the linear model components implemented in limma[24] and empirical Bayes methods were applied in order to attain moderated statistics[25].

Real-time PCR quantification of individual miRNAs

Reverse transcription was performed using the TaqMan miRNA reverse transcription kit and miRNA-specific stem-loop primers (Applied Biosystems). Reverse transcription was performed on a Cyclone PCR system (Peqlab) according to the manufacturer’s protocol. For miRNA assays in plasma samples 3 μl of cDNA was used for reverse transcription because of the low RNA content. Real-time PCR was performed in duplicates and included non-template negative controls. PCR was performed on a ViiA 7 real-time PCR System (Applied Biosystems) following the manufacturer’s protocol. For plasma samples the miRNA concentration levels were normalized to the spiked-in cel-miR-39, for cellular samples the U6 snRNA was used for normalization. Fold changes were calculated using the 2-ΔΔCt method[26].

Statistical analysis

To identify differentially expressed miRNAs in plasma samples of patients prior to radiotherapy and after the second fraction of therapeutic irradiation Wilcoxon test was performed. All miRNAs that were detected in less than 30% of either control samples or irradiated samples were excluded from further analysis. Unsupervised hierarchical clustering was performed using the top distinctive miRNAs applying the parameters maximum distance and Ward’s method. Fisher’s exact test was applied to check for significant clustering of samples from patients prior to treatment and those collected after treatment. P values < 0.05 were considered statistically significant.

Correlation coefficients for miRNA changes in PBMC and plasma as well as in arrays and single assays were calculated using Pearson correlation. All statistical analyses were performed using The R Project for Statistical Computing[23].

Results

MiRNA profiling of the plasma samples from 18 head and neck cancer patients revealed 54 miRNAs with altered expression levels after treatment (Additional file2). Unsupervised hierarchical clustering using the expression levels of the top six significantly deregulated miRNAs revealed two main clusters. Cluster 1 represents ten samples collected prior to treatment and two samples collected after treatment, whereas cluster 2 consists of eight pre-treatment samples and 16 post-treatment samples. Consequently, the miRNA profiles of these six miRNAs differentiated significantly (p = 0.012) between samples collected before therapy and those collected after two days of treatment (Figure 2). For technical validation of the array data eight miRNAs (miR-590-5p, miR-574-3p, miR-425-5p, miR-885-3p, miR-21-5p, miR-28-3p, miR-195-5p, miR-191-5p) were selected from the top list of deregulated miRNAs. Additionally, miR-150-5p and miR-142-3p, both known to play a role in esophageal carcinoma[27, 28], were used for validation. For these ten miRNAs the differences between the two groups were validated using single qRT-PCR assays. The correlation coefficients of the normalized Ct values demonstrated a good correlation between the results of the TaqMan low density array and the TaqMan single assays for most of the tested miRNAs (Additional file3). The single assays of miR-590-5p and miR-885-3p were below the detection level in almost all plasma samples and were consequently excluded from further analyses.
Figure 2

Unsupervised hierarchical cluster analysis of samples collected before therapy and after two days of treatment. Samples collected prior to treatment are marked with blue bars, samples collected after two days of treatment (XRT = radiotherapy) are marked with red bars. Cluster 1 and Cluster 2 indicate the two main clusters emerging from unsupervised cluster analysis using blood plasma concentration levels of the top six distinctive miRNAs.

Many plasma miRNAs are known to originate from peripheral blood mononuclear cells (PBMC)[29]. Therefore, we measured the expression levels of the candidate therapy biomarker miRNAs (miR-574-3p, miR-425-5p, miR-21-5p, miR-28-3p, miR-195-5p, miR-191-5p, miR-150-5p, miR-142-3p) in PBMC samples obtained from the same patients. Correlation coefficients of the normalized Ct values indicated that the changed miRNA levels after treatment do not originate from PBMC (Additional file4).

In order to further investigate whether the plasma miRNA changes following tumor therapy originate from the tumor cells we established an in vitro radiochemotherapy model using primary HNSCC from two different patients, HN1957 and HN2092 (Table 2).

For assessment of the sensitivity of the primary tumor cells to 5-FU and to ionizing radiation dose-response curves (Figure 1a, b) and survival curves (Figure 1c, d) were generated, respectively. HN2092 cells showed higher sensitivity to ionizing radiation but lower sensitivity to 5-FU treatment compared to HN1957 cells. In order to simulate the treatment of the patient samples in the present study cells were irradiated with 2 Gy followed by treatment with 5-FU and on the following day irradiation with a second 2 Gy fraction. Subsequent measurement of the cell viability (24 h after 5-FU treatment) showed no significant effect of 5-FU or irradiation on HN2092 cell cultures, whereas cell viability of HN1957 cell cultures was already reduced by 5-FU treatment alone with a significant additive effect of irradiation on the cell viability (Figure 3).
Figure 3

Cell viability (XTT-assay) of primary HNSCC cell cultures. Absorption was measured 24 h after exclusive 5-FU treatment and 5-FU treatment in combination with 2 x 2 Gy in vitro irradiation relative to DMSO-controls. Error bars represent the standard deviation of three biological replicates (*p < 0.05, **p < 0.01).

MiRNA profiling using Agilent microarrays in the sham-irradiated and DMSO-treated control cells compared to cells treated with in vitro radiochemotherapy revealed several miRNAs that differentiate control cells and treated cells for HN1957 and HN2092 (Additional files5 and6). We selected three miRNAs (miR-425-5p, miR-21-5p, miR-106b-5p) that were altered after in vitro treatment and were also in the top list of therapy-responsive miRNAs in patient plasma samples after in vivo tumor therapy. Additionally, miR-93-5p was included as it shows deregulation after in vitro treatment and in vivo treatment in most of the patients, although it is not among the top significant therapy-responsive miRNAs due to the different directions of regulation (Additional file7). For these miRNAs we conducted TaqMan single assays in control and treated cells to validate the array results. The qRT-PCR assays confirmed an up-regulation of miR-425-5p with a p value close to the significance level (p = 0.552) in HN1957 cell cultures, while for miR-93-5p a significant down-regulation (p < 0.001) in HN2092 cell cultures became apparent (Table 3). Thus, at least two of the four treatment-regulated plasma miRNAs are likely to be associated with the direct response of tumor cells to treatment.
Table 3

MicroRNAs significantly deregulated after radiochemotherapeutic treatment in primary HNSCC cell cultures and in plasma of HNSCC patients

 

Agilent microarray

TaqMan single qRT-PCR assay

miRNA

HN1957

HN2092

HN1957

HN2092

FC (p value)

FC (p value)

FC (p value)

FC (p value)

miR-425-5p

1.27 (0.004)

-

1.28 (0.052)

-

miR-21-5p

1.11 (0.031)

-

1.11 (0.543)

-

miR-106b-5p

0.95 (0.042)

-

1.04 (0.830)

-

miR-93-5p

0.92 (0.023)

0.96 (0.001)

1.03 (0.552)

0.84 (0.000)

FC = fold change.

Discussion

The present study aimed to identify circulating miRNA markers in the plasma of HNSCC patients that are indicative of the efficacy of radiochemotherapy. Our finding that a limited panel of deregulated plasma miRNAs discriminates between pre- and post-radiochemotherapy samples (Figure 2) supports previous evidence[16] for the potential use of plasma miRNAs as radiation-responsive biomarkers. Moreover, several candidates altered after the second fraction of therapeutic irradiation, such as miR-93-5p, miR-142-3p, miR-106b-5p, miR-191-5p and miR-21-5p (Additional file2) have been previously described to be induced by ionizing radiation[30].

Although it is known that circulating miRNAs are present in stable forms, there is still uncertainty about their transport and cellular origin (reviewed in[31]). Since circulating miRNAs may originate from blood cells[29] we tested the hypothesis that the observed therapy-responsive plasma miRNAs originate from peripheral blood mononuclear cells (PBMC) of the same patients. Our analyses of the miRNA profiles of PBMC revealed no significant correlation between the miRNA expression in PBMC and plasma in patients before or after treatment (Additional file4). We conclude that the plasma miRNA changes detected in response to anti-cancer therapy are most likely not originating from circulating PBMC.

It has frequently been hypothesized that plasma miRNAs originate from tumor cells, either based on active secretion[32] or in apoptotic bodies released from dead tumor cells[33]. This is supported by the present study as several candidates among the top significant plasma miRNAs that were altered after radiochemotherapy (Additional file2) have already been reported to be tumor markers before. MiR-195-5p, miR-574-3p and miR-28-3p have been described as being differentially expressed in esophageal squamous cell carcinoma[34, 35]. MiR-191-5p and miR-21-5p play a role in lung cancer diagnosis and prognosis[36], while miR-21-5p in plasma also is known to serve as prognostic marker in esophageal cancer patients[13]. Taken together, these reports indicate that the observed changes in plasma miRNAs in our study might be related to therapy effects on tumor cells.

To further strengthen this hypothesis we investigated the miRNA expression of primary HNSCC cell cultures under simulated radiochemotherapeutic treatment in vitro. The results of the miRNA expression analysis of this in vitro model further supported the assumption that the observed plasma miRNA changes are likely to originate from the tumor cells. We were able to show a significant deregulation of four plasma miRNAs in HN1957 cell cultures (miR-425-5p, miR-21-5p, miR-106b-5p, miR-93-5p) and one plasma miRNA in HN2092 cell cultures (miR-93-5p) by microarray profiling. The observation that the same miRNAs were altered in both tumor cells and plasma leads us to conclude that plasma miRNA changes following radiochemotherapeutic treatment are the result of miRNA release from damaged tumor cells. The observed differences in deregulated miRNAs between the two cell cultures, however, were not comparable. The two primary HNSCC cell cultures had different responses to radiochemotherapeutic treatment. This may be due to differences in the responsiveness of the cell cultures. Indeed, the cell viability assay (XTT) exhibited a significant effect 24 h after irradiation and 5-FU treatment on HN1957 cells, but not on HN2092 cells (Figure 3). This reflects a higher short-term toxicity of the treatment on HN1957 cell cultures compared to HN2092 cell cultures, whereas the long-term effect of irradiation is stronger on HN2092 cell cultures as shown in the colony forming assay (Figure 1c, d). The common deregulation of 43 miRNAs in both primary HNSCC cell cultures (Additional files5 and6) suggests a similar miRNA response to therapeutic treatment independent from the individual tumor. We also observed differences in the miRNA response of both HNSCC cell cultures, which might reflect variations in the individual sensitivity to radiation and chemotherapeutic treatment (Figure 1).

The major finding of this study, i.e. therapy-related miRNA expression in the blood plasma is similar to miRNA changes reported for treated tumor cells, is not only based on global miRNA array screens but was also technically validated by qRT-PCR studies of the most distinctive miRNAs (Table 3). We identified miR-425-5p (HN1957) and miR-93-5p (HN2092) as top candidates of commonly deregulated miRNAs in primary HNSCC cell cultures and blood plasma of HNSCC patients.

In conclusion, the present study compared miRNA data from blood plasma of radiochemotherapy-treated HNSCC patients and from an in vitro model of primary tumor cell radiochemotherapy. This comparison sheds light on the origin and potential use of circulating miRNAs in the blood plasma of HNSCC patients as therapy-responsive biomarkers. We will further integrate these therapy-responsive plasma miRNAs with clinical data of the patients (e.g. overall survival) and seek for independent validation in another patient cohort in order to evaluate their potential for prognostic use. The established HNSCC cell culture model further allows us to investigate the molecular interactome of the deregulated miRNA candidates by adding mRNA expression data and subsequent integrative data analysis.

The current study strongly suggests that alterations of miRNAs following radiochemotherapy in the blood plasma are associated with the tumor response to therapy and therefore might represent novel biomarkers for therapy monitoring.

Notes

Declarations

Acknowledgements

We thank Klaudia Winkler, Aaron Selmeier and Claire Innerlohinger for their excellent technical assistance. We further thank Herbert Braselmann for the graphical design of survival curves as well as Maya Flieger for the collection of clinical data of the patients and the support for patient recruitment. Tissue samples were obtained from the Wales Cancer Bank, which is funded by the Wales Assembly Government and Cancer Research Wales. Other investigators may have received specimens from the same subjects.

Authors’ Affiliations

(1)
Research Unit Radiation Cytogenetics, Helmholtz Center Munich
(2)
Department of Radiation Oncology, University of Munich
(3)
Institute of Radiation Biology, Helmholtz Center Munich
(4)
Clinical Cooperation Group ‘Personalized Radiotherapy of Head and Neck Cancer’, Helmholtz Center Munich

References

  1. May JT, Rao N, Sabater RD, Boutrid H, Caudell JJ, Merchant F, Han G, Padhya TA, Mccaffrey JC, Tanvetyanon T, Deconti R, Kish J, Mccaffrey TV, Trotti A: Intensity-modulated radiation therapy as primary treatment for oropharyngeal squamous cell carcinoma. Head & neck 2013, 35: 1796-1800. 10.1002/hed.23245View ArticleGoogle Scholar
  2. Ho AS, Kraus DH, Ganly I, Lee NY, Shah JP, Morris LG: Decision making in the management of recurrent head and neck cancer. Head & neck 2013, 36: 144-151.View ArticleGoogle Scholar
  3. miRBase: the microRNA database. http://www.miRbase.org
  4. Friedman RC, Farh KK-H, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 2009, 19: 92-105.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Liu X, Jiang L, Wang A, Yu J, Shi F, Zhou X: MicroRNA-138 suppresses invasion and promotes apoptosis in head and neck squamous cell carcinoma cell lines. Cancer Lett 2009, 286: 217-222. 10.1016/j.canlet.2009.05.030View ArticlePubMedPubMed CentralGoogle Scholar
  6. Novello C, Pazzaglia L, Cingolani C, Conti A, Quattrini I, Manara MC, Tognon M, Picci P, Benassi MS: miRNA expression profile in human osteosarcoma: role of miR-1 and miR-133b in proliferation and cell cycle control. Int j of oncol 2013, 42: 667-675.Google Scholar
  7. Wu CW, Dong YJ, Liang QY, He XQ, Ng SS, Chan FK, Sung JJ, Yu J: MicroRNA-18a attenuates DNA damage repair through suppressing the expression of ataxia telangiectasia mutated in Colorectal cancer. PloS one 2013, 8: e57036. 10.1371/journal.pone.0057036View ArticlePubMedPubMed CentralGoogle Scholar
  8. Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer 2006, 6: 857-866. 10.1038/nrc1997View ArticlePubMedGoogle Scholar
  9. Gallego J, Gordon M, Claycomb K, Bhatt M, Lencz T, Malhotra A: In vivo microRNA detection and quantitation in cerebrospinal fluid. J Mol Neurosci 2012, 47: 243-248. 10.1007/s12031-012-9731-7View ArticlePubMedPubMed CentralGoogle Scholar
  10. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, Peterson A, Noteboom J, O’briant KC, Allen A, Lin DW, Urban N, Drescher CW, Knudsen BS, Stirewalt DL, Gentleman R, Vessella RL, Nelson PS, Martin DB, Tewari M: Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci 2008, 105: 10513-10518. 10.1073/pnas.0804549105View ArticlePubMedPubMed CentralGoogle Scholar
  11. Liu R, Chen X, Du Y, Yao W, Shen L, Wang C, Hu Z, Zhuang R, Ning G, Zhang C, Yuan Y, Li Z, Zen K, Ba Y, Zhang C-Y: Serum microRNA expression profile as a biomarker in the diagnosis and prognosis of pancreatic cancer. Clin Chem 2012, 58: 610-618. 10.1373/clinchem.2011.172767View ArticlePubMedGoogle Scholar
  12. Niyazi M, Zehentmayr F, Niemoller O, Eigenbrod S, Kretzschmar H, Schulze-Osthoff K, Tonn J-C, Atkinson M, Mortl S, Belka C: MiRNA expression patterns predict survival in glioblastoma. Radiat oncol 2011, 6: 153. 10.1186/1748-717X-6-153View ArticlePubMedPubMed CentralGoogle Scholar
  13. Komatsu S, Ichikawa D, Takeshita H, Konishi H, Nagata H, Hirajima S, Kawaguchi T, Arita T, Shiozaki A, Fujiwara H, Okamoto K, Otsuji E: Prognostic impact of circulating miR-21 and miR-375 in plasma of patients with esophageal squamous cell carcinoma. Expert Opin on Biol Ther 2012,12(1):S53-59. 10.1517/14712598.2012.636349View ArticleGoogle Scholar
  14. Komatsu S, Ichikawa D, Tsujiura M, Konishi H, Takeshita H, Nagata H, Kawaguchi T, Hirajima S, Arita T, Shiozaki A, Kubota T, Fujiwara H, Okamoto K, Otsuji E: Prognostic impact of circulating miR-21 in the plasma of patients with gastric carcinoma. Anticancer Res 2013, 33: 271-276.PubMedGoogle Scholar
  15. Simone NL, Soule BP, Ly D, Saleh AD, Savage JE, Degraff W, Cook J, Harris CC, Gius D, Mitchell JB: Ionizing radiation-induced oxidative stress alters miRNA expression. PloS one 2009, 4: e6377. 10.1371/journal.pone.0006377View ArticlePubMedPubMed CentralGoogle Scholar
  16. Cui W, Ma J, Wang Y, Biswal S: Plasma miRNA as biomarkers for assessment of total-body radiation exposure dosimetry. PloS one 2011, 6: e22988. 10.1371/journal.pone.0022988View ArticlePubMedPubMed CentralGoogle Scholar
  17. Du L, Pertsemlidis A: microRNA regulation of cell viability and drug sensitivity in lung cancer. Expert Opin on biol Ther 2012, 12: 1221-1239. 10.1517/14712598.2012.697149View ArticleGoogle Scholar
  18. Xiao-Chun W, Wei W, Zhu-Bo Z, Jing Z, Xiao-Gang T, Jian-Chao L: Overexpression of miRNA-21 promotes radiation-resistance of non-small cell lung cancer. Radiat Oncol 2013, 8: 146. 10.1186/1748-717X-8-146View ArticlePubMed CentralGoogle Scholar
  19. Anastasov N, Hofig I, Vasconcellos IG, Rappl K, Braselmann H, Ludyga N, Auer G, Aubele M, Atkinson M: Radiation resistance due to high expression of miR-21 and G2/M checkpoint arrest in breast cancer cells. Radiat Oncol 2012, 7: 206. 10.1186/1748-717X-7-206View ArticlePubMedPubMed CentralGoogle Scholar
  20. Svoboda M, Sana J, Fabian P, Kocakova I, Gombosova J, Nekvindova J, Radova L, Vyzula R, Slaby O: MicroRNA expression profile associated with response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer patients. Radiat Oncol 2012, 7: 195. 10.1186/1748-717X-7-195View ArticlePubMedPubMed CentralGoogle Scholar
  21. Budach V, Stuschke M, Budach W, Baumann M, Geismar D, Grabenbauer G, Lammert I, Jahnke K, Stueben G, Herrmann T, Bamberg M, Wust P, Hinkelbein W, Wernecke KD: Hyperfractionated accelerated chemoradiation with concurrent fluorouracil-mitomycin is more effective than dose-escalated hyperfractionated accelerated radiation therapy alone in locally advanced head and neck cancer: final results of the radiotherapy cooperative clinical trials group of the German cancer society 95–06 prospective randomized trial. J of Clin Oncol: official j of the AmSocof ClinOncol 2005, 23: 1125-1135. 10.1200/JCO.2005.07.010View ArticleGoogle Scholar
  22. Wales cancer bank. http://www.walescancerbank.com
  23. The R project for statistical computing. http://www.r-project.org
  24. Smyth GK: Limma: linear models for microarray data. In Bioinformatics and computational biology solutions using R and bioconductor. Edited by: Gentleman R, Carey V, Huber W, Irizarry R, Dudoit S. New York: Springer; 2005:397-420. Statistics for Biology and HealthView ArticleGoogle Scholar
  25. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl in Genet and Mol Biol 2004, 3: Article3.Google Scholar
  26. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods (San Diego, Calif) 2001, 25: 402-408. 10.1006/meth.2001.1262View ArticleGoogle Scholar
  27. Yokobori T, Suzuki S, Tanaka N, Inose T, Sohda M, Sano A, Sakai M, Nakajima M, Miyazaki T, Kato H, Kuwano H: MiR-150 is associated with poor prognosis in esophageal squamous cell carcinoma via targeting the EMT inducer ZEB1. Cancer Sci 2013, 104: 48-54. 10.1111/cas.12030View ArticlePubMedGoogle Scholar
  28. Lin R-J, Xiao D-W, Liao L-D, Chen T, Xie Z-F, Huang W-Z, Wang W-S, Jiang T-F, Wu B-L, Li E-M, Xu L-Y: MiR-142-3p as a potential prognostic biomarker for esophageal squamous cell carcinoma. J Surg Oncol 2012, 105: 175-182. 10.1002/jso.22066View ArticlePubMedGoogle Scholar
  29. Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, Xiao T, Schafer J, Lee M-LT, Schmittgen TD, Nana-Sinkam SP, Jarjoura D, Marsh CB: Detection of microRNA expression in human peripheral blood microvesicles. PloS one 2008, 3: e3694. 10.1371/journal.pone.0003694View ArticlePubMedPubMed CentralGoogle Scholar
  30. Chaudhry MA, Omaruddin RA, Brumbaugh CD, Tariq MA, Pourmand N: Identification of radiation-induced microRNA transcriptome by next-generation massively parallel sequencing. J Radiat Res 2013, 54: 808-822. 10.1093/jrr/rrt014View ArticlePubMedPubMed CentralGoogle Scholar
  31. Selth LA, Tilley WD, Butler LM: Circulating microRNAs: macro-utility as markers of prostate cancer? Endocr Relat Cancer 2012, 19: R99-R113. 10.1530/ERC-12-0010View ArticlePubMedGoogle Scholar
  32. Hannafon B, Ding W-Q: Intercellular communication by exosome-derived microRNAs in cancer. Int J Mol Sci 2013, 14: 14240-14269. 10.3390/ijms140714240View ArticlePubMedPubMed CentralGoogle Scholar
  33. Zernecke A, Bidzhekov K, Noels H, Shagdarsuren E, Gan L, Denecke B, Hristov M, Koppel T, Jahantigh MN, Lutgens E, Wang S, Olson EN, Schober A, Weber C: Delivery of microRNA-126 by apoptotic bodies induces CXCL12-dependent vascular protection. Sci Signal 2009, 2: ra81.View ArticlePubMedGoogle Scholar
  34. Liu SG, Qin XG, Zhao BS, Qi B, Yao WJ, Wang TY, Li HC, Wu XN: Differential expression of miRNAs in esophageal cancer tissue. Oncology letters 2013, 5: 1639-1642.PubMedPubMed CentralGoogle Scholar
  35. Fu M-G, Li S, Yu T-T, Qian L-J, Cao R-S, Zhu H, Xiao B, Jiao C-H, Tang N-N, Ma J-J, Hua J, Zhang W-F, Zhang H-J, Shi R-H: Differential expression of miR-195 in esophageal squamous cell carcinoma and miR-195 expression inhibits tumor cell proliferation and invasion by targeting of Cdc42. FEBS Lett 2013, 587: 3471-3479. 10.1016/j.febslet.2013.08.036View ArticlePubMedGoogle Scholar
  36. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, Calin GA, Liu C-G, Croce CM, Harris CC: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell 2006, 9: 189-198. 10.1016/j.ccr.2006.01.025View ArticlePubMedGoogle Scholar

Copyright

© Summerer et al.; licensee BioMed Central Ltd. 2013

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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.