Machine learning analysis reveals abnormal functional network hubs in the primary angle-closure glaucoma patients

Abstract

Background: Previous neuroimaging studies demonstrated that primary angle-closure glaucoma (PACG) patients were accompanied by vision and vision-related brain regions changes. However, whether the whole-brain functional network hub changes occur in PACG patients remains unknown. Purpose: The purpose of the study was to investigate the brain function network hubs changes in PACG patients using the voxel-wise degree centrality(DC) method. Materials and methods: Thirty-one PACG patients (21 male and 10 female) and 31 healthy controls (21 male and 10 female) closely matched in age, sex, and education were enrolled in the study. DC method was applied to investigate the brain function network hubs changes in PACG patients. Moreover, the SVM method was applied to distinguishing PACG patient and HC patients. Results: Compared with HC, PACG patients had significantly higher DC values in the right fusiform, left middle temporal gyrus and left Cerebelum_4_5. Meanwhile, PACG patients had significantly lower DC values in the right calcarine, right postcentral gyrus, left precuneus gyrus and left postcentral gyrus. Furthermore, The SVM classification reaches a total accuracy of 72.58% and the ROC curve of the SVM classifier with an AUC value of 0.85.(r=0.25) Conclusion: Our results highlight that .PACG patients showed widespread brain functional network hubs dysfunction relative to visual network, auditory network, default mode network and cerebellum network, which might be shed light on the neural mechanism of optic atrophy in PACG patients.

Introduction

Primary angle-closure glaucoma (PACG) is a serious and irreversible blinding eye disease worldwide. Angle closure is a key pathological process in PACG patients. The main clinical manifestations of PACG are elevated intraocular pressure, eye pain, and vision loss. There are multiple risk factors for PACG, such as genetic predisposition (Wiggs and Pasquale, 2017), hypertension (Lowry and Sanders, 2021), diabetes mellitus (Song et al., 2016), and autonomic dysfunction (Na et al., 2010). Elevated intraocular pressure leads to retinal ganglion cell apoptosis in PACG patients, which is an important pathological feature of PACG. Visual field examination is important for the assessment of visual function in glaucoma patients (De Moraes et al., 2017). However, recent studies have shown that glaucoma causes optic nerve atrophy as well as trans-synaptic degeneration of visual pathways (Lawlor et al., 2018).

Magnetic resonance imaging (MRI) techniques have been successfully applied to identify changes in neural function and structure within the visual cortex of glaucoma patients. Zhou et al. (2017) demonstrated an expanded representation of parafoveal areas in the visual cortex of primary open-angle glaucoma (POAG) patients, compared with healthy controls (HCs). Duncan et al. (2012) demonstrated that POAG patients had differences in cerebral blood flow (CBF) between ventral V1 and dorsal V1. These differences were correlated with visual function in the superior and inferior visual fields. Wang et al. (2017) reported that the POAG patients had decreased functional connectivity within the visual network. Moreover, glaucoma patients exhibited structural changes in the brain. Using diffusion kurtosis imaging, Li et al. (2020) demonstrated that patients with normal-tension glaucoma had microstructural abnormalities in bilateral BA17, BA18, and BA19. Hernowo et al. (2011) found reduced visual pathway volumes in glaucoma patients; affected areas included the optic nerves, optic chiasm, optic tract, and optic radiations. The above studies indicated that patients with glaucoma have functional and structural changes in the visual pathways and visual cortex. To our knowledge, there have been few studies of whole-brain network changes in PACG patients.

The human brain is a complex functional connectome (Sporns, 2011) that plays important roles in physiological processes such as visual function (Wang et al., 2018) and auditory function (Gurtubay-Antolin et al., 2021). Previous neuroimaging studies demonstrated that brain functional network hubs are involved in the integration of information from network elements. Hubs enable the integration of diverse sources of information; they also balance opposing pressures during the evolution of segregated networks. Increased degree centrality (DC) values indicate the enhanced local network information transmission. Meanwhile, the decreased DC values indicate the impaired local network information transmission. The voxel-wise DC method is a graph-based technique that measures the functional relationships of specific voxels throughout the connectivity matrix of the brain, rather than in specific nodes or networks (Zuo et al., 2012). The advantage of the DC method is that the functional network of the whole brain voxel can be calculated without presupposition experimental hypothesis. Machine learning methods combined with the fMRI method have provided a systematic approach for developing an automatic, objective, and sophisticated classification of diseases. The support vector machine (SVM) method is the most commonly used supervised machine learning algorithm for MRI classification that enables individual-level classification and detects biomarkers on the basis of neuroimaging data. Glaucoma is not only associated with functional abnormalities in the visual area but has also been associated with neurodegenerative diseases. However, the effects of optic atrophy on whole-brain network hubs in PACG patients remain unclear. Thus, we hypothesized that PACG patients are accompanied by brain functional network hub changes.

