Introduction

Glycated haemoglobin A1c (HbA1c), which is formed through the non-enzymatic glycation of haemoglobin, is used as a cumulative estimate of mean blood glucose levels from the preceding 5–12 weeks in healthy people and in patients with diabetes [1]. It is a measure of glycaemic control commonly used as a predictor of diabetic microvascular complications in people with type 1 and type 2 diabetes [2]. As HbA1c measurement and patient self-monitoring of blood glucose are the two primary techniques used to assess the effectiveness of diabetes management [3], it is crucial to identify any factors that can bias these values. The recent changes in the ADA criteria for the diagnosis of diabetes [3] contribute to the importance of recognising any variation pertaining to HbA1c levels, especially in those aged 7 years and older, as this is a key age interval for the diagnosis of type 1 diabetes mellitus in paediatric population. Seasonal changes of HbA1c levels have been described in several studies assessing populations of adult patients with type 1 and type 2 diabetes [47]. Such fluctuations were also observed in children with type 1 diabetes [811]. However, most of the studies including paediatric age-groups had short observation periods or included relatively small groups of patients [8, 10], or were not primarily concerned with analysis of seasonal HbA1c fluctuations; these reported relatively small differences in HbA1c between the lowest and peak months [9, 11]. Factors related to potential seasonal variations in HbA1c levels in children have not been previously studied. In a recent study, including adult diabetes patients, Higgins et al. investigated the similarity of mean HbA1c levels with temperature changes in locations spread across the globe, and showed that the amplitude of both variables is strongly correlated [12].

The aim of our study was to evaluate seasonal HbA1c changes over more than 3 years (from August 2006 to October 2009) in a paediatric population of patients with type 1 diabetes. We also evaluated the relation of such changes with objective measures of weather conditions, which may reflect sun exposure and trends in physical activity.

Methods

Study population and laboratory data

