The Effects of Promoter Methylation on Downregulation of DAZAP2 in Multiple Myeloma Cell Lines

Abstract

Our previous studies had shown that DAZAP2 was profoundly downregulated in bone marrow mononuclear cells from multiple myeloma patients. In this report, we analyzed epigenetic changes in multiple myeloma cell lines to understand the molecular mechanisms underlying the downregulation of DAZAP2. Four multiple myeloma cell lines, KM3, MM.1S, OPM-2 and ARH-77, were studied. The results of methylation specific PCR (MSP) showed that the promoter of DAZAP2 was methylated for KM3, MM.1S, OPM-2 and unmethylated for ARH-77. The DAZAP2 promoter region was amplified to obtain a series of different length sequences. All of the amplified sequences were inserted to luciferase reporter vector. The constructs were transfected into COS-7 cells and the luciferase activities were measured to search for the core region of DAZAP2 promoter. Two CpG islands were found in DAZAP2 promoter region. The results of luciferase assay showed that CpG island 1 displayed weak transcriptional activity, whereas CpG island 2 exhibited strong transcriptional activity (273 folds) compared to the control. The sequence that covered both CpG islands 1 and 2 showed higher activity (1,734 folds) compared to the control, suggesting that the two islands had synergistic effect on regulating DAZAP2 expression. We also found that M. Sss I methylase could inhibit the luciferase activity, whereas demethylation using 5-aza-2′-deoxycytidine treatment rescued the expression of DAZAP2 for multiple myeloma cell lines. These data revealed that methylation of DAZAP2 promoter was involved in downregulation of DAZAP2 in multiple myeloma cells.

About the Authors:

Sai-Qun Luo

Contributed equally to this work with: Sai-Qun Luo, Jing-Ping Hu

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Jing-Ping Hu

Contributed equally to this work with: Sai-Qun Luo, Jing-Ping Hu

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Qiang Qu

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Jiang Li

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Wei Ren

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Jia-Ming Zhang

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Yan Zhong

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Wei-Xin Hu

* E-mail: [email protected]

Affiliation: Molecular Biology Research Center, Xiangya School of Medicine, Central South University, Changsha, Hunan, People’s Republic of China

Introduction

Multiple myeloma is a disease with the uncontrolled proliferation and accumulation of malignant plasma cells in the bone marrow. Multiple myeloma cells are characterized by a profound genetic instability resulting in a complex set of numerical and structural chromosomal abnormalities. Pathogenesis of multiple myeloma cells was a multistep process in which plasma cell suffered series of molecular and cellular changes [1], [2].Constitutive genetic alterations of multiple myeloma cells are important determinants of the biological behavior of multiple myeloma cells in their local microenvironment. The interaction of multiple myeloma cells with BMSCs and bone marrow accessory cells upregulated transcripts for cytokines such as IL-6 [3], [4], the anti-apoptotic protein MCL1 [5]; HGF and insulin-like growth factors (IGFs) [6]; and heat shock proteins (HSPs), which regulate the conformation and function of proliferative and anti-apoptotic proteins [7]. This increased autocrine production of cytokines, along with paracrine cytokine production from the BMSCs, stimulates proliferative and anti-apoptotic signaling cascades in multiple myeloma cells. Bone destruction in multiple myeloma is mediated by MIP1α. MIP1α is a potent inducer of osteoclast formation independent of RANKL, and promotes both RANKL-stimulated and IL6-stimulated osteoclast formation [8]. Dickkopf 1 (DKK1) and IL-3 may contribute to the inhibitory effects of multiple myeloma cells on osteoblast differentiation. DKK1 inhibits the canonical Wnt pathway, which mediates the differentiation of osteoblast progenitor cells, DKK1 expression of multiple myeloma cells could inhibit osteoblastogenesis. Increased DKK1 levels in bone marrow plasma and peripheral blood from patients with multiple myeloma are associated with focal bone lesions [9], [10]. This imbalance between bone formation and resorption results in osteolytic lesions, which is a hallmark of multiple myeloma.

The epigenetic alterations are among the earliest molecular abnormalities to occur during tumorigenesis. A variety of genetic changes had been detected in multiple myeloma, including alterations in methylation of p16 gene, a frequent epigenetic event in multiple myeloma patients [11]. CpG island methylation of DKK1 promoter is correlated with in several MM cell lines and in MM cells from advanced MM patients. Demethylation of the DKK1 promoter restores DKK1 expression, which results in inhibition of β-catenin/TCF-mediated gene transcription in MM lines [12]. The chromosomal region at 13q14.1-14.3 is a hotspot of deleting abnormalities, and it appears that some multiple myeloma-specific tumor suppressor genes may reside in this region [13]. Some of immunoglobulin heavy-chain (IgH) translocations were recurrent and had characterized the genomic breakpoints of seven t(4;14) translocations from multiple myeloma patients [14]. The rearrangements of the IgH locus at 14q32 were found in the majority cases of multiple myeloma. p18 was also found to be frequently deleted in multiple myeloma patients [15]. All these abnormalities represented downregulation of tumor suppressor genes.

[Figure omitted. See PDF.]

Table 1. Primers for amplification of DAZAP2 promoter.

https://doi.org/10.1371/journal.pone.0040475.t001

Growth factors and cytokines may also be involved in the regulation of plasma cell proliferation and tumor progression. For example, IL-1β, IL-6, IL-10, IFN-α, IGF-1 and TGF-γ promote proliferation of multiple myeloma cells, whereas IFN-β1 and APO-1/FAS inhibit the growth of multiple myeloma cells [16]–[19]. The cyclin D1 promoter was hypomethylated and hyperacetylated in expressing cell lines and patient samples. RNA polymerase II bound at IgH regulatory sequences can activate the cyclin D1 promoter by either long-range polymerase transferring or tracking [20]. RAS mutations often provide a genetic marker, if not a causal event, in the evolution of MGUS (monoclonal gammopathy of undetermined significance) to multiple myeloma and play a role in the transition from intramedullary to extramedullary tumor [21]. Cell adhesion molecules and cell surface antigens may further contribute to the localization of myeloma cells in the bone marrow [22].