Based on this hypothesis, the present study was performed to determine whether PACG patients exhibited abnormal changes in brain function network hubs. Moreover, the SVM method was used to assess classification power using the DC map as a feature.

Materials and methods

Participants

In total, thirty-one patients with PACG (21 males and 10 females) were enrolled from the Department of Ophthalmology, Jiangxi Provincial People’s Hospital. The diagnostic criteria of PACG were: (1) the intraocular pressure was greater than 21 mmHg in both eyes; (2) optic disc/cup area > 0.6; (3) typical vision field defect (paracentric obscura, arcuate obscura, nasal ladder, fan-shaped field defect, and peripheral field defect); (4) without any other ocular diseases.

The exclusion criteria of PACG were: (1) advanced PACG patients are associated with severe eye pain; (2) PACG patients with a history of surgery; (3) PACG patients with glaucoma-related eye complications, including neovascular glaucoma, high myopia, optic neuritis, cataracts, eye atrophy, and corneal edema; and (4) PACG patients with psychiatric disorders, brain trauma, and other diseases.

Thirty-one HCs (21 males and 10 females) were also recruited for this study. The inclusion criteria were: (1) without any ocular disease with uncorrected visual acuity (VA) > 1.0; 2) no cardiovascular system diseases; and no psychiatric disorders.

Ethical statement

All research methods followed the Declaration of Helsinki and were approved by the Ethical Committee for Medicine of Jiangxi Provincial People’s Hospital. Participants enrolled in the study of their own accord and were informed of the purpose, methods, as well as potential risks before signing an informed consent form.

Clinical evaluation

The VA of all subjects was measured by using the logMAR table and Intraocular pressure was assessed by using automatic intraocular pressure measurement. The best-corrected VA of both eyes was measured in each group.

Magnetic resonance imaging data acquisition

MRI scanning was performed on a 3-Tesla MR scanner (750W GE Healthcare, Milwaukee, WI, United States) with an eight-channel head coil. All participants were required to close their eyes without falling asleep when undergoing MRI scanning. The subjects should be kept calm and not engaged in specific thoughts.

Data pre-processing

All preprocessing was performed using the toolbox for Data Processing & Analysis of Brain Imaging (DPABI)1 (Yan et al., 2016) and briefly following the steps: (1) the first ten volumes of BOLD maps were removed. (2) Slice timing and head motion correction were conducted. (4) Normalized fMRI data were re-sliced with a resolution of 3 × 3 × 3 mm³. (5) Data detrends; (6) Linear regression analysis was applied to regress out several covariates (mean framewise displacement, global brain signal, and averaged signal from cerebrospinal fluid and white matter). (7) Temporal band-pass filtering was performed (0.01–0.08 Hz).

Voxel-wise degree centrality analysis

According to previous studies, we analyzed the different degree of centrality values with a correlation threshold of (r = 0.15, 0.2, 0.25, 0.3) in this study (Zuo et al., 2012). The voxel-wise functional network was generated for each subject, for which we took each voxel as a node and inter-voxel correlations as the edge. Within the default brain mask provided by DPABI (in the MNI-152 standard space with 3 × 3 × 3 mm³ voxel size and resolution of 61 × 73 × 61), we used the preprocessed functional images to perform a voxel-wise correlation analysis. From the n × n Pearson’s correlation coefficient matrix, a map of the degree of the connectivity was computed by counting for each voxel the number of voxels it was correlated to above a threshold of r > 0.25. A high threshold was chosen to eliminate counting voxels that had low temporal correlation attributable to signal noise. The z-score map was smoothed with a 6-mm full-width-half-maximum Gaussian kernel.

Support vector machine analysis

The SVM algorithm was performed by using the Pattern Recognition for Neuroimaging Toolbox (PRoNTo) software Cyclotron Research Centre, University of Liège, Belgium (Schrouff et al., 2013). The following steps were followed: (1) the DC maps (r = 0.15, 0.2, 0.25, 0.3) were used as classification feature. (2) Then, the leave-one-out cross-validation (LOOCV) technique was applied to classifier validation. (3) The total accuracy, specificity, sensitivity, and area under the receiver operating characteristic curve (AUC) were calculated.

