This study examined whether epistasis between single nucleotide polymorphisms (SNPs) within proprotein convertase subtilisin/kexin type 1 (PCSK1) and dopamine β-hydroxylase (DBH) genes is associated with premature ovarian failure (POF). One hundred twenty women with POF and 222 female controls were recruited for this study. To genotype SNPs within PCSK1 and DBH, we used a GoldenGate assay with VeraCode technology, which uses an allele-specific primer extension method. Two SNPs (rs155979 and rs3762986) within PCSK1 and one SNP (rs1611114) within DBH, which were located in the 5' flanking region, were involved in synergistic interactions. The C allele in the rs155979 SNP showed an increased risk of POF in a dominant model when AA genotype in the rs1611114 SNP was present (odds ratio, 3.60; 95% CI, 1.82-7.14; P = 0.00024), whereas the G allele in the rs1611114 SNP showed a reduced risk of POF in a dominant model when at least one C allele at the rs155979 SNP was present (odds ratio, 0.24; 95% CI, 0.11-0.51; P = 0.00018) or one G allele at the rs3762986 SNP was present (odds ratio, 0.33; 95% CI, 0.19-0.60; P = 0.00023). Epistases between SNPs within PCSK1 and DBH genes are significantly associated with susceptibility or resistance to POF.