Impact Factor 3.582

The world's most-cited Neurosciences journals

Original Research ARTICLE

Front. Aging Neurosci., 05 February 2018 | https://doi.org/10.3389/fnagi.2018.00014

Aging Potentiates Lateral but Not Local Inhibition of Orientation Processing in Primary Visual Cortex

  • 1Hefei National Laboratory for Physical Sciences at Microscale and School of Life Science, University of Science and Technology of China, Hefei, China
  • 2Department of Physiology and Pharmacology, Medical School of Ningbo University, Ningbo, China
  • 3Brainnetome Center and National Laboratory of Pattern Recognition, Institute of Automation, and Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Beijing, China
  • 4Department of Basic Medicine, Medical School, Kunming University of Science and Technology, Kunming, China
  • 5School of Computer and Information, Hefei University of Technology, Hefei, China
  • 6State Key Laboratory of Brain and Cognitive Science, Institute of Biophysics, Chinese Academy of Sciences, Beijing, China

Aging-related declines in vision can decrease well-being of the elderly. Concerning early sensory changes as in the primary visual cortex, physiological and behavioral reports seem contradictory. Neurophysiological studies on orientation tuning properties suggested that neuronal changes might come from decreased cortical local inhibition. However, behavioral results either showed no clear deficits in orientation processing in older adults, or proposed stronger surround suppression. Through psychophysical experiments and computational modeling, we resolved these discrepancies by suggesting that lateral inhibition increased in older adults while neuronal orientation tuning widths, related to local inhibition, stayed globally intact across age. We confirmed this later result by re-analyzing published neurophysiological data, which showed no systematic tuning width changes, but instead displayed a higher neuronal noise with aging. These results suggest a stronger lateral inhibition and mixed effects on local inhibition during aging, revealing a more complex picture of age-related effects in the central visual system than people previously thought.

Introduction

Aging has a profound impact on human visual function, where early neuronal processing and perceptual abilities modifications with age have been intensively studied in the last 25 years (Owsley, 2011; Andersen, 2012). Age-related changes in cortical inhibition have been suggested to be involved in visual functions decline. But as the field currently stands, there is no clear unifying account of early stages changes as in spatial vision.

Neurophysiological studies suggested that a reduction in inhibitory function might be related to age-related perceptual dysfunctions in spatial vision (Schmolesky et al., 2000; Leventhal et al., 2003). Nevertheless, some behavioral studies have shown that there are low-level processes of visual functions which are inconsistent with the physiological interpretation of decreased inhibition: (1) psychophysical sensitivities to orientation processing in older populations are not systematically and substantially changed when compared to younger adults (Betts et al., 2007; Delahunt et al., 2008; Govenlock et al., 2009), while the neurophysiological counterpart strongly suggests such a deficit; (2) contrast surround suppression in older adults was found increased, which was interpreted with stronger inhibitory interactions (Karas and McKendrick, 2015). Surround suppression could be corresponding to neurons' response with stimuli at extra classical receptive field (extra-CRF), involving lateral inhibitory interactions between cortical “channels” (Gilbert, 1992; Spillmann and Werner, 1996; Angelucci et al., 2002; Angelucci and Bressloff, 2006), while neurophysiological inference was drawn from classical receptive field (CRF) properties of cells in the visual cortex (area V1). Since the changed cellular properties are commonly attributed to intracortical inhibition (Somers et al., 1995; Shapley et al., 2003), it appears that intracortical local inhibition and lateral inhibition could have underwent heterogeneous modifications with old age.

These two levels of interactions can be accessed non-invasively through the use of psychophysics measures and computational modeling (Gilbert and Wiesel, 1990; Tzvetanov and Womelsdorf, 2008; Tzvetanov, 2012). We used the tilt illusion (TI) paradigm where the presence of an orientated surround stimulus biases the perceived orientation of a simultaneously presented center test. Current opinions propose that tilt repulsion effects can be explained through lateral inhibition between spatially arranged orientation hypercolumns of neurons in V1 (Georgeson, 1973; Wenderoth and Smith, 1999; Kapadia et al., 2000; Figure 1A, red). Specifically, in the classical center-surround tilt repulsion effect (Gilbert and Wiesel, 1990; Angelucci and Bressloff, 2006; Figure 1A), one can infer the local interactions through the indirect measure of the orientation tuning widths of a theoretical population of neurons (Figures 1A,B, top), which is directly related to the local inhibition within one hypercolumn (Somers et al., 1995), but also the lateral interactions strengths between different hypercolumns (blue arrows in Figure 1A). Both local and lateral interactions influence the tilt repulsion curve (Figure 1B, bottom), that describes the misperception of the subject of a vertically oriented center target when the orientation of the surround is varied (Figure 1A, gray stimulus). Through model adjustment to the behavioral data, estimated mean values of these physiological variables of the neural population was accessed for each person.

FIGURE 1
www.frontiersin.org

Figure 1. Illustration of center-surround stimulus, orientation hypercolumns, population tuning curves and behavioral outcome. (A) Illustration of the stimulus, a center Gabor patch surrounded by an annulus of oriented grating, padded with orientation hypercolumns (red circles); the colored lines in right-top circle depict different preferred orientations of the local neuronal population; on the right, illustration of local interaction within an orientation hypercolumn, with three example neurons with different preferred orientations. Two inter-neurons, and their local connections (blue/green are inhibitory/excitatory connections); blue arrows depict inhibitory lateral interactions between hypercolumns. (B) Theoretical population orientation tuning curves (top) together with typical tilt repulsion curve (bottom) describing the orientation misperception of a vertical center stimulus as a function of the orientation of the surround; zero is vertical and positive values are clockwise tilts.

Here, we carried out psychophysical measures in younger and older populations in order to characterize their center-surround tilt repulsion at two spatial frequencies (SFs), one near best perceptual sensitivity and one at higher SF, obtained from the individual contrast sensitivity function (CSF). Then, we performed modeling of subjects' perception assuming it is obtained from decoding these V1 neuronal activities. Contrary to the results from old adult animals suggesting that age-related changes of visual perception in older adults might stem from a reduction in inhibitory functions, the combination of psychophysics measures and modeling showed different aging-related effects on different kinds of inhibition: (1) an increased lateral inhibition with age and (2) maintained local orientation tuning widths, and thus local inhibition (Somers et al., 1995; Carandini and Ringach, 1997; Shapley et al., 2003), of the neuronal population in older adults. To investigate further this later result, we re-analyzed previously published physiological data about orientation processing from our laboratory and extracted the tuning width of neurons. These physiological results are consistent with the behavioral and computational findings of unchanged orientation tuning widths.

Results

Contrast Sensitivity Functions

Contrast sensitivity function (CSF), representing the inverse of the minimum detectable contrast at each spatial frequency (SF), was measured for each subject prior to the TI measure (Figures 2A,C,D). The contrast sensitivity (CS) was significantly modulated by SF [F(10, 380) = 599.56, P < 0.0001] and CS of the elderly group was significantly lower than the young adults group across all spatial frequencies [F(1, 38) = 7.44, P < 0.05], confirming previous reports(Owsley, 2011). The interaction of the SF and groups was also significant [F(10, 380) = 7.23, P < 0.01]. From each individual CSF, we chose two spatial frequencies (SFs), one near the peak sensitivity (low-SF) and one higher (high-SF) (Figure 2C, squares), with the condition that the sensitivity at the high-SF was not too low for allowing stimulus perception at next stage of measures of center-surround TI.

FIGURE 2
www.frontiersin.org

Figure 2. Examples of stimuli used and results of contrast sensitivity function and tilt illusion measures. (A) Schematic illustration of the contrast sensitivity function (CSF) measure; (B) Example of stimuli in tilt illusion measure. (C) Examples of results of CSF measures and fits in the 2-Alternative unforced Choice (2AuFC) design; squares depict the two chosen SFs for the subsequent tilt measures; (D) Averaged values of contrast sensitivities at the 11 different SFs measured for older adults and younger groups; (E) Tilt illusion results, indicated by orientation bias necessary to perceive the center as vertical, as a function of SOs and SFs (low-SF: circles; high-SF: triangles) for older adults (gray) and younger (black) groups; negative (positive) deviations for negative (positive) surround orientations correspond to repulsion; (F) Orientation thresholds around perceived verticality.

Tilt Illusion Measurement

