Altered Functional Connectivity Strength of Primary Visual Cortex in Subjects with Diabetic Retinopathy.

OBJECTIVE: The purpose of the study was to find the differences in intrinsic functional connectivity (FC) patterns of the primary visual area (V1) among diabetic retinopathy (DR), diabetes mellitus (DM), and healthy controls (HCs) applying resting-state functional magnetic resonance imaging (rs-fMRI). PATIENTS AND METHODS: Thirty-five subjects with DR (18 males and 17 females), 22 DM (10 males and 12 females) and 38 HCs (16 males and 22 females) matched for sex, age, and education underwent rs-fMRI scanning. Seed-based FC analysis was performed to find the alterations in the intrinsic FC patterns of V1 in DR compared with DM and HCs.
RESULTS: The study found that DR patients had a significant lower FC between the bilateral calcarine (CAL)/left lingual gyrus (LING) (BA 17/18) and the left V1, and between the bilateral CAL/left LING (BA 17/18) and the right V1 compared with the HCs. Meanwhile, patients with DR exhibited higher FC strength between the left V1 and the bilateral Caudate/Olfactory/Orbital superior frontal gyrus (OSFG), and between the bilateral Caudate/Olfactory/OSFG (BA 3/4/6) and the right V1. Compared with DM group, patients with DR showed increased FC strength between the right CAL (BA 17/18) and the right V1. DM group exhibited lower FC strength between the left fusiform and the left V1, and between the bilateral CAL and the right V1 when compared with HCs. Moreover, DM group was observed to have higher FC strength between the left superior frontal gyrus and the left V1.
CONCLUSION: Our findings indicated that DR patients exhibited FC disruptions between V1 and higher visual regions at rest, which may reflect the aberrant information communication in the V1 area of DR individuals. The findings offer important insights into the neuromechanism of vision disorder in DR patients.

Introduction

Diabetic retinopathy (DR), known as the main cause of adult acquired blindness, is the most common retinal microvascular complication of diabetes mellitus (DM). DR is a progressive disease and can be categorized as an early stage (i.e., nonproliferative DR) and a more advanced stage (i.e., proliferative DR). Multiple risk factors are associated with the progression of DR, such as DM duration, chronic hyperglycemia, hypertension, and dyslipidemia.1,2 With prolonged progression of DM, approximately 50% of affected patients develop DR within 10 years after diagnosis of DM, and this proportion can reach 90% after 25 years.3

Intra-retinal microvascular abnormalities are part of the cause of visual impairment in patients with DR. There is increasing evidence that the vision loss in DR patients is also associated with pathological changes in various parts of the visual pathway, particularly the primary visual area (V1).4,5 V1 (anatomically located in Brodmann’s area 17 [BA 17]), includes the calcarine sulcus of the occipital cortex; this region carries information from the ipsilateral lateral geniculate nucleus to the ventral and dorsal streams.6 A voxel-based morphometry analysis has shown that the gray matter density in the occipital lobe (including V1) was decreased in DR patients.7 Patients with DR caused by type 1 DM (T1DM) showed lower local clustering in the occipital lobe compared with T1DM without retinopathy.8 Moreover, a diffusion-weighted magnetic resonance imaging (MRI) study revealed significantly greater apparent diffusion coefficient values of the visual center, supporting an association between DR and visual functional impairment.9 Notably, DR patients also exhibit functional alterations in the V1. Wang et al. found that patients with DR showed significant alterations in spontaneous brain activity in the bilateral occipital lobe.10 Based on a brain network analysis, Dai and colleagues demonstrated brain functional networks of V1 and secondary visual cortices were altered; these alterations might represent the severity of vision loss in patients with DR.11 The above-mentioned neuroimaging studies confirmed that specific functional and morphological alterations occurred in the V1 in patients with DR. However, few studies focused on whether patients with DR have characteristic changes in functional connectivity (FC) patterns in the V1.

