Introduction

A short treatment with anti-CD3 antibodies at clinical onset of type 1 diabetes has been shown to preserve residual beta cell function for at least 18 months but only in patients with secretory function exceeding 25% of healthy controls as measured by hyperglycaemic clamp [1]. This opens perspectives for prevention trials in the late preclinical disease phase, where beta cell function is even better preserved than at diagnosis [2] and where anti-CD3 antibodies have also been shown to be effective in animal models [3]. However, such an undertaking requires effective identification of individuals at high risk of progression towards hyperglycaemia within a relatively short period, e.g. 5 years, to ensure that potential benefits outbalance the risk of side effects; this would also allow conclusions within a reasonable time frame [4].

Until now, risk assessment for type 1 diabetes has mostly been based on screening for insulin autoantibodies (IAA), autoantibodies against the 65 kDa isoform of GAD, i.e. GAD autoantibodies (GADA) and insulinoma-associated protein-2 (IA-2) autoantibodies (IA-2A) [59]. The presence of IA-2A has been shown to predict impending clinical onset in first-degree relatives of type 1 diabetic patients [6, 9], but allows detection of only about two out of three prediabetic individuals [10]. More recently, antibodies against the related antigen IA-2 beta (IA-2β, antibodies: IA-2βA) were reported to identify a subgroup of first-degree relatives with higher progression rate to diabetes [11, 12]. In 2007, zinc transporter-8 (ZnT8) was identified as a novel diabetes autoantigen [13]. Autoantibodies against the carboxy-terminal intracellular domain of ZnT8 (ZnT8A) have been proposed as additional markers of rapid disease progression [13], but their epitope specificity has been shown to vary according to a non-synonymous single nucleotide polymorphism at codon 325 of SLC30A8, the gene encoding ZnT8 [14].

In preparation of secondary prevention trials with immunointervention, we investigated to what extent screening for individuals at high risk of diabetes based on positivity for IA-2A [6, 9] or the presence of multiple autoantibodies (IAA, GADA and/or IA-2A) [58] could be optimised in terms of sensitivity and progression rate to diabetes by measuring IA-2βA and ZnT8A. A screening for both antibodies was conducted in first-degree relatives who tested positive for IAA, IA-2A or GADA. For the ZnT8A assay we used a radioligand derived from a hybrid ZnT8 cDNA construct generated by fusion of CR and CW (zinc transporter-8 carboxy-terminal constructs carrying 325Arg and 325Trp, respectively) (CW-CR), which is able to recognise antibodies directed against the wild type carboxyterminal domain of ZnT8 carrying 325Arg, as well as antibodies against a polymorphic variant carrying 325Trp [15]. The efficiency of antibody screening was compared before and after omission of relatives protected by HLA-DQ genotype or type of relationship to the diabetic proband [16].

Methods

Participants

Between 30 August 1989 and 17 December 2006, the Belgian Diabetes Registry (BDR) consecutively recruited 6,432 siblings or offspring (under age 40 at entry) of type 1 diabetic probands according to previously defined criteria [10, 16]. The probands are considered representative of the Belgian population of type 1 diabetic patients [17]. After obtaining written informed consent from each relative or their parents, a short questionnaire with demographic, familial and personal information was completed at each visit and blood samples were taken at entry and (as a rule) yearly thereafter. Only relatives (n = 5,635) with two or more contacts during follow-up, the last being at diagnosis in case of prediabetes, were included in this study. This allows the clinical status of relatives at this last time point to be unambiguously ascertained. The study was conducted in accordance with the guidelines in the Declaration of Helsinki as revised in 2008 (www.wma.net/en/30publications/10policies/b3/index.html, accessed 1 June 2009) and approved by the Ethics Committees of the BDR and the participating university hospitals. Blood was sampled at random, divided into aliquots and stored at −80°C until analysed for diabetes-associated autoantibodies and HLA-DQ genotype. Relatives were screened for the presence of IAA, GADA and IA-2A and 409 individuals were positive for at least one of these antibodies at baseline. The median age (interquartile range [IQR]) of these 409 siblings was 12 (6–19) years at baseline; there were 210 males and 199 females (ratio 1.06), 228 siblings, 92 offspring of a diabetic father and 89 offspring of a diabetic mother. Of the 409 siblings, 42% carried susceptible HLA-DQ genotypes and 25% protective ones [18]. During follow-up (median [IQR]: 68 [36–103] months), 86 (21%) antibody-positive relatives developed diabetes; they were identified by BDR as previously described [16].

