Menopausal hormone therapy (MHT) is associated with an elevated risk of breast cancer in postmenopausal women. To identify genetic loci that modify breast cancer risk related to MHT use in postmenopausal women, we conducted a two-stage genome-wide association study (GWAS) with replication. In stage I, we performed a case-only GWAS in 731 invasive breast cancer cases from the German case-control study Mammary Carcinoma Risk Factor Investigation (MARIE). The 1,200 single nucleotide polymorphisms (SNPs) showing the lowest P values for interaction with current MHT use (within 6 months prior to breast cancer diagnosis), were carried forward to stage II, involving pooled case-control analyses including additional MARIE subjects (1,375 cases, 1,974 controls) as well as 795 cases and 764 controls of a Swedish case-control study. A joint P value was calculated for a combined analysis of stages I and II. Replication of the most significant interaction of the combined stage I and II was performed using 5,795 cases and 5,390 controls from nine studies of the Breast Cancer Association Consortium (BCAC). The combined stage I and II yielded five SNPs on chromosomes 2, 7, and 18 with joint P values <6 × 10(-6) for effect modification of current MHT use. The most significant interaction was observed for rs6707272 (P = 3 × 10(-7)) on chromosome 2 but was not replicated in the BCAC studies (P = 0.21). The potentially modifying SNPs are in strong linkage disequilibrium with SNPs in TRIP12 and DNER on chromosome 2 and SETBP1 on chromosome 18, previously linked to carcinogenesis. However, none of the interaction effects reached genome-wide significance. The inability to replicate the top SNP × MHT interaction may be due to limited power of the replication phase. Our study, however, suggests that there are unlikely to be SNPs that interact strongly enough with MHT use to be clinically significant in European women.