DAZAP2 (deleted in azoospermia associated protein 2) had originally been identified as an interacting protein of germ-cell-specific RNA-binding proteins DAZ (deleted in azoospermia) [23].The evolutionary conserved DAZAP2 protein functions as a TCF-4 interacting partner. The knockdown of DAZAP2 could not only reduce the activity of Wnt-signaling as measured by Tcf/−catenin reporters but also additionally alter the expression of Wnt-signaling target genes [24]. DAZAP2 modulates the affinity of TCF-4 for its DNA-recognition motif in chromatin immunoprecipitation studies. Loss of DAZAP2 in embryos resulted in diminished expression of hoxb9 with a concurrent increase in the anterior marker otx2. DAZAP2 is required for FGF dependent posterior patterning. In contrast to FGF activity, DAZAP2 induction of hoxb9 is not blocked by loss of canonical Wnt-signaling, and increasing DAZAP2 level alters neural patterning and induces posterior neural markers [25].

Our previous studies showed that DAZAP2 was the most profoundly downregulated gene in bone marrow mononuclear cells from multiple myeloma patients [26], [27]. The results implicated a role for DAZAP2 as a potential tumor suppressor involved in the origin and development of multiple myeloma. However, the mechanism of its downregulation is so far unknown.

[Figure omitted. See PDF.]

Table 2. Primers for amplification of CpG islands.

https://doi.org/10.1371/journal.pone.0040475.t002

Materials and Methods

Cell Lines

The human multiple myeloma cell line ARH-77, OPM-2, MM.1S and African monkey kidney cell line COS7 were from the American Type Culture Collection (ATCC; Manassas, VA). The human multiple myeloma cell line KM3 was from the Cell Center of Shanghai Chuan-Xiang Biotechnology (Shanghai, China). The cells were cultured in medium RPMI 1640 or DMEM (Gibco) supplemented with 10% heat-inactivated fetal bovine serum (FBS; Gibco), 100 kU/L benzylpenicillin, and 100 mg/L streptomycin at 37°C in a humidified incubator containing 5% CO₂.

Reagents, Polymerases and Restriction Endonucleases, Reportor Gene Vectors and Primers

Taq DNA polymerases, restriction endonuclease such as Xho I, Nhe I, and T4 DNA ligase were from TaKaRa (Dalian, China). CpG methylation enzyme Sss I, AMV Reverse transcriptase and restriction endonuclease BstU I were from the New BioLabs. Agarose gel DNA purification kit was from Ambiogen. Plasmids purification mini kit, DNA Wizard clean-up system, the pGL2 Luciferase Reporter Vectors, Luciferase Assay System, β-galactosidase immuosorbent were from Promega. Lipofectamine™ LTX regent was from Invitrogen. 5-aza-2′-deoxycytidine was from Sigma-Aldrich.

pGL2 vector series containing luciferase reporter gene include of pGL2-Basic plasmid (lacking of promoter and enhancer) as negative control and pGL2-Control plasmid (including SV40 promoter and enhancer) as positive control. pSV-β-gal was used as an internal control to normalize the differences in transfection efficiency.

All primers were synthesised by Invitrogen.

Construction of Recombinant Plasmids

Genomic DNA was extracted from multiple myeloma cell lines. PCR amplification system were as follows: 1–2 µl of genomic DNA (approximately 100 ng),5 µl of 10×Buffer,4 µl of dNTP (2.5 mM), 1 µl of sense (F)/antisense (R) primers (10 µM) (Table 1 and 2), 2 U Taq DNA polymerase, and distilled H₂O was added up to 50 µl. PCR amplification conditions were as follows: hot-started at 94°C for 3 min,30 cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 30 s, followed by one cycle of 72°C for 5 min. PCR products were purified by agarose gel DNA purification kit. The amplifying DNA fragments and pGL2 vector were digested with Xho I and Nhe I respectively, and ligated in the presence of T4 DNA ligase.

Analysis of DAZAP2 Expression Level in Multiple Myeloma Cell Lines

0.5 µg of total RNA extracted from above-mentioned cell lines were mixed with 0.5 µg of oligo(dT)₁₅, then heated for 5 min at 65°C and chilled on ice for 5 min. The reaction system included 4 µl of 5×RT buffer, 20 U RNasin, 4 µl of dNTP (2.5 mM), 10 U AMV Reverse transcriptase (RT) and DEPC treated water was added to 20 µl. The mixture was incubated for 60 min at 42°C. Then the reaction was terminated by heating for 5 min at 95°C.

Fifty µl PCR reaction mixtures included 1.5 mM MgCl₂, 200 µM dNTPs, 1 µl of 10 µM specific primers, 2.5 U Taq DNA polymerase, 5 µl of cDNA template and distilled water was added to 50 µl. Reaction condition was as follows: hot-started for 2 min at 94°C,32 cycles of 94°C for 30 s, 56°C for 30 s and 72°C for 30 s, followed by one cycle of 72°C for 5 min. β-actin was used for normalization of the relative expression level of mRNA, and PCR reaction was stopped before plateau phase. Primers for DAZAP2 are 5′-ACA GCC AAC CTA CCC TGT GCA-3′ (forward) and 5′-CAT GAC TGC AAG CTG AGC AGC-3′ (reverse) Primers for β-actin are 5′-TGA CGG TCA GGT CAT CAT TAT CGG CAA TGA-3′ (forward) and 5′-TTG ATC TTC ATG GTG ATA GGA GCG AGG GCA-3′ (reverse).

[Figure omitted. See PDF.]

Figure 1. DAZAP2 gene expression in multiple myeloma cell lines KM3, MM.1S, OPM-2, and ARH-77.

(A) DAZAP2 gene expression in normal condition. (B) DAZAP2 gene expression after treatment with 5-aza-2′-deoxycytidine. The total RNA was isolated from multiple myeloma cells and RT-PCR was performed with the specific primers to detect the gene expression level of DAZAP2. The PCR products were ran on 1.5% agarose gel. β-actin was used for normalization and verification of sample loading. The length of amplification products for DAZAP2 is 400 bp and for β-actin is 260 bp. M:100 bp DNA ladder plus marker; AZA: 5-aza-2′-deoxycytidine; (−): without treatment of 5-aza-2′-deoxycytidine as control; (+): with treatment of 5-aza-2′-deoxycytidine. Each value represents mean ± S.D. of three independent experiments.

