OBJECTIVE

Islet autoimmunity develops before clinical type 1 diabetes and includes multiple and single autoantibody phenotypes. The objective was to determine age-related risks of islet autoantibodies that reflect etiology and improve screening for presymptomatic type 1 diabetes.

RESEARCH DESIGN AND METHODS

The Environmental Determinants of Diabetes in the Young study prospectively monitored 8,556 genetically at-risk children at 3- to 6-month intervals from birth for the development of islet autoantibodies and type 1 diabetes. The age-related change in the risk of developing islet autoantibodies was determined using landmark and regression models.

RESULTS

The 5-year risk of developing multiple islet autoantibodies was 4.3% (95% CI 3.8–4.7) at 7.5 months of age and declined to 1.1% (95% CI 0.8–1.3) at a landmark age of 6.25 years (P < 0.0001). Risk decline was slight or absent in single insulin and GAD autoantibody phenotypes. The influence of sex, HLA, and other susceptibility genes on risk subsided with increasing age and was abrogated by age 6 years. Highest sensitivity and positive predictive value of multiple islet autoantibody phenotypes for type 1 diabetes was achieved by autoantibody screening at 2 years and again at 5–7 years of age.

CONCLUSIONS

The risk of developing islet autoimmunity declines exponentially with age, and the influence of major genetic factors on this risk is limited to the first few years of life.

Age is a modulator of autoimmune diseases (1). Most autoimmune diseases show periods of high and low incidence, and these differences are likely to have an etiological basis. Type 1 diabetes is an autoimmune disease in which the presymptomatic stage is defined by the presence of autoantibodies against multiple islet autoantigens (2). Most children with multiple islet autoantibodies progress to clinical type 1 diabetes (3). The peak incidence of the development of islet autoantibodies occurs in the first years of life (46). This early peak also indicates that the risk of developing islet autoantibodies attenuates with age and that age influences the child’s risk for developing islet autoantibodies, as demonstrated in relatives of patients with type 1 diabetes (7).

The rate at which the risk of islet autoantibody seroconversion attenuates may also provide insight into the timing of events that cause or protect against autoimmunity. Here, we examined the age-related decline in risk in The Environmental Determinants of Diabetes in the Young (TEDDY) study (8). This prospective study monitored from birth the development of islet autoantibodies in >8,000 children, who were selected based on their genetic predisposition to type 1 diabetes. The TEDDY study includes quarterly measurements of autoantibodies to multiple pancreatic antigens, providing an unprecedented ability to model seroconversion to distinct autoantibody phenotypes throughout childhood. The aim of this study was to determine the rate of change in the risk of type 1 diabetes-associated autoimmunity with increasing age and relative to genetics, demographics, and islet autoantibody phenotypes. The findings were used to identify optimal ages for islet autoantibody screening and to provide etiological models of islet autoimmunity.

Study Design

TEDDY is a prospective cohort study that was designed to identify environmental causes of type 1 diabetes (ClinicalTrials.govNCT00279318). It includes three clinical research centers in the U.S. (Colorado, Georgia/Florida, Washington State) and three in Europe (Finland, Germany, Sweden). The study design and methods were previously published (8). To identify infants at high genetic risk, newborns were screened for HLA DR-DQ genotypes associated with type 1 diabetes (Supplementary Methods). Participants were enrolled in the prospective follow-up before 4.5 months of age. They were followed up every 3 months until age 4 years and at least semiannually thereafter until 15 years of age. For all participants, written informed consent was obtained from a parent or primary caretaker for genetic screening, and separate consent was obtained for participation in the prospective follow-up. The study was approved by the local Institutional Review or Ethics Boards at each site and is monitored by an External Evaluation Committee formed by the National Institutes of Health.

Participants

Between September 2004 and February 2010, the TEDDY study screened 424,788 newborns. Of those, 21,589 had HLA genotypes of interest associated with increased risk. The parents of 8,676 infants consented to the follow-up study. Of these, 120 children, whose eligibility was not confirmed by high-resolution HLA typing, were excluded from the analyses. Therefore, 8,556 children (4,226 girls [49.4%]) were included in the analyses, of which 955 (11.2%) had a first-degree family history of type 1 diabetes (Supplementary Fig. 1 and Supplementary Table 1).

Assessments of Islet Autoantibodies (Study End Points)

Blood samples were obtained at each study visit to determine islet autoantibodies against insulin (IAA), GAD (GADA), and insulinoma antigen-2 (IA-2A). Islet autoantibody positivity (persistent confirmed) was defined as specific autoantibody presence on two consecutive or more visits 3 months apart and confirmed at two laboratories. The date of positivity was the date on which the first sample of the two consecutive samples confirming islet autoantibody positivity was taken. Multiple islet autoantibody positivity (stage 1 type 1 diabetes) was defined as multiple persistent confirmed islet autoantibodies. The date of positivity for multiple islet autoantibodies was the date on which the first sample of the two consecutive samples of the second appearing autoantibody was taken.

Islet Autoantibody Measurements and Single Nucleotide Polymorphism Typing

The levels of IAA, GADA, and IA-2A were measured in two laboratories by radiobinding assays, as previously described (8). In the U.S., all sera were assayed at the Barbara Davis Center for Diabetes at the University of Colorado Denver. In Europe, all sera were assayed at the University of Bristol, Bristol, U.K. Both laboratories have previously shown high sensitivity and specificity as well as excellent concordance (9). Single nucleotide polymorphism genotyping was performed at the Center for Public Health Genomics, University of Virginia, using the Illumina Immunochip (10). Genetic risk scores were calculated using 41 single nucleotide polymorphisms, as previously described (11).

Statistical Analyses

