Trends in High-Risk HLA Susceptibility Genes Among Colorado Youth With Type 1 Diabetes

RESEARCH DESIGN AND METHODS—

Subjects included in this study were identified through two different diabetes registry studies conducted 13 years apart: the Colorado Type 1 Diabetes Registry and the Colorado Site of the SEARCH for Diabetes in Youth Study. The Colorado Type 1 Diabetes Registry was a statewide registry developed to ascertain all new cases of type 1 diabetes in 0- to 17-year-olds in Colorado from 1 January 1978 to 31 December 1988 (10,11). SEARCH is a multicenter, population-based observational registry that ascertained new cases of physician-diagnosed diabetes in youth 0–19 years of age from 1 January 2002 forward (12). A detailed description of this population (2), the Colorado Type 1 Diabetes Registry (10,11), and SEARCH have been published (12). Type 1 diabetes was defined as use of insulin within 2 weeks from diagnosis. Completeness of case ascertainment was assessed by the capture-recapture method (13,14) and estimated to be 96–97% over time (2).

Youth aged 5–17 years identified with incident type 1 diabetes during the period of 1978–1988 and 2002–2004 were eligible to participate in this analysis. In 1978–1988, all Hispanic youth (n = 120) and a similar size random sample of NHW (n = 122) youth identified by the Colorado Type 1 Diabetes Registry were eligible for genetic typing, and 41 Hispanic (34%) and 59 NHW youth (48%) had HLA measured (15). In 2002–2004, all youth over the age of 4 years of Hispanic (n = 87) or NHW (n = 512) origin were eligible for genetic typing. A total of 39 Hispanic (45%) and 225 NHW (44%) youth participated. Age at diagnosis, sex, and family history of diabetes were assessed to determine whether the youth with HLA measured were representative of the overall registered population in each period and racial/ethnic group.

Date of diagnosis and date of birth were obtained from medical records. Family history of diabetes was obtained from medical record abstraction (71% in 1978–1988 and 85% in 2002–2004) or from self-administered questionnaires completed within 5 years from the date of diagnosis (17% in 1978–1988 and 3% in 2002–2004). Children who had a sibling, parent, or grandparent with any type of diabetes were defined as having a positive family history.

To assess whether the sample of youth with available HLA typing was representative of the entire patient population of youth with type 1 diabetes diagnosed in these two time periods, demographic and clinical characteristics of study participants were compared with those of nonparticipants. Table 1 shows these characteristics according to time period and ethnicity. In 1978–1988, there were no significant differences for either NHW or Hispanic youth. In 2002–2004, the only significant difference between those with and without HLA typing was a higher frequency of a positive family history of diabetes in NHW participants with typing available. However, no significant relationship was noted between having the highest-risk HLA genotype (DRB1*03-DQB1*02/DRB1*04-DQB1*03) and having a positive family history of diabetes in either time period (P = 0.96). These data suggest that our samples of NHW and Hispanic youth with measured HLA are likely representative of the larger eligible population with type 1 diabetes in both periods, or that they differ in characteristics not strongly related with the frequency of high-risk HLA genotypes. Therefore, we believe that any observed trends or lack thereof in the frequency and distribution of HLA genotypes over time in our population are not likely to be due to selection bias. All participants provided informed consent and both the Colorado Type 1 Diabetes Registry and SEARCH were approved by relevant institutional review boards.

HLA typing

DRB1 and DQB1 typing was performed using OneLambda PCR assays that were standard for the time period in which they were used (1978–1988 and 2002–2004) (16). The main difference between assays was the number of DQB1 alleles identified that had sequence-specific oligonucleotides (SSOs) available for typing. There were seven DQB1 SSOs identified for typing in 1978–1988, and there were 100 DQB1 SSOs identified for typing in 2002–2004. Due to the difference in the number of DQB1 alleles identified and the availability of SSOs for typing these alleles over the two periods, only the gene locus was used in the analysis to compare genotypes (e.g., DQB1*02 vs. DQB1*0201). Additionally, the current PCR sequence–specific oligonucleotide probes used in 2002–2004 allowed for a larger number of samples to be typed at one time compared with the prior assay (17,18). Previous validation studies based on OneLambda Luminex SSO typing yielded similar results when compared with the earlier SSO assay (19).