Analytical methods

IA-2A [9], GADA [19] and IAA [16] were determined at baseline by liquid-phase radiobinding assays and HLA-DQ polymorphisms by allele-specific oligonucleotide genotyping [18] as described previously. HLA-DQ genotypes were classified as susceptible, protective or neutral on the basis of data from the BDR [18]. Measurements of IA-2βA and ZnT8A CW-CR [12, 14, 15] were also carried out at baseline by liquid-phase radiobinding assay using the same in-house protocol as for IA-2A and GADA, except that 35S-labelled GAD or 35S-labelled IA-2 (intracellular part) were replaced by 35S-labelled IA-2β (intracellular part) or 35S-labelled ZnT8A (intracellular carboxy-terminal part derived from CW-CR constructs). All tracers were purified by ultrafiltration (Amicon Ultra-4 filter units; Millipore, Billerica, MA, USA) and antibody concentrations expressed as percentage of added tracer bound (10,000 cpm/tube). cDNAs for the preparation of radioligands by in vitro transcription-translation were kind gifts of Å. Lernmark (University of Washington, Seattle, WA, USA) for full length 65 kDa GAD, M. Christie (King’s College School of Medicine and Dentistry, London, UK) for IA-2 (intracellular part), V. Lampasona (Instituto San Raffaele, Milano, Italy) for IA-2β (intracellular part; amino acids 662–1033) and J. C. Hutton (Barbara Davis Center for Childhood Diabetes, Aurora, CO, USA) for the dimeric CW-CR ZnT8 construct containing both CR encoding the wild type amino acids 268–369 carrying 325Arg, and CW, a variant carboxy-terminal construct carrying 325Trp. In the Diabetes Antibodies Standardization Program 2009 Workshop diagnostic sensitivity and specificity were respectively 74% and 97% for GADA, 40% and 98% for IAA, 66% and 99% for IA-2A, 53% and 98% for IA-2βA, and 68% and 100% for ZnT8A (CW-CR). Cut-off values for antibody positivity were determined as the 99th percentile of antibody levels in 761 non-diabetic controls and amounted to ≥0.6% for IAA, ≥2.6% for GADA, ≥0.44% for IA-2A, ≥0.39% for IA-2βA and ≥1.20% for ZnT8A. Between-day coefficients of variation determined for the Juvenile Diabetes Research Foundation standard serum were 9% for IAA (n = 413), 10% for GADA (n = 427), 11% for IA-2A (n = 474), 10% for IA-2βA (n = 156) and 8% for ZnT8A (n = 115).

Statistical analysis

Statistical differences between groups were assessed by means of the Mann–Whitney U test for continuous variables and by the χ 2 test, using Yates’ correction or Fisher’s exact test for categorical variables. McNemar’s test was used to assess differences between paired proportions. Kaplan–Meier analysis was used to estimate diabetes-free survival and survival curves were compared using the logrank test. Cox proportional hazards model, performed by forward stepwise method, was used to investigate the independent contributions of risk factors identified by univariate analysis, with calculation of 95% CIs on hazard ratios. In time-to-event analysis, follow-up started at entry and ended at the last contact with the relative or at clinical onset, whichever came first. Except for age and number of antibodies (between 1 and 5), all variables were introduced as categorical variables in Cox regression analysis. Univariate analyses (enter method) were first performed to identify which potential risk factors were significantly associated with progression to diabetes (Electronic supplementary material [ESM] Table 1). The significant variables selected in a first multivariate analysis (ESM Table 1) were then tested in a second model together with combinations of antibody positivity. All statistical tests were performed two-tailed by SPSS for Windows 16.0 (SPSS, Chicago, IL, USA), by EpiInfo version 6 (USD, Stone Mountain, GA, USA) or by GraphPad Prism version 4.00 for Windows (San Diego, CA, USA) and considered significant at p < 0.05.

Results

Progression to diabetes

Of the 409 first-degree relatives selected on the basis of positivity for IAA, GADA and/or IA-2A, 86 (21%) developed diabetes after a follow-up time (IQR) of 39 (15–65) months. Progression occurred more frequently among siblings (55/228; 24%) or offspring of a diabetic father (24/92; 26%) than among offspring of a diabetic mother (7/89; 8%) (p = 0.002). At entry, progressors were younger than non-progressors (median age [IQR]: 9 [5–16] years vs 12 [7–20] years; p = 0.01) and were more often carriers of susceptible HLA-DQ genotypes (72% vs 34%; p < 0.001).

