Introduction

Sensitive periods are epochs during the lifespan within which effects of experience on the brain are particularly strong (Knudsen, 2004). Non-human animal work has established that structural remodelling (Bourgeois, 1997) and the development of local inhibitory neural circuits strongly links to the timing of sensitive periods (Gianfranceschi et al., 2003; Hensch et al., 1998; Hensch & Bilimoria, 2012; Hensch & Fagiolini, 2004; Takesian & Hensch, 2013). Early visual experience has been shown to fine-tune local inhibitory circuits (Benevento et al., 1992; Chattopadhyaya et al., 2004; Gandhi et al., 2008; Toyoizumi et al., 2013), which dynamically control feedforward excitation (Tao & Poo, 2005; Wu et al., 2022). The end of the sensitive period has been proposed to coincide with the maturation of inhibitory neural circuits (Hensch, 2005; Wong-Riley, 2021; H. Zhang et al., 2018). Within this framework, neural circuit stability following the sensitive period is maintained via a balance between excitatory and inhibitory transmission across multiple spatiotemporal scales (Froemke, 2015; Haider et al., 2006; Maffei et al., 2004; Takesian & Hensch, 2013; Wu et al., 2022). Such an excitatory/inhibitory (E/I) ratio has been studied at different organizational levels, including the synaptic and neuronal levels, as well as for neural circuits (Van Vreeswijk & Sompolinsky, 1996; Wu et al., 2022).

The early experience-dependence of local inhibitory circuit tuning is supported by a large body of work in non-human animals. In particular, studies of the mouse visual cortex has demonstrated a disrupted tuning of local inhibitory circuits as a consequence of lacking visual experience at birth (Hensch & Fagiolini, 2004; Levelt & Hübener, 2012). In addition, dark-reared mice have been shown to have spontaneous neural firing in adulthood (Benevento et al., 1992) and a reduced magnitude of inhibition, particularly in layers II/III of the visual cortex (Morales et al., 2002), suggesting an overall higher level of excitation.

Human neuroimaging studies have similarly demonstrated that visual experience during the first weeks and months of life is crucial for the development of visual circuits (Baroncelli et al., 2011; Lewis & Maurer, 2005; Maurer & Hensch, 2012; Röder et al., 2021; Röder & Kekunnaya, 2021; Singh et al., 2018). As studies manipulating visual experience are impossible in human research, much of our understanding of the experience-dependence of visual circuit development comes from patients who underwent a transient period of congenital blindness due to dense bilateral congenital cataracts. If human infants born with dense bilateral cataracts are treated later than a few weeks from birth, they suffer from a permanent reduction of not only visual acuity (Birch et al., 1998; Khanna et al., 2013) and stereovision (Birch et al., 1993; Tytla et al., 1993) but additionally from impairments in higher-level visual functions such as face perception (Le Grand et al., 2001; Putzar et al., 2010; Röder et al., 2013), coherent motion detection (Bottari et al., 2018; Hadad et al., 2012; Maurer & Lewis, 2017), visual temporal processing (Badde et al., 2020) and visual feature binding (McKyton et al., 2015; Putzar et al., 2007). These visual deficits in congenital cataract-reversal individuals have been attributed to altered neural development due to the absence of vision, as individuals who suffered from developmental cataracts do not typically display a comparable severity of impairments (Lewis & Maurer, 2009; Sourav et al., 2020). While extant literature reported correlations between structural changes and behavioral outcomes in congenital cataract-reversal individuals (Feng et al., 2021; Guerreiro et al., 2015; Hölig et al., 2023; Pedersini et al., 2023), functional brain imaging (Heitmann et al., 2023; Raczy et al., 2022) and electrophysiological research (Bottari et al., 2016; Ossandón et al., 2023; Pant et al., 2023; Pitchaimuthu et al., 2021) have started to unravel the neural mechanisms which rely on visual experience during sensitive periods for development.

Resting-state activity measured via fMRI suggested an increased E/I ratio in the visual cortex of congenital cataract-reversal individuals (Raczy et al., 2022): The amplitude of low frequency (<1 Hz) (blood oxygen level-dependent) fluctuations (ALFF) in visual cortex was increased in congenital cataract-reversal individuals compared to normally-sighted controls when they were scanned with their eyes open. Since similar changes were observed in permanently congenitally blind humans, the authors speculated that congenital visual deprivation resulted in increased E/I ratio of neural circuits due to impaired neural tuning, which was not reinstated after sight restoration (Raczy et al., 2022). Other studies measured resting-state electroencephalogram (EEG) activity and analyzed periodic (alpha oscillations) (Bottari et al., 2016; Ossandón et al., 2023; Pant et al., 2023) as well as aperiodic activity (Ossandón et al., 2023). Both measures pointed towards an higher E/I ratio of visual cortex in congenital cataract-reversal individuals (Ossandón et al., 2023). In recent research, authors have interpreted the slope of the aperiodic component of EEG power spectral density as an indication of the relative level of excitation; the flatter the slope, the higher the assumed E/I ratio (Gao et al., 2017; Lombardi et al., 2017; Medel et al., 2020; Muthukumaraswamy & Liley, 2018; Nanda et al., 2023). In fact, an increasing number of prospective studies in children have reported a flattening of this slope that was interpreted as higher levels of excitation with increasing age (Favaro et al., 2023; Hill et al., 2022). Ossandón et al. (2023), however, observed that in addition to the flatter slope of the aperiodic power spectrum in the high frequency range (20-40 Hz), the slope of the low frequency range (1-19 Hz) was steeper in both, congenital cataract-reversal individuals, as well as in permanently congenitally blind humans. The low frequency range has often been associated with inhibition (Jensen & Mazaheri, 2010; Lozano-Soldevilla, 2018; Lozano-Soldevilla et al., 2014). However, it remains unclear how to reconcile EEG resting-state findings for lower and higher frequency ranges.

