Open Access Research

Identification of stable endogenous control genes for transcriptional profiling of photon, proton and carbon-ion irradiated cells

Geeta D Sharungbam12, Christian Schwager12, Sara Chiblak12, Stephan Brons2, Lynn Hlatky3, Thomas Haberer2, Jürgen Debus2 and Amir Abdollahi123*

Author Affiliations

1 Molecular RadioOncology [E210], National Center for Tumor Disease (NCT), German Cancer Research Center (DKFZ), Im Neuenheimer Feld 460, 69120, Heidelberg, Germany

2 Heidelberg Ion Therapy Center (HIT), Heidelberg Institute of Radiation Oncology (HIRO), Department of Radiation Oncology, University of Heidelberg Medical School, Im Neuenheimer Feld 450, 69120, Heidelberg, Germany

3 Center of Cancer Systems Biology, NASA Specialized Center of Research, St. Elizabeth’s Medical Center, Tufts University, 736 Cambridge Street [CBR 1], 02135, Boston, MA, USA

For all author emails, please log on.

Radiation Oncology 2012, 7:70  doi:10.1186/1748-717X-7-70


The electronic version of this article is the complete one and can be found online at: http://www.ro-journal.com/content/7/1/70


Received:1 December 2011
Accepted:17 May 2012
Published:17 May 2012

© 2012 Sharungbam et al.; licensee 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.

Abstract

Background

Quantitative analysis of transcriptional regulation of genes is a prerequisite for a better understanding of the molecular mechanisms of action of different radiation qualities such as photon, proton or carbon ion irradiation. Microarrays and real-time quantitative RT-PCR (qRT-PCR) are considered the two cornerstones of gene expression analysis. In interpreting these results it is critical to normalize the expression levels of the target genes by that of appropriately selected endogenous control genes (ECGs) or housekeeping genes. We sought to systematically investigate common ECG candidates for their stability after different radiation modalities in different human cell lines by qRT-PCR. We aimed to identify the most robust set of ECGs or housekeeping genes for transcriptional analysis in irradiation studies.

Methods

We tested the expression stability of 32 ECGs in three human cancer cell lines. The epidermoid carcinoma cells (A431), the non small cell lung carcinoma cells (A549) and the pancreatic adenocarincoma cells (BxPC3) were irradiated with photon, proton and carbon ions. Expression Heat maps, clustering and statistic algorithms were employed using SUMO software package. The expression stability was evaluated by computing: mean, standard deviation, ANOVA, coefficient of variation and the stability measure (M) given by the geNorm algorithm.

Results

Expression analysis revealed significant cell type specific regulation of 18 out of 32 ECGs (p < 0.05). A549 and A431 cells shared a similar pattern of ECG expression as the function of different radiation qualities as compared to BxPC3. Of note, the ribosomal protein 18S, one of the most frequently used ECG, was differentially regulated as the function of different radiation qualities (p ≤ 0.01). A comprehensive search for the most stable ECGs using the geNorm algorithm identified 3 ECGs for A431 and BxPC3 to be sufficient for normalization. In contrast, 6 ECGs were required to properly normalize expression data in the more variable A549 cells. Considering both variables tested, i.e. cell type and radiation qualities, 5 genes-- RPLP0, UBC, PPIA, TBP and PSMC4-- were identified as the consensus set of stable ECGs.

Conclusions

Caution is warranted when selecting the internal control gene for the qRT-PCR gene expression studies. Here, we provide a template of stable ECGs for investigation of radiation induced gene expression.

Keywords:
Endogenous control genes; Internal control genes; qRT-PCR; Photon; Proton; Carbon-ion; Tumour cells; A431; A549 and BxPC3

Background

In addition to direct, e.g. DNA damaging effect, system level cellular responses to ionizing radiation are attributed to the initiation of intracellular signals and subsequent differential regulation of genes/pathways governing various cellular processes [1]. Therefore, detecting differential regulation of genes is critical for a better understanding of radiation-induced molecular effects. Transcriptional perturbation after cell exposure to different radiation qualities is investigated to unravel the systems biology of cellular response underlying, normal tissue toxicity, carcinogenesis, or anti-cancer effects of irradiation [1-3]. Therefore, these studies have ramification for a broad spectrum of basic and applied sciences ranging from effects of space radiation to carcinogenesis to cancer therapy.

In contrast to conventional photon irradiation the molecular effects of proton or heavier ions (e.g. carbon ions) are less explored yet. However, emerging data indicate molecular differences in transcriptional response of cells to particles as compared to photon irradiation [4,5].

