White matter microstructural alterations in amblyopic adults revealed by diffusion spectrum imaging with systematic tract-based automatic analysis.

BACKGROUND/AIM: We investigated the microstructural changes in white matter of adults with amblyopia using diffusion spectrum imaging with systematic tract-based automatic analysis of the whole brain.
METHODS: Ten adults with amblyopia (six women and four men, 33.6±10.6 years old on average) and 20 age- and sex-matched normal-sighted controls were enrolled. The mean generalised fractional anisotropy (GFA) was measured in 76 white matter tracts and compared between the experimental and control groups using a threshold-free cluster-weighted method and t-test. A 2-percentile cut-off was used to identify segments with the greatest differences between the two groups.
RESULTS: Participants with amblyopia had significantly lower GFA values than the controls in 11 segments located in nine white matter tracts, which included the following: left arcuate fasciculus, left frontal aslant tract, left fornix and left inferior fronto-occipital fasciculus of the association fibres; left thalamic radiations of the auditory nerve and bilateral optic radiations of the projection fibres; and genu and middle temporal gyrus of the callosal fibres. Amblyopic participants had statistically higher GFA values in the bilateral uncinate fasciculus than those of the controls.
CONCLUSION: This preliminary study using whole-brain tractographic analysis of white matter reveals association between abnormal early visual processing and alterations in brain architecture, which may be related to various higher-level deficits, such as audiovisual integration and hand-eye coordination in patients with amblyopia. © Author(s) (or their employer(s)) 2019. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Introduction

Amblyopia, with a prevalence of 2.9%–3.2% in adults,1 is a common cause of uncorrectable visual impairment. This neurodevelopmental disorder of the visual cortex occurs when binocular vision is disrupted during early childhood, leads to lifelong visual deficits affecting both monocular and binocular visual function.2 Common causes of amblyopia include refraction error, anisometropia, strabismus and visual deprivation arising from ptosis or congenital cataract.1 2

Amblyopia reduces corrected visual acuity and is associated with lower-level visual deficits, for example, reduced contrast sensitivity, and higher-level deficits, for example, those involving global processing tasks (such as form perception, contour processing and motion integration) or abnormal processing of information in the ventral and dorsal pathways.2 Abnormal early visual experiences in amblyopia subjects can affect the initial stages of visual–tactile and audiovisual integration and trigger changes in long-term crossmodal plasticity.3–5 Additionally, amblyopia also influences motor functions.6 7

Neuroimaging studies have revealed altered anatomical organisation and neural responses in the lateral geniculate nucleus (LGN) and in striate and extrastriate areas of the visual cortex.8 9 Diffusion-weighted magnetic resonance imaging (MRI) is currently one of the main non-invasive techniques used for in vivo mapping of white matter architecture of the human brain. White matter comprises bundles of axons that connect various cortical or subcortical areas. Diffusion tensor imaging (DTI) has been used to study microstructural changes in white matter of amblyopia subjects, leading to the identification of developmental impairments in the optic radiation (OR) and several fasciculi that are associated with dorsal and ventral visual pathways.10–13 Diffusion spectrum imaging (DSI) is another kind of diffusion-weighted MRI, which can resolve within-voxel heterogeneity of fibre orientations, analyse crossed fibres within white matter, and follow fibre fascicles into the grey matter of cortical and subcortical regions.14 15 The advantage of DSI is the ability to generate anatomically complete results and accurately delineate major neural fibres in complex local environments.14

When diffusion MRI is combined with tractography, the microstructural properties along a specific fibre tract bundle can be determined manually after delineating regions of interest. However, traditional tract-specific analysis is highly operator-dependent and time-consuming, and thus prone to having reconstruction biases and errors. The newly developed whole-brain tract-based automatic analysis (TBAA) has been used for automated whole-brain investigations of over 76 major white matter tracts, using a high-quality DSI template constructed from images acquired from 122 healthy young adults. This comprehensive and standardised method yields results with high geometrical agreement and considerably less measurement variability than manually reconstructed tracts.16 17

