Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Semi-Automatic Segmentation of Optic Radiations and LGN, and Their Relationship to EEG Alpha Waves

  • Emmanuelle Renauld,

    Affiliation Department of Nuclear Medecine and Radiobiology, Faculty of Medicine and Health Science, University of Sherbrooke, Sherbrooke, Qc, Canada

  • Maxime Descoteaux,

    Affiliations Sherbrooke Connectivity Imaging Lab (SCIL), Computer Science Department, Faculty of Science, University of Sherbrooke, Sherbrooke, Qc, Canada, Centre Hospitalier Universitaire de Sherbrooke (CHUS), Sherbrooke, Qc, Canada, Centre d’Imagerie Moléculaire de Sherbrooke (CIMS), Centre de Recherche du CHUS, Sherbrooke, Qc, Canada

  • Michaël Bernier,

    Affiliation Department of Nuclear Medecine and Radiobiology, Faculty of Medicine and Health Science, University of Sherbrooke, Sherbrooke, Qc, Canada

  • Eleftherios Garyfallidis,

    Affiliation Sherbrooke Connectivity Imaging Lab (SCIL), Computer Science Department, Faculty of Science, University of Sherbrooke, Sherbrooke, Qc, Canada

  • Kevin Whittingstall

    kevin.whittingstall@usherbrooke.ca

    Affiliations Department of Nuclear Medecine and Radiobiology, Faculty of Medicine and Health Science, University of Sherbrooke, Sherbrooke, Qc, Canada, Department of Diagnostic Radiology, Faculty of Medicine and Health Science, University of Sherbrooke, Sherbrooke, Qc, Canada, Centre Hospitalier Universitaire de Sherbrooke (CHUS), Sherbrooke, Qc, Canada, Centre d’Imagerie Moléculaire de Sherbrooke (CIMS), Centre de Recherche du CHUS, Sherbrooke, Qc, Canada

Abstract

At rest, healthy human brain activity is characterized by large electroencephalography (EEG) fluctuations in the 8-13 Hz range, commonly referred to as the alpha band. Although it is well known that EEG alpha activity varies across individuals, few studies have investigated how this may be related to underlying morphological variations in brain structure. Specifically, it is generally believed that the lateral geniculate nucleus (LGN) and its efferent fibres (optic radiation, OR) play a key role in alpha activity, yet it is unclear whether their shape or size variations contribute to its inter-subject variability. Given the widespread use of EEG alpha in basic and clinical research, addressing this is important, though difficult given the problems associated with reliably segmenting the LGN and OR. For this, we employed a multi-modal approach and combined diffusion magnetic resonance imaging (dMRI), functional magnetic resonance imaging (fMRI) and EEG in 20 healthy subjects to measure structure and function, respectively. For the former, we developed a new, semi-automated approach for segmenting the OR and LGN, from which we extracted several structural metrics such as volume, position and diffusivity. Although these measures corresponded well with known morphology based on previous post-mortem studies, we nonetheless found that their inter-subject variability was not significantly correlated to alpha power or peak frequency (p >0.05). Our results therefore suggest that alpha variability may be mediated by an alternative structural source and our proposed methodology may in general help in better understanding the influence of anatomy on function such as measured by EEG or fMRI.

Introduction

The relationship between structure and function has long been debated. Is a bigger brain more active? Is the shape of certain structures important? The advent of medical imaging, in particular magnetic resonance imaging (MRI) and even more diffusion (dMRI), has substantially transformed the way researchers address these questions. Structure and function cannot be separated from one another, and we can now investigate the relation between the way the brain is divided into substructures (segregation) and the way these areas communicate with one another (integration) [1].

Technical advances in the field of MRI and dMRI have enabled researchers to investigate how voxel-scale metrics—such as gray (GM) and white (WM) matter density, fractional anisotropy (FA), mean, radial or axial diffusivity and axonal density—are related to behavioral measures. In healthy subjects, for example, this has been addressed by comparing how such structural variations in the same age-group are related to memory, visual or motor skills [24]. However, few studies have analyzed the effect of more macro-structural information such as size, shape or position of WM and GM structures in healthy subjects (for a complete review, see Kanai and Rees [5]). Indeed, this topic is often addressed in healthy aging [6] or in patients with neurological disorders, where changes in shape may be more obvious, for example due to malformations. Here, we consider structural and functional variations across healthy subjects.

A particularly interesting case of functional variability is that derived from human electroencephalography (EEG) recordings. The alpha rhythm, the most prominent EEG rhythm, was first defined by Berger (1930) as EEG oscillations between 8 and 13 Hz, visible from the occipital cortex when the subject is awake with eyes closed and attenuated with eyes opening. In general, EEG alpha variability has been shown to be stronger across subjects (inter-subject variability) compared to repeated measures in the same subject (intra-subject variability) [7]. Yet, despite decades of research, questions regarding their raison d’être and origin are still debated. There have been multiple studies investigating alpha variability in healthy subjects with non-anatomical factors, such as cognitive behaviour, IQ, age and genetics (see Bazanova and Vernon [7] for a review) but, again, few have analyzed the effect of anatomical factors, such as the morphology of GM and WM components in the brain (size, shape, position, etc.).

Structures especially worth investigating are the thalamus, in particular the visual thalamus, and the associated white matter tracts, i.e. the optic tract (OT) and the optic radiation (OR), connecting the retina and visual cortex to the lateral geniculate nucleus (LGN) of the thalamus. Indeed, early microelectrode recordings in animal models showed that neural activity in the LGN is related to cortical alpha activity [8]. Later work in humans showed that spontaneous fluctuations in alpha power are correlated to thalamic blood oxygen level dependent (BOLD) functional MRI (fMRI) [9] and glucose metabolism derived Positron Emission Tomography (PET) (see Liu et al. [10] for a review). Also, knowing that the anatomy of the thalamus and its nuclei influence other phenomenon or diseases, many of which are related to vision (amblyopy, albinism, etc. [11, 12]), it is therefore plausible that such structural variations affect alpha waves, also vision-related. Andrews et al. [13] showed that, although the visual cortex, optic tract and LGN volumes are correlated together and vary between subjects by 100-200%, they are however independent from brain size, which varies by 18%. In terms of the LGN fiber pathways (OR and OT), it has been suggested that it could also influence alpha activity. On the one hand, increasing myelination increases the speed and efficacy of signal conduction, which could in turn promote neural synchrony at its cortical targets, thus influencing EEG measurements at the scalp [14]. On the other hand, variations in OR/OT volume might result in strong volume-conduction effects due to anisotropy [15], while fiber length and/or density could also affect scalp EEG measurements [16]. Yet, little work has analyzed the effect of thalamic structures on EEG alpha waves, which is surprising considering that most of the EEG rhythms generation models using the thalamus use such anatomical parameters [17], and considering that it has been proposed that effect of thalamic and OR variation should be considered [17, 18]. Nunez et al. [19] initially proposed that peak alpha frequency would be related to head size, though this result was recently challenged by Valdés-Hernándes et al. [17], who suggested that FA in sections of thalamic-related fibers was a better alpha correlate. However, given the well-known limitations of the tensor model (from which FA is measured) when quantifying WM, it is unclear whether such results hold when using reconstruction methods that are less sensitive to this limitation, such as High Angular Resolution Diffusion Imaging (HARDI) [20]. To our knowledge, only one study has analyzed the effect of thalamic GM/WM morphology on EEG: Hindriks et al. [21] used DTI connection strength and biophysical models to analyze the importance of the OR. However, connection strength uses the number of streamlines produced by tractography, which might not always accurately describe WM tracks [22]. Also, they could not analyze the effect of the LGN, a satisfactory segmentation of the LGN being difficult to obtain. We believe HARDI-based tractography probabilistic tracking handles the uncertainty in the principal directions extracted from the fiber orientation distribution functions (fODFs) and allows for bundle reconstruction with a higher spatial extent or volumetric coverage. Therefore, in the current study, we combine EEG, fMRI and HARDI-based dMRI to improve quantification of the OR and segmentation of the LGN and to investigate whether structural inter-subject variations in the LGN and its efferent fibers (OR) explain alpha variability in healthy humans.

