Journal of Biomedicine and Biotechnology

Volume 2012 (2012), Article ID 654291, 8 pages

doi:10.1155/2012/654291

Research Article

Analyzing Gene Expression Profile in K562 Cells Exposed to Sodium Valproate Using Microarray Combined with the Connectivity Map Database

Xiang-Zhong Zhang,¹ Ai-Hua Yin,2,3 Dong-Jun Lin,¹ Xiao-Yu Zhu,4,5 Qian Ding,4,6 Chun-Huai Wang,⁴ and Yun-Xian Chen⁴

¹Department of Hematology, The Third Affiliated Hospital of Sun Yat-sen University, Guangdong 510630, Guangzhou, China

²Prenatal Diagnosis Center, The Women and Children Hospital of Guangdong Province, Guangdong 510010, Guangzhou, China

³Key Laboratory of Department of Health of Guangdong Province for Maternal and Children Metabolic-Genetic Disease Diagnosis, Guangdong 510010, Guangzhou, China

⁴Department of Hematology, First Affiliated Hospital of Sun Yat-sen University, Guangdong 510080, Guangzhou, China

⁵Department of Hematology, Anhui Provincial Hospital Affiliated to Anhui Medical University, Anhui 230001, Hefei, China

⁶Department of Hematology, Guizhou Provincial Hospital, Guizhou 550002, Guiyang, China

Received 6 February 2012; Revised 26 March 2012; Accepted 9 April 2012

Academic Editor: Eric W. Lam

Copyright © 2012 Xiang-Zhong Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

To explore the mechanism underlying antileukaemia effect of sodium valproate, the growth and survival of the K562 cell line were investigated. Global profiles of gene expression in K562 cells exposed to sodium valproate were assessed and validated. The differentially expressed genes identified were further used to query the connectivity map database to retrieve a ranked list of compounds that act on the same intracellular targets as sodium valproate. A significant increase in cell apoptosis and a change in gene expression profile were observed in valproate-exposed K562 cells. The significant enrichment analysis of gene ontology terms for the differentially expressed genes showed that these genes were involved in many important biological processes. Eight differentially expressed genes involved in apoptosis were verified by quantitative real-time PCR. The connectivity map analysis showed gene expression profile in K562 cells exposed to sodium valproate was most similar to that of HDACi and PI3K inhibitors, suggesting that sodium valproate might exert antileukaemic action by inhibiting HDAC as well as inhibiting PI3K pathway. In conclusion, our data might provide clues to elucidate the molecular and therapeutic potential of VPA in leukaemia treatment, and the connectivity map is a useful tool for exploring the molecular mechanism of drug action.

1. Introduction

Sodium valproate (VPA) is a well-known antiepileptic drug, also used to treat bipolar disorder, migraine, and neuropathic pain. Recently, VPA has been described as an HDAC inhibitor, resulting in an increased interest for its use in cancer therapy. Studies show that VPA, alone or in combination with other antileukaemic compounds, possesses significant antileukaemic actions on acute myeloid leukemia (AML) cells [1], chronic myeloid leukaemia (CML) cells [2–4], acute lymphoblastic leukemia (ALL) cells [5], and chronic lymphocytic leukemia (CLL) cells [6, 7]. Clinical trials show that VPA therapy is of clinical benefit in patients with acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS) [8–11]. However, the exact molecular mechanisms of VPA action on leukemia treatment remain poorly understood.

Microarray is a high-throughput tool which allows for the analysis of global gene expression profile in a single experiment and has been widely used for exploring molecular mechanisms of pathogenesis and drug treatment. Based on this high-throughput technology, the molecular mechanisms underlying the observed antileukaemic activity of VPA in CLL cells [12] and AML cells [13, 14] have been described. So far, however, there has been no study exploring genomewide gene expression changes in CML.

The Connectivity Map (CMAP) is a collection of genomewide transcriptional expression data from cultured human cells treated with bioactive small molecules. It connects drugs, genes, and diseases together through the transitory feature of common gene-expression changes. By comparing gene-expression signatures, this tool can be used to find connections among small molecules drugs affecting common molecular pathways and putative mechanisms of action of unknown drugs. CMAP has previously been used to discover the mechanisms of drug action [15–18] and disease pathogenesis [19]. In the current study, we first investigated effects of VPA on apoptosis and gene expression profiles of K562 cells, a model for CML, and then mined the CMAP database to explore the molecular mechanism underlying the observed anti-CML effect of VPA.

2. Methods and Materials

2.1. Culture of K562 Cell Line