https://doi.org/10.1371/journal.pone.0040475.g001

[Figure omitted. See PDF.]

Figure 2. MSP analysis of DAZAP2 promoter methylation status in multiple myeloma cell lines KM3, MM.1S, OPM-2, and ARH-77.

MSP was performed with the specific primers to amplify the target sequence and a fragment with length of 214 bp was obtained. M: 100 bp DNA ladder plus marker; U: unmethylated sequence; M: methylated sequence; M+: positive methylation control; U+: negative methylation control.

https://doi.org/10.1371/journal.pone.0040475.g002

Multiple Myeloma Cell Lines were Treated with 5-aza-2′-deoxycytidine

5-aza-2′-deoxycytidine was dissolved in 50% acetic acid before use. Multiple myeloma cell lines KM3, MM.1S, OPM-2, and ARH-77(1×10⁵/ml)were cultured in the medium RPMI 1640 with 10% calf serum for 24 h. The cells were treated with 5 µM 5-aza-2′-deoxycytidine or with solvent as control (0 µM 5-aza-2′-deoxycytidine). The final concentration of acetic acid in the medium was 0.025%. Cells were washed with PBS twice, fresh medium and 5-aza-2′-deoxycytidine were added for every 24 hrs. After triple treatment, cells were collected and total RNA was isolated. RT-PCR was performed to detect the gene expression of DAZAP2 in KM3, MM.1S, OPM-2, and ARH-77 cells.

[Figure omitted. See PDF.]

Figure 3. Amplification of DAZAP2 promoter region.

DAZAP2 promoter region is located upstream of DAZAP2 transcription site. The PCR products were a series of different length sequences. S1∼S15 are the names of DNA fragments corresponding to specific primers in Table 1. (A) the PCR products running on agarose gel; (B) the schematics of amplified DNA fragments in the location of DAZAP2 promoter region.

https://doi.org/10.1371/journal.pone.0040475.g003

[Figure omitted. See PDF.]

Figure 4. The activity analysis of promoter region upstream of DAZAP2 gene.

The amplified DNA fragments in DAZAP2 promoter region were inserted to luciferase reporter gene vector pGL2-Basic and luciferase activity was detected. The empty vector pGL2-Basic was used as negative control. The vector pGL2-Control, containing SV40 promoter and enhancer, can express luciferase with high efficiency and was used as positive control. Correction value = (S_RLU-B_RLU)/(S_β-G-B_β-G); S: recombinant clone sample; B: negative control; RLU (relative light unit): luciferase activity; β-G: β- galactosidase activity. Results are the Mean ± S.D of 3 independent experiments and asterisk symbol indicates P<0.05(*) or P<0.01(**) compared with the control group.

https://doi.org/10.1371/journal.pone.0040475.g004

Analysis of Methylation Specific PCR (MSP)

Fifty µl of genomic DNA (2 µg) was denatured in 5.5 µl of fresh prepared NaOH (3 M) for 15 min at room temperature. The denatured DNA was mixed with freshly prepared 30 µl of 10 mM hydroquinone and 520 µl of 3 M sodium bisulfate. After incubated for 16 h at 50°C in the dark, the treated DNA was purified using DNA Wizard clean-up system and dissolved in 50 µl water. 33 µl of 10 M ammonium acetate and 270 µl anhydrous alcohol were added and then incubated at −80°C so as to precipitate DNA. The solution was centrifuged for 30 min and washed with ice-cold 70% ethanol. The pellet was resuspended in 20 µl distilled water and stored at −20°C.

Two µl of bisulfite-treated genomic DNA was mixed with 2.5 µl of 10×buffer, 2 µl of 2.5 mM dNTPs,0.5 µl of 10 µM primers M1/M2 or U1/U2, 0.3 µl of 1 U/µl Taq DNA polymerase and distilled water was added to 25 µl. Amplification conditions were as follows: hot-started at 94°C for 2 min, 40 cycles of 94°C for 30 s, 53°C for 30 s and 72°C for 30 s, followed by one cycle of 72°C for 5 min. 5 µl of PCR products were run on a 2% agarose gel stained with ethidium bromide to observe the results. Primers (M1/M2) for methylation analysis are 5′-TTT TTT TTA GTT GAG GGA TAT CGG-3′ and 5′-CTT CCT ATT CGA AAC TAC TTC CGT A-3′. Primers (U1/U2) for unmethylation analysis are 5′-TTT TTT TTA GTT GAG GGA TAT TGG-3′ and 5′-CTT CCT ATT CAA AAC TAC TTC CAT A-3′.

[Figure omitted. See PDF.]

Figure 5. Effects of methylated DAZAP2 promoter region on transcription activity.

(A) Methylated and unmethylated recombinant plasmids were digested by BstU I and run on agarose gel. M: 100 bp ladder plus marker; 1, 3, 5: methylated pGL2-DAZAP2/CpG1, CpG2 and DAZAP2; 2, 4, 6: unmethylated pGL2-DAZAP2/CpG1, CpG2 and DAZAP2. (B) Recombinant reporter plasmids bearing methylated DNA fragment were constructed (see materials and methods). COS-7 cells were transfected with recombinant reporter plasmids bearing methylated and unmethylated DNA fragment. pSV-β-gal plasmid was used as an internal control reporter to normalize the differences in transfection efficiency. pGL2-Control vector as positive control and pGL2-Basic vector as negative control were used respectively. Relative luciferase activity was measured in triplicates. Results are the Mean ± S.D of 3 independent experiments and asterisk symbol indicates P<0.05(*) or P<0.01(**) compared with the control group).

https://doi.org/10.1371/journal.pone.0040475.g005

Methylation of CpG Island Fragment and Construction of pGL2–DAZAP2/CpG1(Me), pGL2–DAZAP2/CpG2(Me), pGL2–DAZAP2(Me) Vectors