Background

Segmentation of the thalamus.

The thalamus represents a mixture of white and gray matter, and thus has proven difficult to subdivide from structural images derived from T1 or T2-weighted MR sequences [36, 37], although specific techniques have been proposed [3843]. Thalamic atlases exist, but, given the inter-subject variability of the thalamus and its nuclei [13], this may lead to inaccuracies. Diffusion MRI (dMRI) on the other hand may allow more precision in thalamic segmentation. The founding axiom of the thalamic nuclei segmentation algorithms based on dMRI is that different nuclei have different functions, and thus, different connections to specific parts of the brain. As a result, one could determine a nucleus by clustering the voxels belonging to an associated fiber [44, 45]. Indeed, several studies have used the diffusion tensor (DT) of each voxel to distinguish nuclei from one another without necessarily computing their fibers with tractography [44, 4651], though this approach has drawbacks. There is the problem that DT precision is limited in areas with a high density of gray matter or in crossing fiber regions, such as the thalamus [24]. HARDI could help, but such algorithms have be shown to be less stable as they deal with more information [52]. Also, these clustering techniques make it difficult labelling nuclei automatically and often require visual expertise.

It then seems more promising to use tractography because we have more knowledge about the expected features of the fibers, such as their starting and ending points or their shape. The first articles proposing tractography for the segmentation of the thalamus used the structural connectivity to specific regions of the cortex [36, 45, 5357], labeling the voxels with the name of the associated cortical region. However, first, instead of a thalamus segmentation problem, we then have to face a cortex segmentation one. Also, assessing the quality of the tractography is arduous, though some streamline probability scores have been proposed [32]. That problem will be even more important as HARDI probabilistic tracking becomes more popular, giving more precise streamlines, but along with more false positives [36, 58]. Also, these probabilistic metrics need an additional supervised thresholding step to ensure a valid segmentation. We therefore use the structural connectivity through one well known fiber, the OR, to segment the LGN.

Segmentation of the optic radiation.

The OR has been studied intensively (see Table 1), mainly because of the high risk of vision loss after damages to the OR during temporal lobe epilepsy surgery or tumor resections. The OR originates from the LGN and leads to the occipital lobe after looping anteriorly from the LGN into the temporal lobe, around the temporal horn of the lateral ventricle (Meyer’s loop). It is generally believed that the bulk of the OR leads to the calcarine fissure (in the primary visual cortex, V1), though it has been shown recently that it may reach V2 and V3 as well [59].

OR tractography often depends on a roughly defined LGN from a priori information, a region of interest based on the fractional anisotropy (FA) or on the principal diffusion direction (PDD) of the DT, or manually segmented (see Table 1). However, for this study, we have two objectives: the analysis of the OR and the analysis of the LGN. We are looking for the opposite of what has been published: we want to obtain the LGN from a good segmentation of the OR using tractography. We develop a nearly-automated fiber analysis tool that we believe gives reliable results, based on fiber comparison to a template.

With new information on the size and shape of the thalamus, the optic radiation and the associated nuclei, the LGN, we address the question of the influence of structure on EEG alpha waves. We show that our semi-automatic method using probabilistic tracking allows obtaining good optic radiations. We also show that in this case, structure alone is not sufficient to explain variability in the alpha power or frequency.

Methods

Subjects

Twenty-three healthy subjects were recruited for the study (12 males). All subjects were native French speakers with no psychiatric or neurological symptoms. The study was approved according to the guidelines of the Internal Review Board of the Centre Hospitalier Universitaire de Sherbrooke (CHUS) and all participants provided written consent. Data from two subjects was discarded due to excessive movement during MR acquisition. Another subject showed no EEG alpha peak, and was also excluded.

Data acquisition

MR Imaging data were acquired using a 1.5 T SIEMENS Magnetom (Vision). Noise-reduction headphones and head cushions were used to minimize artifacts. Each session started with an anatomical T1-weighted 1mm isotropic MPRAGE (TR/TE 1860/3.54 msec) acquisition, followed by a fMRI protocol and finally, with a dMRI acquisition. Functional MRI datasets were collected using a standard echo-planar imaging (EPI) sequence: 35 axial image slices, 64 x 64 matrix, TR/TE 2730/40 msec, voxel size 3.438 x 3.438 x 4.2mm. Data were acquired in a box-car format, with subjects alternating between baseline and task conditions via short auditory cues. Subjects performed alternating epochs of 30 sec rest, eyes closed, and 20 sec eyes open, repeated 5 times and ending with a rest epoch, resulting in a total scan time of 4 min and 40 sec. Diffusion MRI datasets were acquired using a single-shot echo-planar (EPI) spin echo sequence of 12 minutes (TR/TE = 11700/98 ms), with b-value of 1000 s/mm2, 64 uniform directions, matrix size of 128x128, 2mm isotropic spatial resolution. To reduce susceptibility distortions, GRAPPA parallel imaging was employed with an acceleration factor of 2.

EEG data were acquired on the same day, for each subject, using a 64-channel EEG system (Brain Products) with sampling rate of 500Hz, according to the 10-20 system, in an eyes closed-eyes open (EC-EO) task guided by an auditory cue (30 seconds EC, 20 seconds EO, for a total recording time of 5 minutes).

