Abstract
Visual working memory (VWM) plays an essential role in many perceptual and higher-order cognitive processes. Despite its reliance on a broad network of brain regions, VWM has a capacity limited to a few objects. This capacity varies substantially across individuals and relates closely to measures of overall cognitive function (Luck and Vogel, 2013). The mechanisms underlying these properties are not completely understood, although the amplitude of neural signal oscillations (Vogel and Machizawa, 2004) and brain activation in specific cortical regions (Todd and Marois, 2004) have been implicated. Variability in VWM performance may also reflect variability in white matter structural properties. However, data based primarily on diffusion tensor imaging approaches remain inconclusive. Here, we investigate the relationship between white matter and VWM capacity in human subjects using an advanced diffusion imaging technique, diffusion kurtosis imaging. Diffusion kurtosis imaging provides several novel quantitative white mater metrics, among them the axonal water fraction (faxon), an index of axonal density and caliber. Our results show that 59% of individual variability in VWM capacity may be explained by variations in faxon within a widely distributed network of white matter tracts. Increased faxon associates with increased VWM capacity. An additional 12% in VWM capacity variance may be explained by diffusion properties of the extra-axonal space. These data demonstrate, for the first time, the key role of white matter in limiting VWM capacity in the healthy adult brain and suggest that white matter may represent an important therapeutic target in disorders of impaired VWM and cognition.
- axonal density
- diffusion kurtosis imaging
- individual differences
- microstructure
- visual working memory capacity
- white matter
Introduction
Visual working memory (VWM; i.e., the short-term storage of visual information) is supported by a distributed neural system, which in humans includes occipital, prefrontal, parietal, and temporal areas (Ungerleider et al., 1998). Functional imaging has been extensively used to identify the cortical areas involved in VWM (Ungerleider et al., 1998; Zimmer et al., 2010). The relationship between VWM and the white matter pathways interconnecting the VWM network nodes has also been studied primarily using diffusion tensor imaging (DTI). DTI measures, such as fractional anisotropy (FA), quantify water molecules' microscopic diffusion and reflect a variety of microstructural features, among them the degree of myelination and fiber density and coherence. Most DTI studies examining the relationship between white matter and working memory (visual or verbal) have focused on pathological groups, such as patients with schizophrenia (Hanlon et al., 2012) or multiple sclerosis (Yu et al., 2012), pathological processes such as aging (Davis et al., 2009), or on early brain development (Klingberg, 2006; Vestergaard et al., 2011; Short et al., 2013; Peters et al., 2014). In healthy adults, modest relationships have been reported between FA and verbal working memory (Karlsgodt et al., 2008; Takeuchi et al., 2011).
A potential limitation of these studies was the use of DTI metrics, which, although highly sensitive to pathological processes, may lack specificity (Walker et al., 2012). Thus, we have used a recently proposed two-compartment model of white matter (Fieremans et al., 2011), based on diffusion kurtosis imaging (DKI), which separates the magnetic resonance signal contribution from intra- and extra-axonal white matter compartments. Several metrics can be obtained using this approach: (1) the axonal water fraction (faxon), which describes the volume of the intra-axonal water relative to the total intra- and extra-axonal volume and relates to axonal density and caliber (Barazany et al., 2009; De Santis et al., 2012); (2) the intra-axonal diffusivity (Daxon), which describes water diffusion within axons, primarily along the axonal axis, and thus reflects microstructural features, such as the size, number, and organization of intra-axonal structures; and (3) the extra-axonal axial (ADextra) and radial (RDextra) diffusivities, which reflect the overall microstructural organization of the extra-axonal space. Recent work has shown that these metrics are sensitive to white matter characteristics that are not described by typical DTI metrics, such as differences in axonal density and caliber in autism (Lazar et al., 2014) and progression from mild cognitive impairment to Alzheimer's disease (Fieremans et al., 2013).
Because kurtosis metrics are more specific to several white matter features that are only indirectly reflected by typical DTI parameters, we hypothesized that they may provide better biomarkers of the relationship between white matter microstructure and physiological markers, such as VWM capacity. Thus, we used a white matter-based voxelwise approach, tract-based spatial statistics (TBSS) (Smith et al., 2006), to investigate correlations between DKI white matter metrics and VWM capacity in a group of 19 healthy male adult participants. The relationship between VWM capacity and DTI-based metrics was also examined.
Materials and Methods
Participants.
Data were acquired from male right-handed healthy controls 30–55 years of age using a protocol approved by our Institutional Review Board. Our age range was selected to minimize the effects of development or aging. Participants were recruited by advertising in the community and through ResearchMatch and Craigslist websites. All participants provided informed consent and participated in a series of cognitive assessments, including VWM. None of the participants revealed a history of brain disorders or head trauma, a history of alcohol or substance abuse in the past 6 months, or substance-induced psychotic disorder or psychotic disorder due to a general medical condition as determined by DSM-IV criteria. The Structured Clinical Interview for DSM-IV Non-Patient Edition (First et al., 2002) was used to assess the potential presence of psychiatric symptoms. Clinical neurological examination of imaging data, performed as part of the typical workflow in our department, was used to assess the presence of any gross brain abnormalities. Subjects with abnormal imaging findings or those presenting psychiatric symptoms were excluded from the study.
A total of 25 subjects were recruited. Two subjects were excluded because of the presence of depression or autism diagnoses. Three additional subjects were excluded because of brain abnormalities on neuroradiological examination, including abnormal presence of white matter hyperintensities and lesions. One subject participated in the neuropsychological testing but failed to participate in the MRI session. The remaining 19 participants were included in analyses.
The age of the included subjects varied between 31 and 55 years (42 ± 6 years). Length of education varied between 12 and 19 years (15.7 ± 1.8 years). Neither education nor the composite global cognitive score (obtained by averaging the z-scores for cognitive assessments) and none of the VWM metrics correlated with age (p > 0.80).
VWM assessment.
VWM capacity was assessed using the Symbol Span subtest of the Wechsler Memory Scale-Fourth Edition (WMS-IV). In Symbol Span, participants are shown a series of abstract symbols on a page for 5 s and then asked to select these symbols from an array of symbols, in the same order they were presented on the previous page. The number of symbols is increased from trial to trial. Each successful trial receives a 2 point score if symbols are remembered in the correct order and 1 point if all symbols are retrieved, but not in the correct order. The test is terminated after three unsuccessful trials. Symbol Span is a measure of VWM capacity as it reflects the longest list of items that a person can repeat back immediately after presentation.
Image acquisition.
Imaging data were acquired on a 3T Trio MRI (Siemens Medical Solutions). Diffusion-weighted images were acquired using Echo Planar Imaging (EPI) for a total of 64 uniformly distributed gradient directions and for two b values (b = 1000 and 2000 s/mm2). Ten nonweighted diffusion images (b = 0 s/mm2) were also collected. A field map image coplanar to the diffusion acquisition was acquired using a pair of non-EPI gradient echo images at two echo times and used to correct for image distortions from B0 field inhomogeneities.
Image processing.
Postacquisition processing was performed with locally developed software and the FSL suite (Jenkinson et al., 2012). Data preprocessing included the following: (1) visual inspection of all images for quality assessment with elimination of any artifactual images; (2) coregistration of diffusion images using the FSL Flirt utility; and (3) correction for B0 field inhomogeneities distortion using the acquired field map data and FSL.
After data correction, diffusion and kurtosis tensors were estimated using a constrained linear least square approach (Tabesh et al., 2011). The encoding gradients' matrix was adjusted to compensate for the effect of image rotation from the motion correction step. The diffusion and kurtosis tensors were used to derive 3D maps of faxon, Daxon, ADextra, and RDextra. Diffusion tensor data were used to derive the DTI metrics FA and axial (AD) and radial (RD) diffusivity at each brain voxel.
TBSS.
Voxelwise analyses were conducted using TBSS (Smith et al., 2006). First, individual FA images of all subjects were aligned to an MNI-152 template using nonlinear registration. Then, a mean FA image was created; a threshold of FA > 0.2 generated an FA skeleton that represents tracts common to all subjects. Individual parametric maps were then projected onto the mean FA skeleton by looking for maximum local values perpendicular to the skeleton. Statistical tests were conducted using the permutation-based nonparametric analysis routine Randomize (Jenkinson et al., 2012) with 5000 iterations. Correlations (positive and negative) were computed to examine the association between each of the diffusion parameters with the Symbol Span Raw Score. As age may influence both microstructural properties and Symbol Span, it was used as a covariate in all analyses to assess relationships between VWM capacity and white matter microstructure independent of potential codependence on age. Primary analyses focused on faxon, Daxon, ADextra, and RDextra metrics. Secondary analyses focused on FA, RD, and AD. Correlations were considered significant when p < 0.05, corrected for multiple comparisons using Threshold Free Cluster Enhancement (Smith and Nichols, 2009).
ROI analyses.
ROIs were obtained for the whole white matter and for tract-specific regions. As we were interested in intrinsic tract properties, we chose to focus on white matter regions of high anisotropy as these regions are less affected by fiber crossings and partial volume averaging. Neither DTI nor DKI white matter metrics adequately describes the individual fiber properties in regions where more than one fiber population is present in the same voxel. To account for this limitation, the following approximations were made: (1) regions of high linear anisotropy primarily contain single-fiber populations and thus better reflect intrinsic tract microstructure; (2) white matter skeleton regions are less affected by partial volume averaging with adjacent tracts or gray matter regions; and (3) high anisotropy white matter skeleton regions are representative of the overall intrinsic tract properties. Thus, ROIs were limited to white matter skeleton areas with linear anisotropy metric, cl (Westin et al., 2002) larger than 0.4. Linear anisotropy was selected for determining the threshold in these analyses as it is larger for prolate tensors, which are more likely to reflect single-fiber populations. Tract-specific masks were generated using the Johns Hopkins University tract atlas overlaid onto the white matter skeleton regions with cl > 0.4. A global white matter mask was also obtained by selecting all skeleton areas that met the same linear anisotropy criterion. DKI and DTI metrics for these ROIs were obtained by averaging all skeleton values covered by the masks.
In a first approach, global- and tract-specific white matter associations with VWM capacity were assessed using two-tailed partial Pearson correlations using age as a covariate. In a second approach, we tested for tract specific associations using two-stage hierarchical multiple linear regression that used Symbol Span as the dependent variable, with age and global white matter value entered as control predictors in the first stage and tract-specific metrics entered as predictors of interest in the second stage. This approach was chosen to identify tract-specific contributions to VWM capacity above global white matter contributions and independent of age. As microstructural properties are highly correlated among various white matter structures (Wahl et al., 2010), only representative tracts were included in the regression analysis to minimize overfitting from highly correlated predictors. Tract selection was performed as follows. First, hierarchical clustering was used to determine groups of tracts with similar properties. Second, forward variable selection was used to identify the best predictors within each group. Hierarchical multiple regression was performed for each of the examined metrics. As different metrics reflect at least in part different white matter properties that may each contribute independently to behavioral performance, a final model was constructed that included contributions from multiple DKI and DTI metrics. Forward selection was used to test whether either global- or tract-specific values of these metrics improved prediction of VWM capacity beyond single metric contributions. Correlations were considered significant if p < 0.01 to account for multiple comparisons. Statistical analyses were conducted using SPSS Statistics for Windows, Version 21.0 (IBM).
Results
TBSS analyses
Voxelwise TBSS results show that VWM capacity positively correlates with faxon in a wide array of white matter tracts (Fig. 1a). These tracts include two prominent pathways of the ventral and dorsal streams, the inferior and superior longitudinal fasciculi, and the optic radiations and posterior thalamus, which are known to play important roles in visual processes. Additionally, faxon of major tracts connecting the prefrontal cortex with the occipital, parietal, and temporal lobes (inferior occipitofrontal fasciculi, cingulum bundles, and uncinate fasciculi) was found to correlate to VWM capacity, which may reflect the key role of the prefrontal cortex in VWM processes. Corpus callosum, another structure prominently highlighted in the correlational maps, underlies the integration of visual information across the right and left visual fields. Aside from faxon, none of the other investigated metrics, whether DTI or DKI based, correlated significantly with VWM capacity in any brain region in TBSS analyses after accounting for multiple comparisons.
ROI analyses
ROI analyses were consistent with the TBSS analyses, with faxon the only metric correlated with VMW capacity globally (p < 0.05; Table 1); however, correlations with global ADextra were also observed at trend level (p = 0.057). VWM capacity correlated significantly with faxon of right cingulum and left superior and inferior longitudinal fasciculi, Daxon of left inferior longitudinal fasciculus, and ADextra of right cingulum. Scatter plots depicting several tract-specific associations between VWM capacity and faxon are presented in Figure 1b–e. Individual tract correlations of the different DTI and DKI diffusion metrics with VWM capacity are presented in Table 1.
Multiple regression analyses were used to assess the relative contributions of age and global- and tract-specific faxon values in predicting VWM capacity. The faxon-based hierarchical clustering of the white matter tracts selected four major clusters, with the first three clusters, including the following: (1) left inferior occipitofrontal fasciculus, anterior thalamic radiations, and genu of corpus callosum; (2) right inferior occipitofrontal fasciculus, right and left inferior longitudinal fasciculi, and splenium of the corpus callosum; and (3) superior longitudinal fasciculi and cingulum bundles (superior segment). A fourth cluster included the corticospinal tracts, parahippocampal cingulum bundles, and uncinate fasciculi. The multistep linear regression model identified faxon values of representative tracts from the first three clusters (the right cingulum, left inferior longitudinal fasciculus, and left inferior occipitofrontal fasciculus) as optimal predictors of VWM capacity controlling for age and global faxon (Table 2). Among these predictors, the individual tracts were the only ones with significant contributions (p < 0.01), with global value not contributing significantly to predicting VWM capacity (Table 2) when the individual tracts were included in the model. The age and global value model explained only 23% of the individual variability in the VWM capacity. The full model accounted for 57% of the interindividual variation in performance. Removal of global faxon, a nonsignificant predictor, increased the adjusted R2 of the overall model to 59%.
Hierarchical regression was also used to assess whether accounting for combined tract contributions in addition to global contributions improves prediction of VWM capacity for metrics other than faxon. No significant relationships (p > 0.01) were found for any of the metrics, whether for models including only age and global values as predictors or more complex models that also included representative tract contributions. Representative results are presented for FA in Table 3.
Evaluation of improvement in the faxon model by considering additional white matter microstructural features described by the other DTI and DKI metrics showed that global ADextra further improved the predictive value of the model (Table 2). The combined faxon-ADextra model explained 72% of the variance in individual VWM capacity.
Discussion
VWM capacity has been modeled in terms of the ability to link distinct cell assemblies across multiple cortical areas (Luck and Vogel, 2013). Our results strongly suggest that the communication routes among cortical areas (i.e., white matter pathways) are a key factor in determining individual VWM capacity. As faxon reflects axonal density and caliber, its relationship to VWM capacity suggests that individuals with denser white matter pathways or larger axons may be able to maintain a larger amount of information “online” for use by various network processors. From this perspective, faxon may index white matter information processing capacity. The rate of information exchange depends on axonal firing rate, which appears to be linearly related to axonal caliber (Perge et al., 2012). Thus, larger axons may contribute to faster information exchange between different network processors and increase the range of synchronized oscillations among distinguishable cell assemblies representing distinct objects stored in memory. Faster axons may also be more proficient in supporting the recurrent feedback loops that are assumed by some models to subserve VWM representations (Luck and Vogel, 2013).
The TBSS analyses delineate the extended network involved in VWM processes. Tract-specific analyses further show that tracts connecting prefrontal, occipital, parietal, and temporal lobe are significant predictors of VWM capacity, contributing beyond global white matter axon properties and independent of age. Thus, these data suggest that VWM storage and retrieval involve collaboration among multiple nodes of the VWM network; the size and density of the axons linking these nodes appear to play an essential role in determining VWM capacity. An additional, but more modest, association was found between VWM capacity and global ADextra. This finding suggests that organization and microarchitecture of extra-axonal space, which include glial cells, such as astrocytes and oligodendrocytes as well as extracellular space, may play a role in transmission of information by white matter. Glia and extracellular space have been long thought to contribute to information transmission in the brain (Fields and Stevens-Graham, 2002). The exact nature of the involved biological substrates and associated diffusion mechanisms, however, is yet to be determined.
Overall, the white matter microstructural features evaluated in this study appear to explain a substantial proportion of the individual variability in VWM capacity (72%). The remaining variance may be explained by features not described by the parameters we measured, such as contributions from gray matter, or differences in overall network organization.
One of the most intriguing properties of VWM capacity is its yet to be explained strong correlation with global cognitive performance (Luck and Vogel, 2013). The extended network of tracts we found related to VWM capacity performance are most likely also involved in more complex cognitive functions. Accordingly, we propose that reliance on a similar white matter network explains, at least in part, the high correlation between VWM capacity and measures of overall cognitive function (Luck and Vogel, 2013). Thus, white matter may represent a promising therapeutic target in disorders of both impaired VWM and impaired general cognitive abilities.
We found no significant correlations between DTI metrics and VWM capacity in our sample. This is consistent with earlier work that found no correlations between FA and cognitive metrics, including VWM capacity, in adulthood (Tamnes et al., 2011; Peters et al., 2014), although modest relationships with verbal working memory have been reported (Takeuchi et al., 2011), sometimes using only analyses restricted to one or two white matter tracts (Karlsgodt et al., 2008). Of note here is that FA correlations with various aspects of working memory (verbal or spatial) have been more widely reported in pathology, such as schizophrenia or multiple sclerosis. Such disorders have well-known impairments in FA; thus, it is expected that working memory in those cases may have different neural mechanisms than in typical healthy brain, with impaired features playing a role in modulating performance.
The DKI-based metrics used in this study are obtained assuming a two-compartment model of white matter microstructure. However, the kurtosis tensor also allows derivation of several model-independent parameters, including mean (MK), radial (RK), and axial kurtosis (AK) (Jensen and Helpern, 2010). RK was shown to be highly correlated with faxon (De Santis et al., 2012). Although we focused our main analyses on the intra- and extra-axonal diffusion metrics due to their relevance to specific white matter features, additional investigation of the relationship with RK in our data revealed, as expected, significant correlations with VWM capacity in the same regions as faxon (data not shown). Thus, the relationships between VWM capacity and white matter microstructural properties are model-independent.
In conclusion, we found that white matter represents an important anatomical substrate of VWM capacity interindividual variability in the adult healthy brain. The DKI-based metric, faxon, was found to be a strong noninvasive imaging predictor of individual visual memory capacity.
Footnotes
This study was supported by National Institutes of Mental Health Awards R21-MH085228 and R03-MH076180. We thank ResearchMatch for their support with advertising our study.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Mariana Lazar, Center for Biomedical Imaging, New York University School of Medicine, 660 First Avenue, Room 423, New York, NY 10016. mariana.lazar{at}nyumc.org