Programmed silencing of γ-globin genes in adult erythropoiesis is mediated by several chromatin remodeling complexes, which determine the stage-specific genome architecture in this region. Identification of cis- or trans-acting mutations contributing to the diverse extent of Hb F might illustrate the underlying mechanism of γ-β globin switching. Here, we recruit a cohort of 1142 β-thalassemia patients and dissect the natural variants in the whole β-globin gene cluster through a targeted next-generation sequencing panel. A previously unreported SNP rs7948668, predicted to disrupt the binding motif of IKAROS as a key component of chromatin remodeling complexes, is identified to be significantly associated with higher levels of Hb F and age at onset. Gene-editing on this SNP leads to elevation of Hb F in both HUDEP-2 and primary CD34+ cells while the extent of elevation is amplified in the context of β-thalassemia mutations, indicating epistasis effects of the SNP in the regulation of Hb F. Finally, we perform ChIP-qPCR and 4C assays to prove that this variant disrupts the binding motif of IKAROS, leading to enhanced competitiveness of HBG promoters to locus control regions. This study highlights the significance of common regulatory SNPs and provides potential targets for treating of β-hemoglobinopathy.