Statistical analysis

The chi-square test was used to calculate the sex and handness and independent sample t-test was used for age, education, and BCVA between the two groups.

The one-sample t-test was conducted to assess the group mean of DC maps (r = 0.15, 0.2, 0.25, 0.3). The two-sample t-test was used to compare the two group differences in the DC maps (r = 0.15, 0.2, 0.25, 0.3) using the Gaussian random field (GRF) method. (two-tailed, voxel-level P < 0.01, GRF correction, cluster-level P < 0.05).

Results

Demographics and disease characteristics

There were no statistically significant differences between the PACG and HC groups in gender, education, or age, but significant differences in BCVA of right eye (p < 0.001), left eye (p < 0.001). The results of these data are listed in Table 1.

TABLE 1

Table 1. Demographic and clinical measurements between patients with PACG and HCs.

Different degree centrality between primary angle-closure glaucoma and healthy control group

The group means of DC maps of the PACG and HC [Figure 1A (r = 0.15) Figure 1B (r = 0.20), Figure 1C (r = 0.25) and Figure 1D (r = 0.3)]. Compared with HC, PACG patients had significantly higher DC values in the left cerebelum_Crus2, left cerebelum_4_5, right para hippocampal, and left thalamus. Meanwhile, PACG patients had significantly lower DC values in the right postcentral, left precuneus, and right postcentral (r = 015) (Table 2 and Figure 2A).

FIGURE 1

Figure 1. The spatial of the DC within the PACG and HC with different correlation thresholds (r = 0.15, 0.2, 0.25, 0.3). The group means DC maps of the PACG and HC [(A) (r = 0.15), (B) (r = 0.20), (C) (r = 0.25), and (D) (r = 0.30)]. PACG patients showed remarkably similar altered degree of centrality brain areas relative to healthy control in the different correlation thresholds (r = 0.15, 0.2, 0.25, 0.3) (FDR correction p < 0.001). DC, degree centrality; PACG, primary angle-closure glaucoma; HC, health control; L, left hemisphere; R, right hemisphere.

TABLE 2

Table 2. Significant differences in the DC between the two groups (r = 0.15).

FIGURE 2

Figure 2. Comparison of the DC between the PACG and HC. (A) (r = 0.15), (B) (r = 0.20), (C) (r = 0.25), (D) (r = 0.30). Significant degree centrality differences were found between the two groups. The blue areas indicate lower degree of centrality values (two-tailed, voxel-level P < 0.01, GRF correction, cluster-level P < 0.05). GRF, Gaussian random field; L, left hemisphere; MOG, middle occipital gyrus; and R, right hemisphere.

Compared with HC, PACG patients had significantly higher DC values in the right fusiform, left middle temporal gyrus, and left Cerebelum_4_5. Meanwhile, PACG patients had significantly lower DC values in the right calcarine, right postcentral gyrus, left precuneus gyrus, and left postcentral gyrus (r = 0.20) (Table 3 and Figure 2B).

TABLE 3

Table 3. Significant differences in the DC between the two groups (r = 0.20).

Compared with HC, PACG patients had significantly higher DC values in the right fusiform, left middle temporal gyrus, and left Cerebelum_4_5. Meanwhile, PACG patients had significantly lower DC values in the right calcarine, right postcentral gyrus, left precuneus gyrus, and left postcentral gyrus (r = 0.25) (Table 4 and Figure 2C).

TABLE 4

Table 4. Significant differences in the DC between the two groups (r = 0.25).

Compared with HC, PACG patients had significantly higher DC values in the right fusiform gyrus, left middle temporal gyrus, left cerebelum_4_5, and right cerebelum_4_5. Meanwhile, PACG patients had significantly lower DC values in the right calcarine, right postcentral, left precuneus, left postcentral and left paracentral_lobule (r = 0.30) (Table 5 and Figure 2D).

TABLE 5

Table 5. Significant differences in the DC between the two groups (r = 0.30).

Support vector machine results

The SVM classification reaches a total accuracy of 72.58% and the ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.15); (Figure 3A). The SVM classification reaches a total accuracy of 75.81% and the ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.20); (Figure 3B). The SVM classification reaches a total accuracy of 74.19% and the ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.25); (Figure 3C). The SVM classification reaches a total accuracy of 77.42% and the ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.30); (Figure 3D).

FIGURE 3