We hypothesised that abnormal early visual experience causes alteration of white matter microstructure in regions involved in visual processing and in areas important for linking visual inputs to the correlated motor tasks. Whole-brain investigations with DSI and TBAA were performed in this study to systematically compare white matter integrity in amblyopic adults with that in normal-sighted controls.

Methods

Participants

The experimental group consisted of subjects aged 20–50 years diagnosed with amblyopia in the ophthalmology clinic at the National Taiwan University Hospital. Amblyopia was defined with a best-corrected visual acuity of less than 20/20 in one or both eyes or a difference of best-corrected visual acuity more than two lines between two eyes. Anisometropic, strabismic, as well as deprivated amblyopia were included in the study. Subjects with neurological or psychiatric disease were excluded. A full-time paediatric ophthalmologist (T-HT), assessed and diagnosed all the subjects. For every amblyopic subject, two age- and sex-matched subjects without strabismus or amblyopia and with normal stereovision were recruited by poster in local communities around the National Taiwan University as controls. All participants underwent a best-corrected visual acuity test, slit lamp and fundus examination, stereotest with random-dot stereogram 300” (NTU 300”) and ocular motility test.

The mean age±SD was 33.6±10.6 years in the amblyopia group and 33.6±10.2 years in the control group (p=0.88). There were six women in the amblyopia group (6/10) and 12 in the control group (12/20). Table 1 summarises the demographic data of the amblyopia group. Subject 1 had history of monocular congenital cataract and underwent surgery at 4 years of age. Subjects 1, 3 and 10 had exotropia in the affected eye. None of the amblyopic participants passed the screening with random-dot stereogram 300”.

Table 1

Clinical characteristics of the participants with amblyopia

No.	Sex	Age, years	AE	Handedness	BCVA	BCVA	Refraction	Refraction	Education	Others
					R	L	R	L
1	M	48	R	L	0.1	1.5	−2.50/–2.00, 15°	−1.75/–0.75, 150°	middle school	RXT, cataract
2	M	35	L	R	1.5	0.1	−4.50/–0.75, 170°	+5.0/–1.25, 180°	graduate school
3	M	23	R	R	0.6	1.0	−1.50/–4.25, 10°	Plano	college	RXT
4	F	28	R	R	0.6	1.2	+4.75	−0.25	college
5	F	29	R	R	0.3	1.0	+2.25/–0.75, 162°	−4.25/–0.50, 180°	college
6	F	30	L	R	1.0	0.6	−3.75/–0.75, 175°	+0.25/–1.75, 180°	high school
7	F	49	L	L	1.0	0.1	−6.0/–0.75, 180°	−1.25/–0.75, 180°	high school
8	F	24	R	R	0.1	0.5	−3.5/–0.50, 160°	0/–2.00, 15°	college
9	F	47	L	L	1.2	0.15	−0.25	+3.5/–0.25, 130°	graduate school
10	M	23	R	R	0.05	1.0	−8.00/–1.00, 15°	−6.75/–1.25, 160°	college	RXT

AE, affected eye; BCVA, best-corrected visual acuity; RXT, right eye exotropia.

This study complied with the Health Insurance Portability and Accountability Act. Written informed consent was obtained prospectively from the study participants after detailed explanations, according to a protocol conforming to the tenets of the Declaration of Helsinki.

MRI acquisition

MR images were acquired using a 3 T MRI system (Trio, Siemens, Erlangen, Germany) with a 32-channel phased-array head coil. The slice orientation of the axial images was defined on the sagittal localiser as the orientation parallel to the line between the anterior and posterior commissures. High-resolution T1-weighted images were acquired using a three-dimensional magnetisation-prepared rapid gradient-echo sequence (parameters: repetition time=2000 ms; echo time=3 ms; flip angle=9°; matrix size=256×192×208; field of view=256×192×208 mm3), resulting in an isotropic spatial resolution of 1 mm3. DSI was performed using a single-shot spin-echo echo-planar imaging sequence with twice-refocused balanced echo for reducing eddy current distortions. Thirty-nine axial images that covered the entire cerebrum were acquired contiguously (parameters: repetition time=9600 ms; echo time=130 ms; matrix size=80×80; field of view=200×200 mm2; slice thickness=2.5 mm). According to a modified optimum sampling scheme, a total of 102 volumes of diffusion-weighted images was acquired with diffusion-encoding gradient directions pointing at the grid points within a half-sphere of the diffusion-encoding space (q-space); the maximum diffusion sensitivity (bmax) was 4000 s/mm2.