pGL2–DAZAP2/CpG1, pGL2–DAZAP2/CpG2 and pGL2–DAZAP2 were digested by Xho I and run on a 0.8% agarose gel with ethidium bromide, then purified by agarose gel DNA purification kit. Methylation assay was incubated for 4 hrs at 37°C as follows: 10 µl of 10×NE buffer II, 10 µl of 10× SAM, 5 µg of digested fragment, 0.5 µl of Sss I methylation enzyme and water was added to 100 µl. Ten µl of methylation mixture was added 0.5 U restriction endonuclease BstU I, then incubated for 1 hr at 37°C. Methylated reaction efficiency was observed by electrophoresis on 0.8% agarose gel.

DNA fragments methylated by Sss I were digested by Nhe I and inserted into linear pGL2-Basic vector which was digested by Nhe I/Xho I. Ligation products were purified by agarose gel DNA purification kit. The new recombinant plasmids with methylated fragments were named as pGL2-DAZAP2/CpG1(Me), pGL2-DAZAP2/CpG2(Me) and pGL2-DAZAP2(Me).

Transfection

COS7 cells in exponential growth phase were rinsed twice with PBS, treated with 0.25% trypsin and added to the DMEM media with 10% fetal bovine serum (FBS). The cells were plated on 24-well plates at density of 0.4×10⁵/well in 500 µl DMEM with 10% FBS. Twenty-four hours later, cells of 80% attachment were transfected with recombinant plasmid. Transfection reagent was diluted with serum and antibiotics-free DMEM medium, and then mixed with 0.4 µg plasmid DNA. The mixture above was incubated at room temperature for 20 min. After moving medium in the plates, cells were washed with PBS and then added into 24-well plates. After 6–8 hrs, medium was changed by fresh DMEM with 10% FBS.

Reporter Gene Activity Assay

After 24–48 hrs of COS7 cells transfection, cell lysate was prepared. DMEM was removed from plates and rinsed with PBS twice. Two hundred µl 1×RLB (reporter lysis buffer) was added per well and the plates were rocked for several times to ensure complete coverage of the cells. Cells were scraped from plates and transferred to 1.5 ml Eppendorf tube on ice. Vortex for 10 s and centrifugation at 12,000×g for 2 min were performed at 4°C. Supernatant was transferred to a new tube and stored at −80°C before use.

Luminometer was programmed to perform a 2-second measurement delay followed by a 10-second measurement read for luciferase activity. One hundred µl of the Luciferase Assay Reagent and 20 µl cell lysate were added into luminometer tube per sample and mixed by pipetting 2–3 times. Tube was placed in the luminometer and results were recorded. The assay was performed in triplicate, and each experiment was repeated at least three times. β- galactosidase activity assay was performed to correct the the value of luciferase activity. Ten µl of cell lysate and 140 µl of 1×RLB were mixed with 150 µl of 2×buffer (ice-cold) and pipetted 2–3 times. Mixture was incubated for 30 min at 37°C. Five hundred µl of 1 M sodium carbonate was added to the tube to terminate the reaction.

Statistical Analysis

All data were expressed as mean ± standard deviation (± SD). Statistical analyses were performed with One-way ANOVA for significance using SPSS11.5 software. P value of less than 0.05 (P<0.05) indicated a significant difference between groups.

Results

DAZAP2 Expression in Multiple Myeloma Cell Lines

Total RNA was extracted from multiple myeloma cell lines KM3, MM.1S, OPM-2 and ARH-77. RT-PCR analysis was performed and the results showed that all cell lines had the band of DAZAP2 (400 bp). The expression level of DAZAP2 in KM3, MM.1S, OPM-2 was lower than that in ARH-77 (Fig.1A).

Effect of 5-aza-2′-deoxycytidine on DAZAP2 Gene Expression in Multiple Myeloma Cell Lines

The above-mentioned multiple myeloma cell lines(1×10⁵/ml)were respectively seeded in plates, cultured in RPMI 1640 medium with 10% FBS. Twenty four hours later, cells were treated with demethylation reagent 5-aza-2′-deoxycytidine at concentration of 0 and 5 µM. The medium was changed and 5-aza-2′-deoxycytidine was added every 24 hrs for 3 times. Cells were collected and total cellular RNA was isolated 24 hrs later at 3rd treatment. RT-PCR was performed so as to evaluate DAZAP2 gene expression level in KM3, MM.1S, OPM-2 and ARH-77 cells. The results showed that the expression level of DAZAP2 in KM3, MM.1S, OPM-2 with treatment of 5-aza-2′-deoxycytidine was much more than that without treatment of 5-aza-2′-deoxycytidine (Fig. 1B). However, there was no change of DAZAP2 expression level in ARH-77 cells after treatment of 5-aza-2′-deoxycytidine.

MSP Analysis of DAZAP2 Promoter

Genomic DNA was extracted from multiple myeloma cell lines. MSP analysis was performed so as to analyze DAZAP2 promoter methylation status in KM3, MM.1S, OPM-2 and ARH-77 cells and results showed that both methylated and unmethylated product could be detected in KM3, MM.1S, OPM-2 cells, while in ARH-77 cells only unmethylated product could be detected. Results indicated that DAZAP2 promoter was partially methylated at least in KM3, MM.1S, OPM-2 cells, while unmethylated in ARH-77 (Fig.2). This could explain that why the expression level of DAZAP2 in KM3, MM.1S, OPM-2 was lower than that in ARH-77.

Amplification of DAZAP2 Promoter Region and Activity Analysis

The DAZAP2 promoter region is located in the upstream of DAZAP2 transcription site and is a GC-rich region (GenBank accession number: JQ425692). The primers were designed to amplified this region (Table 1) and the PCR products were a series of different length sequences (Fig. 3). Bioinformatic analyses showed that there were 2 CpG islands (island 1 and island 2) in DAZAP2 promoter region. To search for the core region of DAZAP2 promoter, all of the amplified sequences were inserted to luciferase reporter gene vector pGL2-Basic. After transfection of COS7 with recombinant plasmids, luciferase activity was detected. Recombinant plasmid containing S10 sequence displayed most active than others (Fig. 4), while island 2 is in it. These results were consistent with bioinformatic analysis of DAZAP2.

Effect of Methylation-treated DAZAP2 Promoter on Transcription Activity

Recombinant plasmids were linearized by Xho I digestion and methylated by Sss I. Methylation efficiency was assessed by BstU I digesting methylated linear fragment. Because BstU I hardly digested methylated cytosine (C^m) sequence, methylated linear fragment was integral. However unmethylated or incompletely methylated fragments were digested into fragments of varying sizes by BstU I (Fig. 5A).