The human brain is a complex dynamic system, which displays signals on blood oxygenation level-dependent (BOLD) imaging during rest. Low-frequency fluctuations (0~0.01 Hz) in BOLD signals at rest reflect spontaneous neuronal activity, as demonstrated by resting state functional MRI (rs-fMRI) studies.6,12 FC algorithm based on rs-fMRI data is a method of the spatiotemporal correlations and synchrony of the BOLD signals among anatomically different regions. The patterns of FC at rest show significant alterations in many diseases, such as Alzheimer’s disease,13 schizophrenia,14 early blindness,15 and glaucoma.12 Patients with primary open angle glaucoma had altered FC within and beyond the visual cortex since the early disease stage, and the findings showed clinical relevance of brain changes by correlating FC values with classical ophthalmologic measures. Another similar study using the FC method found patients with late blindness showed an altered FC between V1 and higher visual areas, and somatosensory cortices during rest. These studies all prove that FC, based on fMRI data, is a useful method for evaluating information transfer between brain regions. To date, no studies have investigated whether patients with DR show FC changes in the V1. Therefore, when compared with HCs and DM patients without retinopathy, we hypothesized that patients with DR would show characteristic alterations in FC patterns in the V1.

Based on this hypothesis, the aim of this rs-fMRI study was to explore differences in the FC strength of V1 between DR patients and control groups. These findings may provide meaningful information to aid in understanding the neuropathologic mechanism in patients with DR.

Patients and Methods

Subjects

In total, 35 DR patients (18 males and 17 females; 53.37±7.56 years), 22 type 2 DM (T2DM, 10 males and 12 females; 55.23±7.73 years) and 38 healthy controls (HCs) (16 males and 22 females; 54.46±8.06 years) were recruited for the present research in the Renmin Hospital of Wuhan University. The study program was carried out in accordance with the Declaration of Helsinki and accepted by the ethical committee of the Renmin Hospital of Wuhan University (No:WDRY2019-K032). These subjects were informed of the research procedures and signed the informed consent agreement prior to the research.

All subjects had to meet these criteria: (1) able to meet the physical demands of the imaging procedure; (2) no cardiovascular diseases; and (3) no serious brain disease (diagnosed by routine brain MRI [including T1WI, and T2WI/FLAIR] assessed by two experienced radiologists).

DR individuals met the following criteria: (1) all patients with DR caused by T2DM according to criteria recommended by the American Diabetes Association;16 (2) proliferative DR with retinal detachment and vitreous hemorrhage were excluded; (3) no additional eye diseases; (4) no ocular surgical history; (5) nonproliferative DR patients without diabetic nephropathy.

All HCs and patients with T2DM in this research had to meet these criteria: (1) normal eye examination and no retinal lesions (visual acuity ≥ 1.0); (2) HCs had no type of DM with normal hemoglobin A1c and blood glucose level; (3) women who were breastfeeding or pregnant during the research were excluded; (4) right-handed.

MRI Methods

The MRI scanning, performed on the 3.0 Tesla MR system (Discovery MR750W system; GE Healthcare, Milwaukee, WI, USA), contained both routine brain MRI and BOLD sequence. The fMRI images were acquired by applying a gradient echo-planar imaging sequence, which was documented in detail in a previous fMRI study.6 The parameters of the above sequences are shown in Table 1.

Table 1

Scanning Parameters of Structural T1-Weighted Images and BOLD Sequence

Sequence	TE(ms)	TR(ms)	Gap(mm)	Thickness(mm)	FOV(mm2)	FA	AM	VS(mm3)
T1WI	3.3	8.5	0	1.0	240×240	12°	256×256	–
BOLD	25	2000	1.2	3.0	240×240	90°	64×64	3.6×3.6×3.6

Abbreviations: T1WI, T1-weighted imaging; BOLD, blood oxygenation level dependent imaging; TR, repetition time; TE, echo time; AM, acquisition matrix; FOV, field of view; FA, flip angle; VS, voxel size.

Functional neuroimaging analysis software, Data Processing & Analysis of Brain Imaging (DPABI-Version 4.1, ), was applied to preprocess the fMRI images. After deleting the first ten images, the remaining 230 volumes for each participant need to be processed in the following steps, as referred to in a recent study:8 time calibration, corrected for head motion, spatial normalization and spatial smoothing.

