Skip to main content

Abnormal intra-network architecture in extra-striate cortices in amblyopia: a resting state fMRI study



Amblyopia (lazy eye) is one of the most common causes of monocular visual impairment. Intensive investigation has shown that amblyopes suffer from a range of deficits not only in the primary visual cortex but also the extra-striate visual cortex. However, amblyopic brain processing deficits in large-scale information networks especially in the visual network remain unclear.


Through resting state functional magnetic resonance imaging (rs-fMRI), we studied the functional connectivity and efficiency of the brain visual processing networks in 18 anisometropic amblyopic patients and 18 healthy controls (HCs).


We found a loss of functional correlation within the higher visual network (HVN) and the visuospatial network (VSN) in amblyopes. Additionally, compared with HCs, amblyopic patients exhibited disruptions in local efficiency in the V3v (third visual cortex, ventral part) and V4 (fourth visual cortex) of the HVN, as well as in the PFt, hIP3 (human intraparietal area 3), and BA7p (Brodmann area 7 posterior) of the VSN. No significant alterations were found in the primary visual network (PVN).


Our results indicate that amblyopia results in an intrinsic decrease of both network functional correlations and local efficiencies in the extra-striate visual networks.


Amblyopia (lazy eye), fundamentally a neurological disorder, is characterized by reduced vision in an otherwise normal eye with the presence of an amblyogenic factor, including early child strabismus (ocular misalignment), anisometropia (difference in refractive error), or ametropia (large symmetric refractive errors) and, more rarely, image deprivation (for review, please see Holmes and Clarke [1]). Extensive neuronal physiological studies have found wide-spread within-cortex neural dysfunctions in amblyopic animal models, including a loss of binocularity in V1 [2], an alteration in the excitatory-inhibitory balance of cortical binocular cells [3], disruption of neuronal receptive-field structures [4], and a degradation of neuronal signals [5]. There is also human imaging evidence that amblyopes have abnormal human middle temporal cortex (hMT) response to pattern motion [6], object-related abnormality in high-order occipitotemporal cortex [7], reduced neural adaptation effects in visual cortices [8] and reduced fidelity of spatial representation to the amblyopic eye’s stimulation [9]. Given the fact that the brain is an extraordinarily complex and highly organized network in which dysfunction can spread easily between linked cortices [10, 11], it is so far unclear how the brain neural network is altered by amblyopia.

The analysis of resting-state functional magnetic resonance imaging (rs-fMRI) provides an effective way to assess the brain’s spontaneous activity and connectivity. Brain regions showing synchronized fluctuations during rs-fMRI form the intrinsic connectivity networks (ICNs), which have been demonstrated to provide the physiological basis for cortical information processing, and to be able to abstract and suspend performance influences combined with various neurological diseases [12]. In recent years, the graph theoretical analysis, which defined a graph as a set of nodes (brain areas) and edges (structural or functional connectivity), provided a powerful tool to examine the topological organization of complex brain networks [13,14,15,16]. Through these approaches, human brain networks have become an optimum small-world and economical topology [17], represented in characteristics of high global and local efficiency of parallel information processing at a low connection cost [18]. Since then, graph theory analysis has been widely and successfully used to explore the brain network architecture in development and neurological diseases, e.g., maturation [19], aging [20], schizophrenia [21], obsessive-compulsive disorder [22], and so on. However, only few rs-fMRI studies investigated amblyopic intrinsic functional connectivity e.g., Ding, Liu [23] found altered connectivity between the primary visual cortex (V1) with the cerebellum and the inferior parietal lobule; Wang, Li [24] have figured out decreased functional connectivity density in the visual ICNs of amblyopic children; and Mendola, Lam [25] have revealed abnormal retinotopically organized functional connectivity of visual areas in amblyopia. It remains unknown whether and how the local efficiency of the brain network evolves from the amblyopes’ abnormal visual experiences.

Here, we measured rs-fMRI to assess 3 ICNs in the visual information processing in 18 healthy volunteers and 18 anisometropic amblyopes: the higher vision network (HVN), the primary visual network (PVN), and the visuospatial network (VSN). Both the intra- and inter-network functional connectivity, as well as the network local efficiency of the visual ICNs were studied. Our results suggest widespread disturbances of functional connectivity and local efficiency in the extra-striate visual networks in amblyopia.

