Gene expression profiling of human dermal fibroblasts exposed to bleomycin sulphate does not differentiate between radiation sensitive and control patients
Radiation Oncology volume 6, Article number: 42 (2011)
Gene expression profiling of the transcriptional response of human dermal fibroblasts to in vitro radiation has shown promise as a predictive test of radiosensitivity. This study tested if treatment with the radiomimetic drug bleomycin sulphate could be used to differentiate radiation sensitive patients and controls in patients who had previously received radiotherapy for early breast cancer.
Eight patients who developed marked late radiation change assessed by photographic breast appearance and 8 matched patients without any change were selected from women entered in a prospective randomised trial of breast radiotherapy fractionation. Gene expression profiling of primary skin fibroblasts exposed in vitro to bleomycin sulphate and mock treated fibroblast controls was performed. 973 genes were up-regulated and 923 down-reguated in bleomycin sulphate treated compared to mock treated control fibroblasts. Gene ontology analysis revealed enriched groups were cellular localisation, apoptosis, cell cycle and DNA damage response for the deregulated genes. No transcriptional differences were identified between fibroblasts from radiation sensitive cases and control patients; subgroup analysis using cases exhibiting severe radiation sensitivity or with high risk alleles present in TGF β1 also showed no difference.
The transcriptional response of human dermal fibroblasts to bleomycin sulphate has been characterised. No differences between clinically radiation sensitive and control patients were detected using this approach.
Gene expression profiling of in vitro cellular responses of human fibroblasts and lymphocytes to radiation has demonstrated that cells undergo complex early transcriptional responses of a wide spectrum of genes from different gene ontologies [1–4]. Microarray studies have demonstrated that the transcriptional response of human cells exposed to radiation in vitro differs between radiation sensitive patients and controls. Therefore this approach has been explored as a predictive test of radiation sensitivity using late normal tissue effects as the endpoint of radiation sensitivity [5–7].
The spectrum of DNA damage caused by bleomycin sulphate is similar but not identical to that caused by ionising radiation, hence its definition as a radiomimetic agent . The molecular and clinical responses after bleomycin sulphate and radiation are similar: both induce post-mitotic differentiation of fibroblasts inducing a senescent phenotype associated with increased collagen production [9–11], activate cascades of profibrotic chemokines and cytokines and cause skin and pulmonary fibrosis in animal models and in the clinic [12–14]. On this basis, the potential of using bleomycin sulphate rather than radiation for predictive testing is here tested in an exploratory study.
Materials and methods
Patients and assessment of late normal tissue injury
Patients with a history of early breast cancer treated with breast conserving surgery and radiotherapy within a clinical trial of radiotherapy fractionation were included. This patient group had prospective scoring of late normal tissue effects . Using photographic scores, cases were identified as marked change in appearance (grade 3) at any assessment or a persistent moderate change (grade 2) for at least 3 consecutive years. Controls had no or minimal change to breast appearance (grade 1) and were matched to cases using defined clinical parameters [see Additional File 1].
One hundred cases with radiation change and 200 matched controls with no change were selected for translational research studies. Primary dermal fibroblasts from 26/100 best matched case control pairs were prepared from explant biopsies of buttock skin as previously described . A subset of 8 case control pairs was selected for this study. Ethical approval was given by the Royal Marsden NHS Foundation Trust Ethics Committee and all patients gave written informed consent.
Cell culture and treatment with bleomycin sulphate
Fibroblasts were seeded at passage 7-9 into T75 cm2 flasks on day 0 and cultured in DMEM/10% FCS. The dose of bleomycin sulphate was previously determined by cell cycle analysis using fluorescence activated cell sorting (FACS) after treatment of fibroblasts with different doses [see Additional File 2]. On day 1, during exponential growth phase, cells were treated with 10 μg/ml bleomycin sulphate or medium alone (mock control) for 6 hours. Cells were then washed in PBS and then cultured in DMEM/10% FCS prior to RNA extraction on day 4.
Cells were washed once in PBS and lysed using the RNeasy Mini Kit (Qiagen). Briefly, cells were lysed in 600 μl extraction buffer and samples were homogenised by centrifugation in a QIAshredder (Qiagen) and stored at -80°C. For RNA isolation, 70% ethanol was added to the sample to bind the RNA to a silica membrane filter, impurities were removed by washing before finally eluting RNA in 30 μl water. RNA was concentrated in a vacuum centrifuge and analysed using an Agilent 2100 Bioanalyser (Agilent Technologies).
RNA processing and hybridisation to Affymetrix chip
Biotinylated target RNA was prepared with minor modifications from the manufacturer's recommendations, on Affymetrix http://www.affymetrix.com/support/technical/manual/expression_manual.affx. Target RNA generated from each sample was then processed using an Affymetrix GeneChip Instrument System and hybridised to HGU133plus2 arrays. Arrays were scanned in an Affymetrix GeneChip Scanner 3000 system. Data pre-processing was carried out using Affymetrix GeneChip Operating software. The data is available on MIAMEVICE http://bioinformatics.picr.man.ac.uk/vice/PublicProjects.vice.
Statistical analysis of Affymetrix data
Subsequent data processing was done using GCRMA package in The R project (R 2.6) http://www.r-project.org/. Data was RMA normalised and mapped to Ensembl gene IDs by Brainarray Custom CDF http://brainarray.mbni.med.umich.edu/Brainarray/Database/CustomCDF/genomic_curated_CDF.asp. Absolute values less than 50 were floored to 50. Log2 ratios were calculated and data was median centred. Genes for which expression values were present in less than 80% of samples were excluded. Significance Analysis of Microarrays (SAM)  was performed to identify genes differentially expressed between sample groups.
Gene Ontology Analysis
To identify enriched gene ontology groups, the functional annotation tool on DAVID Bioinformatics Resources 2007 was used http://david.abcc.ncifcrf.gov/home.jsp. Enriched gene ontologies ranked highly according to statistical significance were identified (EASE score, modified Fisher exact p-value).
Patient characteristics are shown in Table 1. Of the 8 radiation sensitive cases included, 3 cases had grade 3 scores of change in photographic breast appearance at 5 years i.e. marked radiation change. The remainder were scored as cases with moderate change (grade 2) for at least 3 successive years. All control patients had grade 1 scores i.e. no/minimal change in breast appearance.
Transcriptional response of cultured fibroblasts exposed to bleomycin sulphate
To identify the transcriptional response of cells exposed to bleomycin sulphate, paired SAM of bleomycin sulphate and mock treated samples was carried out using all 16 fibroblast cultures. Of the statistically significant differentially expressed genes (false discovery rate = 0), 973 genes were up-regulated and 923 genes were down-regulated in bleomycin sulphate treated compared to mock treated fibroblasts [see Additional File 3].
For fibroblast response to bleomycin sulphate, the highly enriched gene ontologies for up-regulated genes were cellular localization and cell death (Table 2) and for down-regulated genes included regulation of progression through cell cycle, mitotic phase of cell cycle and DNA damage response and repair (Table 3).
Differences between radiation sensitive patients and matched controls
SAM was used to try to identify transcriptional differences between fibroblasts isolated from radiation sensitive patients and controls. One case control pair was excluded from analysis as it was incorrectly matched for radiotherapy fractionation.
Comparisons between radiation sensitive cases and matched controls were made for mock treated fibroblast samples (i.e. not exposed to bleomycin sulphate) and for bleomycin sulphate treated fibroblast samples. No significant differentially expressed genes were identified in either comparison (data not shown).
The next approach taken was to calculate fold induction in transcript levels using gene expression ratios of log2 values for bleomycin sulphate treated compared to mock treated samples. Again, no significant differentially expressed genes were identified (data not shown).
Andreassen et al. reported statistically significant associations for 2 single nucleotide polymorphisms (SNPs) in TGF β1 (positions -509 and codon 10) and the risk of developing late normal tissue effects in the same patient population . Of the 7 case control pairs in the current study, both high risk alleles were present in 5 cases and in 3 controls. The 5 cases with both high risk alleles and 5 matched controls were selected for further analyses. Two out of 5 of these matched controls had both high risk alleles present. SAM analysis was performed on the 5 case control pairs using fold induction values. No significant differentially expressed genes were identified (data not shown). Further analysis using 3 selected cases with marked radiation change and all 7 control patients also did not identify significant differentially expressed genes (data not shown).
In the current study, transcriptional profiling of dermal fibroblasts after exposure to bleomycin sulphate was carried out to determine whether differences in transcriptional response could be identified between patients with late normal tissue radiation effects and matched controls. This was a pilot study to determine if bleomycin sulphate could be used as an alternative to radiation, in this context. No differences were detected between the 2 patient groups. There are a number of possible explanations for this negative finding.
In this study, a score for late normal tissue effects was performed using photographic appearance. Cases had moderate/marked change in breast appearance and matched controls had no/minimal change. The case control selection may be a limitation of the current study. Although confounding factors such as breast size were taken into account, conventional planning techniques were used in this population leading to variations in dosimetry between cases and controls. The limitation of 5 years of follow up may have incorrectly classified radiation sensitive cases into the control group in those patients whose late normal tissue injury became manifest later. However, time to development of late normal tissue injury is a relevant parameter for judging radiation sensitivity. Other than patient 112 for whom data was missing after year 3, the available data showed no apparent injury for 4 control patients at 5 years and for 3 control patients at 10 years after radiation (Table 1).
One of the main limitations of the study was the sample size. This and the related issue of inter-sample variation may have contributed to the negative findings. An additional possible source of variation was that cells were not synchronised prior to treatment. The issue of inter-sample variation was further addressed in the microarray analysis by using fold induction values between drug and mock treated samples. The potential superiority of this approach is supported in another study of predictive testing of radiation sensitivity from the Danish cohort of breast cancer patients . In this study, when cDNA array analysis of basal gene expression was compared between two patient groups, defined as radiation 'sensitive' and radiation 'resistant', only 6 genes were identified as being differentially expressed, suggesting that the difference between untreated fibroblasts from the two groups is likely to be small . The authors selected 17 differentially expressed candidate genes between the two groups, identified in irradiated fibroblast samples, which were further analysed by quantitative real time polymerase chain reaction (Q-RT-PCR) . The study reported that using fold induction values better differentiated radiation 'sensitive' and radiation 'resistant' patients than using either untreated samples or radiation exposed samples alone. Fold induction takes into account background levels (i.e. the transcriptional profile of untreated samples) and thereby controls for genetic variation. However in the current study, a difference between radiation sensitive cases and controls was not detected even with this approach.
Bleomycin sulphate stimulates post-mitotic differentiation of fibroblasts inducing a senescent or 'post-mitotic' phenotype associated with increased collagen production characteristic of the terminal differentiation pathway stimulated by ionising radiation [9, 10, 20]. In the current study, sparse cell cultures were treated for 6 hours with 10 μg/ml bleomycin sulphate on day 1 and analysis was performed at 72 hours to examine the transcriptional response of cells in the post-mitotic state. This dose of bleomycin sulphate has been previously shown to induce post-mitotic differentation in fibroblasts . Under these conditions, transcriptional changes of genes related to the expression of the differentiated phenotype, considered to be relevant to late normal tissue radiation injury, may be seen. However the fibroblast response to bleomycin sulphate did not confirm enrichment of relevant gene ontology categories. For example, for up-regulated genes, the EASE score for extracellular matrix was 3.6E-02 and for response to oxidative stress was 7.3E-02.
Published data report the dose levels of bleomycin sulphate required to produce cell cycle arrest but not cell cytotoxicity. Bleomycin sulphate is known to cause both G1/S and G2/M arrest . At high doses, extensive double strand DNA breaks and apoptosis occur . We aimed to use a dose at which post-mitotic differentiation and cellular response pathways were induced but without causing lethality with predominant death signals. Preliminary experiments carried out in this study using FACS analysis confirmed that the dose applied resulted in predominantly G2/M arrest without significant cell death [see Additional File 2]. Using these experimental conditions, gene expression data did indeed show changes in levels of genes relevant to cell cycle control. Cell death pathway activation was also seen, and this may have contributed to the inability to differentiate radiation sensitive and control groups.
In conclusion, a difference between radiation sensitive cases and matched controls was not detected in this population of breast cancer trial patients who had prospective scoring of late normal tissue effects. This suggests any difference is likely to be small or the variation between patients too great to detect a difference. Limitations of the clinical trial design and the experimental laboratory design could have been contributory.
Dulbecco's modified eagle medium
foetal calf serum
fluorescence activated cell sorting
quantitative real time polymerase chain reaction
significance analysis of microarrays
Rieger KE, Chu G: Portrait of transcriptional responses to ultraviolet and ionizing radiation in human cells. Nucleic Acids Res 2004,32(16):4786-4803. 10.1093/nar/gkh783
Rodningen OK, Overgaard J, Alsner J, Hastie T, Borresen-Dale AL: Microarray analysis of the transcriptional response to single or multiple doses of ionizing radiation in human subcutaneous fibroblasts. Radiother Oncol 2005,77(3):231-240. 10.1016/j.radonc.2005.09.020
Tachiiri S, Katagiri T, Tsunoda T, Oya N, Hiraoka M, Nakamura Y: Analysis of gene-expression profiles after gamma irradiation of normal human fibroblasts. Int J Radiat Oncol Biol Phys 2006,64(1):272-279. 10.1016/j.ijrobp.2005.08.030
Zhou T, Chou JW, Simpson DA, Zhou Y, Mullen TE, Medeiros M, Bushel PR, Paules RS, Yang X, Hurban P, et al.: Profiles of global gene expression in ionizing-radiation-damaged human diploid fibroblasts reveal synchronization behind the G1 checkpoint in a G0-like state of quiescence. Environ Health Perspect 2006,114(4):553-559.
Quarmby S, West C, Magee B, Stewart A, Hunter R, Kumar S: Differential expression of cytokine genes in fibroblasts derived from skin biopsies of patients who developed minimal or severe normal tissue damage after radiotherapy. Radiat Res 2002,157(3):243-248. 10.1667/0033-7587(2002)157[0243:DEOCGI]2.0.CO;2
Rodningen OK, Borresen-Dale AL, Alsner J, Hastie T, Overgaard J: Radiation-induced gene expression in human subcutaneous fibroblasts is predictive of radiation-induced fibrosis. Radiother Oncol 2008,86(3):314-320. 10.1016/j.radonc.2007.09.013
Svensson JP, Stalpers LJ, Esveldt-van Lange RE, Franken NA, Haveman J, Klein B, Turesson I, Vrieling H, Giphart-Gassler M: Analysis of gene expression using gene sets discriminates cancer patients with and without late radiation toxicity. PLoS Med 2006,3(10):e422. 10.1371/journal.pmed.0030422
Povirk LF, Austin MJ: Genotoxicity of bleomycin. Mutat Res 1991,257(2):127-143.
Koslowski R, Morgner J, Seidel D, Knoch KP, Kasper M: Postmitotic differentiation of rat lung fibroblasts: induction by bleomycin and effect on prolyl 4-hydroxylase. Exp Toxicol Pathol 2004,55(6):481-487. 10.1078/0940-2993-00343
Robles SJ, Adami GR: Agents that cause DNA double strand breaks lead to p16INK4a enrichment and the premature senescence of normal fibroblasts. Oncogene 1998,16(9):1113-1123. 10.1038/sj.onc.1201862
Rodemann HP, Binder A, Burger A, Guven N, Loffler H, Bamberg M: The underlying cellular mechanism of fibrosis. Kidney Int Suppl 1996, 54: S32-36.
Cutroneo KR: Evidence for TGF-beta1 and bleomycin intracellular signaling through autocrine regulation of Smad 3 binding to the proximal promoter of the Smad 7 gene. J Cell Biochem 2006,97(5):933-939. 10.1002/jcb.20594
Rubin P, Johnston CJ, Williams JP, McDonald S, Finkelstein JN: A perpetual cascade of cytokines postirradiation leads to pulmonary fibrosis. Int J Radiat Oncol Biol Phys 1995,33(1):99-109. 10.1016/0360-3016(95)00095-G
Yamamoto T, Takagawa S, Katayama I, Yamazaki K, Hamazaki Y, Shinkai H, Nishioka K: Animal model of sclerotic skin. I: Local injections of bleomycin induce sclerotic skin mimicking scleroderma. J Invest Dermatol 1999,112(4):456-462. 10.1046/j.1523-1747.1999.00528.x
Yarnold J, Ashton A, Bliss J, Homewood J, Harper C, Hanson J, Haviland J, Bentzen S, Owen R: Fractionation sensitivity and dose response of late adverse effects in the breast after radiotherapy for early breast cancer: long-term results of a randomised trial. Radiother Oncol 2005,75(1):9-17. 10.1016/j.radonc.2005.01.005
Peacock J, Ashton A, Bliss J, Bush C, Eady J, Jackson C, Owen R, Regan J, Yarnold J: Cellular radiosensitivity and complication risk after curative radiotherapy. Radiother Oncol 2000,55(2):173-178. 10.1016/S0167-8140(00)00173-0
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA 2001,98(9):5116-5121. 10.1073/pnas.091062498
Andreassen CN, Alsner J, Overgaard J, Herskind C, Haviland J, Owen R, Homewood J, Bliss J, Yarnold J: TGFB1 polymorphisms are associated with risk of late normal tissue complications in the breast after radiotherapy for early breast cancer. Radiother Oncol 2005,75(1):18-21. 10.1016/j.radonc.2004.12.012
Alsner J, Rodningen OK, Overgaard J: Differential gene expression before and after ionizing radiation of subcutaneous fibroblasts identifies breast cancer patients resistant to radiation-induced fibrosis. Radiother Oncol 2007,83(3):261-266. 10.1016/j.radonc.2007.05.001
Herskind C, Johansen J, Bentzen SM, Overgaard M, Overgaard J, Bamberg M, Rodemann HP: Fibroblast differentiation in subcutaneous fibrosis after postmastectomy radiotherapy. Acta Oncol 2000,39(3):383-388. 10.1080/028418600750013159
Kaufmann WK, Kies PE: DNA signals for G2 checkpoint response in diploid human fibroblasts. Mutat Res 1998,400(1-2):153-167.
Chen J, Stubbe J: Bleomycins: towards better therapeutics. Nat Rev Cancer 2005,5(2):102-112. 10.1038/nrc1547
We thank Yvonne Hey and Stuart Pepper at the Cancer Research UK Affymetrix Genechip Microarray Service, Paterson Institute for Cancer Research, for performing the Affymetrix Gene Arrays and for technical support and advice and Henrik Edgren at Institute for Molecular Medicine, University of Helsinki, Finland for help with the Affymetrix gene normalisation. We also thank Carsten Herskind at the Department of Radiation Oncology, University Medical Center Mannheim, University of Heidelberg, Germany for helping with the experimental design and Jorge Reis-Filho and Jo Haviland at the Institute of Cancer Research for useful discussion. This work was supported by Cancer Research UK Section of Radiotherapy (CRUK) grant number C46/A2131 and Breakthrough Breast Cancer, the Institute of Cancer Research, UK. We acknowledge the support of the Cancer Research UK Affymetrix Genechip Microarray Service and NHS funding to the NIHR Biomedical Research Centre.
The authors declare that they have no competing interests.
CBW participated in the study design and statistical analysis, carried out cell and RNA preparation and drafted the manuscript. KKS participated in study design, performed statistical analysis and helped to draft the manuscript. ALBD particiapted in aspects of the study design and statistical analysis. JRY was responsible for conceiving the case control design and participated in the study design. CMI participated in the study design and study coordination. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Details of Royal Marsden Hospital/Gloucester Oncology Centre Breast Radiotherapy Fractionation Trial (1986-1998) and scoring of late normal tissue effects. The table includes details of the Royal Marsden Hospital/Gloucester Oncology Centre Breast Radiotherapy Fractionation Trial (1986-1998), methods of assessment of late normal tissue injury and categorisation of patients into cases with late radiation effects and controls with minimal/none. (DOC 34 KB)
Additional file 2:Cell cycle analysis of adult human dermal fibroblasts after treatment with bleomycin sulphate. Cell cycle analysis of adult human skin fibroblasts after treatment with bleomycin sulphate was carried out. Fibroblasts were plated at a density of 2 × 104 cells in a T25 cm2 flask and treated after 1 day in culture with medium alone (control), or with 10 μg/ml or 50 μg/ml bleomycin sulphate for 6 hours or 24 hours as indicated. After 4 days in culture, cells were stained with propidium iodide and analysed by fluorescence activated cell sorting (FACS); a) representative FACS profiles. b) the percentage of cells in each phase of the cell cycle determined using the Watson Pragmatic model. Both doses resulted in accumulation of fibroblasts with 4N-DNA content and therefore the lower dose (10 μg/ml) was used for the treatment of all the experimental samples. (PPT 180 KB)
Additional file 3:Genes differentially regulated between bleomycin sulphate treated and mock treated fibroblasts. Paired SAM analysis of bleomycin sulphate treated and mock treated samples was carried out using 16 fibroblasts cultures (8 radiation sensitive cases and 8 matched controls). 973 genes were up-reulated and 923 genes were down-regulated in bleomycin sulphate treated compared to mock treated fibroblasts. (XLSX 231 KB)
About this article
Cite this article
Westbury, C.B., Sahlberg, K.K., Borresen-Dale, AL. et al. Gene expression profiling of human dermal fibroblasts exposed to bleomycin sulphate does not differentiate between radiation sensitive and control patients. Radiat Oncol 6, 42 (2011). https://doi.org/10.1186/1748-717X-6-42