The V1 area of brain (anatomically equivalent to BA 17) was selected as the region of interest (ROI) in the present study (Figure 1). We defined the center of V1 area as the seed point to calculate FC. The Montreal Neurological Institute coordinates of V1 area were (right V1: 7,-76, 10) and (left V1: −8,-76, 10). The radius of the ROI was set as 5 mm (shown as the red area in Figure 1). The Pearson correlation coefficient between the mean time course of ROI and all other regions’ time courses were calculated. In order to test the normality of the data, Fisher’s z transform analysis was carried out on the Pearson correlation coefficient to acquire an approximate normal distribution for the following statistical analysis.

Figure 1

The center of V1 was chosen as the seed point to calculate FC (the red area).

Ophthalmic Testing

All participants, assessed by two experienced ophthalmologists, received comprehensive eye examinations that included intraocular pressure measurement, best-corrected visual acuity (BCVA) test, slit lamp examination, and retinal examination by fundoscopy. The minimum criteria for the diagnosis of DR was the presence of at least one well-defined retinal microaneurysm in any field of the eye. All patients participating in the study had mild to severe non-proliferative DR and subjects in the DM group did not have retinopathy. Mild non-proliferative DR only shows microaneurysms, and moderate non-proliferative DR shows more than just microaneurysms but less than severe non-proliferative DR. Severe non-proliferative DR has any of the following signs: more than 20 intraretinal hemorrhages in each of 4 quadrants; definite venous beading in 2 quadrants; prominent intraretinal microvascular abnormalities in 1 quadrant and no signs of proliferative DR.

Statistical Analysis

Variables of clinical features of the three groups were analyzed by SPSS version 20.0 software (SPSS Inc, Chicago, Illinois, USA) using Chi-square test and independent-samples t-test. A p value less than 0.05 was considered statistically significant. The sample size was calculated by G-power version 3.1.9.2.

One-sample t-test was used to evaluate intragroup patterns of z-values functional connectivity (zFC) maps by applying SPM8 software (voxel-level: p<0.01). One-way ANOVA was performed to assess the significant zFC maps of V1 differences among three groups. Two-sample two-tailed t-test was conducted to investigate the group differences of the zFC maps and the pairs of zFC maps between V1 and vision-related regions applying the Gaussian random field (GRF) method, which was performed to correct for multiple comparisons and regressed covariates of sex and age applying the SPM 8 software (voxel-level: p<0.01, GRF correction, cluster-level: p<0.05).

Pearson correlation coefficient was conducted to evaluate the relationships between the zFC values of different brain regions and clinical features in DR applying SPSS version 20.0 software.

Results

Demographic and Clinical Characteristics

The demographic and clinical characteristics of DR group (18 males and 17 females, mean age: 53.37±7.56 years), T2DM group (10 males and 12 females; mean age: 55.23±7.73 years) and HC group (16 males and 22 females, mean age: 54.46±8.06 years) are summarized in Table 2. No statistically significant differences were found between the DR and other groups in education, intraocular pressure or age, but significant differences were found in BCVA of the left eye and right eye.

Table 2

Demographics and Visual Measurements Among Three Groups

	DR Group	DM Group	HC Group	DR Group/DM Group		DR Group/HC Group
				T-Values	p-values	T-Values	p-values
Gender (male/female)	18/17	10/12	16/22	N/A	N/A	N/A	N/A
Age (years)	53.37±7.56	55.23±7.73	54.46±8.06	−0.576	0.485	−0.545	0.532
Handedness	35 R	22 R	38 R	N/A	N/A	N/A	N/A
Education (years)	11.85±1.47	12.05±2.17	12.20±1.69	−0.354	0.896	−0.335	0.754
BMI (kg/m2)	23.56±2.43	23.28±2.76	23.02±1.98	1.258	0.187	1.548	0.124
Duration of diabetes (years)	5.52±3.96	5.45±2.65	N/A	1.048	0.654	N/A	N/A
Duration of DR (years)	1.34±0.76	N/A	N/A	N/A	N/A	N/A	N/A
BCVA-OD	0.48±0.28	1.08±0.25	1.16±0.15	−14.963	<0.001	−15.645	<0.001
BCVA-OS	0.43±0.31	1.14±0.18	1.24±0.22	−13.649	<0.001	−12.630	<0.001
Intraocular pressure	15.79±1.37	15.46±2.43	15.12±1.49	1.247	0.573	1.470	0.472
HbA1c (%)	7.25±1.54	7.38±1.68	5.64±0.96	−0.475	0.726	2.690	<0.001
Fasting blood glucose(mmol/L)	7.96±2.45	7.75±2.24	5.75±1.05	−0.568	0.564	3.254	<0.001