Perceived verticality of the target grating under center-surround orientation differences of 0, ±15, ±30, and ±75° for each subject was measured at low- and high-SF (Figure 2B). The amount of orientation misperception for both groups was systematically modulated by surround orientations (SOs) and SFs [Surr. Or.: F(3, 114) = 240.41, P < 0.001; SF: F(1, 38) = 75.14, P < 0.001] (Figure 2E), consistent with previous reports (Georgeson, 1973; Wenderoth and Johnstone, 1988; Smith and Wenderoth, 1999). Importantly, analysis limited to ±15° and ±30° SOs, corresponding to the direct form of the TI (Wenderoth and Johnstone, 1988; Smith and Wenderoth, 1999; Wenderoth and Smith, 1999), repulsion, revealed that perception through older adults' visual system exhibited much stronger tilt repulsion effects across all measured SFs in comparison to the younger adults [F(1, 38) = 6.62, P < 0.05; Figure 2E]. The amount of misperception at ±15° SO was significantly different than that at ±30° [F(1, 38) = 53.79, P < 0.0001]. There were no significant interactions between SOs and groups [F(1, 38) = 0.32, P = 0.57] indicating similar young-elder differences at both SOs. The amount of misperception for low- and high-SFs across groups was significantly modulated [F(1, 38) = 89.72, P < 0.001], but there was no interaction between SFs and groups [F(1, 38) = 1.56, P = 0.22] that suggested same effects of SFs across age groups. There was a significant interaction between SOs and SFs [F(1, 38) = 41.28, P < 0.001], while there were no significant interactions among SOs, SFs and Groups [F(1, 38) = 0.15, P = 0.70], again showing no differences of effects across age.

Orientation discrimination thresholds in elder group were larger than younger adults [F(1, 38) = 7.25, P < 0.05], and they were modulated by SOs [F(3, 114) = 11.86, P < 0.0001] and SFs [F(1, 38) = 49.66, P < 0.0001]. There was also a main interaction between SOs and SFs [F(3, 114) = 7.58, P < 0.001] showing a different trend of threshold variation with surround orientation at different SFs (Figure 2F). All other interactions that included age group were not significant [SOs × Groups: F(3, 114) = 1.73, P = 0.17; SFs × Groups: F(1, 38) = 1.04, P = 0.32; SOs × SFs × Groups: F(3, 114) = 1.73, P = 0.17]. We would like to emphasize that the age effect on thresholds does not allow to draw a direct conclusion that there are different orientation threshold between younger and older groups, due to no adjustment of these threshold values to the individual sensitivity (see Discussion).