Antibody prevalence at baseline

Among the 409 first-degree relatives, 50% were positive for IAA, 62% for GADA, 24% for IA-2A, 14% for IA-2βA and 20% for ZnT8A. The prevalence of IA-2βA and ZnT8A was significantly higher in IA-2A+ than in IA-2A first-degree relatives (Table 1) and in progressors than in non-progressors (IA-2βA 47% vs 6%, p < 0.001; ZnT8A 56% vs 10%, p < 0.001).

Table 1 Prevalence of IA-2βA and ZnT8A at baseline in 409 first-degree relatives according to antibody status for IAA, GADA, and IA-2A

Kaplan–Meier survival analysis

Stratification of relatives according to the status of IA-2A, IA-2βA or ZnT8A (ESM Fig. 1) showed that for each antibody, positivity was associated with a higher progression rate to diabetes than in their absence (p < 0.001). The presence of IA-2βA or ZnT8A significantly decreased diabetes-free survival in IA-2A (Fig. 1c, d), but not in IA-2A+ relatives (Fig. 1a, b). Combined screening for IA-2A and IA-2βA, or IA-2A and ZnT8A identified respectively two or eight prediabetic relatives more (Fig. 1e, f) during the entire follow-up than did screening for IA-2A alone (Fig. 1a, b). Combined testing for IA-2A IA-2βA and ZnT8A detected nine extra prediabetic relatives than testing for IA-2A alone with similarly high overall progression rate (not shown).

Fig. 1
figure 1

Diabetes-free survival (%) in IA-2A-positive first-degree relatives at baseline sampling (n = 97) stratified according to the presence or absence of IA-2βA (a) and ZnT8A (b); in IA-2A-negative first-degree relatives at baseline sampling (n = 312) stratified according to the presence or absence of IA-2βA (c) and ZnT8A (d); and in all antibody-positive first-degree relatives (n = 409) stratified according to the presence of IA-2A or IA-2βA vs their joint absence (e) and the presence of IA-2A or ZnT8A vs their joint absence (f). Dashed lines, positivity; continuous lines, negativity. Number of events (total number of relatives at entry) are indicated next to each arm for each panel. a p = 0.112 (logrank); (b) p = 0.247 (logrank); (c) p = 0.001 (logrank); (d) p < 0.001 (logrank); (e) p < 0.001 (logrank); (f) p < 0001 (logrank)

Table 2 shows that positivity for IA-2A, IA-2βA or ZnT8A tended to be associated with a higher 5-year progression rate to diabetes than positivity for IAA or GADA. Presence of IA-2A and/or ZnT8A identified a subgroup of first-degree relatives containing 77% of all individuals developing diabetes within 5 years, i.e. 10% more than with IA-2A alone (p = 0.031), and with an overall 5-year diabetes risk of 45%, which was similar to the risk observed for IA-2A alone or IA-2βA alone. Testing for IA-2βA in addition to IA-2A or to IA-2A and ZnT8A did not increase screening sensitivity, i.e. the additional IA-2βA-positive relatives did not develop diabetes within 5 years. Including GADA in the test panel improved the sensitivity of detecting prediabetes, but at the expense of an overall lower progression rate (Table 2).

Table 2 Influence of baseline antibody status and protective factors on 5-year progression rate and number of prediabetic first-degree relatives (Kaplan–Meier analysis)

Excluding first-degree relatives with protective HLA-DQ genotypes [18] or relatively protected offspring of diabetic mothers [20] reduced the antibody-positive group by 25% and 22%, respectively, but increased the 5-year progression rate of IA-2A+ or ZnT8A+ individuals by 3% and 7%, respectively (Table 2). Omission of both groups together increased 5-year risk by 10% as compared with that of the entire first-degree relatives group, while the group to be screened was reduced by 40%, without major decrease in the number of prediabetic first-degree relatives identified. In these subgroups of relatives, screening for IA-2A and ZnT8A (but not for IA-2A and IA-2βA) also tended to detect more individuals who developed diabetes within 5 years than screening for IA-2A alone (p = 0.063 to 0.031), but significance was only reached in the group without offspring of a diabetic mother (Table 2). In carriers of HLA-DQ2/DQ8, positivity for IA-2A or ZnT8A defined a group with a 74% 5-year risk, representing 38% of progressors within 5 years (Table 2). In all subgroups, the fraction of rapid progressors identified by screening for IA-2A and ZnT8A was independent of age, but three out of four were between 5 and 30 years (results not shown).

