The role of low-frequency variants in type 1 diabetes (T1D) susceptibility still remains to be clarified. In the present study, we analyzed low-frequency variants of the T1D candidate genes in Japanese. We first screened for protein-changing variants of 24 T1D candidate genes in 96 T1D patients and 96 control subjects, and then the association with T1D was tested in 706 T1D patients and 863 control subjects recruited from the collaborating institutions in Japan. In total, 56 protein-changing variants were discovered; among them, 34 were low-frequency variants (allele frequency <5%). The association analysis of the low-frequency variants revealed that only the A908V variant of GLIS3 was strongly associated with resistance to T1D (Haldane's odds ratio = 0.046, p = 8.21 x 10(-4), and p(c) = 2.22 x 10(-2)). GLIS3 is a zinc finger transcription factor that is highly expressed in pancreatic beta cells, and regulates beta cell development and insulin gene expression. GLIS3 mRNA is also moderately expressed in the human thymus. The precise mechanism responsible for the association is unclear at present, but the A908V variant may affect autoimmunity to the GLIS3 protein itself; the 908V containing epitope may induce central or peripheral tolerance more efficiently than that of 908A. (C) 2013 Elsevier Inc. All rights reserved.