Promoter methylation might shift the balance of Galectin-3 & 12 expression in de novo adult acute myeloid leukemia patients

Acute myeloid leukemia (AML) was reported as the most common type of leukemia among adults. Galectins constitute a family of galactose-binding proteins reported to play a critical role in many malignancies including AML. Galectin-3 and -12 are members of the mammalian galectin family. To understand the contribution of galectin-3 and -12 promoter methylation to their expression, we performed bisulfite methylation-specific (MSP)-PCR and bisulfite genomic sequencing (BGS) of primary leukemic cells in patients with de novo AML before receiving any therapy. Here, we show a significant loss of LGALS12 gene expression in association with promoter methylation. The lowest degree of expression was found in the methylated (M) group while the highest degree was in the unmethylated (U) group and the partially methylated (P) group expression lies in between. This was not the case with galectin-3 in our cohort unless the CpG sites analyzed were outside the frame of the studied fragment. We were also able to identify four CpG sites (CpG number 1, 5, 7& 8) in the promoter region of galectin-12; these sites must be unmethylated so that expression can be induced. As far as the authors know, these findings were not previously concluded in earlier studies.

Introduction

Galectins an evolutionary conserved family are classified into three structural groups. Galectins are characterized by their ability to bind specific carbohydrates involved in a variety of cellular functions including cancer (Verkerke et al., 2022). This family of proteins is encoded by the LGALS genes family in humans and acts as an important recognizing factor towards cancer-associated glycoproteins (Yang et al., 2008). Their expression is very firmly controlled. Altered galectin expression is a hallmark of many cancer cells (Katzenmaier et al., 2017). Galectins can either be tumor promoters or suppressors based on their target cells (Kopitz et al., 2001). Furthermore, these family members had a prognostic influence on various types of malignancies including leukemias (Pena et al., 2014; Chetry et al., 2022). Galectin-12 structurally belongs to the galectins group which contains two homologous carbohydrate recognition domains (CRDs). Galectin-12 is a galectin family member encoded by the LGALS12 gene with a partial expression in leukocytes and adipocytes (Xue et al., 2016). While, galectin-3 belongs to another galectin structural group (chimera group with only one CRD) was reported to play an important role in cancerous'-microenvironments, especially in acute myeloid leukemia (AML) (Burger, 2011; Icard et al., 2014; Evans and Calvi, 2015; Han et al., 2015; Pereira et al., 2015).

AML was reported as a disease with high heterogeneity, is also considered the highest hematologic malignancy in its fatality rate (Chen et al., 2013; Bray et al., 2018; Siegel et al., 2020). Despite the existence of various medications which were approved for patients with AML in the last few years, AML remains a condition with unmet medical needs (Bewersdorf and Abdel-Wahab, 2022).

It was reported that epigenetic dysregulation might contribute to the development of hematological malignancies, such as hypomethylation or increased methylation of the CpG islands in the promoter region of key genes (Gutierrez and Romero-Oliva, 2013). Epigenetic regulation of galectin-1,-3 &-7 was suggested to play a critical role in cancer progression (Chiariotti et al., 1994; Poirier et al., 2001; Margadant et al., 2012; Kim et al., 2013). Galectin-12 was shown to be silenced by DNA methylation in colorectal cancer (CRC) cell lines and primary samples (Xue et al., 2016; Katzenmaier et al., 2017). Moreover, previous studies showed that LGALS3 promoter CpG islands were heavily methylated in the early stages of prostate adenocarcinoma (Ahmed and Bandyopadhyaya, 2015). Interestingly, another study reported that the average methylation degree of five CpG sites in the LGALS3 gene regulatory region was significantly decreased in thyroid cancer tissues (Keller et al., 2013). Subsequently, the present study was conducted to address the analysis of methylation patterns in galectin-3 and -12 promoter regions in patients with de novo AML.

Aim of the study

The present study aims to investigate to what extent the methylation patterns of LGALS 3 & 12 promoter region are correlated with LGALS 3 & 12 gene expression in adult AML patients.

The research basis

This study was based on two previous findings:

First, the expression profiling of eight galectins was previously performed in adult AML. Interestingly, galectin-12 was the only galectin that showed a survival advantage in AML patients when overexpressed in peripheral blood (PB) (El Leithy et al., 2015). In addition, galectin-3 was almost exclusively downregulated in both PB& bone marrow (BM) (El Leithy et al., 2015; Abdelfattah et al., 2021).

Second, alteration of DNA methylation was frequently encountered in AML (Kim et al., 2013). Aberrant methylation of cytosine-5 at CpG sites were clustered in the gene promoter regions. Furthermore, DNA methylation was found to correlate with prognosis in AML (Toyota et al., 2001; Deneberg et al., 2010; Bacigalupo et al., 2013; Chattopadhyaya and Ghosal, 2022).

We conduct this study to find out whether reduced expression of galectins-3 and -12 are associated with methylation of CpG islands in their promoter region.

Materials and methods

Patient cohort and collection of samples

