Introduction

The incidence of type 1 diabetes diagnosed before 15 years of age has increased in most studied populations over the past two to three decades, with an average increase per year of approximately 2.5–3.0% during the 1990s and 2000s [13]. The incidence has tended to rise more steeply, in relative terms, in some of the low prevalence countries, but there are wide variations in the rates of change between countries and between periods [36]. The most recent trends from Sweden and Finland, the countries with the highest incidence, have indicated that the incidence has ceased to increase in the past few years [7, 8].

After Finland and Sweden, Norway is among the countries with the highest incidence of childhood-onset type 1 diabetes in the world [2, 4]. Data published previously showed a stepwise but overall increasing trend during the period 1973–2003, increasing from 19.1 per 100,000 person-years in 1973 to 28.9 per 100,000 person-years in 2001–2003 [2, 9]. Regional variation within Norway has been reported [9, 10].

The aim of this study was to describe the incidence rates of type 1 diabetes with onset before 15 years of age in Norway in the period 1989–2012 and to assess the regional variation within the country for the period 2004–2012. Furthermore, we took advantage of the nationwide Norwegian Prescription Database (NorPD), established in 2004, as a secondary source of cases, allowing estimation of the completeness of ascertainment.

Methods

The Norwegian Childhood Diabetes Registry

In the Norwegian healthcare system all children aged 0–14.9 years with suspected diabetes are referred to a paediatric department. The Norwegian Childhood Diabetes Registry (NCDR) includes all new cases of childhood-onset diabetes, reported from all the paediatric departments in Norway, based on informed consent from the child and/or the child’s parents. Cases are included as type 1 diabetes in the NCDR based on a clinical diagnosis of type 1 diabetes, using the first insulin injection as the date of diagnosis, and residence in Norway at the time of diagnosis, in accordance with the EURODIAB criteria.

NorPD

The NorPD is a national population-based prescription drug registry established at the Norwegian Institute of Public Health in 2004. It captures all prescriptions dispensed at pharmacies to individual patients in open care, and all pharmacies in Norway are obliged by law to forward prescription data to the NorPD [11]. The first dispensation of insulin (insulins or insulin analogues, ATC-code A10A) from a pharmacy was assumed to be near the time of diagnosis of type 1 diabetes for the purpose of estimating completeness of ascertainment during the period 2005–2008 (excluding prevalent users in 2004 to ensure that we only included incident users of insulin from 2005 to 2008).

Statistical analysis

The size of the population in each calendar year, for each sex and 1 year age group, was obtained from Statistics Norway (www.ssb.no) and used to estimate the number of person-years in each subgroup.

The two registers were linked via the 11 digit personal identification code given to each Norwegian citizen; the number of children recorded only in the NCDR, only in the NorPD, or in both, was counted for each calendar year from 2005 to 2008. Completeness of registration in the NCDR was assessed by the capture–recapture method, which requires that the two sources are independent; none of the sources, however, have to be complete [12]. Formulae for ascertainment and 95% CIs are provided in the footnote of Table 1. For the calculation of incidence rates, only cases recorded by the primary register (NCDR) were included.

Table 1 Number of newly diagnosed cases of type 1 diabetes in children aged 0–14.9 years, registered in the NCDR and NorPD

The effects of age, sex and calendar periods were modelled using Poisson regression with Stata software, version 12 (StataCorp LP, College Station, TX, USA). Various models were fitted, including sex, age group, calendar period and county of residence at diagnosis. Age was modelled as a categorical variable grouped in 5-year intervals for most analyses. We also assessed the influence of modelling age in 15 age groups spanning 1 year each. The calendar period was modelled with linear splines to allow estimation of slopes (change per year) and changes in slopes between time segments. Visual inspection of plots of yearly rates clearly suggested three time segments. This was confirmed by plotting expected rates from models with restricted cubic splines using five knots and Harrell’s rule for placement of knots (default option in Stata’s mkspline procedure for restricted cubic splines) (electronic supplementary material [ESM] Fig. 1). The latter model is more flexible. We also explored two-way interactions among the four study variables (using 5-year categories for age). Likelihood ratio tests were used to test the influence of region, sex and age.

The study was approved by the Regional Ethics Committee and the National Data Inspectorate.

Results

Time trend in incidence, 1989–2012