To explore the effects of methylated DAZAP2 promoter region on transcription activity, recombinant reporter plasmids bearing methylated DNA fragment were constructed and COS7 cells were transfected. Compared to pGL2-Basic, luciferase activity of pGL2- DAZAP2/CpG1 was only 7.7 fold (P>0.05) while pGL2-DAZAP2/CpG2 was 273 fold (**P<0.01). Data showed that CpG island 2 exerted an extensive activity, however CpG island 1 displayed inconspicuous activity. Relative luciferase activity of pGL2-DAZAP2 (which covered island 1 and island 2) was 1734 fold of pGL2-Basic (**P<0.01). It was indicated that two CpG islands coalition exerted stronger promoter activity, suggesting that the two islands had synergistic effect on regulating DAZAP2 expression. Relative luciferase activity of pGL2-DAZAP2/CpG1(Me), pGL2-DAZAP2/CpG2(Me) and pGL2-DAZAP2(Me) had no obvious difference compared to pGL2-Basic (P>0.05), which indicated that promoter activity was apparently inhibited attribute to methylated CpG loci blocked the expression of luciferase gene (Fig. 5B).

Discussion

To further understand the molecular mechanism underlying this deregulation, we analyzed epigenetic changes associated with DAZAP2 in multiple myeloma, given that aberrant methylation of promoter contributes to tumorigenesis by inducing transcriptional suppression and tumor suppressor inactivation.

Most malignant features of cancer cells are triggered by activated oncogenes and the loss of tumor suppressors due to mutation or epigenetic inactivation. We previously reported that DAZAP2 was downregulated in newly diagnosed multiple myeloma and this may influence the growth and survival of multiple myeloma cell [26], [27]. The present study was conducted to evaluate the mechanism of DAZAP2 downregulation and determine the methylation status of its promoter in multiple myeloma cells. Two islands with rich GC box in the promoter of DAZAP2 gene were identified. These data indicated that DAZAP2 promoter was hypermethylated and suppressed the expression of DAZAP2 in multiple myeloma cells. We further correlated DAZAP2 expression with normal plasma cells and malignant myeloma cells, as well as the molecular subtypes which the dataset includes 8 genetic subtypes (MY, PR, LB, MS, HP, CD-1, CD-2, and MF) from 351 newly diagnosed myeloma cases [28]. Interestingly, we extended and confirmed our previous discovery that DAZAP2 was significantly downregulated in multiple myeloma cells by using a large uniform dataset (P = 0.004). The low expression of DAZAP2 was especially significant in the subgroups of MY, PR, LB, HP, and CD-1 [28]. This study warrants further investigation of DAZAP2 and its potential role in multiple myeloma. It is indicated that the promoter hypermethylation played a major role in downregulation of DAZAP2 gene expression in multiple myeloma.

DNA methyltransferases (DNMTs) play an important role in maintaining DNA methylation. Abnormal expression of DNMTs and their isoforms have been found in many types of cancer and different DNMT domains are responsible for targeting DNA methylation to specific regions of the genome [29]. An interest in the status of gene methylation has grown in tumor studies in recent years [30], [31]. In 1983, a relationship between the incidence of malignancy and the abnormality of DNA methylation status was confirmed. Some genes in tumor cells exhibit low methylation level, whereas others exhibit high methylation level [32]. Aberrant methylation of genes is an important mechanism for inactivating the genes involved in tumorigenesis [33], [34].

Research has indicated that aberrant gene methylation is an important mechanism of deregulation in tumor cells and methylation of pro-apoptotic genes may accelerate the malignant transformation of tumors. Several genes involved in cell signaling were aberrantly methylated in AML, and the p73 methylation rate was 13.3% [35]. The hypermethylation of the promoter plays a critical role in the inactivation of cancer-associated genes, and some methylated cancer-associated genes may become an independent prognostic factor for their correlated disease. The methylation status may be becoming a new biological marker and risk factor for ALL [36].

Methylation is one kind of reversible genetic variation. Demethylation can be used to reactivate the expression of some important tumor suppressor genes in tumors or precancerous lesions, therefore, it may play a role in tumor therapy or prevention. Demethylation reagents, such as 5-aza-2′-deoxycytidine, are capable of reactivating methylated genes. GSTP1 DNA methylation and protein expression status is correlated with 5-aza-2′-deoxycytidine treatment response in prostate cancer cells [37]. A DNA methylation inhibitor demethylated p73 and activated p21 protein to promote AML cell apoptosis [38]. Recently, 5-aza-2′-deoxycytidine had been used in clinical trials, especially in hematological malignancies. The demethylating drug has been shown to be effective in treating patients with myelodysplastic syndrome (MDS) and blast crisis of chronic myelogenous leukemia (CML-BL) [39], [40]. Recent studies have indicated that the therapeutic effect of 5-aza-2′-deoxycytidine is associated with not only demethylation but also DNA covalent conjugation. Therefore, 5-aza-2′-deoxycytidine induced structural changes in chromatin may be caused by histone acetylation and methylation [41].

The present work provided the direct evidence between methylation of DAZAP2 promoter and downregulation of DAZAP2 in multiple myeloma cell lines. DAZAP2 promoter methylation might be involved in tumorigenesis of multiple myeloma and this would be helpful for the diagnosis and treatment of multiple myeloma patients.

Acknowledgments

We are grateful to Dr. Yi-Wu Shi, who performed the initial experiments with DAZAP2 expression.

Author Contributions

Conceived and designed the experiments: WXH. Performed the experiments: SQL JPH QQ JL WR JMZ YZ WXH. Analyzed the data: SQL JPH QQ JL WXH. Contributed reagents/materials/analysis tools: SQL JL. Wrote the paper: WXH JPH.

Citation: Luo S-Q, Hu J-P, Qu Q, Li J, Ren W, Zhang J-M, et al. (2012) The Effects of Promoter Methylation on Downregulation of DAZAP2 in Multiple Myeloma Cell Lines. PLoS ONE 7(7): e40475. https://doi.org/10.1371/journal.pone.0040475

References