Notes: χ2 test for sex (n). Independent t-test for the other normally distributed continuous data (means ± SD).

Abbreviations: DR, diabetic retinopathy; DM, diabetes mellitus; HC, healthy control; N/A, not applicable; BCVA, best corrected visual acuity; OD, oculus dexter; OS, oculus sinister; Hb, glycosylated hemoglobin; BMI, body mass index.

Imaging of the Retinal Fundus

Figure 2 is a typical fundus photography of a right eye and a left eye, presenting the following signs of DR: retinal punctate hemorrhage, microaneurysms formation, and hard exudation. The minimum criteria for the diagnosis of DR was the presence of at least one well-defined retinal microaneurysm in any field of the eye. All patients enrolled in the study were non-proliferative DR, assessed by two experienced ophthalmologists.

Figure 2

The fundus photos with the typical presentations of DR.

FC Differences

Distribution mean zFC patterns of right V1 and left V1 in group-level of DR, DM and HC participants in the frequency band (0.01~0.08 Hz) are shown in Figure 3A and B. Intra-group patterns of zFC maps of left V1 and right V1 among the three groups are shown in Figures 4 and 5. The study found DR patients had a significantly lower FC between the bilateral calcarine (CAL)/left lingual gyrus (LING) (BA 17/18) and the left V1, and between the bilateral CAL/left LING (BA 17/18) and the right V1 compared with the HCs. Meanwhile, patients with DR exhibited higher FC strength between the left V1 and the bilateral Caudate/Olfactory/Orbital superior frontal gyrus (OSFG), and between the bilateral Caudate/Olfactory/OSFG (BA 3/4/6) and the right V1. Compared with DM group, patients with DR showed increased FC strength between the right CAL (BA 17/18) and the right V1. DM group exhibited lower FC strength between the left fusiform and the left V1, and between the bilateral CAL and the right V1 compared with HCs. Moreover, DM group was observed to have higher FC strength between the left superior frontal gyrus and the left V1 (voxel-level: p<0.01, GRF correction, cluster-level: p<0.05) (Tables 3 and 4).

Figure 3

Distribution group mean zFC among DR, DM and HC group.

Figure 4

Significant zFC maps of left V1 differences among three groups.

Figure 5

Significant zFC maps of right V1 differences among three groups.

Table 3

Significant Differences in the FC of Left V1 Among Three Groups

Brain Areas	MNI Coordinates			Voxels	BA	Peak T Values
	x	y	z
One-way ANOVA
Left OSFG	−24	2	−15	69	-	3.7166
DR vs HC
Bilateral Caudate/Olfactory/OSFG	12	21	−15	939	-	4.5818
Bilateral CAL/ left LING	−3	−72	9	433	17,18	−4.1422
DM vs HC
Left fusiform	−33	−45	−9	140	-	−4.0305
Left superior frontal gyrus	−9	30	36	143	-	3.6485

Notes: The statistical threshold was set at the voxel level with p<0.01 for multiple comparisons using Gaussian random-field theory (two-tailed, voxel-level: p<0.01, GRF correction, cluster-level: p<0.05).

Abbreviations: FC, function connectivity; V1, primary visual cortex; BA, Brodmann area; GRF, Gaussian random field; DR, diabetic retinopathy; DM, diabetes mellitus; HC, healthy control; LING, lingual gyrus; CAL, calcarine; OSFG, orbital superior frontal gyrus; MNI, Montreal Neurological Institute.

Table 4

Significant Differences in the FC of Right V1 Among Three Groups

Brain Areas	MNI Coordinates			Voxels	BA	Peak T Values
	x	y	z