Figure 1 shows the observed incidence rates for each of the years between 1989 and 2012, indicating three distinct segments with different slopes. From 1989 to 1996, the average incidence rate was 22.6 per 100,000 person-years (95% CI 21.4, 23.7), and the estimated change per year was 1.8% (95% CI −0.07, 3.6; p = 0.059). From 1996 to 2004, the average incidence rate was 28.4 per 100,000 person-years (95% CI 27.3, 29.6), with an estimated change per year of 3.4% (95% CI 2.2, 4.7; p < 0.0001). From 2004 to 2012, the average incidence rate was 32.7 per 100,000 person-years (95% CI 31.5, 34.0), with an estimated change per year of 0.3% (95% CI −0.9, 1.6; p = 0.64). Finally, we found that the difference in slope between the first two time segments was not significant (p = 0.24), while the change per year during 2004–2012 was significantly different from the change per year during 1996–2004 (p = 0.006).

Fig. 1
figure 1

Yearly observed type 1 diabetes incidence rates with 95% CIs for children in Norway during 1989–2012. The broken line represents the expected rates from Poisson regression models with linear splines and knots at 1996 and 2004

There was an unexpected single peak in incidence in 2009 (Fig. 1). We explored whether this was consistent over age groups (ESM Fig. 2a), sex (ESM Fig. 2b) and county of residence (not shown) but found no consistent pattern. Slight differences in observed trends among sex, age group and county (none of which were statistically significant) resulted in the observed pattern in Fig. 1. Furthermore, formal modelling of the yearly incidence in the period 2004–2012 (eight degrees of freedom test) was not significant (p = 0.74). When considering 2008, 2009 and 2010 in isolation, 2009 was not significantly different from either of the two neighbouring years (p = 0.09, p = 0.19, respectively). Collectively, these results support the notion that the observed 2009 peak might have been a chance finding.

Variation by age, sex and region, 2004–2012

Because the incidence by age, sex and county for periods up to 2003 has been published previously [2], the remainder of the results are presented for the period 2004–2012 (data in Table 2 and ESM Table 1). The incidence was slightly higher in boys (33.9 per 100,000 person-years; 95% CI 32.2, 35.7) than in girls (31.4 per 100,000 person-years; 95% CI 29.7, 33.2) (Table 2). After adjustment for age and period, the boys had an estimated 8% higher incidence compared with girls (95% CI 0.1, 16; p = 0.047). In both sexes, the highest incidence was at age 10–14 years (ESM Table 1). The age-specific incidence rates during 2004–2012 were 19.2, 35.2 and 43.1 per 100,000 person-years, for the 0–4, 5–9 and 10–14 year olds, respectively (ESM Table 1). The incidence started out very low for those below 1 year of age and increased rapidly to peak at around 10–11 years of age. There was an apparent drop after 11 years of age in girls but not in boys (ESM Fig. 3). No significant two-way interactions among sex, age and period were found in Poisson regression modelling, suggesting that the above-mentioned results for period, sex and age were consistent in subgroups. Also, goodness-of-fit tests (deviance and Pearson) suggested adequate fit of models without interactions.

Table 2 Age- and sex-specific incidence rates of type 1 diabetes in children aged 0–14.9 years in the NCDR, 2004–2012

We observed a significant variation among the 19 counties in Norway (likelihood ratio test with 18 degrees of freedom, after adjustment for age and sex; p < 0.001). The county-wise incidence rates with 95% CIs and change per calendar year are shown in Table 3, and rates are shown in a colour-coded map in ESM Fig. 4. Aust-Agder, in the southern part of Norway, had the highest incidence at 45.7 per 100,000 person-years. Oslo had the lowest incidence at 22.3 per 100,000 person-years. Poisson regression modelling suggested no significant county by calendar period interactions, suggesting that the time trends were relatively homogeneous among counties. When analysing percentage change per year the largest observed increase was in Finnmark (8.1%; p = 0.194), and the largest decrease was observed in Østfold (−6.1%; p < 0.053); neither were significant.

Table 3 Age- and sex-standardised incidence rate of childhood-onset type 1 diabetes in the 19 counties of Norway, 2004–2012

Ascertainment, 2005–2008

During 2005–2008 the NorPD contained data on first-time insulin dispensed for 115 subjects not reported to the NCDR. The completeness of ascertainment in the NCDR was calculated to be 91% for this study period (Table 1). The completeness of ascertainment in the NorPD was calculated to be 98.4%, and the total ascertainment from both sources in the period 2005–2008 was 99.8%.