We performed analyses of data collected until 30 June 2019. The cumulative risks of islet autoantibodies were examined using the Kaplan-Meier method, and between-group comparisons were done using the log-rank test. To evaluate the risk at different ages, we calculated the cumulative islet autoantibody risks for children who remained islet autoantibody negative from the respective landmark to the outcome status. The autoantibody outcomes were as follows: 1) any islet autoantibody if a child was positive at two laboratories in two consecutive samples and 2) multiple islet autoantibodies (stage 1 type 1 diabetes) if a child was positive for two or more islet autoantibodies in two consecutive samples. The outcomes were subclassified into multiple islet autoantibody phenotypes (IAA first to multiple and GADA first to multiple) based on the islet autoantibody present in the first positive sample and as multiple first if the first autoantibody was unknown because the first positive sample had more than one islet autoantibody. IAA and GADA positivity that did not progress to multiple islet autoantibodies were subclassified as single IAA or single GADA if the child had been monitored and tested for ≥3 years after antibody development or if the child developed diabetes without developing multiple islet autoantibodies.

Landmark models (12) were based on the cumulative incidence curves from the respective landmark ages at 7.5 months, 2.125 years, 4.25 years, 6.25 years, and 8.25 years, which were chosen as the time points between scheduled visits. The follow-up time or horizon was set to 5 years from the landmark age. Because we expected a decline in risk (7), one-phase exponential decay functions were used to describe the 5-year horizon risk with increasing age. The curves describe the residual risk of developing islet autoantibodies in the subsequent 5 years for children who were negative at each landmark age. Curves were generated from the 5-year horizon risk at the landmark ages of 7.5 months and at 6-month intervals until 3.625 years, and at 4.25 years and at 12-month intervals until 8.25 years. For single antibody outcomes, the landmark ages >4.25 years were excluded because an additional 3 years of follow-up were required to define these outcomes. Parameters used to describe the risks were the baseline risk at the age of 7.5 months; the plateau risk, which refers to the lowest 5-year horizon risk reached by decay before leveling off; the half-life, which is the time required for the risk between baseline and plateau to halve; the decay constant (λ), which is a measure of the decay rapidity (the larger λ, the faster the decline to the plateau is); and the r value, which represents the goodness of fit of the exponential decay curve. Relative risks between groups across time were calculated as a ratio of the 5-year risks using these curve parameters. GraphPad Prism version 8.0.1 (GraphPad Software, San Diego, CA) was used to generate the exponential decay curves, to calculate the curve parameters, and to compare the values between groups by the extra sum of squares F test. To compare relative risks at different ages, the actual 5-year risks and their CIs of the risks were used to calculate a ratio of relative risks (13). All other analyses were performed using R version 3.6.1 using the package dynpred (version 0.1.2). For all comparisons (log-rank tests, paired Mann-Whitney U tests, and comparison of relative risks over time), a two-tailed P value of <0.05 was considered significant.

In total, 809 of 8,556 children (9.46%, 95% CI 8.85–10.1) enrolled in the TEDDY study developed islet autoantibodies at a median age of 3.2 years (interquartile range [IQR] 1.5–6.5), including 471 who developed multiple islet autoantibodies within a median time of 0.3 years (IQR 0.0–1.0) after becoming islet autoantibody positive (Supplementary Fig. 1). The first autoantibody-positive samples in children who developed multiple islet autoantibodies most frequently contained IAA only (166), GADA only (161), or multiple islet autoantibodies (133). Most of the 338 children with one islet autoantibody had IAA (136) or GADA (190). Most (78.4%) single autoantibody-positive children were monitored for ≥3 years or developed diabetes without testing positive for multiple islet autoantibodies.

The Risk of Islet Autoantibody Seroconversion Declined With Age

The risk of developing islet autoantibodies over a 5-year horizon decreased with increasing landmark age (Fig. 1 and Supplementary Table 2). At the landmark age of 7.5 months, the risk over a 5-year horizon was 6.3% (95% CI 5.7–6.8) for any islet autoantibodies and 4.3% (95% CI 3.8–4.7) for multiple islet autoantibodies compared with 3.2% (95% CI 2.7–3.6) and 1.1% (95% CI 0.8–1.3), respectively, at a landmark age of 6.25 years (both P < 0.0001). This translates to a fourfold reduction in the risk of developing multiple islet autoantibodies (stage 1 type 1 diabetes) among children who are islet autoantibody negative at age 6 years.

Figure 1

Cumulative risk of developing islet autoantibodies at each landmark age. Cumulative risks (line with shaded 95% CI) of developing any (A) or multiple islet autoantibodies (B) in children who were negative for the respective outcome at the ages of 7.5 months (red), 2.125 years (green), 4.25 years (blue), 6.25 years (dark green), and 8.25 years (purple). Exponential decay curves for the 5-year horizon risks for any islet autoantibodies (C) and multiple islet autoantibodies (D). The curve is represented by the equation: Risk = risk at baseline (0.625 years) × exp(λ × age−0.625) − plateau risk. Exponential decay curves are shown for the 5-year horizon risks of developing multiple islet autoantibodies in children with a first-degree relative with type 1 diabetes (black, n = 955) and children in the general population (green, n = 7,601) (E); children with the HLA DR3/4-DQ8 (red, n = 3,339), DR4-DQ8/DR4-DQ8 (black, n = 1,674), DR4-DQ8/DR8 (blue, n = 1,474), and DR3/3 (green, n = 1,791) genotypes (F); and children in the general population with the HLA DR3/DR4-DQ8 or DR4-DQ8/DR4-DQ8 genotypes stratified by their genetic risk score into the upper score quartile (red, n = 1,104), 25th to 75th quartiles (black, n = 2,206), and lower quartile (blue, n = 1,103) (G). H: Also shown is the 5-year risk for multiple islet autoantibodies in children with genetic risk scores in the highest quartile relative to the children in the lowest quartile (blue) and to the children in the 25th to 75th quartiles (black).