Two studies with permanently congenitally blind humans employed Magnetic Resonance Spectroscopy (MRS) to investigate the concentration of both, the inhibitory neurotransmitter Gamma-Aminobutyric Acid (GABA), and the excitatory neurotransmitters Glutamate/Glutamine (Glx) as proxy measures of visual cortex inhibition and excitation, respectively (Coullon et al., 2015; Weaver et al., 2013). Glutamate/Glutamine concentration was significantly increased in the visual cortex of anophthalmic compared to normally-sighted individuals, suggesting increased excitability (Coullon et al., 2015). Preliminary evidence in congenitally permanently blind individuals suggested a decreased GABA concentration in the visual cortex compared to normally-sighted individuals (Weaver et al., 2013). Thus, these MRS studies corroborated the hypothesis that a lack of visual input at birth enhances relative excitation in visual cortex. However, the degree to which neurotransmitter levels recover following sight restoration after a phase of congenital blindness, and how they related to electrophysiological activity, remained unclear.

Here, we filled this gap: we assessed Glutamate/Glutamine (Glx) and Gamma Aminobutyric Acid (GABA+) concentrations using the MEGA-PRESS sequence (Mescher et al., 1998) in individuals whose sight had been restored, on average, after 12 years of congenital blindness. The ratio of Glx/GABA+ concentration was used as a proxy for the ratio of excitatory to inhibitory neurotransmission (Grent-’t-Jong et al., 2022; Liu et al., 2015; Steel et al., 2020; Takei et al., 2016; L. Zhang et al., 2020). Ten congenital cataract-reversal individuals were compared to age-matched, normally sighted controls at rest, both with eyes open and with eyes closed. To link MRS and EEG markers of cortical excitation/inhibition, we subsequently assessed the aperiodic slope of the EEG spectrum in the same subjects across three conditions: at rest with eyes closed, at rest with eyes open, and while viewing visual stimuli which changed in luminance at frequencies ranging from 1 to 30 Hz (Pant et al., 2023). Additionally we calculated the aperiodic intercept, which has previously been linked to broadband neuronal firing (Manning et al., 2009; Musall et al., 2014; Winawer et al., 2013). We predicted an altered visual cortex Glx/GABA+ concentration ratio in the edited MRS signal in congenital cataract-reversal individuals, and a higher intercept as well as altered slope of the EEG aperiodic component. We further hypothesized that neurotransmitter changes would relate to changes in the slope and intercept of the EEG aperiodic activity in the same subjects.

Methods

Participants

We tested two groups of participants. The first group consisted of 10 individuals with a history of dense bilateral congenital cataracts (CC, 1 female, Mean Age = 25.8 years, Range = 11 – 43.5). Participants in this group were all recruited at the LV Prasad Eye Institute (Hyderabad, India) and the presence of dense bilateral cataracts at birth was confirmed by ophthalmologists and optometrists based on a combination of the following criteria: clinical diagnosis of bilateral congenital cataract, drawing of the pre-surgery cataract, occlusion of the fundus, nystagmus, a family history of bilateral congenital cataracts, and a visual acuity of fixating and following light (FFL+) or less prior to surgery, barring cases of absorbed lenses. Absorbed lenses occur specifically in individuals with dense congenital cataracts(Ehrlich, 1948) and were diagnosed based on the morphology of the lens, anterior capsule wrinkling, and plaque or thickness of stroma.

Duration of deprivation was calculated as the age of the participant when cataract removal surgery was performed on the first eye. Two participants were operated within the first year of life (at 3 months and 9 months of age), all other participants underwent cataract removal surgery after the age of 6 years (Mean Age at Surgery = 11.8 years, SD = 9.7, Range = 0.2 – 31.4). All participants were tested at least 1 year after surgery (Mean Time since Surgery = 14 years, SD = 9.1, Range = 1.8 – 30.9) (Table 1).

Clinical and demographic information of the participants with a history of dense bilateral congenital cataracts (CC).

NA indicates that patient’s data for the field were not available. FFL: Fixating and Following Light; CF: Counting Fingers; PL: Perceiving Light. Duration of visual deprivation was calculated by subtracting the date of birth from the date of surgery on the first eye. Time since surgery was calculated by subtracting the date of surgery on the first eye from the date of testing. Visual acuity on the date of testing was measured with the Freiburg Vision Test (FrACT) (Bach 1996).

The second group comprised of 10 normally sighted individuals (SC, 8 males, Mean Age = 26.3 years, Range = 12 – 41.8). Participants across the two groups were individually age matched (± 2 years, t(9) = –0.12, p = 0.91). One additional individual was tested in each group, however, they were excluded from data analysis as their data files were corrupted due to inappropriate file transfer from the scanner. All participants (as well as legal guardians for minors) gave written and informed consent. This study was conducted after approval from the Local Ethical Commission of the LV Prasad Eye Institute (Hyderabad, India) as well as of the Faculty of Psychology and Human Movement, University of Hamburg (Germany).

Data Collection and Analysis

The present study consisted of three data acquisition parts: Magnetic Resonance Spectroscopy (MRS), Electroencephalography (EEG), and visual acuity assessment.

Magnetic Resonance Spectroscopy

Participants underwent MRI and MRS scanning at LUCID Diagnostics in Hyderabad (India) with a 3T GE SIGNA Pioneer MRI machine employing a 24-channel head coil. An attendant was present in the scanning room for the duration of the scan to ensure that participants were comfortable and followed the instructions.

A T1 weighted whole brain image was collected for each participant (Repetition Time (TR) = 14.97 ms, Echo Time (TE) = 6.74 ms, Matrix size = 512 × 512, In-plane resolution = 0.43 × 0.43 mm, Slice thickness = 1.6 mm, Axial slices = 188, Interslice interval= –0.8 mm, Inversion time = 500 ms, Flip angle = 15°). This structural scan enabled registration of every MRS scan to the participants’ anatomical landmarks (Figure 1). For this scan, participants were instructed to keep their eyes closed and stay as still as possible.

Voxel placement for Magnetic Resonance Spectroscopy and electrode placement for Electroencephalography.

a. Position of the frontal cortex (top) and visual cortex (bottom) voxels in a single subject. Figures output from GannetCoRegister. b. Electrode montage according to the 10/20 electrode system with marked occipital electrodes (red) preselected for analysis.

