==== Front
Commun Biol
Commun Biol
Communications Biology
2399-3642
Nature Publishing Group UK London

37386293
4918
10.1038/s42003-023-04918-8
Article
GABA decrease is associated with degraded neural specificity in the visual cortex of glaucoma patients
http://orcid.org/0000-0002-9739-3474
Bang Ji Won JiWon.Bang@nyulangone.org

1
Parra Carlos 1
Yu Kevin 1
Wollstein Gadi 123
http://orcid.org/0000-0002-8885-3766
Schuman Joel S. 1234
http://orcid.org/0000-0003-4012-7084
Chan Kevin C. chuenwing.chan@fulbrightmail.org

12345
1 grid.137628.9 0000 0004 1936 8753 Department of Ophthalmology, NYU Grossman School of Medicine, NYU Langone Health, New York University, New York, New York, 10017 USA
2 grid.137628.9 0000 0004 1936 8753 Center for Neural Science, College of Arts and Science, New York University, New York, New York, 10003 USA
3 grid.137628.9 0000 0004 1936 8753 Department of Biomedical Engineering, Tandon School of Engineering, New York University, New York, New York, 11201 USA
4 grid.137628.9 0000 0004 1936 8753 Neuroscience Institute, NYU Grossman School of Medicine, NYU Langone Health, New York University, New York, New York, 10016 USA
5 grid.137628.9 0000 0004 1936 8753 Department of Radiology, NYU Grossman School of Medicine, NYU Langone Health, New York University, New York, New York, 10016 USA
29 6 2023
29 6 2023
2023
6 67923 11 2022
5 5 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Glaucoma is an age-related neurodegenerative disease of the visual system, affecting both the eye and the brain. Yet its underlying metabolic mechanisms and neurobehavioral relevance remain largely unclear. Here, using proton magnetic resonance spectroscopy and functional magnetic resonance imaging, we investigated the GABAergic and glutamatergic systems in the visual cortex of glaucoma patients, as well as neural specificity, which is shaped by GABA and glutamate signals and underlies efficient sensory and cognitive functions. Our study shows that among the older adults, both GABA and glutamate levels decrease with increasing glaucoma severity regardless of age. Further, our study shows that the reduction of GABA but not glutamate predicts the neural specificity. This association is independent of the impairments on the retina structure, age, and the gray matter volume of the visual cortex. Our results suggest that glaucoma-specific decline of GABA undermines neural specificity in the visual cortex and that targeting GABA could improve the neural specificity in glaucoma.

Reduced GABA in the visual cortex of glaucoma patients is shown to be tightly linked to disruption of stimulus-specific neural representations.

Subject terms

Striate cortex
Sensory processing
https://doi.org/10.13039/100006312 BrightFocus Foundation (BrightFocus) G2019103 Bang Ji Won issue-copyright-statement© Springer Nature Limited 2023
==== Body
pmcIntroduction

Glaucoma is an age-related neurodegenerative disease characterized by the gradual loss of retinal ganglion cells. Such damage disrupts the transmission of the retinal signals to the brain and eventually leads to irreversible blindness. With the increasing aging population, it is estimated that 111.8 million people worldwide will be affected by glaucoma in 2040, suggesting that glaucoma is a major public health problem1. Nevertheless, the exact mechanisms underlying glaucoma remain unclear. Current clinical treatments focus on intraocular pressure (IOP) reduction with laser, medication, or surgery. However, IOP is a major risk factor but not the cause of the disease. Glaucoma progression in certain patients cannot be attenuated even with controlled IOP, indicating that glaucoma cannot be explained by IOP elevation alone. To reduce the prevalence of this irreversible but preventable disease, we need a better understanding of the mechanisms of glaucoma that will lead to improved treatment strategies beyond IOP control.

An increasing amount of studies has pointed out that glaucomatous degeneration occurs along the visual pathway including the optic nerves2,3, optic chiasm3, lateral geniculate nucleus3–9, optic radiation2,10, and visual cortex3,6,7,9,11–13. Further, recent studies have suggested that glaucoma may share some pathogenic mechanisms with Alzheimer’s disease, which is characterized by the accumulation of amyloid β and tau14,15. Critically, glaucoma patients were shown to present abnormal tau proteins in the retina16 and vitreous fluid17. Likewise, in animal models of glaucoma, amyloid β and tau were found in the retina18, lateral geniculate nucleus, and even in the primary visual cortex19, although the amount of deposits was relatively lower in the visual cortex than in the anterior visual pathway. It suggests that the disease may spread trans-synaptically from the anterior to posterior visual pathways19.

Accumulation of amyloid β and tau in Alzheimer’s disease was found to impair the glutamatergic and GABAergic systems, which are the main excitatory and inhibitory systems in the brain. Exposure to amyloid β can increase extracellular glutamate levels by preventing glutamate uptake20,21 or potentiating glutamate release22,23. Elevated amyloid β and tau can also impair synaptic inhibition via downregulation of GABAA receptor24 and induce substantial loss of GABAergic neurons25–27 and GABA-synthetic enzyme25. Of note, loss of GABAergic interneurons co-localizes with tau markers25. If glaucoma patients present amyloid β and tau accumulation in the visual system of the brain, the glutamatergic and/or GABAergic signals are likely to be impaired similar to individuals with Alzheimer’s disease.

In particular, the GABAergic signals play a critical role in shaping the response patterns of neuronal populations28. Recent magnetic resonance spectroscopy (MRS) studies revealed a tight relationship between GABAergic signals and neural specificity at the level of neuronal population. As the cortical GABA levels decrease, neural activity patterns associated with different categories become similar, thus more confusable29,30. Increased confusability of neural activity patterns then could potentially degrade behavioral performance31. Indeed, loss of neural specificity was observed to be associated with declines in memory performance32 and fluid processing abilities such as executive function and speed of processing33. Relatedly, glaucoma patients were found to present impairments on some of the perceptual and cognitive functions where GABAergic signals are thought to be involved. These include abnormalities in visual crowding effect34, binocular rivalry35–38, visual categorization39,40, visual attention41, and cognitive capacity42. Further, glaucoma is thought to accelerate the speed of aging process43, which involves gradual reduction of cortical GABA levels44–49. It suggests that glaucomatous degeneration may involve the reduction of GABA to a greater extent than healthy aging and that this reduction of GABA may affect the neural specificity. Recent hypothesis-independent pathway analysis also implicates GABA and acetyl-CoA metabolism as the novel pathway that is associated with primary open-angle glaucoma50. Nevertheless, other than the limited studies on GABAergic involvements in the retina of experimental glaucoma models51,52, no studies have directly examined glaucoma-specific changes in the GABAergic signals and their relationship with neural specificity in the visual cortex of glaucoma patients.

Therefore, in the current study, we investigated whether the GABA and glutamate levels in the visual cortex are affected by glaucoma and whether these neurochemical changes are associated with the neural specificity independent of the impairments of the retina structure and age. To address these questions, we recruited glaucoma patients and age-matched healthy subjects and conducted functional magnetic resonance imaging (fMRI), MRS of the brain, as well as clinical ophthalmic tests including Humphrey visual field perimetry and optical coherence tomography (OCT) of the retina and optic nerve. Using principal component analyses (PCA) and regression modeling, we observed that among our older adult subjects, GABAergic and glutamatergic signals in the visual cortex gradually reduced with increasing severity of glaucoma regardless of age. Further, we observed that the neural specificity in the visual cortex was tightly related to the glaucoma-specific reduction of GABAergic signals, but not that of glutamatergic signals. This association between GABA and neural specificity was observed independent of age or impairments on the retina. These findings indicate the importance of the GABAergic action in the visual cortex and its implications in sensory encoding in glaucoma.

Results

Forty glaucoma patients and twenty-four age-matched healthy subjects underwent clinical ophthalmic exams and magnetic resonance imaging. Their demographic and clinical characteristics can be found in Table 1. Specifically, we obtained each individual’s structural and functional images of the whole brain, neurochemical profiles of the visual cortex (Fig. 1) as well as clinical ophthalmic measures including peripapillary retinal nerve fiber layer (pRNFL) thickness, macular ganglion cell-inner plexiform layer (mGCIPL) thickness, optic nerve head cup-to-disc (C/D) ratio, and neuroretinal rim (NRR) area from OCT and the visual field mean deviation (MD) from Humphrey standard automated perimetry.Table 1 Demographic and clinical characteristics of the glaucoma and healthy subjects.

	Healthy control