Figure 1

Cumulative risk of developing islet autoantibodies at each landmark age. Cumulative risks (line with shaded 95% CI) of developing any (A) or multiple islet autoantibodies (B) in children who were negative for the respective outcome at the ages of 7.5 months (red), 2.125 years (green), 4.25 years (blue), 6.25 years (dark green), and 8.25 years (purple). Exponential decay curves for the 5-year horizon risks for any islet autoantibodies (C) and multiple islet autoantibodies (D). The curve is represented by the equation: Risk = risk at baseline (0.625 years) × exp(λ × age−0.625) − plateau risk. Exponential decay curves are shown for the 5-year horizon risks of developing multiple islet autoantibodies in children with a first-degree relative with type 1 diabetes (black, n = 955) and children in the general population (green, n = 7,601) (E); children with the HLA DR3/4-DQ8 (red, n = 3,339), DR4-DQ8/DR4-DQ8 (black, n = 1,674), DR4-DQ8/DR8 (blue, n = 1,474), and DR3/3 (green, n = 1,791) genotypes (F); and children in the general population with the HLA DR3/DR4-DQ8 or DR4-DQ8/DR4-DQ8 genotypes stratified by their genetic risk score into the upper score quartile (red, n = 1,104), 25th to 75th quartiles (black, n = 2,206), and lower quartile (blue, n = 1,103) (G). H: Also shown is the 5-year risk for multiple islet autoantibodies in children with genetic risk scores in the highest quartile relative to the children in the lowest quartile (blue) and to the children in the 25th to 75th quartiles (black).

Close modal

The age-related attenuation of risk could be described by a one-phase exponential decay curve (r = 0.992 for any islet autoantibodies and r = 0.998 for multiple islet autoantibodies) (Fig. 1C and D). The half-life, plateau risk, and λ were 0.85 years (95% CI 0.7–1.0), 3.1% (95% CI 3.0–3.3), and 0.82 (95% CI 0.67–0.99), respectively, for any islet autoantibodies and 1.25 years (95% CI 1.1–1.4), 1.0% (95% CI 0.8–1.1), and 0.56 (95% CI 0.50–0.62), respectively, for multiple islet autoantibodies. The risk of developing any islet autoantibodies over a 5-year horizon in children who were negative at their last test was described by risk = 3.3 × exp(−0.82 × [age X − 0.625]) where age X is the age at the preceding negative test. The equation for multiple islet autoantibodies was risk = 3.3 × exp(−0.56 × [age X − 0.625])

The Genetic Influence on Multiple Islet Autoantibody Seroconversion Diminished With Age

The baseline cumulative risk of developing multiple islet autoantibodies varied by sex, first-degree family history status, site, HLA genotype, and genetic risk score (Supplementary Table 1). A decay in the 5-year risk with increasing age was observed in all subgroups examined (Fig. 1, Table 1, and Supplementary Tables 36). The risk decay was similar between children with a first-degree relative with type 1 diabetes (λ, 0.59; 95% CI 0.52–0.68) and children in the general population (λ, 0.53; 95% CI 0.46–0.60; P = 0.44) (Fig. 1E) and between children in European (λ, 0.58; 95% CI 0.50–0.67) and U.S. sites (λ, 0.52; 95% CI 0.40–0.65; P = 0.33) (Supplementary Fig. 2A). In contrast, the characteristics of the risk decay varied by sex (Supplementary Fig. 2B) and, in particular, genetics (Fig. 1F–H). There were marked variations in the baseline risk (P < 0.0001), plateau risk (P < 0.0001), and risk decay (P < 0.0001) between the HLA genotypes (Fig. 1F). Risks differed by HLA genotype at a young age but converged with increasing age so that the remaining risk in autoantibody-negative children at age 6 years was similar between the three DR4-DQ8–containing genotypes. The risk also converged with increasing age when children were categorized by their genetic risk score (Fig. 1G) or INS genotype (Supplementary Fig. 2C). Convergence was confirmed for genetic risk score categories (Fig. 1H and Supplementary Table 7), INS genotype (Supplementary Fig. 2D), and DR4-DQ8–containing genotypes, but not the DR3-containing genotypes (Supplementary Fig. 2E) when the 5-year risk for multiple islet autoantibodies was expressed as relative risks.

Table 1

Characteristics of the one-phase exponential decay curve describing the attenuating risk of developing multiple islet autoantibodies over a 5-year horizon

Outcome and subgroupBaseline risk at 7.5 months landmark, % (95% CI)Half-life, years (95% CI)*Plateau risk, % (95% CI)Decay constant, λ (95% CI)
Multiple islet autoantibody 
 First-degree family history     
  Yes 9.9 (9.5–10.3) 1.2 (1.0–1.4) 1.6 (1.0–1.9) 0.59 (0.52–0.67) 
  No 3.6 (3.5–3.7) 1.3 (1.1–1.5) 0.9 (0.7–1.0) 0.53 (0.46–0.60) 
 Sex     
  Boys 4.6 (4.4–4.9) 1.6 (1.2–2.0) 0.9 (0.5–1.2) 0.45 (0.33–0.57) 
  Girls 3.9 (3.7–4.2) 1.0 (0.8–1.2) 1.0 (0.8–1.2) 0.69 (0.56–0.84) 
 Site     
  European sites 4.9 (4.7–5.1) 1.2 (1.0–1.4) 1.1 (0.9–1.2) 0.58 (0.50–0.67) 
  U.S. sites 3.6 (3.3–3.8) 1.3 (1.1–1.7) 0.8 (0.6–1.0) 0.52 (0.40–0.65) 