1. Hallek M, Bergsagel PL, Anderson KC (1998) Multiple myeloma: increasing evidence for a multistep transformation process. Blood 91: 3–21.M. HallekPL BergsagelKC Anderson1998Multiple myeloma: increasing evidence for a multistep transformation process.Blood91321

2. Kuehl WM, Bergsagel PL (2002) Multiple myeloma: evolving genetic events and host interactions. Nature Reviews Cancer 2: 175–187.WM KuehlPL Bergsagel2002Multiple myeloma: evolving genetic events and host interactions.Nature Reviews Cancer2175187

3. Chauhan D, Uchiyama H, Akbarali Y, Urashima M, Yamamoto K, et al. (1996) Multiple myeloma cell adhesion-induced interleukin-6 expression in bone marrow stromal cells involves activation of NF-κB. Blood 87: 1104–1112.D. ChauhanH. UchiyamaY. AkbaraliM. UrashimaK. Yamamoto1996Multiple myeloma cell adhesion-induced interleukin-6 expression in bone marrow stromal cells involves activation of NF-κB.Blood8711041112

4. Gunn WG, Conley A, Deininger L, Olson SD, Prockop DJ, et al. (2006) A crosstalk between myeloma cells and marrow stromal cells stimulates production of DKK1 and interleukin-6: a potential role in the development of lytic bone disease and tumor progression in multiple myeloma. Stem Cells 24: 986–991.WG GunnA. ConleyL. DeiningerSD OlsonDJ Prockop2006A crosstalk between myeloma cells and marrow stromal cells stimulates production of DKK1 and interleukin-6: a potential role in the development of lytic bone disease and tumor progression in multiple myeloma.Stem Cells24986991

5. Raje N, Kumar S, Hideshima T, Roccaro A, Ishitsuka K, et al. (2005) Seliciclib (CYC202 or R-roscovitine), a small-molecule cyclin-dependent kinase inhibitor, mediates activity via down-regulation of Mcl-1 in multiple myeloma. Blood 106: 1042–1047.N. RajeS. KumarT. HideshimaA. RoccaroK. Ishitsuka2005Seliciclib (CYC202 or R-roscovitine), a small-molecule cyclin-dependent kinase inhibitor, mediates activity via down-regulation of Mcl-1 in multiple myeloma.Blood10610421047

6. Mitsiades CS, Mitsiades NS, McMullan CJ, Poulaki V, Shringarpure R, et al. (2004) Inhibition of the insulin-like growth factor receptor-1 tyrosine kinase activity as a therapeutic strategy for multiple myeloma, other hematologic malignancies, and solid tumors. Cancer Cell 5: 221–230.CS MitsiadesNS MitsiadesCJ McMullanV. PoulakiR. Shringarpure2004Inhibition of the insulin-like growth factor receptor-1 tyrosine kinase activity as a therapeutic strategy for multiple myeloma, other hematologic malignancies, and solid tumors.Cancer Cell5221230

7. Mitsiades CS, Mitsiades NS, McMullan CJ, Poulaki V, Kung AL, et al. (2006) Antimyeloma activity of heat shock protein-90 inhibition. Blood 107: 1092–1100.CS MitsiadesNS MitsiadesCJ McMullanV. PoulakiAL Kung2006Antimyeloma activity of heat shock protein-90 inhibition.Blood10710921100

8. Choi SJ, Oba Y, Gazitt Y, Alsina M, Cruz J, et al. (2001) Antisense inhibition of macrophage inflammatory protein 1-α blocks bone destruction in a model of myeloma bone disease. J Clin Invest 108: 1833–1841.SJ ChoiY. ObaY. GazittM. AlsinaJ. Cruz2001Antisense inhibition of macrophage inflammatory protein 1-α blocks bone destruction in a model of myeloma bone disease.J Clin Invest10818331841

9. Tian E, Zhan FH, Walker R, Rasmussen E, Ma YP, et al. (2003) The role of the Wnt-signaling antagonist DKK1 in the development of osteolytic lesions in multiple myeloma. N Engl J Med 349: 2483–2494.E. TianFH ZhanR. WalkerE. RasmussenYP Ma2003The role of the Wnt-signaling antagonist DKK1 in the development of osteolytic lesions in multiple myeloma.N Engl J Med34924832494

10. Hideshima T, Mitsiades C, Tonon G, Richardson PG, Anderson KC (2007) Understanding multiple myeloma pathogenesis in the bone marrow to identify new therapeutic targets. Nature Reviews Cancer 7: 585–598.T. HideshimaC. MitsiadesG. TononPG RichardsonKC Anderson2007Understanding multiple myeloma pathogenesis in the bone marrow to identify new therapeutic targets.Nature Reviews Cancer7585598

11. Mateos MV, Garcia-Sanz R, Lopez-Perez R, Moro MJ, Ocio E, et al. (2002) Methylation is an inactivating mechanism of the p16 gene in multiple myeloma associated with high plasma cell proliferation and short survival. Br J Haematol 118: 1034–1040.MV MateosR. Garcia-SanzR. Lopez-PerezMJ MoroE. Ocio2002Methylation is an inactivating mechanism of the p16 gene in multiple myeloma associated with high plasma cell proliferation and short survival.Br J Haematol11810341040

12. Kocemba KA, Groen RWJ, van Andel H, Kersten MJ, Mahtouk K, et al. (2012) Transcriptional Silencing of the Wnt-Antagonist DKK1 by Promoter Methylation Is Associated with Enhanced Wnt Signaling in Advanced Multiple Myeloma. PLoS One 7(2): e30359.KA KocembaRWJ GroenH. van AndelMJ KerstenK. Mahtouk2012Transcriptional Silencing of the Wnt-Antagonist DKK1 by Promoter Methylation Is Associated with Enhanced Wnt Signaling in Advanced Multiple Myeloma.PLoS One7(2)e30359

13. Shaughnessy J, Tian E, Sawyer J, Bumm K, Landes R, et al. (2000) High incidence of chromosome 13 deletion in multiple myeloma detected by multiprobe interphase FISH. Blood 96: 1505–1511.J. ShaughnessyE. TianJ. SawyerK. BummR. Landes2000High incidence of chromosome 13 deletion in multiple myeloma detected by multiprobe interphase FISH.Blood9615051511