One reliable and highly sensitive tool that allows rapid and accurate results in gene expression analysis is the qRT-PCR [6,7]. As in any gene expression analysis, selection of a valid normalization or endogenous control to correct for differences in RNA sampling is critical to avoid misinterpretation of results. Inter-sample variation due to sample collection, RNA preparation and quality, inherent sample differences, pipetting errors, different efficacies of the radiation qualities and reverse transcription efficiency are common sources of variability. The ideal endogenous control should have a constant expression level under different experimental conditions and be sufficiently abundant across different samples and cell lines. Although any gene that is stably expressed under a defined experimental condition can be used for normalization, the selection is most commonly made from the constitutively expressed ECGs.

However, the expression levels of the commonly used ECGs may not only vary in different cell lines but also under different experimental treatments or pathological states [8-25]. This necessitates the selection of ECGs which are appropriate for each experimental system. Although, there has been systematic selection of ECGs for various experimental systems, such selection has not been conducted so far for studying the effects of different radiation qualities.

Here, we investigate the expression stability of 32 commonly used ECGs in three human cancer cell lines irradiated with photon, proton and carbon ions. Differential regulation of ECGs was found as the function of both variables, radiation quality and cell type, respectively. Reliable internal control genes for individual cell lines were identified such as PGK1, RPL37A and PSMC4 for A431; RPLPO, UBC, GAPDH, MT-ATP6, CASC3 and PES1 for A549; and RPL37A, RPLPO and CASC3 for BxPC3. A systematic analysis further revealed 5 stable genes among the 32 candidate ECGs tested to normalize gene expression data generated in different cells and after various radiation qualities.

Results

Expression of the 32 ECGs

In this study, 32 ECGs (Additional file 1) were evaluated to identify the most suitable reference genes for gene expression profiling of irradiated cell lines. This collection of genes constitutes frequently used ECGs which were selected based on their relative high abundance and constitutive expression determined by literature search and/or whole genome microarray data. The three prototypic tumour cell lines-- A431, A549 and BxPC3 used in this study are among most commonly investigated model cell lines for each tumor entity. They were irradiated with photon, proton and carbon ions. After total-RNA isolation and quality control using lab-on-chip bioanalyzer, qRT-PCR was performed using Taqman primer and probes.

Additional file 1. List of 32 endogenous control genes used in the study.

Format: PDF Size: 77KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

CT-values and ECG regulation

To get a better overview of the CT-values among all the cell lines, the CT-range along with the minimum and maximum CT-values were listed in Table 1. The variation of CT-values ranged from 8.97 in A431 to 28.58 in A549. The 8.97 CT-value corresponded to 18S indicating its high abundance in the samples, whereas, the CT-value 28.58 corresponded to GADD45A indicating a moderate abundance. Moreover, RPL37A, GAPDH and RPLPO exhibited small CT-range indicating less variation in expression whereas the large CT-range of ACTB, GADD45A, IPO8 showed large variations in their expression.

Table 1. Cycle threshold (CT) values and coefficient of variation (CV) of 32 endogenous control genes across the samples of the cell lines

The CT-range and coefficient of variation for each cell line were listed in additional file 2. Within each individual cell line irradiated with photon, proton and carbon ion, all the ECGs in BxPC3 except for CDKN1A, 18S, POLR2A, PES1 and HMBS with CT-range 1.18, 0.82, 0.80, 0.70 and 0.61, respectively, showed the smallest CT-range, while in A549 all the ECGs except PES1, ACTB and RPS17 showed the largest CT-range (Figure 1A). In other words, as a function of different radiation qualities in BxPC3, the expression of the ECGs was considerably stable, while there was more variation in A549. Moreover, A431 and A549 shared similar expression pattern of the ECGs. The coefficient of variation plotted in Figure 1B also reflected the similar regulation of the 32 ECGs in the three cell lines.

Additional file 2. Cycle threshold range and coefficient of variation (CV) of 32 ECGs in each cell line. The genes are sorted by the coefficient of variation increasing from top to bottom. Four replicates were used in all the three cell lines.

Format: PDF Size: 61KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 1. Radiation induced variation in ECG expression. (A) CT range and (B) coefficient of variation Figure shows the variations in the expression level of each ECG in A431, A549 and BxPC3. BxPC3 showed least variation as a function of different radiation qualities as compared to A431 and A549. A549 showed highest variations.

ECG regulation using heat map

The expression of the 32 candidate ECGs across the samples were also visualized in a Heat map (Figure 2). Direct clustering of the ECGs expression showed that the expression profile of A431 and A549 were more similar as compared to BxPC3. BxPC3 showed least variation whereas A549 showed maximum variation in expression among the samples as well as among the ECGs indicating a cell type specific expression of ECGs. As observed in Table 1, the Heat map also revealed low variation of RPL37A, RPLPO and GAPDH expression, while ACTB, IPO8 and GADD45A depicted high variation in expression. The remaining ECGs were differentially regulated. For example, B2M and ELF1 were up regulated in BxPC3 and A431 but down regulated in A549. On the other hand, MRPL19 and PES1 were down regulated in BxPC3 but up regulated in A431 and A549.

