- Open Access
DNA methylation patterns of β-globin cluster in β-thalassemia patients
Clinical Epigenetics volume 12, Article number: 187 (2020)
Reactivation of fetal hemoglobin (HbF, α2γ2) holds a therapeutic target for β-thalassemia and sickle cell disease. Although many HbF regulators have been identified, the methylation patterns in β-globin cluster driving the fetal-to-adult hemoglobin switch remains to be determined.
Here, we evaluated DNA methylation patterns of the β-globin cluster from peripheral bloods of 105 β0/β0 thalassemia patients and 44 normal controls. We also recruited 15 bone marrows and 4 cord blood samples for further evaluation. We identified that the CpG sites in the locus control region (LCR) DNase I hypersensitive site 4 and 3 (HS4-3) regions, and γ- and β-globin promoters displayed hypomethylation in β0/β0-thalassemia patients, especially for the patients with high HbF level, as compared with normal controls. Furthermore, hypomethylations in most of CpG sites of the HS4-3 core regions were also observed in bone marrows (BM) of β0/β0-patients compared with normal controls; and methylation level of γ-globin promoter -50 and + 17 CpG sites showed lower methylation level in patients with high HbF level compared with those with low HbF level and a negative correlation with HbF level among β0-thalassemia patients. Finally, γ-globin promoter + 17 and + 50 CpG sites also displayed significant hypomethylation in cord blood (CB) tissues compared with BM tissues from normal controls.
Our findings revealed methylation patterns in β-globin cluster associated with β0 thalassemia disease and γ-globin expression, contributed to understand the epigenetic modification in β0 thalassemia patients and provided candidate targets for the therapies of β-hemoglobinopathies.
The human β-globin cluster, spanning a 70-kb region, is composed of five genes (5′-HBE -HBG2 -HBG1 -HBD -HBB-3′; 5′-ε-γG-γA-δ-β-3′) and a distal regulatory element known as LCR , which plays an important role in regulating the expression of all the genes in the cluster. β-thalassemia is characterized by a quantitative defect in the synthesis of β-globin chains underlaid by a marked genotypic heterogeneity of β-globin gene mutations. Reactivating of HbF was reported to be a therapeutic target for β-thalassemia and SCD. Several genetic modulators [2,3,4,5,6,7] and cis-regulatory elements [6, 8,9,10,11,12] involved in the regulation of human fetal hemoglobin (HbF), and concomitant α-thalassemia have been identified as ameliorators of β-thalassemia . DNA methyltransferase (DNMT) inhibitor 5-azacytidine (5-aza) has been applied to hemoglobinopathy patients (mainly in SCD patients) through improvement of HbF expression in the 1980s and 1990s [14, 15], which supported that DNA methylation plays an important role in human globin switching. In 2007, Rodwell Mabaera et al. used three human fetal liver and three BM samples, respectively, to reveal the changes of DNA methylation in β- and γ-globin promoters during fetal-adult globin switching, and proved the negative correlation between promoters DNA methylation level and expression of globin genes in human development stage . Mei Hsu et al. have analyzed the methylation status of regions in the murine β-like globin locus in uncultured primitive and definitive erythroblasts and other cultured primary and transformed cell types . They found a 20-kb domain, extended from the region just past the LCR to before β-major and encompassed the embryonic genes εy, βh1, and βh0, is hypomethylated only in primitive erythroid cells. Even retrotransposons in this region are hypomethylated in primitive erythroid cells.
Although epigenetic-based underlying mechanism driving the fetal-to-adult hemoglobin switch have been studied and reported in mouse and human, the systematic scanning of methylation patterns in β-globin cluster remains to be determined, especially in β-thalassemia patients. Here, we uncovered the methylation patterns of β-globin cluster, and explored its roles in regulating HbF gene expression in PB, BM, or CB samples from human β-thalassemia patients or normal controls.
To evaluate the methylation status of the β-globin clusters in β-thalassemia patients, we first screened the methylation patterns of total 97 CpG sites on the β-globin cluster, including ε-, γ-, δ-, and β-globin promoters, the core sequences of HS4 and HS3 and the 5′ flanking sequences of the LCR region, the CpGs-rich region of the HS2-HS1 intergenic region, the CpG islands (CpGI) in the endogenous retrovirus 9 (ERV-9) located at upstream of HS5, and the predicted binding region of the transcription factor BCL11A on the gene cluster located about 3.5 kb upstream of δ-globin gene  (Fig. 1a), in PB of 105 β0/β0 thalassemia patients and 44 normal controls using bisulfite sequencing method as previously described . We observed that most of CpG sites around the LCR HS4-HS3 regions, and γ- and β-globin promoter displayed significant hypomethylation in β0/β0-thalassemia patients compared with normal controls (Fig. 1b, shown with stars under indicated sites), especially for the TFH patients (Fig. 1c). We didn’t find any significant differences in methylation level in the other examined regions among the three indicated groups (data not shown). Since DNA methylation can be associated with aging progression and significant differences in age between β-thalassemia patients and normal controls, we then included age as a covariate in analysis of covariance (ANCOVA) of methylation difference between β-thalassemia patients and normal controls and observed that hypomethylation in β-thalassemia patients remained significance in most of CpG sites around the LCR HS4 region, and γ- and β-globin promoters (Fig. 1b, shown with stars above indicated sites). We further divided the β0/β0 thalassemia patients into two subgroups: subjects with age < 5 years old (Fig. 1d) and subjects with 5–15 years old (Fig. 1e), and observed that the TFH patients still displayed significant lower methylation levels than the TFL patients, especially for the subjects with 5–15 years old (Fig. 1e).
Given that DNA in PB mainly derived from the lymphocyte cells, which will have some limitations for DNA methylation research with high tissue heterogeneity, we then examined DNA methylation in BM GYPA (Glycophorin A) positive cells and observed that significant hypomethylation at most of CpG sites in core regions of HS4-HS3 and the 5′ flanking sequences, especially in the HS4 region, in β0/β0-patients compared with the normal controls (Fig. 2a). Moreover, we observed that the TFH patients showed significant lower methylation level at the − 50 and + 17 CpG sites of γ-globin promoter than both the TFL patients and normal controls (Fig. 2b). We did not observe any significant differences in β-globin promoter between β0/β0-patients and normal controls. Furthermore, we employed bisulfite clone method to validate methylation pattern of γ-globin promoter in β0/β0-thalassemia patients and normal controls, and observed that hypomethylation at − 50, + 17, and + 50 CpG sites of γ-globin promoter in the TFH patients compared with the TFL patients and normal controls (Fig. 2c).
To further verify the relationship of the methylation status of γ-globin promoter with the HbF level, we performed Pearson correlation analysis of γ-promoter methylation level with HbF level in BM tissues from nine β0/β0-patients and observed that γ-globin promoter—50 (r2 = 0.49, p = 0.03) and + 17 (r2 = 0.91, p < 0.001) CpG sites displayed significantly negative correlations with the HbF level (Fig. 3). More importantly, we observed that γ-globin promoter + 17 and + 50 CpG sites displayed significant hypomethylation in CB tissues compared with BM tissues of normal controls, which is consistent with a higher HbF expression level in CB (HbF: 76.2 ± 7.9%; Fig. 4) than in BM (HbF: 0.3 ± 0.1%). Altogether, these results indicated that the methylation level of LCR HS4-3 regions was associated with thalassemia, and the methylation level of γ-globin promoter was associated with HbF expression.
In this study, we found that the hypomethylation in LCR regions, especially the HS4-3 core regions, was correlated with thalassemia; the methylation level of + 17 CpG site in γ-globin promoter was negatively correlated with HbF expression level, providing an important role of methylation in regulating globin expression in different developmental stage.
It is important to uncover the methylation patterns that account for phenotypic manifestation and HbF expression in β-thalassemia. A previously study has uncovered the methylation patterns in CD34+ primary cells from fetal liver and BM of healthy donors . They found that several differentially methylated CpGs are located near genes with unanticipated roles in red blood cell differentiation and proliferation, and a correlation between the γ-globin promoter DNA methylation and HbF levels in red blood cell progenitors, which was consistent with our results that observed in BM. Some researchers also discovered dynamic DNA methylation of γ-globin promoter in mouse uncultured primitive, definitive erythroblasts, and other cultured primary and transformed cell types, supporting a mechanistic role for DNA methylation in developmental regulation of globin genes . This dynamic DNA methylation pattern was further observed in our study that γ-globin promoter methylation was higher in BM than in CB tissue, elucidating the negative relationship between γ-globin promoter methylation and gene expression with developmental stage.
We further found that the methylation patterns in LCR regions from PB tissue were mostly hypomethylated in β0/β0-thalassemia patients, which was consistent with the case-controls difference patterns observed in BM tissue. However, we found the methylation patterns in γ- globin promoter showed difference between TFH and TFL group of β-thalassemia. The methylation level of most CpG sites in γ-globin promoter were hypomethylated, especially in + 17 and + 50 sites, in the BM tissues of TFH patients, compared with TFL patients, which was consistent with our previous report . The γ-globin promoter methylation, not LCR region DNA methylation, might be the main effect to γ-globin expression, and the methylation alteration in the γ-globin promoter may be a driving factor leading to remote regulation between globin promoter and LCR regions . The hypomethylated patterns in globin promoters might make it easier to access for some transcription factors mediating the looping with globin promoter. In fact, our previous study showed that a SNP rs368698783 was associated with the methylation level on + 17 CpG site in γ-globin promoter , by recruiting DNMT3A to the promoter and further regulating the γ-globin expression, supporting the important role of γ-globin promoter methylation in regulation of γ-globin expression.
Since only 15 BM subjects employed in this study to obtain methylation alterations in LCR regions associated with β0-thalassemia and in γ-globin promoter associated with HbF expression level, it will require larger samples from BM tissues to further verify these results. In addition, given the limitation of ethic approval, we have no evaluated the patterns in fetal liver, which would give a more comprehensive view of the methylation patterns in different developmental stage. Finally, effect of methylation alterations at γ-globin promoter + 17 CpG site on the HbF level should be examined by using site-specific methylation through dCas9-MQ1-sgRNA system in further study.
Our findings revealed methylation patterns in β-globin cluster associated with β0 thalassemia disease or γ-globin expression, contributed to understand the epigenetic modification in β0 thalassemia patients and provided candidate targets for the therapies of β-hemoglobinopathies.
We recruited 105 subjects with β0/β0-thalassemia from southern China. These 105 unrelated Chinese β-thalassemia patients with a large range of HbF expression levels were stratified into low- and high-HbF groups for this study: 52 cases of high HbF (TFH, HbF: 58.9 ± 19.8%) and 53 cases of low HbF (TFL, HbF: 2.0 ± 0.7%). A subset of 44 healthy adult controls (CON, HbF: 0.3 ± 0.2%) was selected from a large random Guangxi population screened as the previously mentioned  (Additional file 1: Tables S1 and S2). DNA extracted from peripheral blood lymphocyte (PB) of above 105 β0/β0 thalassemia patients and 44 healthy donors were used for DNA methylation analysis. Moreover, DNA extracted from GYPA (Glycophorin A) positive cells separated from bone marrows (BM) of nine β0/β0 thalassemia patients including six TFL β0/β0-patients (HbF: 7.46 ± 5.28 g/L) and three TFH β0/β0-patients (HbF: 95.39 ± 9.69 g/L) and six healthy adult controls (HbF: 0.37 ± 0.05 g/L), eight of which were also used in our previous study , and from CB (HbF: 76.2 ± 7.9%) of four healthy adult controls were used for examination of the methylation patterns in different tissues.
Hematological parameters were measured using an automated hematology analyzer (Sysmex F-820; Sysmex Co. Ltd., Kobe, JP). Adult and fetal hemoglobin analysis was carried out using high-performance liquid chromatography (Variant II, Bio-Rad Laboratories, Hercules, CA). DNA was extracted DNA by Phenol–chloroform method. We used XmnI PCR-restriction fragment length polymorphism (RFLP) and high-resolution melting (HRM) assays to genotype KLF1 mutations and three primary HbF-associated single nucleotide polymorphisms (rs7482144, XmnI, in HBG2, rs766432 in BCL11A, and rs9399137 in HBS1L-MYB) as described [4, 12] to unified genetic variants between TFH and TFL group. The study protocol was reviewed and approved by the local medical ethics committee at Nanfang Hospital, China, and conducted in accordance with the Declaration of Helsinki.
DNA methylation assay
We used bisulfite conversion method to examine DNA methylation as previously described . Bisulfite conversion of DNA was performed by denaturing 1 mg genomic DNA with 0.3 M NaOH at 42℃ for 20 min, followed by 95℃ for 3 min and 0℃ for 1 min, and incubating at pH 5.0 with sodium metabisulfite (2.0 M) and hydroquinone (0.5 mM) at 55℃ for 16 h in the dark overlaid with mineral oil. Modified DNA was purified with Promega Wizard DNA Clean Up System (Madison, WI, USA). The eluted DNA was incubated with NaOH (0.3 M) at 37℃ for 15 min and neutralized by 3 M NH4-acetate to pH 7.0. The neutralized DNA was precipitated by 75% ethanol and recovered in 20 ml of TE buffer. We used nested-PCR to amplify the target DNA regions and performed Sanger sequencing. Because unmethylated cytosine (C) was replaced by thymine (T) after PCR reaction, any presence of C at CpG dinucleotides on the sequence represented a methylated C allele in the bisulfite-modified DNA. The percentage methylation of CpG sites was determined based on the C and T allele peaks by means of BioEdit Sequence Alignment Editor v188.8.131.52 (Carlsbad, CA, USA). We further performed bisulfite-treated clone sequencing assay to verify DNA methylation patterns examined by Sanger sequencing. In brief, we cloned BS-PCR products of bisulfite-DNA into the pMD19T cloning vector (D104, TAKARA, Japan), and sequenced 22 clones to calculate the percentage of C as methylation% for each subject. The primers used in this study are listed in Additional file 1: Table S3.
A two-tailed student t-test from SPSS v20 software was used for comparisons between the two indicated groups in this study. ANCOVA analysis with age as a covariate was also used in case–control analysis of methylation difference. P values of < 0.05 were considered as statistically significant.
Availability of data and materials
All the data in this study were available in the figures and tables in the manuscript.
Sickle cell disease
Locus control region
Hypersensitive site 4 and 3
Endogenous retrovirus 9
- Glycophorin A:
Fragment length polymorphism
Levings PP, Bungert J. The human beta-globin locus control region. J Eur J Biochem. 2002;269:1589–99.
Sankaran VG, Menne TF, Xu J, et al. Human fetal hemoglobin expression is regulated by the developmental stage-specific repressor BCL11A. J Science. 2008;322:1839–42.
Zhou D, Liu K, Sun CW, et al. KLF1 regulates BCL11A expression and gamma- to beta-globin gene switching. J Nat Genet. 2010;42:742–4.
Liu D, Zhang X, Yu L, et al. KLF1 mutations are relatively more common in a thalassemia endemic region and ameliorate the severity of beta-thalassemia. J Blood. 2014;124:803–11.
Stadhouders R, Aktuna S, Thongjuea S, et al. HBS1L-MYB intergenic variants modulate fetal hemoglobin via long-range MYB enhancers. J J Clin Invest. 2014;124:1699–710.
Bauer DE, Kamran SC, Lessard S, et al. An erythroid enhancer of BCL11A subject to genetic variation determines fetal hemoglobin level. J Science. 2013;342:253–7.
Grevet JD, Lan XJ, Hamagami N, et al. Domain-focused CRISPR screen identifies HRI as a fetal hemoglobin regulator in human erythroid cells. J Science. 2018;361:285.
Lettre G, Sankaran VG, Bezerra MA, et al. DNA polymorphisms at the BCL11A, HBS1L-MYB, and beta-globin loci associate with fetal hemoglobin levels and pain crises in sickle cell disease. J Proc Natl Acad Sci USA. 2008;105:11869–74.
Sankaran VG, Xu J, Byron R, et al. A functional element necessary for fetal hemoglobin silencing. J N Engl J Med. 2011;365:807–14.
Farrell JJ, Sherva RM, Chen ZY, et al. A 3-bp deletion in the HBS1L-MYB intergenic region on chromosome 6q23 is associated with HbF expression. J Blood. 2011;117:4935–45.
Gilman JG, Huisman TH. DNA sequence variation associated with elevated fetal G gamma globin production. J Blood. 1985;66:783–7.
Chen DY, Zuo YJ, Zhang XH, et al. A genetic variant ameliorates beta-thalassemia severity by epigenetic-mediated elevation of human fetal hemoglobin expression. J Am J Hum Genet. 2017;101:130–8.
Mettananda S, Gibbons RJ, Higgs DR. α-Globin as a molecular target in the treatment of beta-thalassemia. J Blood. 2015;125:3694–701.
Ley TJ, DeSimone J, Anagnou NP, et al. 5-azacytidine selectively increases gamma-globin synthesis in a patient with beta+ thalassemia. J N Engl J Med. 1982;307:1469–75.
Lowrey CH, Nienhuis AW. Brief report: treatment with azacitidine of patients with end-stage beta-thalassemia. J N Engl J Med. 1993;329:845–8.
Mabaera R, Richardson CA, Johnson K, et al. Developmental- and differentiation-specific patterns of human gamma- and beta-globin promoter DNA methylation. J Blood. 2007;110:1343–52.
Hsu M, Mabaera R, Lowrey CH, et al. CpG hypomethylation in a large domain encompassing the embryonic beta-like globin genes in primitive erythrocytes. J Mol Cell Biol. 2007;27:5047–54.
Lessard S, Beaudoin M, Benkirane K, et al. Comparison of DNA methylation profiles in human fetal and adult red blood cell progenitors. J Genome Med. 2015;7:1.
Deng W, Rupon JW, Krivega I, et al. Reactivation of developmentally silenced globin genes by forced chromatin looping. J Cell. 2014;158:849–60.
Xiong F, Sun M, Zhang X, et al. Molecular epidemiological survey of haemoglobinopathies in the Guangxi Zhuang Autonomous Region of southern China. J Clin Genet. 2010;78:139–48.
We thank the patients for their willingness to participate in this study. Prof. Xiangmin Xu for valuable suggestion and Dr. Xinhua Zhang for helping to recruit the participates.
This study was financial supports by the Guangdong Science and Technology Foundation (Grant No. 2019B030316032).
Ethics approval and consent to participate
Informed consent was obtained from all the patients in this study. The study protocol was reviewed and approved by the local medical ethics committee at Nanfang Hospital, China, and conducted in accordance with the Declaration of Helsinki.
Consent for publication
The authors of the manuscript have read and agreed to the consent for publication.
The authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Bao, X., Zuo, Y., Chen, D. et al. DNA methylation patterns of β-globin cluster in β-thalassemia patients. Clin Epigenet 12, 187 (2020). https://doi.org/10.1186/s13148-020-00987-2
- β0 thalassemia patients
- β-globin cluster
- Fetal hemoglobin
- Methylation patterns
- Bone marrow
- Cord blood
- Developmental stage