Data preprocessing

All preprocessing and processing steps are reflected in Fig 1.

thumbnail
Fig 1. Analysis pipeline describing dMRI, fMRI and EEG analysis.

https://doi.org/10.1371/journal.pone.0156436.g001

T1-weighted images.

T1 images were denoised using non-local means (NLM) [60] (Dipy, [61]) and skull-stripped using FSL [62]. T1-weighted images were registered to the upsampled b0 (1x1x1 mm); linear registrations were first performed with FSL FLIRT [63], followed by non-linear registration with ANTS [64]. Non-linear registration was useful even though it is an intra-subject registration to compensate for dMRI deformations.

Diffusion-weighted images.

Diffusion images were denoised using NLM with Rician distribution adaptation [65] available in Dipy. Datasets were upsampled to 1mm isotropic resolution with trilinear interpolation. The b0 image was extracted using AFNI [66] and skulled-stripped using FSL [67]. Tensors, fractional anisotropy (FA) map and fODFs were computed using Dipy.

Functional images.

The description of all fMRI analysis can be found in [68]: We processed the data using slice timing and motion correction, band-pass temporal filtering (0.008 to 0.1 hz) and detrending. The spatial smoothing was replaced by NLM denoising (Dipy). The activation maps were generated by computing the Pearson correlation coefficient between the hemodynamic response function (HRF) convolved stimulus box-car and the voxel time-series. The activation were then converted to z-scores and thresholded at the 96th percentile (equivalent to a z-score of ∼2.5). We then registered the activation maps to the T1-weighted image using the same procedure as described previously (FLIRT+ANTS).

EEG data.

EEG preprocessing was performed using EEGLAB [69]. EEG signals were high-passed filtered at 0.5 Hz. ICA was computed, and bad components were identified as follows. First, the z-score of the power spectral density (PSD) between 5 and 18 Hz was computed. If no point in the alpha band (8-13 Hz) had a z-score higher than 1 (i.e. the component has no alpha peak), the component was removed. Second, the z-score of the weights of the component on each electrode in the mixing matrix was computed. If a component had a z-score higher than 4, it was removed. Thresholds 1 and 4 were chosen after visual analysis of the data. Further details can be found in [70]. For three subjects, one or two channels were interpolated (due to excessive noise) prior to ICA decomposition. Data was rereferenced to an average reference.

Data processing

Brain structures.