14. Sossey-Alaoui K, Kitamura E, Head K, Cowell JK (2002) Characterization of FAM10A4, a member of the ST13 tumor suppressor gene family that maps to the 13q14.3 region associated with B-Cell leukemia, multiple myeloma, and prostate cancer. Genomics 80: 5–7.K. Sossey-AlaouiE. KitamuraK. HeadJK Cowell2002Characterization of FAM10A4, a member of the ST13 tumor suppressor gene family that maps to the 13q14.3 region associated with B-Cell leukemia, multiple myeloma, and prostate cancer.Genomics8057

15. Kulkarni MS, Daggett JL, Bender TP, Kuehl WM, Bergsagel PL, et al. (2002) Frequent inactivation of the cyclin-dependent kinase inhibitor p18 by homozygous deletion in multiple myeloma cell lines: ectopic p18 expression inhibits growth and induces apoptosis. Leukemia 16: 127–134.MS KulkarniJL DaggettTP BenderWM KuehlPL Bergsagel2002Frequent inactivation of the cyclin-dependent kinase inhibitor p18 by homozygous deletion in multiple myeloma cell lines: ectopic p18 expression inhibits growth and induces apoptosis.Leukemia16127134

16. Klein B, Zhang XG, Lu ZY, Bataille R (1995) Interleukin-6 in human multiple myeloma. Blood 85: 863–872.B. KleinXG ZhangZY LuR. Bataille1995Interleukin-6 in human multiple myeloma.Blood85863872

17. Dewan MZ, Watanabe M, Terashima K, Aoki M, Sata T, et al. (2004) Prompt tumor formation and maintenance of constitutive NF-kappaB activity of multiple myeloma cells in NOD/SCID/gammacnull mice. Cancer Sci 95: 564–568.MZ DewanM. WatanabeK. TerashimaM. AokiT. Sata2004Prompt tumor formation and maintenance of constitutive NF-kappaB activity of multiple myeloma cells in NOD/SCID/gammacnull mice.Cancer Sci95564568

18. Anderson KC, Shaughnessy JD Jr, Barlogie B, Harousseau JL, Roodman GD (2002) Multiple myeloma. Hematology 1: 214–240.KC AndersonJD Shaughnessy JrB. BarlogieJL HarousseauGD Roodman2002Multiple myeloma.Hematology1214240

19. Wang MC, Fu XD, Li MX (2006) PI-3K/Akt/GSK-3β signaling cascades stimulated by insulin like growth factor-I contribute to multiple myeloma cells proliferation and survival. Chin Med J 119: 1226–1229.MC WangXD FuMX Li2006PI-3K/Akt/GSK-3β signaling cascades stimulated by insulin like growth factor-I contribute to multiple myeloma cells proliferation and survival.Chin Med J11912261229

20. Liu H, Wang J, Epner EM (2004) Cyclin D1 activation in B-cell malignancy: association with changes in histone acetylation, DNA methylation, and RNA polymerase II binding to both promoter and distal sequences. Blood 104: 2505–2513.H. LiuJ. WangEM Epner2004Cyclin D1 activation in B-cell malignancy: association with changes in histone acetylation, DNA methylation, and RNA polymerase II binding to both promoter and distal sequences.Blood10425052513

21. Rasmussen T, Kuehl M, Lodahl M, Johnsen HE, Dahl IM (2005) Possible roles for activating RAS mutations in the MGUS to MM transition and in the intramedullary to extramedullary transition in some plasma cell tumors. Blood 105: 317–323.T. RasmussenM. KuehlM. LodahlHE JohnsenIM Dahl2005Possible roles for activating RAS mutations in the MGUS to MM transition and in the intramedullary to extramedullary transition in some plasma cell tumors.Blood105317323

22. Hideshima T, Bergsagel PL, Kuehl WM, Anderson KC (2004) Advances in biology of multiple myeloma: clinical applications. Blood 104: 607–618.T. HideshimaPL BergsagelWM KuehlKC Anderson2004Advances in biology of multiple myeloma: clinical applications.Blood104607618

23. Tsui S, Dai T, Roettger S, Schempp W, Salido EC, et al. (2000) Identification of two novel proteins that interact with germ-cell-specific RNA-binding proteins DAZ and DAZL1. Genomics 65: 266–273.S. TsuiT. DaiS. RoettgerW. SchemppEC Salido2000Identification of two novel proteins that interact with germ-cell-specific RNA-binding proteins DAZ and DAZL1.Genomics65266273

24. Lukas J, Mazna P, Valenta T, Doubravska L, Pospichalova V, et al. (2009) Dazap2 modulates transcription driven by the Wnt effector TCF-4. Nucleic Acids Research 37: 3007–3020.J. LukasP. MaznaT. ValentaL. DoubravskaV. Pospichalova2009Dazap2 modulates transcription driven by the Wnt effector TCF-4.Nucleic Acids Research3730073020

25. Roche DD, Liu KJ, Harland RM, Monsoro-Burq AH (2009) Dazap2 is required for FGF-mediated posterior neural patterning, independent of Wnt and Cdx function. Developmental Biology 333: 26–36.DD RocheKJ LiuRM HarlandAH Monsoro-Burq2009Dazap2 is required for FGF-mediated posterior neural patterning, independent of Wnt and Cdx function.Developmental Biology3332636

26. Shi YW, Luo SQ, Peng JB, Huang CH, Tan DR, et al. (2004) The Structure, Expression and Function Prediction of DAZAP2, A Down-Regulated Gene in Multiple Myeloma. Gen Prot Bioinfo 2: 47–54.YW ShiSQ LuoJB PengCH HuangDR Tan2004The Structure, Expression and Function Prediction of DAZAP2, A Down-Regulated Gene in Multiple Myeloma.Gen Prot Bioinfo24754

27. Shi YW, Shen R, Ren W, Tang LJ, Tan DR, et al. (2007) The molecular features and expression of DAZAP2 in human multiple myeloma. Chin Med J 120: 1659–1665.YW ShiR. ShenW. RenLJ TangDR Tan2007The molecular features and expression of DAZAP2 in human multiple myeloma.Chin Med J12016591665