DSI data reconstruction and TBAA

Before reconstruction, the DSI data underwent quality assurance entailing signal drop-out screening, signal-to-noise ratio screening and shaking artefact correction. For each voxel, the 102 diffusion samples were processed using ReMAP estimation, which is a regularised version of Mean Apparent Propagator MRI.18 The generalised fractional anisotropy (GFA) was calculated from the orientation distribution function. GFA represents the diffusion-relevant anisotropy of local architectures, which infers the microstructural integrity of white matter, analogous to the use of fractional anisotropy (FA) in DTI. It ranges from 0 (completely isotropic diffusion) to 1 (diffusion is markedly restricted to one dimension).19 The technical details of TBAA have been described previously.16 17 The position coordinates of the 76 tracts were transformed from the NTU-DSI-122 template to the individual DSI data sets with the deformation maps, and the coordinates were aligned with the proceeding direction of each tract bundle. Each aligned tract was normalised and divided into 100 sampling steps along the sampling coordinates and recorded as the sampling steps in the template space. Linear regression analysis was used to remove the effects of age and sex in each step. The GFA values were sampled in the native space along the sampling steps. Finally, a 3D connectogram (x-axis: steps along sampling coordinates; y-axis: 76 white matter tracts; z-axis: number of each subject) was obtained as the output of TBAA.

Threshold-free Cluster-weighted (TFCW) Scores for group analysis

To compare the differences between the two groups by looking at each sampling step of the tracts and finding the steps with the most significant differences, a total of 76×100 comparisons had to be made. Since no significance of difference would be retained after correcting for this number of multiple comparisons, we adopted a statistical method proposed by Smith and Nichols using TFCW Scores.20 The TFCW Score rescaled the effect size (ES) at each step by considering the contributions of the lengths of continuous segments at different thresholds.

The TFCW Score (Sp) for step p was defined as:

Threshold levels (e) for ES were set from 0.1 to 1 in 10 aliquots, while the Cp(e) was the segment step size at the threshold level e. The value of S was simply ES at step p (ES) multiplied by an area extended from step p under the ES profile (figure 1). This area was equivalent to the length summation of the continuous step segments that contained step p. After calculating the S of all the steps in the connectogram, we plotted a histogram with TFCW Scores (7600 in total) and selected those in the top 2-percentile as the segments with noticeable between-group differences. Figure 1 shows the procedures in the analysis with TFCW Scores. Although the validity of the results seems less strong than that obtained from standard statistics with multiple comparison correction, it is better than those simply using uncorrected p values with a certain threshold of the cluster size. Mean GFA values of the segments that passed the 2-percentile cut-off were calculated and compared by t-test between two groups for further validation.

Figure 1

Threshold-free Cluster-weighted (TFCW) analysis. (1) Tract profiles of the callosal fibres of genu in average generalised fractional anisotropy (GFA) are shown as representative tracts (normal-sighted controls: blue; subjects with amblyopia: magenta). (2) Computing the effect size for each step. (3) Calculating the TFCW Score for each step. (4) Steps with scores in the top 2-percentiles of the histogram were selected. (5) The green segments indicate higher GFA values in controls than in subjects with amblyopia. The red segments indicate lower GFA values in controls than in subjects with amblyopia.

Results

