Overall quantitative results on significant differentially expressed genes obtained by the DNA microarray method, and analysis performed by SAM – Significance Analysis of Microarray (FDR < 5 %), for the comparison irradiated
Glioblastoma multiforme (GBM) is one of the most frequent tumors in the central nervous system and the most malignant tumor among gliomas. In the past two decades, cytogenetic and molecular genetic studies have identified a number of recurrent chromosomal abnormalities and genetic alterations in malignant gliomas, particularly in GBM . It was already described that GBM harbors combinations of the following genetic alterations: loss of heterozygozity of 10q,
The current GBM treatment involves aggressive management including surgery, adjuvant temozolomide-based chemotherapy, and radiotherapy , but GBM patients still present a dismal prognosis, and the median survival is 14.6 months from diagnosis . Although radiotherapy has been found to significantly prolong survival rates for GBM patients, radioresistance is a typical characteristic of this tumor .
Current genome-wide studies and the molecular characterization of GBM have allowed the identification of potential new targets, development of novel therapeutic small molecules and monoclonal antibodies and initiation of clinical trials with these targets [6, 8-10]. However, there is a wide molecular diversity and heterogeneity associated with the aberrantly GBM signaling pathways, culminating in the relative lack of success of these new approaches . Recently, an alternative strategy involves the selective targeting of GBM stem cells, which are resistant to chemo- and radiotherapy. But still, almost all small-molecule inhibitors designed to target these cells failed to demonstrate the effectiveness of this strategy, compared with the conventional therapy .
Considering that most of the treatment protocols are still ineffective, novel approaches are needed towards killing of GBM cells. Transcription machinery, as well as its regulatory elements is also a feasible new target for the application of molecular therapies. Transcription of DNA is dependent on the spatially and temporally coordinated interaction between transcriptional machinery involving RNA polymerase II, transcription factors (TFs)) and transcriptional regulatory components (promoter elements, enhancers, silencers and locus control regions) [12, 13]. The low level of transcription, directed by the general transcription factors associated to RNA polymerase core enzyme, is known as basal transcription . However, there is a rapidly expanding number of ‘context-dependent’ transcription factors that bind DNA and these TFs are capable of positively or negatively regulating the transcription process depending on the context of their binding sites, the complement of protein interactions and other environmental influences .
Postgenomic analyses of major transcription factor families, in both malignant and nonmalignant cell types, have opened new discussions about TF function. The mechanisms by which TFs act in cancer cell systems appear to exhibit a restricted repertoire of skills and plasticity displayed by normal cell systems . The evolution of a restricted malignant transcriptome can be seen clearly in the nuclear receptor superfamily, but is also apparent in the MYC and AP-1 networks . Oncogenic transcriptional rigidity reflects the simultaneous deregulation of target loci such that proliferative and survival signals are enhanced and antimitotic inputs are either limited or lost. Co-repressor proteins significantly contribute with the disruption of these processes . Therefore, understanding mechanisms involved in gene regulation and transcriptional network may lead to a better knowledge about the crucial functions of TFs, providing information to explore possibilities of their application as molecular targets in cancer therapy .
A valuable tool to study the transcription machinery is the DNA microarray technology , which measures the transcript expression of thousands of genes to identify changes in expression profiles at different biological conditions [20-24], thus allowing to compare different cell types under diverse treatment conditions. The influence of
Recently, information on the regulation of gene expression can also be used within the context of functional enrichment tests, and different databases containing TFs binding sites and other regulatory motifs are available, allowing to scan promoter regions of genes to detect the presence of target motifs . This information allow to determining whether a set of pre-selected genes is under control of TFs. FatiGO +  is a web-based tool capable of associating TFs that are common to a gene set used as parameters. This TF prediction method was already applied to a GBM dataset obtained from public repositories of microarray experiments, and the up-regulation of two predicted TFs, E2F1 and E2F4, was validated for several GBM cell lines , demonstrating the suitability of this method.
In the current study, we aimed to identify TFs that could be predicted from significant differentially expressed genes (previously obtained in microarray experiments in irradiated GBM cells) using an
We found few predicted TFs that were common between GBM cell lines, while several exclusive TFs were found for each cell line, indicating that the transcriptional response to ionizing radiation is very particular to each cell line examined in our microarray study, a fact that can be due to the genetic heterogeneity inherent to GBM cells. In spite of this, there was a convergence of biological functions among cell lines; the most relevant processes were related to apoptosis, cell proliferation, cell cycle, DNA repair, oxidative stress, among others. Furthermore, the present results also showed several TFs that were already reported as associated to cancer and stress responses.
2. Materials and methods
2.1. Briefly characterization of the experiment that provided the statistically modulated genes used for TF prediction
2.1.1. Cell culture and irradiation
Human GBM T98G and U87MG cell lines were supplied by the American Type Culture Collection (ATCC) (Rockville, Maryland, USA) and gently donated by Dr. Mari C. Sogayar (Universidade de São Paulo, Brazil). U343MG-a (U343), a cell line established from a primary malignant astrocytoma in an adult , was kindly donated by Dr. James T. Rutka (The Arthur and Sonia Labatt Brain Tumour Research Center, Canada); U251MG cell lines was also purchased from the ATCC (Rockville, MD, USA) and gently donated by Dr. Guido Lenz (Universidade Federal do Rio Grande do Sul, Brazil) . All cell lines grown in the presence of DMEM + HAM F10 medium (Sigma-Aldrich, St. Louis, USA) plus 10% fetal calf serum (Cultilab, Campinas, Brazil), and kept at 37°C and 5% CO2, until they reach semiconfluency. Cells were sub-cultured and 1x106 cells were seeded in 25 cm2 flasks, being incubated at 37°C for 48 h, and irradiated with 8 Gy of gamma-rays (60Co source, dose rate of 2.0 Gy / min., Unit Gammatron S-80, Siemens, 1.25 MeV, HC-FMRP/USP).
2.1.2. cDNA microarrays method and analysis
Two experiments with irradiated and sham-irradiated GBM cells were carried out using a glass slide microarrays containing ~4300 clones of cDNA probe (in replicates) from the human IMAGE Consortium cDNA library ; kindly provided by Dr. Catherine Nguyen (INSERM-CNRS, Marseille, France)], and prepared according to the protocol described by Hegde et al . Microarrays were spotted onto glass slides (Corning, Lowell, MA, USA) by using a Generation III Array Spotter (Amersham Molecular Dynamics, Sunnyvale, USA) according to the manufacturer’s instructions.
Total RNA extraction was performed for all cell lines, 30 min. and 6 h after irradiation, using the Trizol reagent (Invitrogen, Carlsbad, USA) according to manufacturer’s instructions. Each cDNA sample was spotted twice in the slide (duplicate spots). The cDNA complex probes were prepared using the CyScribe Post Labeling Kit (Amersham Biosciences, Buckinghamshire, UK) as previously described . Hybridizations were carried out using an automatic system (Automatic Slide Processor, Amersham Biosciences, UK) and signals were immediately captured after the final wash procedure, using a Generation III laser scanner (Amersham Biosciences, UK). This array platform was already used in several studies [22-26, 36].
2.1.3. Data acquisition and gene expression analysis
The provided microarray data was filtered and normalized [25, 36]. Following the normalization procedure, microarray data was exported to tab-delimited tables in MEV format and analyzed in MEV (v. 3.1) software .
The gene set submitted to SAM (Significance Analysis of Microarray ) were previously obtained by a t-test (α=5%) comparing irradiated (8 Gy)
|U343MG-a (30 min.)||7||116||+1.53 to -2.42|
|U343MG-a (6 h)||3||11||+1.83 to -1.39|
|U87MG (30 min.)||56||73||+1.88 to - 2.95|
|U87MG (6 h)||86||54||+1.68 to -1.95|
|T98G (30 min)||32||0||+2.26 to +1.13|
|T98G (6 h)||16||7||+2.70 to -1.63|
|U251MG (30 min.)||12||69||+1.85 to -1.40|
|U251MG (6 h)||17||20||+2.28 to -2.18|
The list of significantly modulated genes was obtained for a FDR < 5%. U343 cells showed 123 and 14 significantly differentially expressed genes at 30 min and 6 h after irradiation, respectively, whereas U87 showed 129 genes at 30 min and 140 genes at 6 h; T98G cell line displayed 32 and 23 significantly up-regulated genes at 30 min. and 6 h, respectively, whereas U251 showed 81 genes at 30 min. and 37 genes at 6 h (Table 2).
2.2. Transcriptional factor analysis
The analysis of TFs related to the significant differentially expressed genes (SAM) was performed by applying the FatiGO + . This program uses the TRANSFAC , and CisRed  transcription factors database, including their respective binding sites and regulated genes.
FatiGO + analyzes if the pre –selected set of genes (provided after SAM analysis), are under control of the same TF, and search for significant enrichments to each TF that is associated to the gene list compared to the complete reference list, containing ~4300 clones that were spotted onto the microarray slide .
The p-values obtained in the analysis of regulatory elements have been established by the program using the Fisher's exact test for multiple comparisons (unadjusted p-value). The Enrichment Index (EI) calculated for each TF corresponds to the increment obtained regarding the number of genes (%) statistically modulated (SAM) that are associated to a specific TF (List #1) divided by the total number of genes (%) in the array set that were predicted as targets for the same TF (List #2):
EI = % gene List #1/ % gene list #2
The TFs were selected according to unadjusted p-values < 0.05. The genes were submitted to FatiGO + v3.2, using the Gene symbol identifier and the selected gene distance of 10 kb. After selecting the TFs associated to modulated genes (SAM), a search was conducted in PubMed (http://www.ncbi.nlm. Nih.gov / sites / entrez /) looking for biological functions of those TFs.
2.3. Quantitative real-time PCR (qPCR)
We analyzed the transcript expression of HEB, a predict TF that was found associated to 57.7% of up-regulated genes in U87 cells, 30 min. after IR. The reverse transcription step was carried out in the remaining RNA samples from microarray experiments, with the Superscript III Reverse Transcriptase kit (Invitrogen, USA), according to manufacturer’s instructions. The integrity of cDNA samples was validated by the amplification of the endogenous B2M gene and visualization in agarose gel electrophoresis. qPCR was carried out using SYBR green master mix (Applied Biosystems, Foster City, USA) and the expression levels were estimated by the Relative Expression Software Tool (REST) , using 10000 interactions as setup parameter. All primers (Integrated DNA Technologies, Coralville, USA) were designed in Primer3 software  and are displayed on Table 2. The reactions were carried out in the Applied Biosystems 7500 Real-Time PCR System (Applied Biosystems, USA) equipment, using primer sets with an annealing temperature near 60°C and an amplicon of 100–120 bp. The PCR cycle was the following: pre-heating at 50°C for 2 min., 10 min. at 95°C (denaturation step), followed by 40 cycles at 95°C for 15 sec., and at 60°C for 60 sec. The dissociation curves were set up as following: 95°C for 15 sec., 60°C for 20 sec. and 95°C for 15 sec.
||5’- AGGCTATCCAGCGTACTCCA - 3’||112|
||5’ - TCAATGTCGGATGGATGAAA - 3’|
||5’ - CCGCTTGAGTTATCCTCCAC - 3’||116|
||5’ - GTGAGGCAGCAACGTAAGGT - 3’|
2.4. Western Blot (WB)
Protein extraction was performed with the Trizol reagent (Invitrogen, Carlsbad-USA) according to the manufacturer's instructions, using the same samples for RNA extraction. These samples were obtained from U87 cells collected at 30 min. post-irradiation. The expression of HEB was analyzed by Western blot, using ACTB as internal control. Samples were prepared with 30 μg of total protein. After electrophoresis, proteins were transferred from the gel to the membrane Invitrolon PVDF using the XCell IITM Blot Module system (Invitrogen, Carlsbad - USA). The immunodetection and protein visualization were conducted with the WesternBreeze Chromogenic kit (Invitrogen, Carlsbad - USA). The antibodies used in this study were anti-HEB (Santa Cruz, Santa Cruz, USA), and anti-ACTB (Cell Signaling, Danvers, USA), dilution of 1:1000.
We performed densitometric analysis of WB bands using the GelPro Analyzer (MediaCybernetics, Rockville, USA) 4.0, and the relative expression of HEB was calculated relatively to ACTB.
In the FatiGO + analysis, the lists of statistically modulated genes (SAM) were up-loaded in order to find TFs that were significantly associated with up-regulated and down-regulated genes for non-adjusted p-values < 0.05 (Table 3).
A Venn diagram was constructed based on the numbers of predicted TFs from data set previously obtained for each cell line (microarray experiments) (Fig. 1). TFs predicted for 30 min and 6 h were pooled together. Each cell line showed a number of exclusive TFs, but we also observed common TFs between cell lines. Out of 18 exclusive TFs found for U87MG cell line, PEBP (p = 0.008), Bach2 (p = 0.007), Freac-4 (p = 0.003), HLV (p = 0.006), Evi-1 (p = 0.009) displayed the lowest p-values, while PPARG and SEF-1 displayed the highest EI (31.3). U343 presented 9 exclusive TFs; High values of EI were found for MAF (33.1), E2F:DP-1 (22.0), PR (45.5) and STAT3 (38.5), and ARP-1 was the TF presenting the lowest p-value (0.009). T98G cells displayed only 6 exclusive TFs: EBF, Pax, Pbx1b, C/EBP, Poly A downstream element and Pax-9; two of them, EBF and PolyA showed low p-values, 0.005 and 0.007, respectively. Regarding U251MG cells, 13 TFs were predicted, and only APOLYA presented a high EI (27.8) (Table 3).
Interestingly, STAT3 was the common TF found for TP53 wild-type cells; however, this TF was associated with up-regulated genes in U87, and with down-regulated genes in U343. Only one TF (VBP) was common among three cell lines (U87, U343 and U251), being associated with down-regulated genes. Among the
Therefore, our results showed that most of the predicted TFs were exclusive to each cell line and few TFs were common among the GBM cell lines; these results indicate that the transcriptional response to ionizing radiation is very particular to each cell line, and most probably this can be due to the genetic heterogeneity of GBM cells.
By using the real time qPCR method, we confirmed the expression of HEB to validate the
We also look for HEB protein expression by Western Blot; both ACTB and HEB antibodies were used for irradiated and sham-irrradiated U87, 30 min after irradiation (Fig.2B). The relative expression values calculated by densitometric analysis showed that HEB expression was 1.7 higher in irradiated (8 Gy) cells, relatively to the control value (Fig. 2C).
|30 min. (↓)||ARP-1||14.8||6.7||2.2||0.009|
|Imperfect Hogness/Goldberg BOX||2.3||0.2||14.2||0.016|
|Muscle initiator sequence-20||20.5||12.2||1.7||0.031|
|6 h (↑)||MAF||50.0||1.5||33.1||0.031|
|6 h (↓)||ICSBP||40.0||3.7||10.9||0.013|
|30 min. (↑)||C/EBPgamma||85.0||64.4||1.3||0.007|
|30 min. (↓)||Bach2||63.5||43.9||1.4||0.007|
|6 h (↑)||DEC||45.6||31.9||1.4||0.032|
|6 h (↓)||Freac-4||10.5||1.6||6.8||0.003|
|30 min. (↑)||EBF||29.6||10.5||2.8||0.005|
|6 h (↑)||Pbx1b||14.3||1.8||8.2||0.026|
|6 h (↓)||C/EBP||100.0||20.0||5.0||0.002|
|Poly A downstream element||75.0||12.6||6.0||0.007|
|30 min. (↑)||SMAD-4||40.0||10.4||3.9||0.015|
|30 min. (↓)||VBP||21.7||10.4||2.1||0.025|
|6 h (↑)||MAF||15.4||1.5||10.2||0.017|
|6 h (↓)||HES1||45.5||15.6||2.9||0.019|
Recently, genome wide technologies, such as DNA microarrays, provide a huge amount of information about gene expression, but require additional bioinformatics analyses for data interpretation. In order to reduce complex signatures to a small number of activated transcriptional elements, new bioinformatics tools have been developed. To date, genome-wide TF-binding regions and sites were identified using a variety of indirect methods and data sets, revealing abundant binding sites for different TFs in mammalian cells [40-43]. Using lists of differentially expressed genes that were generated by microarray experiments, it is possible to predict TFs that can target common binding sites to a gene set. In the current study, we performed an
4.1. Commonly predicted TFs for two or three GBM cell lines
We found few TFs that were predicted for more than one cell line: MAF, TEF, ATF4, STAT3, VBP and C/EBPGamma. Most of these TFs (STAT3, TEF and VBP) are related to apoptosis, while other biological classes were also found, such as oxidative stress (ATF4), differentiation (MAF) and nucleotide excision repair (C/EBPGamma).
Thyrotroph embryonic factor (TEF) and human hepatic leukemia factor (HLF) are members of the PAR (
Activating transcription factor 4 (
MAF, the v-maf musculoaponeurotic fibrosarcoma oncogene homolog (avian) is a unique subclass of bZIP proteins and was predicted from up-regulated genes in U343 and U251 at 6 h. Depending on the binding site and binding partner, the encoded protein can be a transcriptional activator or repressor. Members of the MAF family appear to play important roles in the regulation of differentiation . MAF was found up-regulated in various cancers, such as colon cancer (but only in tumors that presented high levels of COX-2 expression) , a small subset of myelomas, hairy cell leukemia, T- and NK-cell neoplasms and small cell lymphomas .
C/EBPGamma, a member of the CCAAT/enhancer-binding protein (C/EBP) family of transcription factors was predicted from up-regulated genes in U87 (30 min.) and U251 (6 h). This TF regulates the expression of ERCC5 , and is a participant of DNA repair , particularly in the nucleotide excision repair .
The predicted TFs represent the overall GBM response to irradiation, since they were selected for more than one cell line, and as mentioned above, their functions are directly associated with stress responses involving apoptosis, DNA repair and ER stress. Moreover, as an example of STAT3, which is in clinical trial , predicted TFs may constitute potential targets to be investigated and validated in cancer treatment.
4.2. Exclusively predicted TFs for irradiated U343 cell line
Few predicted TFs (PR, E2F4:DP-1, and ARP-1) associated with statistically significant modulated genes in irradiated U343 cells were found exclusive to this cell line. The functions of these TFs are mainly associated with cell cycle/ tumor growth, being involved in various types of cancer. The overexpression of E2F4 and its binding partner DP-1 revealed a dual function of E2F4, which acts as an activator as well as a repressor, being implicated in positive regulation of the cell cycle . In a previous work, the up-regulation of E2F4 was confirmed for several GBM cell lines , demonstrating the potential of this TF as molecular target in cancer therapy.
Progesterone receptor (PR), a nuclear receptor transcription factor was associated with down-regulated genes in U343 cells (6 h). Steroid hormones participate in several physiological and pathological processes in the brain, including the regulation of tumor cell growth [62-64]. Progesterone exerts many of its effects by the interaction with specific intracellular receptors [62, 65].
ARP-1, also known as orphan nuclear receptor chicken ovalbumin upstream promoter transcription factor II (COUP-TFII) is a member of the steroid/thyroid nuclear hormone receptor superfamily  and can act as transcriptional repressor or activator. ARP-1 plays critical roles in organogenesis [67-70], and is a major angiogenesis regulator within the tumor microenvironment during pancreatic tumor progression and metastasis . Besides, ARP-1 was associated with therapy response in oligodendroglial tumors with 1p/19q loss . In the present study, this TF was predicted from down-regulated genes (30 min. and 6 h), indicating its possible involvement in radiation responses restricted to U343 cell line.
4.3. Exclusively predicted TFs for irradiated U87 cell line
A high number of TFs with low p-value (Bach2, PEBP, Freac-4, HLF and Evi-1) and high EI (PPARG and SEF-1) was associated with significant expressed genes in U87 cell line.
PPARG is a member of the peroxisome proliferator activated receptor (PPAR) family, a subfamily of the nuclear receptor superfamily . The protein level of this receptor has been recently identified as a significant prognostic marker . Interestingly, recent studies have shown that PPARG is expressed in normal and malignant human brain, and the treatment with PPARG agonists induces growth arrest and apoptosis in brain tumor cells
Basic leucine zipper transcription factor 2 (Bach2) is an evolutionarily related member of the BTB-basic region leucine zipper transcription factor family. Bach2 can function as transcriptional activator and repressor . This TF down-regulates cell proliferation of the neuroblastoma cell line N1E-115 and negatively affects their potential to differentiate, being considered as gatekeeper of the differentiated status . Bach2 presents high frequency of loss of heterozygosity of the Bach2 gene in human B cell lymphomas . Consistent with its putative role as a tumor suppressor, Bach2 enhances apoptosis in response to oxidative stress [84, 85].
The transcription Factor PEBP, also called Raf kinase inhibitor protein (RKIP) is a member of the phosphatidylethanolamine-binding protein (PEBP) family. RKIP plays a pivotal modulatory role in several protein kinase signaling cascades. RKIP regulates the activity of the Raf/MEK/ERK, which is responsible for proliferation and differentiation of diverse cell types . It has been reported that RKIP was poorly expressed in primary tumors, being absent in various metastatic cancers; its induction sensitize resistant tumor cells to apoptosis by various chemo- and immunotherapeutic drugs, as well as inhibitors of metastasis .The absence of RKIP is also associated with highly malignant behavior and poor survival of patients .
The forkhead domain is a monomeric DNA binding motif that defines a rapidly growing family of eukaryotic transcriptional regulators. We found Freac-4, also known as Forkhead Box D1 (FOXD1) associated with down-regulated genes in U87 cells (6 h). This gene was found repressed in chemoresistant tumors, as analyzed by microarrays . However, FOXD1 and FOXD2 were highly expressed in prostate cancer and lymph node metastases, among various cancer types . In another study, using kidney-derived cell lines, it was suggested that FOXD1 may be regulated by TP53, WTAR (a mutated form of WT1) and WT1 .
As already mentioned, HLF is a member of PAR bZIP transcription factors, and was associated with down-regulated genes, 6 h after irradiation in U87. PAR bZIP proteins are also involved in apoptosis induction .The fused gene
Activator protein one (AP1) transcription factors are a family of jun and fos proteins, whose subunits present diverse pro/anti-cancer effects, like inhibition or increase in proliferation, inhibition of apoptosis and angiogenesis [93, 94]. AP-1 is one of the genes early activated after radiation in primary human B cells . The inhibition of AP-1 blocks the proliferation of breast tumor cells by suppressing the growth factor signaling . The modulation of AP-1 activity may be a new attempt to reduce the malignant transformation. However, only the function involved with malignancy should be targeted , since AP-1 presents oncogenic and anti-oncogenic properties. This TF was associated to U87 cells (up-regulated genes, 30 min.) and (down-regulated genes, 6 h).
The EVI1 gene encodes a zinc finger transcription factor with important roles in normal development and leukemogenesis. Reports in animal model and findings in
U87 cell line presented several predicted TFs with significant p values and higher enrichment index than other cell lines. Most of the predicted TFs are related to apoptosis (PPARG, Batch2, PEBP, HLF, AP1 and EVI1), but they were associated with up- or down-regulated genes. Overall, the biological functions of these TFs are related to cell proliferation, differentiation and tumor growth, indicating the relevance of their deregulation in cancer development and malignancy, and possibly, in tumor responses to anti-cancer therapies.
4.4. Exclusively predicted TFs for irradiated T98G cell line
Only two TFs were predicted for T98G cells: EBF and C/EBP. The early B-cell factors (EBF) are a family of four highly conserved DNA-binding transcription factors with an atypical zinc-finger and helix-loop-helix motif. Zardo and colleagues found that the EBF3 locus on the human chromosome 10q is deleted or methylated in brain tumors . Functional studies revealed that EBF3 activates genes involved in cell cycle arrest and apoptosis, while in opposite, it represses genes involved in cell survival and proliferation . Therefore, EBFs represent a novel tumor suppressor whose inactivation blocks normal development and contributes to tumorigenesis of diverse types of human cancer .
CEBP is also known as basic leucine zipper transcription factor, CCAAT/enhancer binding protein alpha (CEBPA), which directly interacts with CDK2 and CDK4 and arrests cell proliferation by inhibition of these kinases . CEBPA is crucial for normal granulopoiesis, and dominant-negative mutations of CEBPA gene were found in patients with myeloblastic subtypes (M1 and M2) of acute myeloid leukemia . CEBPA also plays a role in DNA damage response dependent on TP53, as observed in keratinocytes . C/EBPA was found silenced in human squamous cell carcinoma (SCC), and loss of C/EBPA confers susceptibility to UVB-induced skin SCCs involving defective cell cycle arrest in response to UVB . Interestingly, these findings indicate the role of CEBP in DNA damage responses, and possibly, the potential of this TF to be explored as therapeutic molecular target.
4.5. Validation of TF prediction
As a predicted TF associated with up-regulated genes, HEB (p-value = 0.021 and EI = 1.5) was chosen to be studied in terms of expression levels, aiming to validate the
HEB is a member of the class A basic helix-loop-helix (bHLH) family that participates in the nervous system development [107, 108]. According to O'Neil et al. , the repression of E47/HEB has been associated with the induction of leukemia in mice. In another study, it was demonstrated the induction of HEB in gliomas compared with non-neoplastic brain tissue . Moreover, HEB seems to be involved in cell proliferation control of neural stem cells and also progenitor cells, being important to sustain their undifferentiated state during embryonic and adult neurogenesis . Although HEB expression has not yet been correlated with radiation responses in GBM cells, in the present study, we found its association with significant differentially expressed genes at 30 min. following irradiation in U87 cell line; interestingly, we also showed that HEB transcript and protein expression was induced in irradiated U87 cell line, 30 min. after irradiation. This finding, although restrict to one TF, indicates the validity of the TF prediction by
5. Conclusions: Lessons from TF prediction in irradiated GBM cells
The present findings about prediction of TFs associated to differentially expressed genes in GBM cell lines showed that few TFs were shared among different GBM cell lines, while several TFs were found exclusive to each cell line, indicating that the transcriptional response to ionizing radiation is very particular to each cell line, probably due to the genetic heterogeneity, which is characteristic of GBM cell lines. In spite of this observation, several biological functions were similar among cell lines, such as apoptosis, cell proliferation, cell cycle control, DNA repair, ER stress, and differentiation. Furthermore, most of the predicted TFs were already reported as differentially expressed, deleted or mutated in cancer, including GBM. However, apart the similarity of biological functions, different pathways seems to be associated to the predicted TFs. Interestingly, we could not find TP53 as a TF associated to the data set (List #1) analyzed in the present study, even for the GBM cells that were wild-type for
The most intriguing finding refers to apoptotic related TFs. Probably, predicted TFs related to apoptosis control, and found associated with expressed genes at early time (30 minutes) following irradiation, are related to survival in GBM cells; this is supported by reports showing that in general, these cells are very resistant to undergo apoptosis, even under conditions of drug treatment or radiation exposure [112-114]. In fact, GBM cells seem to be capable of activating several pathways to escaping from cell killing by anticancer therapies.
Even considering the relevance of our findings, some methodological limitations should be mentioned regarding
In spite of the limitation mentioned above, in a previous study, we validated the expression of E2F  and the HEB expression was confirmed in the present study, both of them in GBM cell lines. In addition, we selected predicted TFs that were associated with stress response genes, and importantly, the TFs were reported as deregulated or mutated in different cancer types, thus indicating the relevance of further studies to better exploring the role of TFs in the context of therapeutic strategies based on molecular targets.
Chen, J, RM McKay, and LF Parada, Malignant glioma: lessons from genomics, mouse models, and stem cells. Cell, 2012. 149(1): p. 36-47.
Ohgaki, H, Genetic pathways to glioblastomas. Neuropathology, 2005. 25(1): p. 1-7.
McLendon, R et al., Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature, 2008. 455(7216): p. 1061-8.
Parsons, DW, et al., An integrated genomic analysis of human glioblastoma multiforme. Science, 2008. 321(5897): p. 1807-12.
Stupp, R, et al., Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med, 2005. 352(10): p. 987-96.
Ohka, F, A Natsume, and T Wakabayashi, Current trends in targeted therapies for glioblastoma multiforme. Neurol Res Int, 2012. 2012: p. 878425.
Stupp, R, et al., Chemoradiotherapy in malignant glioma: standard of care and future directions. J Clin Oncol, 2007. 25(26): p. 4127-36.
Sathornsumetee, S, et al., Molecularly targeted therapy for malignant glioma. Cancer, 2007. 110(1): p. 13-24.
Perry, J, et al., Novel therapies in glioblastoma. Neurol Res Int, 2012. 2012: p. 428565.
Polivka, J, Jr., et al., New molecularly targeted therapies for glioblastoma multiforme. Anticancer Res, 2012. 32(7): p. 2935-46.
Natsume, A, et al., Glioma-initiating cells and molecular pathology: implications for therapy. Brain Tumor Pathol, 2011. 28(1): p. 1-12.
Thomas, MC and CM Chiang, The general transcription machinery and general cofactors. Crit Rev Biochem Mol Biol, 2006. 41(3): p. 105-78.
Maston, GA, SK Evans, and MR Green, Transcriptional regulatory elements in the human genome. Annu Rev Genomics Hum Genet, 2006. 7: p. 29-59.
Hampsey, M, Molecular genetics of the RNA polymerase II general transcriptional machinery. Microbiol Mol Biol Rev, 1998. 62(2): p. 465-503.
Gaston, K and PS Jayaraman, Transcriptional repression in eukaryotes: repressors and repression mechanisms. Cell Mol Life Sci, 2003. 60(4): p. 721-41.
Battaglia, S, O Maguire, and MJ Campbell, Transcription factor co-repressors in cancer biology: roles and targeting. Int J Cancer, 2010. 126(11): p. 2511-9.
Thorne, JL, MJ Campbell, and BM Turner, Transcription factors, chromatin and cancer. Int J Biochem Cell Biol, 2009. 41(1): p. 164-75.
Mees, C, J Nemunaitis, and N Senzer, Transcription factors: their potential as targets for an individualized therapeutic approach to cancer. Cancer Gene Ther, 2009. 16(2): p. 103-12.
Kel, A, et al., Beyond microarrays: find key transcription factors controlling signal transduction pathways. BMC Bioinformatics, 2006. 7 Suppl 2: p. S13.
Tusher, VG, R Tibshirani, and G Chu, Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A, 2001. 98(9): p. 5116-21.
Sakamoto-Hojo, ET, et al., Gene expression profiles in human cells submitted to genotoxic stress. Mutat Res, 2003. 544(2 3): p. 403-13.
Fachin, AL, et al., Gene expression profiles in human lymphocytes irradiated in vitro with low doses of gamma rays. Radiat Res, 2007. 168(6): p. 650-65.
Fachin, AL, et al., Gene expression profiles in radiation workers occupationally exposed to ionizing radiation. J Radiat Res (Tokyo), 2009. 50(1): p. 61-71.
Carminati, PO, et al., Alterations in gene expression profiles correlated with cisplatin cytotoxicity in the glioma U343 cell line. Genet Mol Biol, 2010. 33(1): p. 159-68.
Godoy, PRDV, et al. Portrait of transcriptional expression profiles displayed by different glioblastoma cell lines. in: M Garami. (ed.) Molecular Targets of CNS Tumors. Rijeka: InTech; 2011. p. 277-300.
da Silva, GN, et al., Expression of genes related to apoptosis, cell cycle and signaling pathways are independent of TP53 status in urinary bladder cancer cells. Mol Biol Rep, 2011. 38(6): p. 4159-70.
Otomo, T, et al., Microarray analysis of temporal gene responses to ionizing radiation in two glioblastoma cell lines: up-regulation of DNA repair genes. J Radiat Res (Tokyo), 2004. 45(1): p. 53-60.
Godoy, PRDV, et al., Glioblastoma cell lines differing in TP53 status (mutant or wild-type) display different molecular signaling responses to gamma-irradiation Submitted, 2013.
Al-Shahrour, F, et al., FatiGO +: a functional profiling tool for genomic data. Integration of functional annotation, regulatory motifs and interaction data with microarray experiments. Nucleic Acids Res, 2007. 35(Web Server issue): p. W91-6.
Al-Shahrour, F, et al., BABELOMICS: a suite of web tools for functional annotation and analysis of groups of genes in high-throughput experiments. Nucleic Acids Res, 2005. 33(Web Server issue): p. W460-4.
Donaires, FS, et al., Transcriptional factors associated to over-expressed genes in glioblastoma multiforme, as predicted by in silico analysis on publicly microarray dataset Submitted, 2013.
Dirks, PB, et al., Retinoic acid and the cyclin dependent kinase inhibitors synergistically alter proliferation and morphology of U343 astrocytoma cells. Oncogene, 1997. 15(17): p. 2037-48.
Filippi-Chiela, EC, et al., Autophagy interplay with apoptosis and cell cycle regulation in the growth inhibiting effect of resveratrol in glioma cells. PLoS ONE, 2011. 6(6): p. e20849.
Lennon, G, et al., The I.M.A.G.E. Consortium: An Integrated Molecular Analysis of Genomes and Their Expression. Genomics, 1996. 33(1): p. 151-152.
Hegde, P, et al., A concise guide to cDNA microarray analysis. Biotechniques, 2000. 29(3): p. 548-50, 552-4, 556 passim.
Mello, SS, et al., Delayed effects of exposure to a moderate radiation dose on transcription profiles in human primary fibroblasts. Environ Mol Mutagen, 2011. 52(2): p. 117-29.
Saeed, AI, et al., TM4: a free, open-source system for microarray data management and analysis. Biotechniques, 2003. 34(2): p. 374-8.
Wingender, E, et al., TRANSFAC: an integrated system for gene expression regulation. Nucleic Acids Res, 2000. 28(1): p. 316-9.
Robertson, G, et al., cisRED: a database system for genome-scale computational discovery of regulatory elements. Nucleic Acids Res, 2006. 34(Database issue): p. D68-73.
Ferretti, V, et al., PReMod: a database of genome-wide mammalian cis-regulatory module predictions. Nucleic Acids Res, 2007. 35(Database issue): p. D122-6.
Bieda, M, et al., Unbiased location analysis of E2F1-binding sites suggests a widespread role for E2F1 in the human genome. Genome Res, 2006. 16(5): p. 595-605.
Birney, E, et al., Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature, 2007. 447(7146): p. 799-816.
Cawley, S, et al., Unbiased mapping of transcription factor binding sites along human chromosomes 21 and 22 points to widespread regulation of noncoding RNAs. Cell, 2004. 116(4): p. 499-509.
Bowman, T, et al., Stat3-mediated Myc expression is required for Src transformation and PDGF-induced mitogenesis. Proc Natl Acad Sci U S A, 2001. 98(13): p. 7319-24.
Brennan, C, et al., Glioblastoma subclasses can be defined by activity among signal transduction pathways and associated genomic alterations. PLoS ONE, 2009. 4(11): p. e7752.
Lo, HW, et al., Constitutively activated STAT3 frequently coexpresses with epidermal growth factor receptor in high-grade gliomas and targeting STAT3 sensitizes them to Iressa and alkylators. Clin Cancer Res, 2008. 14(19): p. 6042-54.
Weissenberger, J, et al., IL-6 is required for glioma development in a mouse model. Oncogene, 2004. 23(19): p. 3308-16.
Doucette, TA, et al., Signal transducer and activator of transcription 3 promotes angiogenesis and drives malignant progression in glioma. Neuro Oncol, 2012.
Sen, M, et al., First-in-Human Trial of a STAT3 Decoy Oligonucleotide in Head and Neck Tumors: Implications for Cancer Therapy. Cancer Discov, 2012. 2(8): p. 694-705.
Gachon, F, et al., The loss of circadian PAR bZip transcription factors results in epilepsy. Genes Dev, 2004. 18(12): p. 1397-412.
Benito, A, et al., A novel role for proline- and acid-rich basic region leucine zipper (PAR bZIP) proteins in the transcriptional regulation of a BH3-only proapoptotic gene. J Biol Chem, 2006. 281(50): p. 38351-7.
De Angelis, R, et al., Functional interaction of the subunit 3 of RNA polymerase II (RPB3) with transcription factor-4 (ATF4). FEBS Lett, 2003. 547(1-3): p. 15-9.
Tian, X, et al., Modulation of CCAAT/enhancer binding protein homologous protein (CHOP)-dependent DR5 expression by nelfinavir sensitizes glioblastoma multiforme cells to tumor necrosis factor-related apoptosis-inducing ligand (TRAIL). J Biol Chem, 2011. 286(33): p. 29408-16.
Blais, JD, et al., Activating transcription factor 4 is translationally regulated by hypoxic stress. Mol Cell Biol, 2004. 24(17): p. 7469-82.
Blank, V and NC Andrews, The Maf transcription factors: regulators of differentiation. Trends Biochem Sci, 1997. 22(11): p. 437-41.
Asting, AG, et al., COX-2 gene expression in colon cancer tissue related to regulating factors and promoter methylation status. BMC Cancer, 2011. 11: p. 238.
Natkunam, Y, et al., Characterization of c-Maf transcription factor in normal and neoplastic hematolymphoid tissue and its relevance in plasma cell neoplasia. Am J Clin Pathol, 2009. 132(3): p. 361-71.
Crawford, EL, et al., CEBPG regulates ERCC5/XPG expression in human bronchial epithelial cells and this regulation is modified by E2F1/YY1 interactions. Carcinogenesis, 2007. 28(12): p. 2552-9.
Mullins, DN, et al., CEBPG transcription factor correlates with antioxidant and DNA repair genes in normal bronchial epithelial cells but not in individuals with bronchogenic carcinoma. BMC Cancer, 2005. 5: p. 141.
O'Donovan, A, et al., Isolation of active recombinant XPG protein, a human DNA repair endonuclease. J Biol Chem, 1994. 269(23): p. 15965-8.
Lee, BK, AA Bhinge, and VR Iyer, Wide-ranging functions of E2F4 in transcriptional activation and repression revealed by genome-wide analysis. Nucleic Acids Res, 2011. 39(9): p. 3558-73.
Graham, JD and CL Clarke, Physiological action of progesterone in target tissues. Endocr Rev, 1997. 18(4): p. 502-19.
Grunberg, SM, et al., Treatment of unresectable meningiomas with the antiprogesterone agent mifepristone. J Neurosurg, 1991. 74(6): p. 861-6.
Schrell, UM, et al., Hormonal dependency of cerebral meningiomas. Part 1: Female sex steroid receptors and their significance as specific markers for adjuvant medical therapy. J Neurosurg, 1990. 73(5): p. 743-9.
Mahesh, VB, DW Brann, and LB Hendry, Diverse modes of action of progesterone and its metabolites. J Steroid Biochem Mol Biol, 1996. 56(1-6 Spec No): p. 209-19.
Tsai, SY and MJ Tsai, Chick ovalbumin upstream promoter-transcription factors (COUP-TFs): coming of age. Endocr Rev, 1997. 18(2): p. 229-40.
Lee, CT, et al., The nuclear orphan receptor COUP-TFII is required for limb and skeletal muscle development. Mol Cell Biol, 2004. 24(24): p. 10835-43.
Takamoto, N, et al., COUP-TFII is essential for radial and anteroposterior patterning of the stomach. Development, 2005. 132(9): p. 2179-89.
You, LR, et al., Suppression of Notch signalling by the COUP-TFII transcription factor regulates vein identity. Nature, 2005. 435(7038): p. 98-104.
You, LR, et al., Mouse lacking COUP-TFII as an animal model of Bochdalek-type congenital diaphragmatic hernia. Proc Natl Acad Sci U S A, 2005. 102(45): p. 16351-6.
Qin, J, et al., COUP-TFII regulates tumor growth and metastasis by modulating tumor angiogenesis. Proc Natl Acad Sci U S A, 2010. 107(8): p. 3687-92.
Shaw, EJ, et al., Gene expression in oligodendroglial tumors. Cell Oncol (Dordr), 2011. 34(4): p. 355-67.
Issemann, I and S Green, Activation of a member of the steroid hormone receptor superfamily by peroxisome proliferators. Nature, 1990. 347(6294): p. 645-50.
Robbins, GT and D Nie, PPAR gamma, bioactive lipids, and cancer progression. Front Biosci, 2012. 17: p. 1816-34.
Cellai, I, et al., Antineoplastic effects of rosiglitazone and PPARgamma transactivation in neuroblastoma cells. Br J Cancer, 2006. 95(7): p. 879-88.
Grommes, C, et al., Inhibition of in vivo glioma growth and invasion by peroxisome proliferator-activated receptor gamma agonist treatment. Mol Pharmacol, 2006. 70(5): p. 1524-33.
Strakova, N, et al., Peroxisome proliferator-activated receptors (PPAR) agonists affect cell viability, apoptosis and expression of cell cycle related proteins in cell lines of glial brain tumors. Neoplasma, 2005. 52(2): p. 126-36.
Chearwae, W and JJ Bright, PPARgamma agonists inhibit growth and expansion of CD133+ brain tumour stem cells. Br J Cancer, 2008. 99(12): p. 2044-53.
Pestereva, E, S Kanakasabai, and JJ Bright, PPARgamma agonists regulate the expression of stemness and differentiation genes in brain tumour stem cells. Br J Cancer, 2012. 106(10): p. 1702-12.
Veliceasa, D, FT Schulze-Hoepfner, and OV Volpert, PPARgamma and Agonists against Cancer: Rational Design of Complementation Treatments. PPAR Res, 2008. 2008: p. 945275.
Oyake, T, et al., Bach proteins belong to a novel family of BTB-basic leucine zipper transcription factors that interact with MafK and regulate transcription through the NF-E2 site. Mol Cell Biol, 1996. 16(11): p. 6083-95.
Shim, KS, et al., Bach2 is involved in neuronal differentiation of N1E-115 neuroblastoma cells. Exp Cell Res, 2006. 312(12): p. 2264-78.
Sasaki, S, et al., Cloning and expression of human B cell-specific transcription factor BACH2 mapped to chromosome 6q15. Oncogene, 2000. 19(33): p. 3739-49.
Hoshino, H, et al., Oxidative stress abolishes leptomycin B-sensitive nuclear export of transcription repressor Bach2 that counteracts activation of Maf recognition element. J Biol Chem, 2000. 275(20): p. 15370-6.
Muto, A, et al., Identification of Bach2 as a B-cell-specific partner for small maf proteins that negatively regulate the immunoglobulin heavy chain gene 3' enhancer. Embo J, 1998. 17(19): p. 5734-43.
Yeung, K, et al., Suppression of Raf-1 kinase activity and MAP kinase signalling by RKIP. Nature, 1999. 401(6749): p. 173-7.
Wu, K and B Bonavida, The activated NF-kappaB-Snail-RKIP circuitry in cancer regulates both the metastatic cascade and resistance to apoptosis by cytotoxic drugs. Crit Rev Immunol, 2009. 29(3): p. 241-54.
Martinho, O, et al., Downregulation of RKIP is associated with poor outcome and malignant progression in gliomas. PLoS ONE, 2012. 7(1): p. e30769.
Ju, W, et al., Identification of genes with differential expression in chemoresistant epithelial ovarian cancer using high-density oligonucleotide microarrays. Oncol Res, 2009. 18(2-3): p. 47-56.
van der Heul-Nieuwenhuijsen, L, NF Dits, and G Jenster, Gene expression of forkhead transcription factors in the normal and diseased human prostate. BJU Int, 2009. 103(11): p. 1574-80.
Pierrou, S, et al., Cloning and characterization of seven human forkhead proteins: binding site specificity and DNA bending. Embo J, 1994. 13(20): p. 5002-12.
Smith, KS, et al., Disrupted differentiation and oncogenic transformation of lymphoid progenitors in E2A-HLF transgenic mice. Mol Cell Biol, 1999. 19(6): p. 4443-51.
Shaulian, E and M Karin, AP 1 in cell proliferation and survival. Oncogene, 2001. 20(19): p. 2390-400.
Schreiber, M, et al., Control of cell cycle progression by c-Jun is p53 dependent. Genes Dev, 1999. 13(5): p. 607-19.
Wilson, RE, et al., Early response gene signalling cascades activated by ionising radiation in primary human B cells. Oncogene, 1993. 8(12): p. 3229-37.
Liu, Y, et al., Inhibition of AP-1 transcription factor causes blockade of multiple signal transduction pathways and inhibits breast cancer growth. Oncogene, 2002. 21(50): p. 7680-9.
Libermann, TA and LF Zerbini, Targeting transcription factors for cancer gene therapy. Curr Gene Ther, 2006. 6(1): p. 17-33.
Wieser, R, The oncogene and developmental regulator EVI1: expression, biochemical properties, and biological functions. Gene, 2007. 396(2): p. 346-57.
Koos, B, et al., The transcription factor evi-1 is overexpressed, promotes proliferation, and is prognostically unfavorable in infratentorial ependymomas. Clin Cancer Res, 2011. 17(11): p. 3631-7.
Zardo, G, et al., Integrated genomic and epigenomic analyses pinpoint biallelic gene inactivation in tumors. Nat Genet, 2002. 32(3): p. 453-8.
Zhao, LY, et al., An EBF3-mediated transcriptional program that induces cell cycle arrest and apoptosis. Cancer Res, 2006. 66(19): p. 9445-52.
Liao, D, Emerging roles of the EBF family of transcription factors in tumor suppression. Mol Cancer Res, 2009. 7(12): p. 1893-901.
Wang, H, et al., C/EBPalpha arrests cell proliferation through direct inhibition of Cdk2 and Cdk4. Mol Cell, 2001. 8(4): p. 817-28.
Pabst, T, et al., Dominant-negative mutations of CEBPA, encoding CCAAT/enhancer binding protein-alpha (C/EBPalpha), in acute myeloid leukemia. Nat Genet, 2001. 27(3): p. 263-70.
Yoon, K and RC Smart, C/EBPalpha is a DNA damage-inducible p53-regulated mediator of the G1 checkpoint in keratinocytes. Mol Cell Biol, 2004. 24(24): p. 10650-60.
Thompson, EA, et al., C/EBPalpha expression is downregulated in human nonmelanoma skin cancers and inactivation of C/EBPalpha confers susceptibility to UVB-induced skin squamous cell carcinomas. J Invest Dermatol, 2011. 131(6): p. 1339-46.
Ravanpay, AC and JM Olson, E protein dosage influences brain development more than family member identity. J Neurosci Res, 2008. 86(7): p. 1472-81.
Uittenbogaard, M and A Chiaramello, Expression of the bHLH transcription factor Tcf12 (ME1) gene is linked to the expansion of precursor cell populations during neurogenesis. Brain Res Gene Expr Patterns, 2002. 1(2): p. 115-21.
O'Neil, J, et al., TAL1/SCL induces leukemia by inhibiting the transcriptional activity of E47/HEB. Cancer Cell, 2004. 5(6): p. 587-96.
Riemenschneider, MJ, TH Koy, and G Reifenberger, Expression of oligodendrocyte lineage genes in oligodendroglial and astrocytic gliomas. Acta Neuropathol, 2004. 107(3): p. 277-82.
Shu, HK, et al., The intrinsic radioresistance of glioblastoma-derived cell lines is associated with a failure of p53 to induce p21(BAX) expression. Proc Natl Acad Sci U S A, 1998. 95(24): p. 14453-8.
Rubel, A, et al., The membrane targeted apoptosis modulators erucylphosphocholine and erucylphosphohomocholine increase the radiation response of human glioblastoma cell lines in vitro. Radiat Oncol, 2006. 1: p. 6.
Honda, N, et al., Radiosensitization by overexpression of the nonphosphorylation form of IkappaB-alpha in human glioma cells. J Radiat Res, 2002. 43(3): p. 283-92.
Yamagishi, N, J Miyakoshi, and H Takebe, Enhanced radiosensitivity by inhibition of nuclear factor kappa B activation in human malignant glioma cells. Int J Radiat Biol, 1997. 72(2): p. 157-62.
Struhl, K, Gene regulation. A paradigm for precision. Science, 2001. 293(5532): p. 1054-5.
Chen, X, et al., Integration of external signaling pathways with the core transcriptional network in embryonic stem cells