thumbnailFigure 2. Direct clustering of 32 ECGs expression. This Heat map represents expression of all 32 ECGs across the three cell lines (A431, A549 and BxPC3) irradiated with photon, proton and carbon ion. Genes were hierarchically clustered by Pearson correlation coefficient using average linkage. Green denotes genes with relatively decreased expression while red denotes genes with relatively increased expression. Scale bar represent log2 expression level of ECGs. C = Carbon, P = Proton, X = Photon and 0 = Control. Expression profile of A431 and A549 are similar as compared to BxPC3 which showed unique expression profile with less variations among the samples.

In Figure 3, we displayed a Heat map generated from a one way ANOVA analysis at p < 0.05 between the three cell lines. It revealed 18 differentially regulated ECGs-- ABL, CDKN1A, PSMC4, EIF2B1, GAPDH, PPIA, TBP, RPS17, UBC, B2M, ELF1, PUM, GADD45A, ACTB, MRPL19, YWHAZ, CASC3 and PES1. These could be cell type specific regulations as their radiation quality variation was minimal. Among them, the first 13 except for PSMC4 were up regulated in BxPC3 and the remaining down regulated. In A431 and A549, 12 of them were down regulated and 6 were up regulated. This supports the finding that the expression levels of ECGs are different in BxPC3 as compared to A431 and A549. BxPC3 showed least variations and A549 showed maximum variations in gene expression as the function of radiation qualities. However, expression levels of ECGs were comparable in A431 and A549 (Figure 1A).

thumbnailFigure 3. Cell type specific regulation of ECGs. This Heat map represents 18 ECGs which are significantly regulated according to ANOVA between the three cell lines, p < 0.05. Relative expression of each gene were normalized to the average intensity of the gene over entire samples (virtual pool). Green denotes genes with relatively decreased expression while red denotes genes with relatively increased expression. Genes are hierarchically clustered by Pearson correlation coefficient using average linkage. Scale bar represent log2 expression level of ECGs. C = Carbon, P = Proton, X = Photon and 0 = Control. Expression levels of the ECGs in BxPC3 is different as compared to A431 and A549. A549 showed maximum variation among the samples.

ANOVA analysis between the three radiation qualities at p < 0.01 revealed one gene, the ribosomal protein 18S, to be differentially regulated after different radiation qualities (Figure 4). This gene is one of the most commonly used internal control genes for normalisation of qRT-PCR based gene expression data. Therefore, caution needs to be practised in using this gene as an internal control gene, in particular when radiation effects are investigated.

thumbnailFigure 4. Radiation induced differential regulation of ECG. Figure displays gene expression levels of 18 S in the three cell lines (A431, A549 and BxPC3) irradiated with photon, proton and carbon ion. Among the 32 ECGs, one gene; the ribosomal protein 18 S was found to be differentially regulated as the function of different radiation qualities (p < 0.01 by ANOVA). Bars indicate mean expression ± standard deviation.

The results in this section corroborated the findings of many other previous studies that the ECGs might be differentially regulated depending upon the experimental set-up and the cell type [8-25]. More importantly, this analysis seems to suggest that the ECGs are differentially regulated by the different cell types and radiation qualities. We attempted to confirm this observation employing a systematic analysis of the expression levels.

Identification of appropriate ECGs

Gene expression levels obtained using PCR should be appropriately normalized by one or more carefully selected stable internal control genes. The geNorm algorithm developed by Vandesompele et al. [26] can determine the expression stability of control genes on the basis of non-normalized expression levels. This measure relies on the principle that the expression ratio of two internal control genes is constant in all samples regardless of the experimental condition or cell type. This algorithm computes a gene expression stability measure (M) for each gene based on the average pairwise expression ratio and then performs a stepwise exclusion of the least stable gene. Then the M values are computed again and stepwise exclusion performed until two genes are left. The genes with the lowest M values are considered to be the most stable across all the samples for each cell line.

Ranking of the 32 ECGs

The M values for all the 32 ECGs in A431, A549 and BxPC3 computed using the geNorm algorithm (integrated into SUMO software) were sorted and ranked in Table 2. This table revealed that the two most stable ECGs irrespective of the radiation qualities were: PGK1-RPL37A for A431, RPLPO-UBC for A549, and RPL37A-RPLPO for BxPC3.

Table 2. Control genes ranked in order of their expression stability*