The MRS scans consisted of single-voxel spectroscopy data that were collected using the MEGA-PRESS sequence, which allows for in-vivo quantification of the low-concentration metabolites GABA and glutamate+glutamine (Glu+Gln) (Mescher et al., 1998; Mullins et al., 2014). Due to the spectral overlap of GABA (3.0 ppm) and Glu/Gln (3.75 ppm) with the higher concentration peaks of N-Acetyl Aspartate (NAA) and Creatine (Cr), accurate quantification of GABA and Glu/Gln is challenging. MEGA-PRESS uses spectral editing to obtain these measurements. Spectroscopy data consist of an edit-ON and an edit-OFF spectrum for each voxel, wherein the “ON” and “OFF” refer to whether the frequency of the editing pulse applied is on– or off-resonance with the signal coupled to the GABA complex (applied at approximately 1.9 ppm). Therefore, subtracting repeated acquisitions of the edit-ON and edit-OFF spectra allows for measurement of the magnitude of signals differing in their response to the editing pulse (e.g. GABA), while cancelling out signals that do not (e.g. Cr) (Mescher et al., 1998). Each MEGA-PRESS scan lasted for 8 minutes and was acquired with the following specifications: TR = 2000 ms, TE = 68 ms, Voxel size = 40 mm x 30 mm x 25mm, 192 averages (each consists of two TRs). Additionally, 8 unsuppressed water averages were acquired, allowing all metabolites to be referenced to the tissue water concentration. Concentrations of GABA and Glu/Gln quantified from these acquisitions are respectively referred to as GABA+, due to the presence of macromolecular contaminants in the signal (Mullins et al., 2014), and Glx, due to the combined quantification of the Glu, Gln and Glutathione peaks.

Two MRS scans were collected from the visual cortex, centered on the calcarine sulcus of every participant (Figure 1). The scans were recorded with the participants’ eyes open (EO condition) and eyes closed (EC condition), respectively. During the scans, participants were instructed to lie as still as possible.

To ensure that we were identifying neurochemical changes specific to visual regions, we selected the frontal cortex as a control region (Figure 1) and collected two scans (EO and EC) from the frontal cortex too. The order of the MRS scans was counterbalanced across individuals for both locations and conditions. Two SC subjects did not complete the frontal cortex scan for the EO condition and were excluded from the statistical comparisons of frontal cortex neurotransmitter concentrations.

MRS Data Analysis

All data analyses were performed in MATLAB (R2018b, The MathWorks Inc.). For MRS data analyses we used Gannet 3.0, a MATLAB based toolbox specialized for the quantification of GABA+ and Glx from edited spectrum data (Edden et al., 2014).

GABA+ and edited Glx concentration values were obtained and corrected using the GannetFit, GannetCoRegister, GannetSegment and GannetQuantify functions (Edden et al., 2014). Briefly, reported concentration values were corrected for the differences in GABA concentration and relaxation times between different kinds of tissue in the voxel (grey matter, white matter and cerebrospinal fluid) (Harris et al., 2015). Gannet uses SPM12 to determine the proportion of grey matter, white matter and cerebrospinal fluid in each individual participant’s voxel (Penny et al., 2007). Note that the tissue fraction values did not differ between groups or conditions (all p’s > 0.19, see Supplementary Material S2). GABA+, Glx and Glx/GABA+ values were compared across groups as proxy measures of inhibition, excitation and E/I ratio respectively. The use of Glx/GABA+ as a proxy measure of E/I neurotransmission is supported by a study that observed a regional balance between Glx and GABA+ at 3T (Steel et al., 2020). Further, the Glx/GABA+ ratio has been employed in prior studies of visual (Takei et al., 2016; L. Zhang et al., 2020), cingulate (Bezalel et al., 2019), frontal (Liu et al., 2015) and auditory cortex (Grent-’t-Jong et al., 2022).

To control for potential unspecified visual cortex changes due to eye pathology, as opposed to genuine changes in neurotransmitter ratio, we compared N-Acetyl Aspartate (NAA) concentrations in the visual cortex of CC vs SC individuals. NAA forms one of the most prominent peaks in the MR spectrum (2.0 ppm chemical shift). NAA has been quantified with high reproducibility in the visual cortex(Brooks et al., 1999) and medial-temporal cortex (Träber et al., 2006) of neuro-typical individuals as well as in various pathologies across visual, frontal and temporal cortex (Paslakis et al., 2014), for example, schizophrenia (Mullins et al., 2003). We did not expect to find differences in NAA concentration between CC and SC individuals as it has not been demonstrated to vary in anophthalmia (Coullon et al., 2015) or permanent early blindness (Weaver et al., 2013) in humans. TARQUIN 4.3.11 was employed to analyze the OFF-spectrum data (Wilson et al., 2011) to assess NAA concentration. FID-A toolbox was used to correct the data for phase errors across acquisitions arising from temporal changes in the magnetic field strength or participant motion (Simpson et al., 2017).

Mean Signal-to-Noise Ratio values for GABA+ and Glx in all groups and conditions were above 19 in the visual cortex and above 8 in the frontal cortex (Supplementary Material S3). A recent study has suggested that an SNR value above 3.8 allows for reliable quantification of GABA+ (Zöllner et al., 2021), in conjunction with considering a given study’s sample size (Mikkelsen et al., 2018). Cramer-Rao Lower Bound (CRLB) values, that is, the theoretical lower limit of estimated error, were 30% or lower for NAA quantification in both groups and conditions (Cavassila et al., 2001). Note that CRLB values above 50% are considered unreliable (Wilson et al., 2019). We confirmed the within-subject stability of metabolite quantification by testing a subset of the sighted controls (n=6) 2-4 weeks apart (Supplementary Material S5).

All reported values are water-normalized.

Prior to in-vivo scanning, we confirmed the GABA+ and GABA+/Glx quantification quality with phantom testing (Henry et al., 2011; Jenkins et al., 2019). Imaging sequences were robust in identifying differences of 0.02 mM in GABA concentration. This 0.02 mM difference was documented by Weaver et al (2013) between the occipital cortices of early blind and sighted individuals. The known vs. measured concentrations of both GABA (r = 0.81, p = 0.004) and GABA/Glx (r = 0.71, p = 0.019) showed significant agreement. The detailed procedure and results are described in the Supplementary Material (Supplementary Material S4). The spectra from all individual subjects are shown in Supplementary Material S10.