(n = 24; mean ± S.E.M.)	Early glaucoma
(n = 14; mean ± S.E.M.)	Advanced glaucoma
(n = 26; mean ± S.E.M.)	Group effect (one-way ANOVA or Chi-squared test)	
Age (year)	64.67 ± 1.56	65.57 ± 2.40	66.19 ± 1.47	F(2,61) = 0.234, P = 0.792	
Sex	11 M, 13F	4M, 10F	13M, 13F	χ2(2) = 1.766, P = 0.414	
Average pRNFL thickness (μm) OD	89.10 ± 2.05	72.00 ± 3.24	62.04 ± 1.95	F(2,56) = 40.923, P < 0.001	
Average pRNFL thickness (μm) OS	91.00 ± 1.84	76.85 ± 3.02	64.62 ± 2.07	F(2,56) = 40.152, P < 0.001	
mGCIPL thickness (μm) OD	79.06 ± 1.44	68.00 ± 3.93	59.52 ± 2.10	F(2,49) = 18.932, P < 0.001	
mGCIPL thickness (μm) OS	78.47 ± 2.00	71.75 ± 2.65	58.57 ± 1.50	F(2,49) = 31.705, P < 0.001	
NRR area (mm2) OD	1.30 ± 0.07	0.86 ± 0.09	0.64 ± 0.04	F(2,56) = 34.963, P < 0.001	
NRR area (mm2) OS	1.35 ± 0.04	0.91 ± 0.11	0.66 ± 0.04	F(2,56) = 41.581, P < 0.001	
Optic nerve head C/D OD	0.42 ± 0.04	0.70 ± 0.04	0.75 ± 0.02	F(2,56) = 30.195, P < 0.001	
Optic nerve head C/D OS	0.41 ± 0.04	0.68 ± 0.03	0.75 ± 0.03	F(2,56) = 35.811, P < 0.001	
Visual field mean deviation (dB) OD	−1.65 ± 0.60	−2.11 ± 0.61	−17.36 ± 1.65	F(2,55) = 48.971, P < 0.001	
Visual field mean deviation (dB) OS	−1.47 ± 0.40	−2.54 ± 0.67	−14.33 ± 1.67	F(2,55) = 31.239, P < 0.001	
pRNFL peripapillary retinal nerve fiber layer, mGCIPL macular ganglion cell-inner plexiform layer, NRR neuroretinal rim, C/D cup-to-disc, OD oculus dexter (right eye), OS oculus sinister (left eye).

Fig. 1 Sample voxel location for proton MRS and representative spectra for GABA and glutamate.

a A 2.2 × 2.2 × 2.2 cm3 voxel (white box) was manually positioned along the calcarine sulci on the occipital lobe. b A sample spectrum from the MEGA-PRESS sequence for GABA. The GABA peak is at 2.8–3.2 ppm. c A sample spectrum from the PRESS sequence for glutamate. The glutamate peak is at 2.2–2.4 ppm.

Glaucoma is accompanied by reduction of GABA in the visual cortex

To examine whether the GABAergic system in the visual cortex alters with disease severity, we first compared the GABA levels across healthy controls, early glaucoma subjects (average MD between eyes better than −6.0 dB) and advanced glaucoma subjects (average MD between eyes worse than −6.0 dB). A one-way ANCOVA while regressing out the effect of age revealed a significant main effect of group (F(2,54) = 6.666, P = 0.003, partial η2 = 0.198; Fig. 2a). Post-hoc tests showed that the levels of GABA are significantly reduced in advanced glaucoma patients compared to healthy controls (early glaucoma vs. healthy controls, Bonferroni-corrected P = 0.330, 95% CI = −0.068 to 0.014, advanced glaucoma vs. healthy controls, Bonferroni-corrected P = 0.002, 95% CI = −0.084 to −0.016, early glaucoma vs. advanced glaucoma, Bonferroni-corrected P = 0.508, 95% CI = −0.018 to 0.064).Fig. 2 GABA levels in the visual cortex.

a The GABA levels (normalized to total creatine levels) were significantly reduced in advanced glaucoma patients compared to healthy controls (Bonferroni-corrected P = 0.002). The distributions are represented using box and whisker plots and the outlier is plotted as a red plus sign. In descending order, the lines in the plots represent: maximum, third quartile, median, first quartile, and minimum. Individual data points are overlaid on the box plots. Healthy control: n = 23, early glaucoma: n = 12, advanced glaucoma: n = 23. b Individual differences in the retinal structure index were significantly associated with individual differences in the GABA levels after controlling for age (r = 0.339, P = 0.020; n = 48).

Next, we also tested whether the impairments on the retina structure, which is a fine-grained, continuous measure of disease severity in the eye are associated with the reduction of GABA levels in the visual cortex. For this, we obtained the marker for retinal structural damage by extracting a common component from each individual’s pRNFL thickness, mGCIPL thickness, C/D ratio, and NRR area of left and right eyes using PCA. The PCA yielded only one component that had an eigenvalue greater than one, explaining 73.42% of variances of all clinical OCT measures. This component, termed as retinal structure index, was significantly different between groups, with the smallest value in the advanced glaucoma subjects and the greatest value in the healthy controls (main effect of group, F(2,48) = 85.042, P < 0.001, partial η2 = 0.780; early glaucoma vs. healthy controls, Bonferroni-corrected P < 0.001, 95% CI = −1.753 to −0.849, advanced glaucoma vs. healthy controls, Bonferroni-corrected P < 0.001, 95% CI = −2.407 to −1.636, early glaucoma vs. advanced glaucoma, Bonferroni-corrected P < 0.001, 95% CI = 0.293 to 1.149). Using this component, we ran a linear regression analysis to estimate the GABA levels. The results showed that the retina structure significantly predicted the GABA levels while controlling for age (T(45) = 2.414, P = 0.020, β = 0.337, partial correlation = 0.339, R2 change = 0.113; Fig. 2b). It indicates that those who have greater impairments on the retina structure (i.e., lower value of the retinal structure index) present lower amounts of GABA in the visual cortex independent of their ages. The effect of age, however, failed to predict the amount of GABA after controlling for the retinal structure index (T(45) = 0.596, P = 0.554, β = 0.083, partial correlation = 0.088, R2 change = 0.007).

Glaucoma is accompanied by reduction of glutamate in the visual cortex

Next, we examined whether the glutamate levels exhibit similar changes as GABA levels in glaucoma. A one-way ANCOVA with a factor of group controlling for age revealed a significant main effect of group (F(2,56) = 5.157, P = 0.009, partial η2 = 0.156; Fig. 3a). The levels of glutamate were significantly lower in advanced glaucoma patients compared to healthy controls and early glaucoma patients (early glaucoma vs. healthy controls, Bonferroni-corrected P = 1.000, 95% CI = −0.110 to 0.153, advanced glaucoma vs. healthy controls, Bonferroni-corrected P = 0.031, 95% CI = −0.236 to −0.009, early glaucoma vs. advanced glaucoma, Bonferroni-corrected P = 0.024, 95% CI = 0.015 to 0.273). Further, the linear regression analysis demonstrated that the retinal structure index significantly predicted the levels of glutamate while the effect of age was controlled (T(46) = 2.654, P = 0.011, β = 0.358, partial correlation=0.364, R2 change = 0.128; Fig. 3b). It suggests that older adults with greater impairments on the retina structure present lower amounts of glutamate in the visual cortex independent of their ages. Nevertheless, age failed to account for the levels of glutamate when the effect of the retina structure was controlled (T(46) = −1.525, P = 0.134, β = −0.206, partial correlation = −0.219, R2 change = 0.042).Fig. 3 Glutamate levels in the visual cortex.

a The glutamate levels (normalized to total creatine levels) were significantly lower in advanced glaucoma patients compared to healthy controls (Bonferroni-corrected P = 0.031) and early glaucoma patients (Bonferroni-corrected P = 0.024). The distributions are represented using box plots and the outliers are plotted as red plus signs. Individual data points are overlaid on the box plots. Healthy control: n = 22, early glaucoma: n = 14, advanced glaucoma: n = 24. b The retinal structure index was positively correlated with the glutamate levels after controlling for age (r = 0.364, P = 0.011; n = 49).

Reduction of GABA, but not glutamate predicts deteriorated neural specificity in the visual cortex

Given the significant reduction of GABA and glutamate in the visual cortex of glaucoma patients, we further tested whether the neural specificity is associated with the levels of GABA and glutamate (see Supplementary Fig. 1 for box plots of the neural specificity). For this, we obtained each individual’s brain activity patterns across V1, V2, ventral V3, and ventral V4, which overlapped with the MRS voxel. During image acquisition, subjects viewed the horizontal and vertical flickering checkerboards in an alternating order. We measured the marker for neural specificity by calculating the fisher z-transformed correlation coefficient between the brain activation patterns representing horizontal and vertical flickering checkerboards. Here, more negative correlation coefficient value indicated a stronger neural specificity.