Genotypes assessed in both time periods were categorized as high risk (DRB1*03-DQB1*02/DRB1*04-DQB1*03); moderate to low risk (DRB1*04-DQB1*03/DRB1*04-DQB1*03, DRB1*04-DQB1*03/X, where X is not DRB1*03-DQB1*02, and DRB1*04-DQB1*03/unknown); low risk (DRB1*03-DQB1*02/DRB1*03-DQB1*02, DRB1*03-DQB1*02/X, where X is not DRB1*04-DQB1*03, and DRB1*03-DQB1*02/unknown); and neutral risk (X/X or X/unknown). The X haplotype denotes all other haplotypes not defined above. The unknown haplotype denotes a haplotype that was unable to be typed.

Statistical analyses

Descriptive univariate analysis was used to determine the frequency of the genotypes for each period. χ² statistics were used to compare frequencies of HLA haplotypes by race and period and genotypes by high-risk (DRB1*03-DQB1*02/DRB1*04-DQB1*03) versus all other genotypes. Multivariate logistic regression was used to adjust for potential confounders (onset age, sex, family history of diabetes, and race/ethnicity) and assess temporal trends. High-risk was compared with all other genotypes using all other genotypes as the referent group. SAS version 9.1 (20) was used for all analyses.

RESULTS—

A total of 364 youth (100 in 1978–1988, and 264 in 2002–2004) with new-onset type 1 diabetes at age 5–17 years had HLA genes measured as described above.

Table 2 shows the frequency of HLA class II DRB1-DQB1 genotypes by time period and ethnicity. More children (both NHW and Hispanic) carried the high-risk genotype (DRB1*03-DQB1*02/DR04-DQB1*03) in 1978–1988 than in 2002–2004 (P = 0.05). Conversely, more children in 2002–2004 carried the moderate to low risk DRB1*04-DQB1*03 (homo- or heterozygous) genotypes than in 1978–1988 (P = 0.04). There was no significant difference by time period or ethnicity in frequency of the neutral risk, DRB1*03-DQB1*02 (homo- or heterozygous), or X/X genotypes. These patterns were similar for NHW and Hispanic youth with type 1 diabetes.

Using multiple logistic regression analysis (Table 2), we assessed the association between having the high-risk genotype (versus all other genotypes) and time period, controlling for age, sex, race/ethnicity, and family history of diabetes. There was a 40% lower odds for carrying the high-risk genotype in 2002–2004 versus 1978–1988 (odds ratio 0.60 [95% CI 0.36–0.99], P = 0.05). A decreasing trend over the two time periods was noted for both NHWs (13%) and Hispanics (6%). Conversely, there was a 70% greater chance of carrying the moderate- to low-risk genotypes in 2002–2004 versus 1978–1988 (1.7 [1.1–2.8], P = 0.03). No significant interactions between ethnicity and time period on the odds of having the highest risk genotype were noted. When stratified according to age-group (Table 3), the most dramatic decreasing trend in the prevalence of the high-risk genotype occurred in those aged 5–9 years (51% [1978–1988] vs. 32% [2002–2004], P = 0.05) with a coincident increasing trend in the proportion with the DRB1*04-DQB1*03 (homo- or heterozygous) genotype in the same age-group (19% [1978–1988] to 40% [2002–2004], P = 0.03). No significant differences over time were noted in the other age-groups.

CONCLUSIONS—

We found that the distribution of HLA genotype changed in youth with type 1 diabetes diagnosed 27 years apart. Having type 1 diabetes at age 5–17 years in 2002–2004 was associated with an 11% difference or 0.6-fold lower odds for carrying the high-risk HLA genotype (DRB1*03-DQB1*02/DRB1*04-DQB1*03) compared with having type 1 diabetes in 1978–1988. The largest proportional difference (19%) over time in the high-risk genotype was found in 5- to 9-year-olds, suggesting that environmental pressures may have a greater impact on type 1 diabetes risk at earlier ages. This is the first documentation of such a trend in a representative biracial sample of youth with type 1 diabetes in the U.S.