The actual stepwise exclusion of the worst-scoring ECG was displayed in Figure 5. In this figure: (1) there was a very steep decrease in the average M value for A431 pointing at two unstably expressed ECGs, (2) the irregular decrease in the average M value for A549 showed that it has aberrantly expressed ECGs, (3) the regular decrease in the average M value for BxPC3 might mean that all the ECGs were stable.

thumbnailFigure 5. Identification of the most stable ECGs in each cell line. Most stable ECGs, i.e. not differentially regulated by different radiation qualities in each cell lines, are identified. Average expression stability M of all remaining control genes after stepwise exclusion of the least stable reference genes in three cell lines are shown. More stably expressed genes are positioned on the right side of the diagram, less stably expressed on the left side. ECGs are ranked in order of their expression stability and presented along x-axis. Stability values (M) determined by geNorm algorithm are presented along y-axis. Low stability value (M) reflects greater stability. For the gene names with their ranking refer Table 2.

Calculation of normalization factor

For each cell line, the normalization factors (NF) were computed, first for the three most stable ECGs, by taking the geometric mean of their expression levels. This is followed by stepwise inclusion of the most stable remaining ECG. Then the pairwise variations Vn(n+1) were calculated for every series of NFn and NFn+1, reflecting the effect of adding an (n + 1)th ECG (Figure 6) [26].

thumbnailFigure 6. Determination of the optimal number of control genes for normalization. Pairwise variation (Vn(n+1)) analysis between the normalization factors (NFn) and (NFn+1) to determine the number of control genes required for accurate normalization. Normalization factors were computed taking 3–11 most stable genes for all the three cell lines. Pairwise variation of 0.15 was taken as a cut off value [26]. For A431 and BxPC3 three ECGs were sufficient for normalization in contrast to six ECGs for A549.

Figure 6 shows that the value of V3/4 was low for A431, implying that the first 3 ECGs (PGK1, RPL37A, PSMC4) were sufficient to be used for normalization. For A549, the low value of V6/7 indicated that the first 6 ECGs (RPLPO, UBC, GAPDH, MT-ATP6, CASC3, PES1) were sufficient for normalization. In BxPC3, the three most stable ECGs (RPL37A, RPLPO, CASC3) were sufficient for normalization purposes.

Validation of the gene-stability measure M

According to Vandesompele et al. [26], three different normalization factors were calculated based on the geometric mean of three genes with, respectively, the smallest M value (NF3(1–3)), the intermediate M value (NF3(11–13)) and the highest M value (NF3(30–32)) as determined by geNorm (Table 2). Further, we determined the average gene-specific variation of the three genes with the most stable expression (i.e., the smallest coefficient of variation) for each normalization factor within each cell line (Figure 7).It is conceivable that the gene-specific variation in all the cell lines were the least when the data are normalized to (NF3(1–3)). This validated that the gene-stability measure effectively identified the ECGs with the most stable expression.

thumbnailFigure 7. Validation of the gene stability measure (M) and the geometric averaging of carefully selected control genes for normalization. The average gene-specific variation (determined as coefficient of variation in percent) for the three control genes with the smallest variation within each cell line after normalization with three different factors calculated as the geometric mean of the three control genes with the lowest (NF3(1–3)), intermediate (NF3(11–13)) and highest (NF3(30–32)) gene-stability values (as determined by geNorm). The data confirmed the stability of the ECGs. The value under the star indicates the normalization factor for each cell line.

Radiation-specific expression of the ECGs within the cell lines

The gene stability measure M value was determined and validated, the candidate ECGs within each cell line were normalized by the appropriate stable ECGs found above and plotted in Figure 8. The variation in regulation of 18S indicated that its expression depended upon the radiation quality (Figure 8A). This observation is in line with other above mentioned analysis performed (Figures 2 and 4). In addition, GADD45A was differentially regulated by radiotherapy in all three cell lines (Figure 8). Together, these data confirm differential regulation of candidate control genes as a function of different radiation qualities.

thumbnailFigure 8. Differential regulation of the ECGs after normalization with stable genes of each cell line. A) A431 normalized by PGK1, RPL37A, and PSMC4; (B) A549 normalized by RPLPO, UBC, GAPDH, MTATP6, CASC3 and PES1; (C) BxPC3 normalized by RPL37A, RPLPO and CASC3. As compared to the raw intensity data the differential regulation of the 32 ECGs are more pronounced after normalization with the identified stable genes for each cell line with correct normalization factor (NFn).

Identifying the consensus set of ECGs for comparative studies across all cell types