To investigate the relationship between neural specificity and the reduced levels of GABA and glutamate, we conducted a series of multiple linear regression analyses while controlling for the effect of age and retinal structure index. Here, we adjusted the effect of the retinal structure index in addition to age because the impairments in the retina may block some of the incoming visual signals, resulting in degraded neural specificity. Our regression analyses showed that those with lower GABA levels present weaker neural specificity in the visual cortex when the retinal structure index and age were adjusted (T(42) = −2.165, P = 0.036, β = −0.317, partial correlation = −0.317, R2 change = 0.088). This association between GABA and neural specificity survived even after we additionally controlled for the levels of glutamate (T(39) = −2.143, P = 0.038, β = −0.328, partial correlation = −0.325, R2 change = 0.091) and the gray matter volume of the corresponding visual areas (T(38) = −2.159, P = 0.037, β = −0.334, partial correlation = −0.331, R2 change = 0.094; Fig. 4a). In contrast, other factors such as the glutamate levels (T(38) = 0.439, P = 0.663, β = 0.070, partial correlation=0.071, R2 change = 0.004; Fig. 4b), age (T(38) = 1.508, P = 0.140, β = 0.223, partial correlation = 0.238, R2 change = 0.046; Supplementary Fig. 2a), retinal structure index (T(38) = −1.690, P = 0.099, β = −0.267, partial correlation = −0.264, R2 change = 0.057; Supplementary Fig. 2b), and the gray matter volume of visual areas (T(38) = 0.684, P = 0.498, β = 0.100, partial correlation = 0.110, R2 change = 0.009; Supplementary Fig. 2c) failed to account for the neural specificity. These results suggest that the most important predictor of neural specificity was the GABA levels in the visual cortex.Fig. 4 The relationship between neural specificity, GABA, and glutamate in the visual cortex for the entire sample.

a The GABA levels were significantly correlated with neural specificity after controlling for the glutamate levels, retinal structure index, age, and the gray matter volume of the visual areas (r = −0.331, P = 0.037; n = 44). b The glutamate levels were not associated with neural specificity after adjusting the effects of the GABA levels, retinal structure index, age, and the gray matter volume of the visual areas (r = 0.110, P = 0.498; n = 44).

In a separate analysis, we conducted a hierarchical regression analysis to test whether the GABA levels explained the significant variance beyond that explained by the retinal structure index, age, glutamate levels, and the gray matter volume of visual areas. The results showed that adding GABA levels to the regression model explained an additional 9.4% of the variations in neural specificity and that this change was significant (F(1,38) = 4.663, P = 0.037).

After confirming that GABA levels significantly accounted for the neural specificity in the entire sample, we further tested whether the same association could be found within the samples of glaucoma patients and healthy controls separately. In the sample of glaucoma patients, we replicated the results of an association between GABA levels and neural specificity when the retinal structure index and age were controlled (T(26) = −2.507, P = 0.019, β = −0.410, partial correlation = −0.416, R2 change = 0.166). This association even became stronger when we additionally adjusted the effects of the levels of glutamate (T(24) = −2.886, P = 0.008, β = −0.517, partial correlation = −0.508, R2 change = 0.219) and the gray matter volume (T(23) = −2.870, P = 0.009, β = −0.524, partial correlation = −0.513, R2 change = 0.223; Fig. 5a). A hierarchical regression analysis also showed that an additional 22.3% of the variations in neural specificity can be explained by adding GABA levels to the regression model (F(1,23) = 8.235, P = 0.009).Fig. 5 The relationship between neural specificity, GABA, glutamate, retinal structure index, and age within the glaucoma group.

a The GABA levels significantly predicted neural specificity after controlling for the glutamate levels, retinal structure index, age, and the gray matter volume (r = −0.513, P = 0.009; n = 29). b The glutamate levels were not associated with neural specificity after controlling for the GABA levels, retinal structure index, age, and the gray matter volume (r = 0.323, P = 0.116; n = 29). c The retinal structure index was not associated with neural specificity after controlling for the GABA levels, glutamate levels, age, and the gray matter volume (r = −0.365, P = 0.072; n = 29). d Age was significantly correlated with neural specificity after adjusting the effects of the GABA levels, glutamate levels, retinal structure index, and the gray matter volume (r = 0.478, P = 0.016; n = 29).

Nevertheless, other factors such as the glutamate levels (T(23) = 1.635, P = 0.116, β = 0.324, partial correlation = 0.323, R2 change = 0.072; Fig. 5b), retina structure index (T(23) = −1.883, P = 0.072, β = −0.341, partial correlation = −0.365, R2 change = 0.096; Fig. 5c), and the gray matter volume of visual areas (T(23) = 0.489, P = 0.630, β = 0.085, partial correlation = 0.101, R2 change = 0.006; Supplementary Fig. 3) failed to predict the neural specificity. This non-significant effect of the glutamate levels, retina structure index, and the gray matter volume of visual areas in the glaucoma group is consistent with the findings from the entire sample including both glaucoma and healthy subjects.

A discrepancy between the sample of glaucoma patients and the entire sample was seen in the effects of age on neural specificity. Glaucoma patients with older ages presented weaker neural specificity while controlling for other confounding effects including retinal structure index, GABA levels, glutamate levels, and the gray matter volume of the visual areas (T(23) = 2.610, P = 0.016, β = 0.496, partial correlation = 0.478, R2 change = 0.185; Fig. 5d). This result indicates that age was a meaningful predictor for neural specificity within the sample of glaucoma patients, but not in the entire sample including both glaucoma and healthy subjects.

Finally, we examined whether the same patterns of relationship could be seen within the healthy control group. The results showed that none of GABA, glutamate, age, or retinal structure index was a significant predictor for the neural specificity in healthy controls after adjusting the remaining factors (GABA: T(9) = −0.594, P = 0.567, β = −0.155, partial correlation = −0.194, R2 change = 0.022; Glutamate: T(9) = 0.310, P = 0.763, β = 0.087, partial correlation = 0.103, R2 change = 0.006; Age: T(9) = −0.455, P = 0.660, β = −0.119, partial correlation = −0.150, R2 change = 0.013; Retinal structure index: T(9) = −0.231, P = 0.823, β = −0.060, partial correlation = −0.077, R2 change = 0.003; Supplementary Fig. 4a–d). However, the gray matter volume of visual areas was a significant predictor for the neural specificity in healthy controls (T(9) = 2.444, P = 0.037, β = 0.659, partial correlation = 0.632, R2 change = 0.375; Supplementary Fig. 4e).

Critically, the overall MRS results yielded by using total creatine for normalization were maintained even when we used N-acetyl-aspartate (NAA), another standard reference resonance for normalization. Using NAA for normalization, we still observed a marginal association between the GABA levels and the neural specificity after adjusting the effects of the retinal structure index, age, the levels of glutamate normalized by NAA, and the gray matter volume in the entire sample (T(38) = −1.890, P = 0.066, β = −0.286, partial correlation = −0.293, R2 change = 0.073). In the sample of glaucoma patients, this association between the GABA level normalized by NAA and the neural specificity became even stronger (T(23) = −2.659, P = 0.014, β = −0.497, partial correlation = −0.485, R2 change=0.192).

Discussion

In the current study, we demonstrated that both GABA and glutamate levels in the visual cortex decline with increasing glaucoma severity. This association between impairments of the retina structure, GABA, and glutamate was independent of age among our older adults. Further, we showed that the reduction of GABA, but not that of glutamate, was associated with the neural specificity in the visual cortex after controlling for confounding effects including age and retinal structural damage. These results provide clear evidence that glaucomatous neurodegeneration involves the reduction of GABA and glutamate levels in the visual cortex and that the reduction of GABA, but not that of glutamate, plays a critical role in degrading neural specificity.

The reduction of GABA levels in the visual cortex of glaucoma patients is consistent with the literature showing that GABA metabolism is a unique pathway associated with glaucoma50. Furthermore, GABA is susceptible to the accumulation of amyloid β and tau24–27,53. Prior studies demonstrated that glaucoma involves accumulation of tau in the vitreous fluid17, retina16,18, lateral geniculate nucleus19 as well as deposits of amyloid β in the lateral geniculate nucleus and the primary visual cortex19. Notably, exposure to amyloid β and tau is known to impair the GABAergic systems24–27,53. In animal models of Alzheimer’s disease, accumulation of amyloid β led to decreases in the number26 and density of GABAergic neurons before loss of pyramidal cells27. Further, amyloid β was shown to weaken synaptic inhibition via downregulating GABAA receptors24. In another animal model expressing high deposits of tau, a substantial decline was observed in the GABA-synthetic enzyme and GABAergic interneurons, which co-localized with tau markers25. This line of studies raises a possibility that glaucoma may share common pathogenic mechanisms with Alzheimer’s disease and that accumulation of amyloid β and tau in the visual pathway may promote loss of GABA neurotransmission in glaucoma patients.

In addition to the possible impact of amyloid β and tau on GABA, several other factors may contribute to GABA reduction. First, impaired visual input may facilitate GABA reduction. Visual deprivation was shown to trigger GABA reduction in the primary visual cortex54. Given that the visual signals are gradually blocked as glaucoma progresses, it is speculated that its impact on GABA would increase as well. Consistent with this, our study observed that the GABA levels in the visual cortex measured by MRS decrease with increasing impairments of the retina structure. Second, accelerated aging process in glaucoma may affect GABA levels. Facilitated aging in glaucoma is supported by the finding that glaucoma patients present elevated amount of advanced glycation end products, which are inevitable compounds of the aging process, as well as upregulation of their receptors in the retina and optic nerve head43. Relatedly, the aging process involves impairments of GABAergic systems55. In healthy animal models and humans, age-related declines were observed in terms of the number of GABAergic neurons56–58, baseline GABA levels, GABA release, GABA receptor binding57, and MRS GABA measures throughout the brain44–49. Therefore, it is possible that the age-related declines of GABA may be accelerated in glaucoma patients. Future studies could address this possibility by obtaining the measures of biological speed of aging from biological assays along with MRS and functional MRI.