HLA genotype     
  HLA DR3/4-DQ8 6.2 (6.0–6.4) 1.6 (1.4–1.9) 0.8 (0.5–1.0) 0.43 (0.37–0.50) 
  HLA DR4-DQ8/DR4-DQ8 4.2 (3.9–4.5) 1.4 (1.0–1.9) 1.7 (1.4–1.9) 0.50 (0.36–0.67) 
  HLA DR4-DQ8/DR8 3.0 (2.7–3.4) 0.4 (0.3–0.6) 1.0 (0.8–1.1) 1.59 (1.09–2.37) 
  HLA DR3/3 1.7 (1.5–1.9) 0.8 (0.6–1.1) 0.3 (0.2–0.4) 0.89 (0.64–1.22) 
 Genetic risk score     
  Upper quartile 11.2 (10.6–11.9) 2.2 (1.9–2.4) 0 (0–0.4) 0.32 (0.29–0.37) 
  25th–75th percentile 5.4 (5.2–5.7) 1.1 (0.9–1.4) 1.7 (1.5–1.9) 0.62 (0.51–0.75) 
  Lower quartile 2.7 (2.1–3.2) 0.5 (0.2–1.3) 1.4 (1.1–1.6) 1.27 (1.12–1.61) 
Autoantibody phenotype     
 IAA first to multiple 1.7 (1.6–1.8) 0.85 (0.8–0.9) 0.1 (0.07–0.14) 0.81 (0.74–0.89) 
 Multiple first 1.2 (1.1–1.3) 0.9 (0.7–1.3) 0.4 (0.4–0.5) 0.74 (0.54–0.98) 
 GADA first to multiple 1.5 (1.4–1.6) 2.5 (1.7–4.0) 0.2 (0.0–0.4) 0.28 (0.17–0.41) 
 Single IAA 0.8 (0.7–0.9) 3.4 (1.1–4.4) 0.1 (0.0–0.4) 0.21 (0.16–0.64) 
 Single GADA 1.0 (0.8–1.2) No decay observed 
Outcome and subgroupBaseline risk at 7.5 months landmark, % (95% CI)Half-life, years (95% CI)*Plateau risk, % (95% CI)Decay constant, λ (95% CI)
Multiple islet autoantibody 
 First-degree family history     
  Yes 9.9 (9.5–10.3) 1.2 (1.0–1.4) 1.6 (1.0–1.9) 0.59 (0.52–0.67) 
  No 3.6 (3.5–3.7) 1.3 (1.1–1.5) 0.9 (0.7–1.0) 0.53 (0.46–0.60) 
 Sex     
  Boys 4.6 (4.4–4.9) 1.6 (1.2–2.0) 0.9 (0.5–1.2) 0.45 (0.33–0.57) 
  Girls 3.9 (3.7–4.2) 1.0 (0.8–1.2) 1.0 (0.8–1.2) 0.69 (0.56–0.84) 
 Site     
  European sites 4.9 (4.7–5.1) 1.2 (1.0–1.4) 1.1 (0.9–1.2) 0.58 (0.50–0.67) 
  U.S. sites 3.6 (3.3–3.8) 1.3 (1.1–1.7) 0.8 (0.6–1.0) 0.52 (0.40–0.65) 
HLA genotype     
  HLA DR3/4-DQ8 6.2 (6.0–6.4) 1.6 (1.4–1.9) 0.8 (0.5–1.0) 0.43 (0.37–0.50) 
  HLA DR4-DQ8/DR4-DQ8 4.2 (3.9–4.5) 1.4 (1.0–1.9) 1.7 (1.4–1.9) 0.50 (0.36–0.67) 
  HLA DR4-DQ8/DR8 3.0 (2.7–3.4) 0.4 (0.3–0.6) 1.0 (0.8–1.1) 1.59 (1.09–2.37) 
  HLA DR3/3 1.7 (1.5–1.9) 0.8 (0.6–1.1) 0.3 (0.2–0.4) 0.89 (0.64–1.22) 
 Genetic risk score     
  Upper quartile 11.2 (10.6–11.9) 2.2 (1.9–2.4) 0 (0–0.4) 0.32 (0.29–0.37) 
  25th–75th percentile 5.4 (5.2–5.7) 1.1 (0.9–1.4) 1.7 (1.5–1.9) 0.62 (0.51–0.75) 
  Lower quartile 2.7 (2.1–3.2) 0.5 (0.2–1.3) 1.4 (1.1–1.6) 1.27 (1.12–1.61) 
Autoantibody phenotype     
 IAA first to multiple 1.7 (1.6–1.8) 0.85 (0.8–0.9) 0.1 (0.07–0.14) 0.81 (0.74–0.89) 
 Multiple first 1.2 (1.1–1.3) 0.9 (0.7–1.3) 0.4 (0.4–0.5) 0.74 (0.54–0.98) 
 GADA first to multiple 1.5 (1.4–1.6) 2.5 (1.7–4.0) 0.2 (0.0–0.4) 0.28 (0.17–0.41) 
 Single IAA 0.8 (0.7–0.9) 3.4 (1.1–4.4) 0.1 (0.0–0.4) 0.21 (0.16–0.64) 
 Single GADA 1.0 (0.8–1.2) No decay observed 
*

The half-life describes the time required for the risk between baseline and plateau to halve.

The plateau risk represents the lowest risk reached by decay.

The decay constant is a measure of the decay rapidity, whereby a larger value represents a faster decline to the plateau risk.

The Risk Decay Was a Feature of Multiple but Not Single Islet Autoantibody Phenotypes