Cox regression analysis

The significant variables selected in a first model (see Methods and ESM Table 1) were tested in a second model additionally including combinations of antibodies. Positivity for HLA-DQ2/DQ8 and IA-2A or ZnT8A were selected as independent predictors of diabetes while having a type 1 diabetic mother proved protective (Table 3).

Table 3 Cox regression analysis in 409 first-degree relatives

Discussion

The present study is the first to examine the added value of measuring both IA-2βA and ZnT8A for prediction of impending diabetes in siblings or offspring of type 1 diabetic patients. It confirms the association of IA-2A [6, 9], IA-2βA [12, 21] and ZnT8A [13, 22] with rapid disease progression and demonstrates that IA-2A and ZnT8A represent the most sensitive combination of two markers to identify relatives with a high progression rate. Exclusion of relatives protected by HLA-DQ genotype [18] or by maternal diabetes [20] reduced the group to be followed by 40% without major loss in screening sensitivity. These findings are useful for screening programmes aiming to enrol high-risk individuals in immunointervention trials.

The strengths of this study are its registry-based nature, the use of a sensitive ZnT8A assay, the broad age range tested and the confirmation of glycaemic status at the last follow-up point for each relative. Criticisms might be that few relatives were followed for more than 10 years and that we did not test all relatives for IA-2βA and ZnT8A. However, the low prevalence of both antibodies in IAA, GADA and IA-2A first-degree relatives (0/441 in our own first-degree relatives, not shown; 0.2% for ZnT8A+ in a previous publication [22]) and of solitary ZnT8A+ in new-onset type 1 diabetes (4% of all patients in a previous study [13]; 2% in >400 patients of the BDR; M. Asanghanwa, unpublished results) supports the validity of our results for all first-degree relatives. Our baseline results cannot rule out the possibility that more first-degree relatives could have developed IA-2βA and/or ZnT8A at a later time point prior to clinical onset.

IA-2A, IA-2βA and ZnT8A target intracellular domains of antigens anchored in the membrane of beta cell secretory vesicles [2325]. These domains may only become visible to the immune system at the outer cell surface in the case of beta cell damage or dysfunction [13, 26, 27]. These antibodies appear generally later than IAA or GADA in prediabetes [13, 2830]. Our finding that IA-2A, IA-2βA and ZnT8A tend to cluster in first-degree relatives and predict, alone or in combination, rapid progression to diabetes is in line with the above observations, indicating that they preferentially mark the later stages of the preclinical phase of type 1 diabetes.

Based on our results, screening for IA-2A and ZnT8A emerge as the most sensitive combination of two markers to identify rapid progressors. Regardless of age, it identified 77% of them, i.e. 10% more than by IA-2A alone, without loss in overall 5-year risk (≥45%). Risk stratification may further be improved by genotyping the ZnT8-encoding SLC30A8 gene [22]. Unlike in the ENDIT study [12], IA-2βA did not improve risk assessment by IA-2A. This may relate to differences in inclusion criteria (registry-based vs preselection on basis of high-titre islet cell antigens) and the lower number of prediabetic participants in the present paper. Inclusion of functional tests may further improve screening specificity [2, 31]. Three out of four rapid progressors were between 5 and 30 years, which is the age group preferentially considered for participation in anti-CD3 trials [32]. Exclusion of genetically or maternally protected relatives increased screening efficiency for rapid progressors as the group to be followed was reduced by 40%, while 5-year risk increased to ≥55% for IA-2A+, IA-2βA+ or ZnT8A+ relatives with minor loss in sensitivity (<10%). There were too few progressors among offspring with a diabetic mother to establish their risk in case of positivity for IA-2A, IA-2βA, ZnT8A and susceptible genotypes, and so any decision on whether to exclude them or not would be ethically based only. Restricting follow-up to carriers of the HLA-DQ2/DQ8 high-risk genotype identified a group with about 75% 5-year risk but at the expense of >50% loss of screening sensitivity. The choice of screening options should be investigated in further studies.

In conclusion, we propose that immunointervention trials in pre-type 1 diabetes should focus on first-degree relatives without protective factors and testing positive for IA-2A and/or ZnT8A.