In this study, MRS measures of GABA cannot distinguish whether the decreased GABA levels in the visual cortex of glaucoma patients reflect reduction of GABAergic neurons, GABA neurotransmitters, or GABA synthesis59. Any combination of these changes could lead to reduced MRS GABA levels. While it is yet unknown which factor contributed the most to our MRS observations, one study using an experimental monkey model observed that after unilateral IOP elevation, the GABAA receptor protein in the primary visual cortex was reduced11, possibly driven by suppression of visually driven activity60. This finding suggests that the declined MRS GABA levels that we observed in the current study are in part reflective of reduced GABA receptor binding. Other possibilities such as reduction of the GABAergic neurons or GABA neurotransmitters need to be further identified in future in vivo, in situ, and ex vivo studies using a finer spatial resolution. Further, it should be noted that in the current study, the GABA-editing technique did not specifically suppress the macromolecules, although macromolecule signals were taken into account in the analysis. Future studies can consider macromolecule-suppressed GABA-edited techniques in order to further reduce the potential contaminations of GABA signals with the macromolecules61.

The decreased GABA signals in the visual cortex of glaucoma patients suggest that the cortical GABAergic function is likely impaired in glaucoma. GABAergic signals are critical for shaping the response of a single neuron as well as the response patterns of a population of neurons28. At the level of a single neuron, GABA-mediated inhibition can sharpen the neuronal response around the preferred stimulus feature62. In line with this, application of GABA agonist increases the neuron’s selectivity to preferred stimulus feature, while application of GABA antagonist decreases the neuron’s selectivity63. Therefore, it is speculated that the selectivity of neurons in the visual cortex may be undermined in glaucoma patients. At the population level, GABAergic signals are related to the specificity of neuronal activity patterns in response to different stimuli. With lower amount of GABA, neuronal activity patterns become less distinct or more confusable29,30. Consistent with these prior findings, we observed that those who had lower GABA levels displayed less specific neural activity patterns in the visual cortex independent of the age and impairments of the retina structure. It suggests that the declines of GABA and neural specificity in the visual cortex may negatively affect the downstream processes along the lower- to higher-order brain regions32,33 in glaucoma.

Alterations of GABA levels and neural specificity in glaucoma patients are likely to have functional influences to various aspects of brain functions. Prior studies showed that GABA levels in healthy subjects are associated with visual-spatial intelligence, visual surround suppression64, motor inhibition45, motor learning65, general sensorimotor functioning49, cognitive performance46,48, experience-dependent plasticity66, and sleep67. Further, neural specificity was found to be related to memory performance32 and fluid processing ability33. Therefore, it is reasonable to suspect that these aspects of brain functions may be undermined if the GABA level and the neural specificity are impaired. In line with this speculation, glaucoma patients were found to present abnormalities in perceptual, cognitive functions and sleep quality. Specifically, these abnormalities include elevated visual crowding effect34, impairments on perceptual switch during binocular rivalry35–37, declines in visual categorization39,40, visual attention41 and cognitive capacity39,40,42 as well as sleep disturbance68,69. Thus, we propose that there could be a functional link between GABA reduction and the reported abnormalities in glaucoma patients. Future studies could address this question by employing both MRS and behavioral measures.

In addition to the GABA levels, we observed that the years of age could predict the neural specificity within the sample of glaucoma patients. Specifically, lower neural specificity was associated with older ages. This association is reasonable given that aging is accompanied by degraded neural distinctiveness/discriminability29–31,70,71. However, in the sample of healthy subjects alone, the relationships between GABA, years of age and neural specificity were not observed. This discrepancy could be partly due to the smaller sample size of the healthy controls or could be related to accelerated biological aging process in glaucoma patients. In order to fully understand this discrepancy, we need future studies with a greater sample size.

Another important finding in this study is that the glutamate levels in the visual cortex decline with the impairments on the retinal structure. This observation may appear inconsistent with the literature showing that accumulation of amyloid β leads to increased glutamate levels in the extracellular space20,21. Nevertheless, it should be noted that glutamate and GABA are closely linked to each other in terms of production, function, and reuptake. At the level of production, GABA is synthesized by decarboxylation of glutamate via glutamate acid decarboxylase. At the functional level, the balance between excitation and inhibition is critical for a neuron to function properly. The neurons become dysfunctional when cortical excitation or inhibition is acutely manipulated72. Further, increases or decreases in the strength of cortical excitation are naturally accompanied by the corresponding proportional changes in the strength of inhibition73,74. At the level of reuptake, glutamate uptake strengthens inhibitory synapses75. Thus, it is possible that reduced glutamate levels in the visual cortex may be at least partially related to the reduced GABA levels. Further, recent studies have demonstrated that amyloid β and glutamatergic signals have multifaceted mechanisms76. In the early stage of Alzheimer’s disease, glutamatergic activity is elevated leading to hyperexcitability of neurons77,78. However, as Alzheimer’s disease progresses, glutamate release and vesicular glutamate transporters become markedly decreased79. Prior MRS studies also observed a decreasing pattern of glutamate levels with severity of Alzheimer’s disease80,81. Thus, our results of decreased glutamate levels in glaucoma may resemble these progressive glutamate changes in the late phases of Alzheimer’s disease.

In the current study, the glutamate levels were not associated with the neural specificity. While it is not clear why the glutamate levels are not linked to the local brain activity, prior studies have consistently shown that the glutamate levels measured during rest are not correlated with the intra-regional brain activity during tasks82–87. Meta-analyses also found evidence in favor of no association between glutamate levels and local brain activity88. This is in sharp contrast with the GABA levels which showed largely negative associations with the local brain activity82,83,85,89,90. Up to date, two explanations have been proposed to interpret the absence of association between glutamate levels and local brain activity. First, the glutamate levels measured during rest may be related to the resting-state brain activity, rather than task-related brain activity91. Second, glutamatergic neurons detected by MRS may influence the brain activity in distal regions more than the local region through their long-range projections88. This interpretation is based on prior studies showing that the glutamate levels are positively correlated with the brain activity outside of the MRS sampling region92–94, but not within the sampling region82–87. However, these explanations remain speculative and thus, need to be further examined in future studies.

While the current study focused on the visual cortex only, it is worth knowing whether the same neurochemical changes occur in the retina of glaucoma patients. A number of patient and experimental animal studies have examined the glutamate-mediated excitotoxicity in the vitreous body and retina in glaucoma95–98. Nevertheless, these studies provided disparate results, making it unclear whether the glutamate in the retina is elevated under glaucoma. Other studies investigating GABA, however, point to a dysfunction of the GABAergic signals in glaucoma. For example, in animal models of glaucoma, the retina showed decreases in the GABA turnover rate, activity of GABA-synthetic enzyme, GABA release and the expression of retinal genes associated with GABAergic systems51,99. Further, studies reported that modulating GABAergic activity in the rat’s retina under chronic glaucomatous conditions may improve retinal ganglion cell viability and function100,101. Future studies are warranted to examine whether GABA plays a role in glaucomatous neuroprotection and functional recovery throughout the visual system of the brain.

To conclude, our findings demonstrate a tight link between declines of GABA and neural specificity in the visual cortex of glaucoma patients. Our results suggest that glaucoma-specific declines of GABA undermine neural specificity in the visual cortex and that targeting GABA could improve the neural specificity in glaucoma patients.

Methods

Subjects

The Institutional Review Board of New York University Grossman School of Medicine approved this study. This study adhered to the tenets of the Declaration of Helsinki. Informed consent was obtained from all subjects prior to participation. Forty glaucoma subjects [age = 65.98 ± 1.26 (mean ± S.E.M.); 42.5% male] and twenty-four healthy controls [age = 64.67 ± 1.56 (mean ± S.E.M.); 45.8% male] were included in the study between August 2018 and May 2022 at New York University Langone Health’s Department of Ophthalmology. All subjects had a best corrected visual acuity (BCVA) of 20/60 or better and showed no past medical history (PMH) or current evidence of retinal or neurological disorders. The demographic and clinical information of those who were included in this study is depicted in Table 1. Subjects in the disease group were clinically diagnosed with primary glaucoma, whereas healthy controls exhibited no clinical evidence of glaucomatous conditions. Subjects were not allowed to participate in the MRI study if they were pregnant or breastfeeding at the time of the study, had any metal parts or fragments in the body with the exception of dental fillings, had conditions such as anxiety or claustrophobia, or had obesity that would hinder placement into the MRI scanner.

Clinical ophthalmic exams

The clinical ophthalmic data for both glaucoma patients and healthy subjects were collected, which included pRNFL thickness, mGCIPL thickness, optic nerve head C/D ratio, and NRR area through an automatic analytic software on board of the Cirrus spectral-domain OCT device (Zeiss, Dublin, CA, USA). Visual field MD was also obtained from the Humphrey Swedish Interactive Thresholding Algorithm (SITA) 24-2 standard (Zeiss, Dublin, CA, USA). The average MD of the left and right eyes combined (OU) were obtained to assign the patients into early or advanced glaucoma group. Early glaucoma was categorized as glaucoma patients with average OU MD better than −6.0 dB and advanced glaucoma was determined as average OU MD worse than −6.0 dB following the previous literature102–105.