Most (97.2%) of the 809 children who developed islet autoantibodies could be placed in one of five antibody profiles: single IAA, single GADA, IAA first with progression to multiple islet autoantibodies, multiple islet autoantibodies in the first positive sample, and GADA first progressing to multiple islet autoantibodies (Supplementary Fig. 1). The risk decay differed markedly among these profiles (P < 0.0001) (Table 1, Fig. 2, and Supplementary Table 8). The risk decay was observed for each of the multiple islet autoantibody phenotypes (Fig. 2A). The IAA-first-to-multiple phenotype showed a steep risk decline (λ, 0.81; 95% CI 0.74–0.89) with a low plateau risk (0.1%). The multiple-antibody-first phenotype also showed a steep risk decline (λ, 0.74; 95% CI 0.54–0.98), but a higher plateau risk (0.4%; P < 0.0001) than the IAA-first-to-multiple phenotype. The risk decline was slower in the GADA-first-to-multiple phenotype (λ, 0.28; 95% CI 0.17–0.41; P < 0.0001). The previously reported association of IAA with HLA DR4-DQ8–containing genotypes (14) was observed for the IAA-first-to-multiple and the multiple-first phenotypes across all landmark ages (Supplementary Fig. 3A and B and Supplementary Table 9). The GADA-first-to-multiple phenotypes was associated with HLA DR3-containing genotypes at the youngest landmark ages (Supplementary Fig. 3C).

Figure 2

Exponential risk decay curves for islet autoantibody phenotypes. Islet autoantibodies were categorized into the multiple (A and B) and single (C) islet autoantibody phenotypes. A: Exponential 5-year horizon risk decay curves for developing the multiple-first (blue), IAA-first-to-multiple (red), and GADA-first-to-multiple (black) phenotypes. B: The 5-year risk for developing the IAA-first-to-multiple or multiple-first phenotype (black line) or in the GADA-first-to-multiple phenotype (red line) in children categorized by their genetic risk scores for genes other than the HLA DR-DQ genotype. The curves show the relative risks for children with scores in the highest quartile (upper 25th percentile) relative to the children with scores in the lower three quartiles (lower 75th percentile). C: Exponential 5-year horizon risk decay curves for developing single IAA (red) and single GADA (black) phenotypes.

Figure 2

Exponential risk decay curves for islet autoantibody phenotypes. Islet autoantibodies were categorized into the multiple (A and B) and single (C) islet autoantibody phenotypes. A: Exponential 5-year horizon risk decay curves for developing the multiple-first (blue), IAA-first-to-multiple (red), and GADA-first-to-multiple (black) phenotypes. B: The 5-year risk for developing the IAA-first-to-multiple or multiple-first phenotype (black line) or in the GADA-first-to-multiple phenotype (red line) in children categorized by their genetic risk scores for genes other than the HLA DR-DQ genotype. The curves show the relative risks for children with scores in the highest quartile (upper 25th percentile) relative to the children with scores in the lower three quartiles (lower 75th percentile). C: Exponential 5-year horizon risk decay curves for developing single IAA (red) and single GADA (black) phenotypes.

Close modal

The similarities in the DR4-DQ8–associated IAA-first-to-multiple and multiple-first phenotypes suggested that these may form one phenotype (IAA/multiple first) that is distinct to the DR3-associated GADA-first phenotype. We, therefore, examined these two outcomes separately and asked how risk declined in relation to susceptibility conferred by genes other than HLA DR-DQ (Fig. 2B and Supplementary Table 7). For both phenotypes, the relative risk in the children with genetic risk scores in the upper quartile compared with the remainder was ∼2 at age 7.5 months, and for the GADA-first phenotype, the relative risk peaked at ∼3 years. Both subsequently declined to parity (1.0).

Unlike the multiple islet autoantibody phenotypes, we found no risk decay with age for the single GADA phenotype and a slow exponential decay for the single IAA phenotype (λ, 0.21; 95% CI 0.16–0.64) (Fig. 2C). The previously reported association of GADA with DR3-DQ2 (15) was observed for single GADA, which was associated with the HLA DR3/3 genotype and also with HLA DR3/4-DQ8 after age 2 years (Supplementary Fig. 3E).

At the landmark age of 7.5 months, the 5-year horizon risk of developing the multiple islet autoantibody phenotype was twofold greater than the risk of developing the single autoantibody phenotype (4.3% vs. 2.0%; P < 0.0001). In contrast, at the landmark age of 6.25 years, the risk of the multiple autoantibody phenotype was twofold less than the risk of the single autoantibody phenotype (1.1% vs. 2.1%; P < 0.0001).

Optimal Islet Autoantibody Screening Ages

A practical aspect of the age-related islet autoantibody risk decay is how it translates into screening for presymptomatic type 1 diabetes. Important features are the ability to identify future cases of type 1 diabetes (i.e., sensitivity) and the positive predictive value (PPV) for developing type 1 diabetes in those identified as positive. Therefore, we examined the sensitivity by using type 1 diabetes by age 12 years (n = 331 children) as the outcome and the 5-year risk of developing type 1 diabetes as the measure of PPV in children who were found to be positive for multiple islet autoantibodies if screening was performed at single time points from age 1 year to 8 years (Fig. 3A). The sensitivity reached a peak if screening was performed at age 4 years (37.2%; 95% CI 32.0–42.6). At least one-third of the cases of type 1 diabetes (sensitivity >33%) could be identified by a single screening visit between 3 and 5 years of age. The PPV was highest when children were identified as positive for multiple antibodies at the age of 1 year (84%; 95% CI 60.7–93.5) or 2 years (68.1%; 95% CI 57.6–76.0). The PPV remained ∼50% when identified from the age of 3 years.

Figure 3

Sensitivity and PPV of multiple islet autoantibodies. Screening for multiple islet autoantibodies was simulated at a single time point (A) and at two time points (B). A: For a single time point, the sensitivity of multiple islet autoantibodies to identify all cases of type 1 diabetes that occurred in TEDDY children by the age of 12 years (blue symbols) and the PPV calculated as the 5-year risk of progressing to type 1 diabetes in children positive for multiple islet autoantibodies (black symbols) are shown for screening at the ages of 1–8 years. Error bars represent the 95% CI. B: For the strategy of screening at two time points, the sensitivity of multiple islet autoantibodies to identify all cases of type 1 diabetes that occurred in TEDDY children by the age of 12 years is given for each combination.