Data from 10 amblyopia subjects and 20 matched controls were analysed. We identified 162 steps that were in the top 2-percentile of the histogram. These steps formed 14 segments in 12 tracts (figures 2–4). Among the association fibres, the tracts included the left arcuate fasciculus (AF), left frontal aslant tract (FAT), left fornix, left inferior fronto-occipital fasciculus (IFOF), and bilateral uncinate fasciculi (figures 2 and 3). Among the callosal fibres (CF), these tracts were the CF of the genu and the CF connecting the middle temporal gyri (figure 3). Among the projection fibres, the tracts were the left frontal striatum (FS) of the orbitofrontal cortex (OFC), left thalamic radiations (TRs) of the auditory nerve and the optic radiations in both hemispheres (figure 4).

Figure 3

The threshold-free cluster-weighted (TFCW) method identified two association fibre tracts and two callosal fibre (CF) tracts that contained segments that were different between the two study groups. Normal-sighted controls, blue; subjects with amblyopia, magenta; black bars under the tract profile indicate the effect size for each step; yellow rectangles indicate the segments that stand out as top 2-percentile differences. GFA, generalised fractional anisotropy.

Figure 4

The threshold-free cluster-weighted (TFCW) method identified four projection fibre tracts that contained segments that were different between the two study groups. Normal-sighted controls, blue; subjects with amblyopia, magenta; black bars under the tract profile indicate the effect size for each step; yellow rectangles indicate the segments that stand out as top 2-percentile differences. FS, frontal striatum; GFA, generalised fractional anisotropy; OFC, orbitofrontal cortex; TR, thalamic radiation.

The threshold-free cluster-weighted (TFCW) method identified four association fibre tracts that contained segments that were different between the two study groups. Normal-sighted controls, blue; subjects with amblyopia, magenta; black bars under the tract profile indicate the effect size for each step; yellow rectangles indicate the segments that stand out as top 2-percentile differences. GFA, generalised fractional anisotropy; IFOF, inferior fronto-occipital fasciculus.

The mean GFA values of these 12 white matter tracts are listed in table 2. The left AF and CF of the genu contained two segments each; therefore, the mean GFA values were determined by averaging the GFA values over the two segments on each tract. In 11 of the 12 tracts, the GFA values were significantly different between the amblyopia and the control groups (p<0.05). Subjects with amblyopia had higher mean GFA values in bilateral uncinate fasciculi and lower mean GFA values in the other nine tracts compared with controls. The GFA value in the segment of the left FS of the OFC was slightly lower in the amblyopia group than in the control, although not statistically significantly (p=0.072).

Table 2

The threshold-free cluster-weighted (TFCW) method revealed 12 white matter tracts with different generalised fractional anisotropy (GFA) values between amblyopic adults and normal-sighted controls

Tract No.	Name of tracts	GFA valuesMean±SD		P value
		Controls	Amblyopia
1	Left optic radiation	0.659±0.038	0.623±0.038	0.0176
2	Right optic radiation	0.654±0.055	0.594±0.044	0.0038
3	Left IFOF	0.717±0.026	0.682±0.047	0.0072
4	Left arcuate fasciculus	0.643+0.026	0.598±0.044	0.0022
5	Left TR of the auditory nerve	0.371±0.030	0.342±0.030	0.0110
6	CF of the middle temporal gyrus	0.161±0.041	0.128±0.018	0.0111
7	Left frontal aslant tract	0.609±0.054	0.554±0.043	0.0063
8	Left fornix	0.560±0.079	0.475±0.088	0.0102
9	CF of the genu	0.592±0.037	0.539±0.042	0.0007
10	Left uncinate fasciculus	0.482±0.050	0.520±0.045	0.0022
11	Right uncinate fasciculus	0.341±0.033	0.390±0.068	0.0056
12	Left FS of the OFC	0.346±0.048	0.311±0.046	0.0720

CF, callosal fibres; FS, frontal striatum; IFOF, inferior fronto-occipital fasciculus; OFC, orbitofrontal cortex; TR, thalamic radiation.

Discussion

Using DSI with TBAA to assess subjects with amblyopia, we found 12 white matter tracts containing segments with differences that were ranked in the top 2-percentile as compared with healthy controls. Results of the whole-brain tractography suggested association between abnormal early visual processing and brain architectural changes.