Materials and methods


A group of adult anisometropic amblyopes (n = 18, mean age: 23.7 ± 1.9 years old) and a group of healthy controls (n = 18, mean age: 25.2 ± 1.8 years old) participated. Anisometropia was defined as refraction differing by 1.0 diopters (D) or more for the two eyes; amblyopia was defined as reduced visual acuity (> 0.1 LogMAR) an otherwise normal eye due to abnormal visual experience early in life. A brief summary of participants’ clinical data is provided in Table 1. A comprehensive eye examination was carried out by a clinician (the author LF) at the First Affiliated Hospital of Anhui Medical University. No strabismic or eye-pathological cases were reported in the participating individuals. Healthy controls had normal or corrected to normal visual acuity (≤ 0.1 logMAR) in both eyes. All subjects were naive as to the purpose of the experiment. A written informed consent was obtained from each participant after explanation of the nature and possible consequences of the study. This study complied with the Declaration of Helsinki and was approved by the University of Science and Technology of China’s Research Ethics Committee.

Table 1 Clinical characteristics of the amblyopes and healthy controls

Image acquisition

The collection of magnetic resonance images was performed using a 3.0 T (Signa HDx; GE Healthcare, Illinois, United States) scanner with an eight-channel head coil. Foam padding and earplugs were used to minimize participants’ head motion and reduce scanner noise. Each scanning session began with an acquisition of high-resolution three-dimensional T1 weighted images using an MP-RAGE sequence (TR = 2300 ms; TE = 2.94 ms; flip angle = 9°; 176 slices; voxel size = 1 × 1 × 1 mm3). The resting state functional scans were T2*-weighted, gradient-echo, planar images (TR = 2000 ms; TE = 28 ms; flip angle = 72°; 40 slices; voxel size = 3 × 3 × 3 mm3). While acquiring resting state images, patients were instructed to keep both eyes closed and to think of nothing in particular.

MRI data analysis

For fMRI data, preprocessing was applied by using Analysis of Functional NeuroImages (AFNI) software tools (Medical College of Wisconsin, Milwaukee, WI, USA; Cox [26]). The first five images of each resting state session were discarded. Each subject’s fMRI data were registered first to his/her anatomical raw data by linear and quadratic registration. This was followed by slice-timing correction, head motion correction, spatial Gaussian smoothing with a kernel width of 6 mm at half maximum (FWHM), as well as temporal detrending respect to the head motion and order 3 polynomial drift correction. Then, all the images were normalized to the Montreal Neurological Institute (MNI) ICBM152 and were resliced by 3.0 × 3.0 × 3.0 mm3 voxels. The motion was assessed and the time points with framewise displacement over 0.2 mm were censored. Data was further filtered with temporal band pass 0.01~0.10 Hz. A regression of motion parameters and their derivatives were applied, and the residual error time series were obtained for further analysis.

Our network nodes were constructed by 19 Regions of interest (ROIs) (Fig. 1, Richiardi, Altmann [27]), part of the Willard 499 ROIs, constituting the primary visual network (PVN), higher visual network (HVN), and visuospatial network (VSN); these ROIs were asymmetrically distributed across the brain. ROI labels were in line with brain anatomy by matching the ROI center coordinates to the AFNI Anatomy Toolbox. Functional connectivity between ROIs were calculated using multivariate distance correlation [28]: For example, suppose areas A and B had t time-points, and vA and vB voxels, respectively. First, z-transfer was applied to each voxel’s time course by its mean and variance; then, the Euclidean distance, dA:t1,t2 and dB:t1,t2, between each pair of time points t1 and t2 was computed for each region:

$$ {d}_{A:t1,t2}=\sqrt{\sum_{v=1}^{v_A}{\left({A}_{v,t1}-{\mathrm{A}}_{v,t2}\right)}^2}\kern0.90em \forall t1,t2=1,\dots, t $$
$$ {d}_{B:t1,t2}=\sqrt{\sum_{v=1}^B{\left({B}_{v,t1}-{B}_{v,t2}\right)}^2}\kern0.90em \forall t1,t2=1,\dots, t $$
Fig. 1

An illustration of the 19 ICN ROIs from the three ICNs (HVN, PVN, and VSN) employed in our study

U-centering was applied to set row and column means to zero.

