Li Gu1†
Yiyao Wang1†Lei Feng1Saiqun Li1Mengwei Zhang1Qingqing Ye1
Yijing Zhuang1
Zhong-Lin Lu2,3,4
Jinrong Li1*Jin Yuan1*- 1State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China
- 2Division of Arts and Sciences, NYU Shanghai, Shanghai, China
- 3Center for Neural Science, Department of Psychology, New York University, New York, NY, United States
- 4NYU-ECNU Institute of Cognitive Neuroscience, NYU Shanghai, Shanghai, China
Previous studies have demonstrated that orientation-specific deprivation in early life can lead to neural deficits of spatial vision in certain space, and can even result in meridional amblyopia (MA). Individuals with astigmatism are the optimal and natural models for exploring this asymmetric development of spatial vision in the human visual system. This study aims to assess the contrast sensitivity function (CSF) and EEG signals along two principal meridians in participants with regular astigmatism when being optimal optical corrected. Twelve participants with astigmatism (AST group, 20 eyes) and thirteen participants with (MA group, 19 eyes) were recruited in the current study. CSFs and spatial sweep visual evoked potentials (sVEP) were measured with vertical and horizontal sinewave gratings along two principal meridians monocularly. Area under log CSF (AULCSF), spatial frequency threshold corresponding to 80% contrast gratings (SF threshold at 80% ctr), and CSF acuity were calculated from CSF test. In addition, sVEP amplitudes and thresholds were calculated with the recursive least square method. Participants with astigmatism exhibited marked vertical-horizontal resolution disparities even after they were corrected with optimal optical corrections. CSF tests showed that AULCSF along weak meridian (measured with horizontal gratings) was lower than that along strong meridian (measured with vertical gratings) in both groups. Significant meridional disparity of CSF acuity was also found in both groups. In addition, the MA group showed larger meridional disparity compared to the AST group. Spatial sVEP thresholds also supported the existence of marked meridional disparity. Our results suggest that meridian-specific partial deprivation in early life might lead to monocularly asymmetric development of spatial vision in the human visual system. In terms of application, we tested the feasibility and reliability of adopting psychophysical and EEG scalings to investigate the asymmetric development of spatial vision related to astigmatism. These paradigms are potentially applicable to reduce and even eliminate the meridional disparity in the primary visual cortex by adopting perceptual learning or other vision-related interventions.
Introduction
The development of visual system is extraordinary sensitive to early visual experience (Berardi et al., 2000; Knudsen, 2004; Kiorpes, 2016). Abnormal visual experience will result in functional and structural deficits in the visual pathway and cerebral cortex, such as deprivation amblyopia (Kiorpes and McKee, 1999; Kiorpes, 2006). In the past decades, psychophysical, electrophysiological, and neuro-imaging techniques make the exploration of neural mechanism underlying human amblyopia feasible and fruitful. However, less attention has been given to meridional amblyopia (MA) resulting from astigmatism. Animal studies have consistently demonstrated that visual deprivation along one orientation would lead to orientation-specific deficits in the visual pathway and cerebral cortex (Blakemore and Cooper, 1970; Hirsch and Spinelli, 1970). However, systematic and comprehensive investigation of meridian-specific deficits in the human visual system is still lacking.
Individuals with astigmatism are the optimal model for investigating the asymmetric meridional development of spatial vision in the human neural system. Astigmatism is a common condition of refractive error in which the eye’s refractive power differs in various meridians, with maximum and minimum powers mutually perpendicular (Read et al., 2014). There are two categories of astigmatism, regular, and irregular astigmatism. The former is normally from abnormal development of corneal curvature, and the latter is from the influence of ocular disease on the components of optical media. Astigmatism has a high prevalence (over 25%) in eastern Asia population most probably due to the different anatomical structure between eyelid, orbit, and eyeball (He et al., 2004, 2007; Rim et al., 2016; Nakamura et al., 2018; Wang et al., 2019).
For regular astigmatism, the parallel rays of light entering the eye are brought to a focus at two distinct focal lines perpendicular to each other rather than to a single focal point (Read et al., 2014). The asymmetrical input of visual signal in perpendicular meridians will affect visual functions, and may even result in MA (Mitchell et al., 1973; Gwiazda et al., 1985). Patients with MA have substandard corrected visual acuity behaviorally, and more notably, astigmatic individuals and amblyopes showed the abnormal development of spatial vision in the visual system (Fiorentini and Maffei, 1973; Freeman and Thibos, 1973; Mitchell et al., 1973; Freeman, 1975). Freeman and his collaborators (Freeman et al., 1972) first demonstrated that astigmatism may contribute to meridian-specific neural deficits in spatial vision, since certain astigmatic participants showed substantial vertical-horizontal resolution differences, even when being fully corrected optically. Since then, researchers further reported on an extensive investigation of the meridional differences in resolution with contrast sensitivity and electrophysiological evidence (Fiorentini and Maffei, 1973; Mitchell et al., 1973; Freeman, 1975). These results indicated that meridian-specific partial deprivation may lead to an analogous modification in the organization of neurons in the human visual system (Freeman et al., 1972). Following previous studies, we used the meridional disparity to represent the resolution differences between vertical and horizontal meridians. However, psychophysical and electrophysiological evidence from a recent study (Yap et al., 2019) demonstrated that non-amblyopic children with and without astigmatism also showed meridional disparity.
Contrast sensitivity and visual evoked potentials (VEPs) have been widely used to investigate the influence of stimulus orientation on the spatial tuning function in healthy subjects and patients (Freeman and Thibos, 1973; Tobimatsu et al., 1993; Arakawa et al., 2000; Tobimatsu and Celesia, 2006; Yap et al., 2019, 2020a). The computerized contrast sensitivity function (CSF) paradigm assesses spatial vision over a wide range of spatial frequencies and contrast levels, and has been demonstrated to be a suitable and applicable paradigm for detecting and diagnosing deficits in spatial vision by a handful of studies (Zhou et al., 2006; Yenice et al., 2007; Huang et al., 2008). Spatial sweep visual evoked potentials (sVEP) has been adopted to measure spatial acuity (Regan, 1977; Tyler, 1979; Norcia and Tyler, 1985; Good and Hou, 2006; Hou et al., 2011; Norcia et al., 2015). The combination of both two scalings would furnish subjective and objective paradigm to future studies related to the improvement of spatial vision, such as diminishing the meridional disparity via perceptual learning in the future study.
Therefore, the goal of the present study was to verify the feasibility and reliability of combining psycho-physical and EEG scalings to build human model of asymmetrical spatial vision development in individuals with astigmatism. To this end, the computerized CSF and sVEP tests were employed. CSF paradigm was adopted to assess spatial vision over a wide range of spatial frequencies and contrast levels (Zhou et al., 2006; Yenice et al., 2007; Huang et al., 2008) and spatial sVEP was adopted to measure spatial acuity (Regan, 1977; Tyler, 1979; Norcia and Tyler, 1985; Good and Hou, 2006; Hou et al., 2011; Norcia et al., 2015). Here, the Area under log CSF (AULCSF), spatial frequency threshold corresponding to 80% contrast gratings (SF threshold at 80% ctr), and CSF acuity were calculated from CSF tests, and the sVEP thresholds were calculated to evaluate meridional disparity in spatial acuity. In these two different paradigms, grating stimuli were presented either along the vertical or the horizontal meridian. Notably, the meridional disparity was quantified under full optical correction with spectacles or contact lenses to eliminate the influence from optical errors. From this, we expect to explore the asymmetrical visual development in human astigmatism in a comprehensive and systematic way.
Materials and Methods
Participants
Twelve participants with astigmatism (AST group, 20 eyes, 6 males, mean age = 9.10 ± 1.37 years, age range:) and thirteen astigmatic participants with substandard corrected visual acuity (VA; defined as meridional amblyopia, MA group, 19 eyes, 7 males, mean age = 13.08 ± 6.75 years) were also recruited in the current study. Their refractive profiles were summarized in Table 1, and clinical demographics details were provided in Supplementary Table 1. Five amblyopes entered this study within a few years (3.89 ± 1.79) of completing conventional amblyopia therapy (i.e., patching therapy). Individuals with regular astigmatism were recruited under the inclusion criteria (Chuck et al., 2018) for eyes in the AST and MA groups. Participants in the study went through careful ocular health examination, VA, autorefraction, manifest subjective refraction, motility examination (including cover-uncover and alternate-cover testing), near stereopsis (marked as Vision Assessment CooperationTM V01, United States), distance stereopsis (Stereo Optical Distance Randot® Stereotest, United States), and corneal topography assessments. Cases of strabismus or micro-strabismus, ocular diseases, and/or abnormalities were excluded.
Table 1. Summary of the refractive profile of MA and AST groups in this study showing the mean refractive error (in DS and DC), power range (in DS and DC), spherical equivalent (in D), the refractive and occlusion history of the participants.
All participants completed the CSF tests, and nine participants in the AST group (13 eyes) and twelve participants in the MA group (16 eyes) completed the EEG scalings. Other participants did not complete the EEG scalings for personal reasons (e.g., limited time, low willing for receiving EEG scalings). This study followed the tenets of the Declaration of Helsinki and was approved by the Zhongshan Ophthalmic Center Ethics Committee. Informed consent was obtained from all participants (or their guardian) prior to data collection. All the participants were wearing 1 month new spectacles or contact lenses under optimal optical corrections prescribed by the same experienced optometrist (author FL) at Zhongshan Ophthalmic Center, in order to avoid the influence of new optical adaptation and previous incorrect optical fitting.
Procedure
Stimuli were displayed on a gamma-corrected AOC G2460PQU/BR LCD computer monitor (120 Hz refresh rate, 1920 × 1080 resolution, and 53.1 cm × 29.8 cm). All experiments were controlled by a PC running MATLAB (Mathworks Inc. Natick, MA, United States) and Psychophysics Toolbox (Brainard, 1997; Pelli, 1997). A special circuit was used to produce 14-bit gray-level resolution (Li et al., 2003; Li and Lu, 2012). The mean background luminance was 27 cd/m2. During the whole experiment, the participant put the head on a chin rest and viewed the stimuli monocularly in a dimly lit room, with the tested eye watching while another eye occluded.
Participants went through a battery of measurements, including EDTRS VA, CSF tests along two perpendicular meridians (i.e., horizontal and vertical meridian, to match the distribution of regular astigmatism), and sVEP tests along the same two perpendicular meridians. Both the CSFs and sVEPs were measured monocularly.
The Tasks
CSF Test
The CSF paradigm was applied to assess the CSF along two meridians. While measuring CSF, the display subtended 4° × 4° at a viewing distance of 1.50 m for participants. To minimize edge effects, a half-Gaussian ramp (σ = 0.5°) was used to blend the gratings into the background. We measured contrast thresholds in a two-interval forced choice grating detection task at seven spatial frequencies (0.5, 1, 2, 4, 8, 12, and 16 c/d) using a three-down one-up staircase procedure that decreased signal contrast by 10% (multiplied the previous value by 0.9) after every three consecutive correct responses and increased signal contrast by 10% after every incorrect response, converging to a performance level of 79.3% correct (Huang et al., 2008; Lu and Dosher, 2014). Following previous studies, each trial started with a 200-ms fixation cross in the center of the display. This was followed by two 200-ms intervals, signaled by a brief tone in the beginning of each and separated by 500 ms. A grating was randomly presented in one of the two intervals. The other interval was blank. Participants indicated the signal interval by pressing one of the keys (F/J keypress) with index finger. The response also initiated the next trial. A reversal results when the staircase changes from increasing to decreasing contrast or vice versa. Usually, one block of three-down one-up staircase with a 10% step size produces about a dozen reversals in 80 trials (Zhang et al., 2019). To make sure getting enough reversals, we set 100 trials to measure the contrast threshold at each spatial frequency. Following the standard practice, we averaged the contrasts of an even number of reversals to estimate the contrast threshold after excluding the first three or four reversals.
Vertical and horizontal sinusoidal gratings were adopted to measure CSF along two meridians, respectively. The seven spatial frequency conditions were randomly mixed, with two meridian blocks completed in random order under each spatial frequency condition. Thus, there were fourteen blocks with 100 trials per block. A 10 min demo program was applied to run before formal data collection to avoid learning effect. A formal testing session of about 1.5 h was required to complete the CSF test in each eye. The sequence of the CSF tests in the two eyes was randomized and counterbalanced across participants.
EEG Data Acquisition and sVEP Test
Participants were seated in a shielded room. The EEG signals were amplified and digitized using an Active Two 64-channel Amplifier with the 64-channel Cap in accordance to the international 10–20 system (Biosemi, Netherland), which provided fast and simple electrode placement. Signals from 64 electrodes were recorded and the impedance of each electrode was kept below 10 kV. Horizontal and vertical electrooculograms (HEOG and VEOG) were also recorded to monitor eye movements. Two reference electrodes were used. The data were sampled at 2048 Hz and filtered with a 0.16–100 Hz band-pass filter.
For spatial sVEP, sweeps of spatial frequency at high contrast were adopted (see Figure 1), and the parameters of stimulus refer to previous research (Hou et al., 2011). An 80% contrast, 6-Hz, phase-reversing cosine grating (shown as a square-wave grating) was swept from 2 to 16 cyc/deg in 10 linear steps. The sweep duration was 10 s. The display subtended 10° × 10° at a viewing distance of 0.90 m for participants. Spatial sVEP test was also monocularly conducted, and the non-viewing eye was occluded with a black eye patch.
Figure 1. Experimental flowchart. We measured the CSF and sVEP signals along two principal meridians (i.e., horizontal meridian and vertical meridian, to match the distribution of with-the-rule astigmatic subjects in the study) monocularly in all participants. N indicates the number of eyes. See also Supplementary Table 1 for clinical details.
Data Analysis
CSF Data
The AULCSF and CSF acuity were calculated (Dorr et al., 2018). Both the spatial frequency and the contrast sensitivity in the logarithmic value were generated. We computed the area under the log CSF (AULCSF) for spatial frequencies ranging from 1.5 to 18 cpd using the trapezoid method. We also computed CSF Acuity, i.e., the intersection of the CSF with the x axis (where contrast threshold is 100%). Additionally, the SF threshold at 80% ctr was calculated and extracted in order to investigate the relationship between CSF and sVEP tests.
sVEP Data
The EEG was analyzed using a customized toolbox (mfeeg. http://sourceforge.net/p/mfeeg) programmed with MATLAB (The Mathworks, Natick, MA, United States). Continuous EEG recordings were band-pass filtered from 1 to 30 Hz.
The stimulus is presented at a given temporal frequency (6 Hz in this study) that drives visual cortical neurons at that frequency and at exact integer multiples of that frequency, as well as the stimulus is in the visible range. The visual response synchronized to the display is sampled with appropriately positioned leads, and the VEP amplitude versus stimulus intensity function is measured as the stimulus-driven response drops into the background EEG noise. In this study, swept parameter presentations were repeated ten times.
To measure the response functions, sVEP recordings for each 10-s trial were divided into 10 sequential epochs that corresponded to the swept stimulus values. For each epoch, a recursive least square (RLS) adaptive filter was used to generate a series of complex valued spectral coefficients representing the amplitude and phase of response components tuned to its multiples of the stimulus frequency (Tang and Norcia, 1995). The signals from three electrodes (Oz, O1, and O2) were averaged for further analysis and we focused on the second harmonics (12 Hz in the current study) of the temporal frequency (Almoqbel et al., 2008, 2011; Hou et al., 2018). These spectral coefficients for each epoch were coherently averaged across trials for each subject and stimulus conditions (Hou et al., 2018). These functions were also used to estimate thresholds for each subject’s individual conditions.
sVEP Threshold Estimation
Response thresholds were estimated by regression of amplitudes from the trial-average epochs for each swept stimulus condition of each subject (Hou et al., 2018). We applied the regression procedure to the sweep response function of each individual subject. For those individuals whose response functions along two meridians both passed the regression criteria, we calculated the resultant thresholds. The regression criteria was adopted from previous research (Norcia and Tyler, 1985), which chose an SNR (signal to noise ratio) of >3:1 (i.e., the amplitude of the peak response signal has to be at least three times larger than the adjacent noise frequency) as a criterion (Norcia and Tyler, 1985).
We also calculated the amplitudes by averaging the sweep response functions of individual participants for a given stimulus condition. In this analysis, each participant contributed equally to all conditions. Error bars in the figures depicting sweep responses are vector standard errors of the vector mean (Hou et al., 2007).
Statistical Analysis
The statistical software package SPSS (version 19; IBM Corp, Armonk, NY, United States) was used for analysis. In order to account for right and left eye-related data (Glynn and Rosner, 2012), linear mixed models (LMMs) analysis was used to investigate the effect of stimulus meridian (strong or weak meridian) and group on the outcome measures. Bonferroni correction was applied to multiple paired comparisons to correct for family-wise error. According to previous research (Huang et al., 2018), monocular data from two eyes was in the same comparison group in the current study design, making F-test of LMMs no better than two sample t-test. Moreover, t-test using the average of two eyes would not perform best considering the small sample size. Thus, comparisons between two meridians were reported with two sample t-test.
Results
CSF Results
Contrast sensitivity function results along two meridians were shown in Table 2, and meridional disparity was present in the AST group (Figure 2A) and MA group (Figure 2B). Individual CSF result for each participant was provided in Supplementary Figure 1. LMMs analysis was used to investigate the effect of stimulus meridian (strong or weak meridian) and group on the AULCSF and Cutoff acuity.
Table 2. Mean values of sVEP threshold, amplitudes at target frequency (12 Hz), and CSF tests along two meridians.
Figure 2. CSF results along two meridians. Participants showed meridional disparity in the AST group (A) and MA group (B). Error bars stand for ± S.E.M.
We first evaluated the meridional disparity of AULCSF. LMMs analysis showed a significant effect of meridian (F1,52.44 = 46.812, p < 0.001), a marginal significant main effect of group (F1,20.22 = 3.176, p = 0.090), and a marginal significant interaction between the two factors (F1,52.44 = 3.461, p = 0.070). F-test of LMMs based on estimated marginal means indicated a significant meridional disparity of AULCSF (Mdiff = 0.223 ± 0.033, F1,52.44 = 46.812, and p < 0.001). Significant meridional disparity of AULCSF was found in the AST (Mdiff = 0.163 ± 0.032, t19 = 5.030, and p < 0.001; Figure 3A) and MA (Mdiff = 0.284 ± 0.028, t18 = 10.209, and p < 0.001; Figure 3B) groups. In addition, the MA group showed larger meridional disparity of AULCSF compared to the AST group (t37 = 2.834, p = 0.007).
Figure 3. Meridional disparity of AULCSF in the AST group (A) and the MA group (B). Error bars represent ± S.E.M. Two asterisks ** indicate a significance level of p < 0.001.
We then evaluated the meridional disparity of CSF acuity. LMMs analysis showed a significant effect of meridian (F1,51.33 = 18.142, p < 0.001), and group (F1,18.58 = 10.793, p = 0.004). The interaction between the two factors was not significant (F1,51.33 = 1.425, p = 0.238). F-test of LMMs based on estimated marginal means indicated a significant meridional disparity of CSF acuity (Mdiff = 0.107 ± 0.025, F1,51.33 = 18.142, and p < 0.001). Significant meridional disparity of CSF acuity was found in the AST (Mdiff = 0.077 ± 0.033, t19 = 2.362, and p = 0.029; Figure 4A) and MA group (Mdiff = 0.138 ± 0.029, t18 = 4.707, and p < 0.001; Figure 4B). Meanwhile, there was no significant difference between meridional disparity of CSF acuity in two groups (t37 = 1.367, p = 0.180).
Figure 4. Meridional disparity of CSF acuity in the AST group (A) and the MA group (B). Error bars represent ± S.E.M. An asterisk * indicates a significance level of p < 0.05. Two asterisks ** indicate a significance level of p < 0.001.
Furthermore, we evaluated the meridional disparity of contrast sensitivity corresponding to individual SF. LMMs analysis showed a significant main effect of meridian (F1,496.69 = 90.249, p < 0.001), a significant main effect of SF (F6,496.69 = 174.645, p < 0.001), a marginal significant main effect of group (F1,20.57 = 3.053, p = 0.096), a significant interaction effect between meridian and group (F1,496.69 = 6.292, p = 0.012), a significant interaction effect between SF and group (F6,496.69 = 8.155, p < 0.001), a significant interaction effect between meridian and SF (F6,496.69 = 3.060, p = 0.006), and but a non-significant interaction between the three factors (F6,496.69 = 0.563, p = 0.759). F-test of LMMs further showed the existence of meridional disparity of contrast sensitivity corresponding to several spatial frequencies (2, 4, 8, 12, and 16 c/d), F1,496.69 = 11.069, p = 0.001; F1,496.69 = 24.160, p < 0.001; F1,496.69 = 22.923, p < 0.001; F1,496.69 = 22.604, p < 0.001; and F1,496.69 = 25.119, p < 0.001, respectively. Post hoc comparisons between two meridians are listed in Table 2.
sVEP Results
For sVEP data, we first evaluated the meridional disparity of sVEP thresholds (Table 2). LMMs analysis showed a significant effect of meridian (F1,54 = 16.234, p < 0.001), and group (F1,54 = 10.354, p = 0.002). The interaction between the two factors was not significant (F1,54 = 0.050, p = 0.825). F-test of LMMs based on estimated marginal means indicated a significant meridional disparity of sVEP thresholds (Mdiff = 2.445 ± 0.607, F1,54 = 16.234, and p < 0.001). Significant meridional disparity of sVEP threshold was found in the AST (Mdiff = 2.580 ± 0.881, t12 = 2.928, and p = 0.013) and MA (Mdiff = 2.310 ± 0.610, t15 = 3.789, and p = 0.002) groups. Meanwhile, the meridional disparity of sVEP threshold in the MA group was no larger than that in the AST group (t27 = 0.259, p = 0.798).
We then examined the meridional disparity of amplitudes for spatial frequencies presented (Figure 5; also see Table 2). LMMs analysis showed a significant main effect of meridian (F1,665.08 = 64.807, p < 0.001), a significant main effect of SF (F9,665.05 = 34.494, p < 0.001), a significant interaction effect between meridian and SF (F9,665.12 = 4.959, p < 0.001, but a non-significant main effect of group (F1,20.75 = 0.241, p = 0.628). No significant interaction between meridian and group (F1,665.08 = 0.043, p = 0.835), between SF and group (F9,665.05 = 1.025, p = 0.418), or between the three factors (F9,665.12 = 0.950, p = 0.481) was found. Multiple comparison analysis further showed the existence of meridional disparity of amplitudes corresponding to medium SF (5.1, 6.7, 8.2, 9.8, and 11.3 c/d), F1,664.87 = 25.053, p < 0.001; F1,664.89 = 32.381, p < 0.001; F1,664.89 = 23.397, p < 0.001; F1,665.01 = 22.570, p < 0.001; and F1,664.95 = 10.252, p = 0.001, respectively. Post hoc comparisons between two meridians are listed in Table 2.
Figure 5. Spatial sVEP voltage response at the second harmonic of the test stimulus reversal rate is plotted as a function of spatial frequency of the test stimulus in the AST group (A) and the MA group (B). Error bars represent ± S.E.M.
Correlation Between CSF and sVEP Tests
Most interestingly, we found measurable correlations between CSF (SF threshold at 80% ctr) and sVEP (sVEP threshold) tests. Results showed that there was a significant correlation between SF threshold at 80% ctr and the sVEP threshold in the AST (r = 0.578, p = 0.002; Figure 6A) and MA (r = 0.671, p < 0.001; Figure 6C) groups. Moreover, the meridional disparity of SF threshold at 80% ctr and meridional disparity of the sVEP threshold was also significantly correlated in the AST (r = 0.569, p = 0.043; Figure 6B) and MA (r = 0.532, p = 0.034; Figure 6D) groups.
Figure 6. A scatter plot of CSF and sVEP measures. (A) SF threshold at 80% ctr and the sVEP threshold in the AST group; (B) meridional disparity of SF threshold at 80% ctr and meridional disparity of the sVEP threshold in the AST group; (C) SF threshold at 80% ctr and the sVEP threshold in the MA group; and (D) meridional disparity of SF threshold at 80% ctr and meridional disparity of the sVEP threshold in the MA group. An asterisk * indicates a significance level of p < 0.05. Two asterisks ** indicate a significance level of p < 0.001.
Correlation Between Cylindrical Refractive Errors and Meridional Disparities of CSF or sVEP Measures
In a final analysis, we examined the relationship between cylindrical refractive errors and meridional disparities of CSF or sVEP measures. There was no significant correlation between cylindrical refractive errors and meridional disparities of any CSF or sVEP measure (AULCSF, CSF acuity, SF threshold at 80% ctr, and sVEP threshold; all p > 0.050). Consistent with previous study (Yap et al., 2020a), this indicated that cylindrical refractive error alone does not constitute the decisive factor for the level of meridional disparity, and other factors may also be related to the asymmetric development in the visual system.
Discussion
The current study investigated the asymmetric development of spatial vision related to astigmatism in the human visual system via combining psycho-physical and EEG scalings. We assessed the CSF and sVEP along two principal meridians in participants with astigmatism when being optimally corrected with spectacles, and confirmed the horizontal and vertical asymmetry of spatial vision in human astigmatism. This finding was in line with animal studies (Blakemore and Cooper, 1970; Hirsch and Spinelli, 1970) and other human studies (Freeman et al., 1972; Fiorentini and Maffei, 1973; Mitchell et al., 1973; Freeman, 1975), suggesting that the meridian-specific partial deprivation in early life can lead to monocularly asymmetric development of spatial vision in the visual system.
Neurons in the primary visual cortex, responsive to the retina mapping projection from optical input, is highly sensitive to visual experience during the critical period (Wiesel, 1963; Morishita and Hensch, 2008). Astigmatism without appropriate optical corrections before the critical period would permanently modify the visual system and result in monocularly and binocularly abnormal visual perception (Freeman et al., 1972). It has been reported that abnormal visual input in early life would affect the establishment of visual perception, and may even lead to MA (Hubel and Wiesel, 1962; Shapley et al., 2003). In this study, converging evidence from CSF and sVEP results suggested that meridional disparity on participants with astigmatism are of neural, rather than optical, origin. In addition, we found that the meridional disparity of AULCSF was more remarkable in MA group. As the AULCSF is a summary measure of spatial vision (Applegate et al., 1998; Oshika et al., 2006; Lesmes et al., 2010), we speculated that this finding indicated a higher level of abnormal spatial vision in meridional amblyopes than astigmats, and that the phenomenon that astigmats had normal corrected VA but showed meridional disparity of spatial vision might be due to the compensation of neural system. It should be noted that higher level of meridional disparity was only found in the AULCSF measurement but not in other measures (CSF acuity, SF threshold at 80% ctr, sVEP threshold, and sVEP amplitudes). Meridional amblyopes would generally be expected to have much greater magnitude of meridional disparities. We speculated that it may be due to low severity of amblyopia or the treatment history of amblyopes. All the amblyopes were already treated via wearing spectacles prior to the measuring of the VEPs, which helped to improve the visual performance (Gao et al., 2018).
The finding that the horizontal meridian was weaker than vertical meridian in the current study, may be attributed to a horizontal effect. Horizontal effect indicated that the horizontal meridian is less sensitive than the rest of meridians. This finding was consistent with previous study showing that horizontal effect existed in both astigmats and non-astigmats (Yap et al., 2019). Therefore, it should be noted that the horizontal effect might be confounding. All the participants included in the current study were with-the-rule astigmatic, which might limit the explanation. Future research could include individuals with other types of regular astigmatism and further clarify the relationship between astigmatism type and meridional disparity.
Consistent with previous study (Yap et al., 2020a), our results demonstrated that astigmatism alone is not the decisive factor on the magnitude of their meridional disparity. We speculated that meridional disparity of spatial vision might result not only from meridional optical blur alone, but also from many alternative mechanisms, such as neural suppression of amblyopia if existed (Hess et al., 2014). Meanwhile, there exists differences between the current study and Yap’s study. They recruited newly diagnosed amblyopic children who have never worn spectacles, and this current study investigated subjects who have worn spectacles for at least a few months prior to testing. Considering that the duration of astigmatic blur (Keech and Kutschke, 1995) and treatment history (Gao et al., 2018) also influenced the magnitude of meridional disparity, the finding that meridional amblyopes demonstrated the horizontal effect in the current study might suggest a consequence of treatment or recovery. Besides, the horizontal effect was only observed in children (3–7 years old) in Yap’s study, and the meridional disparity was observed in subjects of older age range (7–26 years old) in our study. According to recent studies (Yap et al., 2019, 2020b; Yap and Boon, 2020), the horizontal effect is an indicator of normality. It is important to separate the meridian-specific deficit from the normal physiological phenomenon (i.e., horizontal effect), thus future study could include newly diagnosed amblyopes and meridional amblyopes with other types of regular astigmatism to further clarify it.
For the combination of psycho-physical methods and EEG scalings, the spatial sVEP was adopted, instead of contrast sVEP which measured contrast sensitivity, causing the mismatch of two paradigms. We adopted seven spatial frequencies in the CSF test and ten spatial frequencies in the sVEP test, which also caused the mismatching of grating stimuli in two paradigms. Even though the SF thresholds at 80% ctr and the sVEP thresholds were significantly correlated, it is not applicable to explore the relationship between contrast sensitivity and sVEP amplitude of gratings on each spatial frequency. CSF and sVEP tests with the same set of grating stimuli are expected in future research. Though we set an SNR of >3:1 as a criterion to control the overestimation following previous research (Norcia and Tyler, 1985), it is acknowledged that sweep VEPs might overestimate spatial resolution thresholds (Hamilton et al., 2020).
In the current study design, monocular data from two eyes was included in the same group. A larger sample size and data from randomly selected eye may help to clarify these observations. Other limitations of the study include: (1) the age range (7–26 years old) was relatively wide; (2) some amblyopes had VA reduction in both eyes without an interocular acuity difference significantly; (3) several amblyopes with anisometropia were not specifically excluded from this present study, which may contribute to an alternative neural mechanism of amblyopia, such as neural suppression (Hess et al., 2014); and (4) several amblyopes received occlusion treatment for years, which may already compensate for the magnitude of their meridional disparity.
Individuals with astigmatism are the optimal and natural models for exploring the asymmetric development of spatial vision, and further investigations are expected based on our current findings. Firstly, we could carry out more explorations to investigate the occurrence and development of astigmatism and MA. For instance, we could measure CSF in the real-world scenarios, or apply other neuroimaging tools, such as functional magnetic resonance imaging (fMRI) and functional near infrared spectroscopy (fNIRS). Secondly, we could reduce and even eliminate the visual deficit in the primary visual cortex by adopting perceptual learning or other vision-related interventions. Furthermore, we expect individuals with astigmatism to fully restore visual functions even without optical corrections by taking advantage of visual neural plasticity.
Data Availability Statement
The data generated during the current study is available from the corresponding author on reasonable request.
Ethics Statement
The studies involving human participants were reviewed and approved by Zhongshan Ophthalmic Center Ethics Committee. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
JY, JL, and LG designed the research. LG, YW, and LF performed the research. LG analyzed the data and drafted the manuscript. JY, JL, and Z-LL revised the manuscript. All authors commented on and edited the manuscript, and approved the final version of the manuscript.
Funding
This research was supported by a grant from the Key-Area Research and Development Program of Guangdong Province (2019B010152001) to JY, National Key Research & Development Project (2020YFC2003905), National Natural Science Foundation of China (81770954) to JL, grants from the China Postdoctoral Science Foundation Grant (2019M663255), the Natural science foundation of Guangdong province (2020A1515010610), and Fundamental Research Funds for the Central Universities (20ykpy141) to LG.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyg.2021.595536/full#supplementary-material
References
Almoqbel, F. M., Leat, S. J., and Irving, E. (2008). The technique, validity and clinical use of the sweep VEP. Ophthalmic Physiol. Opt. 28, 393–403. doi: 10.1111/j.1475-1313.2008.00591.x
Almoqbel, F. M., Yadav, N. K., Leat, S. J., Head, L. M., and Irving, E. L. (2011). Effects of sweep VEP parameters on visual acuity and contrast thresholds in children and adults. Graefes Arch. Clin. Exp. Ophthalmol. 249, 613–623. doi: 10.1007/s00417-010-1469-8
Applegate, R. A., Howland, H. C., Sharp, R. P., Cottingham, A. J., and Yee, R. W. (1998). Corneal aberrations and visual performance after radial keratotomy. J. Refract. Surg. 14, 397–407. doi: 10.3928/1081-597X-19980701-05
Arakawa, K., Tobimatsu, S., Kurita-Tashima, S., Nakayama, M., Kira, J.-I., and Kato, M. (2000). Effects of stimulus orientation on spatial frequency function of the visual evoked potential. Exp. Brain Res. 131, 121–125. doi: 10.1007/s002219900274
Berardi, N., Pizzorusso, T., and Maffei, L. (2000). Critical periods during sensory development. Curr. Opin. Neurobiol. 10, 138–145.
Blakemore, C., and Cooper, G. F. (1970). Development of the brain depends on the visual environment. Nature 228, 477–478.
Chuck, R. S., Jacobs, D. S., Lee, J. K., Afshari, N. A., Vitale, S., Shen, T. T., et al. (2018). Refractive errors & refractive surgery preferred practice pattern®. Ophthalmology 125, P1–P104. doi: 10.1016/j.ophtha.2017.10.003
Dorr, M., Elze, T., Wang, H., Lu, Z.-L., Bex, P. J., and Lesmes, L. A. (2018). New precision metrics for contrast sensitivity testing. IEEE J. Biomed. Health Inform. 22, 919–925. doi: 10.1109/JBHI.2017.2708745
Fiorentini, A., and Maffei, L. (1973). Evoked potentials in astigmatic subjects. Vis. Res. 13, 1781–1783. doi: 10.1016/0042-6989(73)90095-3
Freeman, R. D. (1975). Contrast sensitivity in meridional amblyopia. Invest. Ophthalmol. Vis. Sci. 14, 78–81.
Freeman, R. D., Mitchell, D. E., and Millodot, M. A. (1972). A neural effect of partial visual deprivation in humans. Science 175, 1384–1386.
Freeman, R. D., and Thibos, L. N. (1973). Electrophysiological evidence that abnormal early visual experience can modify the human brain. Science 180, 876–878. doi: 10.1126/science.180.4088.876
Gao, T. Y., Anstice, N., Babu, R. J., Black, J. M., Bobier, W. R., Dai, S., et al. (2018). Optical treatment of amblyopia in older children and adults is essential prior to enrolment in a clinical trial. Ophthalmic Physiol. Opt. 38, 129–143. doi: 10.1111/opo.12437
Glynn, R. J., and Rosner, B. (2012). Regression methods when the eye is the unit of analysis. Ophthalmic Epidemiol. 19, 159–165. doi: 10.3109/09286586.2012.674614
Good, W. V., and Hou, C. (2006). Sweep visual evoked potential grating acuity thresholds paradoxically improve in low-luminance conditions in children with cortical visual impairment. Invest. Opthalmol. Vis. Sci. 47:3220. doi: 10.1167/iovs.05-1252
Gwiazda, J., Mohindra, I., Brill, S., and Held, R. (1985). Infant astigmatism and meridional amblyopia. Vis. Res. 25, 1269–1276. doi: 10.1016/0042-6989(85)90042-2
Hamilton, R., Bach, M., Heinrich, S. P., Hoffmann, M. B., Odom, J. V., McCulloch, D. L., et al. (2020). VEP estimation of visual acuity: a systematic review. Doc. Ophthalmol. 142, 25–74. doi: 10.1007/s10633-020-09770-3
He, M., Huang, W., Zheng, Y., Huang, L., and Ellwein, L. B. (2007). Refractive error and visual impairment in school children in rural Southern China. Ophthalmology 114, 374–382. doi: 10.1016/j.ophtha.2006.08.020
He, M., Zeng, J., Liu, Y., Xu, J., Pokharel, G. P., and Ellwein, L. B. (2004). Refractive error and visual impairment in urban children in Southern China. Invest. Opthalmol. Vis. Sci. 45:793. doi: 10.1167/iovs.03-1051
Hess, R. F., Thompson, B., and Baker, D. H. (2014). Binocular vision in amblyopia: structure, suppression and plasticity. Ophthalmic Physiol. Opt. 34, 146–162. doi: 10.1111/opo.12123
Hirsch, H. V. B., and Spinelli, D. N. (1970). visual experience modifies distribution of horizontally and vertically oriented receptive fields in cats. Science 168, 869–871. doi: 10.1126/science.168.3933.869
Hou, C., Good, W. V., and Norcia, A. M. (2007). Validation study of VEP vernier acuity in normal-vision and amblyopic adults. Invest. Opthalmol. Vis. Sci. 48:4070. doi: 10.1167/iovs.06-1368
Hou, C., Good, W. V., and Norcia, A. M. (2018). Detection of amblyopia using sweep VEP vernier and grating acuity. Invest. Ophthalmol. Vis. Sci. 59, 1435–1442.
Hou, C., Norcia, A. M., Madan, A., Tith, S., Agarwal, R., and Good, W. V. (2011). Visual cortical function in very low birth weight infants without retinal or cerebral pathology. Invest. Opthalmol. Vis. Sci. 52:9091. doi: 10.1167/iovs.11-7458
Huang, C.-B., Zhou, Y., and Lu, Z.-L. (2008). Broad bandwidth of perceptual learning in the visual system of adults with anisometropic amblyopia. Proc. Natl. Acad. Sci. 105, 4068–4073.
Huang, J., Huang, J., Chen, Y., and Ying, G. (2018). Evaluation of approaches to analyzing continuous correlated eye data when sample size is small. Ophthalmic Epidemiol. 25, 45–54. doi: 10.1080/09286586.2017.1339809
Hubel, D. H., and Wiesel, T. N. (1962). Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex. J. Physiol. 160, 106–154.
Keech, R. V., and Kutschke, P. J. (1995). Upper age limit for the development of amblyopia. J. Pediatr. Ophthalmol. Strabismus 32, 89–93. doi: 10.3928/0191-3913-19950301-07
Kiorpes, L. (2006). Visual processing in amblyopia: animal studies. Strabismus 14, 3–10. doi: 10.1080/09273970500536193
Kiorpes, L. (2016). The puzzle of visual development: behavior and neural limits. J. Neurosci. 36, 11384–11393. doi: 10.1523/JNEUROSCI.2937-16.2016
Kiorpes, L., and McKee, S. P. (1999). Neural mechanisms underlying amblyopia. Curr. Opin. Neurobiol. 9, 480–486.
Knudsen, E. I. (2004). Sensitive periods in the development of the brain and behavior. J. Cogn. Neurosci. 16, 1412–1425. doi: 10.1162/0898929042304796
Lesmes, L. A., Lu, Z.-L., Baek, J., and Albright, T. D. (2010). Bayesian adaptive estimation of the contrast sensitivity function: the quick CSF method. J. Vis. 10, 1–21. doi: 10.1167/10.3.17
Li, X., and Lu, Z. L. (2012). Enabling high grayscale resolution displays and accurate response time measurements on conventional computers. J. Vis. Exp. 60:e3312.
Li, X., Lu, Z. L., Xu, P., Jin, J., and Zhou, Y. (2003). Generating high gray-level resolution monochrome displays with conventional computer graphics cards and color monitors. J. Neurosci. Methods 130, 9–18.
Lu, Z.-L., and Dosher, B. (2014). Visual Psychophysics: From Laboratory to Theory. Cambridge, MA: The MIT Press.
Mitchell, D. E., Freeman, R. D., Millodot, M., and Haegerstrom, G. (1973). Meridional amblyopia: evidence for modification of the human visual system by early visual experience. Vis. Res. 13, 535–558. doi: 10.1016/0042-6989(73)90023-0
Morishita, H., and Hensch, T. (2008). Critical period revisited: impact on vision. Curr. Opin. Neurobiol. 18, 101–107.
Nakamura, Y., Nakamura, Y., Higa, A., Sawaguchi, S., Tomidokoro, A., Iwase, A., et al. (2018). Refractive errors in an elderly rural Japanese population: the Kumejima study. PLoS One 13:e0207180. doi: 10.1371/journal.pone.0207180
Norcia, A. M., Appelbaum, L. G., Ales, J. M., Cottereau, B. R., and Rossion, B. (2015). The steady-state visual evoked potential in vision research: a review. J. Vis. 15:4.
Norcia, A. M., and Tyler, C. W. (1985). Spatial frequency sweep VEP: visual acuity during the first year of life. Vis. Res. 25, 1399– 1408.
Oshika, T., Okamoto, C., Samejima, T., Tokunaga, T., and Miyata, K. (2006). Contrast sensitivity function and ocular higher-order wavefront aberrations in normal human eyes. Ophthalmology 113, 1807–1812. doi: 10.1016/j.ophtha.2006.03.061
Pelli, D. G. (1997). The VideoToolbox software for visual psychophysics: transforming numbers into movies. Spat. Vis. 10, 437–442.
Read, S. A., Vincent, S. J., and Collins, M. J. (2014). The visual and functional impacts of astigmatism and its clinical management. Ophthalmic Physiol. Opt. 34, 267–294. doi: 10.1111/opo.12128
Regan, D. (1977). Speedy assessment of visual acuity in amblyopia by the evoked potential method. Ophthalmologica 175, 159–164.
Rim, T. H., Kim, S.-H., Lim, K. H., Choi, M., Kim, H. Y., Baek, S.-H., et al. (2016). Refractive errors in koreans: the korea national health and nutrition examination survey 2008-2012. Korean J. Ophthalmol. 30:214. doi: 10.3341/kjo.2016.30.3.214
Shapley, R., Hawken, M., and Ringach, D. L. (2003). Dynamics of orientation selectivity in the primary visual cortex and the importance of cortical inhibition. Neuron 38, 689–699. doi: 10.1016/S0896-6273(03)00332-5
Tang, Y., and Norcia, A. M. (1995). An adaptive filter for steady-state evoked responses. Electroencephalogr. Clin. Neurophysiol. 96, 268–277. doi: 10.1016/0168-5597(94)00309-3
Tobimatsu, S., and Celesia, G. G. (2006). Studies of human visual pathophysiology with visual evoked potentials. Clin. Neurophysiol. 117, 1414–1433. doi: 10.1016/j.clinph.2006.01.004
Tobimatsu, S., Kurita-Tashima, S., Nakayama, M., and Kato, M. (1993). Effects of spatial frequency on transient and steady state VEPs: stimulation with checkerboard, square-wave grating and sinusoidal grating patterns. J. Neurol. Sci. 118, 17–24. doi: 10.1016/0022-510X(93)90239-U
Tyler, C. W. (1979). Rapid assessment of visual function: an electronic sweep technique for the pattern visual evoked potential. Invest. Ophthalmol. Vis. Sci. 18, 703–713.
Wang, M., Ma, J., Pan, L., Chen, T., Wang, H. L., Wang, Y. H., et al. (2019). Prevalence of and risk factors for refractive error: a cross-sectional study in Han and Mongolian adults aged 40–80 years in Inner Mongolia, China. Eye 33, 1722–1732. doi: 10.1038/s41433-019-0469-0
Wiesel, T. N. (1963). Single-cell responses in striate cortex of kittens deprived of vision in one eye. J. Neurophysiol. 26, 1002–1017.
Yap, T. P., and Boon, M. Y. (2020). Electrodiagnosis and treatment monitoring of children with refractive amblyopia. Adv. Ophthalmol. Optom. 5, 1–24. doi: 10.1016/j.yaoo.2020.04.001
Yap, T. P., Luu, C. D., Suttle, C., Chia, A., and Boon, M. Y. (2020a). Effect of stimulus orientation on visual function in children with refractive amblyopia. Invest. Opthalmol. Vis. Sci. 61:5. doi: 10.1167/iovs.61.5.5
Yap, T. P., Luu, C. D., Suttle, C. M., Chia, A., and Boon, M. Y. (2019). Electrophysiological and psychophysical studies of meridional anisotropies in children with and without astigmatism. Invest. Opthalmol. Vis. Sci. 60:1906. doi: 10.1167/iovs.18-25924
Yap, T. P., Luu, C. D., Suttle, C. M., Chia, A., and Boon, M. Y. (2020b). Characterising the orientation-specific pattern-onset visual evoked potentials in children with bilateral refractive amblyopia and non-amblyopic controls. Doc. Ophthalmol. 1–15. doi: 10.1007/s10633-020-09794-9
Yenice, O., Onal, S., Incili, B., Temel, A., Afsar, N., and Tanrida Gcaron, T. (2007). Assessment of spatial-contrast function and short-wavelength sensitivity deficits in patients with migraine. Eye 21, 218–223. doi: 10.1038/sj.eye.6702251
Zhang, P., Zhao, Y., Dosher, B. A., and Lu, Z.-L. (2019). Evaluating the performance of the staircase and quick change detection methods in measuring perceptual learning. J. Vis. 19:14. doi: 10.1167/19.7.14
Keywords: contrast sensitivity, sVEP, astigmatism, meridional amblyopia, spatial vision
Citation: Gu L, Wang Y, Feng L, Li S, Zhang M, Ye Q, Zhuang Y, Lu Z-L, Li J and Yuan J (2021) Meridian-Specific and Post-Optical Deficits of Spatial Vision in Human Astigmatism: Evidences From Psycho-Physical and EEG Scalings. Front. Psychol. 12:595536. doi: 10.3389/fpsyg.2021.595536
Received: 03 December 2020; Accepted: 23 February 2021;
Published: 17 March 2021.
Edited by:
Jiawei Zhou, Wenzhou Medical University, ChinaReviewed by:
Marcelo Fernandes Costa, University of São Paulo, BrazilTiong Peng Yap, Paediatric Optometry Centre IGARD, Singapore
Copyright © 2021 Gu, Wang, Feng, Li, Zhang, Ye, Zhuang, Lu, Li and Yuan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jin Yuan, eXVhbmppbmNvcm5lYUAxMjYuY29t; Jinrong Li, bGlqaW5ncjNAbWFpbC5zeXN1LmVkdS5jbg==
†These authors have contributed equally to this work