| BACKGROUND AND OBJECTIVEβ-thalassemia are a group of hereditary blood disorders characterized by anomalies in the synthesis of the beta chains of hemoglobin mainly in tropical and subtropical areas worldwide, including South China. The carrier rate of P-thalassemia is 2.54% in Guangdong province and 6.43% in Guangxi province. When both parents are P-thalassemia carriers there is a 25% risk at each pregnancy of having children with homozygous or compound heterozygous p-thalassemia. The clinical severity of P-thalassemia varies widely from mild to severe forms and genotype-phenotype associations in β-thalassemia have been extensively studied by systematic analysis of defective β-globin genes. However, the same genotype of β-thalassemia can manifest with highly variable patterns of clinical severity ranging from moderate to severe disease due to various genetic modifiers of disease severity. Several genetic modifiers unlinked to the β-globin gene locus have been identified, including genetic variants modulating fetal hemoglobin (HbF) levels, a-globin gene copy numbers/a-thalassemia, alpha hemoglobin stabilizing protein (AHSP), and heme-regulated eIF2alpha kinase (HRI) in either human patients or murine models of the disease. Increased HbF levels and concomitant a-thalassemia are two main modifiers that can ameliorate the clinical and hematological severity of β-thalassemia.The erythroid transcriptional factor KLF1 has recently emerged as one of the key regulators of the γ- to β-globin gene switching. It is now known that the dominant mutation (p.E325K) in KLF1 can cause congenital dyserythropoietic anemia (CDA, type IV) and a compound heterozygous mutation in KLF1 lead to chronic hemolytic anemia but most monoallelic KLF1 mutations give rise to benign phenotypes, including significant increases in levels of HbF and adult hemoglobin A2 (HbA2). However, little is known about the modification of β-thalassemia phenotypes in the context of co-inheritance of KLF1 mutations with β-globin gene mutations. In addition, currently there is little information regarding the population genetics of KLF1 mutations. Here, we investigated the incidence and spectrum of KLF1 mutations in the Chinese population, and show that KLF1 mutations are common in a thalassemia endemic region in southern China and that the hematologic features of β-thalassemia are modulated by different KLF1 mutations.DESIGN AND METHODSPopulation, patients and phenotypic testsWe designed a large population-based cohort study that consisted of a total of 3918 randomly selected unrelated subjects living in southern China. This cohort included 3839 consecutive subjects who were investigated to determine the incidence of KLF1 mutations and 79 healthy subjects with borderline HbA2 (3.3%-4.1%) and/or elevated HbF levels (> 1.5%) to determine KLF1 mutation frequency in this selected population. All subjects in whom KLF1 mutations were identified were available for study of both the spectrum of mutations and of the hematological features of KLF1 variants. The cohort of 3839 consecutive subjects (aged 1—50 years; 1987 males and 1852 females) was grouped into three subgroups based on the β-globin genotypes: 1971 non-thalassemics (cohort A),946 β-thalassemia heterozygotes (cohort B), and 922 β-thalassemia homozygotes or compound heterozygotes (cohort C). Cohort A subjects were randomly drawn from 5789 samples used in a previous study. Consecutive samples in cohort B and cohort C were recruited from November 2008 to June 2013 from southern China. The subset of 79 subjects (cohort D) was also selected from a large random Guangxi population screened for the above two positive indices of suspected KLF1 mutations. All participants are ethnic southern Chinese, whose registered parental origin is from either Guangxi or Guangdong Provinces in southern China where thalassemia is highly prevalent. In addition, we recruited 1190 healthy Chinese subjects from Shandong Province in northern China where thalassemia is much less prevalent. Clinical data was collected retrospectively on all 922 thalassemia patients in cohort C. Diagnosis of thalassemia was made through a comprehensive assessment of clinical presentation, hematological and molecular analysis based on criteria previously described. A summary of clinical and genetic characterizations of 922 patients,568 thalassemia major (TM) and 354 thalassemia intermedia (TI) in cohort C, is presented in the text.Hematological parameters were measured using an automated hematology analyzer (Sysmex F-820, Sysmex Co. Ltd, Kobe, Japan) and hemoglobin analysis was carried out using the high-performance liquid chromatography (Variant Ⅱ, Bio-Rad Laboratories, Hercules, CA, USA). Flow-cytometric analysis of CD44 expression on mature erythrocytes was performed with a FACSCanto Ⅱ flow cytometer (BD Biosciences, San Jose, CA, USA) by immunolabeling of red blood cells (RBCs) with reagents from BD Biosciences-Pharmingen (BD Biosciences) according to standard protocols. The study protocol was reviewed and approved by the local medical ethics committee at Nanfang Hospital, China, and written consent was obtained from parents and/or patients prior to commencing the study. The study was conducted in accordance with the Declaration of Helsinki.Molecular genetic testsGenomic DNA was extracted from peripheral blood using the standard phenol/chloroform method. We used the high-resolution melting (HRM) assay to scan for KLF1 mutations and followed by Sanger sequencing for confirmation of KLF1 variants. A total of nine pairs of primers were designed to amplify the coding region (exons 1-3) and the conserved non-coding regions (promoter region, flanking splice junctions, and 5’- and 3’-untranslated regions) of the KLF1 gene. XmnI PCR-RFLP and HRM assays were used to genotype four primary HbF-associated SNPs (rs7482144, XmnI, in HBG2, rs11886868 and rs766432 in BCL11A, and rs9399137 in HBS1L-MYB. Point mutations and deletions in the α-, β-, and γ-globin genes were detected using previously described methods.Functionally relevant KLF1 mutations were identified on the basis of phenotypic changes (decreased erythrocytes CD44 expression and/or HbF levels> 1%) and functional predictions for non-synonymous substitutions in KLF1 analyzed by SIFT (http://sift.jcvi.org) and PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2). To determine whether the mutations are hereditary or somatic, peripheral blood and DNA samples were available from six parents whose offspring were identified as carriers for KLF1 mutations in cohort C.Statistical analysisStatistical analyses were carried out using SPSS version 19.0 (IBM, Armonk, NY, USA) and P< 0.05 was considered statistically significant. All quantitative variables were compared by either the Mann-Whitney U test or the χ2 test as appropriate. Multivariate analysis was performed using a backward stepwise Cox proportional hazards model in 922 β-thalassemia homozygotes or compound heterozygotes in cohort C, and hazard ratios (HR) with 95% confidence intervals (CI) were obtained to evaluate the associations between putative ameliorating alleles and the patients’ age at first transfusion. All covariates were classified as follows:α-globin gene mutations as 0,1,2 or 3 based on the number of mutated copies of the HBA; For other loci, each variable was defined with values 0,1 or 2 according to the number of copies of putative modifying allele:HBB mutations (β+), rs7482144 (XmnI:+) in HBG2, rs11886868 (C) and rs766432 (C) in BCL11A, rs9399137 (C) in HBS1L-MYB; KLF1 mutations was only scored as 0 or 1 since no homozygous mutation was noted; Gender for male was 1 and for female was 0. The performance of the model was measured by Harrell’s concordance index (C-index) by using the package of rms in R version 3.0.1 (http://www.r-project.org/). Univariate analysis was conducted to compare clinical features by either the Mann-Whitney U test or the X2 test as appropriate in β-thalassemia patients with or without KLFl mutations, and odds ratios (OR) with 95% CI were calculated to assess the risk of receiving systematic transfusion and TI. These two groups had same globin genotype category (β0/β0 and αα/αα) and similar HbF-modulating genetic variants [XmnI (-/-), rs766432 (AA or AC), rs9399137 (TT or CT)]:7 out of 12 KLFl mutation-positive subjects and 362 out of 910 KLF1 mutation-negative subjects in cohort C were included in both case-control and survival analysis. Transfusion-free survival curves were depicted according to the Kaplan-Meier method and compared using the log-rank test.RESULTSThe prevalence of KLF1 mutations is significantly higher in the thalassemia endemic region in southern China.It has been recently reported that KLFl mutations result in elevated HbF levels. As elevated HbF is beneficial in β-thalassemia, we hypothesized that some associations may exist between KLF1 mutations and β-thalassemia in regions with high prevalence of thalassemia. To test this hypothesis, we scanned for mutations of the KLF1 gene in two distinct Chinese populations:a total of 3839 unrelated individuals from three different cohorts in southern China where thalassemia is endemic and a total of 1190 individuals from northern China where thalassemia rarely occurs. We identified 48 subjects with KLF1 heterozygous mutations in the cohort from southern China, giving a heterozygote frequency of 1.25% in this population. In contrast, we identified only one subject with heterozygous missense mutations in the population from north, giving a 0.08% frequency which is significantly lower than found in the population from south (p< 0.001). These findings demonstrate that KLFl mutations have a much higher prevalence in the presence of thalassemia. Comparison of the frequency of KLF1 mutations among the three cohorts from south revealed that the frequency is 1.27%(25/1971) in cohort A,1.16%(11/946) in cohort B and 1.30% (12/922) in cohort C respectively. Statistical analysis showed that there were no statistically significant differences among these three sample groups (p= 0.959), implying that higher KLF1 mutation prevalence is a common phenomenon of a southern Chinese population where thalassemia is endemic.The spectrum of KLF1 mutations in southern China.As described above, we have identified 48 heterozygotes with KLF1 mutations in the southern Chinese population comprising 3839 subjects. To enlarge the size of sample with KLF1 mutations for further analysis, we identified additional 16 subjects with KLF1 heterozygous mutations in 79 subjects (non-thalassemics, cohort D) based on borderline increases in HbA2 (3.3%-4.1%) and/or elevated HbF levels (> 1.5%). Thus a total of 64 mutant alleles were documented in the present study, of which 41 are KLF1 heterozygotes alone,11 are KLFl mutations co-inherited with P-thalassemia heterozygotes, and 12 are KLF1 mutations co-inherited with β-thalassemia homozygotes or compound heterozygotes. We identified four previously reported KLF1 gene mutations (p.Glyl76AlafsX179, p.Ala298Pro, p.Thr334Arg and c.913+1G>A) and three novel mutations. Among these are five missense mutations (p.Glu5Lys, p.Ala298Pro, p.His299Asp, p.Thr334Arg and p.Cys341Tyr), one frameshift mutation (p.Gly176AlafsX179), one splicing mutation (c.913+1G>A) and the number and percentage of mutant alleles that carried these seven different types of KLF1 mutations are listed in the text. We determined that the six mutations (p.Gly176AlafsX179 frameshift, p.Ala298Pro, p.His299Asp, p.Thr334Arg, p.Cys341Tyr and c.913+1G>A) disrupting the zinc-finger (ZF) domain of KLF1 accounted for the vast majority of the molecular changes with two mutations (p.Gly176AlafsX179 and p.His299Asp) being most frequent, accounting for 90.6% of the KLF1 mutant alleles in this population. In addition, we found 13 novel neutral KLF1 variants including 4 missense (tolerated) and 4 synonymous mutations in the coding region and 5 variants in the 3’-UTR in our population. Our findings, together with previously published studies demonstrate the broad spectrum of KLF1 mutations.The effects of KLF1 mutations on hematologic parameters of non-thalassemic subjects.Having demonstrated that KLF1 mutations occur at about the similar frequency in all three cohorts (non-thalassemics, β-thalassemia heterozygotes and β-thalassemia homozygotes/compound heterozygotes), we then examined the effect of these mutations on the hematological parameters of these three cohorts. MCV and MCH values of normal individuals with or without KLF1 mutations. Both MCV and MCH were mildly lower in 41 normal subjects with KLF1 mutations (25 from cohort A and 16 from cohort D), compared to the normal subjects without KLF1 mutations (n= 1946, P< 0.001), demonstrating that mutations of KLF1 result in MCV and MCH values at the lower end of the normal range. In normal subjects without KLF1 mutations, MCV and MCH values were 89.29 ± 4.77 fL (n= 1946, range,79.00-110.10 fL) and 29.46 ± 1.44 pg (n= 1946, range,25.00-36.30 pg), respectively. In contrast, MCV and MCH values were 82.45 ± 5.39 fL (n=41, range,67.50-91.00 fL) and 26.95 ± 2.11 pg (n= 41, range,21.20-30.00 pg), respectively, in normal subjects with the KLF1 mutations. Furthermore, we confirmed borderline increases in HbA2 levels (n= 41,3.28 ± 0.45%, range,2.50-4.87%) and of slightly elevated HbF levels (n= 41,3.44 ± 1.75%, range,1.00-7.40%) in KLFl heterozygous mutants in our normal population, in the two groups (25 from cohort A and 16 from cohort D) that expected to show the same effects of elevated HbA2 and/or HbF.The effects of KLF1 mutations on hematological parameters of β-thalassemia heterozygote subjects.We then examined the effects of KLF1 mutations in modulating hematological parameters of β-thalassemia heterozygotes. For this, we analyzed β0-thalassaemia heterozygotes with or without KLF1 mutations. MCV and MCH values were 63.53 ± 3.94 fL (n= 14, range,57.00-72.00 fL) and 20.26 ± 1.42 pg (n= 14, range,18.20-22.90 pg), for 14 β0-thalassaemia heterozygotes (β0/βN and αα/αα genotypes) with the KLF1 mutations (seven out of 11 from cohort B except for one βE/βN, one β+/βN and two β0/βN+ —/αα and seven from familial cases), while they were 65.11 ± 4.34 fL (n= 217, range,51.00-77.80 fL) and 20.34 ± 1.37 pg (n= 217, range,14.60-27.60 pg), for 217 β0-thalassemia carriers with no co-inherited KLF1 mutations. Surprisingly, in contrast to normal subjects no statistically significant differences between these two groups with respect to these two hematological parameters (P= 0.112 or P= 0.567) were noted. This suggests that mutations of KLF1 gene in the context of β-thalassemia trait does not have significant effect on MCV and MCH. However, we found that co-existence of β-thalassemia trait and KLF1 mutations results in higher levels of HbA2, of over 6.5%(n= 17,6.47 ± 0.81%, range,5.10-8.00%), implying additive effects of these two factors on 5-globin gene expression. In addition, a highly variable increase in HbF levels (n= 17,7.67 ± 4.11%, range,3.60-18.20%) was noted in individuals with co-inheritance of KLF1 mutations in the context of P-thalassemia trait. The basis of reason of this variability needs to be further studied.The effects of KLF1 mutations on P-thalassemia homozygotes/compound heterozygote individuals.To define the role of KLF1 mutations in modulating the clinical and hematological severity, we studied 922 β-thalassemia patients with known genotypes of β-globin and modifier genes including KLF1 mutations. We found it surprising that two KLF1 mutations were selectively represented in 12 TI patients who are all β-thalassemia homozygotes or compound heterozygotes in cohort C:Ten subjects had one frameshift mutation (p.Gly176AlafsX179) and the other two had the missense mutation (p.His299Asp). Both mutations disrupt the KLF1 ZF motif. The clinical features of all 12 affected β-thalassemia individuals are presented in the test. We assessed the influence of different combinations of genetic factors on the entire cohort by multivariate Cox regression analysis, measured by age at first transfusion, and the results showed that the KLF1 mutations had the strongest effect on the clinical severity of β-thalassemia phenotype (HR= 0.213, P< 0.001), followed by β+ mutation (HR= 0.353, P< 0.001), rs7482114 (XmnI,+) in HBG2 (HR= 0.512, P< 0.001), HBA mutations (HR= 0.671, P< 0.001), rs9399137 (C) in HBS1L-MYB (HR = 0.724, P< 0.001) and rs766432 (C) in BCL11A (HR= 0.761, P< 0.001).In these 12 TI patients, seven subjects were β0/β1nd αα/αα genotypes (patients P3, P4, P5, P6, P7, P11 and P12) and three subjects (patients P1, P2 and P8) were (β0/β+ and αα/αα genotypes and two subjects (patients P9 and P10) were β0/β0 co-inherited with the Southeast Asian deletion (SEA) and 3.7kb α-thalassemia deletion respectively. We further tested the effects of KLF1 mutations in modulating the clinical severity of p-thalassemia by univariate analysis of a matched case-control of the above seven subjects (β0/β0 and αα/αα genotypes with KLF1 mutations) and 362 thalassemic subjects without KLF1 mutation in cohort C. Univariate analysis confirmed the links between the KLF1 mutations and decreased clinical severity with remarkably increased hemoglobin concentration (P= 0.006) and HbF levels (P< 0.001) for TI cases with KLF1 mutations, along with delay in mean age at onset of anemia (P= 0.001), decreased transfusion requirements (OR= 0.082, P= 0.007), and a lower risk of exhibiting the features of TM (OR= 0.009, P< 0.001). This difference in HbF levels between 7 vs.362 matched samples is consistent with the results of the HbF levels between the whole TI vs. TM populations, implying a clinically homogeneous group of TI patients in our cohort C whether they had a KLF1 mutation or not. Kaplan-Meier survival curves showed significant differences between the two groups (log-rank test P< 0.001). The median time to the first transfusion in Control Group (362 without KLF1 mutations) was 6 months, compared to a highly significant improvement in transfusion-free survival observed in Case Group (7 with KLF1 mutations) with the median time to the first transfusion at 36 months. Likewise,342 KLF1 mutation-negative TI group received their first transfusion at 24 months, whereas it was initiated at 48 months in 12 KLF1 mutation-positive TI group, although there was no statistically significant difference between two TI subgroups (log-rank test P=0.093).Prevalence of KLF1 mutations is high in individuals with increased HbA2 and HbF levels.To further establish the link between KLF1 mutations and HbA2 and HbF levels, we studied 79 non-thalassemic subjects (cohort D) with elevated levels of HbA2 and HbF. In this sub-group, we found that the prevalence of KLF1 mutations to be 20.25% (16/79), suggesting elevation of HbA2 and HbF levels can potentially be used for screening of KLF1 mutations in the population.DISCUSSIONWe studied the prevalence of KLF1 mutations through a population-based molecular analysis of a large number of samples comprising 3839 participants in southern China (thalassemia endemic region) and 1190 normal individuals in northern China (non-thalassemia endemic region) using a HRM assay we developed. We documented a surprisingly high frequency of KLF1 heterozygous mutations (1.25%) in the population from southern China while it was only 0.08% in northern China. These findings demonstrate that mutations in the KLF1 gene are common in the thalassemia endemic region in southern China. Interestingly, while no functionally relevant KLF1 mutations in the ZF domain were identified in the northern Chinese population such mutations occurred at a very high frequency in southern China, accounting for 98.4% of the KLF1 mutant alleles.We further documented that KLF1 mutation is a genetic modifier ameliorating the clinical severity of β-thalassemia. Moreover, we observed the presence of hypochromic and microcytic red cells (thalassemia-like) or borderline MCV and MCH values in most non-thalassemic individuals as a result of KLF1 mutations. A recent study confirmed that compound heterozygotes for KLF1 mutations are able to cause hypochromic and microcytic anemia. Thus, we hypothesize that the frequent occurrence of the KLF1 mutations that affect protein function (particularly those disrupted the ZF domain of KLF1) might arise due to natural selection, such as in the presence of malaria in the tropical and sub-tropical areas. Our current results imply that KLF1 heterozygous mutations may be one important genetic modifier associated with the inherited hemoglobin disorders, particularly in β-thalassemia, among the southern Chinese population and potentially other ethnic populations.We have identified 64 KLF1 mutant alleles with seven different mutations. The detailed genetic analysis of blood samples from six available families showed germ line origin of KLF1 mutations. The frameshift mutation, p.Gly176AlafsX179, was previously found in Thailand, Vietnamese and Korean populations, and its finding in Chinese populations by us and others suggests that this mutation appears to be predominant in East Asian populations. In our study, only functionally effective KLF1 mutations were scored as positive. We typed all 64 subjects with KLF1 mutations according to the following two major predictors of altered function of KLF1:(1) The presence of reduced CD44 expression on erythrocytes by flow-cytometric analysis, and this was observed in the representative samples from all six types of mutations studied; (2) Determination of elevated HbF levels (>1%), which was the case for all positive samples. All KLF1 variants identified in our study were submitted to the dbSNP (http://www.ncbi.nlm.nih.gov/snp). Decreased CD44 expression due to ZF domain mutations or frameshift mutations disrupting ZF motifs of human KLF1 gene has been previously documented by several studies. Thus, use of a molecular assay to scan for KLF1 mutation-positive in combination with a flow-cytometric analysis for determining decreased CD44 expression will be a valuable tool for establishing functionally relevant KLF1 mutations. In addition to this, the In(Lu) phenotype had been previously used in British and France population to screen for KLF1 mutations to obtain reliable prevalence figures. By ELISA assay using labeled anti-ζ monoclonal antibody, elevated ζ-globin chains expression were found in four of the seven available TI patients (P4, P6, P7 and P9) with KLF1 heterozygous mutations, which has recently been reported in a transfusion-dependent hemolytic anemia caused by compound heterozygosity for KLF1 mutations.It is well known that the clinical manifestations of thalassemia can be modified by several genetic determinants including genetic variants that modulate HbF levels, β-globin mutations and a-globin gene copy numbers. Based on previous studies, which showed that the erythroid transcriptional factor KLF1 was a key regulator of the γ- to β-globin gene switching and KLF1 haploinsufficiency correlated with HPFH caused by decreased BCL11A expression, we decided to investigate the role of KLF1 mutations in modulating the clinical and hematological severity of thalassemia. To the best of our knowledge, this is the first study reporting the clinical features of 12 β-thalassemia patients co-inherited with KLF1 mutations, which had a mean Hb concentration of 82.57±7.81 g/L and a mean HbF level of 39.58±21.80 g/L and indicative of the amelioration of β-thalassemia severity by KLF1 mutations. Multivariate Cox regression analysis showed that the modifier of KLF1 mutations had an effect on the clinical severity of β-thalassemia phenotype and univariate analysis using a matched case-control and transfusion-free survival curve supported results of the multivariate analysis. In our case, two KLF1 mutations, p.Gly176AlafsX179 and p.His299Asp, were both verified to be functional variants of KLF1, which result from generating a truncated protein with loss of the three ZF motifs downstream by frameshift mutation or abnormal structural changes in the ZF1 domain of KLF1 by missense mutation. The clinical features of our KLF1 mutation-positive patients might be ascribed to suppression of BCL11A expression and elevated y-globin levels due to altered DNA-binding domain of KLF1 and impaired transactivation capacity of KLF1, because a similar result due to KLF1 mutations have been previously described. Furthermore, our functional studies in K562 cells using p.His299Asp and p.Gly176AlafsX179 mutant KLF1 constructs support such a possibility (data not shown). In addition, it needs to be noted that three subjects (patients P1, P2 and P8) with the β0/β+genotype and two (patients P9 and P10) who co-inherited with HBA deletions had a very mild, transfusion-independent phenotype highlighting that a combination of KLF1 mutations and known modifiers can greatly modify the clinical phenotype. In conclusion, we have defined a novel and important genetic determinant, KLF1 ZF heterozygous mutations, that contribute to TI. This knowledge should help in clinical diagnoses and genetic counseling, as well as enable designing of appropriate and personalized transfusion program for thalassemia patients with KLF1 mutations. The use of decreased CD44 expression as a surrogate for KLF1 mutations could be a clinically important tool to evaluate individuals with unexplained TI. We suggest that it is important to test for KLF1 mutations to ensure accurate clinical assessment in those TI patients with variable HbF increases since β-thalassemia which benefits from increased HbF can be caused by various genetic factors.This study provides an overview of the important hematological features of individuals with KLF1 heterozygous mutations. Our finding of borderline MCV and MCH values in normal carriers with KLF1 mutations validates previous findings of borderline HbA2 levels and elevated HbF levels resulting from KLF1 mutations. Furthermore, we noted that a heterozygote KLF1 mutation co-inherited with β-globin mutations can lead to a remarkable increase in both HbA2 and HbF levels, compared to those with only KLF1 heterozygous mutations. These two phenotypic outcomes are most likely a result of the cumulative effect of the heterozygous β-thalassemia and KLFl mutations. At the phenotypic level with respect to HbA2 and HbF, there are no differences among subjects with different KLF1 and β-globin mutations. Of interest, in the light of re-analysis of two previously published data, we found that ten KLF1 target genes associated with red blood cell phenotypes might underlie the decreased MCV and MCH values. This could explain at least in part our findings in normal individuals with loss-of-function mutation in KLF1. However, the underlying mechanisms remain to be further investigated. These findings are likely to be useful for the accurate determination of β-thalassemia carriers in screening programs. |