K562 cells, human chronic myelogenic leukemic cell line, were procured from Sun Yat-sen University Cancer Center. The cells were grown in RPMI (Life Technology Corporation, Camarillo, CA, USA), supplemented with 15% fetal bovine serum, 100 U/mL penicillin, and 1 mg/mL streptomycin. Cultures were incubated at 37°C in 5% CO₂.

2.2. Assessment of Apoptosis by Annexin V/PI Dual Staining Method

In order to determine the effect of VPA treatment on apoptosis rate of K562 cells, FACS analysis was carried out as previously described [20].

2.3. Analysis of Gene Expression Profile

To analyze gene expression profile, K562 cells were first cultured for 12 hours with 2 mM VPA or without VPA as a control. Then, the K562 cells were harvested and total RNA was extracted using the RNeasy Mini Kit (QIAGen, USA) according to the manufacturer’s instructions. Total RNA quality was evaluated using formaldehyde agarose gel electrophoresis and was quantified via spectrophotometry (Nanodrop, Wilmington, DE). RNA was amplified and labelled according to a previous protocol [21]. Briefly, 100 ng of total RNA was used to synthesize the double-strand cDNA. RNA was amplified by in vitro transcription using Ambion’s MessageAmp II aRNA Amplification Kits (Life Technologies, Austin, TX, USA). Then, aRNA was reverse-transcribed into cDNA and further labelled with Klenow enzyme. cDNA from VPA-pretreated K562 cells was labelled with Cy3-deoxycytidine triphosphate and cDNA from control K562 cells was labelled with Cy5-deoxycytidine triphosphate. Fluorescent dye-labelled cDNA was hybridised to an Agilent SurePrint G3 Human GE 8 × 60 K Microarray. Hybridisation, scanning, and washing were carried out on Agilent’s Microarray Platform according to Agilent’s standard protocols. The array data were extracted with Agilent Feature Extraction software. After global mean normalization, probes with an intensity <400 were filtered out for further analysis. Differentially expressed genes were further analysed based on a significant enrichment of GO terms using hypergeometric distribution in the R language package software.

2.4. Validation of Differentially Expressed Genes

To explore the mechanism of the observed VPA-inducing apoptosis, eight differentially expressed genes involved in apoptosis were selected for validation of the microarray results by quantitative real-time PCR (qRT-PCR). Total RNAs were first separated from K562 cell line treated with or without VPA. The isolated total RNAs were treated with DNase and exclusively used for validation of the microarray results. Following DNase treatment, 1 μg of DNase-treated RNA was reverse-transcribed cDNA with oligo (dT)15 using M-MLV reverse transcriptase (Life Technologies) in a total volume of 20 μL reaction volume. After reverse transcription reaction, 1 μL of the reaction mixture was used for a qPCR program of 45 cycles consisting of melting (30 s at 94°C), annealing (30 s at 58°C), and extension (30 s at 72°C). The 20 μL reaction mixture contained 1 × PCR Buffer (Mg²⁺ Plus), 0.5 μM of forward primer, 0.5 μM of reverse primer (Table 1), 200 μM of each dNTP, and Eva Green Master Mix in a LightCycler 480 Real-Time PCR System (Roche Applied Science). Data were analysed by the 2−ΔΔCt method [22]. Glyceraldehyde 3-phosphate dehydrogenase (GADPH) was used as a reference gene. All other results were shown as fold change relative to GADPH control.

Table 1: Primer pairs for qRT-PCR validation of differentially expressed genes involved in apoptosis.

Gene symbol
GenBank Accession

Primer (5′3′)

GSN
NM_198252
F
CGAAGAGGTTCCTGGTGAGC

R
GCTGGGTCCGTCTCGATGTA

TP53INP1
NM_033285
F
CTTCTCATTGAACATCCCAGCA

R
TGGGTTGTTCCAGAAAAGTTGTA

TNFSF9
NM_003811
F
GAGCTTTCGCCCGACGAT

R
GCTCCTTCGTGTCCTCTTTGTAG

DDIT4
NM_019058
F
CCGAACTCCCACCCCAGA

R
GCAACGACACCCCATCCA

POLR2G
NM_002696
F
AAGATTGTGGGGACCCGTGT

R
GGATGAGTTGCCTTCCTTGC

LGALS1
NM_033101
F
CAACCCTCGCTTCAACGC

R
GGCGGTTGGGGAACTTGA

BBC3
NM_014417
F
CGGACGACCTCAACGCAC

R
GCACCTAATTGGGCTCCATCT

FOXO3
NM_001455
F
GGGGAACTTCACTGGTGCTAA

R
AACGGTATCACTGTCCACTTGC

PUF60
NM_014281
F
CTGTTTCCCGAGTCAGAGCG

R
CTCCTCTGTCACCTCCCCTTC

