Celiac disease (CD), a gluten intolerance disorder, was implicated to have 57 genetic susceptibility loci for Europeans but not for culturally and geographically distinct ethnic populations like Saudi Arabian CD patients. Therefore, we genotyped Saudi CD patients and healthy controls for three polymorphisms, that is, Phe196Ser in IRAK1, Trp262Arg in SH2B3, and Met518Thr in MMEL1 genes. Single locus analysis identified that carriers of the 518 Thr/Thr (MMEL1) genotype conferred a 1.6-fold increased disease risk compared to the noncarriers (OR = 2.6; 95% CI: 1.22-5.54; P < 0.01). This significance persisted even under allelic (OR = 1.55; 95% CI: 1.05-2.28; P = 0.02) and additive (OR = 0.35; 95% CI: 0.17-0.71; P = 0.03) genetic models. However, frequencies for Trp262Arg (SH2B3) and Phe196Ser (IRAK1) polymorphisms were not significantly different between patients and controls. The overall best MDR model included Met518Thr and Trp262Arg polymorphisms, with a maximal testing accuracy of 64.1% and a maximal cross-validation consistency of 10 out of 10 (P = 0.0156). Allelic distribution of the 518 Thr/Thr polymorphism in MMEL1 primarily suggests its independent and synergistic contribution towards CD susceptibility among Saudi patients. Lack of significant association of IRAK and SH2B3 gene polymorphisms in Saudi patients but their association in European groups suggests the genetic heterogeneity of CD.