MRI data acquisition

Age-matched healthy subjects and glaucoma patients were scanned inside a 3-Tesla MR Prisma scanner (Siemens, Germany) with a 20-channel head coil at the Center for Biomedical Imaging, NYU Langone health, New York University. For anatomical reconstruction, high-resolution T1-weighted MR images were acquired using a multi-echo magnetization-prepared rapid gradient echo sequence, with 256 slices, voxel size = 0.8 × 0.8 × 0.8 mm, 0-mm slice gap, repetition time (TR) = 2400 ms, echo time (TE) = 2.24 ms, flip angle = 8°, field of view = 256 mm, and bandwidth = 210 Hz per pixel.

For MRS acquisition (see Supplementary Table. 1 for details), we manually positioned a 2.2 × 2.2 × 2.2 cm3 voxel along the calcarine sulci in the occipital cortex. The voxel covered parts of the visual cortex including V1, V2, ventral V3, and ventral V4. Shimming was performed by a vendor-provided automated shim tool followed by manual fine adjustment. The shim value defined by the full width at half maximum of the water peak was 13.02 ± 0.13 (mean ± S.E.M.) for gamma-aminobutyric acid (GABA) and 12.93 ± 0.14 (mean ± S.E.M.) for glutamate. The concentration of GABA was measured from a voxel using a MEshcher-GArwood- Point RESolved Spectroscopy (MEGA-PRESS) sequence with double-banded pulses, at TR = 1500 ms, TE = 68 ms, number of averages = 172, and scanning duration = 522 s. The final spectrum was obtained by subtracting the ‘edit-off’ spectrum from ‘edit-on’ spectrum. The concentration of glutamate was measured from the same voxel using a Point RESolved Spectroscopy (PRESS) sequence, with TR = 3000 ms, TE = 30 ms, number of averages = 99, and scanning duration = 300 s. During both MEGA-PRESS and PRESS scans, subjects performed a fixation task. For two subjects, MRS acquisition could not be completed due to technical issues. Additionally, 30 out of the 64 subjects (15 healthy controls and 15 glaucoma patients) underwent water-unsuppressed PRESS scan with number of averages = 16, and scanning duration = 48 s immediately before the PRESS scan in order to perform water scaling.

After MRS acquisition, functional MR images were acquired using a gradient-echo echo-planar imaging (EPI) sequence, with voxel size = 2.3 × 2.3 × 2.3 mm, TR = 1000 ms, TE = 32.60 ms, and scanning duration = 300 s. During fMRI scan, subjects performed a fixation task while a flickering checkerboard pattern was presented at the horizontal vs. vertical meridians. For four subjects (3 healthy controls and 1 glaucoma patient), fMRI images could not be obtained due to technical issues.

Fixation task