RRAGA
NM_006570
F
TGGTCATGTCAGATCCGTCG

R
AAGGCAATTCAGCATTTTCAGA

2.5. Statistics

The microarray data have been deposited in NCBI’s Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) and are accessible through GEO series accession number GSE32189. The significance between groups was analyzed by the Student’s -test (Microsoft Excel, Microsoft Corporation, Seattle, WA). A value of less than 0.05 () was considered statistically significant.

2.6. Comparison of Gene Expression Profiles with the CMAP Database

To further explore the mechanisms underlying the observed effects of VPA, the gene expression profile of VPA treatments was used to query the CMAP database (build 02), which contains more than 7,000 expression profiles representing 1,309 compounds [23]. The similarity between the gene expression profiles of the query signature and a CMAP instance is measured by the connectivity score, ranging from −1 to 1. A high positive connectivity score indicates inducement of the expression of the query signature by the corresponding drug, a high negative, inhibition. Prior to the query in the CMAP database via the Internet (http://www.broad.mit.edu/cmap/), the probe ID of Agilent SurePrint G3 Human GE 8 × 60 K Microarray was transformed into the probe ID defined by the Affymetrix GeneChip Human Genome U133A array according to the probe corresponding gene.

3. Results

3.1. Analysis of Apoptosis Using the Annexin V/PI Dual Staining Method

Using the Annexin V/PI dual staining method just as we previously reported, we measured apoptosis rate of K562 cells. In this study the similar increase in apoptosis was also observed in K562 cells exposed to 2 mM VPA for 48 hours (; , ), compared to the K562 control () [20].

3.2. Differentially Expressed Genes in K562 Cells Pretreated with VPA

Considering that the detection of differentially expressed genes is more easily in cells during the primary response phase of a drug exposure at certain doses that result in a 10–20% decrease in cell viability, we pretreated K562 cells with 2 mM VPA for 12 hours. This treatment resulted in an approximately 12% decrease in cell viability. The transcription profiles of K562 cells pretreated with or without VPA were analyzed using Agilent gene expression microarray. A total of 706 transcripts were identified as being significantly differentially expressed. The significant enrichment analysis of GO terms for the differentially expressed genes using the R language software package showed that apoptosis was the mostly enrichment function genes, and 34 genes involved in apoptosis were listed in Table 2.

Table 2: Differentially expressed genes involved in apoptosis detected by significant enrichment analysis of GO terms using the R language software package.

Genes
GenBank accession
Fold change
Gene description

PDCD7
NM_005707
−2.20
Programmed cell death 7

RRAGA
NM_006570
2.03
Ras-related GTP binding A

DAPK3
NM_001348
2.04
Death-associated protein kinase 3

PUF60
NM_014281
−2.00
Poly-U binding splicing factor 60kda isoform c

FOXO3
NM_001455
2.16
Forkhead box O3A

ARHGEF18
NM_015318
2.20
Rho/Rac guanine nucleotide exchange factor 18

BBC3
NM_014417
−2.18
BCL2 binding component 3 isoform 1

GSN
NM_198252
5.02
Gelsolin isoform c

LGALS1
NM_033101
2.58
Galectin-1

NFKB1
NM_003998
−3.33
Nuclear factor kappa-B, subunit 1 isoform 1

SH3GLB1
NM_016009
2.82
SH3-containing protein SH3GLB1

POLR2G
NM_002696
2.02
DNA directed RNA polymerase II polypeptide G

DDIT4
NM_019058
−2.23
Homo sapiens DNA-damage-inducible transcript 4

STEAP3
NM_182915
−2.52
Dudulin 2 isoform a

DRAM1
NM_018370
2.29
DNA-damage regulated autophagy modulator 1

RNF130
NM_018434
2.07
Ring finger protein 130 precursor

AVEN
NM_020371
−2.62
Cell death regulator Aven

TRIB3
NM_021158
−2.19
Tribbles 3

BAX
NM_138764
2.06
BCL2-associated X protein isoform beta

BID
NM_197966
−2.76
BH3 interacting domain death agonist isoform 1

PERP
NM_022121
9.42
PERP, TP53 apoptosis effector

SPTAN1
NM_003127
2.63
Spectrin, alpha, nonerythrocytic 1

BUB1B
NM_001211
−2.20
Budding uninhibited by benzimidazoles 1 beta

VDAC1
NM_003374
−2.37
Voltage-dependent anion channel 1

BCL2L12
NM_138639
2.19
BCL2-like 12 isoform 1

CUL1
NM_003592
−2.08
Cullin 1

YARS
NM_001040436
−2.30
Tyrosyl-TRNA synthetase

TNFSF9
NM_003811
2.38
Tumor necrosis factor (ligand) superfamily

IER3
NM_003897
2.17
Immediate early response 3

ARHGEF2
NM_004723
−2.10
Rho/Rac guanine nucleotide exchange factor 2

TP53INP1
NM_033285
3.03
Tumor protein p53 inducible nuclear protein 1

ITM2B
NM_021999
2.32
Integral membrane protein 2B

LITAF
NM_004862
2.25
Lipopolysaccharide-induced TNF-alpha factor

ESPL1
NM_012291
−2.11
Separase

−: indicate downregulation.

3.3. Validation of the Differentially Expressed Genes Using qRT-PCR

Based on the microarray results, eight differentially expressed genes involved in apoptosis were selected for qRT-PCR confirmation. Consistent with the results obtained from microarray analysis, PERP, GSN, TP53INP1, NXA4, FI6, ERPINB9, and CDKN2D were upregulated and NFKB1 was downregulated (Figure 1).

Figure 1: Alterations in the expression levels identified by microarray were confirmed using qRT-PCR. Fold changes between K562 cell line treated with and without valproate detected by microarray were compared with those measured by qRT-PCR. In qRT-PCR assay, RNA was isolated from K562 cell line treated with or without valproate exclusively for validation of differentially expressed genes. mRNA levels were normalized with GADPH and fold changes were calculated by dividing mRNA level of each K563 cell line treated with valproate by mean mRNA level from control samples in triplicate. Data are mean values ± standard deviations of the means from three independent experiments performed in triplicate.

3.4. Comparison of Gene Expression Profiles with the CMAP Database

To gain a better understanding of the anticancer and proapoptotic mechanisms of VPA exposure in a CML model, the gene expression profile of differentially expressed genes related to VPA treatment was used to compare with those of 1,309 compounds in the CMAP database [23, 24]. The similarity of the gene expression profile between the query signature and a CMAP instance were measured by the connectivity score, which ranges from −1 to 1. Only drugs with positive connectivity scores, which induce the expression of the query signature, are listed (Table 3). The results showed that the top ten drug hits belong to HDACi PI3K inhibitor, selective estrogen receptor modulator, antipsychotic or antimycobacterial agent, consistent with VPA’s pharmaceutical actions including recently reported action of inhibiting HDAC. Gene expression profile in K562 cells exposed to VPA observed in this study was most similar to that of HDACi and PI3K inhibitor in CMAP database.

Table 3: Top CMAP hits correlated with gene expression signature induced by VPA.

Rank
CMAP drug
P
Mean score
Enrichment
Specificity
Percent nonnull
Function

1
Vorinostat
0.00000
0.918
0.974
0.0201
100
HDACi

2
Trichostatin A
0.00000
0.817
0.957
0.0000
96
HDACi

3
Wortmannin
0.00000
0.335
0.695
0.0065
77
PI3K inhibitor

4
Ly-294002
0.00000
0.325
0.564
0.0201
65
PI3K inhibitor

5
VPA
0.00000
0.370
0.516
0.0000
61
Anti-epileptic; HDACi; Mood stabilizer

6
Raloxifene
0.00002
0.302
0.789
0.0000
85
SERM

7
Trifluoperazine
0.00002
0.303
0.594
0.1154
68
Antipsychotic

8
Scriptaid
0.00004
0.728
0.968
0.0185
100
HDACi

9
Rifabutin
0.00004
0.802
0.967
0.0052
100
Antimycobacterial agent

10
Thioridazine
0.00008
0.260
0.486
0.2831
65
Antipsychotic

The connectivity score, similarity measurement of the gene expression profile between the query signature and that of a CMAP instance, ranges from −1 to 1. Only drugs with positive connectivity scores, which induce the expression of the query signature, are listed. HDACi: histone deacetylase inhibitor. SERM: selective estrogen receptor modulator. VPA: valproic acid.

4. Discussion

VPA is a short-chain, branched fatty acid that has recently been described as a potent HDACi at therapeutic concentrations. Unlike other HDAC inhibitors, which are associated with various toxic side effects, VPA is clinically available. It can be taken orally, can cross the blood-brain barrier, and can be used for extended periods. HDAC inhibition is responsible for the acetylation of histones and nonhistone proteins such as hsp90 [25] and p53 [26]. This event has been associated with the induction of apoptosis in many models of leukaemia [27, 28]. Thus, the use of VPA prior to or concurrently with anticancer drugs may prove to be a beneficial treatment of leukaemia [8, 29]. The K562 cell line, originating from a patient diagnosed with CML in terminal blast crisis, is highly undifferentiated and of the granulocytic series. In this current study, we investigated the apoptosis rate of K562 cells exposed to VPA. The observed apoptosis result was in agreement with our previous reports [20].

Although the molecular mechanisms underlying antileukaemic activity of VPA in CLL [12] and AML [13, 14] models have been reported, genomewide gene expression changes in CML, however, have not been described. CML is a type of myeloproliferative disease associated with a characteristic chromosomal translocation called the Philadelphia chromosome. Conversely, CLL affects B-cell lymphocytes, and AML is characterized by the rapid growth of abnormal white blood cells that accumulate in the bone marrow and interfere with the production of normal blood cells. The observed clinical features and pathology of CML patients, as well as the therapy options, are quite different from those of CLL and AML patients. Consequently, the effect of VPA treatment on CML patients and the underlying molecular mechanisms may be different from effects observed in CLL and AML patients. Based on this information, we investigated the global gene expression profile of K562 cells exposed to VPA using gene expression microarray. Our microarray data demonstrated that VPA exposure resulted in widespread gene expression changes in this cell line. The significant enrichment analysis of GO terms revealed that genes involved in apoptosis were significantly altered by VPA treatment. This alteration was in accordance with the increase in cell apoptosis observed in VPA-pretreated K562 cells using flow cytometry.

To gain a better understanding of the antileukaemic mechanisms of VPA exposure in a CML model, the gene expression profile of differentially expressed genes related to VPA treatment was utilised to compare with those of 1,309 compounds in the CMAP database. CMAP analysis showed that the top ten drug hits included HDACi, PI3K inhibitor and VPA itself. These results were consistent with VPA’s pharmaceutical actions including recently reported action of inhibiting HDAC, indicating the reliability of our result of gene expression profile. In this study, VPA ranked in the top five drug hits. This result not only reflects the reliability of our result of gene expression profile, but also indicates that the connectivity map is a reliable tool for exploring molecular mechanism of drug action. The top four drug hits were vorinostat (suberoylanilide hydroxamic acid, ZolinzaTM), trichostatin A, Wortmannin, and LY-294002. Vorinostat and trichostatin A belong to HDACi. HDACi has a multitude of effects on cancer cells. It can induce cell-cycle arrest, promote differentiation, and stimulate tumor cell death [30]. It activates the apoptosis of cancer cells by the death-receptor (extrinsic) pathway and/or the mitochondrial (intrinsic) pathway [31]. The latter is induced by stress stimuli that disrupt the mitochondrial membrane. Recent studies have showed that HDACi treatment generates reactive oxygen species (ROS) and causes DNA damage in AML cell line, AML patient-derived blasts cultured ex vivo [32], and CML cell line [33], suggesting that ROS inducing intrinsic apoptosis might be one of important mechanisms underlying pro-apoptotic effect of HDACi. The intrinsic apoptosis may be dependent or independent on p53 pathway. We found that expressions of many genes involved in p53-dependent apoptotic pathway were up-regulated by VPA, for example, TP53INP1, PERP, and BAX (Table 2). TP53INP1 promotes p53/TP53 phosphorylation in response to double-strand DNA breaks. PERP plays a role as an effector in the TP53-dependent apoptotic pathway. BAX is regulated by p53 gene. It forms the mitochondrial pore leading to mitochondrial depolarization by oligomerization with BAK [34] and confers an essential gateway for the activation of caspases in the intrinsic apoptosis. Moreover, GO analysis showed that the genes involved in release of cytochrome c from mitochondria (GO: 0001836) were significantly enriched () (data not shown). It has been reported that VPA produces DNA oxidative damage in liver cells [35, 36], increases ROS formation, and induces apoptosis in postimplantation embryos [37]. Taken together, VPA might induce ROS and activate p53-dependent apoptosis of k562.

The expression of the BCR/ABL oncogene is the major pathogenic event in CML and the abnormal BCR/ABL tyrosine kinase promotes leukaemogenesis by inducing the phosphorylation of multiple downstream protein targets that mediate growth promoting and antiapoptotic signals, including MAPK/ERK and PI3K pathways [38–40]. PI3K is required only for the proliferation of BCR/ABL-dependent leukemic cells but not of normal hematopoietic cells [41, 42]. Wortmannin and LY294002 are the earliest and still widely employed inhibitors of PI3K. Studies have shown they are able to block proliferation of Ph1-positive cells [38, 43, 44] and enhance the antileukemia effect of Bcr-Abl tyrosine-kinase inhibitor [45–47], the initial choice of treatment for most CML patients. Based on CMAP analysis, we speculate that VPA might inhibit PI3K pathway too, and such inhibition might be one important mechanism by which VPA exerts antileukaemic action. We found that one member of the forkhead transcription factor O (FOXO) family, FOXO3, was upregulated by VPA (Table 2). FOXO family plays an important role in stress defence and it is negatively regulated by AKT-dependent phosphorylation [48].

In conclusion, we observed that VPA exposure altered mRNA expression in the K562 cell line. The CMAP analysis suggests that these alterations might be associated with inhibiting effect of VPA on HDAC as well as PI3K pathway. The data obtained in this study might provide clues to elucidate the molecular and therapeutic potential of VPA in leukaemia treatment, and the CMAP, by which we confirmed GO analysis and analyze the possible mechanisms of VPA antileukaemia, is a useful tool for exploring molecular mechanisms of drug actions.

Abbreviations

CML: | Chronic myelogenous leukaemia
VPA: | Sodium valproate
GO: | Gene ontology
CMAP: | Connectivity Map.

Conflict of Interests

The authors declare no potential conflict of interests.

Authors’ Contribution

These authors contributed equally to this work X. Z. Zhang and A. H. Yin.

Acknowledgment

This work was supported by Guangdong Province Science and Technology Fund no. 2009B030803033.

References

A. H. Wang, L. Wei, L. Chen et al., “Synergistic effect of bortezomib and valproic acid treatment on the proliferation and apoptosis of acute myeloid leukemia and myelodysplastic syndrome cells,” Annals of Hematology, vol. 90, no. 8, pp. 917–931, 2011.

B. Kircher, P. Schumacher, A. Petzer et al., “Anti-leukemic activity of valproic acid and imatinib mesylate on human Ph+ ALL and CML cells in vitro,” European Journal of Haematology, vol. 83, no. 1, pp. 48–56, 2009.

F. Buchi, R. Pastorelli, G. Ferrari et al., “Acetylome and phosphoproteome modifications in imatinib resistant chronic myeloid leukaemia cells treated with valproic acid,” Leukemia Research, vol. 35, no. 7, pp. 921–931, 2011.

G. Iacomino, M. C. Medici, and G. L. Russo, “Valproic acid sensitizes K562 erythroleukemia cells to TRAIL/Apo2L-induced apoptosis,” Anticancer Research A, vol. 28, no. 2, pp. 855–864, 2008.

E. S. Sung, A. Kim, J. S. Park, J. Chung, M. H. Kwon, and Y. S. Kim, “Histone deacetylase inhibitors synergistically potentiate death receptor 4-mediated apoptotic cell death of human T-cell acute lymphoblastic leukemia cells,” Apoptosis, vol. 15, no. 10, pp. 1256–1269, 2010.

L. Lagneaux, N. Gillet, B. Stamatopoulos et al., “Valproic acid induces apoptosis in chronic lymphocytic leukemia cells through activation of the death receptor pathway and potentiates TRAIL response,” Experimental Hematology, vol. 35, no. 10, pp. 1527–1537, 2007.

J. Stieglmaier, E. Bremer, C. Kellner et al., “Selective induction of apoptosis in leukemic B-lymphoid cells by a CD19-specific TRAIL fusion protein,” Cancer Immunology, Immunotherapy, vol. 57, no. 2, pp. 233–246, 2008.

E. Raffoux, A. Cras, C. Recher, et al., “Phase 2 clinical trial of 5-azacitidine, valproic acid, and all-trans retinoic acid in patients with high-risk acute myeloid leukemia or myelodysplastic syndrome,” Oncotarget, vol. 1, pp. 34–42, 2010.

G. Bug, M. Ritter, B. Wassmann et al., “Clinical trial of valproic acid and all-trans retinoic acid in patients with poor-risk acute myeloid leukemia,” Cancer, vol. 104, no. 12, pp. 2717–2725, 2005.

A. Kuendgen, M. Schmid, R. Schlenk et al., “The histone deacetylase (HDAC) inhibitor valproic acid as monotherapy or in combination with all-trans retinoic acid in patients with acute myeloid leukemia,” Cancer, vol. 106, no. 1, pp. 112–119, 2006.

C. Pilatrino, D. Cilloni, E. Messa et al., “Increase in platelet count in older, poor-risk patients with acute myeloid leukemia or myelodysplastic syndrome treated with valproic acid and all-trans retinoic acid,” Cancer, vol. 104, no. 1, pp. 101–109, 2005.

B. Stamatopoulos, N. Meuleman, C. De Bruyn et al., “Antileukemic activity of valproic acid in chronic lymphocytic leukemia B cells defined by microarray analysis,” Leukemia, vol. 23, no. 12, pp. 2281–2289, 2009.

F. L. Khanim, C. A. Bradbury, J. Arrazi et al., “Elevated FOSB-expression; a potential marker of valproate sensitivity in AML,” British Journal of Haematology, vol. 144, no. 3, pp. 332–341, 2009.

Y. C. Cheng, H. Lin, M. J. Huang, J. M. Chow, S. Lin, and H. E. Liu, “Downregulation of c-Myc is critical for valproic acid-induced growth arrest and myeloid differentiation of acute myeloid leukemia,” Leukemia Research, vol. 31, no. 10, pp. 1403–1411, 2007.

W. Fayad, M. Fryknäs, S. Brnjic, M. H. Olofsson, R. Larsson, and S. Linder, “Identification of a novel topoisomerase inhibitor effective in cells overexpressing drug efflux transporters,” Plos ONE, vol. 4, no. 10, Article ID e7238, 2009.

J. Gheeya, P. Johansson, Q. R. Chen et al., “Expression profiling identifies epoxy anthraquinone derivative as a DNA topoisomerase inhibitor,” Cancer Letters, vol. 293, no. 1, pp. 124–131, 2010.

H. Hieronymus, J. Lamb, K. N. Ross et al., “Gene expression signature-based chemical genomic prediction identifies a novel class of HSP90 pathway modulators,” Cancer Cell, vol. 10, no. 4, pp. 321–330, 2006.

Z. Wen, Z. Wang, S. Wang et al., “Discovery of molecular mechanisms of traditional Chinese medicinal formula Si-Wu-Tang using gene expression Microarray and Connectivity Map,” Plos ONE, vol. 6, no. 3, Article ID e18278, 2011.

D. K. Slonim, K. Koide, K. L. Johnson et al., “Functional genomic analysis of amniotic fluid cell-free mRNA suggests that oxidative stress is significant in Down syndrome fetuses,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 23, pp. 9425–9429, 2009.

X. Z. Zhang, A. H. Yin, X. Y. Zhu, Q. Ding, C. H. Wang, and Y. X. Chen, “Using an exon microarray to identify a global profile of gene expression and alternative splicing in K562 cells exposed to sodium valproate,” Oncology Reports, vol. 27, pp. 1258–1265, 2012.

T. A. Patterson, E. K. Lobenhofer, S. B. Fulmer-Smentek et al., “Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project,” Nature Biotechnology, vol. 24, no. 9, pp. 1140–1150, 2006.

K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method,” Methods, vol. 25, no. 4, pp. 402–408, 2001.

J. Lamb, E. D. Crawford, D. Peck et al., “The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease,” Science, vol. 313, no. 5795, pp. 1929–1935, 2006.

J. Lamb, “The Connectivity Map: a new tool for biomedical research,” Nature Reviews Cancer, vol. 7, no. 1, pp. 54–60, 2007.

J. J. Kovacs, P. J. M. Murphy, S. Gaillard et al., “HDAC6 regulates Hsp90 acetylation and chaperone-dependent activation of glucocorticoid receptor,” Molecular Cell, vol. 18, no. 5, pp. 601–607, 2005.

W. Gu and R. G. Roeder, “Activation of p53 sequence-specific DNA binding by acetylation of the p53 C-terminal domain,” Cell, vol. 90, no. 4, pp. 595–606, 1997.

A. Romanski, B. Bacic, G. Bug et al., “Use of a novel histone deacetylase inhibitor to induce apoptosis in cell lines of acute lymphoblastic leukemia,” Haematologica, vol. 89, no. 4, pp. 419–426, 2004.

A. Morotti, D. Cilioni, F. Messa et al., “Valproate enhances imatinib-induced growth arrest and apoptosis in chronic myeloid leukemia cells,” Cancer, vol. 106, no. 5, pp. 1188–1196, 2006.

H. Fredly, C. Stapnes Bjørnsen, B. T. Gjertsen, and Ø. Bruserud, “Combination of the histone deacetylase inhibitor valproic acid with oral hydroxyurea or 6-mercaptopurin can be safe and effective in patients with advanced acute myeloid leukaemia—a report of five cases,” Hematology, vol. 15, no. 5, pp. 338–343, 2010.

J. S. Carew, F. J. Giles, and S. T. Nawrocki, “Histone deacetylase inhibitors: mechanisms of cell death and promise in combination cancer therapy,” Cancer Letters, vol. 269, no. 1, pp. 7–17, 2008.

R. R. Rosato, J. A. Almenara, Y. Dai, and S. Grant, “Simultaneous activation of the intrinsic and extrinsic pathways by histone deacetylase (HDAC) inhibitors and tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) synergistically induces mitochondrial damage and apoptosis in human leukemia cells,” Molecular Cancer Therapeutics, vol. 2, no. 12, pp. 1273–1284, 2003.