The study was performed in Lodz, an administrative district in central Poland with approximately 2.6 million inhabitants. Data were collected from a central, certified laboratory database, which performs HbA1c testing for more than 90% of paediatric patients with diabetes in the region. Inclusion criteria were: (1) type 1 diabetes diagnosed on basis of clinical symptoms and laboratory results at onset; and (2) age below 18 years. Data from the laboratory database were cross-referenced by two independent researchers (A. Baranowska, W. Fendler) with the clinical database of patients with type 1 diabetes to verify patient status. Over the study period, August 2006 to October 2009, the laboratory method for assessment of HbA1c was consistent. HbA1c assays were performed by ion-exchange HPLC (Variant Hemoglobin A1c Program; Bio-Rad Laboratories, Hercules, CA, USA). The method used has been certified by the National Glycohemoglobin Standardization Program (NGSP) (http://www.ngsp.org/docs/methods.pdf; last accessed 25 November 2010). Reference values for healthy people estimated by the local laboratory were from 4.3 to 5.7%. The within-run CV determined by the manufacturer was 1.05% for people without diabetes and 0.94% for people with diabetes; the between-run CV was 1.61% and 1.16% respectively. Blood samples were collected with a HbA1c capillary collection system (Bio-Rad) and analysed within 2–6 days (according to the manufacturer’s manual, which specified not before the end of 24 h after blood collection to enable complete Schiff-Base removal). Specimens prepared using this procedure are stable for 2 weeks at room temperature or for 4 weeks at 2–8°C. Blood samples were not shipped. The Bio-Rad Variant device was located in the hospital building in a temperature-controlled laboratory room (ambient temperature maintained at 20–25° throughout study period). All HbA1c measurements were labelled by their original date. Laboratory personnel were not aware which results would be used in the study and were instructed not to perform any additional actions concerning HbA1c measurements during the study period.

Meteorological data

Meteorological data were provided by the Institute of Meteorology and Water Management station at Lodz-Lublinek Airport. Three weather-related variables were included in the analysis: (1) mean monthly ambient temperature in centigrade, measured using a meteorological station (MAWS301; Vaisala, Helsinki, Finland); (2) mean monthly duration of bright sunshine (hours of sunshine quantified as h per month); and (3) mean monthly solar irradiance (in J/cm2). The two latter variables were measured using a pyranometer (CM5; Kip&Zonen, Delft, the Netherlands). Both devices were compliant with the World Meteorological Organization and national standards.

Ethical approval

This study protocol was approved by the Local Ethics Committee of the Medical University of Lodz (No RNN/98/10/KE). Guardians of all study participants gave informed consent.

Statistical analysis

Fluctuations of mean HbA1c levels were analysed throughout the study period, divided into months. Following outlier correction procedure (upper and lower 2.5 percentiles of HbA1c), mean HbA1c values were used for seasonal variability evaluation through autocorrelations, a method that is aimed at internal pattern detection and based on mean correlation values between preceding and subsequent values of HbA1c measured in months distant from each other by set intervals. A gap in data for July 2007 (due to the nationwide physicians’ strike, during which no routine metabolic control visits were scheduled within the Clinic) was replaced by interpolated values of four adjacent values for the purpose of uninterrupted time–series analysis. Sun exposure data were also used for autocorrelation analysis. After comparison of autocorrelation patterns, patterns of metabolic control and meteorological data were evaluated using Spearman’s rank correlation with appropriate corrections for multiple comparisons. Statistical significance was set at p ≤ 0.05. Analyses were performed using Statistica 8.0 PL software (Statsoft, Tulsa, OK, USA).

Results

Of the 5,203 measurements performed in type 1 diabetes patients treated at the study centre, 3,935 were available for analysis from 677 children with type 1 diabetes after exclusion of measurements from patients with monogenic diabetes (MODY or permanent neonatal diabetes mellitus, n = 92), patients older than 18 years or those with diabetes duration of less than 12 months (n = 950), measurements done in July 2007 (time of Health Professional Strike) (n = 19) and outliers determined by values falling in the lowest and highest 2.5% (n = 207). Of these 677 children, 88 were of preschool age (<7 years old) and 589 were in primary and secondary school (≥7 years old). Clinical characteristics of the total study group and of the two subgroups are presented in Table 1. The mean (±SD) HbA1c levels for the whole study period was 7.65 ± 1.12%. The highest HbA1c levels were observed during February, November and December, while August and September showed the lowest HbA1c levels. The differences between mean HbA1c levels in late summer and winter months consistently exceeded 0.44%. The maximum difference between means of HbA1c in any two months of the study period was 0.72% (February 2008 vs September 2009). Variables of respective months including number of measurements, means and medians with 25 to 75% ranges are presented in the Electronic supplementary material (ESM), Table 1. A linear decreasing trend of HbA1c values over time was present throughout the study period, but did not reach statistical significance (R = −0.26; p = 0.11) and is shown in ESM Fig. 1. Increases and decreases of mean monthly HbA1c were inversely related to average sun exposure across months as shown in Fig. 1. Periodicity of cycles equalled roughly 12 months, suggestive of a positive correlation between months 1 year apart and a negative correlation for intervals of 6–7 months, reflecting seasonal variability. Autocorrelations of HbA1c declined in strength after two cycles below estimated white noise strength. Sun exposure and ambient temperature data showed a pattern with a periodicity of 11–13 months, with greater amplitude than HbA1c autocorrelations (ESM Fig. 2a–d). Mean HbA1c levels showed a sine-wave pattern of autocorrelations only in schoolchildren (superimposed autocorrelograms for HbA1c with ambient temperature and hours of sunshine per month in schoolchildren are shown in Fig. 2a). Such associations were not observed in preschool children, in whom HbA1c levels showed only white noise variability without any evidence of association with ambient temperature or sun exposure (Fig. 2b). In schoolchildren, correlation analysis showed a statistically significant negative correlation between HbA1c and sun exposure data by hours of sunshine (R = −0.52; p = 0.0007), solar irradiance (R = −0.52; p = 0.0006) and ambient temperature (R = −0.56; p = 0.0002) as shown in Fig. 3a–c. Such associations were not significant in preschoolers, where neither hours of sunshine (R = −0.10; p = 0.52), solar irradiance (R = −0.06; p = 0.69) or ambient temperature (R = −0.18; p = 0.29) correlated with HbA1c levels (Fig. 3d–f).

Table 1 Characteristics of the study group
Fig. 1
figure 1

Fluctuations of HbA1c and sun exposure (expressed as h of sunshine per month) during the study period. Continuous line, HbA1c; dotted line, sun exposure

Fig. 2
figure 2

Autocorrelation spectra of HbA1c (hatched) and h of sunshine per month (grey) in (a) school and (b) preschool children. Dashed lines, 95% CIs of white noise expected for such a sample size

Fig. 3
figure 3

Correlation of meteorological data, i.e. (a, d) h of sunshine per month, (b, e) solar irradiance and (c, f) ambient temperature, with HbA1c levels in schoolchildren (ac) and preschool children (df). a R = −0.52, p = 0.0007; (b) R = −0.52, p = 0.0006; (c) R = −0.56, p = 0.0002; (d) R = −0.10, p = 0.52; (e) R = −0.06, p = 0.69; (f) R = −0.18, p = 0.29

Discussion

This retrospective analysis of HbA1c levels in 677 children and adolescents with type 1 diabetes over more than 3 years revealed a clear seasonal variability with 11–13 months periodicity. Similar results have been described in patients with diabetes by others [411]. In our study, the lowest HbA1c levels were observed in August and September, while peaks were found in November, December and February. Poland, situated on the northern hemisphere, has a moderate climate, and these months correspond to late summer and late autumn/early winter respectively. Mean summer temperatures in Poland are between 16.5 and 20.0°C, mean temperatures in winter are between −6.0 and 0°C. In the present study, winter–summer differences in HbA1c levels for the study interval were from 0.44% in 2007 to 0.67 % in 2009. These values are high, compared with the results obtained by Gerstl et al. in a large cohort of young patients with diabetes observed over several years (mean values for September 7.86% vs January 8.08%) [11], but comparable to the observation of Sakura et al. in adult patients with diabetes, where the differences reached 0.45% [7].

Seasonal variations of HbA1c can be attributed to several factors, which may be classified into two groups: (1) sun exposure/temperature-dependent factors such as physical activity, vitamin D deficiency, serum melatonin concentration; and (2) sun exposure /temperature-independent factors such as school stress, amount of leisure time, mild upper respiratory and other ‘seasonal’ infections, seasonal dietary patterns related to the seasonal availability of ‘healthy food’ and endogenous seasonal hormonal variations (e.g. cortisol) [13]. Herbst et al., in their study including a cohort of 19,143 patients with type 1 diabetes (aged 3–20 years), showed that the frequency of regular physical activity is one of the most important factors influencing HbA1c levels and that HbA1c levels were lower in patients with greater regular weekly physical activity [14]. In countries with moderate climates, both younger (5–6 years old) and older children (10–12 years old) spend more time outdoors in warmer months than in cooler months [15]. Tucker et al. provided summaries of studies published up to 2009, which present evidence that quantitative measures of weather conditions correlate with physical activity of the general population, and particularly that of children [16]. Physical activity of schoolchildren, especially from urban areas, depends on after-school access to school gymnasiums and athletic fields, and on participation in organised activities. Sport facilities, however, are not well maintained, are poorly supported and access to them often costs money. Therefore the main kind of physical activity, especially for boys, are self-organised outdoor activities such as football or basketball, which are highly restricted by daylight hours. In the region covered by this study, about 35% of the population lives in rural areas. In rural areas the traditional ‘manual labour’ has become more mechanised and gradually requires less and less physical effort. Furthermore, on many farms, children with a chronic disease are still stigmatised as weaker, and as a result are often less burdened by physical work. There is also the global problem of increasingly sedentary lifestyles of children, resulting from the ever greater amounts of time spent in front of the computer, video games and television. Thus patterns of physical activity of children and adolescents from rural and urban areas are gradually becoming similar. Based on all these data, daylight hours and ambient temperature are good surrogates of physical activity in Poland, and so we examined their association with HbA1c levels. Our study found a significant negative correlation between HbA1c levels and ‘good’ weather conditions in schoolchildren, but not in preschoolers. Younger children, attending kindergarten or staying at home, are more physically active, both indoors and outdoors, than schoolchildren. During the school year, students spend many hours indoors, due to school work, which affects their daily physical activity patterns. For schoolchildren, weather conditions, day length and the structure of the school year (especially summer holidays) strongly influence outdoor activities, while preschoolers, who are active all day and unrestricted by school-associated workload, seem to be less sensitive to these influences. An association between HbA1c levels and quantified weather conditions, including sun exposure, has not been previously presented for children with diabetes. The data presented in this paper provide a physiological reason both for our findings of an annual pattern of HbA1c levels and for similar results from other groups.

The trend in HbA1c variation is consistent with seasonal variation of type 1 diabetes incidence [17, 18]. An inverse correlation with ultraviolet irradiance and diagnosis of type 1 diabetes has been observed in Newfoundland, Canada [19]. Both phenomena could be dependent on vitamin D synthesis in human skin [20]. The mechanism by which vitamin D possibly protects against type 1 diabetes is still unclear. Although there is some evidence suggesting that vitamin D may protect from cytokine-induced beta cell damage [21], such a mechanism provides little explanation of its association with metabolic control in stable, long-standing type 1 diabetes in children.

On the other hand, melatonin, which is secreted by the pineal gland, synchronises biological rhythms and follows a circadian and circannual rhythm. Melatonin is secreted in darkness and its level increases during the night [22] and in winter (December in the northern hemisphere) compared with daytime hours and summer (July) respectively [23]. Melatonin decreases insulin sensitivity, can decrease insulin secretion and can increase blood glucose levels, while phototherapy sessions (which potentially decrease melatonin levels) have been shown to induce recurrent hypoglycaemia in a patient with type 1 diabetes [2426]. From these observations, one could conclude that higher melatonin levels persisting for a longer time during autumn and winter nights could contribute to higher blood glucose levels. Such an association seems to be even clearer than any relation between vitamin D deficiency and higher blood glucose levels observed in colder months in patients with overt type 1 diabetes.

Other seasonal factors such as different dietary habits (e.g. lower consumption of vegetables and greater fat intake in winter), seasonal changes in body fat content (higher in winter and decreased in summer), more frequent infections and higher cortisol or glucagon levels in autumn and winter than in summer, as well as school stress can also unfavourably influence HbA1c levels through their effect on insulin sensitivity and blood glucose levels [5, 10, 11, 13, 27]. In Poland dietary habits may play an important role. Although food supply is generally unrestricted in the major food markets, due to the regional vegetation pattern prices of fresh vegetables and fruit increase during winter and early spring. In many households, especially those of poorer economic status, this may be compensated by higher intake of food rich in carbohydrates and fat. This can contribute to less favourable glycaemic profiles in young patients with type 1 diabetes who do not adapt their insulin doses appropriately.

Factors that can interfere with the performance of HbA1c assays are: haemoglobin variants and derivatives, shortened erythrocyte survival, vitamin C and E, iron deficiency anaemia, hypertriacylglycerolaemia, hyperbilirubinaemia, uraemia, chronic alcoholism, chronic ingestion of salicylates and opiate addiction (www.ngsp.org/factors.asp; accessed 25 November 2010). Two of these could be speculated to promote a seasonal pattern in a paediatric population, namely vitamin C (as potentially related to seasonal changes in vegetable and fruit intake) and hypertriacylglycerolaemia (related to increased fat intake). However, according to the Bio-Rad Variant Hemoglobin A1c manual (for hypertriacylglycerolaemia) and data from the literature (for pharmacological doses of vitamin C), none of these substances has any impact on HPLC-measured HbA1c levels [28]. In respect to meteorological factors, neither the NGSP (www.ngsp.org/factors.asp; accessed 25 November 2010) nor the manufacturer’s instruction manual mention ambient temperature, humidity or pressure as ‘factors that interfere with HbA1c results’.

Irrespective of their origin, the documented seasonal variations of HbA1c levels should be considered in disease management schedules: adaptation of insulin doses, meal plans, advising and organising physical activity, prevention of physical activity-related fear of hypoglycaemia in young patients with type 1 diabetes and the use of HbA1c as a tool for diabetes diagnosis. This may also be an important consideration for short-time (i.e. running over several months) clinical trials measuring before and after intervention HbA1c levels. In patients using continuous subcutaneous insulin infusion (in many paediatric centres, including ours, more than 50% of patients use this), the basal rate is assessed and adjusted if needed at least every 3 months and could anticipate the observed seasonality of glycaemic control.

There are some limitations that may bias the results of the study. According to ADA recommendations, the HbA1c test should be performed at least twice a year in patients with diabetes who are meeting glycaemic goals and at least four times a year in those who are not meeting treatment goals or whose therapy has changed [3]. The mean number of measurements among our study participants was 1.8 per patient per year, as some patients dropped out upon reaching the age of 18 years and some were evaluated only as inpatients during scheduled metabolic control assessment. However, a small number of measurements per patient reduces the effect of within-patient variability, as patients with better control have fewer observations with less variability and those with higher HbA1c and more frequent testing may weigh the analysis, causing false peaks of bad metabolic control. We would expect this effect to be randomly distributed across the year, thus not biasing any seasonal trend. Another limitation of our analysis is that HbA1c measurements may be affected by factors independent of seasonal variations. HbA1c measurements incorporate the potential biases of a different number of measurements per patient, this is in contrast to the measures of weather variation, which can be observed under consistent conditions. Thus, mean HbA1c levels may depend on the number of measurements each month and sampling of particular patients, while average sun exposure and temperature do not have such variation and allows for a continuous time series composed of the same number of measurements each month. A slight time delay in autocorrelations of HbA1c and weather conditions was also noted. It was probably due to the fact that the HbA1c level is an intrinsically time-lagged variable and depends strongly on the previous 3 months of glycaemic control. Autocorrelations of HbA1c were observed to dwindle faster, due to initial lower strength, which is typical for such autocorrelation effects. Although this may affect the direct relationship of seasonal changes and HbA1c, the established pattern is still valid.

The strength of this study is that it covered a well documented, ethnically homogeneous (white) population with several long-running epidemiological projects such as the developing Nationwide Registry of Pediatric and Adolescent Diabetes and the Polish Registry of Neonatal Diabetes [29, 30]. From a geographic perspective, the region is characterised by the lack of mountains or seaside providing a stable temperate climate with daylight hours of sunshine varying continuously throughout the year and typically no natural disasters.

In conclusion, periodic changes of HbA1c are a biologically significant phenomenon in young patients with type 1 diabetes and should be considered in patient education and diabetes management schedules. They may affect the results of paired, time-dependent comparisons in clinical trials using HbA1c levels as their primary outcome and the HbA1c-based diagnosis of diabetes.