One-way ANOVA
Right OSFG	12	21	−15	263	-	4.3383
Left OSFG	−24	39	−18	182	-	3.8161
DR vs DM
Right CAL	6	−75	12	248	17,18	7.2358
DM vs HC
Bilateral CAL	3	−75	9	1566	17,18	−9.0369
DR vs HC
Bilateral Caudate/Olfactory/OSFG	−24	12	−15	1065	3,4,6	8.798
Bilateral CAL/ left LING	−9	−81	−12	531	17,18	−3.9204

Notes: The statistical threshold was set at the voxel level with p<0.01 for multiple comparisons using Gaussian random-field theory (two-tailed, voxel-level: p<0.01, GRF correction, cluster-level: p<0.05).

Abbreviations: FC, function connectivity; V1, primary visual cortex; BA, Brodmann area; GRF, Gaussian random field; DR, diabetic retinopathy; DM, diabetes mellitus; HC, healthy control; LING, lingual gyrus; CAL, calcarine; OSFG, orbital superior frontal gyrus; MNI, Montreal Neurological Institute.

Correlation Analysis

We did not found obvious correlation between the mean FC values in different brain areas and the clinical characteristics in the DR patients (p>0.05).

Discussion

The FC method, based on rs-fMRI data, provides a novel tool to assess correlation coefficients of BOLD signal time courses among brain regions. To the best of our knowledge, this is the first study to investigate differences in the intrinsic FC patterns of V1 between patients with DR and HCs. We found that DR patients had a significant lower FC between the bilateral calcarine (CAL)/left lingual gyrus (LING) (BA 17/18) and the left V1, and between the bilateral CAL/left LING (BA 17/18) and the right V1 compared with the HCs. Meanwhile, patients with DR exhibited higher FC strength between the left V1 and the bilateral Caudate/Olfactory/Orbital superior frontal gyrus (OSFG), and between the bilateral Caudate/Olfactory/OSFG (BA 3/4/6) and the right V1. Compared with DM group, patients with DR showed increased FC strength between the right CAL (BA 17/18) and the right V1. DM group exhibited lower FC strength between the left fusiform and the left V1, and between the bilateral CAL and the right V1 compared with HCs. Moreover, DM group was observed to have higher FC strength between the left superior frontal gyrus and the left V1. These results may provide some useful information for understanding the neural mechanisms in DR.

The occipital lobe contains most of the anatomical areas of the visual cortex, including the three Brodmann areas (17, 18, and 19) related to vision.17 The CAL is the most reliable anatomic landmark of the medial part of the occipital lobe, which contains the V1 and visual area 2 (V2).18 Lesions of the CAL are reportedly associated with visual field defects and the perception of phosphenes.19 The LING (BA18), located in V2, receives strong feedforward connections from V1. Additionally, the LING is a crucial component of the dorsal visual pathway for visual processing and spatial memory.20 Visual signals received from the retina are transmitted in a stepwise manner along occipito-temporal and occipital-parietal pathways (i.e., ventral and dorsal visual streams). Furthermore, the CAL and LING are both crucial brain regions for the visual pathway.15 Patients with T2DM have a high risk of visual impairment; they have shown reduced FC density in the left CAL and left LING, compared with the HC group.21 Wang and colleagues reported that patients with DR showed abnormal brain activity in the occipital gyrus and LING, relative to HCs; these findings might be related to long-term reductions in stimulation of the visual center.10 Neuroimaging studies have demonstrated that dysfunction of the visual cortex (e.g., V1, CAL, and LING) is present in many patients with advanced visual deprivation, such as retinitis pigmentosa,15 high myopia,22 open globe injury,23 and retinal vein occlusion.24 From the perspective of pathophysiology, the difference between DR and DM is retinal microangiopathy. In our study, we observed reduced FC in the visual cortex between DR and HCs. DM group also showed reduced FC in the bilateral CAL when compared with HC group. These might reflect defects in visual processing in patients with DR. Thus, we speculate that long-term reduction in visual input in patients with DR may result in altered FC in the visual cortex. However, compared with DM group without retinopathy, patients with DR were observed to have an increased FC in the right CAL, which might be an enhanced compensatory connectivity induced by long-term reduction in visual input in patients with DR.