Figure 3

Sensitivity and PPV of multiple islet autoantibodies. Screening for multiple islet autoantibodies was simulated at a single time point (A) and at two time points (B). A: For a single time point, the sensitivity of multiple islet autoantibodies to identify all cases of type 1 diabetes that occurred in TEDDY children by the age of 12 years (blue symbols) and the PPV calculated as the 5-year risk of progressing to type 1 diabetes in children positive for multiple islet autoantibodies (black symbols) are shown for screening at the ages of 1–8 years. Error bars represent the 95% CI. B: For the strategy of screening at two time points, the sensitivity of multiple islet autoantibodies to identify all cases of type 1 diabetes that occurred in TEDDY children by the age of 12 years is given for each combination.

Close modal

Screening for any islet autoantibodies increased sensitivity, especially at younger ages, but was associated with lower PPV than that of multiple islet autoantibodies at all ages (P = 0.008) (Supplementary Table 10). Except for children identified at the age of 1 year, the PPV in children who were identified with a single islet autoantibody was <15%. A two-point screening model for multiple islet autoantibodies was also considered (Fig. 3B). The sensitivity increased to almost 50% if screening for multiple islet autoantibodies was performed at the age of 2 years and again at 5 or 6 years, or at the age of 3 years and again at 6 or 7 years. Of the 331 who developed type 1 diabetes by age 12 years, 43 (13.0%) developed diabetes before age 2 years and 16 (4.8%) developed multiple islet autoantibodies after age 7 years.

We used landmark modeling to assess the influence of age on the risk of developing type 1 diabetes-associated autoantibodies in genetically at-risk children. The risk of developing multiple islet autoantibodies was decreased markedly with age. The influence of HLA and other susceptibility genes on risk also declined and was absent or very weak at the landmark age of 6.25 years. In contrast, the risk of developing one, but not multiple islet autoantibodies, remained constant and was minimally influenced by HLA genotypes.

The strengths of this study include the large number of children monitored from birth for up to 14 years, the high frequency and rigor of islet autoantibody testing, and the relatively large number of outcomes for each of the autoantibody phenotypes analyzed. The attenuated risk observed in our study is consistent with a recent report of children from affected families (7) and in the Diabetes Prevention Trial–Type 1 (DPT-1) study (16). Although the risk did not decrease to zero, both studies showed that the residual risk of developing multiple islet autoantibodies decreased fourfold by 6 years of age, an important finding for families with a family member with type 1 diabetes or a child with genetic features of type 1 diabetes.

A striking feature of the risk decay model was the diminishing influence of factors such as sex, HLA genotype, and other type 1 diabetes susceptibility genes with age. This was particularly evident for the HLA DR3/4-DQ8 genotype, which has the strongest association with the development of islet autoantibodies and type 1 diabetes (14,17). The dominant effect of this genotype over other HLA DR4-DQ8–containing genotypes was absent at the landmark age of 6.25 years for all multiple islet autoantibody phenotypes. We also observed a decreasing influence of the INS gene and the genetic risk score on multiple islet autoantibody risk. These data suggest that a portion of the genetic contribution to the development of multiple islet autoantibodies and, therefore, type 1 diabetes subsides before puberty. The relative reduction in the frequency of HLA DR-DQ risk genotypes in patients with older-onset type 1 diabetes compared with childhood onset (18), the increased contribution of several susceptibility genes to younger-onset type 1 diabetes (19), and the relatively small impact of type 1 diabetes susceptibility genes with progression to diabetes in children positive for multiple islet autoantibodies (20,21) are consistent with this finding.

It is important to note that our findings are relevant to childhood islet autoimmunity and type 1 diabetes onset. Patients who develop type 1 diabetes as adults have some, but not all, of the features of childhood type 1 diabetes. Adults often have fewer islet autoantibodies with an excess of the single GADA phenotype at diabetes onset (22). No study has monitored sufficient adult cases from early age to know when autoimmunity first appears. It is likely, however, that a portion of the adult cases derives from the children who have the single DR3-associated GADA phenotype. This phenotype appeared relatively frequently in late childhood or adolescence compared with other phenotypes.

The risk decay with age is likely to reflect pathogenetic mechanisms and the etiology of type 1 diabetes-associated islet autoimmunity. We suggest three potential models (Supplementary Fig. 4). The stimulus decay model proposes that environmental stimuli for developing multiple islet autoantibodies are strongest or most frequent in the 1st year of life and attenuate thereafter. The early stimulus model suggests that the stimuli are restricted to the 1st year of life and that the likelihood of responding with islet autoantibodies is an exponential function of time. In both models, the potential roles of exposures or the host’s responses to the exposures are strongest in the 1st year of life. The third model, called the response decay model, suggests that the stimuli are present throughout childhood but that the propensity for autoimmune responses decreases exponentially with age. Observations that support a response decay model include the association between type 1 diabetes and early but not later viral infection (23) and the faster rate of progression to diabetes in children who develop multiple islet autoantibodies in the first years of life (3). Understanding which of the three models correctly reflects the pathogenesis has important implications for primary prevention strategies. Models 1 and 2 will likely require strategies to reduce exposure, whereas model 3 will require modulation of host response.

The findings also have implications for when to test children for islet autoantibodies to diagnose presymptomatic type 1 diabetes and recruit into clinical trials (24,25). No single time point captured the majority of cases of type 1 diabetes that occurred by the age of 12 years in the study. The highest sensitivity was achieved by screening between 3 and 5 years of age, similar to the screening age in a general population in the Fr1da study (24). Screening genetically at-risk TEDDY children at two time points increased sensitivity. A combination of testing for islet autoantibodies at the age of 2 years and then at 5–7 years may be advantageous by allowing us to identify at the age of 2 years young children positive for multiple islet autoantibodies with a relatively fast progression to clinical type 1 diabetes. Further testing beyond age 7 years was more likely to identify children who developed persistent single islet autoantibodies than children who developed multiple islet autoantibodies. Further refinement of the single islet autoantibody phenotypes using alternative methods may help us to distinguish single islet autoantibody phenotypes relevant to type 1 diabetes (26).