Figure 3. Classification results using machine learning analysis based on DC values. Function values of two groups (class 1: PACG group; class 2: HC group); The ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.15). (A) Function values of two groups (class 1: PACG group; class 2: HC group); the ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.20). (B) Function values of two groups (class 1: PACG group; class 2: HC group); The ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.25). (C) Function values of two groups (class 1: PACG group; class 2: HC group); The ROC curve of the SVM classifier with an AUC value of 0.85 (r = 0.30) (D).

Discussion

In this study, compared with HCs, PACG patients had significantly higher DC values in the right fusiform gyrus, left middle temporal gyrus, and left cerebellum_4_5. Conversely, PACG patients had significantly lower DC values in the right calcarine sulcus, right postcentral gyrus, left precuneus, and left postcentral gyrus. The SVM classifier had an overall accuracy of 72.58%; the ROC curve of the SVM classifier had an AUC value of 0.85 (r = 0.25). These results demonstrated that PACG patients had widespread brain functional network hub dysfunction that affected the visual, auditory, default mode, and cerebellar networks. The findings might provide insights into the neural mechanism that underlies optic atrophy in PACG patients.

Notably, we found that PACG patients had significantly lower DC values in the right calcarine sulcus. The calcarine sulcus is a core component of the visual cortex; it plays an important role in transmitting visual information to the higher visual cortex. Previous neuroimaging studies showed visual dysfunction and vision-related cortex dysfunction in PACG patients. Colbert et al. (2021). found that glaucomatous mice with elevated intraocular pressure had decreased fractional anisotropy and increased radial diffusivity along the optic nerves and optic tract. Fujishiro et al. (2022) reported that ocular hypertension could lead to reduced cytochrome oxidase expression in the visual cortex. Pankowska et al. (2022) found that normal-tension glaucoma patients had reduced thickness in the right lateral occipital gyrus and left lingual gyrus. Consistent with these findings, we observed significantly lower DC values in the right calcarine sulcus in PACG patients, which suggested abnormal functional connectivity within the visual network.

Another important finding in this study was that PACG patients had significantly lower DC values in the right postcentral gyrus and left postcentral gyrus. The postcentral gyrus has an important role in sensorimotor function. There are two possible reasons for the phenomenon observed in our study. First, PACG patients exhibit visual loss associated with optic atrophy. Using diffusion tensor MRI, Yang et al. (2018) demonstrated that glaucomatous mice had extensive deterioration of visuomotor function by 9 months of age. Zwierko et al. (2022) reported differences in visuomotor task performance among older adults with moderate and advanced glaucoma. Dive et al. (2016) demonstrated visuomotor behavior in patients with glaucoma-related peripheral field reduction during the completion of natural movements. Thus, decreased DC values in the postcentral gyrus might indicate visuomotor dysfunction in PACG patients. Second, the postcentral gyrus plays an important role in nociceptive processing. Glaucoma patients experience painful swelling in both eyes. Dong et al. (2019) demonstrated that patients with eye pain had decreased voxel-mirrored homotopic connectivity in the precentral/postcentral gyrus. Additionally, Pan et al. (2018) reported that eye pain patients had a lower amplitude of low-frequency fluctuation in the left and right precentral/postcentral gyrus and left precuneus. Consistent with these findings, our study revealed that PACG patients had lower DC values in the postcentral gyrus, which might reflect visuomotor dysfunction and eye pain in PACG patients.

Finally, we found that PACG patients had significantly lower DC values in the left precuneus and increased DC values in the left middle temporal gyrus. The precuneus and middle temporal gyrus are core components of the default mode network (DMN). The DMN is an endogenous neural network that shows consistently higher blood oxygenation level-dependent activity during rest. It plays an important role in self-referential thought and introspection; these actions involve various higher cognitive functions such as memory, prospection, and self-processing (Raichle et al., 2001). Wang et al. (2017) demonstrated that POAG patients had decreased functional connectivity in the DMN. Giorgio et al. (2020) reported that ocular hypertension patients had decreased functional connectivity in key cognitive networks [DMN and frontoparietal working memory network (WMN)]. Our finding of significantly lower DC values in the DMN of PACG patients suggests that those patients are experiencing cognitive decline.

In our study, we used the DC maps as a feature. The SVM model was applied to investigate the SVM classification reaches a total accuracy of 72.58%–77.42% and the ROC curve of the SVM classifier with an AUC value of 0.85. The SVM modal showed a high sensitivity for distinguishing the two groups. The DC maps could be sensitive biomarkers for distinguishing patients with PACG from HCs.

