A resting state network in the motor control circuit of the basal ganglia
© Robinson et al; licensee BioMed Central Ltd. 2009
Received: 13 February 2009
Accepted: 23 November 2009
Published: 23 November 2009
In the absence of overt stimuli, the brain shows correlated fluctuations in functionally related brain regions. Approximately ten largely independent resting state networks (RSNs) showing this behaviour have been documented to date. Recent studies have reported the existence of an RSN in the basal ganglia - albeit inconsistently and without the means to interpret its function. Using two large study groups with different resting state conditions and MR protocols, the reproducibility of the network across subjects, behavioural conditions and acquisition parameters is assessed. Independent Component Analysis (ICA), combined with novel analyses of temporal features, is applied to establish the basis of signal fluctuations in the network and its relation to other RSNs. Reference to prior probabilistic diffusion tractography work is used to identify the basal ganglia circuit to which these fluctuations correspond.
An RSN is identified in the basal ganglia and thalamus, comprising the pallidum, putamen, subthalamic nucleus and substantia nigra, with a projection also to the supplementary motor area. Participating nuclei and thalamo-cortical connection probabilities allow this network to be identified as the motor control circuit of the basal ganglia. The network was reproducibly identified across subjects, behavioural conditions (fixation, eyes closed), field strength and echo-planar imaging parameters. It shows a frequency peak at 0.025 ± 0.007 Hz and is most similar in spectral composition to the Default Mode (DM), a network of regions that is more active at rest than during task processing. Frequency features allow the network to be classified as an RSN rather than a physiological artefact. Fluctuations in this RSN are correlated with those in the task-positive fronto-parietal network and anticorrelated with those in the DM, whose hemodynamic response it anticipates.
Although the basal ganglia RSN has not been reported in most ICA-based studies using a similar methodology, we demonstrate that it is reproducible across subjects, common resting state conditions and imaging parameters, and show that it corresponds with the motor control circuit. This characterisation of the basal ganglia network opens a potential means to investigate the motor-related neuropathologies in which the basal ganglia are involved.
A number of studies dating back to 1995 have shown that when subjects are not engaged in processing externally directed tasks or time-varying stimuli - when they are, from a behavioural perspective, at rest - MR images of the brain show correlated, low frequency fluctuations in functionally related areas. This has been interpreted as indicating functional connectivity between regions [1–5]. A number of distinct, largely independent assemblies, or Resting State Networks (RSNs) have been discovered since using semi-exploratory and exploratory analysis methods.
It has recently been shown that RSN fluctuations explain not only inter-trial variation in the BOLD response  but also behaviour  and that some RSNs are disturbed in pathologies such as Alzheimer's disease (e.g. ). This has fuelled efforts to better characterise these networks through behavioural manipulation [9, 10] and by their interdependence on other networks [11, 12], not only to improve experiment design but also to better understand healthy brain function and a range of neurological and psychiatric conditions.
The first RSNs were discovered using functional connectivity analysis, in which correlation is performed between the time course in a seed voxel or region and that in other voxels, in order to reveal areas whose activity is coupled. The discovery that functional connectivity could be observed between ipsilateral and contralateral sensorimotor regions  was rapidly followed by similar observations for visual and auditory areas , the amygdala  and the thalamus and hippocampus . It was later discovered that the group of regions which have come to be known as the Default Mode network, which had been observed to be more active during rest periods than during task processing [15–17], also show fluctuations characteristic of RSNs during rest periods . The development of group Independent Component Analysis methods allowed a fully exploratory approach to identifying RSNs [18–20], and led to the elucidation of other networks in posterior parietal areas, lateralised left and right frontoparietal regions, the anterior temporal lobe, cerebellum and limbic lobe [9, 21–23]. To date, approximately 10 RSNs have been reproducibly identified [9, 23].
There is no a priori model in functional connectivity analysis, but a seed voxel (or ROI) time-course is selected by the experimenter. This process leaves the approach prone to omission unless correlations are computed between a large number of regions (see, e.g., Achard et al. ), and also to weakening by inter-subject variation if seed regions are defined on the basis of template anatomy rather than individual function. Activation results from functional tasks may be used to define seed regions, but this becomes impractical if many networks are to be analysed in the same data. The sensitivity of the analysis is reduced if sub-regions of the same network are separated according to a hypothesis about possible division of function. Alternatively, erroneous conclusions may be drawn about regions functionally connected if seed regions are used that subsume areas which contribute to different networks. In addition, the signal in a seed region comprises many sources, of both neuronal and non-neuronal origin (such as scanner drift). While the inclusion of motion parameters in the general linear model as well as regressors for global, ventricular and white matter signals allows these effects to be reduced [11, 25], it has proved difficult to separate a number of physiological artefacts from RSN-related fluctuations, such as those arising from changes in respiratory rate from the Default Mode . Likewise, the subtraction of global signal to try to mitigate this problem may introduce artificial anticorrelation relationships between component time courses [25, 27, 28], making study of the interrelation between networks problematic. As a genuinely exploratory method, ICA yields potentially interesting components in the context of the other signal fluctuations in the data and is generally capable of separating overlapping signal sources of physiological and neuronal origin . As such is it well suited to an exploratory analysis aimed of RSNs and their interrelations.
Most RSN findings to date relate to the cerebral cortex. There are reports, mostly restricted to the functional connectivity literature, of correlated fluctuations between isolated subcortical structures, including the amygdala  and the thalamus and hippocampus . The involvement of the hippocampus in a (usually subdivided) sensory-auditory RSN  and in the Default Mode has been noted  as has the inclusion of number of thalamic nuclei in the medial visual, auditory and medial frontal RSNs . Jafri et al. also identify elements of the basal ganglia, albeit in a predominantly cortical frontal parietal subcortical network . We recently reported the existence of an RSN in the basal ganglia and thalamus with weak projections to supplementary motor areas . This network overlaps substantially with that identified in a contemporaneous functional connectivity study by Di Martino et al. (corresponding to seeds in the dorsal and caudal putamen) , and a subsequent incidental noting of the same network by Damoiseaux et al. . The network was not present in any other ICA-based study of which we are aware (e.g. [9, 21–23]).
These scant and inconsistent reports leave the consistency and role of this network open to question. A recent dispute indicates that apparently subtle variations in behavioural condition can give rise to the appearance of spontaneous activation . Both the Di Martino et al. study  and our initial report  were based on subjects visually fixating, unlike the majority of RSN studies, in which subjects had their eyes closed, and in which the basal ganglia RSN was not observed. In that light, it seems pertinent to investigate the conditions under which this network is manifest. The origin of these signal fluctuations likewise needs to be established; whether they are neuronal resting-state fluctuations or a physiological artefact. If the origins of these network are neuronal, to which of the many parallel basal ganglia circuits and functions do they correspond, and what relationship do they have to other resting state networks?
The basal ganglia consist of four nuclei; the striatum (which is subdivided into the caudate nucleus and putamen), the globus pallidus or pallidum, the substantia nigra and the subthalamic nucleus . Originally viewed as motor structures, tracing studies in the monkey suggested that the striatum could be divided into two networks, sensorimotor and associative . A more recent model has suggested that the diverse functions in which the basal ganglia are involved are expressed through five parallel segregated circuits; motor, oculomotor, dorsolateral prefrontal, lateral orbitofrontal and anterior cingulate . Each of these receives input from a number of functionally related neocortical regions (e.g. in the case of the motor circuit, the supplementary motor area, arcuate premotor area, motor cortex and somatosensory cortex) and outputs to a single frontal region (e.g. in the motor circuit, the supplementary motor area). In the current view, the basal ganglia are envisaged as comprising just three distinct functional assemblies; the sensorimotor, the associative and the limbic . The fact that each of the corresponding circuits involves distinct regions of striatum, pallidum, substantia nigra, thalamus and cortex offers the possibility - within the limits of activation localisation - of identifying the circuit to which the network corresponds, and with that, its function.
We examine resting-state fluctuations in the thalamus and basal ganglia using two common resting state conditions, two large, independent study groups and a fast EPI protocol optimised for structures with a short T2* [38, 39], combined with high field strength. We apply group ICA tools to identify this network in studies using fixation and eyes closed conditions, and refer to known basal ganglia circuits and probabilistic tractography to identify the function subsumed by this network. Temporal features of independent components are used in a novel classifier to distinguish RSNs from physiological artefacts, and functional network connectivity is applied to elucidate the relationship between the basal ganglia resting state network and the Default Mode.
We find a resting state network involving the thalamus and a large portion of the basal ganglia in groups studied under both the eyes open and fixation conditions, and a temporally coherent network focussed on the caudate in the fixation study only.
A public-access atlas of thalamo-cortical connection probabilities http://www.fmrib.ox.ac.uk/connect/ based on results of a prior study by Behrens et al.  shows that, for the basal ganglia the thalamic activation foci identified connect with high probability to motor areas. For the focus with MNI co-ordinates (13, -20, 0) these probabilities take the values of 0.68 to the primary motor, 0.44 to the sensory and 0.28 to the pre-motor cortex. For the focus at (-13, -20, 0) the respective probabilities are 0.44, 0.38 and 0.43. The probabilities of the primary thalamic connections being to none-motor regions (occipital, prefrontal, posterior parietal and temporal cortices) are below 0.2 for both foci.
The RSNs identified in the MELODIC analysis of the main study group, the names adopted in this work, the relative strength of the components and correspondence to RSNs identified in other studies (the bracketed suffixes L and R indicate that the component identified here corresponds to either the left or right hand side of the network referred to).
Percentage of explained variance
Labeling in [(Calhoun et al., 2008), (Damoiseaux et al., 2006), (Beckmann et al., 2005), (De Luca et al., 2006)] respectively.
Left lateral fronto-parietal
Right lateral fronto-parietal
Anterior temporal lobe
The same resting state networks, including the basal ganglia resting state network, were identified in the GIFT results, which were used to extract single-subject temporal responses for the assessments of temporal and frequency features that follow. The spatial map of the basal ganglia component identified in the GIFT analysis is shown in Additional File 2.
Mean correlation between the basal ganglia component time course and the time courses of other networks, across subjects, listing the most significant mean correlations.
Basal ganglia, Default Mode
8.4 × 10-6
Basal ganglia, right lateral fronto-parietal
5.8 × 10-4
Basal ganglia, left lateral fronto-parietal
9.4 × 10-7
*Default Mode, right lateral fronto-parietal
*Default Mode, left lateral fronto-parietal
1.6 × 10-4
The "Dynamic Range" feature afforded 93% accuracy in distinguishing resting state networks and physiological components, with 0 false negatives (no RSNs incorrectly identified as a physiological signal source) and 1 false positive (1 physiological component labelled incorrectly as an RSN). The basal ganglia component was classified as an RSN.
We detail the structures contributing to a recently reported resting state network in the thalamus and basal ganglia. By using a high field, high BOLD sensitivity experiment design and high resolution analysis, we show that it encompasses not only the thalamus, pallidum, putamen and transverse temporal gyrus - as has been previously noted - but also the substantia nigra and subthalamic nucleus, allowing the basal ganglia circuit to which it corresponds to be identified. Despite its non-observation in most resting-state studies to date [9, 21–23] it was found to be reproducible across subjects and MR measurement parameters and is evident in both the eyes closed and fixation resting state conditions. The network is positively correlated with the left and right lateral fronto-parietal (also called attention or task-positive) networks, and anticorrelated with the Default Mode, whose hemodynamic response it anticipates. The peak frequency and spectral characteristics are similar to those of other RSNs, distinct from physiological fluctuations, and allow it to be classified as an RSN using the unsupervised classifier described. We proceed to identify which of the parallel segregated basal ganglia circuits this network corresponds to based on participating regions.
Because of their historical significance as motor structures, the best studied basal ganglia network is the motor circuit. Cortical input from precentral motor areas, postcentral somatosensory areas, the arcuate premotor area and the supplementary motor areas projects dominantly to the putamen. The putamen sends projections to the interior segment of the pallidum and on to particular thalamic nuclei (the direct pathway) as well as the internal segment of the globus pallidus via the caudolateral substantia nigra to the thalamus (the indirect pathway). Outputs from the thalamus project to the supplementary motor area, the motor cortex and arcuate premotor area, probably in distinct subcircuits . The observed resting state network corresponds well with the motor circuit; the caudate is notably absent, and the putamen constitutes the focus of activation. The activated regions identified as substantia nigra are consistent with coordinates derived from stereotactic electrophysiological studies . The presence of the substantia nigra suggests that the resting state network corresponds either solely to the indirect pathway of the motor circuit, which acts to inhibit movement, or both the indirect and direct pathways. Weak activation in the supplementary motor area is likewise consistent with the motor circuit. Myeloarchitectonics suggest that the thalamic activation corresponds to the ventro-lateral thalamic nucleus, which a prior DTI study has shown to connect dominantly to motor regions . In fact, endogenous BOLD fluctuations in this part of the thalamus have been found to be strongly correlated with motor and premotor areas by Zhang et al. ( - supplementary material) as well as the whole putamen.
As well as there being good agreement between the areas observed in this RSN and the motor circuit, there is disparity between the principle sites of input to the striatum and the other circuits. The caudate provides input to the oculomotor, dorsolateral prefrontal and lateral orbitofrontal circuits and the ventral striatum provides input to the anterior cingulate circuit. The cortical regions to which the circuits send output are the frontal eye fields for the oculomotor circuit and the dorsolateral prefrontal cortex, the lateral orbitofrontal cortex and the anterior cingulate area for those circuits respectively.
No components were identified corresponding to the associative and limbic thalamo-cortical loops . The question arises of whether it is to be expected that an RSN exists for each of the tripartite or Alexander subdivisions of the basal ganglia, and whether basal ganglia RSNs (if present) would be expected to correspond exactly with circuits that have been established from anatomical and afferent projections from the cortex. Looking outside the basal ganglia, we know that whilst RSNs reflect functional task networks, it is not the case that there is an RSN for every functional network defined by task or anatomy. Where there are spontaneous fluctuations in networks at rest, other studies have shown that there is not complete correspondence with the network as defined by anatomy or tasks. The major anatomical afferents from the cortex suggest a functional partitioning of the putamen, with a minor portion classified as associative and limbic, and the majority as being sensorimotor. Zhang et al. found functional connectivity between the motor and premotor cortex and almost the entire putamen , however, as observed here. The temporal cortex also displayed a much weaker correlation with the caudal-ventral putaminal region, a finding not consistent with the classical striate nucleus tripartite subdivision. In addition, when looking at RSNs involving associative cortical area, such as the Default Mode, ventral attention, dorsal attention and executive control networks, only a portion of the caudate nucleus participates, instead of the entire associative subdivision of the striate nucleus.
Our hybrid simulations suggest that ICA is capable of separating circuits which overlap to some extent, but which also have non-overlapping elements and different temporal behaviour (Additional File 3). Further evidence of this is provided by the fact that some brain regions (such as the thalamus) are present in multiple networks. The most likely explanations for not finding RSNs corresponding to the associative and limbic thalamo-cortical circuits is, therefore, that they either do not show spontaneous low-frequency BOLD fluctuations or that these fluctuations are below the sensitivity of this study.
Di Martino et al. have reported a resting state functional connectivity analysis focussing on the striatum  independent of our preliminary reporting of these results . A gender-mixed group of subjects were studied while fixating on the word "Relax" in a study with a sensitivity likely to be equivalent to that here; applying a similar, relatively short TE EPI protocol and comparable number of total image volumes (6895 c.f. 7800 here). The connectivity results obtained for seeds in the dorsal and caudal putamen in the Di Martino et al. study is similar to the RSN observed by us , although the involvement of the substantia nigra (which contributes substantially to the motor circuit attribution) was not reported. Also, the dorsal putamen seeds indicated correlated activity in the anterior cortex cingulate, which relate to executive function, indicating some mixing of fluctuations relating to motor and executive control circuits. The basal ganglia RSN reported here was recently noted as an incidental finding by Damoiseaux et al. , adding independent verification of these results.
The conditions under which this resting state data were acquired were similar to those used in a number of previous studies which also applied group ICA approaches [9, 21–23]. The question then arises as to why this network was not reported in those studies, but has been observed only here and in the most recent publication by Damoiseaux et al. . The basal ganglia resting state network was one of the weakest identified in this study, measured in terms of the percentage of total variance in the data it explains. The sensitivity of this study is likely to be higher than that in the studies cited due to the high field strength (3 T), relatively large number of subjects (26) and short repetition time (1 s). The echo time of 28 ms is quite short and well matched to the T2* of basal ganglia structures at 3 T , yielding optimal BOLD sensitivity. T2* is shorter in the basal ganglia than is typical in the cortex due high iron concentration . While it is likely that, other than reference , previous studies were not sufficiently sensitive to detect this network, it is also possible that it was simply overlooked in the wealth of components.
Previous studies using functional connectivity analyses have demonstrated an anticorrelation relationship between the Default Mode and the dorsal attention network in the resting state [5, 11, 42, 43], which has been interpreted as indicating an interplay between modes thought to reflect stimulus-independent thought and goal-driven activity. These anticorrelation findings have been recently called into question, however, as global signal subtraction performed as a pre-processing step to reduce the influence of physiological noise and scanner drift itself introduces anticorrelations [27, 28]. The approach taken here is not subject to these problems. We calculated the correlation between independent component time courses, with no prior global signal subtraction. The most significant results were as follow. The anticorrelation finding reported previously between the Default Mode and the dorsal attention network [5, 11, 42, 43] was reproduced. The basal ganglia network was also found to be anticorrelated with the Default Mode, and correlated, to a similar degree, with the lateralised attention networks. Although we and others have shown that RSNs possess similar frequency characteristics, it must be the case that they are at most weakly coupled, or they would not be separable in functional connectivity analyses . Correlation values observed here between component time courses are correspondingly low - in the range of 0.19 to 0.27 in magnitude - consistent with those observed in other studies . Despite the fact that they are weak, their consistency across subjects is such that these results are highly significant. There have been suggestions that the thalamus, with involvement of the basal ganglia, may be responsible for instigating the task-independent deactivation of the Default Mode observed when subjects are posed cognitively demanding tasks . This would be consistent with findings by Uddin et al., which have shown that correlations between homologous RSN structures in the cerebral hemispheres are preserved in a patient with complete commissurotomy, indicating that functional connectivity can be mediated by subcortical structures . The hypothesis that the basal ganglia RSN represents the controller of fluctuations in the Default mode is supported by the results of the functional network connectivity analysis. In 18 out the 20 randomly composed groups of 13 subjects, the time course of the basal ganglia independent component preceded that of the Default Mode. Although this demonstrates that the hemodynamic response in the basal ganglia RSN consistently precedes that in the Default Mode, latency differences in the hemodynamic response functions  in both networks would have to be analysed and corrected for  before concluding that activation in the basal ganglia network precedes deactivation in the Default Mode. Possible differences between the precedence of neuronal activation and the measured MR response arising from hemodynamic shape and latency effects  preclude testing order hypotheses with this and other approaches such as Granger Causality Modeling.
Although RSNs show a maximum in frequency spectra in the range 0.01 - 0.04 Hz, this does not accurately reflect the frequency distribution of underlying neuronal fluctuations. The intrinsic autocorrelation of BOLD fMRI data has 1/f behaviour in the frequency domain , and low frequencies are cut off by sampling over a finite duration, leading to a low frequency peak. In fact, when the hemodynamic response function is deconvolved from RSN time courses prior to frequency analysis, the spectra of RSN components are essentially flat up to 0.1 Hz . The analysis of RSN frequencies in this work serves two purposes, neither of which relate to the absolute frequencies observed in the spectra. The first is that differences between the compositions of RSN spectra can be observed. The frequency spectrum of the basal ganglia network is most similar in composition to that of the Default Mode, reinforcing the connection between the two networks. The second is that RSN spectra may be distinguished from components of physiological origin because their frequency distributions reflect the intrinsic convolution with the hemodynamic response function. This is the basis for the classifier used here, which can reliably distinguish RSNs from physiological components using frequency characteristics alone . The basal ganglia RSN clusters clearly with other RSNs.
It is likely that the component which consists solely of the caudate (Figure 3) reflects activity relating to the inhibition of ocular saccades. As such, caudate activity is a task-specific response rather than an RSN. The caudate component is apparent in the fixation study, in which there was a point fixation condition, while no corresponding activity is apparent in the eyes closed study. The role of the caudate in saccadic eye movements is well documented [55–57]. Fluctuations in caudate activity in this capacity (which are a pre-requisite of their identification in an ICA) may arise due to phasic preparation of reflexive saccades and their voluntary inhibition. Ultimately controlled by the superior colliculus, saccades can be generated and inhibited by the input of the caudate nucleus, via efferent projections to the substantia nigra pars reticulata . The oculomotor basal ganglia-thalamocortical circuit, includes, as origins of input to the caudate, frontal eye fields, dorsolateral prefrontal cortex and posterior parietal cortex (Brodmann's areas 8, 9&10 and 7, respectively)  all of which areas are apparent in this component. Frontal eye field neurons are known to fire during passive fixation  as was the condition in the fixation study. The oculomotor network yielded from the dorsal caudate seed in the Di Martino et al. analysis  is similar to the caudate component identified here in the fixation study. Our hypothesis that this is a task-related response to the suppression of ocular saccades is consistent with the use of a fixation condition in that work, and it not being observed in the many previous studies which have used the eyes closed condition [9, 21–23, 32].
Another area of application of this network is as a candidate marker for neuro- and psychopathologies involving the basal ganglia. In a parallel with attempts to use Default Mode activity as a diagnostic marker for Alzheimer's disease , deficits in the basal ganglia RSN may offer a marker for one or more of the diseases in which the basal ganglia are known to play a role. Parkinson's disease is a basal ganglia disorder characterized by the degeneration of dopaminergic neurons in the striatum. This has been shown to have opposing effects on the direct and indirect pathways in the striatum , with hyperkinesia as a result. Functional connectivity in the basal ganglia-cortical circuit has been demonstrated in this condition via synchronous oscillations between local field potentials in the basal ganglia and cortical EEG , the transmission to the basal ganglia of motor cortex electrostimulation in the monkey and rat  and oscillatory high-voltage spindles). Dopaminergic lesion in a rodent model for Parkinson's has been shown to lead to an increase in oscillatory synchronisation in the basal ganglia and increase in frequency and duration of high voltage spindle events). Similarly, excessive synchronisation of subthalamic nucleus neurons has been confirmed as a cause of movement slowing in Parkinsonism . Dopamine-dependent changes in the functional connectivity between the basal ganglia and cortex have likewise been demonstrated . The discovery of the fMRI manifestation of this functional connectivity will allow dysfunction in this system to be probed non-invasively, even if patients are not able to perform motor tasks appropriately.
We report the existence of a resting state in the thalamus and basal ganglia, showing greatest activity in the putamen, pallidum, substantia nigra and subthalamic nucleus, with projections to transverse temporal gyrus and the supplementary motor areas. The network is consistent with the cortico-subcortical motor control circuit of the basal ganglia and is robustly reproducible across subjects, scanning parameters and common behavioural conditions for resting state studies. The frequency spectrum of this component peaks at 0.025 ± 0.007 Hz and is most similar in frequency composition to the Default Mode. Spectral features are similar to those of other resting state networks and distinct from physiological artefacts. Fluctuations in the basal ganglia network precede those in the Default Mode, with which it is anticorrelated. We posit that another independent component focussed in the caudate, which was observed in the fixation condition but not in the eyes closed condition, and which was recently reported as being associated with the resting state of the brain , is related to the suppression of ocular saccades.
The basal ganglia resting state network reported here offers the possibility to improve experiment design and analysis in fMRI studies of the striatum, and a possible window into disorders of the basal ganglia such as Parkinson's disease.
MR data acquisition, fixation study (3 T)
Twenty-six female right-handed subjects with no history of neurological or psychiatric disease and aged between 22 and 41 years (mean 26 ± 5 years) participated in the study, which was approved by the Ethics Committee of the Medical University of Vienna, with informed written consent. Magnetic resonance images were acquired with a 3 T Bruker Medspec S300 scanner (Bruker Biospin, Ettlingen, Germany) using a birdcage head coil. T1-weighted structural images were obtained with a 3D MPRAGE sequence with TE = 8 ms; flip angle = 15°; TA = 13 min. For the resting state run, which was of 5 minutes duration, subjects were asked to fixate on a point projected onto a screen mounted in the scanner bore, which they viewed via a mirror attached to the head resonator. They were instructed not to engage in organised thought and not to sleep. Compliance was established verbally on completion of the experiment. Resting state EPI data was acquired with 18 oblique slices (parallel to a line defined by the anterior and posterior commissures (ACPC)) of 4 mm thickness with a 1 mm nominal gap. Other imaging parameters were as follows: TE/TR = 28/1000 ms, a matrix size of 64 × 64 and a field of view of 21 × 25 cm2 (left-right and anterior-posterior, respectively) yielding 3.3 × 3.9 × 4 mm3 voxels, NR = 300, TA = 5 min. Imaging was prefaced by 10 s of dummy scans to ensure a steady state of longitudinal magnetisation.
MR data acquisition, eyes closed study (4 T)
A second group of subjects was studied in the eyes closed condition. Fifteen subjects (8 males) with no history of neurological or psychiatric disease and aged between 19 and 56 years, mean 36 ± 12 years, participated in the study, which was approved by the Ethics Committee of the University of Trento, with informed written consent. Magnetic resonance images were acquired with a 4 T Bruker Medspec scanner (Bruker Medical, Ettlingen, Germany) using a birdcage-transmit, eight-channel receive head coil (USA Instruments, Inc., Ohio, USA). T1-weighted structural images were obtained using a 3D MPRAGE sequence with TE = 4 ms; flip angle = 7°; iPAT factor 2, TA = 5 min optimized for maximal contrast to noise ratio between grey and white matter at 4 T . Subjects were asked to close their eyes during two resting state runs of 10 minutes duration each, and not engage in organised thought and not sleep (compliance established verbally as in the fixation study). EPI data were acquired along the ACPC line, with 37 oblique slices of 3 mm thickness with a 0.5 mm nominal gap, TE/TR = 33/2200 ms, a matrix size of 64 × 64 and 192 × 192 mm, yielding 3 × 3 × 3 mm voxels. 273 volumes plus 2 preparation scans were acquired in a 10 min session. Resting state runs were prefaced by a point-spread function acquisition to allow the correction of geometric distortions , which has been demonstrated to be effective at high field .
In the eyes closed study, distortion correction of EPI data was performed online using the point-spread function method as implemented in Siemens Distortion Correction WIP Version 2.5 . Functional data acquired in both studies were preprocessed using SPM5 (motion correction, normalisation to the EPI template, no spatial or temporal smoothing). For both the fixation and control studies, Probabilistic Independent Component Analysis (PICA) was performed with MELODIC Version 3.05, part of FSL 4.0 . The multi-session temporal concatenation group ICA approach was employed in order to find common spatial patterns amongst subjects without assuming similarity in temporal responses. Whilst the same approach is employed in the Group ICA of FMRI Toolbox (GIFT) , we used MELODIC to estimate component spatial maps in order both to facilitate comparison with the majority of studies which have looked for RSNs common to healthy subjects [21–23, 32] and because probabilistic approach allows alternative hypothesis testing of the significance of activated voxels . Downsampling of data is common in group ICA to reduce total data volume and to increase SNR. Here we analysed data at two resolutions - downsampled to 4 mm and also with no downsampling, the latter to try to resolve activation in small structures. In the PICA analysis masking was applied to exclude non-brain voxels. The data were de-meaned on a voxel-by-voxel basis, and the voxel-wise variance normalised. The number of components into which the data were decomposed was determined using the Laplace approximation to the Bayesian evidence for a probabilistic principal component analysis model implemented in MELODIC . Estimated component maps were divided by the standard deviation of the residual noise. Statistical significance was attributed by fitting a mixture model to the histogram of intensity values . To establish the reproducibility of responses the same analysis was applied to two equally sized, randomly selected, non-overlapping subgroups of 13 subjects in the fixation study; Subgroup 1, (mean age 25 ± 3 years) and Subgroup 2 (mean age 27 ± 6 years) (age difference not significant in student's t-test; P = 0.37).
Anatomical connectivity determines function . The likelihood of brain regions being connected may be determined from diffusion weighted imaging data in an approach known as probabilistic diffusion tractography , which can provide useful insight into the roles of RSNs . The probabilities of activation foci in the thalamus connecting to particular cortical regions were assessed using an open-access probabilistic diffusion tractography atlas http://www.fmrib.ox.ac.uk/connect/ based on the results of a study by Behrens et al. .
Data in the fixation study were additionally analysed with GIFT . GIFT employs a temporal concatenation approach like that used in this MELODIC analysis, but also back-reconstructs single-subject spatial maps in addition to time courses. In our experience, the temporal responses of single subjects are better separated by GIFT, and we used these as the basis for analysis of temporal features, including characterisation of the frequency distributions of RSNs, the temporal relationships between RSNs (correlations between component time courses) and the attempt to distinguish RSNs and physiological components using temporal features. These further analyses were performed using MATLAB (Mathworks Inc, Natick, MA) routines developed in-house.
Frequency spectra were calculated from component time courses using Welch's averaged, modified periodogram spectral estimation method, using a Hamming window over periods of 64 s, with 50% overlap between segments. The peak frequency was calculated as the mean frequency over subjects at which the spectral power was a maximum.
To investigate correlations between component time courses (which may indicate functional relationships between networks [11, 42]) the correlation coefficients between all possible pairs of component time courses were calculated. To avoid the possibility of introducing anticorrelations between networks [25, 27, 28], no global signal subtraction was performed. Fischer transformation was applied to single-subject correlation values (to normalise their distribution), and the null hypothesis of no correlation between the time courses tested with one-sample t tests at the P < 0.05 level.
An extension of the approach of correlating RSN time courses, termed Functional Network Connectivity , has been developed to allow the study of not only the correlation between networks, but also which leads and which lags in response. The functional network connectivity toolbox http://mialab.mrn.org/ was applied to GIFT results. In keeping with the identification of the peak in spectral power at 0.025 ± 0.007 Hz, correlations were assessed in the frequency range 0.01 to 0.4 Hz. Only lags between components reflecting the Default Mode and the basal ganglia resting state network were assessed, in view of the hypothesised connection between the two networks. A significance threshold of P < 0.05 was set for correlations. As an additional means to assess the reliability of calculated correlations and lags, the analysis was performed for each of 20 randomly-composed groups of 13 subjects from the population of 26 subjects.
where ms,k-1, mi,k-1 are the means of the distributions inferior and superior to the threshold for the previous iteration. This was performed until θk = θk-1. To avoid introducing bias, the basal ganglia component was not included in the development of the classifier or the establishment of discrimination thresholds.
The most salient details of the two study groups and the analysis methods applied to each are shown in a data analysis flow chart, Additional File 4.
This study was supported by the government of the Provincia Autonoma di Trento, Italy, the private foundation Fondazione Cassa di Risparmio di Trento e Rovereto, the University of Trento, Italy, grant 11437 from the Austrian National Bank and by a Siemens Medical Solutions research grant to EM.
- Biswal B, Yetkin FZ, Haughton VM, Hyde JS: Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995, 34 (4): 537-541. 10.1002/mrm.1910340409.View ArticlePubMed
- Biswal BB, Van Kylen J, Hyde JS: Simultaneous assessment of flow and BOLD signals in resting-state functional connectivity maps. NMR Biomed. 1997, 10 (4-5): 165-170. 10.1002/(SICI)1099-1492(199706/08)10:4/5<165::AID-NBM454>3.0.CO;2-7.View ArticlePubMed
- Cordes D, Haughton VM, Arfanakis K, Carew JD, Turski PA, Moritz CH, Quigley MA, Meyerand ME: Frequencies contributing to functional connectivity in the cerebral cortex in "resting-state" data. AJNR Am J Neuroradiol. 2001, 22 (7): 1326-1333.PubMed
- Lowe MJ, Mock BJ, Sorenson JA: Functional connectivity in single and multislice echoplanar imaging using resting-state fluctuations. NeuroImage. 1998, 7 (2): 119-132. 10.1006/nimg.1997.0315.View ArticlePubMed
- Greicius M, Krasnow B, Reiss A, Menon V: Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci USA. 2003, 100 (1): 253-258. 10.1073/pnas.0135058100.PubMed CentralView ArticlePubMed
- Fox MD, Snyder AZ, Zacks JM, Raichle ME: Coherent spontaneous activity accounts for trial-to-trial variability in human evoked brain responses. Nat Neurosci. 2006, 9 (1): 23-25. 10.1038/nn1616.View ArticlePubMed
- Fox MD, Snyder AZ, Vincent JL, Raichle ME: Intrinsic fluctuations within cortical systems account for intertrial variability in human behavior. Neuron. 2007, 56 (1): 171-184. 10.1016/j.neuron.2007.08.023.View ArticlePubMed
- Greicius MD, Srivastava G, Reiss AL, Menon V: Default-mode network activity distinguishes Alzheimer's disease from healthy aging: evidence from functional MRI. Proc Natl Acad Sci USA. 2004, 101 (13): 4637-4642. 10.1073/pnas.0308627101.PubMed CentralView ArticlePubMed
- Calhoun VD, Kiehl KA, Pearlson GD: Modulation of temporally coherent brain networks estimated using ICA at rest and during cognitive tasks. Hum Brain Mapp. 2008, 29 (7): 828-838. 10.1002/hbm.20581.PubMed CentralView ArticlePubMed
- Morcom AM, Fletcher PC: Does the brain have a baseline? Why we should be resisting a rest. Neuroimage. 2007, 37 (4): 1073-1082. 10.1016/j.neuroimage.2006.09.013.View Article
- Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME: The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci USA. 2005, 102 (27): 9673-9678. 10.1073/pnas.0504136102.PubMed CentralView ArticlePubMed
- Jafri MJ, Pearlson GD, Stevens M, Calhoun VD: A method for functional network connectivity among spatially independent resting-state components in schizophrenia. NeuroImage. 2008, 39 (4): 1666-1681. 10.1016/j.neuroimage.2007.11.001.PubMed CentralView ArticlePubMed
- Biswal BB, Yetkin FZ, Haughton VM, Hyde JS: Functional connectivity in the auditory cortex studied with FMRI. NeuroImage. 1996, 3: S305-10.1016/S1053-8119(96)80307-7.View Article
- Stein T, Moritz C, Quigley M, Cordes D, Haughton V, Meyerand E: Functional connectivity in the thalamus and hippocampus studied with functional MR imaging. AJNR Am J Neuroradiol. 2000, 21 (8): 1397-1401.PubMed
- Shulman G, Fiez J, Corbetta M, Buckner R, Miezin F, Raichle M, Petersen S: Common blood flow changes across visual tasks: II. Decreases in cerebral cortex. Journal of Cognitive Neuroscience. 1997, 9 (3): 648-663. 10.1162/jocn.19220.127.116.118.View ArticlePubMed
- Raichle M, MacLeod A, Snyder A, Powers W, Gusnard D, Shulman G: A default mode of brain function. Proc Natl Acad Sci USA. 2001, 98 (2): 676-682. 10.1073/pnas.98.2.676.PubMed CentralView ArticlePubMed
- Mazoyer B, Zago L, Mellet E, Bricogne S, Etard O, Houde O, Crivello F, Joliot M, Petit L, Tzourio-Mazoyer N: Cortical networks for working memory and executive functions sustain the conscious resting state in man. Brain Res Bull. 2001, 54 (3): 287-298. 10.1016/S0361-9230(00)00437-8.View ArticlePubMed
- Beckmann CF, Smith SM: Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans Med Imaging. 2004, 23 (2): 137-152. 10.1109/TMI.2003.822821.View ArticlePubMed
- Beckmann CF, Smith SM: Tensorial extensions of independent component analysis for multisubject FMRI analysis. NeuroImage. 2005, 25 (1): 294-311. 10.1016/j.neuroimage.2004.10.043.View ArticlePubMed
- Calhoun VD, Adali T, Pearlson GD, Pekar JJ: A method for making group inferences from functional MRI data using independent component analysis. Hum Brain Mapp. 2001, 14 (3): 140-151. 10.1002/hbm.1048.View ArticlePubMed
- Beckmann CF, DeLuca M, Devlin JT, Smith SM: Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond B Biol Sci. 2005, 360 (1457): 1001-1013. 10.1098/rstb.2005.1634.PubMed CentralView ArticlePubMed
- De Luca M, Beckmann CF, De Stefano N, Matthews PM, Smith SM: fMRI resting state networks define distinct modes of long-distance interactions in the human brain. NeuroImage. 2006, 29 (4): 1359-1367. 10.1016/j.neuroimage.2005.08.035.View ArticlePubMed
- Damoiseaux JS, Rombouts SA, Barkhof F, Scheltens P, Stam CJ, Smith SM, Beckmann CF: Consistent resting-state networks across healthy subjects. Proc Natl Acad Sci USA. 2006, 103 (37): 13848-13853. 10.1073/pnas.0601417103.PubMed CentralView ArticlePubMed
- Achard S, Salvador R, Whitcher B, Suckling J, Bullmore E: A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J Neurosci. 2006, 26 (1): 63-72. 10.1523/JNEUROSCI.3874-05.2006.View ArticlePubMed
- Weissenbacher A, Kasess C, Gerstl F, Lanzenberger R, Moser E, Windischberger C: Correlations and anticorrelations in resting-state functional connectivity MRI: a quantitative comparison of preprocessing strategies. NeuroImage. 2009, 47 (4): 1408-1416. 10.1016/j.neuroimage.2009.05.005.View ArticlePubMed
- Birn RM, Diamond JB, Smith MA, Bandettini PA: Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. NeuroImage. 2006, 31 (4): 1536-1548. 10.1016/j.neuroimage.2006.02.048.View ArticlePubMed
- Smith S, Niazy R, Beckmann C, Miller K: Resting State Networks - Neither Low Frequency Nor Anticorrelated?. NeuroImage. 2008, 41 (Supplement 1): S181-227.
- Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA: The impact of global signal regression on resting state correlations: are anti-correlated networks introduced?. NeuroImage. 2009, 44 (3): 893-905. 10.1016/j.neuroimage.2008.09.036.PubMed CentralView ArticlePubMed
- Birn RM, Murphy K, Bandettini PA: The effect of respiration variations on independent component analysis results of resting state functional connectivity. Hum Brain Mapp. 2008, 29 (7): 740-750. 10.1002/hbm.20577.PubMed CentralView ArticlePubMed
- Robinson S, Soldati N, Basso G, Sailer U, Jovicich J, Bruzzone Ll, Kryspin-Exner I, Bauer H, Moser E: A Resting State Network in the Basal Ganglia. Proceedings of the Sixteenth Annual Meeting of the ISMRM, Toronto. 2008
- Di Martino A, Scheres A, Margulies DS, Kelly AM, Uddin LQ, Shehzad Z, Biswal B, Walters JR, Castellanos FX, Milham MP: Functional connectivity of human striatum: a resting state FMRI study. Cereb Cortex. 2008, 18 (12): 2735-2747. 10.1093/cercor/bhn041.View ArticlePubMed
- Damoiseaux JS, Beckmann CF, Arigita EJ, Barkhof F, Scheltens P, Stam CJ, Smith SM, Rombouts SA: Reduced resting-state brain activity in the "default network" in normal aging. Cereb Cortex. 2008, 18 (8): 1856-1864. 10.1093/cercor/bhm207.View ArticlePubMed
- Logothetis NK, Murayama Y, Augath M, Steffen T, Werner J, Oeltermann A: How not to study spontaneous activity. NeuroImage. 2009, 45 (4): 1080-1089. 10.1016/j.neuroimage.2009.01.010.View ArticlePubMed
- Kandel ER, Schwartz JH, Jessell TM: Principles of Neural Science. 2000, New York: McGraw-Hill, 4
- Kunzle H: An autoradiographic analysis of the efferent connections from premotor and adjacent prefrontal regions (areas 6 and 9) in macaca fascicularis. Brain, behavior and evolution. 1978, 15 (3): 185-234. 10.1159/000123779.View ArticlePubMed
- Alexander GE, DeLong MR, Strick PL: Parallel organization of functionally segregated circuits linking basal ganglia and cortex. Annu Rev Neurosci. 1986, 9: 357-381. 10.1146/annurev.ne.09.030186.002041.View ArticlePubMed
- Yelnik J: Functional anatomy of the basal ganglia. Mov Disord. 2002, 17 (Suppl 3): S15-21. 10.1002/mds.10138.View ArticlePubMed
- Robinson S, Moser E, Peper M: fMRI of Emotion. fMRI Techniques and Protocols. Edited by: Filippi M. 2009, New Jersey: Humana Press, 41.View Article
- Robinson S, Pripfl J, Bauer H, Moser E: The impact of EPI voxel size on SNR and BOLD sensitivity in the anterior medio-temporal lobe: a comparative group study of deactivation of the Default Mode. MAGMA. 2008, 21 (4): 279-290. 10.1007/s10334-008-0128-0.View ArticlePubMed
- Rombouts SA, Damoiseaux JS, Goekoop R, Barkhof F, Scheltens P, Smith SM, Beckmann CF: Model-free group analysis shows altered BOLD FMRI networks in dementia. Hum Brain Mapp. 2009, 30 (1): 256-266. 10.1002/hbm.20505.View ArticlePubMed
- Behrens TE, Johansen-Berg H, Woolrich MW, Smith SM, Wheeler-Kingshott CA, Boulby PA, Barker GJ, Sillery EL, Sheehan K, Ciccarelli O, et al.: Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat Neurosci. 2003, 6 (7): 750-757. 10.1038/nn1075.View ArticlePubMed
- Fransson P: Spontaneous low-frequency BOLD signal fluctuations: an fMRI investigation of the resting-state default mode of brain function hypothesis. Hum Brain Mapp. 2005, 26 (1): 15-29. 10.1002/hbm.20113.View ArticlePubMed
- Hampson M, Peterson BS, Skudlarski P, Gatenby JC, Gore JC: Detection of functional connectivity using temporal correlations in MR images. Hum Brain Mapp. 2002, 15 (4): 247-262. 10.1002/hbm.10022.View ArticlePubMed
- Slavin KV, Thulborn KR, Wess C, Nersesyan H: Direct visualization of the human subthalamic nucleus with 3T MR imaging. AJNR Am J Neuroradiol. 2006, 27 (1): 80-84.PubMed
- Zhang D, Snyder AZ, Fox MD, Sansbury MW, Shimony JS, Raichle ME: Intrinsic functional relations between human cerebral cortex and thalamus. J Neurophysiol. 2008, 100 (4): 1740-1748. 10.1152/jn.90463.2008.PubMed CentralView ArticlePubMed
- Gelman N, Gorell JM, Barker PB, Savage RM, Spickler EM, Windham JP, Knight RA: MR imaging of human brain at 3.0 T: preliminary report on transverse relaxation rates and relation to estimated iron content. Radiology. 1999, 210 (3): 759-767.View ArticlePubMed
- Calhoun VD, Adali T, Pekar JJ, Pearlson GD: Latency (in)sensitive ICA. Group independent component analysis of fMRI data in the temporal frequency domain. NeuroImage. 2003, 20 (3): 1661-1669. 10.1016/S1053-8119(03)00411-7.View ArticlePubMed
- Gusnard DA, Raichle ME, Raichle ME: Searching for a baseline: functional imaging and the resting human brain. Nat Rev Neurosci. 2001, 2 (10): 685-694. 10.1038/35094500.View ArticlePubMed
- Uddin LQ, Mooshagian E, Zaidel E, Scheres A, Margulies DS, Kelly AM, Shehzad Z, Adelstein JS, Castellanos FX, Biswal BB, et al.: Residual functional connectivity in the split-brain revealed with resting-state functional MRI. Neuroreport. 2008, 19 (7): 703-709. 10.1097/WNR.0b013e3282fb8203.PubMed CentralView ArticlePubMed
- Handwerker DA, Ollinger JM, D'Esposito M: Variation of BOLD hemodynamic responses across subjects and brain regions and their effects on statistical analyses. NeuroImage. 2004, 21 (4): 1639-1651. 10.1016/j.neuroimage.2003.11.029.View ArticlePubMed
- Chang C, Thomason ME, Glover GH: Mapping and correction of vascular hemodynamic latency in the BOLD signal. NeuroImage. 2008, 43 (1): 90-102. 10.1016/j.neuroimage.2008.06.030.PubMed CentralView ArticlePubMed
- Friston K: Causal modelling and brain connectivity in functional magnetic resonance imaging. PLoS biology. 2009, 7 (2): e33-10.1371/journal.pbio.1000033.View ArticlePubMed
- Zarahn E, Aguirre GK, D'Esposito M: Empirical analyses of BOLD fMRI statistics. I. Spatially unsmoothed data collected under null-hypothesis conditions. NeuroImage. 1997, 5 (3): 179-197. 10.1006/nimg.1997.0263.View ArticlePubMed
- Soldati N, Robinson S, Persello C, Jovicich J, Bruzzone L: Automatic classification of brain resting states using fMRI temporal signals. Electronics Letters. 2009, 45 (1): 19-21. 10.1049/el:20092178.View Article
- Scholz VH, Flaherty AW, Kraft E, Keltner JR, Kwong KK, Chen YI, Rosen BR, Jenkins BG: Laterality, somatotopy and reproducibility of the basal ganglia and motor cortex during motor tasks. Brain Res. 2000, 879 (1-2): 204-215. 10.1016/S0006-8993(00)02749-9.View ArticlePubMed
- Hikosaka O, Sakamoto M, Usui S: Functional properties of monkey caudate neurons. I. Activities related to saccadic eye movements. J Neurophysiol. 1989, 61 (4): 780-798.PubMed
- Hikosaka O, Wurtz RH: The basal ganglia. Reviews of oculomotor research. 1989, 3: 257-281.PubMed
- Mohler CW, Goldberg ME, Wurtz RH: Visual receptive fields of frontal eye field neurons. Brain Res. 1973, 61: 385-389. 10.1016/0006-8993(73)90543-X.View ArticlePubMed
- Mallet N, Ballion B, Le Moine C, Gonon F: Cortical inputs and GABA interneurons imbalance projection neurons in the striatum of parkinsonian rats. J Neurosci. 2006, 26 (14): 3875-3884. 10.1523/JNEUROSCI.4439-05.2006.View ArticlePubMed
- Lalo E, Thobois S, Sharott A, Polo G, Mertens P, Pogosyan A, Brown P: Patterns of bidirectional communication between cortex and basal ganglia during movement in patients with Parkinson disease. J Neurosci. 2008, 28 (12): 3008-3016. 10.1523/JNEUROSCI.5295-07.2008.View ArticlePubMed
- Kolomiets BP, Deniau JM, Glowinski J, Thierry AM: Basal ganglia and processing of cortical information: functional interactions between trans-striatal and trans-subthalamic circuits in the substantia nigra pars reticulata. Neuroscience. 2003, 117 (4): 931-938. 10.1016/S0306-4522(02)00824-2.View ArticlePubMed
- Dejean C, Gross CE, Bioulac B, Boraud T: Dynamic changes in the cortex-Basal Ganglia network after dopamine depletion in the rat. J Neurophysiol. 2008, 100 (1): 385-396. 10.1152/jn.90466.2008.View ArticlePubMed
- Chen CC, Litvak V, Gilbertson T, Kuhn A, Lu CS, Lee ST, Tsai CH, Tisch S, Limousin P, Hariz M, et al.: Excessive synchronization of basal ganglia neurons at 20 Hz slows movement in Parkinson's disease. Experimental neurology. 2007, 205 (1): 214-221. 10.1016/j.expneurol.2007.01.027.View ArticlePubMed
- Williams D, Tijssen M, Van Bruggen G, Bosch A, Insola A, Di Lazzaro V, Mazzone P, Oliviero A, Quartarone A, Speelman H, et al.: Dopamine-dependent changes in the functional connectivity between basal ganglia and cerebral cortex in humans. Brain. 2002, 125 (Pt 7): 1558-1569. 10.1093/brain/awf156.View ArticlePubMed
- Papinutto N, Jovicich J: Optimization of brain tissue contrast in structural images at 4T: a computer simulation and validation study. 25th Congress of the European Society for Magnetic Resonance in Medicine and Biology. Valencia, Spain. 2008
- Zaitsev M, Hennig J, Speck O: Point spread function mapping with parallel imaging techniques and high acceleration factors: fast, robust, and flexible method for echo-planar imaging distortion correction. Magn Reson Med. 2004, 52 (5): 1156-1166. 10.1002/mrm.20261.View ArticlePubMed
- Robinson S, Jovicich J: EPI distortion corrections at 4T: multi-channel field mapping and a comparison with the point-spread function method. Proceedings of the Sixteenth Annual Meeting of the ISMRM, Toronto. 2008
- Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, et al.: Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage. 2004, 23 (Supplement 1): S208-219. 10.1016/j.neuroimage.2004.07.051.View ArticlePubMed
- Passingham RE, Stephan KE, Kotter R: The anatomical basis of functional localization in the cortex. Nat Rev Neurosci. 2002, 3 (8): 606-616.View ArticlePubMed
- Oppenheim AVaS RW: Discrete-Time Signal Processing. 1989, Englewood Cliffs, NJ: Prentice-Hall
- Ridler TW, Calvard S: Picture thresholding using an iterative selection method. IEEE Transactions on Systems, Man and Cybernetics. 1978, 8: 630-632. 10.1109/TSMC.1978.4310039.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.