L. A. Petruccelli, D. Dupéré-Richer, F. Pettersson, H. Retrouvey, S. Skoulikas, and W. H. Miller Jr., “Vorinostat induces reactive oxygen species and dna damage in acute myeloid leukemia cells,” Plos ONE, vol. 6, no. 6, Article ID e20987, 2011.

B. Brodska and A. Holoubek, “Generation of reactive oxygen species during apoptosis induced by DNA-damaging agents and/or histone deacetylase inhibitors,” Oxidative Medicine and Cellular Longevity, vol. 2011, Article ID 253529, 7 pages, 2011.

J. M. Adams and S. Cory, “The Bcl-2 apoptotic switch in cancer development and therapy,” Oncogene, vol. 26, no. 9, pp. 1324–1337, 2007.

V. Tong, X. W. Teng, T. K. H. Chang, and F. S. Abbott, “Valproic acid II: effects on oxidative stress, mitochondrial membrane potential, and cytotoxicity in glutathione-depleted rat hepatocytes,” Toxicological Sciences, vol. 86, no. 2, pp. 436–443, 2005.

K. H. Schulpis, C. Lazaropoulou, S. Regoutas et al., “Valproic acid monotherapy induces DNA oxidative damage,” Toxicology, vol. 217, no. 2-3, pp. 228–232, 2006.