Some limitations should be mentioned in the study. First, the sample size of the study is relatively small. Second, the DC values are based on blood oxygen levels dependent (BOLD) on signals, which might be affected by physiological noise. Third, no multimodal MRI methods were used to verify the results of this study.

Conclusion

Our results showed that PACG patients showed widespread brain functional network hub dysfunction relative to the visual network, auditory network, default mode network, and cerebellum network, which might shed new light on the neural mechanism of optic atrophy in PACG patients. Thus, the DC maps could be sensitive biomarkers for distinguishing patients with PACG from HCs. It provides an important imaging reference for clinicians in early diagnosis.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Medical Ethics Committee of the Jiangxi Provincial People’s Hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

R-BC, Y-LZ, HL, and XH contributed to data collection and statistical analyses, wrote the manuscript, designed the protocol, contributed to the MRI analysis, designed the study, oversaw all the clinical aspects of study conduct, and prepared the manuscript. All authors contributed to the article and approved the submitted version.

Acknowledgments

We acknowledge the assistance provided by the Natural Science Foundation of Jiangxi Province (20212BAB216058) and the Jiangxi Provincial Health Technology Project (202210012).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1. ^ http://www.rfmri.org/dpabi

References

Colbert, M. K., Ho, L. C., van der Merwe, Y., Yang, X., McLellan, G. J., Hurley, S. A., et al. (2021). Diffusion tensor imaging of visual pathway abnormalities in five glaucoma animal models. Invest. Ophthalmol. Vis. Sci. 62:21. doi: 10.1167/iovs.62.10.21

De Moraes, C. G., Liebmann, J. M., and Levin, L. A. (2017). Detection and measurement of clinically meaningful visual field progression in clinical trials for glaucoma. Prog. Retin. Eye Res. 56, 107–147. doi: 10.1016/j.preteyeres.2016.10.001

Dive, S., Rouland, J. F., Lenoble, Q., Szaffarczyk, S., McKendrick, A. M., and Boucart, M. (2016). Impact of peripheral field loss on the execution of natural actions: a study with glaucomatous patients and normally sighted people. J. Glaucoma 25, e889–e896. doi: 10.1097/IJG.0000000000000402

Dong, Z. Z., Zhu, F. Y., Shi, W. Q., Shu, Y. Q., Chen, L. L., Yuan, Q., et al. (2019). Abnormalities of interhemispheric functional connectivity in individuals with acute eye pain: a resting-state fMRI study. Int. J. Ophthalmol. 12, 634–639. doi: 10.18240/ijo.2019.04.18

Duncan, R. O., Sample, P. A., Bowd, C., Weinreb, R. N., and Zangwill, L. M. (2012). Arterial spin labeling fMRI measurements of decreased blood flow in primary visual cortex correlates with decreased visual function in human glaucoma. Vision Res. 60, 51–60. doi: 10.1016/j.visres.2012.03.012

Fujishiro, T., Honjo, M., Kawasaki, H., and Aihara, M. (2022). Visual cortex damage in a ferret model of ocular hypertension. Jpn. J. Ophthalmol. 66, 205–212. doi: 10.1007/s10384-022-00901-8

Giorgio, A., Zhang, J., Costantino, F., De Stefano, N., and Frezzotti, P. (2020). Altered large-scale brain functional connectivity in ocular hypertension. Front. Neurosci. 14:146. doi: 10.3389/fnins.2020.00146

Gurtubay-Antolin, A., Battal, C., Maffei, C., Rezk, M., Mattioni, S., Jovicich, J., et al. (2021). Direct structural connections between auditory and visual motion-selective regions in humans. J. Neurosci. 41, 2393–2405. doi: 10.1523/JNEUROSCI.1552-20.2021

Hernowo, A. T., Boucard, C. C., Jansonius, N. M., Hooymans, J. M., and Cornelissen, F. W. (2011). Automated morphometry of the visual pathway in primary open-angle glaucoma. Invest. Ophthalmol. Vis. Sci. 52, 2758–2766. doi: 10.1167/iovs.10-5682

Lawlor, M., Danesh-Meyer, H., Levin, L. A., Davagnanam, I., Vita, E. De, and Plant, G. T. (2018). Glaucoma and the brain: trans-synaptic degeneration, structural change, and implications for neuroprotection. Surv. Ophthalmol. 63, 296–306. doi: 10.1016/j.survophthal.2017.09.010