A limitation of this study is the a priori selection of children based on their HLA genotype, which limits the generalizability of the findings to lower genetic risk strata that represent ∼50% of type 1 diabetes cases throughout life (14). Another limitation is the smaller number of children monitored until teenage years, which prevents us from projecting the residual risk into late adolescence or adulthood and may lead to misclassification of some children positive for a single antibody (27). It is also possible that more frequent sampling would classify some of the multiple-first-islet-autoantibody outcomes as IAA first to multiple or GADA first to multiple. Of note, the multiple-first phenotype shared features with the IAA-first-to-multiple phenotype, and both phenotypes differed from the GADA-first-to-multiple phenotype in the decay constant and HLA associations. Therefore, we suggest that the multiple-first- and the IAA-first-to-multiple are similar phenotypes that may represent a distinct endotype to the GADA-first-to-multiple phenotype (17,28). The association of GADA with diseases that manifest throughout life also supports this hypothesis (19,29). We did not include autoantibodies to zinc transporter 8 in the analysis because these were only measured in children who had developed other islet autoantibodies.

In conclusion, we have shown that the first years of life represent a period of heightened risk of developing type 1 diabetes-relevant autoimmunity and that the risk declined exponentially with age in childhood. The risk decline occurred regardless of the genetic and demographic background but varied among the autoantibody response profiles. The findings have practical implications for implementing population-based screening and have led to models for explaining the development of islet autoimmunity. We suggest that our understanding of the etiological causes of type 1 diabetes will be improved by the findings of future studies aimed at confirming or refuting these models.

Clinical trial reg. no. NCT00279318, clinicaltrials.gov

See accompanying articles, pp. 2189 and 2269.

This article contains supplementary material online at https://doi.org/10.2337/figshare.13643756.

Funding. The TEDDY Study was funded by and National Institute of Diabetes and Digestive Kidney Diseases grants U01 DK63829, U01 DK63861, U01 DK63821, U01 DK63865, U01 DK63863, U01 DK63836, U01 DK63790, UC4 DK63829, UC4 DK63861, UC4 DK63821, UC4 DK63865, UC4 DK63863, UC4 DK63836, UC4 DK95300, UC4 DK100238, UC4 DK106955, UC4 DK112243, UC4 DK117483, and Contract No. HHSN267200700014C, by the National Institute of Allergy and Infectious Diseases, Eunice Kennedy Shriver National Institute of Child Health and Human Development, the National Institute of Environmental Health Sciences, the Centers for Disease Control and Prevention, and the JDRF, and in part by the National Institutes of Health/National Center for Advancing Translational Sciences Clinical and Translational Science Awards to the University of Florida (UL1 TR000064) and the University of Colorado (UL1 TR002535). This work was also supported by funds from the German Federal Ministry of Education and Research (FKZ: 01KX1818) and grants awarded to the German Center for Diabetes Research (DZD e.V.).

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. E.B. and A.-G.Z. designed the study, proposed the analyses, interpreted the findings, and wrote the manuscript. A.W. performed the analyses and drafted and revised the manuscript. C.W. acquired the data, provided input into the analytical plan, and drafted and revised the manuscript. M.H. provided input into the analytical plan and drafted and revised the manuscript. M.J.R., J.T., Å.L., J.-X.S., W.A.H., J.P.K., K.V., D.A.S., and B.A. designed the study and reviewed and edited the manuscript. E.B. and A.-G.Z. are the guarantors of this work, and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

