# GABAergic inhibition in the human visual cortex relates to eye dominance

I. Betina Ip, Uzay Emir, Claudia Lunghi, Andrew Parker, Holly Bridge

► To cite this version:

I. Betina Ip, Uzay Emir, Claudia Lunghi, Andrew Parker, Holly Bridge. GABAergic inhibition in the human visual cortex relates to eye dominance. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-95685-1. doi: 10.1038/s41598-021-95685-1

HAL Id: hal-03357198

https://hal.science/hal-03357198

Submitted on 28 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# GABAergic inhibition in the human visual cortex relates to eye dominance

I. Betina Ip<fnref n="1" />,²✉, Uzay E. Emir³, Claudia Lunghi⁴, Andrew J. Parker²,⁵ & Holly Bridge<fnref n="1" />

Binocular vision is created by fusing the separate inputs arriving from the left and right eyes. ‘Eye dominance’ provides a measure of the perceptual dominance of one eye over the other. Theoretical models suggest that eye dominance is related to reciprocal inhibition between monocular units in the primary visual cortex, the first location where the binocular input is combined. As the specific inhibitory interactions in the binocular visual system critically depend on the presence of visual input, we sought to test the role of inhibition by measuring the inhibitory neurotransmitter GABA during monocular visual stimulation of the dominant and the non-dominant eye. GABA levels were measured in a single volume of interest in the early visual cortex, including V1 from both hemispheres, using a combined functional magnetic resonance imaging and magnetic resonance spectroscopy (combined fMRI-MRS) sequence on a 7-Tesla MRI scanner. Individuals with stronger eye dominance had a greater difference in GABAergic inhibition between the eyes. This relationship was present only when the visual system was actively processing sensory input and was not present at rest. We provide the first evidence that imbalances in GABA levels during ongoing sensory processing are related to eye dominance in the human visual cortex. Our finding supports the view that intracortical inhibition underlies normal eye dominance.

Eye dominance in the healthy visual system is the preference for using one eye over the other in a visual task<fnref n="1" />. Extreme eye dominance is associated with amblyopia, a neurodevelopmental disorder that causes a chronic loss of normal monocular and binocular function with an incidence of 2–4% in the general population²⁻⁴. Input from the two eyes is combined for the first time in the primary visual cortex to serve binocular vision. Therefore, processing at this stage is thought to be decisive in eye preference. Understanding the relationship between eye dominance and the brain is thus an opportunity to investigate the neural mechanisms used by the cerebral cortex to serve perception.

Evidence for a neural mechanism of eye dominance comes from study of the abnormal binocular visual system. In adults, double vision arises when input from the two eyes does not correspond and this is associated with severe discomfort and headaches. In contrast, children who grow up with one deviating eye often develop strabismic amblyopia but perceive only a single image, not two. This is because their visual system has prevented double images from reaching perception by making the input from the deviating eye non-visible and relaying the image from the non-deviating eye. Known as binocular suppression, this type of chronic suppression could be ‘nature’s way out of trouble’⁵. The drawback of suppression is that it is related to a loss of binocular visual function⁶, most notably for stereopsis and binocular summation. Hence, visual cortex suppression has been named as one of the prime causes of the perceptual deficits observed in amblyopia⁷.

Severely imbalanced vision causes both anatomical and functional abnormalities in the primary visual cortex⁸, with functional abnormalities likely maintained by GABAergic inhibition⁹. The strongest evidence in support of this view has come from pharmacological studies in amblyopic animal models. In strabismic amblyopic cats, localized application of the GABA<sub>A</sub> receptor antagonist bicuculline to V1 reduced binocular suppression¹⁰. Similarly, in amblyopic rats, pharmacologically decreasing GABAergic signalling in the visual cortex has been linked to a recovery of normal structure and visual function in the amblyopic eye¹¹,¹². This recovery was blocked by

<footnote n="1">¹Wellcome Centre for Integrative Neuroimaging, FMRIB, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford OX3 9DU, UK. ²Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford OX1 3PT, UK. ³School of Health Sciences, Purdue University, West Lafayette, Indiana 47907, USA. ⁴Laboratoire des systèmes perceptifs, Département d’études cognitives, École normale supérieure, PSL University, CNRS, 75005 Paris, France. ⁵Institut für Biologie, Otto-Von-Guericke Universität, 39120 Magdeburg, Germany. ✉email: betina.ip@ndcn.ox.ac.uk</footnote>
a. Binocular Rivalry

<figure><img src="image_3.png" /></figure>

b. Group Phase Durations

<figure><img src="image_4.png" /></figure>

c. Phase Duration
DE NDE
[%]

Eye dominance

**Figure 1.** (a) Sensory eye dominance was measured in a binocular rivalry experiment outside of the scanner. Stimuli were presented on a mirror stereoscope and participants reported their subjective percept using a continuous key press. (b) Histograms of phase distributions for the dominant (DE, green) and non-dominant (NDE, blue) eye were pooled across participants. Data were normalised by the median dominant eye phase duration. (c) Average phase durations in seconds for DE (green) and NDE (blue) are plotted on the left. On the right are the biases in phase duration quantified for each participant by calculating the difference between the dominant and non-dominant eye phase duration as a percentage of the non-dominant eye phase duration. This index is henceforth referred to as ‘Eye Dominance [%]’ (clear bar, EDI).

cortical delivery of diazepam, a GABA<sub>A</sub> agonist<sup>11</sup>. These studies provide strong support for a role of GABAergic
inhibition in maintaining extreme eye dominance in the adult brain.

Amblyopic suppression may be an extreme form of normal binocular interaction, revealed in a subtle form in normal participants using binocular rivalry<sup>9</sup>. Binocular rivalry is a widely used method to quantify sensory eye dominance in the normal-sighted population<sup>1</sup>. It induces continuous perceptual alternations between incongruent images presented separately to the left or right eye<sup>13,14</sup>. Theoretical models suggest that intracortical inhibition could explain the pattern of perceptual dynamics<sup>15,16</sup>. In recent years, support for this view has come from studies applying magnetic resonance spectroscopy (MRS) to measure GABAergic inhibition in the human brain in the absence of task specific stimulation ('rest'). Greater levels of resting GABA levels in the early visual cortex (EVC) are related to stronger perceptual suppression<sup>17</sup>, and slower, less frequent perceptual switches<sup>18,19</sup>. Pharmacological manipulations of GABA have shown a causal relationship between GABA and binocular rivalry dynamics: increasing GABA reduced perceptual switches and increased percept durations<sup>18</sup> and perceptual suppression<sup>18,20</sup>. In addition, increases in eye dominance duration during binocular rivalry induced by short-term monocular deprivation have been correlated with lower resting GABA levels in the primary visual cortex<sup>21</sup>. These results have highlighted the role of baseline GABAergic inhibition in shaping the competitive interactions in the visual brain, yet the hypothesis that imbalances during active stimulation of each eye are related to eye dominance awaits critical support.

