Abstract
Magnetic resonance imaging (MRI) has greatly extended the exploration of neuroplasticity in behaving animals and humans. Imaging studies recently uncovered structural changes that occur in gray and white matter, mainly after long-term training. A recent diffusion tensor imaging (DTI) study showed that training in a car racing game for 2 h induces changes in the hippocampus and parahippocampal gyri. However, the effect of short-term training on the white matter microstructure is unknown. Here we investigated the influence of short learning tasks on structural plasticity in the white matter, and specifically in the fornix, in humans and rats. Human subjects performed a 2 h spatial learning task, and rats underwent training for 1 d in a Morris water maze. Between tasks, subjects were scanned with DTI, a diffusion MRI framework sensitive to tissue microstructure. Using tract-based spatial statistics, we found changes in diffusivity indices in both humans and rats. In both species, changes in diffusion in the fornix were correlated with diffusion changes in the hippocampus, as well as with behavioral measures of improvement in the learning tasks. These results, which provide the first indication of short-term white matter plasticity in the human brain, suggest that the adult brain white matter preserves dynamic characteristics and can be modified by short-term learning experiences. The extent of change in white matter was correlated with their extent in gray matter, suggesting that all components of the neural network are capable of rapid remodeling in response to cognitive experiences.
Introduction
The nervous system possesses an intrinsic ability to learn and adapt to new experiences throughout life (Pascual-Leone et al., 2005). This neural plasticity is manifested both functionally and structurally. Processes underlying cellular mechanisms of structural plasticity include synapse formation (Knott et al., 2002), expansion and retraction of dendritic spines (Hofer et al., 2009; Xu et al., 2009), axonal remodeling (De Paola et al., 2006), astrocyte modifications (Theodosis et al., 2008), and neurogenesis (Ming and Song, 2005). Magnetic resonance imaging (MRI), by allowing in vivo exploration of whole-brain structural changes, provides invaluable information on neuroplasticity. Morphological changes observed in the adult human brains after prolonged training indicate modifications in gray matter volume and cortical thickness (Maguire et al., 2000; Draganski et al., 2004, 2006). Diffusion tensor imaging (DTI), a diffusion MRI framework, is sensitive to the tissue microstructure and yields quantitative indices reflecting micron-scale tissue density and organization (Assaf and Pasternak, 2008).
Morphometric MRI analyses were complemented by DTI studies of nervous system plasticity showing microstructural differences between nerve fiber tracts of experts and controls (Bengtsson et al., 2005; Imfeld et al., 2009) or before and after long-term training (weeks to months) (Scholz et al., 2009; Takeuchi et al., 2010). However, studies of the short-term training effect on tissue microstructure are scarce. In one recent study, changes in diffusion indices in gray matter were detected after only 2 h of training, pointing to a rapid timescale of structural plasticity (Sagi et al., 2012). However, it is still unknown whether short learning experiences can also modify adult brain white matter. Recent in vitro studies provided evidence for an effect, after only a few hours, of electrical activity on myelination (Wake et al., 2011), as well as for oligodendrocyte development and differentiation into myelinating oligodendrocytes (Stevens et al., 2002).
Here we used DTI to examine structural changes in white matter tracts after performance of short spatial-learning tasks by humans and rats. Human subjects were trained for 2 h in a car racing game, a task found previously to induce changes in diffusion in the hippocampus and parahippocampal gyrus (Sagi et al., 2012). The rat study was based on a short version of Morris water-maze training. Performance in a maze engages the fimbria–fornix (Sutherland and Rodriguez, 1989; de Bruin et al., 2001), and training induces changes in the hippocampal formation (Blumenfeld-Katzir et al., 2011; Lerch et al., 2011). Our analysis was restricted to the fornix, which is the main tract connecting the hippocampus to the medial diencephalon, an area suggested to be involved in memory processes (Gaffan and Gaffan, 1991; Rudebeck et al., 2009). Because the alignment of the whole-brain images used in the above studies was not sufficient to analyze white matter across subjects, we applied tract-based spatial statistics (TBSS), a method used commonly to study white matter (Smith et al., 2006).
Materials and Methods
Human study
Subjects.
The study participants were 70 adult right-handed volunteers (39 women and 41 men; mean ± SD age, 26.6 ± 3.8 years; range, 20–38 years). The subjects were divided into three experimental groups: (1) a learning group (LG; n = 40, 17 males and 23 females; mean age, 26.7 years); (2) a nonspatial learning group that served as an active control (AC; n = 15, 7 males and 8 females; mean age, 26.9 years); and (3) a passive control group (PC; n = 15, 7 males and 8 females; mean age, 26.1 years). This dataset includes subjects from the study by Sagi et al. (2012), as well as a new cohort of subjects. The research protocol was approved by the Institutional Review Board of the Tel Aviv Sourasky Medical Center. Each participant signed a written informed consent form. None of them had a history of neurological disease, psychological disorders, drug or alcohol abuse, or use of neuropsychiatric medication. All had intact vision.
Behavioral paradigm.
Subjects in the LG and AC groups were trained to play a computer car racing game (The Need for Speed; Electronic Arts). The LG training consisted of completing 16 laps (trials) divided into four training sessions (four laps each). The same track was presented in each trial. The aim was to learn the track and achieve better lap times. At the end of each training session, subjects were asked to sketch an outline of the track and were given snapshots of locations along the track that they had to arrange in the correct order. In the AC group, subjects played the same car racing game for the same duration as the LG subjects, but in each lap (trial) the track was different. Subjects in the PC group, like those in the other two groups, underwent two MRI scan, 2 h apart, but unlike the others, they did not participate in any training tasks between the imaging sessions.
Imaging data.
Data acquisition was performed at the Tel Aviv Sourasky Medical Center with a 3 T MRI system (GE Healthcare). Each subject underwent two imaging sessions, ∼2 h apart. For both the LG and AC subjects, these scanning sessions took place before and immediately after performance of their assigned task. The MRI protocol of the first scanning session included both conventional anatomy sequences and DTI acquired with an eight-channel head coil. In the second scanning session, only DTI scans were administered.
The DTI protocol included acquisition of whole-brain diffusion-weighted images with the following parameters: 2.1 × 2.1 × 2.1 mm3 resolution reconstructed to 1.58 × 1.58 × 2.1 mm3, up to 70 axial slices, Δ/δ = 33/26 ms, b = 1000 s/mm2, with 19 gradient directions and an additional image with no diffusion weighting (b0 image). To increase the signal/noise ratio, the DTI scan was repeated three times.
To correct for head motion, we used a least-squares algorithm and six-parameter (rigid-body) transformations using SPM (version 2; University College London). DTI was analyzed by the DiVa (diffusion imaging visualization and analysis) software (Pasternak et al., 2008, 2009) implemented in MATLAB 7.3.0 (MathWorks), from which maps of fractional anisotropy (FA), mean diffusivity (MD), principal (axial) diffusivity (λ1), and radial diffusivity (RD) were computed.
Analysis of DTI data.
The data were subjected to voxelwise statistical analysis by TBSS (Smith et al., 2006), part of the FSL [for FMRIB Software Library (in which FMRIB is Functional MRI of the Brain, Oxford University, Oxford, UK)] package (Smith et al., 2004). Briefly, FA images from both scans for each subject were first nonlinearly aligned into to a 1 × 1 × 1 mm common space using the nonlinear image registration tool of FMRIB. Two LG subjects were excluded because their ventricles were enlarged, causing misalignment of the fornix tract after nonlinear normalization. The transformed maps were averaged to generate a mean FA image, which was then “thinned” to create a white matter FA skeleton at a threshold FA value of 0.2. The aligned FA image of each subject was then projected onto the mean FA skeleton. The MD, λ1, and RD maps were analyzed by application of the same spatial transformation parameters as those estimated in the initial FA analysis. The course of the fornix in the ventricles makes it prone to contamination by CSF. To minimize artifacts resulting from partial volume effects, we excluded the data from subjects in whom >10% of the voxels in the fornix skeleton had MD values, indicating CSF contamination (MD > 2). Accordingly, five subjects were excluded from the LG, two from the AC, and six from the PC groups.
In view of our hypothesis that the fornix plays a role in this short-term learning, we used anatomical atlases to create a mean fornix mask for statistical analysis. Our focus on analyzing the fornix tract instead of the whole-brain skeleton also helped to reduce multiple comparisons issues and to limit false-positive errors.
Rat study
Behavioral paradigm.
For this part of the study, 44 male Wistar rats (4 months old) underwent two MRI scans during the light phase, 1 d apart. All rats were maintained on a 12 h light/dark cycle with access to food and water ad libitum. Before scanning, the rats were anesthetized with 1–2% isoflurane in oxygen.
The rats were divided into three groups: (1) a learning group (L; n = 24); (2) a cued (visible platform) swim group (C; n = 10); and (3) a passive (naive) control group (N; n = 10). This dataset was included in the study by Sagi et al. (2012). Between scans, the rats in group L were trained to perform a short water-maze task in a 120-cm-diameter pool. Training consisted of three repeats of four trials, with ∼40 min between trials. At each trial, the rat was placed in a different quadrant of the pool and had to find, based on spatial cues, the location of a hidden platform. An HVS video-tracking system (HVS Image) was used to record the time taken to reach the platform, the time spent in each quadrant, and the swimming speed. Rats in group C were given the same task, but the platform was made visible by an object placed on top of it. The rats in group N remained in their cages between the two scans.
Imaging data.
MRI scanning was performed using a 7 T MRI scanner (Bruker) with a 30 cm bore and a gradient strength of up to 400 mT/m. The MRI protocol included (DTI) acquisition with a diffusion-weighted spin-echo echo-planar-imaging pulse sequence. Acquired volumes were 21 slices, each 0.75 mm thick with the following parameters: resolution of 0.2 × 0.2 mm2 (matrix size, 128 × 128; field of view, 25.6 mm2), repetition time of 4000 ms, echo time of 25 ms, Δ/δ = 10/4.5 ms, four echo-planar-imaging segments, and 15 noncollinear gradient directions with a single b-value shell at 1000 s/mm2 and one image with a b value of 0 s/mm2 (b0 image). The imaging protocol was repeated three times for signal averaging and to compensate for acquisition in which significant head motion was observed. Each DTI acquisition took 4.5 min, and the entire MRI protocol lasted ∼20 min.
Because excessive breathing during DTI acquisition can lead to motion artifacts in the sampled slice, images of every slice at each gradient direction were screened automatically for motion artifacts. If the signal in the outer borders of the brain was significant, the slice was omitted from analysis. DTI acquisition was repeated three times to compensate for the extracted slices. DTI analysis was performed with in-house software, implemented on MATLAB, and FA and MD maps were produced.
Image normalization was performed using SPM2. Each rat was normalized to a rat atlas template that was created from one representative rat and included b0 and FA images registered with a digitized version of the Paxinos and Watson (2005) stereotaxic atlas. The normalization procedure initially included bias correction of all b0 images and coregistration of each b0 image with the b0 template using a six-parameter rigid-body transformation. The parameters of the coregistration were applied on the different FA and MD maps. The coregistered FA maps were then normalized to the FA template using a 12-parameter affine nonlinear transformation and 0.2 mm smoothing. Normalization parameters were applied on all DTI-indexed maps (FA, MD, λ1, and RD). After the normalization, data were reexamined, and some rats were excluded because of motion artifacts, asymmetrical ventricles, or Nyquist ghost artifacts in relevant brain areas. Five rats were omitted from group L, three from group C, and four from group N.
A mean FA map was created from the normalized images, and a skeleton was generated from it using the TBSS procedure, with a threshold at FA ≥ 0.1. Rat FA data were projected onto the skeleton. Projection of MD, λ1, and RD values to the skeleton was also applied. Statistical analysis was restricted to the fornix by use of a mean fornix mask that was created using the Paxinos and Watson (2005) stereotaxic atlas.
Statistical analyses.
To ascertain fair comparison across groups and to allow the replication of results, the learning groups in both the human and rat studies were divided into two subgroups. Division of the learning groups in this way helped to reduce the difference in sample size between the learning groups and the control and provided an additional cohort of subjects that could be independently analyzed while also serving as a replication group for the learning effect. Such replication strengthens the inference derived from the results of the statistical analysis.
The learning group in the human study was divided on the basis of the analysis described by Sagi et al. (2012). The first learning cohort (LG1) was compared with both control groups (LG1, n = 14; AC, n = 13; PC, n = 9). The second learning cohort (termed LG2 in this study; n = 19) comprised the group used by Sagi et al. (2012) for replication with additional new subjects who underwent the same experimental procedure.
The learning group in the rats study was randomly divided into two. One group (L1) was compared with the C and N groups (L1, n = 10; C, n = 7; N, n = 6), and the other was used for replication (L2, n = 9). We report on changes in diffusion indices that pass p < 0.01 in the across groups analyses (learning vs controls), as well as in the replication group.
Data both in the human and rat study were statistically analyzed using the Randomize program (part of FSL), carrying voxel-by-voxel nonparametric permutation test (5000 permutations) and thresholded using threshold-free cluster enhancement, which obviates the use of an arbitrary threshold for the initial cluster formation (Smith and Nichols, 2009). General linear model designs were applied for analysis of group contrasts × time interaction, by paired t test analysis for the replication groups, and by linear correlations with changes in gray matter (hippocampus) and behavioral measurements.
To examine whether changes in the fornix were correlated with changes in the hippocampus, we conducted paired t tests on the normalized MD maps of the learning groups of both humans and rats, using SPM2 (for the normalization process on the human images, see Sagi et al., 2012). Mean values of MD changes in the cluster of results found in the hippocampus were extracted for each subject, and the percentage of change in MD was calculated. These values were applied in a voxel-based correlation analysis with the differences (post-task vs pre-task) in diffusion indices in the fornix.
Results
Behavioral results
In the human study, improvement in the learning task was demonstrated by all LG subjects, as indicated by a reduction of 22.5 ± 3% (mean ± SEM, p < 0.0001; n = 33) in their normalized lap times and an improvement in their scores for picture arrangement (p < 0.0001; Fig. 1A). In contrast, the AC group showed no improvement in their normalized lap time (Fig. 1A). In the rat study, latency time of the L group rats in the water maze were significantly reduced (57.4 ± 9.4%, mean ± SEM, p < 0.0001; n = 19; Fig. 1B).
Short-term structural changes in the fornix
To assess the effect of learning on structural changes in the fornix, we examined group (learning vs two control groups) × time (scan time) interaction. We did this by looking for clusters in which changes in diffusion indices were significantly different in the LG than in both of the control groups. Significant differences in MD were found in two clusters in the body of the fornix (p < 0.01 cluster ≥ 5, uncorrected; at most significant voxel, p = 0.001). In these clusters, MD in LG1 was decreased by an average ± SEM of 6.75 ± 1.64% (Fig. 2A,C). Reduction in MD was also found in LG2 in the column of the fornix (p < 0.05 cluster ≥ 5, uncorrected; at most significant voxel, p = 0.006; Fig. 2A). To confirm that voxels showing significant effect were indeed located within the tract, significant clusters were de-projected (using an FSL de-project tool) back to each subject's standard-space (nonlinearly registered) FA images.
A decrease in FA was found in the analysis across the three groups. However, the same change was not observed at the same statistical threshold in the LG2, in which a decrease was found only at p < 0.05. However, additional analysis revealed that the observed reduction in MD derived from reduction in the λ1 (p < 0.01 cluster ≥ 5; Fig. 2B) but not from changes in the RD. Reduction in λ1 was detected in the same regions as those showing reductions in MD.
In the rat study, analyses of the three groups (L vs C and N) × time interaction disclosed a similar pattern of change in MD as was found in the human study. An average ± SEM decrease of 5 ± 1.08% in MD was found in group L1 (Fig. 2D,F). This reduction was also observed in the replication group (L2) (p < 0.01 cluster ≥ 4, uncorrected; at most significant voxel, p = 0.002; Fig. 2D). Additional analysis revealed that the decrease in MD originates from a reduction in λ1 (p < 0.01 cluster ≥ 3; Fig. 2E). In both analyses, a decrease in FA was found at less strict p values of p < 0.05.
Correlation with structural changes in the gray matter of the hippocampus
Previous whole-brain analyses of the learning groups in both human and rat studies revealed significant reductions in MD in the hippocampus (Sagi et al., 2012). In the present study, we found that this change in gray matter in the hippocampus is correlated with the difference between pre-task and post-task diffusion indices in the fornix. In the human study, changes in MD values in the hippocampus was correlated with the difference in MD values in the fornix (Fig. 3A,B; p < 0.05, corrected; r = 0.62), as well as with the difference in pre-task and post-task λ1 values in the fornix (Fig. 3A,C,D; p < 0.05, corrected; r = 0.87). Correlation between change in MD in the hippocampus and change in λ1 in the fornix was found in the left and right crus of the fornix, as well as in the body of the fornix in which we found a significant reduction in MD. In the rat study, change in MD in the hippocampus was also correlated with the differences between MD values in the fornix (Fig. 4A,B; p < 0.05, corrected; r = 0.9) and with the difference between RD values in the fornix (Fig. 4C,D; p < 0.05, corrected; r = 0.87).
Correlation with behavioral data
We also analyzed the correlation between behavioral measures and differences in diffusion indices. The rate of improvement in the car racing game was found to correlate with reductions in MD (p < 0.005 cluster ≥ 5; r = 0.55) and FA (Fig. 5A,B; p < 0.005 cluster ≥ 5; r = −0.74) in the column and right crus of the fornix. Significant correlations were also found between the averaged picture scores and differences in MD (Fig. 5C,D; p < 0.005 cluster ≥ 5; r = 0.52) and λ1 (p < 0.005 cluster ≥ 5; r = 0.53) in the body of the fornix.
One measure of improvement in the water-maze paradigm is a reduction in the time required for the rat to reach the platform. This learning measure was found to be correlated with changes in MD in the fornix (p < 0.005 cluster ≥ 5; r = 0.58), as well as with RD differences in the fornix (Fig. 5E,F; p < 0.005 cluster ≥ 5; r = 0.78).
Discussion
Our results provide the first indication of short-term white matter plasticity in the human brain. Using diffusion MRI, we showed that short periods of spatial learning induce changes in diffusion indices in the fornix. We found that these changes exist across species and are related to coextending changes in diffusion in gray matter and correlated with behavioral learning progression and improvement.
The hypothesis put forward in this study was that microstructural changes occur in the white matter after short learning experiences and that the rapidly achieved structural plasticity can be detected by DTI. This imaging technique yields quantitative measurements of white matter microstructure and can therefore be used to study plasticity in vivo (Scholz et al., 2009; Blumenfeld-Katzir et al., 2011; Lerch et al., 2011). Studies have indeed demonstrated anatomical differences in the white matter of groups of experts, such as musicians (Imfeld et al., 2009), as well as correlations between fiber tract organization and the amount of piano practice (Bengtsson et al., 2005), and changes in diffusion after training in juggling (Scholz et al., 2009) or after meditating (Tang et al., 2010). Our study further affirms the utility of diffusion imaging in neuroplasticity research and presents new evidence for the temporal progression of structural plasticity in white matter within short periods of time. In agreement with studies such as those mentioned above, we suggest that white matter in the adult brain preserves dynamic characteristics. It seems clear that biochemical and electrophysiological changes (see below and Wake et al., 2011) can induce rapid structural remodeling in the white matter in the form and content of myelin and other cellular structures.
The findings on changes in diffusion in the fornix support and complement a previous research on short-term plasticity that revealed a decrease in MD in the human hippocampus and parahippocampal gyrus after the same spatial learning task as the one used in the present work. We showed here that those gray matter modifications correlate with structural changes in the fornix, the main hippocampal projection. A similar strong correlation between changes in MD in the fornix and changes in the hippocampus was found in our rat study. Together, these findings raise the possibility that microstructural alterations occur in related gray and white matter and that the extent of change they undergo is associated. These correlations further point to a connected physiology of the neural network, in which all of its components interact and become modified by new experiences. Although MRI does not reveal the starting point of these changes or indicate whether changes in one component generate changes in others, it does suggest that a comprehensive study of neuroplasticity should incorporate scrutiny of both gray and white matter. Furthermore, changes in diffusion indices in the fornix, in both humans and rats, were found to correlate with improvement in the learning procedures. Therefore, it seems that our aptitude to learn is connected to the biological processes underling learning progression. In addition, the observed correlation between the mean picture scores that subjects achieved in the car racing task and the changes in diffusion in the fornix supports a role for the fornix in memory and recollection processes (Rudebeck et al., 2009).
It should be noted that the results showing a decrease in MD and λ1 in the fornix in learning task participants did not pass multiple comparisons. Structural changes in white matter induced over a short periods (hours) would not be expected to produce extreme modifications in the tissue and would be expected, as a consequence, to be buried in noise. However, the same results showing a decrease in MD emerged from the two separate analyses that were conducted in both the human and rat studies. In addition, results of the correlation analyses further support a role for the fornix in learning-induced plasticity. Aside from demonstrating rapid attainment of plasticity of the white matter, this work underscores the advantages of using a tract-based approach in the investigation of white matter. Whole-brain voxel-based analysis without track registration failed to provide adequate alignment of the relevant tracts to enable proper voxel-based comparisons within and between subjects.
Interpretation of the changes in diffusion in white matter is not a straightforward matter. MD and FA are generally believed to depict tissue density and fiber organization, respectively. Reduction in MD might be attributed to changes in extracellular volume (Ransom et al., 1985; Sykova, 1997), swelling of cells (Le Bihan, 2007), or an increase in glia cell volume (Kleim et al., 2007; Theodosis et al., 2008; Markham et al., 2009). Electrical activity might also induce structural changes in axons. Swelling of axons was found in association with action potentials (Iwasa and Tasaki, 1980; Iwasa et al., 1980; Tasaki and Byrne, 1990). Activated axons transmit signals to neighboring glial cells, thereby promoting myelination (Ishibashi et al., 2006; Wake et al., 2011). It could be postulated that this axon-to-glia signaling promote changes in glial cell morphology. Remodeling of astrocytic processes was found to occur after training in the water maze (Blumenfeld-Katzir et al., 2011; Lerch et al., 2011; Sagi et al., 2012). In addition, oligodendrocyte maturation and morphology were found to be altered as a result of electrical activity in axons (Stevens et al., 2002). Changes in cell size and the numerous cellular processes in glia may result in a decrease in λ1. Additional research is needed to reveal the biological processes that underlie learning-induced plasticity in white matter and their relationship to diffusion parameters. Although DTI represents a state-of-the-art method for studying white matter, alternative techniques using diffusion MRI were recently suggested. Methods based on models such as CHARMED (for composite hindered and restricted model of diffusion) (Assaf et al., 2004), AxCaliber (Assaf et al., 2008), and ActiveAx (Alexander et al., 2010) have the advantage of quantifying the contribution of different tissue components to the signal (e.g., axonal compartment and axonal composition). These methods allow direct examination of specific tissue components. This may be helpful in determining the biological origins of white matter plasticity and should be applied in future studies.
To summarize, this study revealed changes in diffusion in the fornix of adult humans and rats after only 2 h of training in a learning task, as well as marked correlations between behavior and structural changes in the fornix origin, the hippocampus. Because diffusion indices are indicative of tissues microstructure, the results suggest that short-term learning can induce structural plasticity in the white matter. Additional research is needed to understand the cellular processes that underlie these changes, as well as the time course of their development. This study suggests that the multiple components of the neural network constitute a dynamic system that is able to achieve rapid remodeling as a result of cognitive experiences. An additional striking finding is the ability of MRI to map and quantify micron-scale processes in the neural tissue, offering the opportunity to study the physiology of the human brain in different timescales and under various spatial conditions.
Footnotes
This work was supported by Israel Science Foundation Grant 994/08 and the Future and Emerging Technologies Programme within the Seventh Framework Programme for Research of the European Commission (FET-Open “CONNECT” Project Grant 238292).
The authors declare no competing financial interests.
- Correspondence should be addressed to Prof. Yaniv Assaf, Department of Neurobiology, George S. Wise Faculty of Life Science, Tel Aviv University, 69978 Tel Aviv, Israel. assafyan{at}post.tau.ac.il