Early visual experience is essential in building the neural network of the visual system. Although white matter architecture can change with age and learning experience, neural plasticity is at its highest in young children during a period termed ‘the critical period of vision development’. Disruption of visual input during this period may result in amblyopia, and functional recovery is difficult when neuroplasticity slows down at the end of the critical period. Amblyopia had been found to be associated with changes of functional brain organisation in primary visual cortex and LGN, volume loss in visual cortical regions and alteration of white matter fasciculi associated with visual processing.8–13

On comprehensive evaluation of major white matter fibre architecture in amblyopic adults, we found that the amblyopia group had lower mean GFA values in fascicles involved in visual processing. The OR connects the visual cortex and the LGN, and the IFOF is associated with the ventral stream of visual pathway. A previous study that used DTI to investigate children with anisometropic amblyopia showed a bilateral decrease in the FA in the OR.10 12 In another DTI study, altered FA values were found in the right OR, left inferior longitudinal fasciculus and IFOF.11

Interestingly, recent studies showed that amblyopia is associated with altered ‘audiovisual integration,’ which manifests as a reduced McGurk effect. This defect of fusing incongruent auditory and visual syllables may persist into adulthood in patients with amblyopia.4 5 Some of the altered white matter tracts found in our study are known to be involved in language and auditory functions. The AF is one of the four subcomponents of the superior longitudinal fasciculus, and is regarded as a part of the dorsal stream for language which plays a role in mapping acoustic sounds to articulatory representations.21 The left TR of the auditory nerve connects the thalamus to Heschl's gyrus, which processes incoming auditory information. The CF connecting the middle temporal gyri may be involved in language processing.22

Some authors reported that children with amblyopia were impaired in fine motor skill and hand−eye coordination tasks, particularly manual dexterity tasks requiring speed and accuracy.6 7 The frontal aslant tract (FAT), which connects the inferior frontal gyrus with the supplementary and presupplementary motor areas, is regarded as part of a 'motor stream' involved in visuomotor processing. The bilateral FAT is associated with the closing grip phase of hand kinematics, which determines efficient and stable grasping of a target object.23 Altered FAT might be involved in the neuromechanism underlying defective upper limb coordination in patients with amblyopia.

In a resting-state functional MRI study, Lin et al found that adults with anisometropic amblyopia had reduced regional homogeneity in the prefrontal cortices and the OFC.24 Interestingly, our study found that the GFA value was lower for a segment of CFs of the genu, which connect the bilateral prefrontal cortices and OFCs. The uncinate fasciculi are bidirectional, long-range white matter tract that connects lateral OFCs to the anterior temporal lobes, and they were the only white matter tracts of the amblyopia group with higher mean GFA values than those of the controls in the present study. We speculated that it is the result of overcompensation during development as the uncinate fasciculus transmits information from the visual ventral stream to the OFC for visual discrimination.

To our knowledge, DSI and TBAA have not previously been used to investigate white matter organisation in adults with amblyopia. Compared with previous studies of human amblyopia using white matter tractography, the major advantage of our study is that we were able to screen 76 major fibre tracts comprehensively. However, without assessing behavioural parameters and structural MRI data in the same patients, any attempt to link the present data to behaviour must remain speculative.

The primary limitation of our study is its small sample size. A larger study with participants from different age groups and categorisation of amblyopia is warranted, and the correlation between the effect of eye dominance, visual acuity and the level of stereopsis on white matter structure should be analysed to gain a better understanding of neural mechanisms underlying amblyopia.

In conclusion, in this preliminary study, we found reduced GFA in 11 segments of nine white matter fibre tracts in adults with amblyopia. Our results suggest that abnormal early visual experience may induce structural alterations in white matter of the brain, which may be associated with various deficits, such as visual processing, audiovisual integration and hand–eye coordination. TBAA and DSI are potentially promising tools for evaluating microstructural integrity in individuals with amblyopia and for elucidating the possible neural mechanisms underlying the disease symptoms.