In the celiac disease analysis, we found 5 regions (TAGAP, IL18R1, RGS21, PLEK, and CCR9) with SNPs that had p-values <10−4: rs1054091/TAGAP (HR = 1.59; p = 5.8x10-6); rs4851575/IL18R1 (HR = 1.45; p = 5.7x10-5); rs1936670/RGS21 (HR = 2.23; p = 5.1x10-5); rs114569351/PLEK (HR = 2.64; p = 4.2x10-5); and rs12493471/CCR9 (HR = 1.40; p = 6.4x10-5) (Table 1, S2 Fig). Here, RGS21 is linked to celiac disease.