How genetic variations in apoptosis pathway interact with environmental factors to contribute to esophageal adenocarcinoma (EA) risk has not been comprehensively investigated. We conducted a case-only analysis in 335 Caucasian EA patients that were genotyped for 242 single nucleotide polymorphisms (SNPs) in 43 apoptotic genes. Gene-environment interactions were assessed using a two-step approach. First, random forest algorithm was used to screen for the potential interacting markers. Next, we used case-only logistic regression model to estimate the effects of gene-environment interactions on EA risk. Four SNPs (PERP rs648802; PIK3CA rs4855094, rs7644468 and TNFRSF1A rs4149579) had significant interaction with gastroesophageal reflux disease (GERD). The presence of variant alleles in TP53BP1 rs560191, CASP7 rs7907519 or BCL2 rs12454712 enhanced the risk of smoking by 2.08-2.58 times [interaction odds ratio (ORi)=2.08-2.58, adjusted P-value (Padj)=0.02-0.04]. Compared with patients carrying ≤1 risk genotype, the risk of GERD on EA was increased in persons with two (ORi=1.89, Padj=0.016) or ≥3 (ORi=4.30, Padj<0.0001) risk genotypes. Compared with cases with ≤1 risk genotype, smoking-associated EA risk increased by 3.15 times when ≥2 risk genotypes were present (ORi=3.15, Padj<0.0001). In conclusion, interactions among apoptotic SNPs and GERD or smoking play an important role in EA development.