Two-locus haplotype analysis identified DRB1*03:01∼DQB1*02:01 [0.44 vs. 0.18, OR (95% CI) = 3.44 (2.33–5.1), Pc = 3.48 × 10−10]; DRB1*04:02∼DQB1*03:02 [0.077 vs. 0.014, OR = 6.06 (2.03–24.37), Pc = 2.3 × 10−3] and DRB1*04:05∼DQB1*03:02 [0.060 vs. 0.010, OR = 6.24 (1.79–33.34), Pc = 0.011] as positively associated, and DRB1*16:02∼DQB1*05:02 [0.024 vs. 0.075, OR = 0.3 (0.11–0.74), Pc = 0.041] as negatively associated with T1D, after applying Bonferroni correction. Here, HLA-DRB1 is linked to type 1 diabetes mellitus.