It is possible that the different high-SF used in younger and elder groups contributed to the differences in misperception between them. We found, as in previous work (Smith and Wenderoth, 1999), that the amount of orientation bias was significantly related to the SF under both ±15° and ±30° SOs conditions for both elder and younger groups (Spearman rank's correlation; the elder: r = 0.58, P = 8.5*10−5 at ±15°, r = 0.39, P = 0.013 at ±30°; the younger adults: r = 0.74, P = 3.6*10−8 at ±15°, r = 0.63, P = 1.2*10−5 at ±30°) (Figures 3A,B). This suggested that adjusting for the variation of TI with SF (Smith and Wenderoth, 1999) still showed that older adults exhibited stronger misperception. Additionally, the old persons group had globally lower high-SF [t(38) = 5.37, P < 0.0001, unpaired t-test] than the younger group, showing that their stronger repulsion effects cannot be explained by this experimental manipulation.

FIGURE 3
www.frontiersin.org

Figure 3. Relation between tilt repulsion, SF and CS at low-SF. Relations between tilt repulsion and SFs under ±15° (triangles) and ±30° (circles) SOs across all measured subjects in older adults (A) and younger (B) groups; Relations between tilt repulsion and CS at low-SF in the elderly group [±15° in (C); ±30° in (D)] and younger (±15° in (E); ±30° in (F)] groups. Error bars are bootstrapped SE.

Relation between Tilt Illusion Bias and Contrast Sensitivity

The elderly group exhibited decreased contrast sensitivity and increased amount of TI compared to the younger adults. To reveal any possible relation between tilt repulsion and peak CS, we correlated tilt and CS measured at low-SF condition, where the SF values for both groups are similar. The results showed that there was no significant relation between bias and CS within each age group (the elder: r = −0.25, P = 0.29 at ±15°, r = −0.36, P = 0.12 at ±30°; the younger adults: r = −0.21, P = 0.38 at ±15°, r = −0.02, P = 0.92 at ±30°) (Figures 3C–F).

Modeling

To further our understanding of early visual processing and the plausible underlying network changes with aging, based on the learned knowledge that neuronal activity in V1 could be the substrate of orientation identification and detection (Georgeson, 1973; Gilbert and Wiesel, 1990; Kapadia et al., 2000; Tzvetanov, 2012), we modeled tilt repulsion effects and contrast sensitivity by including the spatial lateral inhibition in V1. Psychophysical modeling based on the V1 neuronal responses relates each behavioral variable to different tuning characteristics: orientation misperception is obtained through orientation tuning and center-surround inhibitory interactions (Figures 4A–C), while CSF could be obtained through the SF and contrast tuning of the neurons (Figures 4D–G).

FIGURE 4
www.frontiersin.org

Figure 4. V1 model illustration. (A–C) Example and prediction for orientation coding and decoding. (A) Uniform orientation tuning of the neuronal population. (B) Response of the neuronal population to center of 0°. and two different surround orientations of ±30° (top) and ±15° (bottom). (C) Orientation prediction of the model from the population responses for various surround orientations. Examples for three different set of parameters in red, black and blue. (D–G) Example and prediction for (SF) and contrast tuning coding and decoding. (D) SF tuning examples, with the characteristic tuning width decrease with increasing preferred SF. (E) Example of the relation between the minimum contrast semi-saturation constant and preferred SF. (F) Examples of CSF prediction for two sets of model parameters. (G) Examples of contrast response functions in the model for the minimum semi-saturation constant at few preferred SFs (arrows depict half-amplitude constant).

To obtain an insight into the differences between older adults and younger adults (for instance, lateral inhibition, orientation tuning, contrast tuning or a combination of them), we modeled a two-layer feed-forward surround-to-center inhibitory network of V1 neurons. In this model, the TI and contrast sensitivity function (Figures 4C,F) are related through a single parameter, the amount of surround-to-center inhibition (defined as Iinh). The amplitude and shape of the orientation misperception are dependent on Iinh and orientation tuning width (defined as σθ) (Figures 4A–C). The CSF is dependent on Iinh (because the grating stimulus excites center and surround mechanisms) together with the smallest contrast semi-saturation constant (defined as cmin) and its relation to the SF tuning (Figures 4D–G).

We fitted the model first to tilt perception data (two free parameters: Iinh, σθ), and then to CSF data by using lateral inhibition as a fixed parameter (for CSF prediction, there were three free parameters: cmin, a, b; with cmin corresponding to the best neuronal contrast sensitivity across all SFs; e.g., Figure 4E, peak value; see Methods for fitting details). The model was fit to the tilt repulsion data of ±15 and ±30° that provided independent estimates of orientation tuning widths and surround to center inhibitory amplitudes at both low- and high-SFs for each subject (Figures 5A,B). Age differences in tilt repulsion near the peak of perceptual sensitivity (low-SF) were explained by stronger lateral inhibition in the aged group [t(38) = 2.95, P = 0.0054, log-transformed values], while no difference in orientation tuning widths were observed [t(38) = 1.25, P = 0.22]. Additionally, both parameters co-varied with the SF at which it was measured across all population: lateral inhibition increased with higher SFs (r = 0.65, P = 4.4*10−11, n = 80), while neuronal tuning widths tended to decrease with higher SFs (r = −0.31, P = 0.0053, n = 80). This last finding resonates with the established physiological findings of sharper orientation tuning widths for cells with higher preferred SFs (Tolhurst and Thompson, 1981; De Valois et al., 1982), and provides strong support for our modeling approach. The model fit to the CSF data provided the minimum semi-saturation constant of the contrast response function near the peak of the CSF. Age group did not show differences in this optimal neuronal contrast sensitivity [t(38) = −0.19, P = 0.85, log-transformed values; Figure 5C], but instead strongly correlated with the peak contrast sensitivity of the subjects (young group: r = −0.67, P = 0.001; old group: r = −0.51, P = 0.022). On the other hand, lateral inhibition was not correlated with peak contrast sensitivity (young: r = −0.27, P = 0.24; old: r = −0.37, P = 0.11; Figure 5D).

FIGURE 5
www.frontiersin.org

Figure 5. Model fitting results. Summary plots of subjects model parameters (Iinh, σθ, cmin) obtained from tilt repulsion data (A,B) and contrast sensitivity data (C) and relations to the SF (A,B) or peak contrast sensitivity (C,D). Insets depict histograms of the corresponding variable on the ordinate with significance values for the low-SF case.

This combined modeling of orientation identification and contrast detection shows that neuronal contrast sensitivity and lateral interactions are not straightforwardly associated to the perceptual outcome and need simultaneously to be considered. Aging affected the strength of lateral interactions while individual subjects differences in neuronal contrast sensitivities still reliably represented their individual perceptual sensitivities.

Neurophysiological Data Re-Analysis

The model fitting results showed that orientation tuning widths are not systematically changed with aging, as proposed by previous behavioral reports (Betts et al., 2007; Delahunt et al., 2008; Govenlock et al., 2009), which we sought could be confirmed by orientation response properties of V1 neurons recorded from younger and older animals. Our published neurophysiological evidences from monkeys and cats showed that orientation processing in V1 is substantially affected by aging according to changes in orientation and direction bias (OB and DB) calculated by vector summation (Leventhal et al., 2003; Hua et al., 2006; Fu et al., 2010, 2013). It is difficult to directly compare orientation tuning widths from model results and neurophysiological data from these previous published analyses because these later did not consider the orientation tuning width of neurons (see Supplementary Material-5). Consequently, we re-analyzed our neurophysiological data through curve fitting, which allowed to separate effects onto the parameters of tuning: tuning width, minimum firing rate, and amplitude. This provided a clearer picture about age-related changes in orientation tuning properties.

Here, we re-analyzed one data set (Fu et al., 2013) by fitting direction/orientation tuning functions to each neuronal data set (Figures 6A,B). This data set contained 264 cells from 10 senescent monkeys and 140 cells from 4 young monkeys. We first tested whether a given neuron's data can be considered as tuned to orientation/direction of motion with an F-test (Mazurek et al., 2014; Tzvetanov, 2016). Then, we extracted three important parameters of the tuning curve from the tuned cells: the minimum firing rate in the tuning curve (r0), amplitude of firing (A), and tuning width as half-width at half-amplitude (HWHA).

FIGURE 6
www.frontiersin.org

Figure 6. Neurophysiological data re-analysis of orientation tuning. Examples of fitting results for orientation (A) and direction (B) of motion selective cells; OB, orientation bias index, DB, direction bias index. Distribution of minimum firing rate (C), amplitude (D), HWHA (E), and maximum firing rate (F) for the young and old cells. Arrows depict mean values.

Re-analysis results showed that in older monkeys, 143/264 (54.17%) neurons were found as tuned, a proportion which significantly differed (χ2 = 5.014, P < 0.05) from the 92/140 (65.71%) tuned cells in the younger animals. This confirmed the previous reports of smaller number of orientation tuned cells in V1 of senescent animals (Schmolesky et al., 2000; Leventhal et al., 2003; Hua et al., 2006; Fu et al., 2013). Importantly, r0, which here we consider as neural noise, was substantially increased [young: 32.83 ± 1.80, old: 52.11 ± 2.68; t(226.51) = 5.28, P < 0.001], but amplitude A (young: 79.47 ± 4.60; old: 77.66 ± 4.29) and tuning width HWHA (young: 35.73 ± 1.37; old: 33.85 ± 1.00) were similar [A: t(233) = −0.28, P = 0.78; HWHA: t(233) = −1.13, P = 0.26] in the two age groups of neurons (Figures 6C–E). A control analysis on a subset of the data showed that these results were not due to possible confounding factors such as mixing different single units together (see Supplementary Material-4). Additionally, the maximum firing rate of the cells (r0+A) in the two populations (young: 112.30 ± 4.86; old: 129.77 ± 4.88) differed [t(233) = 2.42, P < 0.05; Figure 6F], which was mainly attributed to the change in r0. The key results described above were also confirmed, out of the amplitude effect, by re-analyzing a second data set (Fu et al., 2010) from our laboratory about orientation/direction tuning and surround effects onto CRF (see Supplementary Material-3).

Discussion

Neurophysiological studies on orientation processing in the primary visual cortex of animals have proposed that inhibitory function generally declines with age, and thus healthy older humans might have worse early perceptual abilities than younger persons. However, various psychophysical reports on aging have shown no changes in perceptual capabilities of orientation processing or even opposite effects, where interpretation leads to stronger inhibition. Firstly, we investigated these issues through a combination of psychophysical and modeling approaches. It was found that the amount of center-surround tilt repulsion, attributed to lateral inhibitory interactions and local orientation tuning widths, in the older group was higher compared to younger observers, and that contrast sensitivity of older adults was also globally lower than younger adults. We found a common explanation of all these phenomena in a single model of V1 that dissociated lateral and local inhibitory effects on the perceptual outcome of tilt perception. Our behavioral results and computational modeling provide important evidence that low-level processing deficits in the visual system of elders could be attributed to stronger lateral inhibition. Additionally, the modeling predicted that orientation tuning width in V1 globally should not change with aging. Therefore, secondly, we re-analyzed previous physiological data published by our laboratory and found that the neuronal data followed that prediction. That is, orientation tuning widths seemed globally stable across age, while the minimum firing rate of the tuning curve had a strong increase. Thus, our psychophysical modeling and re-analysis of neurophysiological results consistently revealed a detail picture of age-related changes in orientation processing, which we think solved some contradictions between neurophysiological and behavioral reports, and uncovered a differential age-related effects on local inhibition.

Lateral inhibition between neural mechanisms tuned to different orientations in V1 was proposed as an explanation of the repulsive TI already more than 40 years ago (Blakemor et al., 1970; Blakemore et al., 1973; Georgeson, 1973). An overall increase in illusory bias of tilt observed in the elderly group suggested an increase in lateral inhibition within V1 during aging of early visual function, which was confirmed by the model results. Our results are in line with recent reports that demonstrated an increase in perceptual surround suppression of contrast, orientation and other features in the older adults when compared to younger observers (Karas and McKendrick, 2015; Nguyen and McKendrick, 2016). We reproduced these previous contrast surround suppression results with our own stimuli by conducting a center-surround contrast discrimination task for all participants in our study, which correlated to the tilt biases of the subjects (see Supplementary Material-1). Thus, the similarity in results and interpretation of these two indexes of tilt repulsion and contrast suppression support the idea that they allow to measure a common neural mechanism, and we propose that both could be interpreted as stronger lateral inhibition. This stronger lateral inhibition of older adults should be dissociated from the standard interpretation of reduced local inhibition.

The model, based on the tilt repulsion data, also predicted that the orientation tuning widths of neurons (corresponding to local inhibition; Somers et al., 1995; Shapley et al., 2003) were similar between older adults and the younger population. These findings account for psychophysical results that sensitivities for orientation processing of older adults are not systematically and substantially changed when compared to younger adults (Betts et al., 2007; Delahunt et al., 2008; Govenlock et al., 2009). In one study (Govenlock et al., 2009), using a typical orientation masking paradigm that relates to the underlying orientation tuning, the authors did not find differences in tuning widths between groups of young and old persons, as in our results. Although our raw orientation thresholds were different between the two aged groups, these raw results should not be misinterpreted as evidence of differences, due to no adjustment of these values to the individual sensitivity. Indeed, the previous reports (Delahunt et al., 2008) showed that orientation discrimination thresholds are similar across age groups once they are adjusted for inter-subject variability of detection sensitivity, and thus age, which we could not perform due to lack of measurement in our experimental design. Instead, our methods allowed us to infer the underlying theoretical population tuning widths through the modeling of the tilt repulsion effect (Tzvetanov, 2012).

The computational prediction that orientation tuning width of neurons were unchanged in older adults when compared to younger populations was confirmed by re-analyzing our previous neurophysiological data collected in young and senescent macaques. These neurophysiological data was collected from classical receptive field measures, which relates to local inhibition in our model. The results showed an increased minimum firing rate of the tuning curve and spontaneous activities, which we consider here as neural noise, and maintained orientation tuning widths with aging. Additionally, the agreement of results with a second data set leads to a consistent and deeper understanding of neuronal tuning changes and aging effects onto orientation processing (see Supplementary Material-3). These outcomes suggest that the two levels of local inhibition, an untuned inhibition affecting neural noise and a tuned inhibition influencing orientation tuning width (Somers et al., 1995; Shapley et al., 2003), were changed by aging differentially, with a weaker untuned local inhibition but a maintained tuned inhibition. These findings revealed a more complex picture of age-related effects in the local inhibitory circuits than previously thought.

At the neurophysiological level, previous studies from different laboratories reported an age-related increase in spontaneous neuronal activity in macaque V1 (Leventhal et al., 2003; Wang et al., 2005; Zhang et al., 2008), cat A17 (Hua et al., 2006; Wang et al., 2014), rat A17 (Wang et al., 2006), and rat primary auditory area A1 (Hughes et al., 2010). The higher spontaneous activities in older animals cannot be accounted for with possible technical issue of differential effects of anesthesia upon cortical function in young and old animals. This possibility has been examined previously in both NHP (Schmolesky et al., 2000; Leventhal et al., 2003; Nguyen and McKendrick, 2016) and cats (Hua et al., 2006). The properties of individual cells, in both old and young animals, were recorded while systematically varying anesthesia and paralytic levels. It was found that giving as much as four times the minimum level of general anesthesia or paralytic required to anesthetize or paralyze both old and young animals does not alter the degree of selectivity for orientation and direction exhibited by V1 cells. Strong increase in anesthesia decreased neuronal responsiveness similarly in both groups. In addition, a greater sensitivity of old animal to anesthesia is hard to reconcile with the findings that old cells exhibit higher spontaneous firing than do young cells. Thus, differences of anesthesia effects in old animals does not seem of a strong concern. One mechanism hypothesized to underlie such age-related increase in spontaneous neuronal activity was a reduction of neuronal local inhibition due to diminished GABAergic neurotransmission (Leventhal et al., 2003; Wang et al., 2005; Zhang et al., 2008). Since weaker orientation tuned inhibition should strongly broaden the tuning widths (Somers et al., 1995; Shapley et al., 2003), which we did not observe, we propose that a strong weakening of untuned inhibition with aging provides an interpretation of the previously observed changes, while the tuned inhibition remains intact.

However, lateral inhibition, which involves long-range interactions among neurons, seems to increase with aging. The weak GABA mediated local inhibition might result in much higher background firing rates. Hence, a speculative possibility is that an increase in the lateral inhibition is created to counterbalance the decline in local inhibition. However, this “chain of effects” hypothesis needs further research.

In view of previous reports showing a decrease in surround suppression in older adults during motion discrimination task (Betts et al., 2005), the finding of older adults displaying a stronger center-surround inhibition in TI supports the proposal that there is an inhibitory process in the “static” domain of V1 that is different from the mechanism of dynamic motion integration in MT. Physiological studies on primate motion processing reported that center-surround receptive fields are particularly common in MT (Born and Bradley, 2005; Anton-Erxleben et al., 2009) except the input layer IV (Raiguel et al., 1995), which is one of the evidences indicating that surround inhibition observed in motion direction discrimination paradigm involved area MT and might not be inherited from its feed-forward inputs. Some findings (Tadin, 2015) also showed that there seem to be a direct relationship between area MT/V5 and spatial suppression in motion direction discrimination paradigm. Additionally, they found that subjects' performance in motion discrimination task did not change when V1 processing was disrupted through TMS pulses. On the other hand, lateral inhibition in V1 between neural mechanisms tuned to different orientations was suggested to underlie the repulsive TI (Blakemor et al., 1970; Georgeson, 1973). Thus, these two psychophysical measures seem to estimate different aspects of cortical function and seem to reflect independent neuronal mechanisms.

Last, it is noteworthy to mention that the current study on aging brings together different methods and results obtained from different species, humans and non-human primates. It would be an important step forward to be able to demonstrate in a single observer the link between the neurophysiological changes and their putative perceptual consequences. Unfortunately, there is no study about the behavioral changes on orientation discrimination abilities of senescent NHP subjects yet. We know of one recent study demonstrating visual decline in old NHP, but animals performed a simple fixation task with biologically irrelevant stimuli (Csete et al., 2015). At our best knowledge, such a work in elder studies is missing, we think because of the vulnerability of elder animals. Currently, we are only able to compare human behavioral data and physiological analyses from other species.

In conclusion, our behavioral, computational and physiological findings provide a new perspective on aging of the visual system. We unveiled two different types of aging-related changes at the primary visual cortex, increase in lateral inhibition and unchanged local orientation coding capacities, but associated to higher neuronal noise. This suggested that local and lateral inhibition were differently affected by aging, and explained disparate results among previous behavioral and neurophysiological studies.

Materials and Methods

Subjects

The experiment included two groups: 21 younger adults (19 young adults of 20–30 years, plus co-author TT who was 41 years old) and 21 older adults (65–79 years). One young subject data was excluded due to no clear staircase convergence in the TI task (very high thresholds). One elder's data was excluded because of incorrect task performance in the replication of the contrast discrimination task (staircases without convergence properties). Subjects younger than 30 years old were students of the University of Science and Technology of China, while others were recruited from local communities. All subjects were naive to the purpose of the experiments (except the author), and their informed consents were obtained before participation. Examination (MMSE) was performed on older adults to exclude probable dementia. Alcoholism, stroke and depression were also exclusion criteria by questionnaires before conducting experiments. Participants also provided information about their general health, to exclude people with systemic conditions known to affect visual function (for example, diabetes, migraine, schizophrenia, and epilepsy) or who were taking medications known to affect visual function (e.g., anti-anxiety or anti-depressant medications). All participants were measured with normal or corrected-to-normal visual acuity ]younger = 0.016 ± 0.050, elder = 0.013 ± 0.071 (logMAR); mean ± std]. This research has been approved by the ethics committee of University of Science and Technology of China and followed the guidelines of the Declaration of Helsinki. Subjects were paid for their participation hourly.

Set-Up

The experiment was conducted in a dimly illuminated room. Stimuli were displayed on a 40.0 × 30.0 cm CRT monitor (Sony G520; 85 Hz, resolution of 1,600 × 1,200 pixels) with self-programmed Matlab functions (Mathworks Inc.) using the Psychophysics toolbox (Brainard, 1997). Stimuli were displayed using an NVIDIA Quadro K600 graphics system and viewed binocularly. To avoid local cues for vertical/horizontal and position, the screen was delimited by a 30.0 cm diameter circular window cut in a black cardboard (Tzvetanov, 2012). Luminance values were obtained with the help of the contrast box switcher (Li et al., 2003), that allowed to extend luminance range digitization above 10 bits, and thus provided the necessary minimum contrast step for CSF measure. Calibration was performed each day.

The eye-to-screen distance was maintained with a chin rest and fixed at 4 meters for CSF measure and the experiments for TI test with high-SF, and 2 m for TI test at low-SF. Experiments were initiated by subjects with a keyboard press. Subjects were requested to fixate on a small black square on the screen center and stimuli would be presented 200 ms after fixation point disappearance. Subjects responded by pressing corresponding keyboard keys.

Stimuli

Stimuli used for CSF measure were vertical sine-wave gratings with different SFs (0.71, 1, 1.41, 2, 2.83, 4, 5.66, 8, 11.31, 16, and 22.63 cycle/degree) (Figure 1A). All stimuli subtended 2° aperture, and were presented on the center of monitor with 40 cd/m2 background. A border-mask blended the stimuli to the background to minimize edge effects. There were 275 trials in total (25 trials/SF). Each trial consisted of two 100 ms intervals that separated by a 500 ms inter-stimulus interval (ISI), and each interval was announced by a prior short beep. The target sine-wave grating was presented in one of the two intervals.

Stimuli used in the tilt illusion (TI) experiments were a central Gabor patch (test grating) surrounded by an annulus of sine-wave grating (inducer) with different surround orientation (SO: 0°, ±15°, ±30°, ±75°; angle was defined with respect to the orientation of the center test) (Figure 1B). The orientation of the central Gabor patch was varied around vertical from trial to trial to measure subject's perceived verticality. The Gabor patch was defined through the following equation:

L(x, y)= Lo+LoC exp ((x2+y2)/σ2)                 ×cos(2πf(xcosθ+ysinθ))    (1)

with L0 the background luminance of the screen, C the Gabor patch contrast, and f and θ its spatial frequency and angle relative to vertical. The size of the stimuli were scaled with SF, keeping the central test stimulus fixed at 4 cycles. Surround annulus width was equal to center test diameter. Both test and inducer had same spatial frequency (SF) and constant contrast (90%) and were presented on a mean background luminance of 35 cd/m2. Specifically, TI was measured under two SF values, a low-SF located around peak contrast sensitivity, a high-SF which smaller than cutoff SF but higher than low-SF. Each participant had his/her own low-SF and high-SF values, obtained from CSF curve individually. Low-SF and high-SF conditions were conducted in two separate blocks and each block consisted of 420 trials (60 trials/SO). The stimulus in each trial was presented for nine video frames (~100 ms) and no feedback was provided.

Procedure

Before formal measure, each subject received a short practice session (about 80 trials) for each task. A CSF measure was conducted prior to TI measure to extract low-SF and high-SF values. TI measure for high- and low-SF were conducted sequentially in two separate blocks.

A 3-key response design (Kaernbach, 2001) was used for CSF measure, with two keys for responding in which interval the subjects perceived the signal and a third key for undecided decision. The contrast of the target sine-wave grating was varied in multiple of 10% according to the up-down staircase procedures, with correct responses followed by decreasing contrast with 4.5/1.5 steps and incorrect responses by increasing contrast of 8/7 steps for SFs of 2n (with n from ½ to 4½ in steps of 1, and from 1 to 4 in steps of one respectively). Undecided key presses were randomly drawn as correct/incorrect. Starting points were contrasts of 0.5,192 0.005, 0.5, 0.005, 0.5, 0.005, 0.5, 0.005, 0.7 0.05, 0.8 for the successive 11 SFs, respectively, and each staircase “down” step-size was additionally 3 times bigger for the first 4 trials (example results in Figure 1C). Feedback was provided by different sound for correct and incorrect responses and mute for the undecided key.

The 3-key design was also used for TI measure where the observer was required to report whether the orientation of the center test was clockwise (CW) or counter clockwise(CCW) from his/her internal vertical standard by pressing two predefined keys. The third key was allowed if the subject felt he/she did not perceive central target grating (due to surround suppression; especially at high SFs), and “not seen” cases were randomly drawn as CW or CCW for the staircase procedure. The center orientation was varied according to the weighted up-down rule with steps up/down of 5/2 and 2/5 degrees with base step of 1°. Staircases started at the opposite side of the convergence point allowing rapid measures within the transition region of the psychometric function. No feedback for any responses.

Data Analysis

We used Bayesian fitting to adjust theoretical psychometric functions to the CSF and TI data (Treutwein and Strasburger, 1999).

Contrast Sensitivity Function

A 2D psychometric function was fit to the 2D contrast-SF (c, f) data, with the probability of correct response defined as:

P(c, f)=γ+​ 1γλ1+exp(log(214)(log(c)log(1S(f)))/σ    (2)

with parameters γ and λ being subject's “guess rate” (see below) and lapse rate, and 2σ defining the spread between 16 and 84% of the function in the range γ to 1-λ (assuming constant spread, σ, across SFs). S(f) is the standard 3-parameters sensitivity function (Rohaly and Owsley, 1993):

S(f)=Mfaexp(-fb)    (3)

used in previous studies to define the CSF shape in the SF dimension. The 3-key design data was processed following Garcia-Perez (García-Pérez, 2010), which in the event the subject followed the 3rd key instructions allows a decrease in measurement variability. The Likelihood function is then:

logML=(iyilog(p(ci,fi))+(1yi)log(1p(ci,fi)))              +(k1/2log(p(ck,fk))+1/2log(1p(ck,fk)))    (4)

with the first sum (trials with index i) running over all responses 1st or 2nd interval and the second sum (trials with index k) running over the “undecided” 3rd key responses. The variables in the above equation are: (ci, k, fi,k)–contrast & SF pair presented at trial i or k, yi–correct/incorrect (1/0) response of subject at contrast & SF levels (ci, fi), p (ci, k, fi, k)-theoretical probability of correct response (psychometric function). In this ML equation, the first term is the standard log-ML term for fitting binomial data; the second “Fechner” term is simply log (p(1-p)), the logarithm of the binomial variability at the stimulus levels for which the subject pressed “undecided”/“not seen”; it is maximized when p(ck, fk) is 0.5, i.e., when the subject is totally uncertain about the interval of signal presentation, and thus provides a firm theoretical ground for introducing “undecided”/“not seen” responses into the 2AFC design. The lapse rate was fixed at 1% for all but one subject where it was zero. The “guess rate” γ was 0.5. An example CSF fit for one subject and his responses to each presented stimulus is displayed in Figure 1C.

Tilt repulsion. We fit a 1D psychometric function to the orientation discrimination data for each surround orientation, with probability of CW responses to target orientation θ given by:

P(θ)=λ+12λ1+exp(log(214)(θa)/σ    (5)

where λ is subject's lapse rate, and a and σ being the perceived vertical orientation (also called “bias”) for the given surround and the threshold of the subject for perceiving a deviation from verticality, respectively. Because of the symmetry in the experimental design (symmetric surrounds of ±15, ±30, and ±75°), for the fitting we imposed that thresholds of opposite surround orientations (e.g., −30 and +30°) are the same. The lapse rate was fixed at 0% for low-SF condition, and 1% for high-SF. The data was processed by eliminating any datum with 3rd key responses (subject did not see the target), and we computed the amount of surround suppression as the proportion of 3rd key presses. Bias values were computed as the half-difference between two opposite surround orientations.

Statistical Analysis

For contrast sensitivity, a two factors between-within ANOVA (age group × spatial frequency) was conducted. For TI, a 3-way between-within subject ANOVA (age group × surround orientation × spatial frequency) was conducted on all data. For center-surround contrast discrimination task, a t-test was used to analyze the suppression ratio data (paired or unpaired, as appropriate). All statistical levels use Geisser-Greenhouse epsilon-hat adjusted values, where appropriate. Spearman rank correlation was used. Data were expressed as mean ± SE.

Model

The model was developed for another study in our laboratory, concerning Amblyopia deficits (Huang et al., 2017), and a detailed description was given in the accompanying manuscript (Huang et al., 2017). We provide the model description as in the original work for consistent methodological description in the current work.

Simple Model of V1 Surround-to-Centre Interactions

We assume, as in many previous studies, that simple feature perception as local orientation and contrast can be explained through the decoding of primary visual cortex neuronal activities. Therefore, we investigated a simple V1 model of two-layer neurons coding the main features of interest in the study: orientation, contrast, SF, and space. The model consists of orientation hyper-columns arranged into a hexagonal structure, with each hyper-column containing neurons responding to various contrasts and SFs. First layer neurons can be thought of simple cells whose responses are as follows:

rijk(θ,f,c)=A×T(θ;θi)×F(f;fj)×C(c;ck)    (6)

with “preferred” features (θi,fj,ck) and the three normalized tuning functions to orientation, SF and contrast are described as wrapped-Gaussian (Swindale, 1998), log-Gaussian (Yuan et al., 2014), and hyperbolic ratio (Albrecht and Hamilton, 1982), respectively (A is the maximum amplitude of firing of the neuron). They are:

T(θ;θi)=e-12(θ-θiσθ)2    (7)
F(f;fj)=e-12(f-fjσSF)2    (8)
C(c;ck)=cnckn+cn    (9)

Remark: for the contrast tuning, ck is the semi-saturation constant and can be called the “preferred” contrast of the neuron, since for contrasts around ck the neuron is the most informative above the input contrast (Chirimuuta and Tolhurst, 2005; May and Solomon, 2015a,b), and away from ck it asymptotes and provides no information about contrast input.

These simple cells feed the second layer of neurons through a spatial (excitatory center)–(inhibitory surround) connectivity structure, whose responses Rijk follow the conductance-based model (Grossberg, 1988; Piëch et al., 2013):

Rijk(θ,f,c)=h(vijk)    (10)
τdvijkdt=-vijk+(ve-vijk)ge+(vi-vijk)gi    (11)
ge=mnoωijk,mnocenrmno    (12)
gi=mno,(x,y)ωijk,mno,(x,y)surrmno,(x,y)    (13)

with h() a transducer (rectifying) function transforming voltage to firing rate and feature weights ω's defined as:

h(v)=m×max(0,v-T)    (14)
ωijk,mnocen=IcGimcGjnc    (15)
ωijk,mno,(x,y)sur=IsGimsGjnsGx,y    (16)

and the various parameters are: T is the voltage threshold of firing, m is the slope of voltage-to-firing rate relation, τ is the cell time constant, ve and vi are the excitatory and inhibitory equilibrium voltage potentials, ge and gi are respectively the excitatory and inhibitory conductances feeding the corresponding neuron through a weighted sum of first layer activities (ge sum is within hyper-column; gi sum is over all surrounding hyper-columns), Gim, jn are Gaussian tuned feature weights [respectively within orientation and within SF; with possible different tuning widths indexed (c,s)], Gx, y is a spatial weight function summing surrounding hyper-columns activity, and Ic, s are the center/surround excitatory/inhibitory input strengths, respectively. Here, it is assumed that the weights are independent across features and iso-feature tuned (peaking at the receiving neuron preferred value (i,j), i.e., iso-orientation and iso-SF).

In the feed-forward model equation (11) can be analytically solved, giving:

vijk=vege+vigi1+ge+gi    (17)
Rijk=h(vijk)    (18)

Using all relations above and an input with uniform surround (all surrounding hyper-columns are stimulated with the same stimulus of orientation θs and contrast cs), assuming the center stimulus (θc,cc) excites the center hyper-column, the center input excitatory/inhibitory conductances can be analytically computed:

ge,i=Ic,sAc,sσθc,sσSFc,sθc,s BC,S×e-12(θi-θc,s)2(σθc,s)2+σθ2×e-12(f-fj)2(σSFc,s)2+σSF2    (19)

The various constants are defined as:

BC,S=((SFc,s)2(log(2))22+Kc,slog(2))SFc,se    (20)
(θc,s)-2=(σθc,s)-2+(σθ)-2    (21)
(SFc,s)-2=(σSFc,s)-2+(σSF)-2    (22)
Kc,s(f;fj)=(σSFc,s)2log f+σSF2log fj(σSFc,s)2+σSF2    (23)

where σθc,s and σSFc,s are the orientation and SF tuning widths of the weight functions Gim and Gjn for center-center and surround-to-center connections, respectively, Ac, s are the contrast-weighted amplitude of firing of the input neurons rijk for center/surround respectively, Ic, s are the excitatory and inhibitory weight amplitudes, and

Is=ns×Iinh    (24)

defines the total inhibitory input from all surrounding hyper-columns with mean inhibitory strength per hyper-column Iinh, respectively; ns is a “mean” number of surround hyper-columns influencing the center.

Fixed Parameters

Summary of the set of parameters that were fixed in the model (following a “normalized” conductance based subtractive inhibition model, e.g., Piëch et al., 2013): ns = 6, n = 2, A = 2, σSF = 1, m = 1, T = 1, Ic = 1, ve = 14/3, v I = −2/3. SF dimension sampling was every ¼ octaves from ½ up to 64 c/d, orientation feature sampling every 2°. The contrast tuning relation was kept normalized by multiplying its amplitude by a factor (1 + ckn).

Model Based Fitting of CSF and Tilt Perception Data

Here we investigate how based on the output activity Rijk of the network we can predict the perception of the subject in the two main features of interest, contrast detection for predicting the contrast sensitivity function and orientation identification for predicting the tilt repulsion effect. It is assumed that perception is based on decoding of the center hyper-column activities as described below.

Modeling the Contrast Sensitivity Function (CSF in 2D)

In this experiment, the target stimulus is a vertical and uniform sine-wave grating limited in a circular spatial window, whose strength (contrast) is varied in order to measure the perception threshold across all SFs.

Given the uniform input stimulus, all neurons rijk have exactly the same input and thus their activation across orientation, SF, and contrast have the same profiles and peaks (for any i,j,k: rijkc=rijks). Second, given the task of detecting always a vertically oriented grating, by assuming that subjects disregard other orientated activities through an unspecified attentional mechanism, we simplified the term over orientations in equation (19) into a value of one (in practice, this simplification can be thought of pooling these orientation neuronal activities into the constants Ic, s or Ac, s) and modeled only one orientation of network activities.

In the above model description, one important function in predicting the contrast sensitivity at a given SF is the hyperbolic ratio of the neuronal population (Chirimuuta and Tolhurst, 2005; May and Solomon, 2015a,b). We assume that contrast detection across all SFs is performed by decoding the activities of the neurons with the smallest semi-saturation constant ckmin(SF)=min(ck|SF). In our model, to predict the CSF across SFs we additionally need to properly describe the relation between ckmin and SF together with SF tuning width versus preferred SF. Based on previous neurophysiological reports (see Supplementary Material-2), we fixed:

ckmin(fj)=cmin×(ab)ae-afjae-fjb    (25)
σSF=e-fj-20.2    (26)

where fj is the preferred SF of the neuron and parameters (cmin,a,b) define the neuronal sensitivity function across SFs. We fixed σSFc,s=σSF/2 because of center/surround tuning widths entanglement in predicting behavioral results, Ic = 1; ns = 6, and neurons with ck > 1,000 were pruned.

Last step, for predicting the contrast sensitivity function across SFs, we used the standard signal detection theory, defining the psychometric function:

P(c,f)=γ+(1-γ-λ)Pth    (27)
Pth=2(0+RjR0Var(Rj)+Var(R0)dRj0.5)    (28)

with the lapse and guess rates obtained as described in the above section Data analyses, and R0 being the activity of the neurons for no signal input.

From the above model description, there are only four free model parameters that need to be adjusted for predicting the full CSF: cmin, Iinh, a, b. It was done by replacing equation (2) with equation (27) and following all remaining steps.

Modeling Orientation Identification (Tilt Illusion)

For this feature, we make a different set of simplifications in the model feature space based on the experimental design for tilt perception. Here, center and surround hyper-columns are stimulated with varying orientations while the contrast of the center and surround stimuli are kept constant and equal. Therefore, we can describe the two-layer neuronal network activities through the above mathematical development. But here, we fixed the input layer contrast activation at A = 2 and c = 1 (near maximum contrast, e.g., Piëch et al., 2013), and given the task of identifying the orientation of the center stimulus for a fixed SF (Georgeson, 1973), we simplified into equation (19) the term over SFs by assuming subjects disregard other SF neuronal activities (eg through an unspecified attentional mechanism), giving:

ge,i=Ic,sAc,sσθc,sθc,s(e-12(θi-θc,s)2(σθc,s)2+σθ2)    (29)

with all constants defined in model description section. We fixed the center and surround feed-forward orientation tuning widths to σθc,s=σθ /2.

Last, to predict the orientation psychometric function (tilt data), we decoded the perceived orientation of the stimulus [a in Equation (5)] as the vector average orientation of the center hyper-column activities (Gilbert and Wiesel, 1990; Seung and Sompolinsky, 1993; Tzvetanov, 2012). For this feature, there are only two free parameters, Iinh and σθ, that are sufficient to provide tilt misperception description (Gilbert and Wiesel, 1990; Tzvetanov, 2012). Lapse rate was fixed at 1%.

General Fitting Procedure

First, each individual data set was fit with the tilt prediction part only for SO of ±15 and ±30° in order to extract Iinh at both high and low SFs. For a given SF, a single discrimination threshold for all SOs was set as free parameter. Then, the parameter Iinh was used as fixed in the CSF fitting procedure in order to find the best three parameters (cmin, a, b) that described subject's CSF data.

Neurophysiological Data Re-Analysis

The direction of motion tuning curves were fit with a wrapped von-Mises function of the form:

r(Δθ1)=b+a1exp(k(cos(2Δθ1))-1)r(Δθ2)=b+(a1-a2)exp(-2k)               +a2exp(k(cos(2Δθ2)-1))    (30)

where r is the response of the neuron, and the parameters b, a1, a2, and k allow to compute the physiological parameters of interest, minimum firing rate, the two amplitudes to opposite direction of motions, and half-width at half-amplitude, as:

                        r0=b+a1exp(-2k)            A(θpref)=a1(1-exp(-2k))A(θpref+180)=a2(1-exp(-2k))               HWHA=0.5cos1(1+log1+exp(2k)2)/k)    (31)

The two variables Δθ are defined as Δθ1 = θ-θpref and Δθ2 = θ-θpref-180, each defined in the range ±90°, and the additional term for r(Δθ2) allows continuity at the boundary. The fitting was done by minimizing the χ2 between the data and the model function, and because standard errors were not present, it was assumed sei = (FANO × yi/5)½ with FANO = 1.5, and for the few cases where yi = 0, sei = 0.5477. We performed three successive fits: (1) an orientation tuned function fixing a1 = a2 (4 free parameters: b, a, k, θpref); (2) a direction tuned function leaving parameters a1 and a2 independent; and (3) a flat-topped version of the equation (Albrecht and Hamilton, 1982) allowing for more peaked or flatter types of tuning curves, that included one more parameter ν [in equation (30), Δθ1,2 is replaced by Δθ1,2/2 + ν sin(Δθ1,2)]. Multiple refits were done with random starting points for finding the best set of parameters. Each fitted curve was then used to test, with a nested F-test (Tzvetanov, 2016), whether it described the data better than the global mean of the data, i.e., whether the cell could be described as tuned. Then, if two or more functions described the data better than the mean, we further used a nested F-test to test whether the more complex functions (with additional parameters) described the data better than the simpler one, and selected the corresponding model. In the models with different amplitudes to opposite motion direction, the amplitude of firing was defined as the higher of the two. In the flat-top model, the half width of half amplitude (HWHA) was directly searched in a discretized direction space (0.1° steps). At the end, we additionally eliminated any fitted curve with too narrow of tuning (HWHA < 15), which discarded 8 old cells and 7 young cells. Two example fits are presented on Figures 6A,B for one orientation and one direction tuned cells, together with their orientation and direction bias values. Here, we were not interested of the overall goodness-of-fit of the model to the data. Nevertheless, in Supplementary Material-4, we show that there are no essential differences in fit quality between old and young cell types once variability is properly accounted for.

Neurophysiological Experimental Materials and Methods

Please see Fu et al. (2013) for detailed information regarding electrophysiogical experiments in monkeys. Briefly:

Animals

Subjects for this study were 14 male rhesus monkeys (Macaca mulatta) classified into 2 groups: the young (YA) group consisted of 4 monkeys who were 5.8 ± 0.5 (mean ± standard deviation) years old; the old group included early senescent (ES) group (7 monkeys at ages of 21.6 ± 1.5 years) and late senescent (LS) group (3 monkeys at ages of 29.3 ± 2.3 years). Monkeys were examined ophthalmoscopically and had no indication of optical or retinal problems that would impair visual function. Retinal blood vessels, lens clarity, and the maculae were all within the normal range.

Animal Preparation and Recording

Subjects were sedated with Ketamine HCl (10–15 mg/kg, i.m., ParkeDavis, Morris Plains, NJ, USA) and then anesthetized with 3–5% halothane (Halocarton Laboratories, River Edge, NJ, USA) in a 70:30 mixture of N2O: O2. Intravenous and tracheal cannulae were inserted. Animals were then placed in a stereotaxic frame. All pressure points and incisions were infiltrated with a long-acting local anesthetic (2% lidocaine HCl; Copley Pharmaceuticals, Canton, MA, USA). A craniotomy was performed above the area in V1 corresponding to the representation of central vision. After surgical procedures were finished and sufficient level of anesthesia was assured, a mixture of d-tubocurarine (0.4 mg/kg/h; Sigma, St. Louis, MO, USA) and gallamine triethiodide (7 mg/kg/h; Sigma) was infused intravenously to induce and maintain paralysis. During recording, animals were artificially ventilated. End tidal CO2 partial pressure was monitored and maintained at ~4–4.5%. Body temperature was maintained at 38–38.5°C by an automatically regulated heating pad. The anesthesia was maintained with a mixture of N2O (70%), O2 (30%), and halothane (0.25–1.0%, as needed). Electrocardiography was continuously monitored. The level of anesthesia was assessed by periodically applying nociceptive stimulation (toe pinch) and monitoring the changes in the heart rate.

The nictitating membrane was prepared with neosynephrine. The pupils were dilated with atropine, and the eyes were covered with contact lenses for protection from desiccation. The optics of the eyes and retinal vasculature were monitored throughout the experiment. No visible deterioration in optics occurred during the experimental period in any of the animals. Action potentials of isolated units were recorded extracellularly by glass microelectrode (filled with 2-M NaCl, 1–3-MΩ impedance at 1000 Hz), which was advanced by a hydraulic microdrive (David Kopf Instruments, Tojunga, CA, USA).

Visual Stimulation

For each isolated single unit, the eye preference was determined and all subsequent stimuli were presented monocularly to the dominant eye. Cell's receptive field was carefully plotted on a tangent screen by hand with the use of an ophthalmoscope. After that, a 17″ CRT color monitor (85-Hz refresh rate; Sony, Tokyo, Japan) was positioned at the distance of 57 cm in front of the animal's eyes and centered on the receptive field of the cell. The program to generate the visual stimuli was written in MATLAB (MathWorks, Natick, MA, USA), using the extensions provided by the high-level Psychophysics Toolbox (Brainard, 1997) and low-level Video Toolbox (Pelli, 1997). To quantify the orientation and direction selectivity of V1 cell, drifting bars were used, whose width, length, and moving speed were adjusted to elicit strongest response from the recorded cell. The direction of motion of each presented bar was orthogonal to its orientation. We used moving bars at 24 randomly chosen movement directions, ranging from 0° to 360° in steps of 15° to compile the tuning curves for the cells studied. Each direction was presented 10 times. The inter-trial interval was 2–5 s. The luminance of the stimuli used was 39 cd/m2 for white and 0.95 cd/m2 for black.

Data Collection

Signals recorded from the microelectrode were amplified (1,000×), band-pass filtered (300–10,000 Hz), and then digitized (sampling frequency of 20 kHz) by using an acquisition board (National Instruments, Austin, TX, USA) controlled by IGOR software (WaveMetrics, Portland, OR, USA). Such original signals were stored in a computer for offline analysis. The responses to moving bars were defined as the maximal value in the peristimulus time histogram (PSTH, bin width of 10 ms) during the stimulation period. Before each drifting bar was presented, the spontaneous (baseline) activity was recorded during a period of 0.5–0.7 s while the screen with average luminance was presented.

Ethics Statement

This study was carried out in accordance with the recommendations of the ethics committee of University of Science and Technology of China with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the the ethics committee of University of Science and Technology of China.

Authors Contributors

ZW and TT designed and analyzed behavioral experiments; ZW performed experiments; ZW and TT performed data analysis and neurophysiological re-analysis; SY and YF provided physiological data and discussed its re-analysis; TT performed modeling; ZW, YZ, and TT provided project supervision and funds; ZW and TT wrote the paper; all authors discussed and commented on the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (31230032 to YZ), General Financial Grant from the China Postdoctoral Science Foundation (2015M571940 to ZW), the Fundamental Research Funds for the Central Universities (to TT), the Hundred-Talent Program of the Chinese Academy of Sciences (to SY) and the K. C. Wong Magna Fund in Ningbo University (to ZW).

Conflict of Interest Statement

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/fnagi.2018.00014/full#supplementary-material

References

Albrecht, D. G., and Hamilton, D. B. (1982). Striate cortex of monkey and cat: contrast response function. J. Neurophysiol. 48, 217–237. doi: 10.1152/jn.1982.48.1.217

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersen, G. J. (2012). Aging and vision: changes in function and performance from optics to perception. Wiley Interdiscip. Rev. Cogn. Sci. 3, 403–410. doi: 10.1002/wcs.1167

PubMed Abstract | CrossRef Full Text | Google Scholar

Angelucci, A., and Bressloff, P. C. (2006). Contribution of feedforward, lateral and feedback connections to the classical receptive field center and extra-classical receptive field surround of primate V1 neurons. Prog. Brain Res. 154, 93–120. doi: 10.1016/S0079-6123(06)54005-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Angelucci, A., Levitt, J. B., Walton, E. J., Hupe, J. M., Bullier, J., and Lund, J. S. (2002). Circuits for local and global signal integration in primary visual cortex. J. Neurosci. 22, 8633–8646.

PubMed Abstract | Google Scholar

Anton-Erxleben, K., Stephan, V. M., and Treue, S. (2009). Attention reshapes center-surround receptive field structure in macaque cortical area MT. Cereb. Cortex 19, 2466–2478. doi: 10.1093/cercor/bhp002

PubMed Abstract | CrossRef Full Text | Google Scholar

Betts, L. R., Sekuler, A. B., and Bennett, P. J. (2007). The effects of aging on orientation discrimination. Vision Res. 47, 1769–1780. doi: 10.1016/j.visres.2007.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Betts, L. R., Taylor, C. P., Sekuler, A. B., and Bennett, P. J. (2005). Aging reduces center-surround antagonism in visual motion processing. Neuron 45, 361–366. doi: 10.1016/j.neuron.2004.12.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Blakemor, C., Carpente, R. H., and Georgeso, M. A. (1970). Lateral inhibition between orientation detectors in human visual system. Nature 228, 37–39. doi: 10.1038/228037a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Blakemore, C., Muncey, J. P., and Ridley, R. M. (1973). Stimulus specificity in the human visual system. Vision Res. 13, 1915–1931. doi: 10.1016/0042-6989(73)90063-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Born, R. T., and Bradley, D. C. (2005). Structure and function of visual area MT. Annu. Rev. Neurosci. 28, 157–189. doi: 10.1146/annurev.neuro.26.041002.131052

PubMed Abstract | CrossRef Full Text | Google Scholar

Brainard, D. H. (1997). The psychophysics toolbox. Spat. Vis. 10, 433–436. doi: 10.1163/156856897X00357

PubMed Abstract | CrossRef Full Text | Google Scholar

Carandini, M., and Ringach, D. L. (1997). Predictions of a recurrent model of orientation selectivity. Vision Res. 37, 3061–3071. doi: 10.1016/S0042-6989(97)00100-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Chirimuuta, M., and Tolhurst, D. J. (2005). Does a Bayesian model of V1 contrast coding offer a neurophysiological account of human contrast discrimination? Vis. Res. 45, 2943–2959. doi: 10.1016/j.visres.2005.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Csete, G., Bognár, A., Csibri, P., Kaposvári, P., and Sáry, G. (2015). Aging alters visual processing of objects and shapes in inferotemporal cortex in monkeys. Brain Res. Bull. 110, 76–83. doi: 10.1016/j.brainresbull.2014.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Delahunt, P. B., Hardy, J. L., and Werner, J. S. (2008). The effect of senescence on orientation discrimination and mechanism tuning. J. Vis. 8, 5 1–9. doi: 10.1167/8.3.5

PubMed Abstract | CrossRef Full Text | Google Scholar

De Valois, R. L., Albrecht, D. G., and Thorell, L. G. (1982). Spatial frequency selectivity of cells in macaque visual cortex. Vision Res. 22, 545–559. doi: 10.1016/0042-6989(82)90113-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y., Wang, X. S., Wang, Y. C., Zhang, J., Liang, Z., Zhou, Y. F., et al. (2010). The effects of aging on the strength of surround suppression of receptive field of V1 cells in monkeys. Neuroscience 169, 874–881. doi: 10.1016/j.neuroscience.2010.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y., Yu, S., Ma, Y., Wang, Y., and Zhou, Y. (2013). Functional degradation of the primary visual cortex during early senescence in rhesus monkeys. Cereb. Cortex 23, 2923–2931. doi: 10.1093/cercor/bhs282

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Pérez, M. A. (2010). Denoising forced-choice detection data. Br. J. Math. Stat. Psychol. 63, 75–100. doi: 10.1348/000711009X424057

PubMed Abstract | CrossRef Full Text | Google Scholar

Georgeson, M. A. (1973). Spatial frequency selectivity of a visual tilt illusion. Nature 245, 43–45. doi: 10.1038/245043a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilbert, C. D. (1992). Horizontal integration and cortical dynamics. Neuron 9, 1–13. doi: 10.1016/0896-6273(92)90215-Y

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilbert, C. D., and Wiesel, T. N. (1990). The influence of contextual stimuli on the orientation selectivity of cells in primary visual cortex of the cat. Vision Res. 30, 1689–1701. doi: 10.1016/0042-6989(90)90153-C

PubMed Abstract | CrossRef Full Text | Google Scholar

Govenlock, S. W., Taylor, C. P., Sekuler, A. B., and Bennett, P. J. (2009). The effect of aging on the orientational selectivity of the human visual system. Vision Res. 49, 164–172. doi: 10.1016/j.visres.2008.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossberg, S. (1988). Nonlinear neural networks: principles, mechanisms, and architectures. Neural Networks 1, 17–61. doi: 10.1016/0893-6080(88)90021-4

CrossRef Full Text | Google Scholar

Hua, T., Li, X., He, L., Zhou, Y., Wang, Y., and Leventhal, A. G. (2006). Functional degradation of visual cortical cells in old cats. Neurobiol. Aging 27, 155–162. doi: 10.1016/j.neurobiolaging.2004.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Zhou, Y., Liu, C., Liu, Z., Luan, C., and Tzvetanov, T. (2017). The neural basis of spatial vision losses in the dysfunctional visual system. Sci Rep. 7:11376.

PubMed Abstract | Google Scholar

Hughes, L. F., Turner, J. G., Parrish, J. L., and Caspary, D. M. (2010). Processing of broadband stimuli across A1 layers in young and aged rats. Hear. Res. 264, 79–85. doi: 10.1016/j.heares.2009.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaernbach, C. (2001). Adaptive threshold estimation with unforced-choice tasks. Percept. Psychophys. 63, 1377–1388. doi: 10.3758/BF03194549

PubMed Abstract | CrossRef Full Text | Google Scholar

Kapadia, M. K., Westheimer, G., and Gilbert, C. D. (2000). Spatial distribution of contextual interactions in primary visual cortex and in visual perception. J. Neurophysiol. 84, 2048–2062. doi: 10.1152/jn.2000.84.4.2048

PubMed Abstract | CrossRef Full Text | Google Scholar

Karas, R., and McKendrick, A. M. (2015). Contrast and stimulus duration dependence of perceptual surround suppression in older adults. Vision Res. 110(Pt A), 7–14. doi: 10.1016/j.visres.2015.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Leventhal, A. G., Wang, Y., Pu, M., Zhou, Y., and Ma, Y. (2003). GABA and its agonists improved visual cortical function in senescent monkeys. Science 300, 812–815. doi: 10.1126/science.1082874

PubMed Abstract | CrossRef Full Text | Google Scholar

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. doi: 10.1016/S0165-0270(03)00174-2

PubMed Abstract | CrossRef Full Text | Google Scholar

May, K. A., and Solomon, J. A. (2015a). Connecting psychophysical performance to neuronal response properties I: discrimination of suprathreshold stimuli. J. Vis. 15:8. doi: 10.1167/15.6.8

PubMed Abstract | CrossRef Full Text | Google Scholar

May, K. A., and Solomon, J. A. (2015b). Connecting psychophysical performance to neuronal response properties II: contrast decoding and detection. J. Vis. 15:9. doi: 10.1167/15.6.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazurek, M., Kager, M., and Van Hooser, S. D. (2014). Robust quantification of orientation selectivity and direction selectivity. Front. Neural Circuits 8:92. doi: 10.3389/fncir.2014.00092

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, B. N., and McKendrick, A. M. (2016). Visual contextual effects of orientation, contrast, flicker, and luminance: all are affected by normal aging. Front. Aging Neurosci. 8:79. doi: 10.3389/fnagi.2016.00079

PubMed Abstract | CrossRef Full Text | Google Scholar

Owsley, C. (2011). Aging and vision. Vision Res. 51, 1610–1622. doi: 10.1016/j.visres.2010.10.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelli, D. G. (1997). The videotoolbox software for visual psychophysics: transforming numbers into movies. Spat. Vis. 10, 437–442. doi: 10.1163/156856897X00366

PubMed Abstract | CrossRef Full Text | Google Scholar

Piëch, V., Li, W., Reeke, G. N., and Gilbert, C. D. (2013). Network model of top-down influences on local gain and contextual interactions in visual cortex. Proc. Natl. Acad. Sci. U.S.A. 110, E4108–E4117. doi: 10.1073/pnas.1317019110

PubMed Abstract | CrossRef Full Text | Google Scholar

Raiguel, S., Van Hulle, M. M., Xiao, D. K., Marcar, V. L., and Orban, G. A. (1995). Shape and spatial distribution of receptive fields and antagonistic motion surrounds in the middle temporal area (V5) of the macaque. Eur. J. Neurosci. 7, 2064–2082. doi: 10.1111/j.1460-9568.1995.tb00629.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rohaly, A. M., and Owsley, C. (1993). Modeling the contrast-sensitivity functions of older adults. J. Opt. Soc. Am. A 10, 1591–1599. doi: 10.1364/JOSAA.10.001591

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmolesky, M. T., Wang, Y., Pu, M., and Leventhal, A. G. (2000). Degradation of stimulus selectivity of visual cortical cells in senescent rhesus monkeys. Nat. Neurosci. 3, 384–390. doi: 10.1038/73957

PubMed Abstract | CrossRef Full Text | Google Scholar

Seung, H. S., and Sompolinsky, H. (1993). Simple models for reading neuronal population codes. Proc. Natl. Acad. Sci. U.S.A. 90, 10749–10753. doi: 10.1073/pnas.90.22.10749

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S., and Wenderoth, P. (1999). Large repulsion, but not attraction, tilt illusions occur when stimulus parameters selectively favour either transient (M-like) or sustained (P-like) mechanisms. Vision Res. 39, 4113–4121. doi: 10.1016/S0042-6989(99)00118-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Somers, D. C., Nelson, S. B., and Sur, M. (1995). An emergent model of orientation selectivity in cat visual cortical simple cells. J. Neurosci. 15, 5448–5465.

PubMed Abstract | Google Scholar

Spillmann, L., and Werner, J. S. (1996). Long-range interactions in visual perception. Trends Neurosci. 19, 428–434. doi: 10.1016/0166-2236(96)10038-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Swindale, N. V. (1998). Orientation tuning curves: empirical description and estimation of parameters. Biol. Cybern. 78, 45–56. doi: 10.1007/s004220050411

PubMed Abstract | CrossRef Full Text | Google Scholar

Tadin, D. (2015). Suppressive mechanisms in visual motion processing: from perception to intelligence. Vision Res. 115(Pt A), 58–70. doi: 10.1016/j.visres.2015.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Tolhurst, D. J., and Thompson, I. D. (1981). On the variety of spatial frequency selectivities shown by neurons in area 17 of the cat. Proc. R. Soc. Lond. B. Biol. Sci. 213, 183–199. doi: 10.1098/rspb.1981.0061

PubMed Abstract | CrossRef Full Text | Google Scholar

Treutwein, B., and Strasburger, H. (1999). Fitting the psychometric function. Percept. Psychophys. 61, 87–106. doi: 10.3758/BF03211951

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzvetanov, T. (2012). A single theoretical framework for circular features processing in humans: orientation and direction of motion compared. Front. Comput. Neurosci. 6:28. doi: 10.3389/fncom.2012.00028

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzvetanov, T. (2016). Commentary: robust quantification of orientation selectivity and direction selectivity. Front. Neural Circuits 10:25. doi: 10.3389/fncir.2016.00025

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzvetanov, T., and Womelsdorf, T. (2008). Predicting human perceptual decisions by decoding neuronal information profiles. Biol. Cybern. 98, 397–411. doi: 10.1007/s00422-008-0226-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Xie, X., Li, X., Chen, B., and Zhou, Y. (2006). Functional degradation of visual cortical cells in aged rats. Brain Res. 1122, 93–98. doi: 10.1016/j.brainres.2006.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhou, Y., Ma, Y., and Leventhal, A. G. (2005). Degradation of signal timing in cortical areas V1 and V2 of senescent monkeys. Cereb. Cortex 15, 403–408. doi: 10.1093/cercor/bhh143

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Yao, Z., Yuan, N., Liang, Z., Li, G., and Zhou, Y. (2014). Declined contrast sensitivity of neurons along the visual pathway in aging cats. Front. Aging Neurosci. 6:163. doi: 10.3389/fnagi.2014.00163

PubMed Abstract | CrossRef Full Text | Google Scholar

Wenderoth, P., and Johnstone, S. (1988). The Different mechanisms of the direct and indirect tilt illusions. Vision Res. 28, 301–312. doi: 10.1016/0042-6989(88)90158-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wenderoth, P., and Smith, S. (1999). Neural substrates of the tilt illusion. Aust. N. Z. J. Ophthalmol. 27, 271–274. doi: 10.1046/j.1440-1606.1999.00191.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, N., Liang, Z., Yang, Y., Li, G., and Zhou, Y. (2014). Changes of spatial and temporal frequency tuning properties of neurons in the middle temporal area of aged rhesus monkeys. Eur. J. Neurosci. 40, 2652–2661. doi: 10.1111/ejn.12634

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Wang, X., Wang, Y., Fu, Y., Liang, Z., Ma, Y., et al. (2008). Spatial and temporal sensitivity degradation of primary visual cortical cells in senescent rhesus monkeys. Eur. J. Neurosci. 28, 201–207. doi: 10.1111/j.1460-9568.2008.06300.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: behavioral measure, computational modeling, neurophysiology, early visual processing, senescence

Citation: Wang Z, Yu S, Fu Y, Tzvetanov T and Zhou Y (2018) Aging Potentiates Lateral but Not Local Inhibition of Orientation Processing in Primary Visual Cortex. Front. Aging Neurosci. 10:14. doi: 10.3389/fnagi.2018.00014

Received: 27 May 2017; Accepted: 12 January 2018;
Published: 05 February 2018.

Edited by:

Ying Xu, University at Buffalo, United States

Reviewed by:

Guilan Shi, Old Dominion University, United States
Tun Li, University of California, San Francisco, United States

Copyright © 2018 Wang, Yu, Fu, Tzvetanov and Zhou. 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 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: Tzvetomir Tzvetanov, tzvetan@ustc.edu.cn
Yifeng Zhou, zhouy@ustc.edu.cn