MRS Statistical Analysis

All statistical analyses were performed using MATLAB R2018b.

We compared the visual cortex concentrations of 3 neurochemicals (GABA+ and Glx from the DIFF spectrum, NAA from the edit-OFF spectrum) between the two groups. For each metabolite we submitted the concentration values from the visual cortices of CC and SC individuals to a group (2 Levels: CC, SC)-by-condition (2 Levels: EO, EC) ANOVA model. To compare the Glx/GABA+ ratio between groups, we additionally submitted this ratio value to a group-by-condition ANOVA. Identical analyses were performed for the corresponding frontal cortex neurotransmitter values. Wherever necessary, post-hoc comparisons were performed using t-tests.

Electrophysiological recordings

EEG data were collected to investigate aperiodic activity in the same participants and on the same day. Data were acquired in three conditions: at rest with eyes open (EO, 3 minutes), at rest with eyes closed (EC, 3 minutes) and during visual stimulation with stimuli that changed in luminance (LU) (Pant et al., 2023). We used the slope of the aperiodic (1/f) component of the EEG spectrum as an estimate of E/I ratio (Gao et al., 2017; Medel et al., 2020; Muthukumaraswamy & Liley, 2018). Further, we compared the intercept of the aperiodic component in the human EEG between groups, as an estimate of broadband neuronal firing activity (Haller et al., 2018; Manning et al., 2009; Miller, 2010).