This cohort study included 171 samples; 73 BM and 98 PB samples from AML patients, 97 males (56.7%) and 74 females (43.3%) with a mean age of 38.5 years (SD 12.5) ranging from 18 to 63 years. All of these patients were presented to the inpatient clinic at the National Cancer Institute (NCI), Cairo University (CU), diagnosed between July 2012 to December 2017. Patients underwent routine laboratory investigations and imaging diagnoses and were classified according to the standard morphological and immunophenotyping (IPT) criteria. All samples were collected before treatment; patients were treated intensively with the standard protocol. Patients with acute promyelocytic leukemia (APL) were given All-trans retinoic acid (ATRA). Other FAB subtypes were given the 3 + 7 treatment protocol. Response to induction therapy was assessed between days 14 and 28 after induction therapy and none of them received hypomethylating therapy. Written informed consents were obtained from the patients or their legal guardians, and this study was approved by the ethical committee of NCI, CU, Egypt, and was in accordance with the 2011 Declaration of Helsinki (IRP Approval No. 201902012.4). The age and sex-matched group consists of 15 PB samples from healthy donors from the same hospital, and eight BM samples from volunteers for BM transplantation were selected as a control group.

Methods

The present study sample size was estimated according to sample size estimation using the G*Power program (University of Düsseldorf, Düsseldorf, Germany) which is related to our previous work (El Leithy et al., 2015; Abdelfattah et al., 2021). The molecular assays were done on BM and PB whole white blood cell pellets; the initial blast cells median was 54% (ranging from 30% to 98%). It should be mentioned that the patient number is not the same in each technique depending on the availability of the samples in the lab. However, all samples were selected out of consecutive cohorts which meet the study eligibility criteria. That’s to say methylation analysis and gene expression were done on the same patients, after normalization of the genes expression to healthy donors, while the methylation analysis was assessed only in the patient cohort. This was done to evaluate the effect of the genes promotor methylation on the corresponding gene expression in the same samples.

RNA extraction and quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR)

