Introduction

In addition to the well-known and well-described microvascular complications of diabetes, i.e. nephropathy, retinopathy and peripheral neuropathy, a condition that may be tentatively referred to as ‘diabetic encephalopathy’ is attracting more attention. This term was introduced in 1965 [1] to describe the pathological changes in the central nervous system of young patients with long-term type 1 diabetes mellitus. The relationship between the development of this encephalopathy and exposure to chronic hyperglycaemia, leading to microangiopathy, is not well established. A recent meta-analysis concluded that lowered cognitive performance of patients with type 1 diabetes appears not to be associated with recurrent episodes of severe hypoglycaemia, as often speculated, but rather with microvascular complications such as retinopathy [2]. A longitudinal study [3] showed that cognitive efficiency may decline over time in diabetic adults, and that this neurocognitive change may be linked, at least in part, to the occurrence of microvascular complications such as proliferative retinopathy and to elevated blood pressure. Other evidence for a damaging effect of chronic hyperglycaemia on the brain comes from the work of Ferguson et al. [4]. Subjects who had background retinopathy performed less well across most cognitive domains examined. The authors also found that background retinopathy was associated with small punctate white matter lesions corresponding to enlarged perivascular spaces.

To date, the few neuroradiological studies that have assessed structural brain abnormalities in patients with type 1 diabetes have demonstrated conflicting results [48]. The majority of these studies have been based on manual or semi-automated region-of-interest-guided measurements, which may be inherently biased.

This prompted us to employ a fully automated, sensitive, standardised, whole-brain technique known as voxel-based morphometry (VBM) [9], applying it to magnetic resonance images to study structural differences in cortical grey matter density (GMD) in patients with type 1 diabetes. Our aim was to establish whether long-term hyperglycaemia resulting in advanced retinopathy was associated with reduced GMD in the brain. We used advanced retinopathy as a marker of microvascular disease because the retinal vessels share common anatomical, embryological and physiological characteristics with the cerebral microcirculation [10, 11]. Abnormalities of the retinal arterioles (e.g. microaneurysms, retinal haemorrhages and arteriolar narrowing) may be markers of concomitant cerebral small vessel diseases [12].

We therefore included patients with type 1 diabetes and proliferative retinopathy, and compared these patients with type 1 diabetic patients who did not have retinopathy in order to assess the relation between microvascular changes and GMD. Both patient groups were compared with healthy controls to assess the effect of diabetes in itself.

We hypothesised that patients with retinopathy have reduced GMD compared with patients without retinopathy and healthy controls.

Subjects, materials and methods

Subjects

Three age- and sex-matched groups of participants were included: 18 patients with type 1 diabetes (World Health Organization criteria, 1999) without microvascular complications (non-diabetic retinopathy [NDRP], maximum of three microaneurysms), 13 patients with type 1 diabetes with a microvascular complication, i.e. diabetic proliferative retinopathy (DRP) (grade 5 DRP according to the EURODIAB classification [13]) and 21 healthy controls (Table 1). All subjects were right-handed and normotensive (<140/90 mmHg). None of the NDRP patients had a micro- or macrovascular diabetic complication. Four DRP patients were known to have nephropathy and two other patients were known to have nephropathy and neuropathy. Exclusion criteria were: previous alcohol or drug abuse; a history of psychiatric disease/treatment; a history of severe head trauma accompanied by loss of consciousness; stroke; epilepsy; and a history of severe, recurrent hypoglycaemia (defined as more than five episodes that required external assistance for recovery) [14]. Written consent was obtained from all subjects and the study was approved by the local medical ethics advisory committee. The investigations were carried out in accordance with the principles of the Declaration of Helsinki as revised in 2000.

Table 1 Characteristics of participants

Magnetic resonance image acquisition

Imaging was performed on a 1.5 T Siemens Sonata (Siemens, Erlangen, Germany) scanner using a standard circularly polarised head coil, with foam padding to restrict head motion. A localiser scan was first performed for positioning of the image planes, followed by an automated shim procedure to improve magnetic field homogeneity. Scans were obtained as whole-brain T1-weighted MPRAGE (magnetisation prepared rapid acquisition gradient echo) volumes and were acquired in the coronal plane (inversion time 950 ms, repetition time 2,700 ms; echo time 5.15 ms; flip angle 8°; 160 slices; voxel size 1×1×1.5 mm).

This imaging protocol was part of a functional magnetic resonance imaging (MRI) study which did not include a fluid-attenuating inversion recovery sequence. For this reason, white matter lesions could not be scored.

Statistical analysis