The EEG was recorded using Ag/AgCl electrodes attached according the 10/20 system(Homan et al., 1987) to an elastic cap (EASYCAP GmbH, Herrsching, Germany) (Figure 1). We recorded 32 channel EEG using the BrainAmp amplifier, with a bandwidth of 0.01–200 Hz, sampling rate of 5 kHz and a time constant of 0.016 Hz /(10 s) (http://www.brainproducts.com/). All scalp recordings were performed against a left ear lobe reference.

Participants were asked to sit as still as possible while EEG was being recorded. First, resting-state EEG data were collected. During the EO condition, participants were asked to fixate on a blank screen. During the EC condition, participants were instructed to keep their eyes closed. The order of conditions was randomized across participants. The EEG data sets reported here were part of data published earlier (Ossandón et al., 2023; Pant et al., 2023).

Subsequently, EEG data were recorded during 100 trials of a target detection task with stimuli that changed in luminance (LU). Stimuli were presented with a Dell laptop, on a Dell 22 inch LCD monitor with a refresh rate of 60 Hz. They were created with MATLAB r2018b (The MathWorks, Inc., Natick, MA) and the Psychtoolbox 3 toolbox (Brainard, 1997; Kleiner et al., 2007). On each trial, participants observed a circle at the center of a black screen, subtending a visual angle of 17 degrees. The circle appeared for 6.25 s and changed in luminance with equal power at all frequencies (0-30 Hz). At the end of every trial, participants had to indicate whether a target square, subtending a visual angle of 6 degrees, appeared on that trial. The experiment was performed in a darkened room (for further details, see (Pant et al., 2023)).

EEG Data Analysis

Data analysis was performed using the EEGLab toolbox on MATLAB 2018b(Delorme & Makeig, 2004). All EEG datasets were filtered using a Hamming windowed sinc FIR filter, with a high-pass cutoff at 1 Hz and a low-pass cutoff at 45 Hz. Eye movement artifacts were detected in the EEG datasets via independent component analysis using the runica.m function’s Infomax algorithm in EEGLab. Components corresponding to horizontal or vertical eye movements were identified via visual inspection and removed (Plöchl et al., 2012).

The two 3 minutes long resting-state recordings (EC, EO) were divided into epochs of 1 s. Epochs with signals exceeding ±120 μV were rejected for all electrodes. We then calculated the power spectral density of the EO and EC resting-state data using the pwelch function (window length = 1000, overlap = 0).

Datasets collected while participants viewed visual stimuli that changed in luminance (LU) were down-sampled to 60 Hz (antialiasing filtering performed by EEGLab’s pop_resample function) to match the stimulation rate. The datasets were divided into 6.25 s long epochs corresponding to each trial. Subsequently, baseline removal was conducted by subtracting the mean activity across the length of an epoch from every data point. After baseline removal, epochs with signals exceeding a threshold of ±120 μV were rejected in order to exclude potential artifacts. Finally, we calculated the power spectral density of the LU data using the pwelch function (window length = 60 samples, overlap = 0).

We calculated the aperiodic (1/f) component of the power spectrum for the EO, EC and LU conditions (Donoghue, Haller, et al., 2020; Schaworonkow & Voytek, 2021). First, we fit the 1/f distribution function to the frequency spectrum of each participant, separately for each electrode. The 1/f distribution was fit to the normalized spectrum converted to log-log scale (range = 1-20 Hz) (Donoghue, Dominguez, et al., 2020; Gyurkovics et al., 2021; Schaworonkow & Voytek, 2021). We excluded the alpha range (8-14 Hz) for this fit to avoid biasing the results due to documented differences in alpha activity between CC and SC individuals (Bottari et al., 2016; Ossandón et al., 2023; Pant et al., 2023). This 1/f fit resulted in a value of the aperiodic slope, an aperiodic intercept value corresponding to the broadband power of 1-20 Hz, and a fit error value for the spectrum of every participant, individually for all electrodes. The visual cortex aperiodic slope and intercept values were obtained by averaging across pre-selected occipital electrodes O1 and O2, giving one value of broadband slope and intercept per participant for the EO, EC and LU conditions (Figure 1).

EEG Statistical Analysis

We compared the average visual cortex aperiodic slope and intercept in separate group (2 Levels: CC, SC) by condition (3 levels: EC, EO, LU) ANOVA models.

Visual acuity

Visual acuity was measured for every participant on the date of testing, using the Freiburg Visual Acuity Test (FrACT) (Bach 1996, Bach 2007, https://michaelbach.de/fract/). Visual acuity is reported as the logarithm of the mean angle of resolution (logMAR, Table 1), wherein higher values indicate worse vision(Elliott, 2016). As in previous studies, we ran a number of exploratory correlation analyses between GABA+, Glx and Glx/GABA+ concentrations, and visual acuity at the date of testing, duration of visual deprivation, and time since surgery respectively in the CC group (Birch et al., 2009; Guerreiro et al., 2015; Kalia et al., 2014; Rajendran et al., 2020). We additionally tested the correlation between the aforementioned metrics and chronological age across the CC and SC groups.

Correlation analyses between MRS and EEG measures

Exploratory correlation analyses between EEG and MRS measures were run separately for CC and SC individuals. We calculated correlations between the aperiodic intercept and GABA+, Glx and Glx/GABA+ concentrations. Further, correlations between aperiodic slope, and the concentrations of GABA+, Glx and Glx/GABA+ were assessed. MRS measures collected at rest with eyes open (EO) and eyes closed (EC) were correlated with the corresponding resting-state EEG conditions (EO, EC). The correlation between EEG data collected while participants viewed flickering stimuli (LU) calculated with GABA+, Glx and Glx/GABA+ concentration measured while participants’ eyes were open at rest. We did not have prior hypotheses as to the best of our knowledge no extant literature has tested the correlation between aperiodic EEG activity and MRS measures of GABA+,Glx and Glx/GABA+. Therefore, we corrected for multiple comparisons using the Bonferroni correction (6 comparisons).

Results

Transient visual deprivation lowered Glx/GABA+ concentration in the visual cortex

The water-normalized Glx/GABA+ concentration ratio was significantly lower in the visual cortex of congenital cataract-reversal (CC) than age-matched, normally sighted control (SC) individuals (main effect of group: F(1,39) = 5.80, p = 0.021) (Figure 2). This effect did not vary with eye closure (main effect of condition: F(1,39) = 2.29, p = 0.139, group-by-condition interaction: F(1,39) = 1.15, p = 0.290). As a control for unspecified effects of surgery or visual deprivation on neurochemistry, the frontal cortex Glx/GABA+ concentration was compared between groups. There was no difference between CC and SC individuals in their frontal cortex Glx/GABA+ concentration (main effect of group: F(1,37) = 0.05, p = 0.82, main effect of condition: F(1,37) = 2.98, p = 0.093, group-by-condition interaction: F(1,37) = 0.09, p = 0.76) (Figure 2).

Edited spectra obtained from Magnetic Resonance Spectroscopy (MRS).

a. Average edited spectra showing GABA+ and edited Glx peaks in the visual cortices of normally sighted individuals (SC, green) and individuals with reversed congenital cataracts (CC, red) are shown. Edited MRS DIFF spectra are separately displayed for the eyes open (EO) and eyes closed (EC) conditions using dashed and solid lines respectively. The standard error of the mean is shaded. Water-normalized GABA+, Glx and Glx/GABA+ concentration distributions for each group and condition are depicted as violin plots on the right. The solid black lines indicate mean values and dotted lines indicate median values. The coloured lines connect values of individual participants across conditions. b. Corresponding average edited MRS spectra and water-normalized GABA+, Glx and Glx/GABA+ concentration distributions measured from the frontal cortex are displayed. c. Correlations between visual cortex Glx/GABA+ concentrations in the visual cortex of CC individuals and visual acuity in logMAR units are depicted for the eyes closed (EC, left) and eyes open (EO, right) conditions.

When separately comparing CC and SC individuals’ GABA+ and Glx concentrations in the visual cortex, we did not find any group difference (GABA+ main effect of group: F(1,39) = 2.5, p = 0.12, main effect of condition: F(1,39) = 0.6, p = 0.43, group-by-condition interaction: F(1,39) = 0.03, p = 0.86; Glx main effect of group: F(1,39) = 2.8, p = 0.103, main effect of condition: F(1,39) = 1.8, p = 0.19, group-by-condition interaction: F(1,39) = 1.27, p = 0.27) (Figure 2). In the frontal cortex, GABA+ and Glx concentrations did not vary either with group or condition (all p values > 0.19) (Figure 2).

Glx/GABA+ concentration measured when CC individuals’ eyes were closed correlated positively with visual acuity on the logMAR scale (r = 0.65, p = 0.044), indicating that CC individuals with higher Glx/GABA+ values had worse visual acuity (Figure 2C, Supplementary Material S5). The same correlation was not significant for the eyes opened condition (r = –0.042, p = 0.908) (Figure 2C). Duration of deprivation and time since surgery did not significantly predict Glx/GABA+, GABA+ or Glx concentrations in the CC group (all p values > 0.088, Supplementary Material S6).

No difference in NAA concentration between CC and SC individuals’ visual cortices

As a control measure to ensure that between-group differences were specific to hypothesized changes in Glx and GABA+ concentrations, we compared the NAA concentration between CC and SC individuals. The NAA concentration did not significantly differ between groups, neither in their visual (main effect of group: F(1,39) = 0.03, p = 0.87, main effect of condition: F(1,39) = 0.31, p = 0.58, group-by-condition interaction: F(1,39) = 0.09, p = 0.76) nor their frontal cortices (main effect of group: F(1,37) = 1.1, p = 0.297, main effect of condition: F(1,37) = 0.14, p = 0.71, group-by-condition interaction: F(1,37) = 0.03, p = 0.86) (Supplementary Material S1, Figure S1).

Transient visual deprivation resulted in a steeper aperiodic slope and higher aperiodic intercept at occipital sites

The aperiodic slope (1-20 Hz), measured via EEG as an electrophysiological estimate of the E/I ratio(Gao et al., 2017; Muthukumaraswamy & Liley, 2018), was compared between CC and SC individuals. The aperiodic slope was significantly steeper i.e. more negative, at occipital electrodes in CC than in SC individuals (F(1,59) = 13.1, p < 0.001) (Figure 3). Eye closure and visual stimulation did not affect the steepness of the aperiodic slope (F(2,59) = 0.78, p = 0.465, group-by-condition interaction: F(2,59) = 0.12, p = 0.885).

Full spectrum and aperiodic activity of the electroencephalogram (EEG).

a. EEG spectra across O1 and O2 with the corresponding aperiodic (1/f) fits for normally sighted individuals (SC, blue, left) and individuals with reversed congenital cataracts (CC, red, right). Spectra of EEG recordings are displayed for the eyes closed (EC) and eyes opened (EO) conditions, as well as while viewing stimuli that changed in luminance (LU). Shaded regions represent the standard error of the mean. b. Aperiodic intercept (top) and slope (bottom) value distributions for each group and condition are displayed as violin plots. Solid black lines indicate mean values, dotted black lines indicate median values. Coloured lines connect values of individual participants across conditions.

The aperiodic intercept (1-20 Hz) was compared between CC and SC individual to estimate group differences in broadband neural activity(Manning et al., 2009; Musall et al., 2014; Winawer et al., 2013) and was found to be significantly larger at occipital electrodes in CC than SC individuals (main effect of group: F(1,59) = 5.2, p = 0.026) (Figure 3). Eye closure did not affect the magnitude of the aperiodic intercept in either group(main effect of condition: F(2,59) = 0.16, p = 0.848, group-by-condition interaction: F(2,59) = 0.11, p = 0.892).

Within the CC group, visual acuity, time since surgery and duration of blindness did not significantly correlate with the aperiodic slope or the intercept (all p’s > 0.083, Supplementary Material S7). Age negatively correlated with the aperiodic intercept across CC and SC individuals, that is, a flattening of the intercept was observed with age. Similar effects of chronological age have been previously observed (Hill et al., 2022; Voytek et al., 2015) (Supplementary Material S6, Supplementary Material S9).

Glx concentration predicted the aperiodic intercept in CC individuals’ visual cortices during ambient and flickering visual stimulation

We exploratorily tested the relationship between Glx, GABA+ and Glx/GABA+ measured at rest and the EEG aperiodic intercept measured at rest and during flickering visual stimulation, separately for the CC and the SC group.

Visual cortex Glx concentration in CC individuals was positively correlated with the aperiodic intercept either when participants had their eyes open during rest (r = 0.91, p = 0.001, Bonferroni corrected) or when they viewed flickering stimuli (r = 0.90, p < 0.001, Bonferroni corrected). Corresponding correlations were not significant for Glx concentrations in the eyes closed condition (r = 0.341, p > 0.99, Bonferroni corrected). By contrast, in SC individuals no significant correlation was observed between visual cortex Glx concentration and aperiodic intercept in any condition (all p’s > 0.99, Bonferroni corrected) (Figure 4).

Exploratory correlation analyses between the aperiodic intercept (1-20 Hz) and glutamate/glutamine (Glx) concentration in the visual cortex.

Correlations between water-normalized Glx concentration and aperiodic intercept are shown for the eyes closed (EC, left), eyes open (EO, middle) and visual stimulation (LU, right) conditions for sighted controls (SC, green, top) and individuals with reversed congenital cataracts (CC, red, bottom). Displayed p values are Bonferroni corrected for multiple comparisons.

A negative correlation between the aperiodic slope and Glx concentration in CC individuals (i.e., steeper slopes with increasing Glx concentration) was observed during visual stimulation, but did not survive correction for multiple comparisons (Supplementary Material S8). No such correlation was observed between Glx concentration and aperiodic slope in the eyes open or closed conditions. Visual cortex GABA+ concentration and Glx/GABA+ concentration ratios did not significantly correlate with the aperiodic intercept or slope in either CC or SC individuals, during any experimental condition (Supplementary Material S8).

Discussion

Research in non-human animals has provided convincing evidence that the ratio of excitation to inhibition (E/I) in visual cortex is reliant on early visual experience (Froemke, 2015; Haider et al., 2006; Hensch et al., 1998; Takesian & Hensch, 2013; Wu et al., 2022). In parallel, research in humans who were born blind due to dense bilateral cataracts, and received delayed sight restoration surgery in childhood or as adults, has found limited recovery of both basic visual and higher order visual functions (Birch et al., 2009; Röder & Kekunnaya, 2021). The present study tested whether neurotransmitter concentrations and electrophysiological markers of cortical E/I ratio depend on early visual experience in humans, and how possible changes in visual cortex E/I ratio relate to visual recovery. First, we employed Magnetic Resonance Spectroscopy (MRS) and assessed Glutamate/Glutamine (Glx) and Gamma-Aminobutyric Acid (GABA+) concentrations, as well as their ratio, in the visual cortex (Shibata et al., 2017; Steel et al., 2020; Takei et al., 2016). Second, the slope and intercept of the aperiodic spectrum of EEG resting-state activity during eye opening and eye closure (Gao et al., 2017; Muthukumaraswamy & Liley, 2018; Ossandón et al., 2023), as well as during flickering stimulation (Pant et al., 2023), were measured over the occipital cortex in the same individuals. The EEG measures allowed us to relate neurotransmitter changes to neural activity changes in congenital cataract-reversal individuals.

We found a lower Glx/GABA+ concentration ratio in the visual cortex of congenital cataract-reversal (CC) individuals as compared to normally sighted controls (SC). Additionally, the slope of the aperiodic EEG power spectrum was steeper for the low frequency range (1-20 Hz), and its intercept was higher in CC than SC individuals. In the CC group, Glx concentration correlated with the intercept of the aperiodic component during ambient and flickering stimulation. Glx/GABA+ concentration ratio during eye closure predicted visual acuity of CC individuals. Together, the present results provide strong evidence for experience-dependent development of the E/I ratio in the human visual cortex, with consequences for behavior.

Altered Glx/GABA+ ratio after sight restoration in congenitally blind humans

Previous MRS studies in the visual cortex of permanently congenitally blind humans reported higher Glx concentrations (Coullon et al., 2015) in five anophthalmic humans, and numerically lower GABA concentrations in congenitally blind humans (Weaver et al., 2013) (n = 9), as compared to normally sighted individuals. These results were interpreted as suggesting a higher E/I ratio in the visual cortex of permanently congenitally blind humans, which would be consistent with the extant literature on higher BOLD activity in the visual cortices of the same population (Bedny, 2017; Röder & Kekunnaya, 2022). We observed a lower Glx/GABA+ ratio in CC individuals, which suggests a lower rather than higher E/I ratio. Our results imply a change in neurotransmitter concentrations as a consequence of restoring vision following congenital blindness. Here, we speculate that due to limited structural plasticity after a phase of congenital blindness, the neural circuits of CC individuals, which had adapted to blindness after birth, employ available, likely predominantly physiological plasticity mechanisms (Knudsen, 1998; Mower et al., 1985; Röder et al., 2021), in order to re-adapt to the newly available visual excitation following sight restoration.

Structural remodeling for typical E/I balance requires visual experience following birth (Hensch & Fagiolini, 2005; Takesian & Hensch, 2013; H. Zhang et al., 2018) and is linked to a sensitive period (Desai et al., 2002; Hensch & Fagiolini, 2005). A repeatedly documented finding in permanently congenitally blind humans is the increased thickness of visual cortex (Andelin et al., 2019; Guerreiro et al., 2015; Hölig et al., 2023; Jiang et al., 2009), which was not observed in late-blind humans (Andelin et al., 2019). These structural changes in permanently congenitally blind individuals were interpreted as a lack of experience dependent pruning of exuberant synapses and/or reduced myelination typically shifting of the grey-white matter boundary (Natu et al., 2019). In parallel, it was observed that the overproduction of synapses during the initial phase of brain development in non-human primates was independent of experience, but that synaptic pruning, predominantly of excitatory synapses, depended on visual experience (Bourgeois, 1996; Bourgeois et al., 1989). This lack of excitatory pruning has been though to underlie the observed higher excitability of visual cortex due to congenital visual deprivation (Benevento et al., 1992; Huang et al., 2015; Morales et al., 2002), which was interpreted as a homeostatic adjustment of the E/I ratio resulting from the lack of visual feedforward excitation (Huang et al., 2015; Turrigiano & Nelson, 2004). Indeed, sight restoration in non-human primates after several months of congenital bilateral lid suture resulted in higher spontaneous firing in visual association cortex (Hyvärinen et al., 1981). Further, neuronal spiking was increased in dark reared mice when exposed to ambient light (Benevento et al., 1992). Crucially, increased visual cortex thickness (Feng et al., 2021; Guerreiro et al., 2015; Hölig et al., 2023) and higher BOLD activity during rest with the eyes open (Raczy et al., 2022) have been observed for CC individuals as well, suggesting incomplete recovery of cortical structure and function after sight restoration in humans. Thus, the restored feedforward drive to visual cortex after surgery might reach a visual cortex with a lower threshold for excitation, possible due to a relatively higher number of excitatory synapses (Bourgeois, 1996; Morales et al., 2002).

Studies in non-human animals have demonstrated that excitation and inhibition appear to go hand-in-hand (Froemke, 2015; Haider et al., 2006; Isaacson & Scanziani, 2011; Tao & Poo, 2005). Outside of the sensitive period, the visual cortices of CC individuals to maintain neural circuit stability (Lee & Kirkwood, 2019; Turrigiano & Nelson, 2004). An overall reduction in Glx/GABA ratio would counteract the aforementioned adaptations to congenital blindness, e.g. a lower threshold for excitation, which might come with the risk of runaway excitation in the presence of restored visually-elicited excitation. Previous MRS studies have observed a reduction of GABA concentrations in visual cortex (Lunghi et al., 2015) and an increase in the BOLD response (Binda et al., 2018) following monocular blindfolding in adulthood, which was interpreted as indicative of adult homeostatic plasticity. Further, studies in adult mice have provided support for a homeostatic adjustment of the E/I ratio following prolonged changes in neural activity (Chen et al., 2022; Goel & Lee, 2007; Keck et al., 2017; Whitt et al., 2013). For example, a long period of decreased activity following enucleation in adult mice commensurately decreased inhibitory drive (Keck et al., 2011), primarily onto excitatory neurons (Barnes et al., 2015). In line with the lowered Glx/GABA+ ratio being a compensatory measure, the observed correlation of the Glx/GABA+ ratio during eye closure and visual acuity suggests that the more successful this downregulation of the E/I ratio, the better the visual recovery.

The correlation of a lower Glx/GABA+ ratio with better visual acuity is reminiscent of a previously identified correlation in a larger group of CC individuals between decreased visual cortex thickness and better visual acuity (Hölig et al., 2023). Hence, CC individuals with more advanced structural normalization appear to have a better starting point for functional recovery mediated by physiological plasticity.

Glx correlated with the aperiodic intercept of EEG resting-state activity during visual stimulation in congenital cataract-reversal individuals

An increased intercept of the aperiodic component of EEG activity was observed in the same CC individuals who underwent MRS assessment, irrespective of eye opening or eye closure, as well as during flickering stimulation. The intercept of the aperiodic component has been linked to overall neuronal spiking activity (Manning et al., 2009; Musall et al., 2014) and fMRI BOLD activity (Winawer et al., 2013).

Moreover, the slope of the aperiodic component for the low frequency range (1-20 Hz) was steeper in CC individuals. By contrast, the slope of the higher (20-40 Hz) frequency range was flatter in CC than SC individuals (Ossandón et al., 2023). Higher frequencies (such as 20-40 Hz) have been predominantly associated with local circuit activity and feedforward signaling (Bastos et al., 2018; Van Kerkoerle et al., 2014); the increased 20-40 Hz slope may therefore signal increased spontaneous spiking activity in local networks. We speculate that the steeper slope of the aperiodic activity for the lower frequency range (1-20 Hz) in CC individuals reflects the concomitant increase in inhibition.

Interestingly, in CC individuals, the intercept of the aperiodic activity was highly correlated with the Glx concentration during rest with eyes open, and during flickering stimulation (also see Supplementary Material S11). Based on the assumption that the aperiodic intercept reflects broadband firing (Manning et al., 2009; Winawer et al., 2013), this suggests that the Glx concentration might be related to broadband firing in CC individuals during active and passive visual stimulation. As this is the first study testing a relationship between aperiodic EEG parameters and MRS parameters in humans, we interpret this finding based on results from non-human animals. In typically developed adults, sparse neuronal coding (Toyoizumi et al., 2013) after synaptic pruning reflected decorrelated activity (Chini et al., 2022; Sompolinsky et al., 2001; Trägenap et al., 2023; Vinje & Gallant, 2000). Given that visual deprivation resulted in interrupted pruning of predominantly excitatory synapses in monkeys’ visual cortices (Bourgeois, 1996, 1997), we speculate that in the absence of adequate inhibitory tuning for incoming visual input, restored feedforward drive might result in more correlated spiking activity in CC individuals. Indeed, poorly tuned visual responses were evidenced in prior work with CC individuals, including larger (population) receptive fields (Heitmann et al., 2023) and a broader scalp topography of the first cortical visual event related potential (Sourav et al., 2018), as well as higher gamma band activity (Ossandón et al., 2023), which has been linked to firing in feedforward circuits (Van Kerkoerle et al., 2014). Electrophysiological studies with CC individuals have additionally demonstrated a reduced stimulus-selectivity of visual association cortex (Le Grand et al., 2001; Röder et al., 2013; Segalowitz et al., 2017), consistent with imprecise neural representations in this population.

Why might an enhanced inhibitory drive, as indicated by the lower Glx/GABA ratio and, possibly, the 1-20 Hz range of aperiodic EEG activity, not allow for complete recovery of neural tuning? As mentioned above, E/I balance is reflective of multiple interacting mechanisms, and typically maintained by both feedforward inhibition and feedback-mediated inhibition (Froemke, 2015; Isaacson & Scanziani, 2011; Keck et al., 2017; Wu et al., 2022). Feedback connectivity typically emerges later than feedforward connectivity (Batardière et al., 2002; Burkhalter, 1993), and thus is likely more susceptible to damage from early sensory deprivation (Magrou et al., 2018; Yusuf et al., 2022). The markedly reduced alpha oscillatory activity during rest (Ossandón et al., 2023) and in response to white-noise stimulation (Pant et al., 2023) has been interpreted as indirect evidence for deficits in recurrent connectivity within the visual system of CC individuals. Moreover, deficits in the higher harmonic responses and intermodulation frequency responses in steady state visual-evoked responses (Pitchaimuthu et al., 2021), both indicative of bidirectional processing in visual cortex (Kim et al., 2011), have been reported as additional evidence for altered recurrent circuitry in CC individuals. Therefore, the downregulated E/I ratio in CC individuals reflected in the present study likely differs from fine-tuned inhibition within recurrent networks necessary for typical visual processing.

Aperiodic EEG activity correlated with chronological age in sight recovery individuals and normally sighted controls

In both the CC and SC group, the intercept of the aperiodic component of EEG activity decreased with chronological age, irrespective of condition. This chronological age effect replicates corresponding earlier reports in healthy populations (Hill et al., 2022; Voytek et al., 2015). The effect of chronological age in the present study is, again, similar to a correspondingly preserved effect of chronological age on the cortical thickness in CC individuals in cross-sectionally assessed MRI data (Hölig et al., 2023). The reduction of visual cortex thickness during childhood, after an initial increase, is a well-documented trend in brain development (Gilmore et al., 2020; Natu et al., 2019). Thus, we speculate that some typical developmental changes emerge despite aberrant visual experience.

Limitations

To the best of our knowledge, the present study is the first assessment of neurotransmitter concentrations within the visual cortex of sight recovery individuals with a history of congenital blindness by employing non-invasive MRS.

We are aware that MRS has a low spatial specificity. Moreover, MRS measures do not allow us to distinguish between presynaptic, postsynaptic and vesicular neurotransmitter concentrations. However, previous work has validated the link between MRS measures of GABA and Glutamate and the activity of inhibitory and excitatory neuronal assemblies, respectively, both in rodents and non-human primates (Bielicki et al., 2004; Takado et al., 2022). Further, our phantom testing showed high correlations between the experimentally varied metabolite concentrations and the extracted GABA and Glx concentrations, validating the employed assessment and analysis pipelines (Figure S3). Finally, all reported group differences in MRS parameters were specific to visual cortex, and were not found for the frontal control voxel.

The sample size of the present study is relatively high for the rare population, but undoubtedly, overall, rather small. We nevertheless think that our results are valid. Our findings neurochemically (Glx andGABA+ concentration), and anatomically (visual cortex) specific. The MRS parameters varied with parameters of the aperiodic EEG activity and visual acuity. The group differences for the EEG assessments corresponded to those of a larger sample of CC individuals (n=38) (Ossandón et al., 2023), and effects of chronological age were as expected from the literature.

Conclusion

The present study in sight recovery individuals with a history of congenital blindness indicates that E/I balance is a result of early experience, and crucial for human behavior. We provide evidence that E/I balance in sight recovery individuals is altered even years after surgery, which might result from the previous adaptation to congenital blindness.

Conflict of interest

Dr. Sunitha Lingareddy is the Managing Director Radiology at Lucid Medical Diagnostics, Hyderabad, India. All other authors have no conflicts to declare.

Author contributions

RP, KP, JO, JF and BR conceptualized the study. RP, KP, IS and PR collected the data. IS, PR and RK diagnosed, recruited and provided clinical assessments of participants. RP, JO and KP analyzed the data. BR and JF supervised data analysis and methodological decisions. BR, RK and SL provided infrastructure, resources and funding. RP and BR wrote the original draft of the manuscript, and all authors provided edits and reviews on the final draft of the manuscript.

Acknowledgements

We thank the technical staff of the Lucid Medical Diagnostics Center, Banjara Hills, Hyderabad, India, in particular Mr. Balakrishna Vaddepally, for technical assistance during collection of MRS/MRI data. We would like to acknowledge Dr. Suddha Sourav for technical support, and Ms. Prativa Regmi for assistance with phantom testing and data collection. We are grateful to D. Balasubramanian of the L.V. Prasad Eye Institute for initiating and supporting our research. The study was funded by the German Research Foundation (DFG Ro 2625/10-1 and SFB 936-178316478-B11) to Brigitte Röder. RP was supported by a PhD student fellowship from the Hector Fellow Academy GmbH.