Li, T., Qu, X., Chen, W., Wang, Q., Wang, H., Wang, Y., et al. (2020). Altered information flow and microstructure abnormalities of visual cortex in normal-tension glaucoma: evidence from resting-state fMRI and DKI. Brain Res. 1741:146874. doi: 10.1016/j.brainres.2020.146874

Lowry, E. A., and Sanders, D. S. (2021). Hypertension management and glaucoma: hypothesizing causes in correlational data. Ophthalmology 128, 401–402. doi: 10.1016/j.ophtha.2020.10.020

Na, K. S., Lee, N. Y., Park, S. H., and Park, C. K. (2010). Autonomic dysfunction in normal tension glaucoma: the short-term heart rate variability analysis. J. Glaucoma 19, 377–381. doi: 10.1097/IJG.0b013e3181c4ae58

Pan, Z. M., Li, H. J., Bao, J., Jiang, N., Yuan, Q., Freeberg, S., et al. (2018). Altered intrinsic brain activities in patients with acute eye pain using amplitude of low-frequency fluctuation: a resting-state fMRI study. Neuropsychiatr. Dis. Treat. 14, 251–257. doi: 10.2147/NDT.S150051

Pankowska, A., Matwiejczuk, S., Koziol, P., Zarnowski, T., Pietura, R., and Kosior-Jarecka, E. (2022). Visual tract degradation in bilateral normal-tension glaucoma-cortical thickness maps and volumetric study of visual pathway areas. J. Clin. Med. 11:1907. doi: 10.3390/jcm11071907

Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682. doi: 10.1073/pnas.98.2.676

Schrouff, J., Rosa, M. J., Rondina, J. M., Marquand, A. F., Chu, C., Ashburner, J., et al. (2013). PRoNTo: pattern recognition for neuroimaging toolbox. Neuroinformatics 11, 319–337. doi: 10.1007/s12021-013-9178-1

Song, B. J., Aiello, L. P., and Pasquale, L. R. (2016). Presence and risk factors for glaucoma in patients with diabetes. Curr. Diab. Rep. 16:124. doi: 10.1007/s11892-016-0815-6

Sporns, O. (2011). The human connectome: a complex network. Ann. N.Y. Acad. Sci. 1224, 109–125. doi: 10.1111/j.1749-6632.2010.05888.x

Wang, J., Li, T., Zhou, P., Wang, N., Xian, J., and He, H. (2017). Altered functional connectivity within and between the default model network and the visual network in primary open-angle glaucoma: a resting-state fMRI study. Brain Imaging Behav. 11, 1154–1163. doi: 10.1007/s11682-016-9597-3

Wang, Z., Zeljic, K., Jiang, Q., Gu, Y., Wang, W., and Wang, Z. (2018). Dynamic network communication in the human functional connectome predicts perceptual variability in visual illusion. Cereb. Cortex 28, 48–62. doi: 10.1093/cercor/bhw347

Wiggs, J. L., and Pasquale, L. R. (2017). Genetics of glaucoma. Hum. Mol. Genet. 26, R21–R27. doi: 10.1093/hmg/ddx184

Yan, C. G., Wang, X. D., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4

Yang, X. L., van der Merwe, Y., Sims, J., Parra, C., Ho, L. C., Schuman, J. S., et al. (2018). Age-related changes in eye, brain and visuomotor behavior in the dba/2j mouse model of chronic glaucoma. Sci. Rep. 8:4643. doi: 10.1038/s41598-018-22850-4

Zhou, W., Muir, E. R., Nagi, K. S., Chalfin, S., Rodriguez, P., and Duong, T. Q. (2017). Retinotopic fmri reveals visual dysfunction and functional reorganization in the visual cortex of mild to moderate glaucoma patients. J. Glaucoma 26, 430–437. doi: 10.1097/IJG.0000000000000641

Zuo, X. N., Ehmke, R., Mennes, M., Imperati, D., Castellanos, F. X., Sporns, O., et al. (2012). Network centrality in the human functional connectome. Cereb. Cortex 22, 1862–1875. doi: 10.1093/cercor/bhr269

Zwierko, T., Jedziniak, W., Florkiewicz, B., Lesiakowski, P., Sliwiak, M., Kirkiewicz, M., et al. (2022). Physical activity is associated with improved visuomotor processing in older adults with moderate and advanced glaucomatous visual field defect: a cross-sectional study. Int. J. Environ. Res. Public Health 19:1760. doi: 10.3390/ijerph19031760