Tumor Biol. DOI 10.1007/s13277-015-3513-0

RESEARCH ARTICLE

Expression profiling of cancer-related galectins in acute myeloid leukemia Asmaa A. El Leithy 1 & Reham Helwa 2 & Magda M. Assem 1 & Nagwa H. A. Hassan 2

Received: 13 January 2015 / Accepted: 27 April 2015 # International Society of Oncology and BioMarkers (ISOBM) 2015

Abstract Acute myeloid leukemia (AML) is the most common type of leukemia in adults with the lowest survival rate of all the leukemias. It is a heterogeneous disease in which a variety of cytogenetic and molecular alterations have been identified. Some galectins were previously reported to have important roles in cancer-like neoplastic transformation, tumor cell survival, angiogenesis, and tumor metastasis. Previous studies have showed that some galectin family members play a role in various types of leukemia. The present study aims at evaluating and clarifying the diagnostic and prognostic value of the expression of cancer-related galectins in relation to the clinicopathological characters of AML patients. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect expression profile of eight galectin family members (galectin-1, -2, -3, -4, -8, -9, -12, and -13) in 53 newly diagnosed de novo AML patients. The samples were collected from the inpatient clinic at National Cancer Institute (NCI), Cairo University (CU), diagnosed between July 2012 and May 2013. Our results show that patients with lower LGALS12 gene expression have a lower overall survival than those with higher expression (P value 6.94 PLT (*10^9/L) PLT (*10^9/L) ≤27 >27 Initial BM blast (%) Initial BM blast ≤54 % >54 % Initial BM promyelocyte Present (0 % blast)

n

%

36 11

76.5 23.5

*47

*7 4 3

57 43

25 5 1

81 16 3

15 17 5 16

28.3 32.1 9.4 30.2

*31

Median

Range

29.75

1.75–414

26 27 Mean 6.94

49.1 50.9 SD ±1.99

27 26 Median 27

50.9 49.1 Range 3–302

27 26 Median 54 %

50.9 49.1 Range 0–98

27 26

50.9 49.1

5

9.4

AML acute myeloid leukemia, FAB French–American–British classification of acute myeloid leukemia

LGALS9, LGALS12, and LGALS13 by qRT-PCR. Peripheral blood (PB) samples from both AML and healthy donors were collected on Trizol and frozen at −80 °C until use. The expression of LGALS1, 2, 3, 4, 8, 9, 12, and 13 was determined before treatment. Among the 53 patients, 45 (84.9 %) received standard induction chemotherapy 3+7 type regimen (doxorubicin 45 mg/m2 per day on days 1 to 3 and cytarabine 100 mg/

Tumor Biol. Table2 Immunophenotyping (IPT) of AML patients

IPT*

n

%

CD 117 present CD34 present MHC CLASSII present CD7 present CD14 present

20 20 25 5 9

37.7 37.7 47.2 9,4 17

2 2 1 5 4 3

3.8 3.8 1.9 9,4 7.5 5.7

CD2 present HLA DR present CD64 present CD4 present CD11c present CD19 present

m2 per day on days 1 to 7). Twenty-six out of 45 patients (57.78 %) reached a complete remission (CR), 3 patients (5.7 %) received 2+5 type regimen (doxorubicin 45 mg/m2 per day on days 1 to 2 and cytarabine 100 mg/m2 per day on days 1 to 5). Two out of the three patients (66.7 %) reached CR. The remaining 5 patients received AIDA protocol (ATRA (45 mg/m2/day) associated with intravenous idarubicin (IDA) (12 mg/m2/day)×3 days), and 2 out of 5 patients (40 %) reached CR.

RNA isolation and qRT-PCR The total cellular RNA from peripheral blood samples was purified using Biozol Total RNA Extraction Reagent (BioFlux) and reversely transcribed to cDNA. The qRT-PCR was performed using SYBR Green Real Time PCR master mix (EXPRESS SYBR® GreenER™ qPCR Supermix, universal (Invitrogen™)) according to the manufacturer’s instruction on Mx3000P QPCR system (Agilent Technologies, California, USA). The sequences of the forward and reverse primers for LGALS1, 2, 3, 4, 8, 9, 12, 13 and GAPDH are shown in Table 3. The primers were retrieved from PrimerBank, Massachusetts General Hospital, The Center for Computational and Integrative Biology, Harvard Medical School [11–14]. The CT values were obtained and normalized to GAPDH. Then, the fold changes were calculated using 2−ΔΔCT method. All of the CT values were in the linear range of detection.

Heatmap visualization The mRNA expression as fold changes was clustered using Multi Experiment Viewer [15]. Expression levels of all genes were visualized in the same green–black–red scheme in order to evaluate relative expression levels in all samples.

Table 3 The sequences of the forward and reverse primers for LGALS1, 2, 3, 4, 8, 9, 12, 13 & GAPDH Gene

Primer sequence

LGALS1 Primer forward Primer reverse LGALS2

TCGCCAGCAACCTGAATCTC GCACGAAGCTCTTAGCGTCA

Primer forward Primer reverse LGALS3 Primer forward Primer reverse LGALS4 Primer forward Primer reverse LGALS8 Primer forward Primer reverse LGALS9 Primer forward Primer reverse LGALS12 Primer forward Primer reverse LGALS13

ATGACGGGGGAACTTGAGGT CAGGTTCAGCTTGTCTGTCC

Primer forward Primer reverse GAPDH Primer forward Primer reverse

TATTGCCTTCCGTTTCCGAGT GGCACGTAGTCTGTTGTCTCC

ATGGCAGACAATTTTTCGCTCC GCCTGTCCAGGATAAGCCC CGACGCTGCCTTACTACCAG CCAACCACAAAGTTCACGAAGA ATACGTGGGCATGTTCCTAGT CGGCCCTTTTGAAACGAGG GGACGGACTTCAGATCACTGT CCATCTTCAAACCGAGGGTTG GCCTGGGCAGGTCATCATAG GAGTTCTGTCTGCGAAGGAGG

CTGGGCTACACTGAGCACC AAGTGGTCGTTGAGGGCAATG

Statistical analysis Data was analyzed using SPSS win (Statistical Package for Social Science), version 21. Numerical data was expressed as mean and standard deviation (SD), median, and range as appropriate. Qualitative data was expressed as frequency and percentage. Chi-square (Fisher’s exact) test was used to examine the relation between qualitative variables as appropriate. Comparison between groups with respect to numerical variables was made using Student’s t test for not normally distributed variables while Mann–Whitney U test has been used for not normally distributed variables. Survival analysis was made using Kaplan–Meier method. Comparison between two survival curves used log rank test, P value ≤0.05 was considered significant, and all test was two-tailed. P value was corrected using Benferroni adjustment to avoid hyperinflation of type 1 error resulting from multiple comparisons. Overall survival (OS) was calculated from date of diagnosis till date of death or last follow-up. Disease-free survival (DFS)

Tumor Biol.

The profiling of LGALS1, 2, 3, 4, 8, 9, 12, and 13 in PB cells of the AML patients exhibited differential expression comparable to normal controls. The fold changes of the studied galectins were then analyzed to define lower and higher expression groups as shown later in the results. The expression profiling of LGALS1, 2, 3, 4, 8, 9, 12, and 13 is shown in a descending order of downregulation in Table 4. The fold changes were then subjected to clustering analysis as shown in Fig. 1.

expression profiling and response to chemotherapy showed no statistically significant association in all types of galectins. The median follow-up time was 1.25 months (ranging from 0.5 to 13.5 months); OS of AML patients is shown in Fig. 2. There is no statistically significant association between AML patient OS and clinicolaboratory data of studied group as shown in Table 5. Relation of OS and gene expression is shown in Table 6; patients with lower LGALS12 gene expression have a lower OS than those with higher expression (P value 29.75 Hb (g/dL) ≤6.94 >6.94 PLT (*10^9/L) ≤27

27 47 36 11

0.396

2.204

0.351 0.551

1.678 NR

15 17 5 16

0.350 0.441 0.800 0.304

1.743 2.204 NR 1.678

26 27

0.505 0.300

6.776 1.678

27 26

0.412 0.382

2.204 1.743

27

0.464

5.395

>27 Initial BM blast ≤54 % >54 % Initial BM promyelocyte Present (0 % blast) Absent (0 % promyelocyte)

26

0.340

1.678

27 26

0.497 0.313

5.395 1.217

5 48

0.800 0.368

NR 1.743

*NR means not reached OS

0.952

0.318

0.497

0.390

0.871

0.531

0.100

0.211

Tumor Biol. Table 6 Relation of overall survival (OS) and gene expression of studied group

Gene expression

n

Proportion surviving at 6 months

26 27

0.383 0.438

1.678 2.204

15

0.583

10.658

Downregulation LGALS3 Upregulation Downregulation LGALS4 Upregulation Downregulation LGALS8 Upregulation Downregulation LGALS9 Upregulation Downregulation LGALS12 Upregulation Downregulation LGALS13 Upregulation

38

0.334

1.743

10 43

0.549 0.359

NR 1.678

8 45

0.750 0.355

NR 1.678

19 34

0.237 0.489

1.118 2.368

25 28

0.303 0.481

1.184 2.368

6 47

1.711 0.461

NR 1.678

11

0.273

1.678

Downregulation

42

0.428

2.368

LGALS1 Upregulation Downregulation LGALS2 Upregulation

value 0.993

0.195

0.130

0.144

0.116

0.239

0.026

0.591

mammary cancer cells [23] and alter the cell cycle when added to human T lymphocytes [24]. On the other hand, it was also found that galectin-1 was associated with disease progression. It is also secreted by myeloid cells and was suggested to influence leukemic B cell biology [25]. In our study, galectin-1 was downregulated in 50.9 % of AML samples. According to the previous publication, the downregulation of galectin-1 could be a favor for cancer cells to skip apoptosis/cell cycle arrest. Fig. 3 OS and its relation to LGALS12 expression

Median survival estimate (months)

Also, in previous studies, galectin-3 was reported to be downregulated during progression of different types of solid tumors [26–29]. On the other side, others have reported higher expression levels of galectin-3 in tumors of the pancreas [30], thyroid [31], colon [32], and gastric tissues [33]. According to previous studies, alteration in galectin-3 gene expression can lead to promotion or inhibition of tumor proliferation, differentiation, or metastasis, depending on the cell types and tumors [26]. However, galectin-3 is highly expressed by chronic

Tumor Biol. Fig. 4 Disease-free survival of all AML patients

Table 7 Relation of disease-free survival (DFS) and clinicolaboratory data of studied group

Variable

n

Proportion surviving at 6 months

Median survival estimate (months)

All Gender Male Female Age ≤40

30

0.774

8.882

16 14

0.856 0.682

8.882 NR

15

0.857

8.882

>40 FLT3-ITD W/W W/ITD FAB classification M0&M1 M2 M3 M4 Initial CBC TLC (*10^9/L) ≤29.75 >29.75 Hb (g/dL) ≤6.94 >6.94 PLT (*10^9/L) ≤27

15 27 21 6

0.727

NR

0.774 0.667

8.882 NR

6 10 3 11

0.800 0.750 NR 0.771

NR NR NR NR

16 14

0.685 0.909

8.388 10.855

15 15

0.755 0.779

NR 8.882

16

0.848

8.882

>27 Initial BM blast ≤54 % >54 % Initial BM promyelocyte Present (0 % blast) Absent (0 % promyelocyte)

14

0.682

10.855

18 12

0.779 0.379

NR 7.895

3 27

NR 0.749

NR NR

P value − 0.732

0.666

0.846

0.559

0.307

0.917

0.905

0.316

0.251

**Number of patients achieved CR is 30

Tumor Biol. Table 8 Relation of disease-free survival (DFS) and gene expression of studied group

Gene expression

n

Proportion surviving at 6 months

Median survival estimate (months)

14 16

0.786 0.800

8.882 10.855

10

0.889

10.855

Downregulation LGALS3 Upregulation Downregulation LGALS4 Upregulation Downregulation LGALS8 Upregulation Downregulation LGALS9 Upregulation Downregulation LGALS12 Upregulation Downregulation LGALS13 Upregulation

20

0.728

8.882

8 22

0.833 0.745

NR 8.388

6 24

0.500 0.733

10.855 8.388

8 22

0.857 0.752

8.882 NR

12 18

0.875 0.728

8.882 8.388

6 24

0.667 0.699

NR 8.882

4

0.667

8.882

Downregulation

26

0.782

10.855

LGALS1 Upregulation Downregulation LGALS2 Upregulation

P value 0.744

0.618

0.412

0.391

0.546

0.532

0.175

0.809

myelogenous leukemia (CML) cells in the bone marrow. Interestingly, in a previous study, it was found that coculturing bone marrow stromal cells with CML cell lines induced the expression of galectin-3 and caused drug resistance in vitro. Similarly, it was also reported that galectin-3 upregulation in CML cells increases their proliferative and chemotactic capacity and their resistance to therapeutic drugs. In addition to this, galectin-3 facilitates CML cell migration and long-term residence in the bone marrow [34]. In AML, upregulation of galectin-3 in bone marrow was reported to be an independent poor prognostic factor, irrespective of age, leukocyte counts, or karyotype, and might serve as a new biomarker in AML [7]. This was not the case in the present study where galectin-3 was downregulated in the peripheral blood in 81.1 % of AML in the present cohort. This difference could be attributed to three reasons, first the fact that galectin-3 expression in the bone marrow is different than in the peripheral blood. Second, the difference could be due to the difference in ethnicity as our cases were Egyptians and there were Taiwan, or due to the different AML biology. Galectin-8 binds to integrins, inhibits cell adhesion, and induces apoptosis [35]. Moreover, Nagy et al. reached that galectin-8 plays a significant role in the biology of human colon cancers where a malignant colon tissue exhibits a

significantly lower galectin-8 concentration than normal or benign tissue within colon cancers [36]. Also, galectin-8 was identified as a major target for tumor-specific monoclonal antibodies raised against prostate and lung cancers [37, 38]. Downregulation of galectin-8 in 64.2 % of our samples could reflect one possible mechanism of evading apoptosis in AML patients. Rechreche et al. identified the human galectin-4 as a downregulated gene during colorectal carcinogenesis [39], where it was upregulated in hepatocellular carcinomas [40] and in gastric cancer cells with increased metastatic potential [41], as compared to the low level in the corresponding normal tissues. Hereby, ours is the first study to report down-expression of galectin-4 in the peripheral blood of AML (in 84.9 % of 53 patients). Moreover, by comparing the expression profile of the studied galectins and the clinical characteristics of AML patients, the expression of galectin-4 was found to be statistically associated with patient age. The expression of galectin-4 was higher in younger age group of AML (adjusted P value

Expression profiling of cancer-related galectins in acute myeloid leukemia.

Acute myeloid leukemia (AML) is the most common type of leukemia in adults with the lowest survival rate of all the leukemias. It is a heterogeneous d...
994KB Sizes 2 Downloads 10 Views