The caudate nucleus is a C-shaped structure in the basal ganglia striatum. On account of the extensive afferent and efferent connections, as well as cortico-striato-thalamic loops, the caudate nucleus is responsible for both motor and cognition function.25 There is an increasing number of reports that functional and structural abnormalities in the caudate nucleus are present in patients with hyperkinetic and hypokinetic movement disorders (e.g., Parkinson’s and Huntington’s diseases).26,27 The coordination of motor and visual cortices plays an important role in visuomotor function.28 Additionally, a recent analysis revealed that patients with DR had abnormal intranetwork FCs in the bilateral thalamus and right caudate nucleus of the basal ganglia network, suggesting that DR patients may exhibit movement control dysfunction.29 Dai et al. found that characteristic path length was enhanced and degree distribution was reduced (indicating greater connectivity) in the caudate nucleus and most frontal gyri of patients with DR, which may relate to compensatory recruitment.11 Consistent with the prior findings, our results showed that the DR group had enhanced FC between the bilateral V1 and bilateral caudate nucleus, which may reflect enhanced visuomotor function in patients with DR. We presume that the reorganization of the V1 and motor areas might contribute to the enhanced compensatory connectivity between the caudate and V1.

The olfactory region is present in the foremost part of the frontal lobe. This region performs an important role in olfactory function; thus, it exhibits high levels of functional and biological plasticity.30 There is increasing evidence that the olfactory region is closely linked to several key physiological functions, such as visual memory performance and olfactory modulation of vision.31,32 Olfactory dysfunction is also an early biomarker for neurodegenerative diseases, including Alzheimer’s disease.33 van Rijzingen et al. found that ablation of the olfactory bulb could impair visuospatial learning in rats.34 Some studies have investigated the olfactory modulation of visual perception; the results suggest that olfactory inputs can modulate visual information.32,35 Gouveri et al. reported that, in T2DM patients, the presence of olfactory dysfunction and retinopathy was significantly associated with lower olfactory scores.36 Additionally, evidence from a neuroimaging study suggests that patients with DR exhibit enhanced connectivity between the visual cortex and olfactory region.37 Accordingly, our results showed that patients with DR had enhanced FC between the bilateral olfactory region and bilateral V1, which might compensate for visual impairment in patients with DR. Therefore, we presume that reduction in visual input may promote the olfactory region to perform more robust visual modulation by enhanced connectivity with the V1.

Finally, DR individuals showed significantly enhanced FC between the bilateral V1 and bilateral OSFGs, relative to FC in the HC group. The OSFG, located in the prefrontal lobe, receives neural signals from the dorsolateral medial thalamus, temporal gyrus, olfactory system, and amygdala.4 Notably, the OSFG is responsible for regulation and modulation of mood and vision.8 The alteration of apparent diffusion coefficient values investigated applying the diffusion-weighted imaging in the brains of patients with DR was similar in both the visual area and orbitofrontal cortex.38 A previous fMRI study revealed enhancement of FC between the left orbitofrontal gyrus and the temporal gyrus in T1DM patients, which suggested a compensatory or adaptive response to impaired functional areas in those patients.8 Moreover, dysfunction in the OSFG has been observed in patients with DR.11,39 In our study, enhanced FC between the bilateral V1 and OSFG indicated improved functional reorganization, which may relate to OSFG compensatory regulation or recruitment of visual resources from other regions.

The main limitation of this research was its relatively small sample size. Moreover, the study did not include assessments of retinal neurodegeneration and cognitive impairment in DR patients. We will expand the sample size and collect more complete clinical data in our future study.

Conclusion

In conclusion, our results indicated that patients with DR exhibited FC disruptions between V1 and higher visual regions at rest, which may reflect the aberrant information communication in the brain of DR individuals, and imply the presence of impaired vision in these patients. The findings may be helpful for understanding the neuromechanism of vision disorder in patients with DR. Based on our findings, systemic use of drugs to improve microcirculation and neurotrophic drugs might be beneficial to patients with DR.