Regarding the 51 cases identified exclusively from the NorPD, 25 (49%) of the cases were clustered at the beginning of 2005. This might reflect that these patients had their clinical diagnosis in 2004. Regarding the 29 cases exclusively identified from the NorPD in 2008, 16 (55%) were clustered at the end of 2008. The cases exclusively found in the NorPD (n = 115) had a mean age of 10.7 years (SD 3.6) and 52% were boys. They were spread among 17 of 19 counties, but 63% (n = 72) were clustered in five of 19 counties: Hordaland, Rogaland, Akershus, Vest-Agder and Sør-Trøndelag. This might reflect local difficulties as a delay in registration of new cases of type 1 diabetes in the NCDR.

The cases exclusively registered in the NCDR (n = 20) had a mean age of 9.6 years (SD 4.3) and 30% were boys. They came from eight out of 19 counties, and 40% (n = 8) were clustered in one county (Telemark).

Discussion

Our results show that the incidence rate among Norwegian children was stable at about 33 per 100,000 person-years during 2004–2012, and that the time trend in this period was significantly different from that during 1996–2004. Our nationwide registry has registered incident cases prospectively since 1989 using standard criteria and methods. Other notable results include significant regional variation, discussed in more detail below. We have also estimated ascertainment in the NCDR using the nationwide NorPD as an independent source for the period 2005–2008.

Regional variation has previously been described in Norway [9, 10]. The current data show a slightly different pattern than in previous years. While the southernmost county in Norway, Vest-Agder, previously had the highest incidence, and the neighbouring Aust-Agder county had an incidence near the country average [10], the current pattern is reversed. Oslo had the lowest incidence; in the time period 1989–1998, however, Oslo had a low, but not the lowest, incidence [10]. The relatively small population sizes make it difficult to study potential causes in any county. The reason for the low incidence in Oslo is also unknown. Compared with other counties, however, Oslo has a relatively high proportion of ethnic groups known from other studies to have a lower incidence of type 1 diabetes. Finnmark, the northernmost county that had the lowest incidence in 1989–1998, is still among the counties with the lowest incidence. Our present data are based on county of residence at the time of diagnosis. We were not able to take into account whether or when some individuals might have changed county of residence before diagnosis; however, most people in Norway tend to move within counties when they move house. Future analyses should investigate small area variations and the potential roles of factors such as population density, urban/rural setting, latitude and climate.

The strength of this study is the nationwide, population-based design. The linking of two independent population-based registries ensured the capture of cases with new-onset type 1 diabetes.

The study disclosed that some of the cases in the NCDR were registered with an incorrect personal identification number, and clerical errors may lie behind the less than 100% estimated completeness of ascertainment. This provides an important impetus for us to continue to work with quality assurance at all levels of our registry. It is also possible that a small number of children with prevalent type 1 diabetes might have emigrated from another country to Norway and incorrectly been counted as incident type 1 diabetes in the NorPD. The latter would thus have biased the estimated completeness of the NCDR downwards. In future studies, we hope to be able to correct any such bias.

In this study we did not include children with type 2 or other forms of diabetes in the NCDR, although a few cases in the NorPD might have been included. This limitation of not removing the small number of cases in the NorPD where insulin was being given for conditions other than type 1 diabetes (e.g. type 2 diabetes, cystic fibrosis or steroid-induced diabetes) might have led to an underestimation of the completeness of the NCDR.

There are clear indications of geographic differences in the trends in incidence of childhood-onset type 1 diabetes, but the overall annual increase is estimated to be around 2.5–3% [3, 13, 14]. Several European studies have suggested, in relative terms, that increases are greatest among children aged 0–4 years [5, 1517]. This has never been the case in Norway, where the highest increase has always been in the age group 10–14 years.

Our result is in accordance with the most recent observations from Sweden [7] and Finland [8], as well as from some lower incidence countries such as the Czech Republic [18].

This 9-year period of no apparent increase in incidence suggests that time trends in exposure to unknown factors conferring an increased risk of type 1 diabetes have levelled off, or that factors conferring protection have ceased decreasing in Norway. As non-genetic factors able to explain changing trends in the epidemiology of type 1 diabetes have not been established, we find it difficult to speculate about the explanations, but we hope that these observations can stimulate generation of new hypotheses to be tested in future aetiological studies.