To maintain subjects’ attention and vigilance levels across the scans, we provided subjects with the fixation task during MEGA-PRESS and PRESS scans. The fixation point was presented at the center of the gray background and subjects were asked to fixate their eyes at the point. The fixation point changed its color unpredictably from white ([R, G, B] = [255, 255, 255] to red ([R, G, B] = [255, 127, 127]) and returned to white 1.5 s later. When the button was pressed within a 1.5-s time window, this response was recorded as a hit. If not, the response was recorded as a miss. The mean accuracy (±S.E.M.) was 98.87 ± 0.58%. The accuracy did not differ across healthy controls, and early and advanced glaucoma patients (F(2,32) = 0.263, P = 0.770, partial η2 = 0.016).

Neural specificity task

During fMRI scans (1 run=300 s), subjects performed the same fixation task, while the background was filled with flickering checkerboard patterns. The checkerboard patterns were presented at either horizontal or vertical meridians in a gray background. The purpose of adding the flickering checkerboard patterns in the background was to obtain the neural specificity for horizontal and vertical meridians. Each of the horizontal and vertical meridians was presented for 8 s in alternation, with 18 trials for horizontal meridians and 18 trials for vertical meridians in total. At the first and the last 6-s periods, the fixation point was presented only without any checkerboard patterns. The mean accuracy of the fixation task (±S.E.M.) was 94.10 ± 1.11% and was comparable across healthy controls, and early and advanced glaucoma patients (F(2,42) = 1.468, P = 0.242, partial η2 = 0.065).

MRS data analysis

MRS data was fitted in the frequency domain using the LCModel software106. We used basis functions that include models of macromolecular spectra to adjust the baseline generated by macromolecular and lipid components. We assessed the quality of fitting by visually inspecting LCmodel fitted spectra and examining their Cramer-Rao lower bounds (CRLB) and spectral signal-to-noise (S/N) ratio. Four spectra from the MEGA-PRESS and two spectra from the PRESS scans were excluded from further analyses because of poor fitting (CRLB > 20%, S/N < 8). The CRLB was 8.76 ± 0.15% (mean ± S.E.M.) for GABA and 8.65 ± 0.40% (mean ± S.E.M.) for glutamate. The S/N ratio was 23.84 ± 0.38 (mean ± S.E.M.) for GABA and 26.6 ± 1.36 (mean ± S.E.M.) for glutamate. We normalized the concentrations of GABA and glutamate using the amount of total creatine, which is commonly used as a standard reference resonance106,107. For complementary support, we also normalized the GABA and glutamate concentrations by the amount of N-acetyl-aspartate (NAA). In order to examine whether the concentrations of total creatine and NAA were stable regardless of glaucoma, we also normalized the total creatine and NAA to the water signal using the water-unsuppressed PRESS scans collected from 30 out of the 64 subjects. The absolute values of total creatine and NAA did not differ between groups (total creatine, F(2,27) = 0.465, P = 0.633, partial η2 = 0.033; NAA, (F(2,27) = 1.431, P = 0.257, partial η2 = 0.096, one-way ANOVA) and did not have any association with the retina structure index (total creatine, T(23) = −1.152, P = 0.261, β=−0.234, R2 = 0.055; NAA, T(23) = 1.174, P = 0.253, β = 0.238, R2 = 0.056, linear regression). These findings indicate that the total creatine and NAA were not substantially affected by glaucoma, and can be used as reliable references.

fMRI data analysis

We analyzed fMRI data using Freesurfer software version 7.2 (http://surfer.nmr.mgh.harvard.edu/) and MATLAB. The fMRI data was preprocessed with motion correction but not with spatial or temporal smoothing. Then the functional data was registered to the individual’s structural template. We extracted the blood-oxygenation-level-dependent (BOLD) signals from the individual’s region-of-interest (ROI) masks generated from Freesurfer’s cortical reconstruction process. The ROI masks that we used were four cytoarchitectonic areas in the occipital lobe (hOc1-hOc4v) which spatially overlapped with the location of the MRS voxel. The extracted BOLD signals were then shifted by 6 s to account for the hemodynamic delay. We excluded voxels with spikes greater than 10 standard deviations from the mean and removed the linear trend in the BOLD time course. Then we normalized the BOLD signals using z-score for each voxel. All volumes were included in the analysis.

Calculation of neural specificity

We defined the neural specificity as the dissimilarity in brain responses between different visual conditions. For this, we averaged the z-normalized BOLD signals of each voxel across 8 volumes (8 s) which corresponded to the duration of horizontal/vertical meridian presentation. These BOLD signals of each voxel were then averaged across 18 trials of horizontal meridians or 18 trials of vertical meridians separately. Given that these two patterns of BOLD signals represented the horizontal or vertical meridians, we computed the Fisher z-transformed correlation coefficient between these two representations to obtain multi-voxel pattern similarity. Finally, we averaged the correlation coefficients across V1, V2, ventral V3, and ventral V4.

Statistics and reproducibility

The sample size was not predetermined but similar or even greater to those reported in prior studies2,3,10,12,13,108. The investigators were not blinded to disease diagnosis during data collection. For all statistical comparisons, we conducted two-tailed parametric tests with P < 0.05 as the criterion for statistical significance. For ANOVAs, we assessed the assumption of homogeneity of variances using Levene’s test. The following post-hoc tests were conducted with Bonferroni corrections. For multiple regression analyses, we confirmed that the assumptions of independence of observations and non-multicollinearity were not violated using the Durbin-Watson statistic and tolerance/VIF values. Further, we extracted a common component from clinical ophthalmic measures including pRNFL thickness, mGCIPL thickness, optic nerve head C/D ratio, and NRR area for each individual using PCA. Then we entered this component, a measure of retina structure in the regression model rather than entering all clinical ophthalmic measures because the clinical measures were highly correlated to each other, violating the non-multicollinearity. We set the criteria of selecting the component in the PCA as the eigenvalue higher than one. To verify that the overall MRS results hold regardless of the method of normalization, we used NAA in addition to total creatine for normalization.

Apparatus

We created all visual stimuli in MATLAB using Psychophysics Toolbox 3109. The stimuli were presented via an MRI-compatible projector (1024 × 768 resolution, 60 Hz refresh rate).

Inclusion and ethics statement

The authors followed the recommendations set out in the Global Code of Conduct for Research in Resource-Poor Settings. The research was determined and conducted by local researchers. The roles and responsibilities were agreed amongst collaborators ahead of the research. This research did not involve any of the health, safety, security or other risk to researchers, and did not have any severe restrictions in the setting of the researchers. This study did not result in any of the stigmatization, incrimination, discrimination or personal risk. The authors took local and regional research relevant to the current study into account in citations.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Reporting Summary

Supplementary information

The online version contains supplementary material available at 10.1038/s42003-023-04918-8.

Acknowledgements

The authors would like to thank Ms Tonya Robins, Zena Moore, Jamika Singleton-Garvin and members of the Neuroimaging and Visual Science Laboratory at New York University Grossman School of Medicine for their help with subject recruitment and technical support. This work is supported in part by the National Institutes of Health R01-EY028125, R01-EY013178, and P41-EB017183 (Bethesda, Maryland), BrightFocus Foundation G2016030, G2019103, and G2021001F (Clarksburg, Maryland), and an unrestricted grant from Research to Prevent Blindness to NYU Langone Health Department of Ophthalmology (New York, NY).

Author contributions

J.W.B., C.P., G.W., J.S.S., and K.C.C. conceived and designed experiments. J.W.B. and C.P. carried out experiments. J.W.B., C.P., and K.Y. analyzed the data. G.W., J.S.S., and K.C.C. supervised the experiments. J.W.B., C.P., K.Y., G.W., J.S.S., and K.C.C. wrote the article.

Peer review

Peer review information

Communications Biology thanks Molly Simmonite and Jihong Wu for their contribution to the peer review of this work. Primary Handling Editors: Joao Valente.

Data availability

All data underlying the results and figures are freely available in OSF with the identifier doi: 10.17605/OSF.IO/85CDS (https://osf.io/85cds/)110

Code availability

All statistical codes for SPM are freely available in OSF with the identifier doi: 10.17605/OSF.IO/85CDS (https://osf.io/85cds/)110

Competing interests

J.S.S. declares the following competing interests: Royalties e Zeiss, Dublin, CA.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Tham YC Global prevalence of glaucoma and projections of glaucoma burden through 2040: a systematic review and meta-analysis Ophthalmology 2014 121 2081 2090 10.1016/j.ophtha.2014.05.013 24974815
2. Garaci FG Optic nerve and optic radiation neurodegeneration in patients with glaucoma: in vivo analysis with 3-T diffusion-tensor MR imaging Radiology 2009 252 496 501 10.1148/radiol.2522081240 19435941
3. Zikou AK Voxel-based morphometry and diffusion tensor imaging of the optic pathway in primary open-angle glaucoma: a preliminary study AJNR Am. J. Neuroradiol. 2012 33 128 134 10.3174/ajnr.A2714 22116110
4. Weber AJ Chen H Hubbard WC Kaufman PL Experimental glaucoma and cell size, density, and number in the primate lateral geniculate nucleus Invest Ophthalmol. Vis. Sci. 2000 41 1370 1379 10798652
5. Yucel YH Zhang Q Gupta N Kaufman PL Weinreb RN Loss of neurons in magnocellular and parvocellular layers of the lateral geniculate nucleus in glaucoma Arch. Ophthalmol. 2000 118 378 384 10.1001/archopht.118.3.378 10721961
6. Crawford ML Harwerth RS Smith EL 3rd Shen F Carter-Dawson L Glaucoma in primates: cytochrome oxidase reactivity in parvo- and magnocellular pathways Invest Ophthalmol. Vis. Sci. 2000 41 1791 1802 10845600
7. Imamura K Molecular imaging reveals unique degenerative changes in experimental glaucoma Neuroreport 2009 20 139 144 10.1097/WNR.0b013e32831d7f82 19057418
8. Chaturvedi N Hedley-Whyte ET Dreyer EB Lateral geniculate nucleus in glaucoma Am. J. Ophthalmol. 1993 116 182 188 10.1016/S0002-9394(14)71283-8 8352303
9. Gupta N Ang LC Noel de Tilly L Bidaisee L Yucel YH Human glaucoma and neural degeneration in intracranial optic nerve, lateral geniculate nucleus, and visual cortex Br. J. Ophthalmol. 2006 90 674 678 10.1136/bjo.2005.086769 16464969
10. Engelhorn T Diffusion tensor imaging detects rarefaction of optic radiation in glaucoma patients Acad. Radio. 2011 18 764 769 10.1016/j.acra.2011.01.014
11. Lam DY Kaufman PL Gabelt BT To EC Matsubara JA Neurochemical correlates of cortical plasticity after unilateral elevated intraocular pressure in a primate model of glaucoma Invest. Ophthalmol. Vis. Sci. 2003 44 2573 2581 10.1167/iovs.02-0779 12766059
12. Boucard CC Changes in cortical grey matter density associated with long-standing retinal visual field defects Brain 2009 132 1898 1906 10.1093/brain/awp119 19467992
13. Murphy MC Retinal Structures and Visual Cortex Activity are Impaired Prior to Clinical Vision Loss in Glaucoma Sci. Rep. 2016 6 31464 10.1038/srep31464 27510406
14. Ghiso JA Doudevski I Ritch R Rostagno AA Alzheimer’s disease and glaucoma: mechanistic similarities and differences J. Glaucoma 2013 22 S36 S38 10.1097/IJG.0b013e3182934af6 23733125
15. Sen S Saxena R Tripathi M Vibha D Dhiman R Neurodegeneration in Alzheimer’s disease and glaucoma: overlaps and missing links Eye (Lond.) 2020 34 1546 1553 10.1038/s41433-020-0836-x 32152519
16. Gupta N Fong J Ang LC Yucel YH Retinal tau pathology in human glaucomas Can. J. Ophthalmol. 2008 43 53 60 10.3129/i07-185 18219347
17. Yoneda S Vitreous fluid levels of beta-amyloid((1-42)) and tau in patients with retinal diseases Jpn. J. Ophthalmol. 2005 49 106 108 10.1007/s10384-004-0156-x 15838725
18. Chiasseu M Tau accumulation, altered phosphorylation, and missorting promote neurodegeneration in glaucoma J. Neurosci. 2016 36 5785 5798 10.1523/JNEUROSCI.3986-15.2016 27225768
19. Yan Z Elevated intraocular pressure induces amyloid-beta deposition and tauopathy in the lateral geniculate nucleus in a monkey model of glaucoma Invest. Ophthalmol. Vis. Sci. 2017 58 5434 5443 10.1167/iovs.17-22312 29059309
20. Fernandez-Tome P Brera B Arevalo MA de Ceballos ML Beta-amyloid25-35 inhibits glutamate uptake in cultured neurons and astrocytes: modulation of uptake as a survival mechanism Neurobiol. Dis. 2004 15 580 589 10.1016/j.nbd.2003.12.006 15056466
21. Harris ME Amyloid beta peptide (25-35) inhibits Na+-dependent glutamate uptake in rat hippocampal astrocyte cultures J. Neurochem. 1996 67 277 286 10.1046/j.1471-4159.1996.67010277.x 8667003
22. Chin JH Ma L MacTavish D Jhamandas JH Amyloid beta protein modulates glutamate-mediated neurotransmission in the rat basal forebrain: involvement of presynaptic neuronal nicotinic acetylcholine and metabotropic glutamate receptors J. Neurosci. 2007 27 9262 9269 10.1523/JNEUROSCI.1843-07.2007 17728440
23. Kabogo D Rauw G Amritraj A Baker G Kar S ss-amyloid-related peptides potentiate K+-evoked glutamate release from adult rat hippocampal slices Neurobiol. Aging 2010 31 1164 1172 10.1016/j.neurobiolaging.2008.08.009 18819729
24. Ulrich D Amyloid-beta impairs synaptic inhibition via GABA(A) Receptor Endocytosis J. Neurosci. 2015 35 9205 9210 10.1523/JNEUROSCI.0950-15.2015 26085642
25. Levenga J Tau pathology induces loss of GABAergic interneurons leading to altered synaptic plasticity and behavioral impairments Acta Neuropathol. Commun. 2013 1 34 10.1186/2051-5960-1-34 24252661
26. Krantic S Hippocampal GABAergic neurons are susceptible to amyloid-beta toxicity in vitro and are decreased in number in the Alzheimer’s disease TgCRND8 mouse model J. Alzheimers Dis. 2012 29 293 308 10.3233/JAD-2011-110830 22232004
27. Ramos B Early neuropathology of somatostatin/NPY GABAergic cells in the hippocampus of a PS1xAPP transgenic model of Alzheimer’s disease Neurobiol. Aging 2006 27 1658 1672 10.1016/j.neurobiolaging.2005.09.022 16271420
28. Isaacson JS Scanziani M How inhibition shapes cortical activity Neuron 2011 72 231 243 10.1016/j.neuron.2011.09.027 22017986
29. Chamberlain JD GABA levels in ventral visual cortex decline with age and are associated with neural distinctiveness Neurobiol. Aging 2021 102 170 177 10.1016/j.neurobiolaging.2021.02.013 33770531
30. Lalwani P Neural distinctiveness declines with age in auditory cortex and is associated with auditory GABA levels Neuroimage 2019 201 116033 10.1016/j.neuroimage.2019.116033 31326572
31. Simmonite M Polk TA Age-related declines in neural distinctiveness correlate across brain areas and result from both decreased reliability and increased confusability Neuropsychol. Dev. Cogn. B: Aging Neuropsychol. Cogn. 2022 29 483 499 10.1080/13825585.2021.1999383 34757860
32. Fandakova Y Leckey S Driver CC Bunge SA Ghetti S Neural specificity of scene representations is related to memory performance in childhood Neuroimage 2019 199 105 113 10.1016/j.neuroimage.2019.05.050 31121295
33. Park J Carp J Hebrank A Park DC Polk TA Neural specificity predicts fluid processing ability in older adults J. Neurosci. 2010 30 9253 9259 10.1523/JNEUROSCI.0853-10.2010 20610760
34. Ogata NG Visual crowding in glaucoma Invest Ophthalmol. Vis. Sci. 2019 60 538 543 10.1167/iovs.18-25150 30716149
35. Tarita-Nistor L Samet S Trope GE Gonzalez EG Intra- and inter-hemispheric processing during binocular rivalry in mild glaucoma PLoS ONE 2020 15 e0229168 10.1371/journal.pone.0229168 32097443
36. Issashar Leibovitzh G Trope GE Buys YM Tarita-Nistor L Perceptual grouping during binocular rivalry in mild glaucoma Front Aging Neurosci. 2022 14 833150 10.3389/fnagi.2022.833150 35693345
37. Joao CAR Scanferla L Jansonius NM Binocular interactions in glaucoma patients with nonoverlapping visual field defects: contrast summation, rivalry, and phase combination Invest. Ophthalmol. Vis. Sci. 2021 62 9 10.1167/iovs.62.12.9 34505864
38. Tarita-Nistor L Samet S Trope GE Gonzalez EG Dominance wave propagation during binocular rivalry in mild glaucoma Vis. Res. 2019 165 64 71 10.1016/j.visres.2019.10.006 31678616
39. Lenoble Q Lek JJ McKendrick AM Visual object categorisation in people with glaucoma Br. J. Ophthalmol. 2016 100 1585 1590 10.1136/bjophthalmol-2015-308289 27503395
40. Roux-Sibilon A Scene and human face recognition in the central vision of patients with glaucoma PLoS ONE 2018 13 e0193465 10.1371/journal.pone.0193465 29481572
41. Swenor BK Impact of the ability to divide attention on reading performance in glaucoma Invest Ophthalmol. Vis. Sci. 2017 58 2456 2462 10.1167/iovs.17-21520 28460047
42. Gangeddula V Ranchet M Akinwuntan AE Bollinger K Devos H Effect of cognitive demand on functional visual field performance in senior drivers with glaucoma Front. Aging Neurosci. 2017 9 286 10.3389/fnagi.2017.00286 28912712
43. Tezel G Luo C Yang X Accelerated aging in glaucoma: immunohistochemical assessment of advanced glycation end products in the human retina and optic nerve head Invest. Ophthalmol. Vis. Sci. 2007 48 1201 1211 10.1167/iovs.06-0737 17325164
44. Chalavi S The neurochemical basis of the contextual interference effect Neurobiol. Aging 2018 66 85 96 10.1016/j.neurobiolaging.2018.02.014 29549874
45. Hermans L Brain GABA levels are associated with inhibitory control deficits in older adults J. Neurosci. 2018 38 7844 7851 10.1523/JNEUROSCI.0760-18.2018 30064995
46. Simmonite M Age-related declines in occipital GABA are associated with reduced fluid processing ability Acad. Radio. 2019 26 1053 1061 10.1016/j.acra.2018.07.024
47. Gao F Edited magnetic resonance spectroscopy detects an age-related decline in brain GABA levels Neuroimage 2013 78 75 82 10.1016/j.neuroimage.2013.04.012 23587685
48. Porges EC Frontal gamma-aminobutyric acid concentrations are associated with cognitive performance in older adults Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2017 2 38 44 28217759
49. Cassady K Sensorimotor network segregation declines with age and is linked to GABA and to sensorimotor performance Neuroimage 2019 186 234 244 10.1016/j.neuroimage.2018.11.008 30414983
50. Bailey JN Hypothesis-independent pathway analysis implicates GABA and acetyl-CoA metabolism in primary open-angle glaucoma and normal-pressure glaucoma Hum. Genet. 2014 133 1319 1330 10.1007/s00439-014-1468-7 25037249
51. Moreno MC Effect of ocular hypertension on retinal GABAergic activity Neurochem. Int. 2008 52 675 682 10.1016/j.neuint.2007.08.014 17928106
52. Zhou X Differential modulation of GABAA and NMDA receptors by an alpha7-nicotinic acetylcholine receptor agonist in chronic glaucoma Front. Mol. Neurosci. 2017 10 422 10.3389/fnmol.2017.00422 29326549
53. Najm R Jones EA Huang Y Apolipoprotein E4, inhibitory network dysfunction, and Alzheimer’s disease Mol. Neurodegener. 2019 14 24 10.1186/s13024-019-0324-6 31186040
54. Lunghi C Emir UE Morrone MC Bridge H Short-term monocular deprivation alters GABA in the adult human visual cortex Curr. Biol. 2015 25 1496 1501 10.1016/j.cub.2015.04.021 26004760
55. Loerch PM Evolution of the aging brain transcriptome and synaptic regulation PLoS ONE 2008 3 e3329 10.1371/journal.pone.0003329 18830410
56. Stanley DP Shetty AK Aging in the rat hippocampus is associated with widespread reductions in the number of glutamate decarboxylase-67 positive interneurons but not interneuron degeneration J. Neurochem. 2004 89 204 216 10.1111/j.1471-4159.2004.02318.x 15030405
57. Caspary DM Milbrandt JC Helfert RH Central auditory aging: GABA changes in the inferior colliculus Exp. Gerontol. 1995 30 349 360 10.1016/0531-5565(94)00052-5 7556513
58. Hua T Kao C Sun Q Li X Zhou Y Decreased proportion of GABA neurons accompanies age-related degradation of neuronal function in cat striate cortex Brain Res. Bull. 2008 75 119 125 10.1016/j.brainresbull.2007.08.001 18158105
59. Stagg CJ Bachtiar V Johansen-Berg H What are we measuring with GABA magnetic resonance spectroscopy? Commun. Integr. Biol. 2011 4 573 575 10.4161/cib.16213 22046466
60. Hendry SH Fuchs J deBlas AL Jones EG Distribution and plasticity of immunocytochemically localized GABAA receptors in adult monkey visual cortex J. Neurosci. 1990 10 2438 2450 10.1523/JNEUROSCI.10-07-02438.1990 2165524
61. Gu M GABA editing with macromolecule suppression using an improved MEGA-SPECIAL sequence Magn. Reson. Med. 2018 79 41 47 10.1002/mrm.26691 28370458
62. Wu GK Arbuckle R Liu BH Tao HW Zhang LI Lateral sharpening of cortical frequency tuning by approximately balanced inhibition Neuron 2008 58 132 143 10.1016/j.neuron.2008.01.035 18400169
63. Leventhal AG Wang Y Pu M Zhou Y Ma Y GABA and its agonists improved visual cortical function in senescent monkeys Science 2003 300 812 815 10.1126/science.1082874 12730605
64. Cook E Hammett ST Larsson J GABA predicts visual intelligence Neurosci. Lett. 2016 632 50 54 10.1016/j.neulet.2016.07.053 27495012
65. Stagg CJ Bachtiar V Johansen-Berg H The role of GABA in human motor learning Curr. Biol. 2011 21 480 484 10.1016/j.cub.2011.01.069 21376596
66. Boroojerdi B Battaglia F Muellbacher W Cohen LG Mechanisms underlying rapid experience-dependent plasticity in the human visual cortex Proc. Natl Acad. Sci. USA 2001 98 14698 14701 10.1073/pnas.251357198 11734655
67. Park S Shorter sleep duration is associated with lower GABA levels in the anterior cingulate cortex Sleep. Med. 2020 71 1 7 10.1016/j.sleep.2020.02.018 32447224
68. Qiu M Ramulu PY Boland MV Association between sleep parameters and glaucoma in the United States population: national health and nutrition examination survey J. Glaucoma 2019 28 97 104 10.1097/IJG.0000000000001169 30585942
69. Wang H Zhang Y Ding J Wang N Changes in the circadian rhythm in patients with primary glaucoma PLoS ONE 2013 8 e62841 10.1371/journal.pone.0062841 23658653
70. Carp J Park J Polk TA Park DC Age differences in neural distinctiveness revealed by multi-voxel pattern analysis Neuroimage 2011 56 736 743 10.1016/j.neuroimage.2010.04.267 20451629
71. Bowman CR Chamberlain JD Dennis NA Sensory representations supporting memory specificity: age effects on behavioral and neural discriminability J. Neurosci. 2019 39 2265 2275 10.1523/JNEUROSCI.2022-18.2019 30655350
72. Dudek FE Sutula TP Epileptogenesis in the dentate gyrus: a critical perspective Prog. Brain Res. 2007 163 755 773 10.1016/S0079-6123(07)63041-6 17765749
73. Atallah BV Scanziani M Instantaneous modulation of gamma oscillation frequency by balancing excitation with inhibition Neuron 2009 62 566 577 10.1016/j.neuron.2009.04.027 19477157
74. Haider B Duque A Hasenstaub AR McCormick DA Neocortical network activity in vivo is generated through a dynamic balance of excitation and inhibition J. Neurosci. 2006 26 4535 4545 10.1523/JNEUROSCI.5297-05.2006 16641233
75. Mathews GC Diamond JS Neuronal glutamate uptake Contributes to GABA synthesis and inhibitory synaptic strength J. Neurosci. 2003 23 2040 2048 10.1523/JNEUROSCI.23-06-02040.2003 12657662
76. Findley CA Bartke A Hascup KN Hascup ER Amyloid beta-related alterations to glutamate signaling dynamics during Alzheimer’s disease progression ASN Neuro 2019 11 1759091419855541 10.1177/1759091419855541 31213067
77. Siskova Z Dendritic structural degeneration is functionally linked to cellular hyperexcitability in a mouse model of Alzheimer’s disease Neuron 2014 84 1023 1033 10.1016/j.neuron.2014.10.024 25456500
78. Hascup KN Hascup ER Altered neurotransmission prior to cognitive decline in AbetaPP/PS1 mice, a model of Alzheimer’s disease J. Alzheimers Dis. 2015 44 771 776 10.3233/JAD-142160 25374106
79. Minkeviciene R Age-related decrease in stimulated glutamate release and vesicular glutamate transporters in APP/PS1 transgenic and wild-type mice J. Neurochem. 2008 105 584 594 10.1111/j.1471-4159.2007.05147.x 18042177
80. Rupsingh R Borrie M Smith M Wells JL Bartha R Reduced hippocampal glutamate in Alzheimer disease Neurobiol. Aging 2011 32 802 810 10.1016/j.neurobiolaging.2009.05.002 19501936
81. Kantarci K Proton MR spectroscopy in mild cognitive impairment and Alzheimer disease: comparison of 1.5 and 3 T AJNR Am. J. Neuroradiol. 2003 24 843 849 12748083
82. Muthukumaraswamy SD Edden RA Jones DK Swettenham JB Singh KD Resting GABA concentration predicts peak gamma frequency and fMRI amplitude in response to visual stimulation in humans Proc. Natl Acad. Sci. USA 2009 106 8356 8361 10.1073/pnas.0900728106 19416820
83. Violante IR GABA deficit in the visual cortex of patients with neurofibromatosis type 1: genotype-phenotype correlations and functional impact Brain 2013 136 918 925 10.1093/brain/aws368 23404336
84. Bridge H Altered neurochemical coupling in the occipital cortex in migraine with visual aura Cephalalgia 2015 35 1025 1030 10.1177/0333102414566860 25631169
85. Bednarik P Neurochemical and BOLD responses during neuronal activation measured in the human visual cortex at 7 Tesla J. Cereb. Blood Flow. Metab. 2015 35 601 610 10.1038/jcbfm.2014.233 25564236
86. Ip IB Combined fMRI-MRS acquires simultaneous glutamate and BOLD-fMRI signals in the human brain Neuroimage 2017 155 113 119 10.1016/j.neuroimage.2017.04.030 28433623
87. Schallmo MP Glutamatergic facilitation of neural responses in MT enhances motion perception in humans Neuroimage 2019 184 925 931 10.1016/j.neuroimage.2018.10.001 30312807
88. Kiemes A GABA, glutamate and neural activity: a systematic review with meta-analysis of multimodal (1)H-MRS-fMRI studies Front. Psychiatry 2021 12 644315 10.3389/fpsyt.2021.644315 33762983
89. Donahue MJ Near J Blicher JU Jezzard P Baseline GABA concentration and fMRI response Neuroimage 2010 53 392 398 10.1016/j.neuroimage.2010.07.017 20633664
90. Muthukumaraswamy SD Evans CJ Edden RA Wise RG Singh KD Individual variability in the shape and amplitude of the BOLD-HRF correlates with endogenous GABAergic inhibition Hum. Brain Mapp. 2012 33 455 465 10.1002/hbm.21223 21416560
91. Enzi B Glutamate modulates resting state activity in the perigenual anterior cingulate cortex - a combined fMRI-MRS study Neuroscience 2012 227 102 109 10.1016/j.neuroscience.2012.09.039 23022216
92. Duncan NW Negative childhood experiences alter a prefrontal-insular-motor cortical network in healthy adults: A preliminary multimodal rsfMRI-fMRI-MRS-dMRI study Hum. Brain Mapp. 2015 36 4622 4637 10.1002/hbm.22941 26287448
93. Duncan NW Enzi B Wiebking C Northoff G Involvement of glutamate in rest-stimulus interaction between perigenual and supragenual anterior cingulate cortex: a combined fMRI-MRS study Hum. Brain Mapp. 2011 32 2172 2182 10.1002/hbm.21179 21305662
94. Cadena EJ A longitudinal multimodal neuroimaging study to examine relationships between resting state glutamate and task related BOLD response in schizophrenia Front Psychiatry 2018 9 632 10.3389/fpsyt.2018.00632 30555359
95. Dreyer EB Zurakowski D Schumer RA Podos SM Lipton SA Elevated glutamate levels in the vitreous body of humans and monkeys with glaucoma Arch. Ophthalmol. 1996 114 299 305 10.1001/archopht.1996.01100130295012 8600890
96. Brooks DE Garcia GA Dreyer EB Zurakowski D Franco-Bourland RE Vitreous body glutamate concentration in dogs with glaucoma Am. J. Vet. Res. 1997 58 864 867 9256971
97. Carter-Dawson L Vitreal glutamate concentration in monkeys with experimental glaucoma Invest. Ophthalmol. Vis. Sci. 2002 43 2633 2637 12147596
98. Honkanen RA Vitreous amino acid concentrations in patients with glaucoma undergoing vitrectomy Arch. Ophthalmol. 2003 121 183 188 10.1001/archopht.121.2.183 12583783
99. Gramlich OW Godwin CR Wadkins D Elwood BW Kuehn MH Early functional impairment in experimental glaucoma is accompanied by disruption of the GABAergic system and inceptive neuroinflammation Int. J. Mol. Sci. 2021 22 7581 10.3390/ijms22147581 34299211
100. Zhou X Zhang R Zhang S Wu J Sun X Activation of 5-HT1A receptors promotes retinal ganglion cell function by inhibiting the cAMP-PKA pathway to modulate presynaptic GABA release in chronic glaucoma J. Neurosci. 2019 39 1484 1504 10.1523/JNEUROSCI.1685-18.2018 30541912
101. Zhou X Alpha7 nicotinic acetylcholine receptor agonist promotes retinal ganglion cell function via modulating GABAergic presynaptic activity in a chronic glaucomatous model Sci. Rep. 2017 7 1734 10.1038/s41598-017-02092-6 28496108
102. Cello KE Nelson-Quigg JM Johnson CA Frequency doubling technology perimetry for detection of glaucomatous visual field loss Am. J. Ophthalmol. 2000 129 314 322 10.1016/S0002-9394(99)00414-6 10704546
103. Burgansky-Eliash Z Optical coherence tomography machine learning classifiers for glaucoma detection: a preliminary study Invest. Ophthalmol. Vis. Sci. 2005 46 4147 4152 10.1167/iovs.05-0366 16249492
104. Mohammadzadeh V Longitudinal macular structure-function relationships in glaucoma Ophthalmology 2020 127 888 900 10.1016/j.ophtha.2020.01.023 32173112
105. Nouri-Mahdavi K Prediction of visual field progression from OCT structural measures in moderate to advanced glaucoma Am. J. Ophthalmol. 2021 226 172 181 10.1016/j.ajo.2021.01.023 33529590
106. Provencher SW Automatic quantitation of localized in vivo 1H spectra with LCModel NMR Biomed. 2001 14 260 264 10.1002/nbm.698 11410943
107. Provencher SW Estimation of metabolite concentrations from localized in vivo proton NMR spectra Magn. Reson. Med. 1993 30 672 679 10.1002/mrm.1910300604 8139448
108. Trivedi V Widespread brain reorganization perturbs visuomotor coordination in early glaucoma Sci. Rep. 2019 9 14168 10.1038/s41598-019-50793-x 31578409
109. Brainard DH The psychophysics toolbox Spat. Vis. 1997 10 433 436 10.1163/156856897X00357 9176952
110. Bang, J. W. Data sets for GABA reduction in the visual cortex of glaucoma patients. OSF 10.17605/OSF.IO/85CDS (2022)