E. W. Tung and L. M. Winn, “Valproic acid increases formation of reactive oxygen species and induces apoptosis in postimplantation embryos: a role for oxidative stress in valproic acid-induced neural tube defects,” Molecular Pharmacology, vol. 80, pp. 979–987, 2011.

T. Skorski, P. Kanakaraj, M. Nieborowska-Skorska et al., “Phosphatidylinositol-3 kinase activity is regulated by BCR/ABL and is required for the growth of Philadelphia chromosome-positive cells,” Blood, vol. 86, no. 2, pp. 726–736, 1995.

J. M. Goldman and J. V. Melo, “Mechanisms of disease: chronic myeloid leukemia—advances in biology and new approaches to treatment,” The New England Journal of Medicine, vol. 349, no. 15, pp. 1451–1464, 2003.

X. Zhu, L. Wang, B. Zhang, J. Li, X. Dou, and R. C. Zhao, “TGF-β1-induced PI3K/Akt/NF-κB/MMP9 signalling pathway is activated in Philadelphia chromosome-positive chronic myeloid leukaemia hemangioblasts,” Journal of Biochemistry, vol. 149, no. 4, pp. 405–414, 2011.

A. M. Martelli, M. Nyåkern, G. Tabellini et al., “Phosphoinositide 3-kinase/ Akt signaling pathway and its therapeutical implications for human acute myeloid leukemia,” Leukemia, vol. 20, no. 6, pp. 911–928, 2006.