28. Zhan F, Huang Y, Colla S, Stewart JP, Hanamura I, et al. (2006) The molecular classification of multiple myeloma. Blood 108: 2020–2028.F. ZhanY. HuangS. CollaJP StewartI. Hanamura2006The molecular classification of multiple myeloma.Blood10820202028

29. Choi SH, Heo K, Byun HM, An WJ, Lu WG, et al. (2011) Identification of preferential target sites for human DNA methyltransferases. Nucleic Acids Research 39: 104–118.SH ChoiK. HeoHM ByunWJ AnWG Lu2011Identification of preferential target sites for human DNA methyltransferases.Nucleic Acids Research39104118

30. Dejeux E, Rønneberg JA, Solvang H, Bukholm I, Geisler S, et al. (2010) DNA methylation profiling in doxorubicin treated primary locally advanced breast tumours identifies novel genes associated with survival and treatment response. Mol Cancer 9: 68.E. DejeuxJA RønnebergH. SolvangI. BukholmS. Geisler2010DNA methylation profiling in doxorubicin treated primary locally advanced breast tumours identifies novel genes associated with survival and treatment response.Mol Cancer968

31. Corn PG, Kuerbitz SJ, van Noesel MM, Esteller M, Compitello N, et al. (1999) Transcriptional silencing of the p73 gene in acute lymphoblastic leukemia and Burkitt’s lymphoma is associated with 5′ CpG island methylation. Cancer Res 59: 3352–3356.PG CornSJ KuerbitzMM van NoeselM. EstellerN. Compitello1999Transcriptional silencing of the p73 gene in acute lymphoblastic leukemia and Burkitt’s lymphoma is associated with 5′ CpG island methylation.Cancer Res5933523356

32. Feinberg AP, Tycko B (2004) Timeline - The history of cancer epigenetics. Nature Reviews Cancer 4: 143–153.AP FeinbergB. Tycko2004Timeline - The history of cancer epigenetics.Nature Reviews Cancer4143153

33. Siu LLP, Chan JK, Wong KF, Choy C, Kwong YL (2003) Aberrant promoter CpG methylation as a molecular marker for disease monitoring in natural killer cell lymphomas. Br J Haematol 122: 70–77.LLP SiuJK ChanKF WongC. ChoyYL Kwong2003Aberrant promoter CpG methylation as a molecular marker for disease monitoring in natural killer cell lymphomas.Br J Haematol1227077

34. Pei JH, Luo SQ, Zhong Y, Chen JH, Xiao HW, et al. (2011) The association between non-Hodgkin lymphoma and methylation of p73. Tumor Biology 32: 1133–1138.JH PeiSQ LuoY. ZhongJH ChenHW Xiao2011The association between non-Hodgkin lymphoma and methylation of p73.Tumor Biology3211331138

35. Galm O, Wilop S, Luders C, Jost E, Gehbauer G, et al. (2005) Clinical implications of aberrant DNA methylation patterns in acute myelogenous leukemia. Annals of Hematology. pp. 39–46.O. GalmS. WilopC. LudersE. JostG. Gehbauer2005Clinical implications of aberrant DNA methylation patterns in acute myelogenous leukemia.Annals of Hematology3946

36. Roman-Gomez J, Agirre X, Jiménez-Velasco A, Arqueros V, Vilas-Zornoza A, et al. (2009) Epigenetic regulation of microRNAs in acute lymphoblastic leukemia. J Clin Oncol 27: 1316–1322.J. Roman-GomezX. AgirreA. Jiménez-VelascoV. ArquerosA. Vilas-Zornoza2009Epigenetic regulation of microRNAs in acute lymphoblastic leukemia.J Clin Oncol2713161322

37. Chiam K, Centenera MM, Butler LM, Tilley WD, Bianco-Miotto T (2011) GSTP1 DNA methylation and expression status is indicative of 5-aza-2′-deoxycytidine efficacy in human prostate cancer cells. PLoS One 6(9): e25634.K. ChiamMM CenteneraLM ButlerWD TilleyT. Bianco-Miotto2011GSTP1 DNA methylation and expression status is indicative of 5-aza-2′-deoxycytidine efficacy in human prostate cancer cells.PLoS One6(9)e25634

38. Schmelz K, Wagner M, Dörken B, Tamm I (2005) 5-Aza-2′-deoxycytidine induce p21(WAF) express by demethylation of p73 leading to p53-independent apoptosis in myeloid leukemia. Int J Cancer 114: 683–695.K. SchmelzM. WagnerB. DörkenI. Tamm20055-Aza-2′-deoxycytidine induce p21(WAF) express by demethylation of p73 leading to p53-independent apoptosis in myeloid leukemia.Int J Cancer114683695

39. Wijermans PW, Lübbert M, Verhoef G, Klimek V, Bosly A (2005) An epigenetic approach to the treatment of advanced MDS; the experience with the DNA demethylating agent 5-aza-2′-deoxycytidine (decitabine) in 177 patients. Ann Hematol. pp. 9–17.PW WijermansM. LübbertG. VerhoefV. KlimekA. Bosly2005An epigenetic approach to the treatment of advanced MDS; the experience with the DNA demethylating agent 5-aza-2′-deoxycytidine (decitabine) in 177 patients.Ann Hematol917

40. Yang AS, Doshi KD, Choi SW, Mason JB, Mannari RK, et al. (2006) DNA methylation changes after 5-aza-2′-deoxycytidine therapy in patients with leukemia. Cancer Res 66: 5495–5503.AS YangKD DoshiSW ChoiJB MasonRK Mannari2006DNA methylation changes after 5-aza-2′-deoxycytidine therapy in patients with leukemia.Cancer Res6654955503

41. Link PA, Baer MR, James SR, Jones DA, Karpf AR (2008) p53-inducible ribonucleotide reductase (p53R2/RRM2B) is a DNA hypomethylation-independent decitabine gene target that correlates with clinical response in myelodysplastic syndrome/acute myelogenous leukemia. Cancer Res 68: 9358–9366.PA LinkMR BaerSR JamesDA JonesAR Karpf2008p53-inducible ribonucleotide reductase (p53R2/RRM2B) is a DNA hypomethylation-independent decitabine gene target that correlates with clinical response in myelodysplastic syndrome/acute myelogenous leukemia.Cancer Res6893589366