$$ {D}_{A:t1,t2}=\left\{\begin{array}{c}{d}_{A:t1,t2}-\frac{1}{t-2}{\sum}_{p=1}^t{d}_{A:t1,p}-\frac{1}{t-2}{\sum}_{q=1}^t{d}_{A:q,t2}+\frac{1}{\left(t-1\right)\left(t-2\right)}{\sum}_{p,q=1}^t{d}_{A:q,p},\mathrm{t}1\ne t2\\ {}0,t1=t2\end{array}\right. $$

The distance correlation, dCor, was then computed as follows:

$$ dCor\left(A,B\right)=\left\{\begin{array}{c}\sqrt{dCov\left(A,B\right)/\sqrt{dVar(A) dVar(B)}}, dCov\left(A,B\right)>0\\ {}0, dCov\left(A,B\right)\le 0\end{array}\right. $$

where dCov was distance covariance and dVar was distance variance.

$$ dCov\left(A,B\right)=1/t\left(t-3\right){\sum}_{t1,t2=1}^t{D}_{A:t1,t2}{D}_{B:t1,t2} $$
$$ dVar(A)=1/t\left(t-3\right){\sum}_{t1,t2=1}^t{D}_{A:t1,t2}^2 $$

The method is similar to the well-established univariate functional connectivity analysis [28], but allows inference based on multivoxel information within each ROI rather than the averaged global BOLD time series. The distance correlation, a metric of multivariate dependence of high dimensional vectors [29], is more reliable and robust than univariate methods [30, 31].

Network edges were obtained by Fisher-transformed distance correlation (z = 0.5ln[(1 + dCor)/(1- dCor)], where dCor is the distance correlation between the time series of each ROI, resulting in a 19*19 functional correlation matrix for each subject. These matrices were further used for the network analysis and for graph theoretical analysis using the GRETNA toolbox [32]. The nodal local efficiency is defined as the harmonic mean of the inverse of the l, which is the minimum value of the sum of weights over all possible paths between the immediate neighborhood nodes of the node:

$$ {E}_{local}=\frac{1}{N_{G_i}\left({N}_{G_i}-1\right)}{\sum}_{j,k\in {G}_i}\frac{1}{l_{j,k}} $$

where the subgraph Gi is defined as the set of nodes that are directly connected by a single edge to the ith node, and NGi is the number of nodes in the Gi [18].

The functional connectivity analysis was conducted by the programs in the MATLAB (MathWorks, Natick, MA). The correlation between network node and other node within one network (HVN, PVN, or VSN) is defined as intra-network connectivity, whereas the correlation between a node of a given network and that of another network is defined as the inter-network connectivity. Fisher transformation was applied to intra-and inter-network matrices of each subject to yield mean Fisher-transformed correlation values. The group differences in average intra-network or inter-network were assessed by the repeated-measured analysis of variance (ANOVA), and the false-discovery-rate (FDR) corrected t-test, with p-value corrected according to the Algorithm 2 by Storey [33]. Linear correlation analysis were also applied to assess the relationship between the visual acuity and the network connectivity of the amblyopes.

To characterize network efficiency, the local efficiency (LE) of each visual ICN node was computed as a function of the minimum path length between regions [18]. A series of sparsity threshold (0.2 ≤ sparsity ≤ 0.8, interval = 0.05) were applied to measure the individual correlation matrices, for there was no gold standard for selecting a proper single sparsity threshold. The LE at each sparsity was calculated and the area under the curve (AUC) for LE was employed to be a summarized scalar [14]. Group differences in AUC of LE (aLE) of each network node were reported after FDR-corrected t-test ([FDR-corrected] P < 0.05) separately.


Functional connectivity analysis

Figure 2 shows the results of the functional connectivity analysis for amblyopes (Fig. 2a) and healthy controls (HC; Fig. 2b). As expected, both the HC and the amblyopic matrices showed more positive correlations within each network than those between networks. The amblyopic matrix (Fig. 2a) showed generally reduced correlations compared with the HC matrix (Fig. 2b). This difference was apparent in correlation difference matrix (Amblyopia minus HC) shown by Fig. 2c.

Fig. 2

ICN nodes connectivity and group differences. 19–19 matrices were computed in all ROIs for all ICNs in amblyopes (a) and healthy controls (b). The nodes are grouped by ICNs. The intra-network connectivity is plotted as diagonal colored blocks and the inter-network connectivity is plotted as off-diagonal blocks. The group differences are plotted with 19–19 matrices (c) and bar graphs (d/e). Error bars represent standard deviations; *: P < 0.05, FDR corrected

For the intra-network connectivity, we used a mixed-design repeated-measured ANOVA, with group (amblyopes vs. HCs) as the between-subject factor and with intra-network connectivity (HVN, PVN, and VSN) as the within-subject factor. ICN nodes connectivity was significantly different between these two groups (F (1,34) = 4.21, P = 0.048); such group difference existed in all the three intra-network conditions, as the interaction between the group and the intra-network was not significant (F (2,68) = 0.15, P = 0.86). Two-sample t-test of each ICN further showed that the connectivity within the VSN (t (34) = − 2.64, Uncorrected [FDR corrected] P = 0.0124 [0.037]) and the HVN (t (34) = − 2.29, Uncorrected [FDR corrected] P = 0.0284 [0.043]) were significantly decreased in amblyopes (Fig. 2d).

Secondly, we performed a mixed-design repeated-measured ANOVA, with group (amblyopes vs. HCs) as the between-subject factor and with inter-network connectivity (HVN-PVN, HVN-VSN, and PVN-VSN) as the within-subject factor. Connectivity was significantly different across the inter-network pairs (F (2,68) = 31.2, P < 0.001), while neither the between-group effect (F (1,34) = 1.17, P = 0.29) nor the interaction effect (F (2,68) = 1.21, P = 0.31) was significant. No significant alteration of between network connectivity was observed in amblyopes (Fig. 2e) after the FDR corrected t-test. We then applied Pearson’s correlation analysis and found that neither intra-network connectivity nor inter-network connectivity showed any significant correlation with the corrected visual acuity in the amblyopes (P > 0.5).

The local efficiency analysis

To further investigate the effects of amblyopia within the visual ICNs, we conducted a local efficiency analysis. In Fig. 3, we plotted the averaged visual ICNs of amblyopes (Fig. 3a) and the HCs (Fig. 3b). The visual ICNs of both groups demonstrated small-world network architecture and the local efficiencies at the extra-striate cortices were significantly decreased in amblyopes, evidenced by the smaller node size illustrated in Fig. 3. A mixed ANOVA, with group (amblyopes vs. HCs) as the between-subject factor and with network nodes as the within-subject factor also showed that there was a significant difference between groups (F (1,34) = 6.27, P = 0.017) and nodes (F (18,612) = 2.6, P < 0.001). We further conducted a series of t-tests (amblyopes vs. healthy controls) across all the 19 visual ICN nodes. The t-test analysis showed that the aLE at the lPFt, lhIP3, lBA7p, rhIP3, lV3v, rV3v, and rV4 were significantly smaller in amblyopes than that in healthy controls (Table 2).

Fig. 3

A widespread of extra-striate cortices showed significant decrease of aLE in the amblyopia. The LE was calculated at a series of sparsity threshold (0.2 ≤ sparsity ≤ 0.8, interval = 0.05), and the area under the curve (AUC) for LE (aLE) was obtained as a summarized scalar. FDR-corrected t-test showed a significant decrease of aLE in lPFt, lhIP3, lBA7p, lV3v, rhIP3, rV3v, and rV4 (labelled with arrows) of the amblyopia (smaller node sizes) (a) and healthy controls (b). Note that a sparsity of 0.35 was used here for illustration, with the sizes of the nodes proportional to the aLE of each node and ICNs labeled with different colors

Table 2 MNI coordinates, cortical regions, and the effect of amblyopia on the aLE of each node


Our principal finding is that anisometropic amblyopes suffer from a decrease of intra-network functional connectivity and local efficiency within the brain extra-striate cortices. To our limited knowledge, this study is the first demonstration of an intrinsic alteration of the brain extra-striate visual networks in adult amblyopes, which suggests an underlying pathological process engaged in amblyopia.

Our analysis included distance correlation to assess functional connectivity [29]. This multivariate method was similar to the standard univariate functional connectivity method in obtaining correlations between brain ROIs from resting state fMRI data [30]. Furthermore, the distance correlation, by using multivariate patterns to measure the dependences between two brain regions, could effectively detect the non-linearity while avoiding any within ROI signal averaging. The method is capable of encoding information of associations between brain areas that was lost by averaging [31, 34]. We have also employed a pre-computed group-level brain network parcellation, which represents the functional organization of the brain, and is integrally correlated with genes linked to synaptic function [27].

An important step in understanding how the abnormal visual experience of amblyopia influenced visual neural network is the generation of a map of the connectivity architecture of the brain. The application of network science and graph theory has enabled detailed descriptions of how disease affects the brain [12,13,14,15,16]. Through diffusion tensor imaging (DTI), amblyopic brain structural connectivity studies have found increased mean diffusivity (MD) in thalamo-cortical visual pathways [35] and vertical occipital fasciculus [36], as well as decreased fractional anisotropy (FA) in the optic radiation, inferior longitudinal fasciculus/inferior fronto-occipital fasciculus and superior longitudinal fasciculus [37]. Previous rs-fMRI works have reported altered connectivity between the primary visual cortex (V1) with the cerebellum and the inferior parietal lobule [23], decreased functional connectivity density in the visual ICNs in amblyopic children [24], as well as disrupted retinotopically functional connectivity of visual areas in amblyopes [25]. Through network analysis, our present work further delineated the amblyopic deficits in visual network architectures. The observed reduction in the present work suggests that amblyopes have a less efficient visual network compared with that of healthy controls.

A previous study has observed impaired visual functional connectivity in amblyopia while processing the visual information from the amblyopic eye [38]. By using rs-fMRI, we were able to extend the observation to the intrinsic functional connectivity, i.e., no visual inputs. We demonstrated reduced intra-network correlations within the HVN. The deleterious effects of amblyopia on HVN could also be localized in terms of reduced local efficiency of V3v and V4. Since the local efficiency shows how efficient the processing is between the immediate neighbors of a node when the node is removed, it reveals the degree of fault tolerance of the system [39]. Thus, the results suggest that the V3v and V4 were intrinsically less fault tolerant in amblyopes and can be interpreted to have a more fragile visual system intolerant to fault or conflicting information inputs [40, 41].

Furthermore, our results suggest a loss of functional connectivity within the VSN of amblyopes, as well as a reduction of local efficiency of the VSN nodes (hIP3, PFt and BA7p). The VSN are cortices that deal with processing of spatial working memory, visually guided action, eye movements and navigation [42]. The hIP3 has been found to be highly structurally and functionally connected to the visual cortex and plays an important role in attentional selection between peripherally presented stimuli [43]. The PFt participates in the action observation and imitation network [44], and the BA7p is a key hub of the VSN bridging to the executive network [45]. Wang, Crewther [46] have found that when amblyopes viewed visual motion stimulus through amblyopic eyes, both the activation and the functional connectivity of VSN were weaker compared to that while viewing through their fellow eyes. Through intrinsic functional network analysis, our results suggested that the amblyopic deficits reflected impaired neural synchronizations within the visuospatial network nodes. This is consistent with a recent study of pathological perturbations to widespread white matter fiber tracts in amblyopia [47]. Our results of reorganization of the visuospatial network that is remote from the primary visual cortex suggest functional pathological cascades encompassing large swathes of the visuospatial system in amblyopia. However, the question of how architecture alterations of the visual networks are linked to amblyopic clinical deficits requires further investigation.


In summary, we compared the visual ICNs of amblyopes with those of normal observers and found decreased intra-network functional connectivity and local efficiency in some brain areas within the visual ICNs. These findings suggest that amblyopes suffer from a reduction of both internal neural functional connectivity and local efficiency within extra-striates and visuospatial networks.

Availability of data and materials

The code supporting the findings of this study are available from the corresponding author upon request.



Analysis of variance


Area under curve


Brodmann area 44


Brodmann area 45


Brodmann area 7 anterior


Brodmann area 7 posterior


Blood oxygenation level dependent


Fractional anisotropy


False discovery rate


Front eye field


Functional magnetic resonance imaging


Healthy controls


Human intraparietal area 3


Human middle temporal cortex


Higher visual network


Intrinsic connectivity network


Inferior temporal cortex


Mean diffusivity


Primary visual network


Resting state functional magnetic resonance imaging


Echo time


Repetition time


Primary visual cortex


Secondary visual cortex


Third visual cortex, ventral part


Fourth visual cortex


Vertical occipital fasciculus


Visuospatial network


  1. 1.

    Holmes JM, Clarke MP. Amblyopia. Lancet. 2006;367(9519):1343–51.

    Article  Google Scholar 

  2. 2.

    Hubel DH, Wiesel TN. Binocular interaction in striate cortex of kittens reared with artificial squint. J Neurophysiol. 1965;28(6):1041–59.

    CAS  Article  Google Scholar 

  3. 3.

    Sengpiel F, Blakemore C. The neural basis of suppression and amblyopia in strabismus. Eye. 1996;10(Pt 2):250–8.

    Article  Google Scholar 

  4. 4.

    Tao X, Zhang B, Shen G, Wensveen J, Smith EL 3rd, Nishimoto S, et al. Early monocular defocus disrupts the normal development of receptive-field structure in V2 neurons of macaque monkeys. J Neurosci. 2014;34(41):13840–54.

    Article  Google Scholar 

  5. 5.

    Shooner C, Hallum LE, Kumbhani RD, Ziemba CM, Garcia-Marin V, Kelly JG, et al. Population representation of visual information in areas V1 and V2 of amblyopic macaques. Vis Res. 2015;114:56–67.

    Article  Google Scholar 

  6. 6.

    Thompson B, Villeneuve MY, Casanova C, Hess RF. Abnormal cortical processing of pattern motion in amblyopia: evidence from fMRI. NeuroImage. 2012;60(2):1307–15.

    CAS  Article  Google Scholar 

  7. 7.

    Lerner Y, Pianka P, Azmon B, Leiba H, Stolovitch C, Loewenstein A, et al. Area-specific amblyopic effects in human occipitotemporal object representations. Neuron. 2003;40(5):1023–9.

    CAS  Article  Google Scholar 

  8. 8.

    Li X, Coyle D, Maguire L, McGinnity TM, Hess RF. Long timescale fMRI neuronal adaptation effects in human amblyopic cortex. PLoS One. 2011;6(10):e26562.

    CAS  Article  Google Scholar 

  9. 9.

    Farivar R, Zhou J, Huang Y, Feng L, Zhou Y, Hess RF. Two cortical deficits underlie amblyopia: a multifocal fMRI analysis. NeuroImage. 2019;190:232–41.

    Article  Google Scholar 

  10. 10.

    Tononi G, Sporns O, Edelman GM. A measure for brain complexity: relating functional segregation and integration in the nervous system. Proc Natl Acad Sci U S A. 1994;91(11):5033–7.

    CAS  Article  Google Scholar 

  11. 11.

    Fornito A, Zalesky A, Breakspear M. The connectomics of brain disorders. Nat Rev Neurosci. 2015;16(3):159–72.

    CAS  Article  Google Scholar 

  12. 12.

    Ewers M, Sperling RA, Klunk WE, Weiner MW, Hampel H. Neuroimaging markers for the prediction and early diagnosis of Alzheimer's disease dementia. Trends Neurosci. 2011;34(8):430–42.

    CAS  Article  Google Scholar 

  13. 13.

    Bullmore E, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. 2009;10(3):186–98.

    CAS  Article  Google Scholar 

  14. 14.

    He Y, Evans A. Graph theoretical modeling of brain connectivity. Curr Opin Neurol. 2010;23(4):341–50.

    PubMed  Google Scholar 

  15. 15.

    Bassett DS, Bullmore E, Verchinski BA, Mattay VS, Weinberger DR, Meyer-Lindenberg A. Hierarchical organization of human cortical networks in health and schizophrenia. J Neurosci. 2008;28(37):9239–48.

    CAS  Article  Google Scholar 

  16. 16.

    Salvador R, Suckling J, Coleman MR, Pickard JD, Menon D, Bullmore E. Neurophysiological architecture of functional magnetic resonance images of human brain. Cereb Cortex. 2005;15(9):1332–42.

    Article  Google Scholar 

  17. 17.

    Watts DJ, Strogatz SH. Collective dynamics of 'small-world' networks. Nature. 1998;393(6684):440–2.

    CAS  Article  Google Scholar 

  18. 18.

    Achard S, Bullmore E. Efficiency and cost of economical brain functional networks. PLoS Comput Biol. 2007;3(2):e17.

    Article  Google Scholar 

  19. 19.

    Zhong S, He Y, Shu H, Gong G. Developmental changes in topological asymmetry between hemispheric brain white matter networks from adolescence to young adulthood. Cereb Cortex. 2017;27(4):2560–70.

    PubMed  Google Scholar 

  20. 20.

    Cotier FA, Zhang R, Lee TMC. A longitudinal study of the effect of short-term meditation training on functional network organization of the aging brain. Sci Rep. 2017;7(1):598.

    Article  Google Scholar 

  21. 21.

    Zhu J, Wang C, Liu F, Qin W, Li J, Zhuo C. Alterations of functional and structural networks in schizophrenia patients with auditory verbal hallucinations. Front Hum Neurosci. 2016;10:114.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Jung WH, Yucel M, Yun JY, Yoon YB, Cho KI, Parkes L, et al. Altered functional network architecture in orbitofronto-striato-thalamic circuit of unmedicated patients with obsessive-compulsive disorder. Hum Brain Mapp. 2017;38(1):109–19.

    Article  Google Scholar 

  23. 23.

    Ding K, Liu Y, Yan X, Lin X, Jiang T. Altered functional connectivity of the primary visual cortex in subjects with amblyopia. Neural Plast. 2013;2013:612086.

    Article  Google Scholar 

  24. 24.

    Wang T, Li Q, Guo M, Peng Y, Li Q, Qin W, et al. Abnormal functional connectivity density in children with anisometropic amblyopia at resting-state. Brain Res. 2014;1563:41–51.

    CAS  Article  Google Scholar 

  25. 25.

    Mendola JD, Lam J, Rosenstein M, Lewis LB, Shmuel A. Partial correlation analysis reveals abnormal retinotopically organized functional connectivity of visual areas in amblyopia. Neuroimage Clin. 2018;18:192–201.

    CAS  Article  Google Scholar 

  26. 26.

    Cox RW. AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res. 1996;29(3):162–73.

    CAS  Article  Google Scholar 

  27. 27.

    Richiardi J, Altmann A, Milazzo AC, Chang C, Chakravarty MM, Banaschewski T, et al. Correlated gene expression supports synchronous activity in brain networks. Science. 2015;348(6240):1241–4.

    CAS  Article  Google Scholar 

  28. 28.

    Geerligs L, Cam-Can, Henson RN. Functional connectivity and structural covariance between regions of interest can be measured more accurately using multivariate distance correlation. Neuroimage. 2016;135:16–31.

    Article  Google Scholar 

  29. 29.

    Székely GJ, Rizzo ML, Bakirov NK. Measuring and testing dependence by correlation of distances. Ann Stat. 2007;35(6):2769–94.

    Article  Google Scholar 

  30. 30.

    Geerligs L, Tsvetanov KA, Cam-Can, Henson RN. Challenges in measuring individual differences in functional connectivity using fMRI: the case of healthy aging. Hum Brain Mapp. 2017;38(8):4125–56.

    Article  Google Scholar 

  31. 31.

    Yoo K, Rosenberg MD, Noble S, Scheinost D, Constable RT, Chun MM. Multivariate approaches improve the reliability and validity of functional connectivity and prediction of individual behaviors. Neuroimage. 2019;197:212–23.

    Article  Google Scholar 

  32. 32.

    Wang J, Wang X, Xia M, Liao X, Evans A, He Y. GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front Hum Neurosci. 2015;9:386.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Storey JD. A direct approach to false discovery rates. J R Stat Soc Ser B (Stat Methodol). 2002;64(3):479–98.

    Article  Google Scholar 

  34. 34.

    Anzellotti S, Coutanche MN. Beyond functional connectivity: investigating networks of multivariate representations. Trends Cogn Sci. 2018;22(3):258–69.

    Article  Google Scholar 

  35. 35.

    Allen B, Spiegel DP, Thompson B, Pestilli F, Rokers B. Altered white matter in early visual pathways of humans with amblyopia. Vis Res. 2015;114:48–55.

    Article  Google Scholar 

  36. 36.

    Duan Y, Norcia AM, Yeatman JD, Mezer A. The structural properties of major white matter tracts in strabismic amblyopia. Invest Ophthalmol Vis Sci. 2015;56(9):5152–60.

    Article  Google Scholar 

  37. 37.

    Li Q, Jiang Q, Guo M, Li Q, Cai C, Yin X. Grey and white matter changes in children with monocular amblyopia: voxel-based morphometry and diffusion tensor imaging study. Br J Ophthalmol. 2013;97(4):524–9.

    Article  Google Scholar 

  38. 38.

    Li X, Mullen KT, Thompson B, Hess RF. Effective connectivity anomalies in human amblyopia. Neuroimage. 2011;54(1):505–16.

    Article  Google Scholar 

  39. 39.

    Latora V, Marchiori M. Efficient behavior of small-world networks. Phys Rev Lett. 2001;87(19):198701.

    CAS  Article  Google Scholar 

  40. 40.

    Zhou JW, Liu R, Feng LX, Zhou YF, Hess RF. Deficient binocular combination of second-order stimuli in amblyopia. Invest Ophthalmol Vis Sci. 2016;57(4):1635–42.

    Article  Google Scholar 

  41. 41.

    Huang CB, Zhou J, Lu ZL, Feng L, Zhou Y. Binocular combination in anisometropic amblyopia. J Vision. 2009;9(3):17.1–16.

    Article  Google Scholar 

  42. 42.

    Kravitz DJ, Saleem KS, Baker CI, Mishkin M. A new neural framework for visuospatial processing. Nat Rev Neurosci. 2011;12(4):217–30.

    CAS  Article  Google Scholar 

  43. 43.

    Gillebert CR, Mantini D, Peeters R, Dupont P, Vandenberghe R. Cytoarchitectonic mapping of attentional selection and reorienting in parietal cortex. NeuroImage. 2013;67:257–72.

    Article  Google Scholar 

  44. 44.

    Caspers S, Zilles K, Laird AR, Eickhoff SB. ALE meta-analysis of action observation and imitation in the human brain. Neuroimage. 2010;50(3):1148–67.

    Article  Google Scholar 

  45. 45.

    Humphreys GF, Lambon Ralph MA. Fusion and fission of cognitive functions in the human parietal cortex. Cereb Cortex. 2015;25(10):3547–60.

    Article  Google Scholar 

  46. 46.

    Wang H, Crewther SG, Liang M, Laycock R, Yu T, Alexander B, et al. Impaired activation of visual attention network for motion salience is accompanied by reduced functional connectivity between frontal eye fields and visual cortex in strabismic amblyopia. Front Hum Neurosci. 2017;11:195.

    CAS  Article  Google Scholar 

  47. 47.

    Tsai TH, Su HT, Hsu YC, Shih YC, Chen CC, Hu FR, et al. White matter microstructural alterations in amblyopic adults revealed by diffusion spectrum imaging with systematic tract-based automatic analysis. Br J Ophthalmol. 2019;103(4):511–6.

    Article  Google Scholar 

Download references


Not applicable


The author(s) disclose receipt of the following financial support for the research, authorship, and/or publication of this article: this study was supported by the National Natural Science Foundation of China (grant numbers 81501942, 81701665, 81500754) and by the Fundamental Research Funds for the Central Universities (grant number WK2100230016).

Author information




Conception and design of the experiments: JZ, XW. Acquisitions of MRI and subject data: QL, YH, HW, LF. Analysis of the data: ZL, QL, YW, XW. Original drafting of the paper: ZL, YH, XW. Critical revisions for the paper and the approval of the final version: BN, YW, GL, YZ, YH, BQ, JZ. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jiawei Zhou or Xiaoxiao Wang.

Ethics declarations

Ethics approval and consent to participate

This study complied with the Declaration of Helsinki and was approved by the University of Science and Technology of China’s Research Ethics Committee (ID:USTCEC201200006). A written informed consent was obtained from each participant after explanation of the nature and possible consequences of the study.

Consent for publication

All study subjects gave informant consent.

Competing interests

The author(s) declare no potential competing interest with respect to the research, authorship, and/or publication of this article.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lu, Z., Huang, Y., Lu, Q. et al. Abnormal intra-network architecture in extra-striate cortices in amblyopia: a resting state fMRI study. Eye and Vis 6, 20 (2019).

Download citation


  • Amblyopia
  • fMRI
  • Resting state
  • Extra-striate cortices
  • Visuospatial network
  • Graph analysis