To compare the ECGs expression levels across A431, A549 and BxPC3 cells, first a consensus set of ECGs was identified for normalization of expression data using the algorithm suggested by Vandesompele et al. [26]. RPLPOUBCPPIATBP and PSMC4 were selected by eliminating the ECGs with high M value to normalize and compare the cell type specific gene-expressions (Figure 9). Although the overall abundance of most ECGs among different cell lines was relatively similar, cell-line specific gene-expression were identified for some candidate ECGs such as, 18SB2MYWHAZPGK1CDKN1A and GADD45A. In contrast, ECGs with a relatively constant expression included GAPDHRPLPORPL30APPIAUBC etc. In A431, a 422-fold expression difference was observed between the most stable gene (PGK1) and the least stable gene (18S) whereas a 530 and 375 fold difference in expression was found in A549 and BxPC3, respectively.

thumbnailFigure 9. Logarithmic histogram of the expression levels of 32 ECGs in all the cell lines. The 32 ECGs were normalized to the geometric mean of five control genes (RPLPO, UBC, TBP, PPIA and PSMC4). In A431, a 422-fold expression difference is observed between the most stable gene (PGK1) and the least stable gene (18 S), while in A549 and BxPC3, a fold difference of 530 and 375 respectively were observed between them. In addition it shows the cell line specific differences in expression levels of particular genes (e.g. YWHAZ).

Discussion

The emergence of a growing number of particle therapy facilities worldwide will stimulate comparative studies aiming to decipher the molecular mechanisms underlying differential biological effects of these novel radiation qualities. Comparative investigation of gene regulation on transcriptional level as the function of radiation treatment constitutes a cornerstone of these studies. Quantitative real time PCR (qRT-PCR) is considered the most sensitive method for detection of gene expression level. One limitation of this method is the need for proper endogenous control gene. To generate relative expression levels, the expression of the reference gene/s needs minimally alter among different types of cells or treatments. The goal of this study was to identify such ECGs.

We analysed the expression levels of 32 ECGs using the clustering, statistical methods such as ANOVA and the geNorm algorithm. Global analysis lead to the finding that gene expression profile in pancreatic cancer cells (BxPC3) is different as compared to the two other epithelial cancer cells tested i.e. epidermoid and lung carcinoma cells (A431 and A549). The ECGs in BxPC3 showed least variation in expression whereas A549 showed maximum variation in expression as the function of radiation qualities. Among the three cell lines, the ECGs were more stable in BxPC3. From the point of view of selecting appropriate ECGs this feature might be advantageous. On the other hand, it means that the ECGs in BxPC3 are less regulated by different radiation qualities.

ANOVA analysis of non normalized data revealed that 18 out of the 32 ECGs showed cell type specific differential regulation as shown by the differences in the expression profile between A431, A549 and BxPC3. In addition, significant radiation quality specific regulation was shown by one gene the ribosomal protein 18S. Of note, this gene is one of the most commonly used internal control genes. Therefore, our data suggest validation of this gene prior to its use as internal control in radiation biology studies.

Although the clustering and ANOVA analysis of the raw data provided an overall overview and information about ECGs regulation, they do not allow the selection of appropriate stable ECGs for normalization of the qRT-PCR data. Next, the geNorm algorithm was employed to test for the stability of the 32 candidate ECGs as reference genes as well as selection of the optimal number of genes for normalization of gene expression.

Using the geNorm algorithm the two most stable ECGs for each cell line were determined- PGK1-RPL37A in A431, RPLPO-UBC in A549, and RPL37A-RPLPO in BxPC3. Further the optimal number of ECGs for the normalization of gene expression in each cell line was determined and validated: three ECGs each for A431 (PGK1, RPL37A, and PSMC4) and BxPC3 (RPL37A, RPLPO and CASC3) were recommended. In contrast, six ECGs for A549 data (RPLPO, UBC, GAPDH, MTATP6, CASC3 and PES1) were required.

The normalized gene expression for each cell line in Figure 8 showed that 18S and B2M are unstable genes under different radiation qualities. This is in contrast to earlier studies using 18S and B2M as reference genes for proton therapy [4]. Besides, PPIAACTB and UBC for particle therapy using 0.5 Gy 4He ions in normal human lung fibroblasts [27] and ACTB for A549 [5] were reported as reference genes. However, Table 2 showed that these genes have intermediate stability within the cell line examined here.

Figure 8 demonstrates the regulation of particular genes treated with different radiation qualities. For instance, Figure 8B showed that in A549, CDKN1A-- a gene downstream of p53 pathway which is also implicated in regulation of cell growth and cell response to DNA damage-- is up regulated under all the radiation qualities, with maximum up regulation under photon. The gene YWHAZ-- involved in signal transduction by binding to phosphorylated serine residues on a variety of signaling molecule-- is up regulated in A431- photon while minimally regulated in A549-photon. GADD45A-- which binds to proliferating cell nuclear antigen, stimulates DNA excision repair in vitro and inhibits entry of cells into S phase-- was down regulated in A431 but up regulated in BxPC3. 18S-- a component of the ribosome, the protein manufacturing machinery of all living cells-- is seen to be up regulated in carbon ion while its regulation varies for proton and photon in different cell lines.