1.
Winkler
C
,
Jolink
M
,
Knopff
A
, et al
.
Age, HLA, and sex define a marked risk of organ-specific autoimmunity in first-degree relatives of patients with type 1 diabetes
.
Diabetes Care
2019
;
42
:
1684
1691
2.
Insel
RA
,
Dunne
JL
,
Atkinson
MA
, et al
.
Staging presymptomatic type 1 diabetes: a scientific statement of JDRF, the Endocrine Society, and the American Diabetes Association
.
Diabetes Care
2015
;
38
:
1964
1974
3.
Ziegler
AG
,
Rewers
M
,
Simell
O
, et al
.
Seroconversion to multiple islet autoantibodies and risk of progression to diabetes in children
.
JAMA
2013
;
309
:
2473
2479
4.
Ziegler
AG
;
BABYDIAB-BABYDIET Study Group
.
Age-related islet autoantibody incidence in offspring of patients with type 1 diabetes
.
Diabetologia
2012
;
55
:
1937
1943
5.
Parikka
V
,
Näntö-Salonen
K
,
Saarinen
M
, et al
.
Early seroconversion and rapidly increasing autoantibody concentrations predict prepubertal manifestation of type 1 diabetes in children at genetic risk
.
Diabetologia
2012
;
55
:
1926
1936
6.
Krischer
JP
,
Lynch
KF
,
Schatz
DA
, et al.;
TEDDY Study Group
.
The 6 year incidence of diabetes-associated autoantibodies in genetically at-risk children: the TEDDY study
.
Diabetologia
2015
;
58
:
980
987
7.
Hoffmann
VS
,
Weiß
A
,
Winkler
C
, et al
.
Landmark models to define the age-adjusted risk of developing stage 1 type 1 diabetes across childhood and adolescence
.
BMC Med
2019
;
17
:
125
8.
TEDDY Study Group
.
The Environmental Determinants of Diabetes in the Young (TEDDY) study: study design
.
Pediatr Diabetes
2007
;
8
:
286
298
journal
9.
Bonifacio
E
,
Yu
L
,
Williams
AK
, et al
.
Harmonization of glutamic acid decarboxylase and islet antigen-2 autoantibody assays for national institute of diabetes and digestive and kidney diseases consortia
.
J Clin Endocrinol Metab
2010
;
95
:
3360
3367
10.
Törn
C
,
Hadley
D
,
Lee
HS
, et al.;
TEDDY Study Group
.
Role of type 1 diabetes-associated SNPs on risk of autoantibody positivity in the TEDDY study
.
Diabetes
2015
;
64
:
1818
1829
11.
Bonifacio
E
,
Beyerlein
A
,
Hippich
M
, et al.;
TEDDY Study Group
.
Genetic scores to stratify risk of developing multiple islet autoantibodies and type 1 diabetes: a prospective study in children
.
PLoS Med
2018
;
15
:
e1002548
12.
Putter
H
,
van Houwelingen
HC
.
Understanding landmarking and its relation with time-dependent Cox regression
.
Stat Biosci
2017
;
9
:
489
503
13.
Altman
DG
,
Bland
JM
.
Interaction revisited: the difference between two estimates
.
BMJ
2003
;
326
:
219
14.
Erlich
H
,
Valdes
AM
,
Noble
J
, et al.;
Type 1 Diabetes Genetics Consortium
.
HLA DR-DQ haplotypes and genotypes and type 1 diabetes risk: analysis of the type 1 diabetes genetics consortium families
.
Diabetes
2008
;
57
:
1084
1092
15.
Hagopian
WA
,
Sanjeevi
CB
,
Kockum
I
, et al
.
Glutamate decarboxylase-, insulin-, and islet cell-antibodies and HLA typing to detect diabetes in a general population-based study of Swedish children
.
J Clin Invest
1995
;
95
:
1505
1511
16.
Vehik
K
,
Haller
MJ
,
Beam
CA
, et al.;
DPT-1 Study Group
.
Islet autoantibody seroconversion in the DPT-1 study: justification for repeat screening throughout childhood
.
Diabetes Care
2011
;
34
:
358
362
17.
Krischer
JP
,
Liu
X
,
Vehik
K
, et al.;
TEDDY Study Group
.
Predicting islet cell autoimmunity and type 1 diabetes: an 8-year TEDDY study progress report
.
Diabetes Care
2019
;
42
:
1051
1060
18.
Caillat-Zucman
S
,
Garchon
HJ
,
Timsit
J
, et al
.
Age-dependent HLA genetic heterogeneity of type 1 insulin-dependent diabetes mellitus
.
J Clin Invest
1992
;
90
:
2242
2250
19.
Inshaw
JRJ
,
Cutler
AJ
,
Crouch
DJM
,
Wicker
LS
,
Todd
JA
.
Genetic variants predisposing most strongly to type 1 diabetes diagnosed under age 7 years lie near candidate genes that function in the immune system and in pancreatic β-cells
.
Diabetes Care
2020
;
43
:
169
177
20.
Beyerlein
A
,
Bonifacio
E
,
Vehik
K
, et al.;
TEDDY Study Group
.
Progression from islet autoimmunity to clinical type 1 diabetes is influenced by genetic factors: results from the prospective TEDDY study
.
J Med Genet
2019
;
56
:
602
605
21.
Redondo
MJ
,
Geyer
S
,
Steck
AK
, et al.;
Type 1 Diabetes TrialNet Study Group
.
A type 1 diabetes genetic risk score predicts progression of islet autoimmunity and development of type 1 diabetes in individuals at risk
.
Diabetes Care
2018
;
41
:
1887
1894
22.
Tuomi
T
,
Carlsson
A
,
Li
H
, et al
.
Clinical and genetic characteristics of type 2 diabetes with and without GAD antibodies
.
Diabetes
1999
;
48
:
150
157
23.
Beyerlein
A
,
Donnachie
E
,
Jergens
S
,
Ziegler
AG
.
Infections in early life and development of type 1 diabetes
.
JAMA
2016
;
315
:
1899
1901
24.
Ziegler
AG
,
Kick
K
,
Bonifacio
E
, et al.;
Fr1da Study Group
.
Yield of a public health screening of children for islet autoantibodies in Bavaria, Germany
.
JAMA
2020
;
323
:
339
351
25.
McQueen
RB
,
Geno Rasmussen
C
,
Waugh
K
, et al
.
Cost and cost-effectiveness of large-scale screening for type 1 diabetes in Colorado
.
Diabetes Care
2020
;
43
:
1496
1503
26.
Fouts
A
,
Pyle
L
,
Yu
L
, et al.;
Type 1 Diabetes TrialNet Study Group
.
Do electrochemiluminescence assays improve prediction of time to type 1 diabetes in autoantibody-positive TrialNet subjects?
Diabetes Care
2016
;
39
:
1738
1744
27.
Bosi
E
,
Boulware
DC
,
Becker
DJ
, et al.;
Type 1 Diabetes TrialNet Study Group
.
Impact of age and antibody type on progression from single to multiple autoantibodies in type 1 diabetes relatives
.
J Clin Endocrinol Metab
2017
;
102
:
2881
2886
28.
Battaglia
M
,
Ahmed
S
,
Anderson
MS
, et al
.
Introducing the endotype concept to address the challenge of disease heterogeneity in type 1 diabetes
.
Diabetes Care
2020
;
43
:
5
12
29.
Solimena
M
,
Folli
F
,
Denis-Donini
S
, et al
.
Autoantibodies to glutamic acid decarboxylase in a patient with stiff-man syndrome, epilepsy, and type I diabetes mellitus
.
N Engl J Med
1988
;
318
:
1012
1020
Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at https://www.diabetesjournals.org/content/license.