The size of each brain was measured as proposed by FSL-FAST (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FAST). Thalamic masks and occipital lobe masks were obtained using FSL FIRST [71] and FreeSurfer (http://freesurfer.net/), respectively. Thalamic masks were further modified by constraining them to voxels whose FA was in the range [0.1, 0.5] in order to remove parts of the cortico-spinal tract (CST) that were incorrectly included in the original segmentation, as proposed by Mang et al. [56].

Tracking.

Tracking was performed separately for the left and right thalami. Seedpoints were placed in the posterior voxels of the thalamus, found automatically as the 3/4 posterior portion between the most anterior point and the most posterior point. Both deterministic and probabilistic tracking were performed with streamtrack from MRtrix [72], with default values, on a FA > 0.1 mask, giving 50n unidirectional tracks where n is the number of voxels in the seeding region, excluding tracks crossing the anterior part of the thalamus. In most of the articles working on the segmentation of the OR, DT deterministic is used (see Table 1). However, the use of HARDI probabilistic tracking has proven useful, allowing reconstructing a larger extent of the OR, namely a longer Meyer’s loop [27]. Tracking results were then filtered, keeping only streamlines touching the whole occipital lobe mask. Streamlines reaching the opposite hemisphere or the anterior part of the brain were discarded with exclusion boxes.

OR and splenium templates.

We observed that the tracking results included part of the splenium of the corpus callosum (SCC). To correct for this, we manually created templates for each bundle—the OR and the SCC—by visual inspection of subject 1’s tractography result. Subject 1 was chosen randomly because his tractogram was standard, i.e average quality not no apparent erroneous results. Here, we kept tracks that broadly resembled either the OR and SCC (see Fig 2, on the left) such that they could be quickly identified in the remaining subjects. The bundles were subsampled using QuickBundles, a clustering algorithm from Dipy, to keep only a few representative streamlines, acting as the templates, as indicated in Fig 3.

thumbnail
Fig 2. Templates of OR (green) and splenium (red).

The difference between deterministic + DTI tracking (paler) and HARDI + probabilistic (darker) was major, mainly in the Meyer’s loop. Probabilistic template was used.

https://doi.org/10.1371/journal.pone.0156436.g002

thumbnail
Fig 3. Templates (on the left) and fiber comparison with MDF (on the right).

Each streamline of the subject (white) is kept if it resembles the OR template (green), rejected otherwise (cases were it resembles the splenium template (red)).

https://doi.org/10.1371/journal.pone.0156436.g003

OR segmentation and metrics.

The OR bundle of each subject was automatically segmented by comparing each streamline of the subject’s tracking result to the OR and SCC templates. Templates were first registered to the diffusion space of that subject by a linear registration. The registration matrix was defined by FLIRT, using a mask combining the two thalami and the occipital lobe. We chose this mask because it was only important to have the posterior part of the brain aligned with the subject’s space, since it contains the OR. The comparison between the streamlines and the templates was performed with the Mean Direct-Flipped (MDF) metric from Dipy. If MDF(streamline, OR) was smaller than MDF(streamline, SCC), the track was kept. Otherwise, it was removed (see Fig 3).

Mean FA and apparent fiber density (AFD) were measured [73, 74], as well as distances between the anterior tip of the temporal horn and the anterior tip of Meyer’s loop in ICBM space (TH-ML distances). TH was selected visually (y = 86) and ML was defined as the most anterior (y) point with at least 5 streamlines. This criteria was added to avoid remaining spurious tracks. Streamline count was computed at each step of the process (initial tracking, after filtering the streamlines reaching the occipital lobe only, after deletion of the SCC). However, to characterize the final OR, streamline count was not used used to quantify the OR since this measure is known to be easily biased [75, 76]. Rather, we computed the mean cross-sectional area (CSA), which is a more robust tract-based metric as it takes into consideration the shape (size, length) of the fiber bundle. This technique is already used in histology [13] and is becoming popular with tractography [77]. We selected 18 equally spaced plans perpendicular to the bundle (sections). We discarded the three first ones since they were too close to the thalamus (see Fig 4). For each section, the CSA was defined as the number of voxels in that plane where passed at least 5 streamlines. The mean CSA was computed for sections 4-15, as well as the minimum CSA, acting as a measure of the stem of the OR. Statistics on the three last sections (16-18), where the OR bundle significantly spreads out, were computed separately to measure the spanning of the OR in the occipital lobe: the maximum and the mean CSA were computed. If a voxel was included in more than one plane, it was counted only once, in the plane closest to the occipital cortex. Voxel count was prefered to other measures such as radius to prevent artificially increased CSAs in spanning areas (ex. planes 13 and 14 in Fig 4).

thumbnail
Fig 4. Superior view of the OR.

Cross section areas (CSAs) were defined as the number of voxels where passed at least 5 streamlines in each section. Mean and minimum CSA of sections 4-15, and mean and maximum CSA of sections 16-18 were measured.

https://doi.org/10.1371/journal.pone.0156436.g004

LGN segmentation.

For each subject, the number of OR streamlines in each voxel was transformed into a density map with MRtrix. Voxels outside the thalamus in the OR density map of each subject were set to zero, and the z-scores of the non-zero voxels were calculated. We defined the LGN of each hemisphere as the voxels with z-score higher than 4. This threshold was chosen to obtain LGN sizes corresponding to the sizes presented in the literature (115mm3 on the left and 131mm3 on the right [13]). The sizes with threshold 4 or with other thresholds were highly correlated, and thus, the threshold value does not influence the final correlations with EEG alpha metrics. LGNs were also registered on ICBM452 standard space [78]. An occurrence map (in percentage of subjects) was obtained in the ICBM space with, for each voxel’s value, the number of subjects whose LGN touched this voxel.

EEG metrics.

EO and EC windows were separated, and the PSD was computed for each window separately. PSD estimates were then averaged separately for the EO and EC conditions. Data were averaged over all P, PO and O electrodes. EC maximum and average power in the alpha band (8-13 Hz) over these electrodes was measured, and same in the spectrum obtained by subtracting the EO spectrum from the EC spectrum, such as was done by Posthuma et al. [79]. The peak frequency, or alpha peak (Pα), was defined as the frequency of the maximum EC power. The same results were measured for the left occipital electrodes and for the right occipital electrodes separately.

fMRI metrics.

The activation map of the eyes open-closed task was obtained by computing the mean across each subject’s activation map and thresholding it for z-scores 2.5 and higher.

Results

Thalami masks

The thalami mask had a mean size of 9450 mm3 ± 1158 in the left hemisphere and 9110 mm3 ± 1044 in the right, representing 81.4% and 85.1% of the mask given by FSL, respectively (see Fig 5).

thumbnail
Fig 5. Thalamus mask modification.

Blue: thalamic masks as defined by FSL. Green: voxels where FA ∈[0.1, 0.5]. Background: the FA map.

https://doi.org/10.1371/journal.pone.0156436.g005

Tractography, templates and OR segmentation

Tracking results showed a large Meyer’s loop and a wide spanning of streamlines in the occipital lobe (see Figs 2 and 6). Number of streamlines touching the occipital area represented, in the left and right hemisphere, 3.3% ± 0.9 and 2.4% ± 0.7 of the initial tracking. The final number of streamlines in the OR (SCC deleted) then represented 70% ± 13 and 48% ± 17 of the occipital-touching ones. Results were similar with deterministic tracking. Cross-sectional areas (CSAs) are reported in Table 2 and represented in Fig 4. Although the fiber count indicated a much bigger left OR than right (2.7 times bigger, p < 0.1e−6), the mean CSA was more similar (left was approximately 1.3 times larger than right, p = 0.02). Also, left and right OR CSAs were correlated (p < 0.0005), while fiber count was not (p = 0.7). The average FA and AFD within the OR are shown in Table 3. TH-ML distances were of −6 ± 3 mm on the left and −4 ± 3 mm on the right (negative results are anterior to the tip of the temporal horn). TH-ML distance with DTI and deterministic tracking were, respectively, of −2 ± 4 mm and 2 ± 5 mm.

LGN segmentation

LGN sizes obtained from the thresholded density map (see Fig 7) were of 116 ± 18mm3 on the left, and 100 ± 26mm3 on the right. The position of the centroid of the LGN (the point with the strongest density, i.e. the most streamlines), varied very little across subjects, as shown in Table 4.

thumbnail
Fig 7. Effect of the thresholding on the ORs.

In the thalamus: in blue, the final OR, voxels crossed by at least 5 fibers and in red, with a threshold on the z-score of 4. In the optic radiation: density map, in number of streamlines. On the left: z = 64. On the right: z = 69.

https://doi.org/10.1371/journal.pone.0156436.g007

thumbnail
Table 4. Position of the strongest density point, in ICBM space.

https://doi.org/10.1371/journal.pone.0156436.t004

EEG alpha waves measurements

Visual inspection confirmed that no subject presented a peak outside the 8-13 Hz range. All of these measures were also similar across hemispheres, as shown in Table 5 (p >0.05 in all cases).

Relation with fMRI results

To assess the quality of our LGN position, we visually compared the occurence maps of the OR in the ICBM space with the maps of the mean fMRI activation of all subjects. Thalamic voxels with highest occurence score were localized near the fMRI activation sites, as shown in Fig 8.

thumbnail
Fig 8. Validation of LGN segmentation.

(left) The red/yellow overlay represents fMRI activity during the EC-EO task. Note the widespread activation in the visual cortex and more focal bilateral activation in the thalamus (LGN). (right) An enhanced view shows that our reconstructed streamlines converge nicely onto the fMRI activation sites (outline in white).

https://doi.org/10.1371/journal.pone.0156436.g008

Link between structure and EEG alpha waves

Overall, we computed 68 Pearson correlations and associated p-values between EEG alpha and dMRI measures in each hemisphere (see Table 6). All scatter plots were visually inspected to ensure that these results were not driven by outliers, and the smallest p-value was 0.06.

thumbnail
Table 6. p-values for all combinations between EEG data and structural data.

Smallest value was 0.06.

https://doi.org/10.1371/journal.pone.0156436.t006

Discussion

The main goal of this study was to investigate the relationship between structure and function in the brain by considering the influence of thalamus and optic radiation (OR) morphology on EEG alpha waves. To do so, we developed a novel approach for quantifying subject-to-subject variations in different measures of thalamic and OR structures with a nearly-automated fiber analysis tool. The selection of the ROIs to perform the tractography is automatic. The sorting of the OR streamlines uses templates that can be defined by a non-expert. The comparison of the streamlines to the templates is automatic and does not need thresholds. Our method allowed us to obtain precise OR to isolate the LGN in multiple subjects, proving its efficiency. Nevertheless, no significant correlation with EEG was observed.

Structure

Optic radiation.

The OR is a well-studied fiber bundle. In humans, it is known to emerge from the LGN, to display a relatively large loop that partially extends into the anterior temporal lobe (Meyer’s loop), and to end as the inferior sagittal stratum in the visual cortex. Although it is made of one single sheet, it is constituted of three parts, each of them having a different curvature [26].

Our reconstructed tracts also follow this pattern (Fig 6). To ensure recovery of the whole OR, we used few anatomical priors, placing a large ROI for the seeding points (3/4 of the thalamus, with no condition on the initial direction), and using the whole occipital cortex as target region for the streamlines. Probabilistic tracking allowed the streamlines to successfully follow the OR, even in the highly curved sections, without additional ROI. The fact that the occipital-touching streamlines represented only around 2-3% of the initial streamlines highlights the difficulty in recovering the OR with tractography. This number would be smaller with whole brain tractography, and an extremely large number of initial seeds would be needed before obtaining a satisfactory OR. For this reason, tracking from the thalamus allows better results, though it may include half of the splenium if the initial thalamus segmentation is overly conservative. Our fiber comparison to a template allowed differentiating the OR bundle from the SCC.

Distances between the anterior tip of Meyer’s loop and the anterior tip of the temporal horn were comparable to other literature reports. Indeed, in Ebeling and Reulen [80], TH-ML distances in cadavers were of -5 mm in average, which is similar to our results. Few in vivo studies have reproduced this [27, 81]. It might be due to the use of probabilistic tracking, as discussed by Lilja et al. [27], but also to the use of HARDI instead of DTI. Indeed, these distances with our DTI and deterministic tracking were smaller (i.e. the Meyer’s loops were less anterior). Also, as reported by others [23, 25, 27], left Meyer’s loop (ML) was more anterior, as revealed with the TH-ML measure, and left OR was bigger, in the streamlines count, but also, to a smaller degree, in the average cross-sectional area (CSA). Our methodology permits different measures, in addition to the cross-sections (e.g. circumference), though this is outside the main scope of our study and remains an active field of research for bundle-based and tract-based metrics. In summary, our OR reconstruction exhibits properties that closely resemble its known anatomy despite the fact that the diffusion images were acquired at relatively low field strength (1.5T). Increasing the spatial resolution is unlikely to affect our results as we purposely focused on more macroscopic measures of tractography.

Density maps and LGN.

Our track density maps appeared to peak near the LGN, which was then confirmed using fMRI (Fig 8). The peak density location was very stable across subjects (see Fig 6), suggesting that our approach can not only accurately localize the LGN, but also that these 3D coordinates may be included as an anatomical reference in a standardized space, such as the ICBM.

However, it should be noted that the nearby pulvinar, whose main projections mainly target visual areas V3 to V5, also show connections the V1 and V2 [82]. It is therefore possible that some of our reconstructed fibers may also arise from the pulvinar and it is perhaps more accurate to label our purported LGN area as “visual thalamus”. We are currently working toward a more accurate segmentation of the LGN and pulvinar.

We should also note that LGN size depends most probably on the initial thalamus size, in itself a research subject [37, 83].

Function

For all subjects, the peak alpha frequency was, as expected [7], situated between 8 and 13 Hz, and all subjects showed a robust change in power during EO-EC in occipital/parietal electrodes, revealing a good control of the signal to noise ratio. For our analysis, we used the maximum alpha power. For the average power, it has been suggested that a subject-driven band instead of the usual 8-13 Hz is more appropriate [7]. However, the use of average power was unnecessary here because it was highly correlated to the maximum and brought no additional information on the EEG variability across subjects.

Structure vs function

No significant correlation was found, for any pair of metrics between structural metrics and EEG metrics (all p-values > 0.05), suggesting that OR and LGN structures have little influence on EEG alpha inter-individual variability. This idea is further supported by the fact that OR segmentation results showed lateralization, but not EEG. The reason behind this dissociation is unclear. On the one hand, our results are in contrast to those reported by Hernandez, who showed a significant correlation between tensor-derived FA of the posterior and superior corona radiata, probably associated to interactions between thalamus and cortex and in the posterior commissural fibers of the Corpus Callosum (splenium), and alpha peak frequency [17]. However, as stated in the introduction, such FA metrics are strongly biased by fibre crossings, which might explain our lack of correlation when using HARDI-based reconstructions. Other studies using voxel based morphology (VBM) have also reported low correlations (r = 0.1—0.3) between WM volume and alpha power [6, 84, 85], which is more in line with our findings, though it should be noted that VBM- and dMRI-derived measures of WM are based on very different phenomena. More work in how these two modalities are related could help in reconciling these discrepancies. For instance, it possible that certain microstructural properties of the LGN/OR/OT, such as myelin content [86], are better correlated to EEG than the WM metrics identified in the current study.

In summary, we here report on a novel approach for identifying the visual thalamus as well as its efferent fibre tracts. This method could be easily extended for accurately segmenting other thalamic nuclei and fiber bundles in different cortical areas, using templates in bundle atlases, for improving our understanding of structure-function relationships in humans.

Acknowledgments

Many thanks to Gabriel Girard and Maxime Chamberland for their precious listening and their valued suggestions. Thanks to the Fonds Québécois de Recherche en Nature et Technologie (FQRNT), to the Faculty of Medecine of the University of Sherbrooke and to the The Natural Sciences and Engineering Research Council of Canada (NSERC) for their financial support.

Author Contributions

Conceived and designed the experiments: KW. Performed the experiments: ER MB. Analyzed the data: ER MB. Contributed reagents/materials/analysis tools: MD EG. Wrote the paper: ER KW MD. Helped with the conception of figures: ER MB.

References

  1. 1. Sporns O. Network attributes for segregation and integration in the human brain. Current Opinion in Neurobiology. 2013;23(2):162–171. pmid:23294553
  2. 2. Sasson E, Doniger GM, Pasternak O, Assaf Y. Structural correlates of memory performance with diffusion tensor imaging. NeuroImage. 2010;50:1231–1242. pmid:20045476
  3. 3. Steele CJ, Scholz J, Douaud G, Johansen-Berg H, Penhune VB. Structural correlates of skilled performance on a motor sequence task. Frontiers in Human Neuroscience. 2012;6. pmid:23125826
  4. 4. Moretti DV, Paternicò D, Binetti G, Zanetti O, Frisoni GB. EEG markers are associated to gray matter changes in thalamus and basal ganglia in subjects with mild cognitive impairment. NeuroImage. 2012 Mar;60(1):489–96. pmid:22166796
  5. 5. Kanai R, Rees G. The structural basis of inter-individual differences in human behaviour and cognition. Nature Reviews Neuroscience. 2011;12(4):231–242. pmid:21407245
  6. 6. Whitford TJ, Rennie CJ, Grieve SM, Clark CR, Gordon E, Williams LM. Brain maturation in adolescence: Concurrent changes in neuroanatomy and neurophysiology. Human Brain Mapping. 2007;28:228–237. pmid:16767769
  7. 7. Bazanova O, Vernon D. Interpreting EEG alpha activity. Neuroscience & Biobehavioral Reviews. 2013 May;.
  8. 8. Lopes da Silva FH, Vos JE, Mooibroek J, van Rotterdam A. Relative contributions of intracortical and thalamo-cortical processes in the generation of alpha rhythms, revealed by partial Coherence analysis. Electroencephalography and clinical neurophysiology. 1980;50:449–456. pmid:6160987
  9. 9. Goldman RI, Stern JM, Engel J Jr, Cohen MS. Simultaneous EEG and fMRI of the alpha rhythm. Neuroreport. 2002;13(18):2487–2492. pmid:12499854
  10. 10. Liu Z, de Zwart Ja, Yao B, van Gelderen P, Kuo LW, Duyn JH. Finding thalamic BOLD correlates to posterior alpha EEG. NeuroImage. 2012 Nov;63(3):1060–9. pmid:22986355
  11. 11. Barnes GR, Li X, Thompson B, Singh KD, Dumoulin SO, Hess RF. Decreased gray matter concentration in the lateral geniculate nuclei in human amblyopes. Investigative ophthalmology & visual science. 2010;51:1432–1438.
  12. 12. Mcketton L, Kelly KR, Schneider Ka. Abnormal lateral geniculate nucleus and optic chiasm in human albinism. The Journal of comparative neurology. 2014 Aug;522(11):2680–2687. pmid:24639208
  13. 13. Andrews TJ, Halpern SD, Purves D. Correlated size variations in human visual cortex, lateral geniculate nucleus, and optic tract. The Journal of neuroscience: the official journal of the Society for Neuroscience. 1997 Apr;17(8):2859–2868.
  14. 14. Musall S, Von Pföstl V, Rauch A, Logothetis NK, Whittingstall K. Effects of neural synchrony on surface EEG. Cerebral Cortex. 2014;24:1045–1053. pmid:23236202
  15. 15. Haueisen J, Tuch D, Ramon C. The influence of brain tissue anisotropy on human EEG and MEG. Neuroimage. 2002;15:159–166. pmid:11771984
  16. 16. Nunez PL. Implications of white matter correlates of EEG standing and traveling waves. NeuroImage. 2011 Aug;57(4):1293–9. pmid:20382232
  17. 17. Valdés-Hernández Pa, Ojeda-González A, Martínez-Montes E, Lage-Castellanos A, Virués-Alba T, Valdés-Urrutia L, et al. White matter architecture rather than cortical surface area correlates with the EEG alpha rhythm. NeuroImage. 2010 Feb;49(3):2328–2339. pmid:19850139
  18. 18. Bürgel U, Schormann T, Schleicher A, Zilles K. Mapping of histologically identified long fiber tracts in human cerebral hemispheres to the MRI volume of a reference brain: position and spatial variability of the optic radiation. NeuroImage. 1999 Nov;10(5):489–499. pmid:10547327
  19. 19. Nunez PL, Reid L, Bickford RG. The relationship of head size to alpha frequency with implications to a brain wave model. Electroencephalography and Clinical Neurophysiology. 1978;44:344–352. pmid:76540
  20. 20. Tuch DS, Reese TG, Wiegell MR, Makris N, Belliveau JW, Wedeen VJ. High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Magnetic Resonance in Medicine. 2002;48(4):577–582. pmid:12353272
  21. 21. Hindriks R, Woolrich M, Luckhoo H, Joensson M, Mohseni H, Kringelbach ML, et al. Role of white-matter pathways in coordinating alpha oscillations in resting visual cortex. NeuroImage. 2014 Oct;106:328–339. pmid:25449741
  22. 22. Calamante F, Smith RE, Tournier JD, Raffelt D, Connelly A. Quantification of voxel-wise total fibre density: Investigating the problems associated with track-count mapping. NeuroImage. 2015;117:284–293. pmid:26037054
  23. 23. Kammen A, Law M, Tjan BS, Toga AW, Shi Y. Automated retinofugal visual pathway reconstruction with multi-shell HARDI and FOD-based analysis. NeuroImage. 2015;125:767–779. pmid:26551261
  24. 24. Martínez-Heras E, Varriano F, Prčkovska V, Laredo C, Andorrà M, Martínez-Lapiscina EH, et al. Improved Framework for Tractography Reconstruction of the Optic Radiation. Plos One. 2015;10(9):e0137064. Available from: http://dx.plos.org/10.1371/journal.pone.0137064. pmid:26376179
  25. 25. Dreessen de Gervai P, Sboto-Frankenstein UN, Bolster RB, Thind S, Gruwel MLH, Smith SD, et al. Tractography of Meyer’s Loop asymmetries. Epilepsy research. 2014 Jul;108(5):872–882. pmid:24725809
  26. 26. Benjamin CFA, Singh JM, Prabhu SP, Warfield SK. Optimization of tractography of the optic radiations. Human Brain Mapping. 2014;35:683–697. pmid:23225566
  27. 27. Lilja Y, Ljungberg M, Starck G, Malmgren K, Rydenhag B, Nilsson DT. Visualizing Meyer’s loop: A comparison of deterministic and probabilistic tractography. Epilepsy research. 2014 Mar;108(3):481–490. pmid:24559840
  28. 28. Winston GP, Mancini L, Stretton J, Ashmore J, Symms MR, Duncan JS, et al. Diffusion tensor imaging tractography of the optic radiation for epilepsy surgical planning: A comparison of two methods. Epilepsy Research. 2011;97:124–132. pmid:21885257
  29. 29. Stieglitz LHL, Lüdemann WO, Giordano M, Raabe A, Fahlbusch R, Samii M. Optic radiation fiber tracking using anteriorly angulated diffusion tensor imaging: a tested algorithm for quick application. Neurosurgery. 2011 May;68(5):1239–51. pmid:21273926
  30. 30. Chen X, Weigel D, Ganslandt O, Buchfelder M, Nimsky C. Prediction of visual field deficits by diffusion tensor imaging in temporal lobe epilepsy surgery. NeuroImage. 2009 Apr;45(2):286–97. pmid:19135156
  31. 31. Yogarajah M, Focke NK, Bonelli S, Cercignani M, Acheson J, Parker GJM, et al. Defining Meyers loop—temporal lobe resections, visual field deficits and diffusion tensor tractography. Brain. 2009;132:1656–1668. pmid:19460796
  32. 32. Sherbondy AJ, Dougherty RF, Napel S, Wandell BA. Identifying the human optic radiation using diffusion imaging and fiber tractography. Journal of vision. 2008;8(10):12. pmid:19146354
  33. 33. Yamamoto A, Miki Y, Urayama S, Fushimi Y, Okada T, Hanakawa T, et al. Diffusion tensor fiber tractography of the optic radiation: analysis with 6-, 12-, 40-, and 81-directional motion-probing gradients, a preliminary study. AJNR American journal of neuroradiology. 2007;28:92–96. pmid:17213432
  34. 34. Okada T, Miki Y, Kikuta K, Mikuni N, Urayama S, Fushimi Y, et al. Diffusion tensor fiber tractography for arteriovenous malformations: quantitative analyses to evaluate the corticospinal tract and optic radiation. AJNR American journal of neuroradiology. 2007;28(6):1107–13. pmid:17569969
  35. 35. Kikuta KI, Takagi Y, Nozaki K, Hanakawa T, Okada T, Miki Y, et al. Early experience with 3-T magnetic resonance tractography in the surgery of cerebral arteriovenous malformations in and around the visual pathway. Neurosurgery. 2006;58(2):331–337. pmid:16462487
  36. 36. Traynor C, Heckemann RA, Hammers A, O’Muircheartaigh J, Crum WR, Barker GJ, et al. Reproducibility of thalamic segmentation based on probabilistic tractography. NeuroImage. 2010;52:69–85. pmid:20398772
  37. 37. Keller SS, Gerdes JS, Mohammadi S, Kellinghaus C, Kugel H, Deppe K, et al. Volume estimation of the thalamus using freesurfer and stereology: consistency between methods. Neuroinformatics. 2012 Oct;10(4):341–50. pmid:22481382
  38. 38. Magnotta Va, Gold S, Andreasen NC, Ehrhardt JC, Yuh WT. Visualization of subthalamic nuclei with cortex attenuated inversion recovery MR imaging. NeuroImage. 2000 Apr;11(4):341–346. pmid:10725190
  39. 39. Fujita N, Tanaka H, Takanashi M, Hirabuki N, Abe K, Yoshimura H, et al. Lateral geniculate nucleus: anatomic and functional identification by use of MR imaging. AJNR American journal of neuroradiology. 2001 Oct;22(9):1719–1726. pmid:11673167
  40. 40. Deoni SCL, Rutt BK, Parrent AG, Peters TM. Segmentation of thalamic nuclei using a modified k-means clustering algorithm and high-resolution quantitative magnetic resonance imaging at 1.5 T. NeuroImage. 2007;34:117–126. pmid:17070073
  41. 41. Calamante F, Oh SH, Tournier JD, Park SY, Son YD, Chung JY, et al. Super-resolution track-density imaging of thalamic substructures: comparison with high-resolution anatomical magnetic resonance imaging at 7.0T. Human brain mapping. 2013 Oct;34(10):2538–2548. pmid:23151892
  42. 42. Traynor CR, Barker GJ, Crum WR, Williams SCR, Richardson MP. Segmentation of the thalamus in MRI based on T1 and T2. NeuroImage. 2011;56:939–950. pmid:21310246
  43. 43. Tourdias T, Saranathan M, Levesque IR, Su J, Rutt BK. Visualization of intra-thalamic nuclei with optimized white-matter-nulled MPRAGE at 7T. NeuroImage. 2014 Jan;84:534–545. pmid:24018302
  44. 44. Wiegell MR, Tuch DS, Larsson HBW, Wedeen VJ. Automatic segmentation of thalamic nuclei from diffusion tensor magnetic resonance imaging. NeuroImage. 2003 Jun;19(2):391–401. pmid:12814588
  45. 45. Behrens TEJ, Johansen-Berg H, Woolrich MW, Smith SM, Wheeler-Kingshott CAM, Boulby Pa, et al. Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nature neuroscience. 2003 Jul;6(7):750–757. pmid:12808459
  46. 46. Ziyan U, Tuch D, Westin CF. Segmentation of thalamic nuclei from DTI using spectral clustering. Medical Image Computing and Computer-Assisted Intervention—MICCAI 2006. 2006 Jan;p. 807–814.
  47. 47. Duan Y, Li X, Xi Y. Thalamus segmentation from diffusion tensor magnetic resonance imaging. International journal of biomedical imaging. 2007 Jan;2007(2):7.
  48. 48. Jonasson L, Hagmann P, Pollo C, Bresson X, Richero Wilson C, Meuli R, et al. A level set method for segmentation of the thalamus and its nuclei in DT-MRI. Signal Processing. 2007;87:309–321.
  49. 49. Unrath A, Klose U, Grodd W, Ludolph AC, Kassubek J. Directional colour encoding of the human thalamus by diffusion tensor imaging. Neuroscience letters. 2008 Apr;434(3):322–327. pmid:18325671
  50. 50. Ziyan U, Westin CF. Joint segmentation of thalamic nuclei from a population of diffusion tensor MR images. Medical Image Computing and Computer-Assisted Intervention—MICCAI 2008. 2008;Springer B:279–286.
  51. 51. Rittner L, Lotufo RA, Campbell J, Pike GB. Segmentation of thalamic nuclei based on tensorial morphological gradient of diffusion tensor fields. Biomedical Imaging: From Nano to Macro, 2010 IEEE International Symposium on. 2010;p. 1173–1176.
  52. 52. Schultz T. Segmenting thalamic nuclei: what can we gain from HARDI? Medical Image Computing and Computer-Assisted Intervention—MICCAI 2011. 2011 Jan;Springer B:141–148.
  53. 53. Johansen-Berg H, Behrens TEJ, Sillery E, Ciccarelli O, Thompson AJ, Smith SM, et al. Functional-anatomical validation and individual variation of diffusion tractography-based segmentation of the human thalamus. Cerebral Cortex. 2005;15(1):31–39. pmid:15238447
  54. 54. Jbabdi S, Woolrich MW, Andersson JLR, Behrens TEJ. A Bayesian framework for global tractography. NeuroImage. 2007;37:116–129. pmid:17543543
  55. 55. Zhang D, Snyder AZ, Shimony JS, Fox MD, Raichle ME. Noninvasive functional and structural connectivity mapping of the human thalamocortical system. Cerebral Cortex. 2010 May;20(5):1187–1194. pmid:19729393
  56. 56. Mang SC, Busza A, Reiterer S, Grodd W, Klose U. Thalamus segmentation based on the local diffusion direction: a group study. Magnetic resonance in medicine: official journal of the Society of Magnetic Resonance in Medicine / Society of Magnetic Resonance in Medicine. 2011 Jan;67(1):118–126.
  57. 57. Elias WJ, Zheng ZA, Domer P, Quigg M, Pouratian N. Validation of connectivity-based thalamic segmentation with direct electrophysiologic recordings from human sensory thalamus. NeuroImage. 2012;59:2025–2034. pmid:22036683
  58. 58. Côté MA, Girard G, Boré A, Garyfallidis E, Houde JC, Descoteaux M. Tractometer: Towards validation of tractography pipelines. Medical Image Analysis. 2013;17:844–857. pmid:23706753
  59. 59. Alvarez I, Schwarzkopf DS, Clark CA. Extrastriate projections in human optic radiation revealed by fMRI-informed tractography. Brain structure & function. 2014 Jun;p. 1–14.
  60. 60. Coupé P, Yger P, Prima S, Hellier P, Kervrann C, Barillot C. An Optimized Blockwise Non Local Means Denoising Filter for 3D Magnetic Resonance Images. IEEE Transactions on Medical Imaging. 2008;27:425–441. pmid:18390341
  61. 61. Garyfallidis E, Brett M. Dipy, a library for the analysis of diffusion MRI data. Frontiers in neuroinformatics. 2014;8:1–14.
  62. 62. Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. FSL. NeuroImage. 2012;62:782–790. pmid:21979382
  63. 63. Jenkinson M, Bannister P, Brady M, Smith S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. NeuroImage. 2002;17:825–841. pmid:12377157
  64. 64. Avants B, Tustison N, Song G. Advanced Normalization Tools (ANTS). Insight Journal. 2009;p. 1–35.
  65. 65. Descoteaux M, Wiest-Daesslé N, Prima S, Barillot C, Deriche R. Impact of Rician adapted non-local means filtering on HARDI. Medical Image Computing and Computer-Assisted Intervention—MICCAI 2008. 2008;Springer B:122–130.
  66. 66. Cox RW. AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Computers and biomedical research, an international journal. 1996;29:162–173. pmid:8812068
  67. 67. Smith SM. Fast robust automated brain extraction. Human Brain Mapping. 2002;17:143–155. pmid:12391568
  68. 68. Bernier M, Chamberland M, Houde JC, Descoteaux M, Whittingstall K. Using fMRI non-local means denoising to uncover activation in sub-cortical structures at 1.5 T for guided HARDI tractography. Frontiers in Human Neuroscience. 2014 Sep;8:1–11.
  69. 69. Delorme A, Makeig S. EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. Journal of Neuroscience Methods. 2004;134:9–21. pmid:15102499
  70. 70. Whittingstall K, Bartels A, Singh V, Kwon S, Logothetis NK. Integration of EEG source imaging and fMRI during continuous viewing of natural movies. Magn Reson Imaging. 2010;2010/06/29. pmid:20579829
  71. 71. Patenaude B, Smith SM, Kennedy DN, Jenkinson M. A Bayesian model of shape and appearance for subcortical brain segmentation. NeuroImage. 2011;56:907–922. pmid:21352927
  72. 72. Tournier JD, Calamante F, Connelly A. MRtrix: Diffusion tractography in crossing fiber regions. International Journal of Imaging Systems and Technology. 2012;22:53–66.
  73. 73. Raffelt D, Tournier JD, Rose S, Ridgway GR, Henderson R, Crozier S, et al. Apparent Fibre Density: A novel measure for the analysis of diffusion-weighted magnetic resonance images. NeuroImage. 2012;59:3976–3994. pmid:22036682
  74. 74. Dell’Acqua F, Simmons A, Williams SCR, Catani M. Can spherical deconvolution provide more information than fiber orientations? Hindrance modulated orientational anisotropy, a true-tract specific index to characterize white matter diffusion. Human Brain Mapping. 2013;34:2464–2483. pmid:22488973
  75. 75. Jones DK, Knösche TR, Turner R. White matter integrity, fiber count, and other fallacies: The do’s and don’ts of diffusion MRI. NeuroImage. 2013;73:239–254. pmid:22846632
  76. 76. Willats L, Raffelt D, Smith RE, Tournier JD, Connelly A, Calamante F. Quantification of track-weighted imaging (TWI): Characterisation of within-subject reproducibility and between-subject variability. NeuroImage. 2014;87:18–31. pmid:24246491
  77. 77. Groeschel S, Tournier JD, Northam GB, Baldeweg T, Wyatt J, Vollmer B, et al. Identification and interpretation of microstructural abnormalities in motor pathways in adolescents born preterm. NeuroImage. 2014 Feb;87:209–219. pmid:24185027
  78. 78. Rex DE, Ma JQ, Toga AW. The LONI Pipeline Processing Environment. NeuroImage. 2003;19:1033–1048. pmid:12880830
  79. 79. Posthuma D, Neale MC, Boomsma DI, de Geus EJ. Are smarter brains running faster? Heritability of alpha peak frequency, IQ, and their interrelation. Behavior genetics. 2001 Nov;31(6):567–579. pmid:11838534
  80. 80. Ebeling U, Reulen HJ. Neurosurgical topography of the optic radiation in the temporal lobe. Acta neurochirurgica. 1988;92:29–36. pmid:3407471
  81. 81. Wu W, Rigolo L, O’Donnell LJ, Norton I, Shriver S, Golby AJ. Visual pathway study using in vivo diffusion tensor imaging tractography to complement classic anatomy. Neurosurgery. 2013;70:1–19.
  82. 82. Leh SE, Chakravarty MM, Ptito A. The connectivity of the human pulvinar: A diffusion tensor imaging tractography study. International Journal of Biomedical Imaging. 2008;2008. pmid:18274667
  83. 83. Amini L, Soltanian-Zadeh H, Lucas C, Gity M. Automatic Segmentation of Thalamus from Brain MRI Integrating Fuzzy Clustering and Dynamic Contours. IEEE Transactions on Biomedical Engineering. 2004;51:800–811. pmid:15132506
  84. 84. Babiloni C, Frisoni G, Steriade M, Bresciani L, Binetti G, Del Percio C, et al. Frontal white matter volume and delta EEG sources negatively correlate in awake subjects with mild cognitive impairment and Alzheimer’s disease. Clinical Neurophysiology. 2006;117:1113–1129. pmid:16564740
  85. 85. Smit DJa, Boomsma DI, Schnack HG, Hulshoff Pol HE, de Geus EJC. Individual differences in EEG spectral power reflect genetic variance in gray and white matter volumes. Twin research and human genetics. 2012 Jun;15(3):384–392. pmid:22856372
  86. 86. Nunez PL, Srinivasan R, Fields RD. EEG functional connectivity, axon delays and white matter disease. Clinical neurophysiology: official journal of the International Federation of Clinical Neurophysiology. 2014 Apr;.