F. J. Giles and M. Albitar, “Mammalian target of rapamycin as a therapeutic target in leukemia,” Current Molecular Medicine, vol. 5, no. 7, pp. 653–661, 2005.

Q. Xu, S. E. Simpson, T. J. Scialla, A. Bagg, and M. Carroll, “Survival of acute myeloid leukemia cells requires PI3 kinase activation,” Blood, vol. 102, no. 3, pp. 972–980, 2003.

B. Calabretta and T. Skorski, “BCR/ABL regulation of PI-3 kinase activity,” Leukemia and Lymphoma, vol. 23, no. 5-6, pp. 473–476, 1996.

A. Klejman, L. Rushen, A. Morrione, A. Slupianek, and T. Skorski, “Phosphatidylinositol-3 kinase inhibitors enhance the anti-leukemia effect of STI571,” Oncogene, vol. 21, no. 38, pp. 5868–5876, 2002.

S. B. Marley, J. L. Lewis, H. Schneider, C. E. Rudd, and M. Y. Gordon, “Phosphatidylinositol-3 kinase inhibitors reproduce the selective antiproliferative effects of imatinib on chronic myeloid leukaemia progenitor cells,” British Journal of Haematology, vol. 125, no. 4, pp. 500–511, 2004.

E. Weisberg, L. Banerji, R. D. Wright et al., “Potentiation of antileukemic therapies by the dual PI3K PDK-1 inhibitor, BAG956: effects on BCR-ABL and mutant FLT3-expressing cells,” Blood, vol. 111, no. 7, pp. 3723–3734, 2008.

H. Lei and F. W. Quelle, “FOXO transcription factors enforce cell cycle checkpoints and promote survival of hematopoietic cells after DNA damage,” Molecular Cancer Research, vol. 7, no. 8, pp. 1294–1303, 2009.