To test the long-standing prediction that imbalances in cortical inhibition during monocular visual response are related to eye dominance⁹, we set out to measure GABA during functional stimulation. Specifically, we wanted to quantify the ‘interocular’ difference in GABA during monocular stimulation of either the stronger or the weaker eye. Conventionally, MRS experiments measure GABA levels during the resting state, when participants have their eyes closed, or simply watch a movie with both eyes open. In absence of a specific task, such measurements at rest are thought to reflect the overall inhibitory tone in the brain area under study²². In contrast, our study takes advantage of recent advances in ultra-high field functional MRS methodology, that allow measurement of metabolite levels during task performance²³⁻²⁶. For all participants, we determined eye dominance both using binocular rivalry and the magnitude of interocular inhibition, manifest in the difference in GABA levels during stimulation of the stronger and weaker eyes respectively. Our results provide the first evidence that eye dominance in the human visual cortex is related to neural mechanisms of active interocular inhibition.

Results

**Binocular rivalry quantifies eye dominance.** For each participant, sensory eye dominance was measured individually using binocular rivalry stimuli presented by means of a Wheatstone stereoscope (Fig. 1a). Binocular rivalry phase distributions pooled across all observers were not normally distributed (Shapiro-Wilk Normality Test, W-stats = 0.64, *p* = 0.0001) as demonstrated by the typical skewedness²⁷ of bistable percept phase durations (Fig. 1b). The median phase durations in seconds were used as a measure of central tendency²⁸. Median phase durations across 12 observers were normally distributed (W-stats = 0.95, *p* = 0.337), hence parametric statistics were used for significance testing with an uncorrected alpha of 0.05. After confirming that there was no overall difference in right vs left eye median phase durations (two-tailed paired samples *t*-test, *t*(11) = 0.228, *p* = 0.824), eye dominance was assigned based on the length of the median phase duration (Fig. 1c, left; 5 observers = left eye dominant, 7 observers = right eye dominant, left vs right median). Not surprisingly, median phase durations were greater for the DE compared to the NDE (two-tailed paired samples *t*-test, *t*(11) = 2.960, *p* = 0.013). To obtain a singular measure of the increase in eye dominance, an eye dominance index (Fig. 1c right;
a

Visual Stimulation fMRS

<figure><img src="image_5.png" /></figure>

b

<figure><img src="image_6.png" /></figure>

**Figure 2.** The combined fMRI-MRS monocular visual stimulation paradigm. (a) 7-Tesla combined fMRI-MRS data were collected from the early visual cortex during viewing of visual stimuli through the dominant or non-dominant eye. Visual input comprised of a constant fixation task, while a 50% contrast flashing checkerboard was cycled on ('checkerboard') and off ('fixation') in the background. (b) A representative MRS spectrum during viewing of visual stimuli through one eye, created by averaging across checkerboard and fixation periods. From bottom to top: blue spectrum shows LCModel fit to raw data. Black spectrum shows raw data. Immediately above are the residuals of the LCModel fit. Topmost line shows the LCModel fit to the GABA signal. For illustrative purposes, the GABA signal was multiplied by 10.

'EDI', see "Methods" section, Eq. (1)) was calculated. The EDI measures the increase in duration of the DE phase as a percentage of the NDE phase (range 0.64–23.1%, mean: 7.9 ± 8.2%).

**Monocular visual stimulation protocol during combined fMRI-MRS acquisition.** To measure cortical inhibition during monocular visual stimulation in the visual cortex (Fig. 2a), a 2×2×2 cm³ MRS volume-of-interest (VOI) was symmetrically placed along the midline of the calcarine sulcus to cover the early visual cortex (EVC) using a combined fMRI-MRS sequence (see 'Supplementary Information'). We obtained high quality MRS spectra from the EVC volume (Fig. 2b), including detection of inhibitory neurotransmitter GABA. Across participants, an average of 85 ± 13.7% of the MRS VOI overlapped with the fMRI activation map (checkerboard > fixation), confirming that measurement of the neurochemical levels was measured from visually-stimulated tissue. No difference in percentage overlap was found between monocular stimulation conditions (dominant eye: 85 ± 13.7%, non-dominant eye: 86 ± 11.3%, Wilcoxon Rank Sum Test, *p* = 0.72).