Statistical analysis was performed using SPSS version 11.0 (SPSS, Chicago, IL, USA). Demographic data of patients and healthy controls were analysed using one-way ANOVA for continuous variables and the χ 2 test for binomial variables.

MRI data analysis

Data analysis was performed using statistical parametric mapping software (SPM2; Wellcome Department of Cognitive Neurology, Institute of Neurology, London, UK, http://www.fil.ion.ucl.ac.uk/spm/) with the VBM tool Jena script http://dbm.neuro.uni-jena.de/vbm.html; both websites accessed on 9 January 2006) in MATLAB version 6 (Mathworks, Natick, MA, USA).

We used optimised VBM with a custom template and priors because this method improves the plausibility of intergroup comparisons, presumably because of improved segmentation and spatial normalisation [15, 16].

Briefly, the processing steps outlined by this protocol are: (1) the creation of a customised anatomical T1-weighted template and prior probability images separately for each group, by normalising the brain images to the default SPM T1-weighted template, segmenting, averaging and smoothing the averaged brains in each group; and (2) the normalisation of the structural brain images in each group using these customised templates, segmenting and cleaning the original T1-weighted images. This procedure has been described in detail previously [16].

Statistical analysis used the general linear model and is based on the random Gaussian field theory [17]. We performed analyses on the smoothed grey matter images only. Differences in GMD between the three groups were examined using one-way ANOVA and post hoc comparisons. We calculated the positive contrasts (controls > NDRP + DRP; controls > NDRP; controls > DRP and NDRP > DRP) and the negative/inverse contrasts. Results were thresholded at p<0.05 corrected for multiple comparisons at the cluster level.

This analysis was based on a strict p value (p<0.05, cluster-corrected). In a second analysis, we restricted the analysis to regions of interest which followed from one of the contrasts (see Results), using less stringent statistics (p<0.001 uncorrected). Thus, this procedure limited statistical comparisons in SPM to the desired anatomical region and increased sensitivity (lower p value) while controlling for false positives (restricted to certain regions).

Additionally, multiple regression analyses between GMD and HbA1c, diabetes duration, onset of diabetes and blood pressure were performed, adjusting for age, to identify clusters of voxels whose GMD related to these parameters. A threshold of p<0.05 was used, corrected for multiple comparisons at the cluster level.

Results

Controls vs NDRP and DRP groups

There were no significant differences in GMD between healthy volunteers and type 1 diabetes patients (NDRP and DRP groups combined).

NDRP group vs DRP group

Reduced GMD was observed in the DRP group compared with the NDRP group in the left middle frontal gyrus (Z=3.86), right inferior frontal gyrus (Z=4.43), right occipital lobe (Z=4.36) and left cerebellum (Z=4.08) (Fig. 1). Table 2 lists the stereotaxic coordinates of the most significant clusters of voxels.

Fig. 1
figure 1

Voxel-based morphometry comparisons between 13 patients with type 1 diabetes and proliferative retinopathy (DRP) and 18 patients with type 1 diabetes without retinopathy (NDRP). Right and left are inverted in the image, so the right side of the brain appears on the left side in the image. Changes in GMD are shown on a mean normalised structural image in axial planes. Significantly reduced GMD is seen in the right frontal gyrus (marked A), right occipital lobe (marked B), the left cerebellum (marked C) and left middle frontal gyrus (marked D) in the DRP group compared with the NDRP group (p<0.05 corrected for multiple comparisons)

Table 2 Areas of grey matter density reduction in the DRP group vs the NDRP group

Differences between healthy volunteers and the NDRP group and the DRP group in these same regions were not statistically significant, although GMDs tended to be reduced in the DRP group relative to controls.

In the second analysis we calculated the following contrasts: controls > NDRP (mask NDRP > DRP) and controls > DRP (mask NDRP > DRP). The former contrast did not show significant differences in GMD between the groups. The latter, however, displayed significant reductions at the voxel level (p<0.001) of grey matter in the DRP group compared with controls in the right inferior frontal gyrus (Z=3.19) and right occipital lobe (Z=3.77). Examination of the inverse contrasts, i.e. reduced GMD of controls compared with DRP and NDRP and reduced GMD of NDRP compared with DRP, did not yield any significant results for the analyses mentioned above.

We observed no significant association between atrophy and HbA1c, diabetes duration, onset of diabetes or blood pressure within the type 1 diabetes group. These variables did not account for the differences in GMD between the DRP and the NDRP groups.

Discussion

We found that the DRP patients, compared with the NDRP patients, showed reduced GMD in the right inferior frontal gyrus, right occipital lobe, left cerebellum and left middle frontal gyrus. Two of these regions, the right inferior frontal gyrus and the right occipital lobe, also showed reduced GMD in comparison with the healthy participants. The NDRP group did not show reduced GMD compared with healthy participants.

Central and peripheral atrophy in type 1 diabetes [7, 8, 18] and subcortical white matter lesions [4, 5] have been reported in the literature. The methods used in these studies are semiquantitative. In the study by Ferguson and colleagues [4], a relationship with retinopathy was observed. However, they used several scoring systems to capture the intensity, distribution and appearance of white matter lesions, because no single scale was judged to be an adequate summary. Moreover, small punctate white matter lesions, which represent enlarged perivascular spaces, were scored on a scale of 0–3, because these lesions were not accounted for by any of the above rating scales.

The strength of our study lies in the inclusion of matched non-diabetic controls and the use of a sensitive, fully automated whole-brain technique that detects regionally specific differences in brain tissue composition on a voxel-by-voxel basis. At its simplest, it involves a voxelwise comparison of the local concentration of grey matter between groups of subjects [9].

Recently, an interesting paper on the effects of type 1 diabetes on GMD as measured by VBM was published [19]. This study investigated whether lower grey matter densities in patients with type 1 diabetes were present, and if so, whether they were associated with glycaemic control and/or a history of severe hypoglycaemic events. It was found that, compared with healthy controls, patients with type 1 diabetes showed lower GMD in several brain areas. One or more severe hypoglycaemic events was associated with less GMD in the left cerebellar posterior lobe. Moreover, the higher the HbA1c level (suggesting compromised glycaemic control) the lower the GMD concentration, including the left posterior cingulate, right parahippocampal gyrus, left hippocampus, left superior temporal gyrus, right occipital cuneus and left thalamus. In an additional (post hoc) analysis, the authors found that a greater degree of retinopathy was associated with greater GMD loss in brain regions used for cognition (frontal and temporal regions). It was concluded that persistent hyperglycaemia may be associated with lower levels of GMD in patients with type 1 diabetes and that the superior temporal gyrus is particularly vulnerable to the effects of diabetes. Although this study and our study found a relation between retinopathy and the degree of density reduction, there are important differences in design and patient population that complicate comparison of results. Whereas our study was designed as a case–control study to test the effect of microangiopathy on the brain, the study by Musen et al. [19] investigated whether GMD reduction occurred in type 1 diabetic patients, and, if so, whether hypoglycaemia and/or longstanding hyperglycaemia were associated with reduced GMD. Overall, the effects we found were less significant, probably because our initial analysis was based on a stricter statistical threshold (p<0.05 corrected for multiple comparisons). Lowering our threshold would probably reveal more significant reductions in GMD. We did not indicate the magnitude of the GMD reduction, but comparing Z-scores indicates that the effects found are in fact of the same magnitude (i.e. both studies found subtle effects).

Our results support the notion that type 1 diabetes is associated with modest regional cortical atrophy, but these changes are limited to patients with advanced microvascular complications of diabetes. From the results of the regression analysis (no significant association between atrophy and HbA1c, diabetes duration, onset of diabetes or blood pressure within the type 1 diabetes group and between the DRP and NDRP groups), one might conclude that the GMD differences found in patients with type 1 diabetes are associated with long-term hyperglycaemic end-organ damage instead of more short-term indices of hyperglycaemia itself, such as HbA1c. The question arises whether some individuals are more sensitive to hyperglycaemic damage (and will therefore develop retinopathy and brain atrophy), or that the GMD reductions we see are indeed the result of prolonged exposure to hyperglycaemia in the past. To answer this question, a study with a longitudinal design is warranted.

Our results indicate that DRP patients, compared with NDRP patients, showed reduced GMD in four brain areas (left middle frontal gyrus, Z=3.86; right inferior frontal gyrus, Z=4.43; right occipital lobe, Z=4.36; left cerebellum, Z=4.08).

The inclusion of a matched healthy control group is an important strength of our study, allowing us to test whether reduced GMD is associated with diabetes in itself or with the presence of microvascular disease. If reduction of GMD occurs as a result of diabetes (other than microvascular disease), one would expect that the differences in GMD reduction between the DRP and the NDRP group would be less pronounced than those between the DRP group and healthy controls. In a second analysis we investigated whether the two patient groups differed from the healthy controls, applying a less stringent threshold of p<0.001 (uncorrected for multiple comparisons; a widely accepted and used statistical threshold) to increase sensitivity to group differences, while restricting our analysis to predefined regions of interest. Only voxels that were found to be ‘active’ in a previous contrast were queried (in this case, we used the left middle frontal gyrus, the right inferior frontal gyrus, the right occipital lobe and the left cerebellum as regions of interest for further analysis). We found significant reductions at the voxel level (p<0.001) of GMD in the DRP group compared with the healthy participants in the right inferior frontal gyrus (Z=3.19) and the right occipital lobe (Z=3.77). We did not find differences between the healthy participants and the NDRP group in this analysis. The Z-scores in the NDRP–DRP comparison (initial analysis) are somewhat higher than the Z-scores in the comparison of the DRP group with the healthy participants (the second, less stringent analysis), The Z-scores of this latter comparison did not reach significance with the initial, stringent thresholding, but were not statistically lower than the Z-scores of the DRP–NDRP comparison (which did reach significance in the initial analysis). Hence, we can conclude that the DRP group showed significantly reduced GMD compared with the NDRP group and compared with the healthy controls, and that the differences between the DRP and the NDRP groups were not greater than the differences between the DRP group and the healthy participants. Furthermore, because we did not find differences between the NDRP group and the healthy participants, we conclude that the reduction in GMD is not related to diabetes in itself, but rather has a microvascular origin.

The mechanism behind the association between retinopathy and reduced GMD is unknown, but hyperglycaemia is likely to be an important mediator in both. Hyperglycaemia may lead to an accumulation of potentially toxic glucose metabolites, oxidative stress, accelerated formation of advanced glycation end-products and microvascular changes in the brain, analogous to those observed in peripheral neuropathy [20]. Considering vascular changes in the brain in diabetes, it is known that structural abnormalities at the microvascular level include thickening of capillary basement membranes and decreased capillary density [21]. During ageing, brain capillaries undergo progressive degeneration, as evidenced by amyloid deposition, thickened basement membranes and reduced vessel elasticity [22, 23]. Considering the similarities of these events, it is possible that diabetes and ageing share the same pathogenic processes. Since only the DRP group showed a reduction in GMD, it is likely that hyperglycaemia-induced changes contribute to accelerated ageing of the brain, as reflected in increased cortical atrophy. However, there is evidence that accelerated ageing is associated with global grey matter volume decreases rather than only with local areas of accelerated loss [15, 24]. The fact that we did not find a global ageing effect can probably be attributed to below-threshold effects, but this needs to be clarified in future research.

Another possible contributing mechanism underlying the diminished GMD in the DRP group may be (recurrent) microinfarction in the grey matter, resulting in its loss. Microinfarction has been reported as the basis of dementia [25, 26]; it is associated with a history of transient ischaemic attacks [27] and contributes significantly to the progression of cognitive deficits in brain ageing [28]. The inferior frontal gyrus and the occipital lobe, showing decreased GMD in the DRP group compared with the NDRP group and the healthy controls, are the watershed areas of the medial cerebral artery and the posterior cerebral artery. Watershed areas are located peripherally to primary distribution areas of the vascular systems. Because they are located at the extreme ends of the arterial distribution, they are particularly vulnerable to ischaemia and infarction in people with circulation problems. It is well established that vascular watershed zones are the first to be deprived of sufficient blood flow in the event of cerebral hypoperfusion [29]. Therefore, we suggest that the restriction of cortical microinfarcts in our population to watershed cortical areas points to cerebral hypoperfusion (resulting from hyperglycaemia-induced vascular abnormalities) as the determining factor in their pathogenesis. It is possible that diabetes patients with a microvascular complication have an additional risk of developing a reduction in GMD, since both ageing and hypoperfusion can affect the integrity of the grey matter (Fig. 2).

Fig. 2
figure 2

Possible mechanism underlying reduced grey matter density in the DRP patient group

Further information about the occurrence and severity of watershed cortical infarcts, and hypoperfusion as a possible determinant of these infarcts with respect to the progression and severity of diabetes is needed.

The possible impact of isolated microvascular lesions on cognition remains controversial, and data regarding this issue are very scarce. Although some studies point to a possible causal relationship with dementia [3033], others challenge this point of view [34]. Shrinkage and atrophy of the brain are generally considered to be a risk factor for cognitive decline [3539].

Hypertension has been associated with changes in brain structure and cognitive function [4043]. It is known that hypertension accelerates the decline in cerebral perfusion [44, 45]. Altered regional cerebral bloodflow in hypertensive subjects in the middle/posterior watershed areas has been found [45]. Moreover, memory performance in this group was related to this diminished response, so hypertension represents a risk for generating watershed cortical infarcts. Therefore, it is of importance to control for hypertension when studying brain structure, as we did.

In summary, our data show that patients with type 1 diabetes who had developed advanced retinopathy also had modest focal cortical atrophy on MRI of the brain. Prospective studies are warranted to establish whether reduced GMD in patients with diabetes is progressive, and how it affects various domains of cognitive functioning.