Our study supports previous findings from Finland (9) and the U.K (8), both including Caucasian adults with childhood-onset type 1 diabetes diagnosed up to 80 years apart. The U.K. study showed that the frequency of the high-risk genotype was 12% lower in adults diagnosed up to 15 years of age in 1985–2002 compared with those diagnosed between 1922–1946 and 21% lower in those diagnosed under 5 years of age (8). The Finnish study showed that the frequency of the high-risk genotype was 7.1% lower and that of protective genotypes was 7.2% higher in adults with childhood-onset type 1 diabetes diagnosed after 1990 compared with those diagnosed before 1965 (9). Two older studies conducted in Finland in the 1980s reported a decrease in HLA-DR3 by 25% for those diagnosed with type 1 diabetes in the 1980s compared with those diagnosed in the 1960s (6,7). However, the major limitation of these studies was the potential for selection bias. Neither of these studies provided information on the size or representativeness of the referent population. More importantly, two studies (8,9) measured HLA in adults who were diagnosed with type 1 diabetes as children in the first half of the 20th century, when early mortality associated with type 1 diabetes was high. Therefore, the frequency and distribution of HLA genotypes in this sample may be significantly influenced by factors associated with survival of type 1 diabetes. An underestimate of the proportional distribution of high-risk HLA genotypes in earlier time periods (due to a potential association of high-risk HLA genes with increased mortality) may have resulted in underestimating the true temporal trends in the proportion of type 1 diabetes cases with high-risk HLA genotypes in these earlier studies. In contrast, our study genotyped youth with type 1 diabetes relatively close to their disease onset and thus is not affected by survivor bias. In addition, we were able to provide evidence that the sample included in the analysis was representative of the Colorado population of youth with type 1 diabetes in both time periods.

Nevertheless, our study has several limitations. Importantly, although the same method was used for typing, due to the specific time periods in which they were performed, two slightly different assays were used, namely, OneLambda LabType Luminex SSO in 2002–2004 and OneLambda SSO in 1978–1988. We were unable to validate the results obtained with the earlier SSO assay due to the lack of stored DNA samples and inability to recontact participants. However, previous studies (19,21) suggest comparable agreement in serologic specificity with an 85% allele agreement. Due to this and the difference in the number of allele SSOs available for typing over time, only the gene locus was used in this analysis. By not using expanded allele information, we potentially introduced some misclassification, likely increasing the proportion with the high-risk genotype. However, misclassification would be similar for both time periods and therefore unlikely to influence our analysis of trends. Because for 2002–2004 we did have allele information on 100 DQB1 SSOs, we estimated the magnitude of the potential misclassification, and only four (1.5%) cases were misclassified as having the highest-risk genotype in 2002–2004 using the expanded allele information only. In addition, there were 12 youth from 1978–1988 that did not have complete genotypes. These youth were included in the study and classified as haplotype unknown. This categorization could have led to some misclassification. Specifically, if the unknown haplotype was DRB1*03-DQB1*02 or DRB1*04-DQB1*03, youth would have been classified as high risk based on the typed haplotype or as heterozygous for the above haplotypes. For example, in the extreme situation where all youth with the unknown haplotype in the first time period would have been in fact at “high-risk,” the proportion of youth with the highest risk genotype in 1978–1988 would have been 51% (instead of 39%). This would have made our trend estimates even more substantial (i.e., a 23% decrease in the proportion with high-risk genotype). Consequently, our estimated change of 11% is conservative.

There are several hypotheses proposed to explain the increase in type 1 diabetes incidence in youth and the earlier age at onset. The “accelerator” and “overload” hypotheses suggest that an increase in body mass in genetically susceptible individuals accelerates β-cell decline or stresses (overloads) the β-cells, leading to an early age at type 1 diabetes diagnosis (22,23). The “hygiene hypothesis” postulates that the decreasing early life exposure to infectious agents in Westernized societies has led to an impairment in the maturation of the immune system, thus permitting an increased occurrence of immune-mediated disorders, such as asthma and type 1 diabetes (24). Further testing is required for all of these hypotheses.

Given that there is an increasing incidence of type 1 diabetes in Colorado youth (2), our findings support the hypothesis that environmental pressures may be increasing the disease penetrance despite a shift in the distribution of HLA genotypes from high to moderate-low risk. Although this and previous studies (6–9) evaluating the distribution and frequency of HLA susceptibility genes over time have provided evidence in support of this hypothesis, the nature of environmental factors responsible for the increasing incidence of type 1 diabetes remains unknown.