In addition, five most stable ECGs (RPLPOUBCTBPPPIA and PSMC4) in three cell lines were selected as internal control genes for the normalisation of the gene expression independent of radiation qualities and cell type. This selection was based on the guideline of Vandesompele et al. [26].

The gene expression in each cell lines normalized by the selected five stable ECGs was shown in Figure 9. The expression of ACTB showed 2.8-fold difference between the highest and lowest expression levels, whereas YWHAZ, 18S, GADD45A showed 5.8, 3.14 and 4-fold difference between the highest and lowest expression levels. The expression Heat map (Figure 3) also illustrate the regulation of YWHAZ, GADD45A and 18S.

Selecting the ECGs for normalization across all the cell lines is a subtle issue. Although in [26], an algorithm to select the best ECGs within each specific cell line is presented, a clear method of selecting the best ECGs for all the cell lines is not given. More precisely, among the 5 selected ECGs (Table 3) -- TBP, UBC, RPLPO, PPIA, PSMC4-- one ECG could be stable in one cell line, while it could be of intermediate stability in the other. However, our selection is supported by the fact that, in Table 2, all these genes are of intermediate stability in each of the cell lines.

Table 3. List of ECGs qualified as internal control genes across the cell and radiation qualities

Conclusions

Careful selection and validation of ECGs prior to conducting radiation biology experiment is warranted. We report that different radiation qualities induced differential regulation of a number of ECGs among the candidate 32 “housekeeping genes”. Additional cell type specific gene expression was observed. Identification of the best internal control gene is a prerequisite for a successful quantitative measure of gene expression via RT-PCR. In this paper we provide a template for the identification of appropriate ECGs for the radiation induced gene expression studies. We identified reliable genes for individual expression profiling of the cell lines, the normalization of A431 may be done by PGK1, RPL37A and PSMC4; A549 by RPLPO, UBC, GAPDH, MT-ATP6, CASC3 and PES1; and BxPC3 by RPL37A, RPLPO and CASC3. However, the 5 ECGs-- TBP, UBC, RPLPO, PPIA, PSMC4-- can be taken as the most suitable candidate reference genes for radiation response expression profiling in the tumor models studied. Moreover, this robust set of the most suitable candidate ECGs for radiation experiment may be applied and validated for the clinicopathological analysis of cancer specimens of epithelial tumors, non-small cell lung cancer and pancreatic adenocarcinoma.

Methods

Cell lines

The three different human tumour cell lines, i.e., the lung carcinoma cells (A549), the epidermoid carcinoma cells (A431) and pancreatic cancer cells (BxPC3) were used for the study. A549 and A431 cell lines were obtained from Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH (DSMZ) and BxPC3 from the American Type Culture Collection. The A549 and A431 cell lines were grown in 5 ml Dulbeccos Modified Eagle's Medium (DMEM) (Biochrom), BxPC3 was grown in 5 ml RPMI 1640 medium (GIBCO Invitrogen) supplemented with 10.0% FCS in T25 flasks (Becton Dickinson). Cells were cultured under standard conditions in a fully humidified incubator with 5.0% CO2 at 37.0°C.

Irradiation

Cells were irradiated in T25 flasks with 2Gy of photon, 2Gy of proton and 1Gy of carbon ion. Photon was delivered by a linear accelerator at 6 Mev (Mevatron, Siemens, Erlangen, Germany). Particle irradiation with proton and carbon ion was done using a pencil beam in a spread out Bragg peak with 1.5 cm width equivalent to a depth of 14.0 cm in water, at the Heidelberg Ion Therapy Center (HIT) [28]. After irradiation, the cells were incubated for 12 h at 37.0°C. Control cells were treated identically but without irradiation. Cells were scrapped using the cell scraper after adding 300.0 μl TRIzol (Invitrogen) and collected in 1.5 ml Eppendorf tubes and subsequently stored at −20.0°C.

RNA isolation and cDNA synthesis

RNA was isolated in phase lock tubes using TRIzol (Invitrogen) according to the manufacturer’s protocol. To avoid genomic DNA contamination RNA was treated with Dnase I (Ambion). Purified RNA was eluted in 20.0μL of nuclease-free water and stored at −20.0°C. RNA concentration and purity was assessed using a Nanodrop ND-1000 spectrophotometer (Peqlab). Integrity and concentration of RNA samples were determined by using RNA 6000 Nano Lab Chip kits and a 2100 Bioanalyzer (Agilent). RNA (2.0 μg) was subjected to reverse transcription reaction using the high-capacity cDNA reverse transcription kit (Applied Biosystems) according to the manufacturer’s protocol.