**Interocular GABA from early visual cortex correlates with eye dominance.** Mutual inhibition between neural inputs from each eye to the primary visual cortex is thought to drive eye dominance. To reveal differences in inhibitory responses between the dominant (DE) and non-dominant eye (NDE), participants alternately viewed visual stimuli with either their DE or NDE while the non-stimulated eye was covered with a semi-transparent occluder. Across participants, dominant eye GABA tended to be higher compared to NDE viewing, however the difference was not statistically significant (GABA:H₂O, paired *t*-test, *t*(11) = 2.13, *p* = 0.06; GABA:tCr, paired *t*-test, *t*(11) = 1.72, *p* = 0.11). The difference in GABA levels between eyes was then quantified as a single metric, the 'interocular' GABA (Fig. 3a). We related the behaviourally derived eye dominance index (EDI), measured in the binocular rivalry experiment, to interocular GABA. EDI was significantly correlated with interocular GABA when using GABA scaled to the internal reference signal of unsuppressed water (Fig. 3b, GABA:H₂O, Spearman's Rank Correlation *r* = 0.62, uncorrected *p* = 0.037). To control for the possibility that the metabolite reference method influenced the results, we also scaled GABA to the summed signal of creatine and phosphocreatine ('GABA:tCr'). The correlation remained significant (GABA:tCr, *r* = 0.59, uncorrected *p* = 0.049). In an exploratory analysis of this relationship, we show that there is no correlation between EDI and dominant eye GABA (Fig. 3c: Spearman's Rank Correlation, GABA:H₂O, *r* = -0.11, *p* = 0.73; GABA:tCr, *r* = -0.06, *p* = 0.83), but a negative association between EDI with non-dominant eye GABA (Fig. 3d: Spearman's Rank Correlation, GABA:H₂O, *r* = -0.59, *p* = 0.048; GABA:tCr, *r* = -0.61, *p* = 0.037). We conclude that stronger eye dominance corresponds to a greater difference in GABAergic inhibition during dominant compared to non-dominant eye stimulation. Investigating the correlation within conditions suggested that the effect is partly accounted for by a failure of the NDE to inhibit the DE, as stronger eye dominance was associated with lower GABA during stimulation of the weaker eye.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 3.** (a) The interocular GABA metric was calculated to show the difference in GABAergic inhibition between activation of the DE or NDE eye. (b) Eye dominance was highly correlated with interocular GABA. Further exploratory analysis showed that (c) eye dominance was not correlated with DE GABA:H₂O, but that it was correlated with NDE GABA:H₂O (d). Solid line represents best fitting linear model, hatched lines represent 95% confidence interval of the regression line. Bold typed **r**= Spearman's Rank Correlation Coefficient for GABA scaled to water, ‘r’ Spearman's Rank Correlation Coefficient for GABA scaled to total Creatine, \*uncorrected *p* < 0.05.</figcaption></figure>

Resting MRS data was also acquired while participants kept their eyes closed and no stimuli were delivered. We found no relationship between EDI and resting GABA (GABA:H₂O, *r* = -0.14, *p* = 0.67; GABA:tCr, *r* = -0.09, *p* = 0.77). This control analysis shows that monocular visual stimulation was necessary to reveal inhibitory interactions relevant to eye dominance in the early visual cortex. Exploratory correlations between other metabolites and EDI are reported in the 'Supplementary Information' section in Table S1. We also tested whether the interocular difference in GABA levels (ΔGABA) was related to eye dominance. ΔGABA was calculated by comparing checkerboard > fixation prior to calculating the interocular ΔGABA metric. We found no correlation between eye dominance and interocular ΔGABA (*r* = 0.24, *p* = 0.44) suggesting that the addition of a checkerboard to fixation did not reveal a relationship to eye dominance.

**Eye dominance was not correlated with interocular BOLD-signal to flashing checkerboards.** The combined fMRI-MRS sequence measures fMRI and MRS data in the same TR. The next step was to investigate whether the BOLD-fMRI response from the same region in the early visual cortex also related to eye dominance. The complementary measures reflect different aspects of binocular function and require different analysis approaches. The fMRS analysis revealed how inhibitory signals differed between DE and NDE irrespective of checkerboard or fixation condition. On the other hand, the fMRI-signal reveals differences in BOLD signal (checkerboard > fixation) between DE and NDE (Fig. 4a). The fMRI depends on increased activation to a flashing checkerboard relative to fixation, before it is compared across DE and NDE eyes, and therefore cannot provide any information about activation to the fixation task alone. We found no difference in the checkerboard-evoked BOLD-signal between DE and NDE (Fig. 4b: Wilcoxon's Rank Sum, *z* = -0.77; *p* = 0.43). We then calculated the interocular BOLD-signal (for 'interocular BOLD', see "Methods" section Eq. (3)). There was no correlation between EDI and interocular BOLD (Fig. 4c, Spearman's *r* = 0.251, *p* = 0.430).
a

<figure><img src="image_11.png" /></figure>

b

<figure><img src="image_12.png" /></figure>

c

<figure><img src="image_13.png" /></figure>

**Figure 4.** (a) %BOLD-change was obtained by comparing checkerboard > fixation prior to calculating the interocular BOLD metric. (b) %BOLD-signal change to a flashing checkerboard did not differ between dominant and non-dominant eye. (c) There was no relationship between eye dominance and the interocular %BOLD to the flashing checkerboard. Errors are 1 ± SEM. *r* Spearman's Rank Correlation. Solid line represents best fitting linear model, hatched lines represent 95% confidence interval of the regression line.

a

Perceptual suppression =

<figure><img src="image_14.png" /></figure>

b

<figure><img src="image_15.png" /></figure>

**Figure 5.** (a) To compare with previous data, we calculated an alternative measure of binocular rivalry dynamics, the 'perceptual suppression' rate (Robertson et al.¹⁷). (b) Resting GABA:tCr levels was correlated with binocular rivalry perceptual suppression [dominant eye/(dominant eye + mixed percept)]. *r* = Pearson's linear correlation co-efficient. *p* < 0.05, uncorrected. Solid line represents best fitting linear model, hatched lines represent 95% confidence interval of the regression line.

**Resting GABA relates to binocular rivalry suppression.** Our dataset provided an opportunity to replicate a prior finding on binocular dynamics in the early visual cortex, obtained using different methods. Robertson et al., found a relationship between GABAergic inhibition and the proportion of perceptual suppression¹⁷. In their study, participants viewed stimuli with both eyes and performed a simple task at fixation, a state referred to as 'resting MRS'. We used our resting GABA, when participants had their eyes closed and no stimulation was delivered, to attempt to replicate their findings. Using the same data from which we derived the eye dominance metric, we calculated binocular rivalry suppression, the proportion of the time seeing the dominant percept, divided by the sum of the dominant and mixed percept proportions (Fig. 5a). There was a significant correlation between perceptual suppression and resting GABA:tCr (Fig. 5b, *r* = 0.66, *p* = 0.02) although the correlation was not significant for GABA:H₂O (*r* = 0.48, *p* = 0.12). This result suggests that the link between GABAergic inhibition and perceptual suppression can be replicated at ultra-high field strength, with the combined fMRI-MRS sequences and a different definition of 'resting MRS'. The finding that GABA:tCr was correlated with binocular suppression while GABA:H₂O was not, suggests that GABA levels may not be wholly independent from influences of Creatine energy metabolism²⁹.
Discussion

Eye dominance in the normal visual system may be regarded as a window into how the brain selects and com-
bines information from the two eyes. Extreme eye dominance has a pathological form manifest in 'amblyopia',
which impairs normal visual function. We demonstrate a link between intracortical inhibition and normal eye
dominance in the healthy visual system, based on a novel approach to measurement of inhibition in the cortex.
We applied monocular visual stimulation during combined fMRI-MRS at 7-Tesla. There was no link between
eye dominance and inhibition in the simultaneously acquired BOLD-signal, or in resting GABA levels.

We show that GABA-levels measured during visual stimulation in the early visual cortex are linked to eye dominance in the healthy human brain. Our results are in agreement with the view that the neural mechanism of eye dominance is mediated by GABAergic inhibition<sup>9,30</sup>. With regard to the difference in GABA between DE and NDE eye ('interocular GABA'), individuals with stronger eye dominance showed a greater difference in GABA. In an exploratory analysis, we further show that this difference is partially driven by a failure to inhibit the stronger eye during activation of the weaker eye. Eye dominance may therefore be due to a systematically weaker ability of the non-dominant eye to inhibit the stronger eye during active viewing. In confirmation, no relationship to eye dominance was present when GABA levels were measured during the resting state, during which participants had their eyes closed. Our control analyses ('Supplementary Information', Table S2) showed that spectral quality measures from our data (SNR, tCr line width), did not correlate with interocular GABA. In addition, interocular GABA related to eye dominance independently of the reference method. These results suggest that our main result was not driven by differences in spectral quality or the metabolite reference method.

What might the visually driven GABA signals represent? While the coarse spatial and temporal resolution of single-voxel MRS precludes a definite assignment of MRS-visible GABA to a particular cellular function or compartment, our functional paradigm associated GABA with neural activity during monocular visual stimulation. Monocular stimulation would have increased local metabolic demand and neurotransmitter release across neurons involved in visual processing and task-performance. Our analysis shows that the majority of the MRS VOI overlaps with functional activation and primary visual cortex (V1). While around 15-20% of the neurons in the cerebral cortex are inhibitory, a high percentage (~70-80%) of GABAergic neurons in V1 are parvalbumin expressing (PV +) interneurons<sup>31,32</sup>. The network of PV + interneurons in V1 has been extensively linked to early ocular dominance plasticity<sup>33</sup>. Thus, visually driven GABA signals observed in our study could represent intra-cortical inhibition from GABAergic interneurons that bias perception in favour of the dominant eye. Our findings here emphasize that revealing of inhibitory interactions in the binocular visual system depends on the status of visual input.

There was no relationship between eye dominance and the interocular BOLD-signal. It is important to point
out that the BOLD-signal was dependent the contrast between checkerboard and fixation, whereas the fMRS
GABA measures were calculated across both stimulus conditions. Hence, the two measures cannot be directly
compared. Only a handful of studies have investigated the relationship between eye dominance and visual cortex
BOLD-signal change and it is fair to say that the results have so far been inconclusive, partly due to differences
in how eye dominance was assigned. Greater BOLD-activation in V1 has been found during DE versus NDE
stimulation<sup>34,35</sup>, whereas others find no difference<sup>36</sup>. In these studies, eye dominance was assigned using other
approaches, such as 'sighting dominance'. Different methods of assigning eye dominance often do not agree
with one another (for a review see<sup>37</sup>). We measured sighting dominance for all our subjects and found that it
corresponded with assignments of sensory eye dominance assignments in only 7 out of 12 participants. A better
comparison with our study may be found from a recent study using binocular rivalry: here they reported that
increases in eye dominance induced by monocular deprivation are linked to increases in deprived eye BOLD-
activation<sup>38</sup>. The lack of any correlation in our BOLD-signal data is consistent with the view that the BOLD-
signal is not sensitive to interocular inhibition in early visual cortex<sup>39</sup>, whereas MRS has significant sensitivity
to this signal<sup>17,18,21</sup>.

While we measured eye dominance using binocular rivalry in participants with normal stereo vision, a signifi-
cant minority of the general population do not have normal binocular functions. Using a meta-analysis, Chopin
et al.<sup>40</sup> estimated that ~7% of the general population fall into this category, which can include those with binocular
vision pathologies. More specifically, stereo-anomalous can be defined as individuals who fail stereo-vision tests,
but unlike amblyopes, they do not exhibit a difference in visual acuity between eyes. In such individuals, mean
binocular rivalry eye dominance durations were shown to be twice as long as those of stereo-normals<sup>41</sup>. It is
well established that amblyopes have abnormal eye dominance<sup>42</sup> that can approximate normal temporal patterns
when vision in the fellow eye is attenuated using a neutral density filter<sup>43</sup>. More recently, sensory eye dominance
has been used as an index for enhancement of binocular vision in amblyopes, with eye dominance modulated
by visual perceptual training regimes<sup>44,45</sup>. Our result showing that eye dominance in the normal observer is cor-
related with GABAergic signalling could provide a basis for understanding the vision of individuals with greater
asymmetry between eyes such as stereo-anomalous individuals and amblyopes.

Our study specifically set out to evaluate differences in interocular inhibition. However, it included a rela-
tively small cohort (N = 12) and statistical results were not corrected for multiple comparisons. Nevertheless,
our methods successfully replicated a recent result linking an alternate measure of interocular suppression and
resting GABA levels, despite using a smaller sample size than the original study¹⁷. Additionally, our MRS-voxel
size did not permit the spatial resolution to target the primary visual cortex alone. While a large percentage of
the MRS voxel in early visual cortex co-localised to the primary visual cortex (47%), portions of V2 and V3 were
also included. We therefore cannot rule out a contribution from visual areas beyond V1. Whilst we are confi-
dent that our protocol dissociated dominant from non-dominant eye stimulation at the input level, a long-term
improvement in the spatial resolution of MRS imaging is required to resolve these signals at the circuit level in
ocular dominance columns.
## Conclusion

In conclusion, these results have shown that GABAergic responses in the healthy human visual cortex relate to eye dominance. This relationship is specific to the conditions of visual stimulation and not present during GABA measured during rest. Although our study was conducted using a relatively small cohort of participants and should be interpreted with caution, we provide compelling evidence supporting a role of GABAergic inhibition in sensory eye dominance⁹,¹⁵,¹⁶. The extent to which these relationships can be related to clinical conditions of extreme eye dominance remains to be determined.

## Methods

**Participants.** Thirteen participants (7 females, 29.2 ± 6.0 years, age range 21–42 years) including two of the authors took part in the study. All data presented were collected as part of a behavioural and MRI-data set, of which a subset has been published²⁵. One participant was identified as an outlier due to a large percentage (68.5%) of mixed periods in the binocular rivalry experiment. Mixed periods are time points in which the participant reports seeing a piecemeal version of left and right eyes' percept. Participants with a high percentage of mixed periods are excluded because they would not contribute enough data for analysis, may have an undiagnosed visual condition or an error in performing the task. Participants had normal or corrected-to-normal vision, no neurological impairments and normal stereo-acuity (< 120 arc s, TNO Stereo test, Lameris, Utrecht). All subjects were involved in a 1-h psychophysical session to measure binocular and monocular visual function, as well as eye dominance, and took part in a 1.5-h MRI session to measure interleaved changes in neurochemical levels and BOLD-activity. All subjects gave written informed consent. Approval for the study was obtained by the University of Oxford Research Ethics Committee (MSD-IDREC-C1-2014-146). The study was carried out in accordance to the Declaration of Helsinki.

**Behavioural protocol and procedure.** Participants' eye dominance was measured in a separate psychophysical session outside of the scanner (Fig. 1a). Sighting eye dominance was measured using the Miles Test⁴⁶, which identifies the dominant eye as the one used for sighting of a distant target when viewed through an aperture. Sensory eye dominance was measured using binocular rivalry. While both sighting and sensory eye dominance tests were collected, only sensory eye dominance was used in subsequent analysis to provide a quantitative measure of eye dominance. Head position was stabilised with a chin and headrest. Stimuli were displayed on two gamma-linearised CRT monitors (viewing distance 57 cm) viewed through a Wheatstone mirror stereoscope. Stimuli were generated using MATLAB (The MathWorks, Natick, MA) with Psychophysics Toolbox⁴⁷ running on an Apple Mac-mini. Binocular rivalry stimuli were two achromatic gratings (orientation: ± 45 deg off vertical, spatial frequency: 6 cycle/deg, contrast: 100%, diameter: 3.2 deg) presented on a uniform mid-grey background in the centre of vision. After successfully fusing a Nonius fusion target, participants self-initiated the trial and reported the perceived orientation of the tilted grating using the left (counter-clockwise) and right (clockwise) arrows of a computer keyboard; the upward arrow was pressed when a mixed percept was perceived. The orientations of the members of the grating pair were randomised across runs between the left and right eyes. After a practice run, all participants took part in three binocular rivalry runs, each lasting 180 s.

**Binocular rivalry analysis.** Data were first pre-processed: missing data before perceptual transitions were assigned to the subsequent percept. Responses of < 200 ms, mixed responses or missing data with no response were removed from the analysis. Because of the skewedness of phase durations (Fig. 1b), the median phase durations were calculated²⁸. To quantify sensory eye dominance, an 'eye dominance index' (EDI) was calculated as:

$$EDI = [(d_{DE} - d_{NDE})/(d_{NDE})] \times 100, \quad (1)$$

where $d_{DE}$ is the median phase duration through the dominant eye and $d_{NDE}$ is the median phase duration through the non-dominant eye. The EDI quantifies the percentage increase in median phase duration of the dominant over the non-dominant eye.

**Visual stimulation in MRI scanner.** The visual stimulation paradigm continuously engaged neural pathways involved in visual processing (Fig. 3a). Identical visual stimulation was delivered to the dominant eye or the non-dominant eye. The non-viewing eye was occluded with a translucent, form-depriving occluder²¹. In order to minimise head motion, participants practised switching the occluder to the non-viewing eye prior to scanning and were encouraged not to give any verbal responses while in the scanner. To estimate GABA levels, data were analysed as a sustained visual stimulation design across the entire scan (128 spectral averages) during which a behavioural task at fixation was continuously present. Visual stimuli were comprised of the 'checkerboard' blocks, showing high contrast flashing checkerboards, and a 'fixation' block, showing an active fixation task presented on a black background. Participants performed an attention demanding task at fixation, during which they monitored the appearance of a red target and pressed a button on a button box as soon as possible. The task was performed continuously and irrespective of checkerboard or fixation blocks. Blocks with checkerboards present or absent were treated as the same for the MRS analysis, because the continuous fixation task alone is a form of visual stimulation that demands attention⁴⁸ from monocular units. As a control analysis, we also performed the neurochemical analysis by subdividing the run into 'checkerboard' and 'fixation' blocks, and calculating the difference between block types.

Stimuli were generated on a Mac Book Pro laptop using Matlab and Psychtoolbox-3⁴⁷ and displayed on a gamma-linearised Eiki LC-XL100 projector (resolution: 1024 × 768, refresh rate: 60 Hz). Participants viewed the stimuli through 45° angled mirrors on a back-projection screen (viewing distance: 60 cm). Visual stimuli
were full-field checkerboards, contrast reversing at 8 Hz (stimulus size = 19.82° × 14.25°, 8 Hz flicker, mean luminance = 385 cd/m²; 50% contrast). The fixation condition was a uniform black screen (2.33 cd/m²) with a fixation dot task. Each run consisted of four alternations of fixation (64 s) followed by flashing checkerboards (64 s). A central fixation dot was displayed (white with black border, size = 0.75°) throughout the experiment. Participants were instructed to maintain fixation and press any button on a MRI-safe button box when the marker turned red (500 ms, ~ once in every 3 s).

**MR protocol.** Magnetic Resonance data were acquired using a 7 T whole body (Siemens, Erlangen) MR-scanner with a Nova Medical head coil (single transmit, 32 receive channels). A T1-weighted structural scan was acquired for each participant with a 1-mm isotropic resolution (MPRAGE, repetition time TR = 2.2 s, inversion time $T_1$ = 1.05 s, echo time TE = 2.82 ms, FOV = 192 × 192 × 176 mm, flip angle = 7°, total acquisition time = 171 s) to permit placement of the visual cortex voxel-of-interest (VOI). A 2 × 2 × 2 cm (8 cm³) MRS VOI was positioned to target the early visual cortex (EVC) in the occipital lobe. The EVC VOI was centred on the calcarine sulcus and at the midline to include V1 in both hemispheres. fMRI and MRS data were acquired using a combined fMRI-MRS sequence²⁵, yielding a slab of EPI (3D EPI, resolution = 4.3 × 4.3 × 4.3 mm; flip angle = 5°, repetition time TR<sub>epi</sub> = 40 ms, TE = 25 ms, FOV = 240 mm, 16 slices) and a MR spectrum upon each TR. The measurements used a single scan approach and were not phase-cycled, although there was a phase cycling step to minimize potential unwanted signals from outside of the VOI⁴⁹. For each experimental run, 128 spectral averages were acquired using short-echo semi-localisation by adiabatic selective refocusing (semi-LASER) pulse sequence (TE = 36 ms, TR<sub>mrs</sub> = 4 s) with VAPOR water suppression and outer volume suppression⁵⁰,⁵¹. Data were collected as single free induction decay per excitation. A dielectric pad measuring 110 × 110 × 5 mm³ containing a suspension of Barium Titanate (BaTiO₃) and deuterated water (mass-mass ratio of 3:1) was placed behind the occiput of each participant⁵²,⁵³ to increase the effective transmit field efficiency close to the pad⁵⁴.

**fMRI analysis.** fMRI data were analysed using FEAT (FMRI Expert Analysis Tool) v.6.00, part of the FSL software distribution (FMRIB's Software Library, www.fmrib.ox.ac.uk/fsl; RRID:SCR_002823). Pre-processing included motion correction MCFLIRT⁵⁵; non-brain tissue extraction⁵⁶; 5 mm smoothing, grand-mean intensity normalization and high pass temporal filtering (Gaussian-weighted least squares straight line fitting, main experiment = 132 s). Registration of EPI to an initial 2-mm structural image used 6 DOF, followed by registration to the 1-mm isotropic T1-weighted structural image using boundary-based registration (BBR) in FLIRT⁵⁵,⁵⁷. BOLD-change in the MRS-voxel was calculated using Featquery.

**fMRS analysis.** Pre-processing for raw MRS data was performed in MRspa (https://www.cmrr.umn.edu/downloads/mrspa/) and included eddy current correction, frequency alignment to the tNAA singlet at 2.01 ppm and phase correction using a least-square algorithm prior to signal averaging. For each participant, 128 spectral averages were collected for each experimental condition. Metabolites within the chemical shift range of 0.5 to 4.5 were analysed using LCModel⁵⁸,⁵⁹. Group averages were calculated by averaging across LCModel metabolite estimates from each individual. Metabolite levels were estimated using a basis set of alanine (Ala); ascorbate/vitamin C (Asc); aspartate (Asp); glycerophosphorylcholine (GPC); phosphorylcholine (PCho); creatine (Cr); phosphocreatine (PCr); γ-aminobutyricacid (GABA); glucose (Glc); glutamine (Gln); glutamate (Glu); glutathione (GSH); inositol (Ins); lactate (Lac); phosphoeth anolamine (PE); scylo-inositol (sIns); taurine (Tau); N-acetyl-aspartate multiplet (mNAA); N-acetyl-aspartate singlet (sNAA); Acetyl moiety of N-acetylaspartylglutamate (sNAAG); aspartyl moiety of NAAG (mNAAG); glutamate moiety of NAAG (gNAAG). Consistent with our previous studies²⁵,⁶⁰, the basis set allowed splitting of singlets and multiplets of NAA and NAAG spectra to account for differences in transverse relaxation and J-modulation at moderate echo times⁶¹,⁶². Macromolecule inclusion procedures were performed as in Bednarik et al.²³. Macromolecular (MM) spectra were added in the LCModel basis set. MM spectra were measured from the occipital cortex of 3 participants using an inversion recovery sequence (TR = 3 s, TE = 36 ms, inversion time TI = 0.685 s), and were included in the model spectra. The residual signal of the methylene resonance of tCr at 3.93 ppm was removed by post processing and high-frequency noise was suppressed using a Gaussian filter ($\sigma$ = 0.05 s) before including the MM spectrum into the LCModel basis set. A correlation coefficient is calculated to quantify the independence of metabolite estimates from each other. Two metabolites that have an absolute correlation coefficient > 0.5 cannot be separated from each other. In these cases, the sum rather than the pair is reported. Correlation coefficients across all metabolites were determined from the LCModel fitting of semi-LASER spectra. For example, the LCModel correlation coefficients were more negative than – 0.5 for the following pair of metabolites: Cr and PCr (r = – 0.83).

The amount of cerebro-spinal fluid (CSF, 6.2 ± 2.3%), white matter (WM, 50 ± 3.9%) and grey matter (GM, 43.8 ± 2.9%) inside the EVC MRS voxel were estimated using the brain-extracted T1-weighted high resolution anatomy scan. Tissue fractions were determined by using FSL anatomical processing script fsl_anat and automated tissue segmentation (FSL v6.0 FAST⁶³) with 5 mm Full-Width Half Maximum bias field correction. Percentage of tissue types were quantified by using the FSL command line tool fslstats. We identified the impact of BOLD-effects in metabolite spectra by estimating the total Creatine line width at 3.03 ppm (tCrLW) during dominant and non-dominant eye runs⁶⁴. No differences in BOLD-effects between DE and NDE spectra were found (p = 0.82). Metabolites referenced to water were corrected for the amount of CSF using the equation: $[M_{corr}] = [M] \times (1/[1 - f_{CSF}])$, where $M_{corr}$ is the corrected metabolite level, M the metabolite value from LCModel and $f_{CSF}$ the CSF fraction²¹,⁶⁵. GABA was scaled to the unsuppressed tissue water spectrum collected from the same MRS volume and at the start of each MRI session as a metabolite reference⁶⁶. To reduce dependency of the GABA measure on tissue fraction, we applied the α-correction method⁶⁷⁻⁶⁹ by Harris et al.⁶⁷, that first calculates
GABA for a hypothetical all gray matter voxel, and then normalizes the values relative to the group tissue fraction using the equation:

$$c_{GMWMcorr} = \frac{c_{meas}}{(f_{GM} + \alpha f_{WM})} \frac{\mu_{GM} + \alpha \mu_{WM}}{\mu_{GM} + \mu_{WM}}, \quad (2)$$

where $c_{GMWMcorr}$ is the corrected value, $c_{meas}$ is the uncorrected value, $f_{GM}$ and $f_{WM}$ are the participant's GM and WM fractions, and $\alpha$ is the ratio of the metabolite estimate in grey and white matter, set to 0.5, and $\mu_{GM}$ and $\mu_{WM}$ are the GM and WM fractions across the group. An $\alpha$-value of 0.5 has been chosen as white matter is assumed to have half the GABA level as gray matter⁶⁷. Since the estimated GABA levels were not further corrected for tissue relaxation, we report metabolite levels scaled to unsuppressed water as 'GABA:H₂O'. As an alternative approach, GABA was also scaled to the sum of creatine and phosphocreatine ('tCr') signals acquired in the same voxel at the same time ('GABA:tCr')⁷⁰. Referencing to tCr effectively controls for variations in tissue composition and CSF proportion across participants. GABA:tCr results are presented alongside results using GABA:H₂O.

The quality of the MRS fits was quantified by the LCModel Cramér-Rao lower bounds (CRLB), and the criterion of 30% across participants (128 spectra/participant, N = 12) was chosen to exclude metabolites that were less detectable from noise. The averaged tissue corrected metabolite levels across 12 participants (128 spectra/participant) for GABA:H₂O was 0.91 ± 0.12 (GABA:tCr, 0.81 ± 0.10, CRLB = 28.1 ± 3.7%). To assess the stability of GABA measurements over run time, we calculated the intra-subject test-retest reliability of GABA values using the coefficient of variation (CoV) (Supplementary Information, Fig S2). The intra-subject CoV was calculated by dividing the standard deviation of the measurements by their mean level⁷¹. GABA CoV was not affected by run time (Fig. S2b), supporting a stable measurement of GABA over acquisition time, nor was it affected by viewing condition (Fig. S2c).

The interocular difference metric in neural response is calculated for MRS measured GABA and for the %BOLD-signal change using the same equation:

$$\text{Interocular N} = (N_{DE} - N_{NDE})/N_{NDE}, \quad (3)$$

where $N_{DE}$ is the neural response value during dominant eye viewing, $N_{NDE}$ is the neural response value during non-dominant eye viewing.

To quantify the percentage of the MRS VOI in respect to specific visual areas, we measured the percentage of the VOI that overlapped with cytoarchitectonic maps of human post-mortem data from bilateral V1, V2²⁷, V3v and hV4³³ using the FSL command line tool atlasquery (Supplementary Information, S1). No further thresholding was performed on the probabilistic atlases.

**Statistical analysis.** Data were assessed for normality using the Shapiro-Wilk Test. If normality was rejected, (SW-test, $p < 0.05$) non-parametric statistics were used. The Wilcoxon Signed Rank test was used to test for significant differences in the median between two repeated measures. Spearman's Rank correlation coefficients were calculated to evaluate the relationship between two independent variables of interest. The $p$-value was calculated using the exact distribution test. If normality was not refuted, paired t-tests were applied to test for differences between group means. Pearson's Linear Correlation coefficients were computed to assess the relationship between two independent variables. In all cases, the significance value alpha was set to 0.05, uncorrected. Due to the low number of participants, it is possible that extreme data points influenced the data. However, no outliers were identified in metabolite or eye dominance measures using the function boxplot.m in Matlab.

Received: 8 June 2021; Accepted: 27 July 2021

Published online: 23 August 2021

## References

1. Ooi, T. L. & He, Z. J. Sensory eye dominance: Relationship between eye and brain. *Eye Brain* **12**, 25–31. https://doi.org/10.2147/EB.S176931 (2020).

2. Brown, S. A. et al. Prevalence of amblyopia and associated refractive errors in an adult population in Victoria, Australia. *Ophthalm. Epidemiol.* **7**, 249–258 (2000).

3. Friedman, D. S. et al. Prevalence of amblyopia and strabismus in white and African American children aged 6 through 71 months the Baltimore Pediatric Eye Disease Study. *Ophthalmology* **116**(2128–2134), e2121–e2122. https://doi.org/10.1016/j.ophtha.2009.04.034 (2009).

4. Attebo, K. et al. Prevalence and causes of amblyopia in an adult population. *Ophthalmology* **105**, 154–159. https://doi.org/10.1016/S0161-6420(98)91862-0 (1998).

5. Burian, H. M. Symposium; strabismus; adaptive mechanisms. *Trans. Am. Acad. Ophthalmol. Otolaryngol.* **57**, 131–144 (1953).

6. Sireteanu, R. & Fronius, M. Naso-temporal asymmetries in human amblyopia consequence of long-term interocular suppression. *Vis. Res.* **21**, 1055–1063. https://doi.org/10.1016/0042-6989(81)90010-9 (1981).

7. Harrad, R. Psychophysics of suppression. *Eye (Lond.)* **10**(Pt 2), 270–273. https://doi.org/10.1038/eye.1996.57 (1996).

8. Wiesel, T. N. & Hubel, D. H. Single-cell responses in striate cortex of kittens deprived of vision in one eye. *J. Neurophysiol.* **26**, 1003–1017. https://doi.org/10.1152/jn.1963.26.6.1003 (1963).

9. Sengpiel, F., Blakemore, C., Kind, P. C. & Harrad, R. Interocular suppression in the visual cortex of strabismic cats. *J. Neurosci.* **14**, 6855–6871 (1994).

10. Sengpiel, F., Jirmann, K. U., Vorobyov, V. & Eysel, U. T. Strabismic suppression is mediated by inhibitory interactions in the primary visual cortex. *Cereb. Cortex* **16**, 1750–1758. https://doi.org/10.1093/cercor/bhj110 (2006).

11. Maya Vetencourt, J. F. et al. The antidepressant fluoxetine restores plasticity in the adult visual cortex. *Science* **320**, 385–388. https://doi.org/10.1126/science.1150516 (2008).
12. Harauzov, A. et al. Reducing intracortical inhibition in the adult visual cortex promotes ocular dominance plasticity. *J. Neurosci.* **30**, 361–371. https://doi.org/10.1523/JNEUROSCI.2233-09.2010 (2010).

13. Levelt, W. On Binocular Rivalry (Institution for Perception, 1965).

14. Alais, D. B. R. *Binocular Rivalry* (MIT Press, 2005).

15. Blake, R. A neural theory of binocular rivalry. *Psychol. Rev.* **96**, 145–167. https://doi.org/10.1037/0033-295x.96.1.145 (1989).

16. Seely, J. & Chow, C. C. Role of mutual inhibition in binocular rivalry. *J. Neurophysiol.* **106**, 2136–2150. https://doi.org/10.1152/jn.00228.2011 (2011).

17. Robertson, C. E., Ratai, E. M. & Kanwisher, N. Reduced GABAergic action in the autistic brain. *Curr. Biol.* **26**, 80–85. https://doi.org/10.1016/j.cub.2015.11.019 (2016).

18. van Loon, A. M. et al. GABA shapes the dynamics of bistable perception. *Curr. Biol.* **23**, 823–827. https://doi.org/10.1016/j.cub.2013.03.067 (2013).

19. Pitchaimuthu, K. et al. Occipital GABA levels in older adults and their relationship to visual perceptual suppression. *Sci. Rep.* **7**, 14231. https://doi.org/10.1038/s41598-017-14577-5 (2017).

20. Mentch, J., Spiegel, A., Ricciardi, C. & Robertson, C. E. GABAergic inhibition gates perceptual awareness during binocular rivalry. *J. Neurosci.* **39**, 8398–8407. https://doi.org/10.1523/JNEUROSCI.0836-19.2019 (2019).

21. Lunghi, C., Emir, U. E., Morrone, M. C. & Bridge, H. Short-term monocular deprivation alters GABA in the adult human visual cortex. *Curr. Biol.* **25**, 1496–1501. https://doi.org/10.1016/j.cub.2015.04.021 (2015).

22. Stagg, C. J. et al. Relationship between physiological measures of excitability and levels of glutamate and GABA in the human motor cortex. *J. Physiol.* **589**, 5845–5855. https://doi.org/10.1113/jphysiol.2011.216978 (2011).

23. Bednarik, P. et al. Neurochemical and BOLD responses during neuronal activation measured in the human visual cortex at 7 Tesla. *J. Cereb. Blood Flow Metab.* **35**, 601–610. https://doi.org/10.1038/jcbfm.2014.233 (2015).

24. Apsvalka, D., Gadie, A., Clemence, M. & Mullins, P. G. Event-related dynamics of glutamate and BOLD effects measured using functional magnetic resonance spectroscopy (fMRS) at 3T in a repetition suppression paradigm. *Neuroimage* **118**, 292–300. https://doi.org/10.1016/j.neuroimage.2015.06.015 (2015).

25. Ip, I. B. et al. Combined fMRI-MRS acquires simultaneous glutamate and BOLD-fMRI signals in the human brain. *Neuroimage* **155**, 113–119. https://doi.org/10.1016/j.neuroimage.2017.04.030 (2017).

28. Frangou, P., Emir, U.E., Karlaftis, V.M. et al. Learning to optimize perceptual decisions through suppressive interactions in the human brain. *Nat Commun* **10**, 474. https://doi.org/10.1038/s41467-019-08313-y (2019).

27. Levelt, W. J. Note on the distribution of dominance times in binocular rivalry. *Br. J. Psychol.* **58**, 143–145. https://doi.org/10.1111/j.2044-8295.1967.tb01068.x (1967).

28. Carmel, D., Arcaro, M., Kastner, S. & Hasson, U. How to create and use binocular rivalry. *J. Vis. Exp.* https://doi.org/10.3791/2030 (2010).

29. Rackayova, V., Cudalbu, C., Pouwels, P. J. W. & Braissant, O. Creatine in the central nervous system: From magnetic resonance spectroscopy to creatine deficiencies. *Anal. Biochem.* **529**, 144–157. https://doi.org/10.1016/j.ab.2016.11.007 (2017).

30. Sengpiel, F. Motion perception is learned, not innate. *Nat. Neurosci.* **9**, 591–592. https://doi.org/10.1038/nn0506-591 (2006).

31. Kelly, J. G., Garcia-Marin, V., Rudy, B. & Hawken, M. J. Densities and laminar distributions of Kv3.1b-, PV-, GABA-, and SMI-32-immunoreactive neurons in macaque area V1. *Cereb. Cortex* **29**, 1921–1937. https://doi.org/10.1093/cercor/bhy072 (2019).

32. Vanbrederode, J. F. M., Mulligan, K. A. & Hendrickson, A. E. Calcium-binding proteins as markers for subpopulations of gabaergic neurons in monkey striate cortex. *J. Comp. Neurol.* **298**, 1–22. https://doi.org/10.1002/cne.902980102 (1990).

33. Hensch, T. K. Critical period regulation. *Annu. Rev. Neurosci.* **27**, 549–579. https://doi.org/10.1146/annurev.neuro.27.070203.144327 (2004).

34. Mendola, J. D. & Conner, I. P. Eye dominance predicts fMRI signals in human retinotopic cortex. *Neurosci. Lett.* **414**, 30–34. https://doi.org/10.1016/j.neulet.2006.12.012 (2007).

35. Rombouts, S. A., Barkhof, F., Sprenger, M., Valk, J. & Scheltens, P. The functional basis of ocular dominance: Functional MRI (fMRI) findings. *Neurosci. Lett.* **221**, 1–4. https://doi.org/10.1016/s0304-3940(96)13260-2 (1996).

36. de Best, P. B., Raz, N., Dumoulin, S. O. & Levin, N. How ocular dominance and binocularity are reflected by the population receptive field properties. *Investig. Ophthalmol. Vis. Sci.* **59**, 5301–5311. https://doi.org/10.1167/iovs.18-24161 (2018).

37. Mapp, A. P., Ono, H. & Barbeito, R. What does the dominant eye dominate? A brief and somewhat contentious review. *Percept. Psychophys.* **65**, 310–317. https://doi.org/10.3758/bf03194802 (2003).

38. Binda, P. et al. Response to short-term deprivation of the human adult visual cortex measured with 7T BOLD. *Elife* https://doi.org/10.7554/eLife.40014 (2018).

39. Logothetis, N. K. What we can do and what we cannot do with fMRI. *Nature* **453**, 869–878. https://doi.org/10.1038/nature06976 (2008).

40. Chopin, A., Bavelier, D. & Levi, D. M. The prevalence and diagnosis of 'stereoblindness' in adults less than 60 years of age: A best evidence synthesis. *Ophthalm. Physiol. Opt.* **39**, 66–85. https://doi.org/10.1111/opa.12607 (2019).

41. Cooper, P. R. & Mendola, J. D. Abnormal sensory eye dominance in stereoanomalous subjects. *J. Vis. Exp.* **19**, 14. https://doi.org/10.1167/19.13.14 (2019).

42. Porac, C. & Coren, S. Suppressive processes in binocular vision—Ocular dominance and amblyopia. *Am. J. Optom. Phys. Opt.* **52**, 651–657 (1975).

43. Leonards, U. & Sireteanu, R. Interocular suppression in normal and amblyopic subjects: The effect of unilateral attenuation with neutral density filters. *Percept. Psychophys.* **54**, 65–74. https://doi.org/10.3758/bf03206938 (1993).

44. Ooi, T. L., Su, Y. R., Natale, D. M. & He, Z. J. A push-pull treatment for strengthening the 'lazy eye' in amblyopia. *Curr. Biol.* **23**, R309–R310. https://doi.org/10.1016/j.cub.2013.03.004 (2013).

45. Lunghi, C. et al. A new counterintuitive training for adult amblyopia. *Ann. Clin. Transl. Neurol.* **6**, 274–284. https://doi.org/10.1002/acn.698 (2019).

46. Miles, W. R. Ocular dominance: Methods and results. *Psychol. Bull.* **25**, 155–156 (1928).

47. Brainard, D. H. The psychophysics toolbox. *Spat. Vis.* **10**, 433–436 (1997).

48. Schwartz, S. et al. Attentional load and sensory competition in human vision: Modulation of fMRI responses by load at fixation during task-irrelevant stimulation in the peripheral visual field. *Cereb. Cortex* **15**, 770–786. https://doi.org/10.1093/cercor/bhh178 (2005).

49. Tkac, I. & Gruetter, R. Methodology of H-1 NMR spectroscopy of the human brain at very high magnetic fields. *Appl. Magn. Reson.* **29**, 139–157. https://doi.org/10.1007/Bf03166960 (2005).

50. van de Bank, B. L. et al. Multi-center reproducibility of neurochemical profiles in the human brain at 7 T. *NMR Biomed.* **28**, 306–316. https://doi.org/10.1002/nbm.3252 (2015).

51. Oz, G. & Tkac, I. Short-echo, single-shot, full-intensity proton magnetic resonance spectroscopy for neurochemical profiling at 4 T: Validation in the cerebellum and brainstem. *Magn. Reson. Med.* **65**, 901–910. https://doi.org/10.1002/mrm.22708 (2011).

52. Lemke, C. et al. Two-voxel spectroscopy with dynamic B0 shimming and flip angle adjustment at 7 T in the human motor cortex. *NMR Biomed.* **28**, 852–860. https://doi.org/10.1002/nbm.3328 (2015).

53. Luo, W. et al. Permittivity and performance of dielectric pads with sintered ceramic beads in MRI: Early experiments and simulations at 3 T. *Magn. Reson. Med.* **70**, 269–275. https://doi.org/10.1002/mrm.24433 (2013).
54. Teeuwisse, W. M., Brink, W. M. & Webb, A. G. Quantitative assessment of the effects of high-permittivity pads in 7 Tesla MRI of the brain. Magn. Reson. Med. 67, 1285-1293. https://doi.org/10.1002/mrm.23108 (2012).

55. Jenkinson, M., Bannister, P., Brady, M. & Smith, S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825-841 (2002).

56. Smith, S. M. Fast robust automated brain extraction. *Hum. Brain Mapp.* **17**, 143–155. https://doi.org/10.1002/hbm.10062 (2002).

57. Jenkinson, M. & Smith, S. A global optimisation method for robust affine registration of brain images. Med. Image Anal 5, 143–156 (2001).

58. Provencher, S. W. Estimation of metabolite concentrations from localized in vivo proton NMR spectra. Magn. Reson. Med. 30, 672–679 (1993).

59. Provencher, S. W. Automatic quantitation of localized in vivo 1H spectra with LCModel. NMR Biomed. 14, 260–264 (2001).

60. Ip, I. B., Emir, U. E., Parker, A. J., Campbell, J. & Bridge, H. Comparison of neurochemical and BOLD signal contrast response functions in the human visual cortex. J. Neurosci. 39, 7968-7975. https://doi.org/10.1523/Jneurosci.3021-18.2019 (2019).

61. Deelchand, D. K., Henry, P. G., Ugurbil, K. & Marjanska, M. Measurement of transverse relaxation times of J-t coupled metabolites in the human visual cortex at 4 T. Magn. Reson. Med. 67, 891–897. https://doi.org/10.1002/mrm.23080 (2012).

62. Marjanska, M. et al. Localized H-1 NMR spectroscopy in different regions of human brain in vivo at 7T: T-2 relaxation times and concentrations of cerebral metabolites. NMR Biomed. 25, 332–339. https://doi.org/10.1002/nbm.1754 (2012).

63. Zhang, Y., Brady, M. & Smith, S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imaging 20, 45-57. https://doi.org/10.1109/42.906424 (2001).

64. Zhu, X. H. & Chen, W. Observed BOLD effects on cerebral metabolite resonances in human visual cortex during visual stimulation: A functional (1)H MRS study at 4 T. Magn. Reson. Med. 46, 841–847 (2001).

65. Deelchand, D. K. et al. Two-site reproducibility of cerebellar and brainstem neurochemical profiles with short-echo, single-voxel MRS at 3T. Magn. Reson. Med. 73, 1718-1725. https://doi.org/10.1002/mrm.25295 (2015).

66. Gasparovic, C. et al. Use of tissue water as a concentration reference for proton spectroscopic imaging. Magn. Reson. Med. 55, 1219-1226. https://doi.org/10.1002/mrm.20901 (2006).

67. Harris, A. D., Puts, N. A. & Edden, R. A. Tissue correction for GABA-edited MRS: Considerations of voxel composition, tissue segmentation, and tissue relaxations. J. Magn. Reson. Imaging 42, 1431-1440. https://doi.org/10.1002/jmri.24903 (2015).

68. Porges, E. C. et al. Impact of tissue correction strategy on GABA-edited MRS findings. Neuroimage 162, 249–256. https://doi.org/10.1016/j.neuroimage.2017.08.073 (2017).

69. Puts, N. A. J. et al. Reduced GABA and altered somatosensory function in children with autism spectrum disorder. Autism Res. 10, 608–619. https://doi.org/10.1002/aur.1691 (2017).

70. De Graaf, R. A. In-Vivo NMR Spectroscopy: Principles and Techniques (Wiley, 2007).

71. Terpstra, M. et al. Test-retest reproducibility of neurochemical profiles with short-echo, single-voxel MR spectroscopy at 3T and 7T. Magn. Reson. Med. 76, 1083-1091. https://doi.org/10.1002/mrm.26022 (2016).

72. Amunts, K., Malikovic, A., Mohlberg, H., Schormann, T. & Zilles, K. Brodmann's areas 17 and 18 brought into stereotaxic space-where and how variable? Neuroimage 11, 66–84. https://doi.org/10.1006/nimg.1999.0516 (2000).

73. Rottschy, C. et al. Ventral visual cortex in humans: Cytoarchitectonic mapping of two extrastriate areas. Hum. Brain Mapp. 28, 1045–1059. https://doi.org/10.1002/hbm.20348 (2007).

## Acknowledgements

This work was supported by the Medical Research Council (MR/K014382/1), The Royal Society (Dorothy Hodgkin Research Fellowship to BI, University Research Fellowship to HB), and the French National Research Agency (ANR-19-CE28-0008, PlaStiC). We would like to thank Dr William Clarke for advice on MRS quantification. This research was funded in whole, or in part, by the Wellcome Trust [Grant number 203139/Z/16/Z]. For the purpose of open access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

## Author contributions

Conceptualization, I.B.I., C.L. and H.B.; Methodology, I.B.I., U.E.E., C.L. and H.B.; Investigation, I.B.I. and U.E.E.; Writing, I.B.I; Review and editing, I.B.I, U.E.E., C.L., H.B. and A.J.P.; Funding Acquisition, A.J.P. and H.B.; Resources, C.L. and U.E.E.; Supervision, H.B.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-95685-1.

**Correspondence** and requests for materials should be addressed to I.B.I.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_16.png" /></figure>

**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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021