The total cellular RNA from the number of BM and PB blood samples was purified to profile and associate the expression of LGALS 3 & 12. This was done using Invitrogen™ TRIzol™ Reagent (Invitrogen™, Thermo Fisher Scientific). Afterwards it was reversely transcribed to cDNA using Applied Biosystems™ (Thermo Fisher Scientific). The qRT-PCR was performed using Applied Biosystems PowerUp™ SYBR™ Green Master Mix (Thermo Fisher Scientific) according to the manufacturer’s instruction on Applied Biosystems Step One™ Real-Time PCR System (Thermo Fisher Scientific). The sequences of the forward and reverse primers related to LGALS3, LGALS12, and GAPDH are shown in Table 1. The primers were designed using the Primer-BLAST tool available on the NCBI website (https://www.ncbi.nlm.nih.gov/tools/primer-blast/).

TABLE 1

TABLE 1. The sequences of the forward and reverse primers for LGALS 3, 12, & GAPDH.

Differential DNA methylation analysis

Genomic DNA extraction and bisulfite conversion: Genomic DNA was isolated from 52 BM samples using G-spin™ Total DNA Extraction Mini Kit (LiliF Diagnostic Products). Then genomic DNA was bisulfite-modified using Thermo Scientific™ EpiJET™ Bisulfite Conversion Kit (Thermo Fisher Scientific) according to the manufacturer’s instruction. Subsequently, to investigate the effect of the DNA methylation on galectin-3&12 expressions in AML patients, methylation was analyzed by two methods on two separate cohorts as follows: 1- methylation-specific (MSP)-PCR for LGALS3 & 12 in 24 cases and 2- bisulfite genomic sequencing (BGS) for LGALS12 in 28 cases.

Methylation-specific (MSP)-PCR for LGALS3 & 12

The methylation pattern of LGALS3 & 12 genes promoter region was carried out on 24 subjects. For performing MSP, four pairs of primers were applied and specified for the methylated and un-methylated targeted sequences as shown in Figures 1A, B. Primer sequences are shown in Table 2. All of the MSP primers were designed using the MethPrimer design tool (http://www.urogene.org/methprimer/). Then the specificity of the primers was applied using the “by search” tool.

FIGURE 1

FIGURE 1. LGALS3 & 12 promoter regions: (A) Promoter region of LGALS3: Chromosome 14: 55128400-55132801 & (B) Promoter region of LGALS12: Chromosome 11: 63490635-63492346.

TABLE 2

TABLE 2. Profile of primers used in MSP.

Amplification was carried out in a thermocycler. For controlling and optimizing the MSP reactions, the EpiTect® PCR Control DNA Kit (Qiagen, Hilden, Germany) was used according to the manufacturer’s instructions. Methylation in this case is considered based on the amplification shown by samples with M primer as well as with both M and U primers. Samples showing amplification with only U were considered as unmethylated. Bands corresponding to methylated partners determined complete methylation. Partial methylation is determined by bands corresponding with both M and U primers.

Genomic DNA methylation sequencing of galectin-12 as a validation method for its methylation pattern

The second cohort included twenty-eight BM samples at diagnosis; they were enrolled for targeted bisulfite sequencing. MethPrimer software, (http://www.urogene.org/cgi-bin/methprimer/methprimer.cgi), was used for designing a specific set of primer pairs (M13-tailed PCR and LGALS12 primers) shown in Table 3 that binds only to bisulfite-modified DNA. This analysis was focused on a genomic region rich in CpG islands (by using DBCAT software, http://dbcat.cgm.ntu.edu.tw), containing 11 CpGs at nt −445 to nt −213 upstream of exon 1 in the predicted promoter region of LGALS12. This region also includes binding site for the transcription factor well known as SP1 that binds to GC-rich motives of many promoters (Katzenmaier et al., 2017).

TABLE 3

TABLE 3. Profile of primer used for CpG methylation analysis (used in BGS).

PCR reaction was carried out on Applied Biosystems™ Veriti™ 96-Well Fast Thermal Cycler (Thermo Fisher Scientific). Sequencing was carried out using BigDye™ Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific) according to the manufacturer’s instructions after the cycle sequencing. Products were purified using The BigDye® XTerminator™ Purification Kit (Thermo Fisher Scientific). Applied Biosystems™ 3500 XL Genetic Analyzer (Thermo Fisher Scientific) carried out targeted automatic bisulfite sequencing reaction and sequencing data analysis.

Statistical analysis

The data were analyzed using a statistical package for social science ‘IBM SPSS Statistics for Windows, version 26 (IBM Corp., Armonk, N.Y., United States). Continuous normally distributed variables were represented as mean ± SD with a 95% confidence interval, while non-normal variables were summarized as median with 25 and 75 percentiles, and using the frequencies and percentage for categorical variables; a p-value < 0.05 was considered statistically significant. To compare the means of normally distributed variables between groups, the student’s t-test was performed, while the Mann-Whitney U test was used in non-normal variables. χ² test was used to determine the distribution of categorical variables between groups. Pearson correlation was done to measure if there was any linear association between LGALS3 & 12 gene expression. For MSP-PCR results, univariate analysis was conducted to determine the prognostic performance of each studied biomarker. In addition, Survival analysis was performed using a Kaplan-Meier test.

Results

LGALS3 & 12 expressions in AML patients

The profiling of LGALS3 & 12 in both PB and BM of the AML patients showed differential expression as shown in (Figures 2A,B). The LGALS3 gene expression in PB and BM showed an association in the gene down-regulation in BM more than in PB (82.6% vs. 66.3% p-value = 0.044). Where, LGALS12 gene expression exhibited more downregulation in PB than in BM (63.3% vs. 43.8% p-value = 0.012) (Table 4; Figures 2C,D).

FIGURE 2

FIGURE 2. LGALS3 & 12 expressions in AML patients: (A) Box plot showing differential expression of the LGALS3 gene in PB & BM; the X-axis represents the LGALS3 gene and the Y-axis shows BM and PB cohorts. (B) Box plot showing differential expression of the LGALS12 gene in PB & BM; the X-axis represents the LGALS12 gene and the Y-axis shows BM and PB cohorts. (C) Cluster bar showing association between the LGALS3 gene expression in PB & BM; the X-axis represents the percent and the Y-axis shows the LGALS3 gene expression in BM and PB cohorts. (D) Cluster bar showing association between LGALS12 gene expression in PB & BM; the X-axis represents the percent and the Y-axis shows the LGALS12 gene expression in BM and PB cohorts. (E) Scatter plot showing the correlation between LGALS3&12 genes expression within BM. The X-axis represents the LGALS3 gene expression and the Y-axis shows the LGALS12 gene expression. (F) Scatter plot showing the correlation between LGALS3&12 genes expression within PB; the X-axis represents the LGALS3 gene expression and the Y-axis shows the LGALS12 gene expression. Abbreviations: BM (Bone marrow); PB (Peripheral blood).

TABLE 4

TABLE 4. Expression profile of LGALS3 & 12 in both PB and BM of the AML patients.

By studying the correlation between LGALS3 & 12 genes expression in PB and BM we found that; no statistically significant correlation was found between both genes’ expression in BM (r = −0.036 and p-value = 0.812). Conversely in the PB, a statistically significant moderate positive correlation was found (r = 0.5 and p-value < 0.001), Figures 2E,F.

The methylation pattern of LGALS3 & 12 gene promoter region

(MSP)-PCR for LGALS3 exhibited that one (4%) out of 24 examined subjects were completely methylated (M), 5 (21%) were partially methylated (P) and 18 (75%) were unmethylated (U). While LGALS12 promoter region showed that 5 (21%) were completely M out of 24 examined subjects, 12 (50%) were P and 7 (29%) were U.

Because most of the studied group by MSP-PCR exhibited an unmethylated pattern (75%) for the LGALS3 promoter region, BGS as a validation step was carried out only for LGALS12.

Relation of MSP of LGALS 12 with clinicolaboratory data

In this study, we compared all the clinical characteristics of patients with LGALS12 methylation patterns (Table 5). The results showed that there was a positive association with the adverse genetic risk in P and M cases when compared with the U group (p-value = 0.05 and 0.001, respectively). Furthermore, there was a positive risk for M with Odd Ratio (OR) = 4 (p-value = 0.01). Likewise, there was a positive association with the adverse genetic risk in M cases when compared with the P group (p-value < 0.001). In addition, there was a positive risk for the M group with OR = 2.00 (p-value = 0.05).

TABLE 5

TABLE 5. Relation of the methylation pattern of galectin 12 and clinicolaboratory data of the studied group.

On the other hand, there was a negative association with intermediate genetic risk in M cases when compared with the U group (p-value = 0.01) with protection from the M with OR = 0.33 (p-value = 0.01). Similar, results were given when the M group was compared with the P group (p-value = 0.05), OR = 0.25 (p-value < 0.001).

Regarding the hepatomegaly and splenomegaly associations, the results showed that there was a positive association with M in the cases that had organomegaly of the liver and spleen when compared with the U and/or P groups (p-value = 0.01 and 0.01, respectively).

In addition, the results showed that there was a positive association with lymph nodes (LNs) involvement in P cases and M cases when compared with the U group (p-value = 0.001 and 0.001, respectively). Likewise, there was a positive association with lymph node-positive in M cases when compared with the P group (p-value < 0.001). Moreover, there was a positive risk for M with OR = 7.50 (p-value = 0.01).

Regarding the immunophenotyping (IPT), the results showed that there was a positive association with the myelomono IPT in M cases when compared with the P group (p-value = 0.05). We found that there was a negative association with mono IPT in M cases when compared with the U group (p-value = 0.01) with protection from the M (OR = 0.50 (p-value = 0.02)). Similar results were given when the M group was compared with the P group, OR = 0.29 (p-value < 0.001).

Besides, the results showed a positive association with the mortality rate in M cases when compared with the P group (p-value = 0.04).

Impact of methylation pattern of LGALS12 on its gene expression and patients’ overall survival

The association results showed a positive association between LGALS-12 downregulation and an increasing incidence of methylation patterns. LGALS-12 was downregulated in P or M cases when compared to U cases in PB samples (p-value = 0.03 and 0.01, respectively) and BM samples (p-value = 0.003 and 0.01, respectively). Moreover LGALS-12 was downregulated in M cases when compared to P cases in PB samples (p-value = 0.01) and BM samples (p-value = 0.01), (Table 6; Figure 3).

TABLE 6

TABLE 6. The association between the methylation pattern of LGALS12 and its gene expression.

FIGURE 3

FIGURE 3. Log2 fold change of the LGALS12 gene expression associated with the methylation pattern of its promoter region. Abbreviations: M (methylated); P (partially methylated), U (unmethylated), BM (Bone marrow); PB (Peripheral blood). (A) Box plot showing Log2 fold change of the LGALS12 gene expression associated with the methylation pattern of its promoter region in PB; the X-axis represents the methylation pattern of LGALS12 gene and the Y-axis shows Log2 fold change of the LGALS12 gene expression in PB cohort. (B) Box plot showing Log2 fold change of the LGALS12 gene expression associated with the methylation pattern of its promoter region in BM; the X-axis represents the methylation pattern of LGALS12 gene and the Y-axis shows Log2 fold change of the LGALS12 gene expression in BM cohort.

The linear regression analysis results confirmed the abovementioned association results. As the OR of the P or M cases was compared with U cases in PB samples = 0.88 (p-value = 0.04), 0.03 (p-value < 0.001) and OR of BM samples = 0.76 (p-value = 0.04), and 0.78 (p-value = 0.01), respectively. In addition, OR of M cases was compared with P cases in PB samples = 0.38 (p-value < 0.001) and BM samples = 0.81 (p-value = 0.01), Table 6.

Follow-up of cases was done for 20 months. The median follow-up time was 3.75 months (ranging from 0.07 to 19.28 months). The OS of AML patients was measured from the date of diagnosis until the date of death or censoring for patients alive at the last follow-up. Studying the relation of OS and LGALS-12 methylation pattern showed a very close time to death in M cases when compared with U cases (p-value < 0.001). The same results were given, in the case of M group compared with P (p-value < 0.001), Figure 4.

FIGURE 4

FIGURE 4. Overall survival (OS) and its relation to methylation pattern of LGALS12. Abbreviations: M (methylated); P (partially methylated), U (unmethylated), BM (Bone marrow); PB (Peripheral blood).

Galectin-12 methylated CpG islands

The total number of CpG sites analyzed was 308 in 28 AML patients (Figure 5). Patients were divided based on their expression pattern into 2 groups. In the first non-expressed group (23 patients), the total number of CpG sites studied in the non-expressed group was 253 mostly methylated 203; the unmethylated sites were only 50/353 (19.8%). Second, expressed group (five patients), the total number in this expressed group was 55 mostly unmethylated 39/55 (70.1%), and 16 sites were methylated.

FIGURE 5

FIGURE 5. Bisulfite Genomic Sequencing of the promoter region of LGALS12. Methylation status of CpGs of Gelactin-12 promoter. CpGs are shown as either black (methylated) or white (unmethylated) circles and are numbered by Roman numerals.

Analysis of methylation pattern per patient revealed that:

In the non-expressed group (23 patients) five patients had all eleven (11/11) CpG sites methylated. Three patients had 10/11 CpG sites methylated. Two patients had 9/11 sites methylated. Nine patients had 8/11 CpG sites methylated and four patients had 7/11 CpG sites methylated.

In the expressed group (five patients) one patient had 2/11 CpG sites methylated, two patients had 3/11 CpG sites methylated and two patients had 4/11 CpG sites methylated. None of the studied AML patients showed 5/11 or 6/11 methylated sites.

Analysis of methylation pattern per CpG locus revealed that:

In the unexpressed group; CpG number 1 9/23 were unmethylated, CpG number 2 &5 8/23 were unmethylated, CpG number 3 & 6 6/23 were unmethylated, CpG number. 4 5/23 were unmethylated, CpG number 7& 10 1/23 were unmethylated, CpG number 8 4/23 were unmethylated and CpG number nine all were methylated. CpG number 11 2/23 were unmethylated.

In the expressed group; CpG number 1, 5, 7 & 8 were all unmethylated, CpG number 4 & 6 were unmethylated in 4/5, CpG number 2 &11 were unmethylated in 3/5, CpG number 9&10 were unmethylated in 2/5 and CpG number 3 was unmethylated in 1/5.

Discussion

The stimulator to leukemogenesis is a result of coordinated alterations in epigenetic regulation including methylation. Both gene-specific and global methylation patterns could predict outcomes in patients in AML (Deneberg et al., 2010). Previous studies showed that abnormal expression of some galectins correlates with tumor growth, cell migration, invasion, tumor aggressiveness, metastasis, recurrence, and poor prognosis in various cancers including leukemia (Ahmed et al., 2009; Bacigalupo et al., 2013; Ruvolo, 2019).

In the present study, despite the fact there was a statistically significant association between LGALS3 &12 genes expression in both PB and BM, there was no significant correlation between the expression in BM & PB of these two galectins in our previous study (Abdelfattah et al., 2021). This novelty could be due to the relatively large sample size studied here. In the current study, LGALS3 & 12 were both downregulated (66.3% & 63.3% respectively) in PB, and LGALS3 was even more downregulated in the BM (82.6%) while LGALS12 was less downregulated in the BM (43.8%). This could possibly be attributed to the fact that its expression in the BM is not only in leukocytes but also is in adipocytes (Xue et al., 2016). Galectin-12 was found to be preferentially expressed by human adipocytes and functions as an intrinsic negative regulator of lipolysis. In addition to its important function as an intracellular regulator of sebocyte proliferation (Tsao et al., 2022).

The study shows that, LGALS3 gene expression in PB and BM was mostly downregulated but contrary to our findings, Cheng et al. (2013) demonstrated that a higher bone marrow LGALS3 protein expression was an independent unfavorable prognostic factor for OS in patients with AML in Taiwan. This difference could be attributed to racial disparity and ethnic variation (Balan et al., 2008). Unlike our finding Ruebel et al. (2005) found that LGALS3 gene expression is decreased upon methylation in its promoter region in some pituitary tumors, however, their study was only restricted to cell lines for different types of cancer.

In fact, galectin-3 functionality depends on its subcellular localization, whether nuclear, cytoplasmic, cell surface, or circulating. Whereby cytoplasmic and circulating galectin-3 provide the most cell growth promotion (Newlaczyl and Yu, 2011). Ruvolo (2019) elaborated on the survival advantage induced by galectin-3 in the leukemic niche. He & Farhad et al. (2018) showed the role of MSC-derived galectin 3 in the AML microenvironment (Farhad et al., 2018; Ruvolo, 2019).

Herein, this study examined the methylation pattern of galectin-3 using MSP-PCR which resulted in most of the cases being unmethylated (18/24), six cases were partially methylated and only one case was completely methylated, despite the predominant low expression pattern of LGALS3. The role of LGALS3 expression pattern in carcinogenesis was extensively investigated in previous studies (Ebrahim et al., 2014). Silencing of galectin-3 expression by methylation of its promoter was associated with early stages of prostate cancer (Cummings et al., 2022). Our study outcome however did not give the same results in AML adult patients, unless, methylated CpG sites might be present outside the studied region.

Regarding LGALS12, we analyzed the methylation pattern of its promoter region by 2 methods. In the first cohort (MSP-PCR), most of our cases 12/24 (50%) were in the P category. This could be attributed to contamination by normal cells as reported previously (Quesnel et al., 1998), or by the fact that methylation in the CpG islands was not consistent in all AML samples (Galm et al., 2005). Our validation cohort supports the second notion since the percentage of methylation in the 11 CpG sites varied among patients ranging from 100% (11/11) to 7/11 but not less than seven sites methylation in patients who did not express galectin-12. At least seven CpG loci out of the eleven were methylated in the non-expressed group. In addition, we identified four CpG sites (1, 5, 7& 8) in the promoter region of galectin-12. All four must be unmethylated so that LGALS12 expression can be induced. To the authors’ knowledge; it is the first time to report such novelty in AML patients, an assumption that needs to be proven.

Regarding the clinical data of the patients, in the same institution, a previous study showed a significant association between splenomegaly and a relatively higher LGALS3 expression (Abdelfattah et al., 2021). Galectin-3 is known to be a powerful chemoattractant for monocytes, macrophages, and dendritic cells (Sano et al., 2000; Hsu et al., 2009). Thus we can postulate that its relatively higher expression in the bone marrow might attract the cells in numerous numbers which will be successfully drained into the spleen causing its expansion.

Here, our results showed a significant association between the promoter methylation status of galectin-12 & splenomegaly, hepatomegaly and lymphadenopathy. All cases with lymph nodes (LNs) enlarged were methylated 5/5 either partially or completely methylated (p-value < 0.001) for each. None of LNs enlargement was in the unmethylated group (p-value < 0.001). Higher expression of LGALS12 was shown to cause cell cycle arrest and apoptosis; probably causing shrinkage of the spleen and lymph nodes (Yang et al., 2001). Interestingly all 3 cases in the adverse genetic risk group were either partially (one case) or completely methylated (2 cases) in galectin-12 (p-value < 0.001) while none were in the unmethylated group. Also, it is in accordance with Farzaneh et al (2022), they showed that methylation as a biological process influences gene expression by affecting the promoter activity in colorectal cancer (Ghadiri Moghaddam et al., 2022). We found LGALS12 expression in the bone marrow only is border line significantly associated with AML-M4 compared to the other AML subtypes (p-value = 0.05) (data not shown), however, our methylation analysis showed a statistically significant association between complete methylation and unmethylation in the mono subtype only (p-value = 0.01). Moreover, mortality rate was increased in the methylated group. This finding is consistent with our previous finding that, patients with higher LGALS12 expression have the better overall survival (El Leithy et al., 2015).

Conclusion

The methylation pattern of the promoter region affects the expression only in galectin-12 but not in galectin-3. Our findings identify that hypermethylation of galectin-12 promoter is a common event in de novo adult AML. The abnormally hypomethylated and over-expressed galectin-12 cases had a relatively overall survival advantage. Galectin-3 downregulation is not a consequence of promoter methylation. However, intergenic and out of studied fragment DNA methylation cannot be excluded.

Recommendation and future prospective

Galectin-12 promoter hypomethylation and relative over-expression showed an overall survival advantage in AML patients. The present study findings confirmation in other cohorts in the same and different populations is recommended. Consequently, future research for specifically targeting a hypomethylating therapy agent for methylated galectin-12 promoter region could be an advance in the treatment of AML. Furthermore, the prospective evaluation of the methylation status of the galectin-12 promoter region in AML patients will be highly recommended for adjusting the patient treatment protocol.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

Ethics statement

The studies involving human participants were reviewed and approved by Written informed consents were obtained from the patients or their legal guardians, and this study was approved by the ethical committee of NCI, CU, Egypt, and was following the 2011 Declaration of Helsinki (IRP Approval No. 201902012.4). The patients/participants provided their written informed consent to participate in this study.

Author contributions

Conceptualization: MA and AE; Data curation: MA, AE, RE, RN, MK, and AA; Formal Analysis: MA, AE, RE, RN, MK, and AA;Funding acquisition: AA, AB, MA, MA, RE, MA and AE; Methodology: AE, RE, RN, MK, and MA; Software:RE, HG, AA, and AE; Supervision: MA, AE, and AA; Visualization: AE; Writing–review and editing: MA, AA, and AE. All authors approved the final version of the manuscript.

Acknowledgments

The authors thank Omar Ellethy, MSc. Neuroscience, VU-Amsterdam-RU Nijmegen, for his support during manuscript final revision.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abdelfattah, M. M., Assem, M. M., El Leithy, A. A., Mahmoud, R., Allam, N. M., and Helwa, R. (2021). Galectins in acute myeloid leukemia: Gene expression profiling with clinical emphasis. JBUON 5, 6–13.

Google Scholar

Ahmed, H., and Bandyopadhyaya, G. (2015). Examination of the regulation of galectin-3 expression in cancer. Methods Mol. Biol. 1207, 355–369. doi:10.1007/978-1-4939-1396-1_24

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmed, H., Cappello, F., Rodolico, V., and Vasta, G. R. (2009). Evidence of heavy methylation in the galectin 3 promoter in early stages of prostate adenocarcinoma: Development and validation of a methylated marker for early diagnosis of prostate cancer. Transl. Oncol. 2, 146–156. doi:10.1593/tlo.09118

PubMed Abstract | CrossRef Full Text | Google Scholar

Bacigalupo, M. L., Manzi, M., Rabinovich, G. A., and Troncoso, M. F. (2013). Hierarchical and selective roles of galectins in hepatocarcinogenesis, liver fibrosis and inflammation of hepatocellular carcinoma. World J. Gastroenterol. 19, 8831–8849. doi:10.3748/wjg.v19.i47.8831

PubMed Abstract | CrossRef Full Text | Google Scholar

Balan, V., Nangia-Makker, P., Schwartz, A. G., Jung, Y. S., Tait, L., Hogan, V., et al. (2008). Racial disparity in breast cancer and functional germ line mutation in galectin-3 (rs4644): A pilot study. Cancer Res. 68, 10045–10050. doi:10.1158/0008-5472.CAN-08-3224

PubMed Abstract | CrossRef Full Text | Google Scholar

Bewersdorf, J. P., and Abdel-Wahab, O. (2022). Translating recent advances in the pathogenesis of acute myeloid leukemia to the clinic. Genes. Dev. 36, 259–277. doi:10.1101/gad.349368.122

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi:10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

Burger, J. A. (2011). Nurture versus nature: The microenvironment in chronic lymphocytic leukemia. Hematol. Am. Soc. Hematol. Educ. Program 2011, 96–103. doi:10.1182/asheducation-2011.1.96

CrossRef Full Text | Google Scholar

Chattopadhyaya, S., and Ghosal, S. (2022). DNA methylation: A saga of genome maintenance in hematological perspective. Hum. Cell. 35, 448–461. doi:10.1007/s13577-022-00674-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Zheng, R., Zhang, S., Zhao, P., Li, G., Wu, L., et al. (2013). Report of incidence and mortality in China cancer registries, 2009. Chin. J. Cancer Res. 25, 10–21. doi:10.3978/j.issn.1000-9604.2012.12.04

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, C.-L., Hou, H.-A., Lee, M.-C., Liu, C.-Y., Jhuang, J.-Y., Lai, Y.-J., et al. (2013). Higher bone marrow LGALS3 expression is an independent unfavorable prognostic factor for overall survival in patients with acute myeloid leukemia. Blood 121, 3172–3180. doi:10.1182/blood-2012-07-443762

PubMed Abstract | CrossRef Full Text | Google Scholar

Chetry, M., Bhandari, A., Feng, R., Song, X., Wang, P., and Lin, J. (2022). Overexpression of galectin2 (LGALS2) predicts a better prognosis in human breast cancer. Am. J. Transl. Res. 14, 2301–2316.

PubMed Abstract | Google Scholar

Chiariotti, L., Benvenuto, G., Zarrilli, R., Rossi, E., Salvatore, P., Colantuoni, V., et al. (1994). Activation of the galectin-1 (L-14-I) gene from nonexpressing differentiated cells by fusion with undifferentiated and tumorigenic cells. Cell. Growth Differ. 5, 769–775.

PubMed Abstract | Google Scholar

Cummings, R. D., Liu, F.-T., Rabinovich, G. A., Stowell, S. R., and Vasta, G. R. (2022). Galectins essentials of glycobiology [internet]. 4th edition.

Google Scholar

Deneberg, S., Grovdal, M., Karimi, M., Jansson, M., Nahi, H., Corbacioglu, A., et al. (2010). Gene-specific and global methylation patterns predict outcome in patients with acute myeloid leukemia. Leukemia 24, 932–941. doi:10.1038/leu.2010.41

PubMed Abstract | CrossRef Full Text | Google Scholar

Ebrahim, A. H., Alalawi, Z., Mirandola, L., Rakhshanda, R., Dahlbeck, S., Nguyen, D., et al. (2014). Galectins in cancer: Carcinogenesis, diagnosis and therapy. Ann. Transl. Med. 2, 88. doi:10.3978/j.issn.2305-5839.2014.09.12

PubMed Abstract | CrossRef Full Text | Google Scholar

El Leithy, A. A., Helwa, R., Assem, M. M., and Hassan, N. H. (2015). Expression profiling of cancer-related galectins in acute myeloid leukemia. Tumour Biol. 36, 7929–7939. doi:10.1007/s13277-015-3513-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, A. G., and Calvi, L. M. (2015). Notch signaling in the malignant bone marrow microenvironment: Implications for a niche-based model of oncogenesis. Ann. N. Y. Acad. Sci. 1335, 63–77. doi:10.1111/nyas.12562

PubMed Abstract | CrossRef Full Text | Google Scholar

Farhad, M., Rolig, A. S., and Redmond, W. L. (2018). The role of Galectin-3 in modulating tumor growth and immunosuppression within the tumor microenvironment. Oncoimmunology 7, e1434467. doi:10.1080/2162402X.2018.1434467

PubMed Abstract | CrossRef Full Text | Google Scholar

Galm, O., Wilop, S., Lüders, C., Jost, E., Gehbauer, G., Herman, J. G., et al. (2005). Clinical implications of aberrant DNA methylation patterns in acute myelogenous leukemia. Ann. Hematol. 84, 39–46. doi:10.1007/s00277-005-0005-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghadiri Moghaddam, F., Farajnia, S., Karbalaei-Mahdi, M., and Monir, L. (2022). Epigenetic insights in the diagnosis, prognosis, and treatment selection in CRC, an updated review. Mol. Biol. Rep. 49, 1–10. doi:10.1007/s11033-022-07569-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutierrez, S. E., and Romero-Oliva, F. A. (2013). Epigenetic changes: A common theme in acute myelogenous leukemogenesis. J. Hematol. Oncol. 6, 57. doi:10.1186/1756-8722-6-57

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Zhang, Y., Jia, T., and Sun, Y. (2015). Molecular mechanism underlying the tumor-promoting functions of carcinoma-associated fibroblasts. Tumour Biol. 36, 1385–1394. doi:10.1007/s13277-015-3230-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, D. K., Chernyavsky, A. I., Chen, H.-Y., Yu, L., Grando, S. A., and Liu, F. T. (2009). Endogenous galectin-3 is localized in membrane lipid rafts and regulates migration of dendritic cells. J. Invest. Dermatol 129, 573–583. doi:10.1038/jid.2008.276

PubMed Abstract | CrossRef Full Text | Google Scholar

Icard, P., Kafara, P., Steyaert, J. M., Schwartz, L., and Lincet, H. (2014). The metabolic cooperation between cells in solid cancer tumors. Biochimica Biophysica Acta (BBA)-Reviews Cancer 1846, 216–225. doi:10.1016/j.bbcan.2014.06.002

CrossRef Full Text | Google Scholar

Katzenmaier, E. M., Kloor, M., Gabius, H. J., Gebert, J., and Kopitz, J. (2017). Analyzing epigenetic control of galectin expression indicates silencing of galectin-12 by promoter methylation in colorectal cancer. IUBMB Life 69, 962–970. doi:10.1002/iub.1690

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, S., Angrisano, T., Florio, E., Pero, R., Decaussin-Petrucci, M., Troncone, G., et al. (2013). DNA methylation state of the galectin-3 gene represents a potential new marker of thyroid malignancy. Oncol. Lett. 6, 86–90. doi:10.3892/ol.2013.1312

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. J., Hwang, J. A., Ro, J. Y., Lee, Y. S., and Chun, K. H. (2013). Galectin- 7 is epigenetically regulated tumor suppressor in gastric cancer. Oncotarget 4, 1461–1471. doi:10.18632/oncotarget.1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Kopitz, J., von Reitzenstein, C., Andre, S., Kaltner, H., Uhl, J., Ehemann, V., et al. (2001). Negative regulation of neuroblastoma cell growth by carbohydrate-dependent surface binding of galectin-1 and functional divergence from galectin-3. J. Biol. Chem. 276, 35917–35923. doi:10.1074/jbc.M105135200

PubMed Abstract | CrossRef Full Text | Google Scholar

Margadant, C., van den Bout, I., van Boxtel, A. L., Thijssen, V. L., and Sonnenberg, A. (2012). Epigenetic regulation of galectin-3 expression by β1 integrins promotes cell adhesion and migration. J. Biol. Chem. 287, 44684–44693. doi:10.1074/jbc.M112.426445

PubMed Abstract | CrossRef Full Text | Google Scholar

Newlaczyl, A. U., and Yu, L. G. (2011). Galectin-3--a jack-of-all-trades in cancer. Cancer Lett. 313, 123–128. doi:10.1016/j.canlet.2011.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Pena, C., Figueroa, J. A., Hosiriluck, N., Suvorava, N., Trotter, K., Reidy, A., et al. (2014). Galectins as therapeutic targets for hematological malignancies: A hopeful sweetness. Ann. Transl. Med. 2, 87. doi:10.3978/j.issn.2305-5839.2014.09.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereira, E. R., Jones, D., Jung, K., and Padera, T. P. (2015). The lymph node microenvironment and its role in the progression of metastatic cancer. Semin. Cell. Dev. Biol. 38, 98–105. doi:10.1016/j.semcdb.2015.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirier, F., Bourin, P., Bladier, D., Joubert-Caron, R., and Caron, M. (2001). Effect of 5-azacytidine and galectin-1 on growth and differentiation of the human b lymphoma cell line bl36. Cancer Cell. Int. 1, 2. doi:10.1186/1475-2867-1-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Quesnel, B., Guillerm, G., Vereecque, R., Wattel, E., Preudhomme, C., Bauters, F., et al. (1998). Methylation of the p15(INK4b) gene in myelodysplastic syndromes is frequent and acquired during disease progression. Blood 91, 2985–2990. doi:10.1182/blood.v91.8.2985.2985_2985_2990

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruebel, K. H., Jin, L., Qian, X., Scheithauer, B. W., Kovacs, K., Nakamura, N., et al. (2005). Effects of DNA methylation on galectin-3 expression in pituitary tumors. Cancer Res. 65, 1136–1140. doi:10.1158/0008-5472.CAN-04-3578

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruvolo, P. P. (2019). Galectins as regulators of cell survival in the leukemia niche. Adv. Biol. Regul. 71, 41–54. doi:10.1016/j.jbior.2018.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sano, H., Hsu, D. K., Yu, L., Apgar, J. R., Kuwabara, I., Yamanaka, T., et al. (2000). Human galectin-3 is a novel chemoattractant for monocytes and macrophages. J. Immunol. 165, 2156–2164. doi:10.4049/jimmunol.165.4.2156

PubMed Abstract | CrossRef Full Text | Google Scholar

Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer statistics, 2020. CA Cancer J. Clin. 70, 7–30. doi:10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

Toyota, M., Kopecky, K. J., Toyota, M. O., Jair, K. W., Willman, C. L., and Issa, J. P. (2001). Methylation profiling in acute myeloid leukemia. Blood 97, 2823–2829. doi:10.1182/blood.v97.9.2823

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsao, C.-H., Hsieh, W.-C., Yang, R.-Y., Lo, Y.-H., Tu, T.-J., Ke, L.-Y., et al. (2022). Galectin-12 modulates sebocyte proliferation and cell cycle progression by regulating cyclin A1 and CDK2. Glycobiology 32, 73–82. doi:10.1093/glycob/cwab100

PubMed Abstract | CrossRef Full Text | Google Scholar

Verkerke, H., Dias-Baruffi, M., Cummings, R. D., Arthur, C. M., and Stowell, S. R. (2022). Galectins: An ancient family of carbohydrate binding proteins with modern functions. Methods Mol. Biol. 2442, 1–40. doi:10.1007/978-1-0716-2055-7_1

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, H., Yang, R. Y., Tai, G., and Liu, F. T. (2016). Galectin-12 inhibits granulocytic differentiation of human NB4 promyelocytic leukemia cells while promoting lipogenesis. J. Leukoc. Biol. 100, 657–664. doi:10.1189/jlb.1HI0316-134R

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, R.-Y., Hsu, D. K., Yu, L., Ni, J., and Liu, F. T. (2001). Cell cycle regulation by galectin-12, a new member of the galectin superfamily. J. Biol. Chem. 276, 20252–20260. doi:10.1074/jbc.M010914200

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, R. Y., Rabinovich, G. A., and Liu, F. T. (2008). Galectins: Structure, function and therapeutic potential. Expert Rev. Mol. Med. 10, e17. doi:10.1017/S1462399408000719

PubMed Abstract | CrossRef Full Text | Google Scholar