Association and interaction of the IL4R, IL4, and IL13 loci with type 1 diabetes among Filipinos

Abstract

In the search for genes involved in type 1 diabetes (T1D), other than the well-established risk alleles at the human leukocyte antigen loci, we have investigated the association and interaction of polymorphisms in genes involved in the IL4/IL13 pathway in a sample of 90 Filipino patients with T1D and 94 controls. Ten single-nucleotide polymorphisms (SNPs), including two promoter SNPs in the IL4R locus on chromosome 16p11, one promoter SNP in the IL4 locus on chromosome 5q31, and four SNPs--including two promoter SNPs--in the IL13 locus on chromosome 5q31 were examined for association, linkage disequilibrium, and interaction. We found that both individual SNPs (IL4R L389L; odds ratio [OR] 0.34; 95% confidence interval [CI] 0.17-0.67; P=.001) and specific haplotypes both in IL4R (OR 0.10; 95% CI 0-0.5; P=.001) and for the five linked IL4 and IL13 SNPs (OR 3.47; P=.004) were strongly associated with susceptibility to T1D. Since IL4 and IL13 both serve as ligands for a receptor composed, in part, of the IL4R alpha chain, we looked for potential epistasis between polymorphisms in the IL4R locus on chromosome 16p11 and the five SNPs in the IL4 and IL13 loci on chromosome 5q31 and found, through use of a logistic-regression model, significant gene-gene interactions (P=.045, corrected for multiple comparisons by permutation analysis). Our data suggest that the risk for T1D is determined, in part, by polymorphisms within the IL4R locus, including promoter and coding-sequence variants, and by specific combinations of genotypes at the IL4R and the IL4 and IL13 loci.

Figure 1

Gene-gene interaction. ORs and patient and control counts of specific IL4R and IL4-IL13 genotypes are shown. A, IL4R E375A (Glu encoded by the A allele, Ala by the C allele) with IL4 T-524C. B, IL4R Q551R (Gln encoded by the A allele, Arg by the G allele) with IL13 A-1512C. C, IL4R C-3223T and IL13 A-1512C. D, IL4R C-3223T and IL13 C-1112T. The P value for the association of each genotype combination is shown above each OR bar.