Real-time PCR

To assess the expression of Human Endogenous Control gene set, real-time quantitative reverse transcription PCR (qRT-PCR) was performed on 32 candidate genes using TaqMan® chemistry (Applied Biosystems). Experiments were performed at least in triplicates for the three cell lines-- A431, A549 and BxPC3. Plates were run on a 7900HT Fast Real-Time PCR Systems (Applied Biosystems) using Fast 96-well blocks, Standard Fast PCR cycling protocol with 10.0μL reaction volumes. Cycling conditions used were-- 1 cycle initiation at 50.0°C for 2 min and 1 cycle at 95.0°C for 10 min, followed by amplification for 40 cycles at 95.0°C for 15 s and 60.0°C for 1 min. Amplification data were collected via Sequence Detection Systems 2.3 software (Applied Biosystems). The CT-values were computed with RQ Manager 2.xx (Applied Biosystems).

Statistical analyses

Statistical analysis of data was performed using SUMO software package (http://www.oncoexpress.de/software/sumo webcite). ANOVA was used to detect variation in the expression of the ECGs across the samples according to the radiation qualities and cell lines respectively. The average expression stability measure values (M) were computed using the geNorm algorithm suggested by Vandesompele et al. [26] (also incorporated in the SUMO program package).

Competing interests

The authors declare that they have no competing interests.

Authors contributions

GDS and AA designed the experiment, performed research, analysed data and wrote the manuscript. CS developed software and performed data analysis and statistics. SC, LH and JD analysed data and edited the manuscript. SB and TH performed the heavy ion irradiation planning and treatment. All authors read and approved the final manuscript.

Acknowledgements

We thank Claudia Rittmüller, Christiane Rutenberg and Barbara Schwager for the excellent technical assistance. This work was supported in part by the German Krebshilfe (Deutsche Krebshilfe, Max-Eder 108876), DFG National Priority Research Program: the Tumour-Vessel Interface “SPP1190”, NASA Specialized Center of Research NNJ04HJ12G, and the German Federal Ministry of Research and Technology (Bundesministerium für Bildung und Forschung – BMBF 03NUK004C).

References

  1. Abdollahi A, Folkman J: Evading tumour evasion: current concepts and perspectives of anti-angiogenic cancer therapy.

    Drug Resist Updat 2010, 13:16-28. OpenURL

  2. Abdollahi A, Li M, Ping G, Plathow C, Domhan S, Kiessling F, Lee LB, McMahon G, Groene HJ, Lipson KE, Huber PE: Inhibition of platelet-derived growth factor signaling attenuates pulmonary fibrosis.

    J Exp Med 2005, 201:925-935. OpenURL

  3. Abdollahi A, Schwager C, Kleef J, Esposito I, Domhan S, Peschke P, Hauser K, Hahnfeldt P, Hlatky L, Debus J, Peters JM, Friess H, Folkmann JM , Huber PE: Transcriptional network governing the angiogenetic switch in human pancreatic cancer.

    Proc Natl Acad Sci 2007.

    doi :

    Publisher Full Text OpenURL

  4. Suzuki Y, Nakano T, Ohno T, Oka K: Comparison of the radiobiological effect of carbon ion beam therapy and conventional radiation therapy on cervical cancer.

    J Radiat Res 2007, 49:473-479. OpenURL

  5. Akino Y, Teshima T, Kihara A, Kodera-Suzumoto Y, Inaoka M, Higashiyama S, Furusawa Y, Matsuura N: Carbon-ion beam irradiation effectively suppresses migration and invasion of human non small-cell lung cancer cells.

    Int J Radiat Oncol 2009, 75:475-481. OpenURL

  6. Domhan S, Muschal S, Schwager C, Morath C, Wirkner U, Ansorge W, Maercker C, Zeier M, Huber PE, Abdollahi A: Molecular mechanisms of the antiangiogenic and antitumor effects of mycophenolic acid.

    Mol Cancer Ther 2008, 7:1656-1668. OpenURL

  7. Almog N, Ma L, Raychowdhury R, Schwager C, Erber R, Short S, Hlatky L, Vajkoczy P, Huber PE, Folkman J, Abdollahi A: Transcriptional switch of dormant tumors to fast-growing angiogenic phenotype.

    Cancer Res 2009, 69:836-844. OpenURL

  8. Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR.

    Biotechniques 2004, 37:112-119. OpenURL

  9. Haberhausen G, Pinsl J, Kuhn C, Markert-Hahn C: Comparative study of different standardization concepts in quantitative competitive reverse transcription-PCR assays.

    J Clin Microbiol 1998, 36:628-633. OpenURL

  10. Hamalainen H, Tubman J, Vikman S, Kyrola T, Ylikoski E, Warrington JA, Lahesmaa R: Identification and validation of endogenous reference genes for expression profiling of T helper cell differentiation by quantitative real-time RT-PCR.

    Anal Biochem 2001, 299:63-70. OpenURL

  11. Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations.

    Genes Immun 2005, 6:279-284. OpenURL

  12. Mahoney D, Carey K, Fu M, Snow R, Cameron-Smith D, Gianni Parise G, Tarnopolsky MA: Real-time RT-PCR analysis of housekeeping genes in human skeletal muscle following acute exercise.

    Physiol Genomics 2004, 18:226-231. OpenURL

  13. Barnard G, Staniunas R, Bao S, Mafune K, Steele GD, John JL, Gollan JL, Chen LB: Increased expression of human ribosomal phosphoprotein P0 messenger RNA in hepatocellular carcinoma and colon carcinoma.

    Cancer Res 1992, 52:3067-3072. OpenURL

  14. Henry J, Coggin D, King C: High-level expression of the ribosomal protein L19 in human breast tumors that overexpress erbB-2.

    Cancer Res 1993, 53:1403-1408. OpenURL

  15. Suzuki T, Higgins P, Crawford D: Control selection for RNA quantitation.

    Biotechniques 2000, 29:332-337. OpenURL

  16. Murphy R, Watt K, Cameron-Smith D, Gibbons CJ, Snow RJ: Effects of creatine supplementation on housekeeping genes in human skeletal muscle using real-time RT- PCR.

    Physiol Genomics 2003, 12:163-174. OpenURL

  17. Ross D, Scherf U, Eisen M, Perou CM, Rees C, Spellman P, Iyer V, Jeffrey S, de Rijn MV, Waltham M, Pergamenschikov A, Jeffrey CF, Lee JC, Lashkari D, Shalon D, Myers TG, Weinstein JN, Botstein D, Brown PO: Systematic variation in gene expression patterns in human cancer cell lines.

    Nat Genet 2000, 24:227-235. OpenURL

  18. Schmidt U, Fuessel S, Koch R, Baretton GB, Lohse A, Tomasetti S, Unversucht S, Froehner M, Wirth MP, Meye A: Quantitative multi-gene expression profiling of primary prostate cancer.

    Prostate 2006, 66:1521-1534. OpenURL

  19. Khimani A, Mhashilkar A, Mikulskis A, O’Malley M, Liao J, Golenko EE, Mayer P, Chada S, Killian JB, Lott ST: Housekeeping genes in cancer: normalization of array data.

    Biotechniques 2005, 38:739-745. OpenURL

  20. Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits.

    Biotechnol 1999, 75:291-295. OpenURL

  21. Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, Pazzagli M, Bustin SA, Orlando C: Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies.

    Anal Biochem 2002, 309:293-300. OpenURL

  22. Tsuji N, Kamagata C, Furuya M, Kobayashi D, Yagihashi A, Morita T, Horita S, Watanabe N: Selection of an internal control gene for quantitation of mRNA in colonic tissues.

    Anticancer Res 2002, 22:4173-4178. OpenURL

  23. Vaarala M, Porvari KS, Kyllonen AP, Mustonen MV, Lukkarinen O, Vihko PT: Several genes encoding ribosomal proteins are over-expressed in prostate-cancer cell lines: confirmation of L7a and L37 over-expression in prostate-cancer tissue samples.

    Int J Cancer 1998, 78:27-32. OpenURL

  24. Viraj P, Melissa AH, Philippa H, Hillyer P, Navarro MB, Rabin RL: Systematic method for determining an ideal housekeeping gene for real-time PCR analysis.

    J Biomol Tech 2008, 19:342-347. OpenURL

  25. Akerman G, Rosenzweig B, Domon O, Tsai C, Bishop M, McGarrity LJ, James T: Alterations in gene expression profiles and the DNA-damage response in ionizing radiation-exposed TK6 cells.

    Environ Mol Mutagen 2005, 45:188-205. OpenURL

  26. Vandesompele J, De Preter K, Pattyn F, Poppe B, Roy NV, Paepe LJ, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.

    Genome Biol 2002., 3

    research0034.1-0034.11

    OpenURL

  27. Ghandhi SA, Yaghoubian B, Amundson SA: Global gene expression analyses of bystander and alpha particle irradiated normal human lung fibroblasts: Synchronous and differential responses.

    BMC Medical Genomics 2008, 1:63.

    doi

    BioMed Central Full Text OpenURL

  28. Combs S, Jakel O, Haberer T, Debus J: Particle therapy at the Heidelberg Ion Therapy Center (HIT) - Integrated research-driven university-hospital-based radiation oncology service in Heidelberg, Germany.

    Radiother Oncol 2010, 95:41-44. OpenURL