We have investigated, in 282 multiplex Caucasian families (the Human Biological Data Interchange Repository), the association of type 1 diabetes with polymorphisms in the IL4R gene. IL4R encodes a subunit of the interleukin-4 receptor, a molecule critical to T-helper cell development. By genotyping eight different IL4R single-nucleotide polymorphisms (SNPs) and identifying haplotypes (complex alleles) in the multiplex type 1 diabetic families who were stratified for HLA genotype, we have observed significant evidence of linkage and association of the IL4R gene to type 1 diabetes. In particular, we have identified a specific haplotype that appears to be protective and observed that this protective effect is strongest among individuals not carrying the HLA DR3/DR4 genotype (which confers the strongest genetic risk for type 1 diabetes). These findings suggest an important role for the IL4R gene in immune-related disease susceptibility and illustrate the value of using multi-SNP haplotype information in association studies.
Type 1 diabetes, an early-onset autoimmune disease of glucose homeostasis, is caused by destruction of insulin-producing cells of the pancreas and involves a cellular T-helper response characteristic of the Th1 (cell-mediated) pathway (1, 2). One molecule critical to T-helper cell differentiation and commitment is the interleukin-4 (IL4) receptor heterodimer, and in particular, the subunit encoded by the IL4R gene located at 16p11.2–12.1. Many of the known single-nucleotide polymorphisms (SNPs) in IL4R change the protein primary structure (reviewed in 3). Changes to IL4R amino acid residues, particularly at Ile50Val, Ser478Pro, and Gln551Arg, are suggested to be functionally important and impact IL4 receptor signaling (ref. 4 and citations therein). IL4R SNPs have been associated with susceptibility to immune-related diseases, including (but not limited to) asthma and atopy (4–7). However, a consistent pattern of allele or genotype associations with specific disease-related phenotypes has not emerged. Furthermore, the genotyping of single IL4R SNP alleles appears insufficient to capture the genetic diversity of the IL4R locus (8–13), and this may in part explain the discrepancies cited. Here, by genotyping eight different IL4R SNPs and identifying haplotypes (complex alleles) in multiplex type 1 diabetic families stratified with respect to high-risk HLA genotypes, we have observed significant evidence of linkage and association of the IL4R gene to type 1 diabetes.
We genotyped 282 type 1 diabetic multiplex families at eight SNPs in the IL4R gene (Table 1) using a multiplex PCR assay coupled with sequence-specific hybridization to an immobilized probe array (Fig. 1). The 282 Human Biological Data Interchange (HBDI) families used in this study are nuclear families of Caucasian ethnicity, with unaffected parents and two ascertained affected siblings (14). The two-generation family structure and the absence of recombination among SNPs enabled us to deduce IL4R haplotypes for all individuals and assess transmission patterns. Herein, we will use the term “haplotype” to represent the ordered set of IL4R SNP alleles present on a single chromosome. The “haplotype” notation is representative of the primary sequence of the protein and is consistent with terminology used by others (8, 11–13).
There is substantial evidence that many genetic loci, as well as environmental factors, are involved in predisposition to type 1 diabetes (15). However, the genetic locus with the largest contribution by far to familial clustering is the HLA region on chromosome 6p21.3 (IDDM1). This region is thought to contribute at least 50% of the overall genetic risk (16). Multiple loci within the HLA region (e.g., HLA-DRB1, -DQA1, -DQB1, -DPB1, and -A), and specific HLA locus genotypes and haplotypes, all contribute to the genetic risk (15–18). Some HLA alleles are positively associated with type 1 diabetes (conferring increased risk), while others are negatively associated (conferring reduced risk). One heterozygous genotype (DRB1*0301-DQB1*0201/DRB1* 0401-DQB1*0302) confers a significantly increased risk, with an odds ratio (OR) of 20–30. The identification of type 1 diabetes risk loci outside the DR-DQ region in some cases has required stratification of DR3/4 and non-DR3/4 genotypes (17).
In our analyses of the IL4R data, we stratified the families into two groups: those in which either or both of the affected sibs are DR3/4 versus those in which neither of the sibs is DR3/4 (see research design and methods). Since the HLA DR3/4 genotype confers such a large genetic risk for type 1 diabetes, our approach was to look in families without the high-risk HLA DR-DQ genotypes so that other less strong linkage and association signals might be detected. Since below we present results together for the full dataset and for the two subgroups, we apply a conservative Bonferroni adjustment to the unadjusted P values to take account of these multiple comparisons.
The transmission disequilibrium test (TDT) (19) was performed on the IL4R genotype data to test for the association of the individual alleles of the eight IL4R SNPs to type 1 diabetes. The results of the SNP locus TDT are shown in Table 2. Nearly statistically significant evidence of excess transmission is only seen at SNP 6 (S478P) and only in “neither sib DR3/4” families (unadjusted P = 0.03, adjusted P = 0.10), although among all families, the data show a strong trend for SNP 6. The transmissions of the “ reference” allele (Table 1) of SNPs 3, 4, 5, and 7 are elevated to the same extent (∼58%) as for SNP 6. We note that the comparisons of the eight SNPs presented are not independent statistical tests, as the SNPs are in very strong linkage disequilibrium (LD) (Table 1).
Evidence for linkage of IL4R to type 1 diabetes was assessed by the haplotype sharing method (16), following inference of haplotypes and DR3/4 stratification (see research design and methods). The distribution of haplotype identity-by-descent (IBD) sharing is shown in the online appendix. A nearly statistically significant deviation from the expected distribution of IBD sharing (unadjusted P = 0.03, adjusted P = 0.08) was observed in the “neither sib DR3/4” subgroup (online appendix), providing evidence for linkage of the IL4R locus to type 1 diabetes, but only in that subgroup. A λs = 1.16 for IL4R can be calculated for the complete dataset (20).
Genome-wide linkage scans have identified many potential regions of linkage of type 1 diabetes in addition to the HLA region (IDDM1) and INS (IDDM2). Whereas evidence exists for a linked locus at 16q22-q24 (logarithm of odds [LOD] = 3.93) in U.K. and U.S. families, support for linkage has also been reported at 16p11-q12 (LOD = 1.74), where IL4R resides (21). Our IBD sharing data (online appendix) provide evidence consistent with those reports for linkage of IL4R with type 1 diabetes in the “neither sib DR3/4” group. One explanation for the lower LOD score near IL4R in those other genome-wide scans is that they may not have had sufficient power to detect loci with a λs < 1.3.
Association of IL4R haplotypes with type 1 diabetes was assessed using a modified TDT haplotype association analysis (Table 3). We evaluated the nine most frequent eight-locus haplotypes using a two-sided modified TDT for the null hypothesis of no association and no linkage (see research design and methods). The TDT analyses (Table 3) revealed significant evidence of undertransmission of a single eight-locus IL4R haplotype (labeled “21222221” or H-3) in the full cohort (unadjusted P = 0.007, adjusted P = 0.021) and the “neither sib DR3/4” group (unadjusted P = 0.005, adjusted P = 0.015). In the full cohort, the H-3 haplotype has an OR of 0.63 (unadjusted 95% CI = 0.44–0.89), and in the “neither sib DR3/4” group, the OR = 0.49 (unadjusted 95% CI = 0.28–0.81) (Table 3). The haplotype H-3 can thus be regarded as conferring protection to type 1 diabetes.
The ratio of transmissions to nontransmissions of H-3 to affected sibs is 28:44 from H-1/H-3 parents and 17:25 from H-2/H-3 parents, among unambiguously determined haplotypes (data not shown). These transmission proportions of ∼40% differ from the observed 50% transmission of H-1 or H-2 from H-1/H-2 parents expected under the null hypothesis (data not shown). The effect of H-3 is thus specific and not due to an apparent overtransmission of any other haplotype.
The use of multiplex families such as the HBDI cohort increases the power to detect a genetic effect precisely because sibships are more likely to share a genetic factor. Because of this ascertainment method, the observed association of IL4R in type 1 diabetes might apply only to a subset of the general population.
A recent report by Reimsnider et al. (7) investigated the association of IL4R and type 1 diabetes using the same HBDI family cohort (among others), but concluded that there was none. There are several possible explanations for the differences between the results of that report and ours. First, our IBD sharing results revealed marginal evidence of linkage only when DR3/4 stratification was applied. Second, in that previous study, only a single IL4R SNP (no. 7, Q551R) was genotyped (7), and it was not the SNP (no. 6, S478P) for which we detected the best evidence for a transmission bias (Table 2). More importantly, the observation of eight-SNP haplotypes increases the informativeness of the locus by converting biallelic SNPs into one multiallelic locus (Table 3).
The association of haplotype H-3 with type 1 diabetes may reflect the combined effect of the eight SNPs that compose the haplotype. Alternatively, the eight IL4R SNPs studied here may not capture all the variation defining the type 1 diabetes-protective haplotype, and other SNPs on or in LD with haplotype H-3 may be involved. We have found that haplotype H-3 is associated with protection from type 1 diabetes, with a similar OR, in a Filipino population (T. Bugawan, D.B.M., A.M.V., A. Panelo, C. Solfelix, P. Pozzilli, R. Buzetti, H.A.B., manuscript in preparation) and favor the hypothesis that this specific variant of the IL4R protein confers protection to type 1 diabetes. An alternative hypothesis, in which H-3 acts as a marker in LD with other causal variation(s), seems less likely because the LD patterns and haplotype distributions among Filipinos are different from those of the Caucasian HBDI population.
Our results suggest that a specific combination of the eight SNPs, rather than any individual SNP, is responsible for the observed protection from type 1 diabetes (Table 3). The evidence for the negative association of the H-3 haplotype with type 1 diabetes is stronger than that seen with any of the individual SNPs (33 vs. 42% transmission in the “neither DR3/4” group; unadjusted P = 0.005 vs. 0.03, respectively). The effect of individual SNPs in a haplotype seems to depend on the nature of the other SNPs, which is biologically plausible, as these SNPs encode allelic variants of the same IL4R protein. For example, haplotypes H-1 and -2 differ only at SNP 1 and have the same OR (Table 3); H-3 and -5 also differ only at SNP1, although H-5 does not appear to influence susceptibility to type 1 diabetes, whereas H-3 does. The role of SNP 1 in susceptibility to type 1 diabetes is therefore context dependent, highlighting the value of using haplotype data.
How might haplotype H-3 confer protection to type 1 diabetes given the known role of IL4R variation in Th2-related disease? Although numerous studies report associations of individual IL4R SNP alleles and genotypes with asthma and atopy, the results are not consistent. In some studies, the variant allele is associated with a disease-related phenotype, while in others it is the reference allele. If, in fact, the variant alleles (or haplotypes bearing them, such as H-3) are associated with a Th2 phenotype, then the protection to type 1 diabetes reported here may reflect a shift of the balance between Th1 and Th2 responses toward Th1. If, on the other hand, the reference alleles (or haplotypes bearing them) are associated with asthma and atopy, then this simple model of Th1 versus Th2 balance does not fit the data. Perhaps, therefore, variations in IL4R influence the regulation of both Th1 and Th2 responses, as elaborated further below.
The SNPs we have used to define IL4R haplotypes do not fully coincide with those used by others in assessing haplotypic association to asthma and atopy phenotypes. Nonetheless, what Ober et al. (11) called haplotype “222122” and Hackstein et al. (8) called IL4RA*03 both correspond to our haplotype H-3. IL4RA*03 was reported to be associated with lower total IgE levels in nonatopics than were more common haplotypes (equivalent to our haplotypes H-1, -2, and -4) (8). Also, haplotype “222122” was observed by TDT to be significantly undertransmitted to asthmatics in an isolated Hutterite population, although in that study the TDT was used as a test of linkage (11). These observations suggest that H-3 might be associated with a dampened T-helper cell response in either the Th1 or Th2 direction. In vitro and in vivo studies comparing the functions of the receptor molecules expressed from fully defined IL4R haplotypes will be instructive. Whatever the mechanism, allelic variants of the IL4R protein, one of which is encoded by haplotype H-3, are likely to be modifiers of immune disease susceptibility.
RESEARCH DESIGN AND METHODS
Subjects.
The family-based samples were provided as purified genomic DNA from the HBDI, a repository for cell lines from families affected with type 1 diabetes (14). All of the 282 HBDI families used in this study are nuclear families with unaffected parents and at least two affected siblings. In addition to the 564 offspring (2 sibs in each of 282 families in the affected sib pairs on which ascertainment was based), there were 26 other affected children and 270 unaffected offspring among these families. Further description of these samples and details of HLA class II genotyping can be found elsewhere (16,17).
Definition of DR3/4, stratification.
Families were stratified into two groups based on the HLA genotype of the two ascertained affected children. A child possesses the high-risk HLA-DR3/DR4 genotype if he/she carries a DRB1* 0301-DQB1*02 haplotype and a DR4 allele that is not DRB1*0403 (a protective allele) in a haplotype with DQB1*02 or DQB1*0302 (but not *0301 or *0304). In “either/both sib DR3/4” families (n = 135), one or two children possess the DR3/DR4 genotype. In “neither sib DR3/4” families (n = 147), neither child possesses the DR3/DR4 genotype.
IL4R genotyping
List of SNPs detected.
The eight SNPs in the human IL4R gene are listed in Table 1. The SNPs will subsequently be referred to by the SNP number listed in Table 1. The “reference” allele for each SNP is that present in GenBank accession X52425.1 and is denoted by a “1” in subsequent tables; the variant allele for each SNP is denoted by a “2.”
PCR primers and amplification.
Primers were purchased from Operon (Alameda, CA) and were modified at the 5′ phosphate by conjugation with biotin. Pairs of primers were designed and used to amplify the regions encompassing the eight SNPs. Other SNPs not relevant to this report were genotyped concurrently by coamplification within the same multiplex PCR and hybridization to an immobilized probe array (data not shown).
PCR amplification was performed in a total reaction volume of 50–100 μl containing the following reagents: 70 mmol/l KCl, 12 mmol/l Tris-HCl (pH 8.3 at 25°C), 3 mmol/l MgCl2, 0.2 ng/μl purified human genomic DNA, 0.2 μmol/l each primer, 200 μmol/l each of dATP, dTTP, dCTP, and dGTP, and 0.25 units/μl AmpliTaq Gold DNA polymerase. Amplification was carried out in a GeneAmp PCR System 9600 thermal cycler (PE Biosystems, Foster City, CA) using the following specific temperature cycling profile: an initial hold at 94°C for 12.5 min; 33 cycles of 95°C for 45 s, 61°C for 30 s, and 72°C for 45 s; and a final extension step of 72°C for 7 min.
Allele-specific detection by probe hybridization.
The detection of allelic variants at each SNP is performed by stringent hybridization of the biotinylated PCR products to the immobilized sequence-specific probes, followed by visualization using an enzyme-catalyzed color development step. The procedure is essentially as described in Cheng et al. (22), except that hybridization and washing were performed at 55°C. An example of linear array results is shown in Fig. 1. Haplotypes were determined as described below.
Inference of haplotypes.
We inferred IL4R haplotypes by a combination of two methods: 1) using the GeneHunter program (Falling Rain Genomics, Palo Alto, CA), which calculates haplotypes from genotype data from pedigrees; and 2) inspecting each HBDI family pedigree individually using the Cyrillic program (Cherwell Scientific Publishing, Palo Alto, CA) to resolve any ambiguous or unsupported haplotype assignments. Haplotypes could be confidently and unambiguously assigned in all but 6 of the 282 families; in the remainder, the phase of SNP 1 could not be determined.
Statistical analysis, methods, and algorithms.
Evidence for linkage of IL4R to type 1 diabetes was assessed by the haplotype sharing method. IBD values of parental IL4R haplotypes in the affected sibs could be determined unambiguously in 256 families. The locus-specific sibling relative risk was calculated using the formula λ s = 0.25/(proportion of families IBD = 0) (20). CIs were calculated using the standard error of the IBD proportions. The TDT was performed for individual SNPs by the method of Spielman (19). Association of IL4R haplotypes with type 1 diabetes was assessed using the modified TDT haplotype association analysis (described below).
LD among the eight IL4R SNPs was assessed using the pool of AFBAC haplotypes deduced for the entire cohort. The AFBAC pool represents an unbiased estimate of the population haplotype frequencies (23). We calculated the values of two distinct but related metrics for LD, namely D’ and Δ, by standard methodology.
Modified TDT haplotype association analysis.
We evaluated the nine most frequent eight-locus haplotypes using a two-sided modified TDT for the null hypothesis of no association and no linkage. As obtaining haplotypes for parents based on data from the children can introduce bias (24), we used a statistic that does not rely on this information. The statistic is the proportion of times a haplotype was transmitted among informative transmissions where it was available. When both parents and a child are heterozygous at n ≥ 1 loci, the statistic gives weight = (0.5)n to each haplotype arrangement consistent with the observed diplotypes. We used 2,000 randomly permuted datasets to assess significance (25). A child’s complementary diplotype is defined as the entire set of alleles remaining when the child’s observed diplotype is subtracted from the combined diplotypes of the two parents. The permutation procedure randomly and independently switches the roles of each child’s observed and complementary diplotypes to create 2,000 permuted datasets. This procedure is valid when phase is unknown, even if recombination is present. Empirical CIs were based on 2,000 bootstrap samples, created by resampling the observed families with replacement. Intervals for the proportion transmitted were converted to equivalent intervals for the reported ratio, which can be interpreted as an OR for observed segregation relative to Mendelian segregation. The haplotype analysis was performed using the S-PLUS software package (Insightful, Seattle, WA), and S-PLUS bootstrap functions were used for both resampling methods.
SNP . | Amino acid variation . | IL4R SNP position in X52425.1 . | IL4R SNP position in AC004525.1 . | SNP . | . | . | . | . | . | . | . | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | . | . | 1 . | 2 . | 3 . | 4 . | 5 . | 6 . | 7 . | 8 . | |||||||
1 | I50V | 398A>G | 94272T>C | — | 1.00 | 0.37 | 0.37 | 0.33 | 0.52 | 0.36 | 1.00 | |||||||
2 | N142N | 676C>T | 92548A>G | 0.35 | — | −1.00 | −1.00 | −1.00 | 0.41 | 0.38 | 1.00 | |||||||
3 | E375A | 1374A>C | 76608T>G | 0.16 | −0.15 | — | 1.00 | 1.00 | 0.97 | 1.00 | −1.00 | |||||||
4 | L389L | 1417G>T | 76565C>A | 0.16 | −0.15 | 1.00 | — | 1.00 | 0.97 | 1.00 | −1.00 | |||||||
5 | C406R | 1466T>C | 76516A>G | 0.13 | −0.14 | 0.96 | 0.96 | — | 1.00 | 1.00 | −1.00 | |||||||
6 | S478P | 1682T>C | 76300A>G | 0.26 | 0.28 | 0.79 | 0.79 | 0.78 | — | 1.00 | 1.00 | |||||||
7 | Q551R | 1902A>G | 76080T>C | 0.21 | 0.23 | 0.72 | 0.72 | 0.69 | 0.88 | — | 1.00 | |||||||
8 | S761P | 2531T>C | 75451A>G | 0.10 | 0.29 | −0.04 | −0.04 | −0.04 | 0.20 | 0.18 | — |
SNP . | Amino acid variation . | IL4R SNP position in X52425.1 . | IL4R SNP position in AC004525.1 . | SNP . | . | . | . | . | . | . | . | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | . | . | 1 . | 2 . | 3 . | 4 . | 5 . | 6 . | 7 . | 8 . | |||||||
1 | I50V | 398A>G | 94272T>C | — | 1.00 | 0.37 | 0.37 | 0.33 | 0.52 | 0.36 | 1.00 | |||||||
2 | N142N | 676C>T | 92548A>G | 0.35 | — | −1.00 | −1.00 | −1.00 | 0.41 | 0.38 | 1.00 | |||||||
3 | E375A | 1374A>C | 76608T>G | 0.16 | −0.15 | — | 1.00 | 1.00 | 0.97 | 1.00 | −1.00 | |||||||
4 | L389L | 1417G>T | 76565C>A | 0.16 | −0.15 | 1.00 | — | 1.00 | 0.97 | 1.00 | −1.00 | |||||||
5 | C406R | 1466T>C | 76516A>G | 0.13 | −0.14 | 0.96 | 0.96 | — | 1.00 | 1.00 | −1.00 | |||||||
6 | S478P | 1682T>C | 76300A>G | 0.26 | 0.28 | 0.79 | 0.79 | 0.78 | — | 1.00 | 1.00 | |||||||
7 | Q551R | 1902A>G | 76080T>C | 0.21 | 0.23 | 0.72 | 0.72 | 0.69 | 0.88 | — | 1.00 | |||||||
8 | S761P | 2531T>C | 75451A>G | 0.10 | 0.29 | −0.04 | −0.04 | −0.04 | 0.20 | 0.18 | — |
Information identifying the IL4R SNPs is presented in the first four columns of the table. Each SNP is described by its position in the GenBank cDNA reference accession sequence X52425.1, and in the reference genomic sequence for chromosome 16 BAC clone AC004525.1 (these two sequences are antiparallel). The amino acid variation numbering refers to that of the mature IL4R protein (4,11). The remaining eight columns present two measures of linkage disequilibrium (LD) between all pairs of IL4R SNPs used in this study. D’ values are shown in the upper-right triangle; Δ values in the lower-left triangle. Bold text indicates LD measures whose absolute values are near 1.0, indicating “complete” or “absolute” pairwise LD.
SNP . | All families (n = 282) . | . | . | . | Either/both sib DR3/4 (n = 135) . | . | . | . | Neither sib DR3/4 (n = 147) . | . | . | . | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | Number of informative transmission events . | % Transmission of reference allele to affecteds . | McNemar χ2 statistic . | P . | Number of informative transmission events . | % Transmission of reference allele to affecteds . | McNemar χ2 statistic . | P . | Number of informative transmission events . | % Transmission of reference allele to affecteds . | McNemar χ2 statistic . | P . | |||||||||
1 | 564 | 47.3% | 1.60 | 0.21 | 256 | 51.6% | 0.25 | 0.62 | 308 | 53.6% | 1.57 | 0.21 | |||||||||
2 | 188 | 54.8% | 1.72 | 0.19 | 88 | 52.3% | 0.18 | 0.67 | 100 | 57.0% | 1.96 | 0.16 | |||||||||
3 | 216 | 56.0% | 3.13 | 0.08 | 104 | 53.8% | 0.62 | 0.43 | 112 | 58.0% | 2.89 | 0.09 | |||||||||
4 | 216 | 56.0% | 3.13 | 0.08 | 104 | 53.8% | 0.62 | 0.43 | 112 | 58.0% | 2.89 | 0.09 | |||||||||
5 | 214 | 55.1% | 2.26 | 0.13 | 104 | 53.8% | 0.62 | 0.43 | 110 | 56.4% | 1.78 | 0.18 | |||||||||
6 | 314 | 55.4% | 3.68 | 0.06 | 138 | 52.2% | 0.26 | 0.61 | 176 | 58.0% | 4.45 | 0.03 | |||||||||
7 | 368 | 54.1% | 2.45 | 0.12 | 170 | 52.9% | 0.59 | 0.44 | 198 | 55.1% | 2.02 | 0.16 | |||||||||
8 | 24 | 50.0% | 0.00 | 1.00 | 14 | 42.9% | 0.29 | 0.59 | 10 | 60.0% | 0.40 | 0.53 |
SNP . | All families (n = 282) . | . | . | . | Either/both sib DR3/4 (n = 135) . | . | . | . | Neither sib DR3/4 (n = 147) . | . | . | . | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | Number of informative transmission events . | % Transmission of reference allele to affecteds . | McNemar χ2 statistic . | P . | Number of informative transmission events . | % Transmission of reference allele to affecteds . | McNemar χ2 statistic . | P . | Number of informative transmission events . | % Transmission of reference allele to affecteds . | McNemar χ2 statistic . | P . | |||||||||
1 | 564 | 47.3% | 1.60 | 0.21 | 256 | 51.6% | 0.25 | 0.62 | 308 | 53.6% | 1.57 | 0.21 | |||||||||
2 | 188 | 54.8% | 1.72 | 0.19 | 88 | 52.3% | 0.18 | 0.67 | 100 | 57.0% | 1.96 | 0.16 | |||||||||
3 | 216 | 56.0% | 3.13 | 0.08 | 104 | 53.8% | 0.62 | 0.43 | 112 | 58.0% | 2.89 | 0.09 | |||||||||
4 | 216 | 56.0% | 3.13 | 0.08 | 104 | 53.8% | 0.62 | 0.43 | 112 | 58.0% | 2.89 | 0.09 | |||||||||
5 | 214 | 55.1% | 2.26 | 0.13 | 104 | 53.8% | 0.62 | 0.43 | 110 | 56.4% | 1.78 | 0.18 | |||||||||
6 | 314 | 55.4% | 3.68 | 0.06 | 138 | 52.2% | 0.26 | 0.61 | 176 | 58.0% | 4.45 | 0.03 | |||||||||
7 | 368 | 54.1% | 2.45 | 0.12 | 170 | 52.9% | 0.59 | 0.44 | 198 | 55.1% | 2.02 | 0.16 | |||||||||
8 | 24 | 50.0% | 0.00 | 1.00 | 14 | 42.9% | 0.29 | 0.59 | 10 | 60.0% | 0.40 | 0.53 |
The table lists the observed transmissions of the reference allele (as defined in Table 1) at each of the eight IL4R SNP loci individually. The TDT analysis was performed on the complete HBDI cohort and after stratification for HLA-DR3/4 status, as described in the text. The number of informative transmissions for each SNP and the percent transmission are listed, as are the associated McNemar χ2 statistic and the significance level. P values are shown before adjustment.
Haplotype ID . | Haplotype . | All families . | . | . | . | Either/both sib DR3/4 . | . | . | . | Neither sib DR3/4 . | . | . | . | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | Number of informative transmission events . | % transmitted . | OR . | 95% CI . | Number of informative transmission events . | % transmitted . | OR . | 95% CI . | Number of informative transmission events . | % transmitted . | OR . | 95% CI . | |||||||||
H-1 | 11111111 | 550 | 52% | 1.09 | 0.92–1.30 | 269 | 52% | 1.08 | 0.85–1.38 | 281 | 53% | 1.11 | 0.86–1.40 | |||||||||
H-2 | 21111111 | 498 | 51% | 1.06 | 0.90–1.25 | 226 | 50% | 1.01 | 0.80–1.28 | 272 | 52% | 1.10 | 0.86–1.42 | |||||||||
H-3 | 21222221 | 122 | 39%* | 0.63 | 0.44–0.89 | 61 | 44% | 0.80 | 0.50–1.27 | 62 | 33%† | 0.49 | 0.28–0.81 | |||||||||
H-4 | 22111111 | 93 | 47% | 0.89 | 0.61–1.30 | 50 | 45% | 0.81 | 0.47–1.36 | 43 | 50% | 0.99 | 0.56–1.71 | |||||||||
H-5 | 11222221 | 86 | 49% | 0.95 | 0.59–1.50 | 37 | 40% | 0.68 | 0.30–1.42 | 49 | 55% | 1.24 | 0.67–2.36 | |||||||||
H-6 | 22111221 | 77 | 44% | 0.78 | 0.45–1.24 | 27 | 52% | 1.07 | 0.49–2.40 | 50 | 40% | 0.65 | 0.31–1.23 | |||||||||
H-7 | 11111121 | 77 | 57% | 1.33 | 0.90–2.06 | 43 | 55% | 1.23 | 0.71–2.20 | 33 | 60% | 1.48 | 0.75–3.19 | |||||||||
H-8 | 21111121 | 27 | 43% | 0.75 | 0.38–1.60 | 12 | 38% | 0.61 | 0.18–2.22 | 15 | 47% | 0.88 | 0.38–3.09 | |||||||||
H-9 | 22111222 | 16 | 45% | 0.83 | 0.25–2.68 | 10 | 50% | 1.00 | —‡ | 6 | 38% | 0.60 | —‡ |
Haplotype ID . | Haplotype . | All families . | . | . | . | Either/both sib DR3/4 . | . | . | . | Neither sib DR3/4 . | . | . | . | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | Number of informative transmission events . | % transmitted . | OR . | 95% CI . | Number of informative transmission events . | % transmitted . | OR . | 95% CI . | Number of informative transmission events . | % transmitted . | OR . | 95% CI . | |||||||||
H-1 | 11111111 | 550 | 52% | 1.09 | 0.92–1.30 | 269 | 52% | 1.08 | 0.85–1.38 | 281 | 53% | 1.11 | 0.86–1.40 | |||||||||
H-2 | 21111111 | 498 | 51% | 1.06 | 0.90–1.25 | 226 | 50% | 1.01 | 0.80–1.28 | 272 | 52% | 1.10 | 0.86–1.42 | |||||||||
H-3 | 21222221 | 122 | 39%* | 0.63 | 0.44–0.89 | 61 | 44% | 0.80 | 0.50–1.27 | 62 | 33%† | 0.49 | 0.28–0.81 | |||||||||
H-4 | 22111111 | 93 | 47% | 0.89 | 0.61–1.30 | 50 | 45% | 0.81 | 0.47–1.36 | 43 | 50% | 0.99 | 0.56–1.71 | |||||||||
H-5 | 11222221 | 86 | 49% | 0.95 | 0.59–1.50 | 37 | 40% | 0.68 | 0.30–1.42 | 49 | 55% | 1.24 | 0.67–2.36 | |||||||||
H-6 | 22111221 | 77 | 44% | 0.78 | 0.45–1.24 | 27 | 52% | 1.07 | 0.49–2.40 | 50 | 40% | 0.65 | 0.31–1.23 | |||||||||
H-7 | 11111121 | 77 | 57% | 1.33 | 0.90–2.06 | 43 | 55% | 1.23 | 0.71–2.20 | 33 | 60% | 1.48 | 0.75–3.19 | |||||||||
H-8 | 21111121 | 27 | 43% | 0.75 | 0.38–1.60 | 12 | 38% | 0.61 | 0.18–2.22 | 15 | 47% | 0.88 | 0.38–3.09 | |||||||||
H-9 | 22111222 | 16 | 45% | 0.83 | 0.25–2.68 | 10 | 50% | 1.00 | —‡ | 6 | 38% | 0.60 | —‡ |
Each haplotype is listed by the ordered set of alleles present at each of the eight IL4R SNPs. A “1” denotes the presence of the reference allele, and a “2” indicates the presence of the variant allele for each SNP, as described in Table 1. Haplotypes are also given arbitrary haplotype ID numbers (H-1 through -9). The table shows only haplotypes that were observed in at least 15 effective transmission events in the affected children (representing 1,546/1,658 = 93% of observed haplotypes in informative transmission events). For each haplotype, percent transmitted is the ratio of actual transmissions to possible transmissions to affected children of informative parents. The ORs compare actual with Mendelian segregation ratios. The 95% CI column lists the upper and lower bounds of the 95% CI of the OR, calculated based on the number of families (not shown) informative for that haplotype. P values are based on permutation.
Significant after adjustment at P = 0.021;
significant after adjustment at P = 0.015;
CIs are not reported in instances where 10 or fewer families were informative for the haplotype.
Article Information
This work was funded in part by a Career Development Award to J.A.N. from the American Diabetes Association and by National Institutes of Health grants R01 AI29042-09A1 (to H.E.) and R01 GM62628-02 and 5R01 GM53275-08 (to L.C.L.).
We thank Amy White for technical assistance and our colleagues at Roche Molecular Systems for comments and discussion.
REFERENCES
Address correspondence and reprint requests to Henry A. Erlich, Roche Molecular Systems, 1145 Atlantic Ave., Alameda CA 94501. E-mail: [email protected].
Received for publication 8 February 2002 and accepted in revised form 26 July 2002.
HBDI, Human Biological Data Interchange; IBD, identity by descent; IL4, interleukin-4; LD, linkage disequilibrium; LOD, logarithm of odds; OR, odds ratio; SNP, single-nucleotide polymorphism; TDT, transmission disequilibrium test.