In vivo anatomical and functional identification of V5/MT using high-resolution MRI: a technique for relating structure and function in the human cerebral cortex
Dr. Stephanie Barnes, Prof. Gary Egan & Prof. John Watson
Wednesday, April 1st, 2015
Previous in vivo neuroimaging studies have clearly demonstrated the functional specialisation of the human cerebral cortex. However, precise anatomical localisation of functionally defined cortical areas is an ongoing challenge due to the poor spatial resolution of functional imaging techniques and significant inter-individual differences in the complex morphological structure of the human cortex. The present study used high-resolution MRI to identify V5/MT in three subjects based on its distinctive MR-visible myeloarchitectonic structure. Consistent with previous studies, V5/MT was localised to the junction of the ascending limb of the inferior temporal sulcus and the lateral occipital sulcus. This anatomically defined location of V5/MT was shown to correspond with its functionally defined location, identified using fMR in one subject. Structural MR images with high spatial resolution were acquired in this study by combining increased MR field strength, a multi-channel phased-array head coil for image acquisition and signal averaging across a series of T1-weighted images. This study thus confirmed that MR contrast can be used to resolve intracortical lamination known to be present on a histological level, enabling cortical substructure to be visualised in vivo. It provided proof of concept in a single human subject; therefore, further validation of this novel technique for identification of V5/MT and other functionally defined cortical areas is required. Application of this methodology in its own right, or integrated with other MR-based neuroanatomical mapping techniques, will facilitate structure-function correlation throughout neocortex in living human subjects.
Functional specialisation of human cerebral cortex has been demonstrated in vivo using positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) [e.g. 1]. Current research is aimed at precise anatomical localisation of functionally defined areas. However, accurate localisation of active regions is an ongoing challenge given the poor spatial resolution of functional imaging and inter-individual variability in the complex morphological structure of cerebral cortex [2,3]. Consequently, location of cortical areas must be examined within an individual to obtain precise structure-function information.
Various techniques have been developed to indirectly correlate structure and function in cerebral cortex. Post-mortem human brain analysis following functional imaging would be ideal for precise anatomical localisation of functional regions within an individual; however, this is not readily available. Consequently, novel techniques were required to enable structure-function correlation in living humans.
Early research relied on presumed homology between human and non-human primate brains. Unfortunately, comparative variations in brain size, complexity, orientation and potential rearrangement of cortical areas between species makes comparison of functional anatomy difficult .
Brain atlases developed using traditional histological techniques to provide standardised coordinates for neuroanatomical landmarks have also been widely used. However, significant inter-individual variation in brain topography limits their usefulness in precisely identifying functional regions. For example, Talairach and Tournoux’s  atlas is based on one brain’s structure, over which Brodmann’s cytoarchitectonic map  was projected. Given that functionally defined areas can vary in location by centimetres, such an atlas can provide only gross localisation. Additionally, Brodmann’s two-dimensional map contained no data about the intrasulcal surface, two-thirds of cortical surface, so we can only estimate these borders [7,8].
Humans who have had a stroke, tumour, or traumatic brain injury also provide localisation information. Correlation of lesion location and subsequent neurological deficits provides information about the damaged region’s role. However, this technique’s usefulness is limited because (1) lesions are often extensive so accurately locating the area responsible for the missing function is difficult and (2) rarely is only a single function lost, because lesions tend to incorporate areas performing a range of roles .
Recent studies have revealed a successful new method of achieving precise anatomical localisation of functionally identified areas. Co-registration of structural MR with functional images from the same individual enables functionally defined regions to be mapped onto the specific morphological structure of each subject [10,11,12,13]. This approach to structure-function correlation relies on distinct anatomical features, in particular MR-visible cortical myeloarchitecture, e.g. the densely myelinated stria of Gennari, which demarcates primary visual cortex [14,15]. These techniques were used here to identify the visual motion area, V5/MT.
V5/MT is readily identifiable histologically because of its characteristic myelination. Clarke and Miklossy  first identified putative V5/MT in post-mortem brains at the occipito-temporal junction featuring distinctive myelination. These MR-visible myelin bands in layers I, IV and V (the latter two, the external and internal bands of Baillarger) and radial fibres crossing layers VI to IV help anatomically localise putative V5/MT at the junction of ascending limb of inferior temporal sulcus (ALITS) and lateral occipital sulcus (LOS) in parieto-temporo-occipital cortex [11,13]. This is largely consistent with functional results, showing that this functionally defined area also has anatomical identifiers.
V5/MT has been studied in non-human and human primates. Its location and role in non-human primates is similar, although not identical, to humans [17,18], as expected given the limitations of such studies, discussed above. Research in brain-damaged humans also supports V5/MT’s role in visual motion detection [19,20,21,22]; however, lack of imaging or post-mortem analysis has prevented accurate lesion localisation. Even using transcranial magnetic stimulation to mimic deficits experimentally [23,24] is limited by the extensive area such lesions encompass, likely responsible for a range of functions. Thus, structure-function correlation based on MR imaging in healthy humans represents significant progress in the field.
Recent studies have produced a robust non-invasive method of identifying functionally defined cortical areas in living humans [e.g. 11]. However, acquiring MR images with sufficient spatial resolution to precisely characterise underlying microarchitecture is an ongoing challenge. Our approach focused on improving resolution by applying advanced technology and analysis now available.
Precise anatomical localisation of functionally defined cortical areas relies on correlation of functional and structural images. Visualisation of cortical lamination requires a minimum MR resolution of 200-300m, since the thickest myelin band is 250m thick . Standard T1-weighted images are generally acquired at 1-1.5mm3 resolution – it would take several hours to acquire a single T1-weighted image of 200-300m3 with equivalent signal intensity and signal-to-noise ratio (SNR). It is clearly unfeasible to expect subjects to remain stationary within the scanner for this time.
Several ways of increasing spatial resolution of T1-weighted images without relying on long scan times were used here. Firstly, the scanner’s magnetic field strength was increased from 1.5 to 3 Tesla (T), nearly doubling available signal, producing a corresponding SNR improvement , thus enhancing resolution. Secondly, a multi-channel phased-array head coil was used for image acquisition, increasing SNR using radiofrequency coils and combining multiple coils with individual receiver channels into an array covering the same volume as a larger coil with slight sensitive volume overlap, producing signal of equivalent amplitude but greatly reduced noise . This improvement in SNR is particularly evident in surface regions , such as V5/MT. Thirdly, averaging signal across a series of T1-weighted images significantly reduced noise, improving SNR and enhancing visibility of fine architectonic detail , a technique validated by Walters et al. . Additionally, voxel oversampling during averaging due to jitter from small inter-scan head movements increased signal and reduced partial volume error of single scans, improving neuroanatomical detail .
The current study thus applied a well-established methodology of non-invasive in vivo structural identification of functionally defined cortical areas developed by Walters et al.  to the precise anatomical localisation of V5/MT. We aimed to (1) validate this technique and (2) use technological advances to enhance MR images and improve microarchitecture detection. We hypothesised that this would produce improved spatial resolution of structural MR images, enhancing visualisation of cortical lamination within V5/MT.
Experiment 1: Identification of cortical microarchitecture using high-resolution MRI
Ethical approval was obtained. Three healthy male subjects (mean age 52) participated with informed consent.
Twelve to sixteen high-resolution T1-weighted part brain images and three to four whole brain images were acquired over several sessions on a Siemens Trio 3T scanner using a 32-channel phased-array head coil (Siemens AG, Germany). Part brain parameters: three-dimensional magnetisation prepared rapid gradient echo (3D MP-RAGE) sequence: slices=144; thickness=0.75mm; field of view (FOV)=220mm; in-plane resolution=0.5×0.5mm2; echo time (TE)=3.41ms; repeat time (TR)=1800ms; inversion time (TI)=900ms; flip angle (FA)=9º; number of excitations (NEX)=1. Whole brain parameters: 3D MP-RAGE sequence: slices=256; thickness=0.60 mm; FOV=265 mm; in-plane resolution=0.6×0.6mm2; TE=2.81ms; TR=1900ms; TI=900ms; FA=9º; NEX=1. Raw images were transferred via a DICOM client program (Digital Jacket, Hewlett-Packard, CA), composed into contiguous volumes, and saved in Analyze (Radiological) format (Biomedical Imaging Resource, Mayo Foundation, MN).
Images were analysed using tools from Oxford Centre for fMRI of the Brain (FMRIB) Software library [29,30] and Mricro . T1-weighted images were cropped at rostral spinal cord using Mricro  and automatically segmented to remove non-brain tissue using the FMRIB Brain Extraction Tool (BET) . Each image was resampled at half the acquired voxel dimensions, producing volumes with voxel dimensions of 0.25×0.38×0.25mm3 for part and 0.30×0.30×0.30mm3 for whole brain images.
One part brain with minimal motion artefact was made the template for each scanning session. All images obtained in that session were registered to the template using a rigid body model with six degrees of freedom with the FMRIB Linear Image Registration Tool . Template images for each session were then registered to the template for session one. Transformation matrices were concatenated and applied to each image. Transformed images were averaged using fslmaths , producing a mean high-resolution image.
A single whole brain T1-weighted image was acquired in each scanning session. Each of these images was registered to that acquired in session one. Transformation matrices were applied to each and transformed images averaged using fslmaths, producing a mean high-resolution image. SNR was calculated prior to and following image co-registration with regional intensity measured using ImageJ version 1.45 .
V5/MT’s site was estimated based on its postulated location at the ALITS and LOS intersection. Slices through this region were identified and two-dimensional cortical lamination analysis conducted. Intensity line profiles were manually generated using ImageJ version 1.45 . The number, intensity and relative location of each stationary point or point of inflection between cortical surface and grey-white matter boundary were calculated. These measurements were used to generate a cortical lamination map and enabled comparison of lamination between subjects.
Experiment 2: In vivo structural identification of V5/MT using high-resolution MRI
Subject 2 was previously involved in Walters et al. .
fMR data for subject 2 was obtained from Walters et al. . Images were acquired on a 1.5T scanner (Signa Echospeed, General Electric). Subjects observed a moving checkerboard stimulus . Further details available in Walters et al. .
Functional analysis was carried out by Walters et al.  using FLIRT  and SPM99 [35,36]. High-resolution T1-weighted anatomical image obtained for subject 2 was aligned with the average greyscale-normalised surface coil T1-weighted image acquired by Walters et al. . The functional activation map was overlaid and used to identify functionally defined V5/MT for comparison with the location of anatomically defined putative V5/MT from Experiment 1.
Experiment 1: Identification of cortical microarchitecture using high-resolution MRI
Figure 1 shows three T1-weighted MR slices for subject 2. Figure 1(i) shows the raw T1-weighted MR image, prior to de-skulling. The second panel (Fig 1(ii)) is from a single T1-weighted MR image; the third panel (Figure 1(iii)) shows the effect of co-registering multiple T1-weighted MR images within and across scanning sessions. The significant increase in spatial resolution produced by averaging is clearly evident. Quantitatively, this reflects a 34% improvement in SNR due to co-registration.
Slices through putative anatomically defined V5/MT in the co-registered high-resolution structural T1-weighted MR images for all subjects are shown in Figure 2 (left). Visual examination of the areas of interest (Figure 2, middle) showed two horizontally oriented bands within the cortical ribbon. Intensity line profile analysis through putative V5/MT at this point (AB) is also shown in Figure 2 (right). This enabled quantification of the location of these bands. The first was close to the cortical surface, at 30% grey matter depth while the second was at 65% cortical thickness. The intensity line profile for subject 2 (Figure 2 (vi)) also suggested a third band near the grey-white matter junction (90% depth), consistently identified as either a local maximum or point of inflection.
Intensity line profiles of putative V5/MT are clearly different to surrounding cortex. The opposite sulcal bank is characterised by a single, wide peak at 45% cortical thickness (Figure 3).
Experiment 2: In vivo structural identification of V5/MT using high-resolution MRI
Functionally-defined V5/MT for subject 2 was identified on the high-resolution structural T1-weighted MR image acquired in this study using functional data from Walters et al. . This region sits at the junction of ALITS and LOS. Other areas of activation correspond largely with lower order visual areas; there are also areas of parietal activation. Figure 4(i) shows a high-resolution slice through subject 2’s brain with overlaid foci of functional activation. Figure 4(ii) is an enlarged view of the region, with functional activation overlying the location of putative V5/MT identified in Experiment 1. These functional results correlate strongly with the spatial location of putative V5/MT. Characteristic intensity line profiles are thus effective anatomical identifiers for V5/MT.
This study employed a well-established methodology of non-invasive in vivo identification of functionally defined cortical areas to determine V5/MT’s precise location. Efficacy of the technique pioneered by Watson et al.  and extended to fMRI by Walters et al.  was confirmed. High-resolution structural MR images were successfully acquired and intensity line profiles drawn through putative V5/MT. Co-registration of functional data from Walters et al.  with new high-resolution structural data confirmed the putative anatomical location of V5/MT at the junction of ALITS and LOS in subject 2, consistent with previous studies [10,11,13].
Intensity line profile analysis has been used previously [37,11] to quantify cortical lamination. Current results are consistent with previous findings, demonstrating light-coloured bands at 30% and 65% cortical depth. Subject 2’s data suggest a third band near the grey-white matter junction. The first two likely correspond to the heavily myelinated internal and external bands of Baillarger identified in putative V5/MT in post-mortem brains . The origin of the third band may correspond with radial fibres traversing lower cortical layers, also described in post-mortem specimens. These results thus confirm V5/MT’s characteristic T1-weighted MR appearance. Intensity line profiles distinguish putative V5/MT from surrounding cortex. Profiles of the opposite sulcal bank are characterised by a single, wide peak at 45% depth (Figure 3), consistent with lamination described by Walters et al. .
There is a strong correlation between putative V5/MT’s spatial location and functional results for subject 2. This confirms that the characteristic intensity line profiles from this region are distinctive anatomical identifiers for human V5/MT, verifying Walters et al. . Quantitative comparison of the intensity of line profile maxima through V5/MT here and in Walters et al. was not possible due to different techniques used to derive cortical lamination maps. Qualitatively, however, the peaks of local maxima (Figure 2) are less well defined than those previously obtained. This reflects reduced SNR with the present paradigm.
On a whole brain level, high-resolution T1-weighted MR images acquired in this study clearly have superior spatial resolution compared to those in Walters et al.  However, this overall improvement in resolution is at the expense of reduced SNR at the site of interest i.e. putative V5/MT. This is reflected qualitatively as less well-defined maxima on intensity line profiles. In comparison, images obtained previously have high resolution at V5/MT, but lower resolution across remaining cortex. This represents a significant benefit of using surface coils placed directly over the cortex compared to the current paradigm. The potential for concurrent use of these techniques warrants exploration.
This novel approach to improving MR spatial resolution combined technological advances with innovative image analysis. MR scanning hardware is improving rapidly with improved access to high magnetic field strength scanners and multi-channel phased-array head coils. The technique of signal averaging across images from successive sessions is also relatively recent. We demonstrated that these new tools can be successfully applied to produce T1-weighted MR images with high SNR, enabling detection of cortical lamination.
The application of this new approach to in vivo structure-function correlation in the human cortex requires further validation. This study provides proof of concept in a single subject. High-resolution structural MR images obtained for two others clearly demonstrate the efficacy of this approach in enabling visualisation of cortical microarchitecture; unfortunately, functional identification of V5/MT in these subjects was not possible. Future extension of this study will help resolve this issue.
Imaging techniques used here significantly reduced analysis complexity. Twin surface coil arrangements used by Walters et al.  with narrow fields of view resulted in complicated signal attenuation requiring correction prior to registration. This processing step was not required here. Further, use of a multi-channel phased-array head coil also enabled examination of subcortical microarchitecture since high SNR is not restricted to cortical surface. Acquisition of images with high SNR over a larger field of view represents significant progress.
As suggested by Walters et al. , this technique can be broadly applied to identification of cortical architecture in living humans. Its usefulness has not yet been fully explored with studies restricted to V1 and V5/MT [14,11,12]. Given the strong contribution of myeloarchitecture to MR signal, further studies investigating non-visual cortex should focus on Flechsig’s37 fields of increased myelination to maximise initial success.
Development of imaging and analysis techniques enabling visualisation of cortical lamination opens up new research areas. For example, if functionally active areas are well-characterised microanatomically in vivo using high-resolution MR, major input and output cortical layers can be identified. This would require new task paradigms with multiple conditions activating a functional area in different ways, thus involving distinct pathways . Additionally, anatomical localisation of functionally defined areas could guide medical therapy, like that achieved with deep brain stimulation in Parkinson’s disease . Further research in these areas is required.
This study confirmed that MR contrast can resolve intracortical lamination present histologically, enabling visualisation of cortical substructure in vivo. It employed improved MR hardware and analysis to validate Walters et al. , including V5/MT’s characteristic MR profile, and identified further microarchitectonic detail. Further optimisation of techniques to improve lamination detection is required to maximise results. Application of this methodology alone, or integrated with other MR-based mapping, will facilitate structure-function correlation throughout the neocortex in living humans.
I am very grateful for the training in MR analysis provided by Jennifer Leech, Ayaka Ando, Alicia Dymowski and Dr Marcus Gray, technical support from Sina Sadeghi and Dr Neil Killeen, and assistance with MR imaging from Michael Kean. This research was supported by funds from the Florey Neuroscience Institutes, Melbourne. I was also grateful to receive the Dr Charlotte Last Memorial Scholarship for Honours in Medicine from the Sydney Medical School.
Conflict of Interest
S Barnes: firstname.lastname@example.org
 Zeki, S., Watson, J.D., Lueck, C.J., Friston, K.J., Kennard, C., & Frackowiak, R.S. (1991). A direct demonstration of functional specialization in human visual cortex. The Journal of Neuroscience, 11(3), 641-649.
 Iaria, G., & Petrides, M. (2007). Occipital sulci of the human brain: variability and probability maps. The Journal of Comparative Neurology, 501, 243-259.
 White, L. E., Andrews, T. J., Hulette, C., Richards, A., Groelle, M., Paydarfar, J., & Purves, D. (1997). Structure of the human sensorimotor system. I: Morphology and cytoarchitecture of the central sulcus. Cerebral cortex, 7, 18-30.
 Sereno, M. I., & Tootell, R. B. H. (2005). From monkeys to humans: what do we now know about brain homologies? Current Opinion in Neurobiology, 15, 135-144.
 Talairach, J., & Tournoux, P. (1988). Co-planar stereotaxic atlas of the human brain: 3-dimensional proportional system: an approach to cerebral imaging. New York: Thieme Medical Publishers.
 Brodmann, K. (1909). Vergleichende Lokalisationslehre der GroBhirnrinde in ihren Prinzipien dargestellt auf Grund des Zellenbaues. Leipzig: Barth.
 Zilles, K., & Amunts, K. (2010). Centenary of Brodmann’s map – conception and fate. Nature Reviews Neuroscience, 11.
 Zilles, K., Armstrong, E., Schleicher, A., & Kretschmann, H.-J. (1988). The human pattern of gyrification in the cerebral cortex. Anatomy and Embryology, 179, 173-179.
 Bridge, H., & Clare, S. (2006). High-resolution MRI: in vivo histology? Philosophical Transactions of the Royal Society B: Biological Sciences, 361, 137-146.
 Dumoulin, S., Bittar, R., Kabani, N., Baker, C., Le Goualher, G., Pike, G., & Evans, A. (2000). A New Anatomical Landmark for Reliable Identification of Human Area V5/MT: a Quantitative Analysis of Sulcal Patterning. Cerebral Cortex, 10, 454-463.
 Walters, N., Egan, G., Kril, J., Kean, M., Waley, P., Jenkinson, M., & Watson, J. (2003). In vivo identification of human cortical areas using high-resolution MRI: An approach to cerebral structure-function correlation. Proceedings of the National Academy of Sciences, 100, 2981-2986.
 Walters, N., Eickhoff, S., Schleicher, A., Zilles, K., Amunts, K., Egan, G., & Watson, J. (2007). Observer-Independent Analysis of High-Resolution MR Images of the Human Cerebral Cortex: In Vivo Delineation of Cortical Areas. Human Brain Mapping, 28, 1-8.
 Watson, J., Myers, R., Frackowiak, R., Hajnal, J., Woods, R., Mazziotta, J., Shipp, S., & Zeki, S. (1993). Area V5 of the Human Brain: Evidence from a Combined Study Using Positron Emission Tomography and Magnetic Resonance Imaging Cerebral Cortex, 3, 79-94.
 Clark, V., Courchesne, E., & Grafe, M. (1992). In vivo Myeloarchitectonic Analysis of Human Striate and Extrastriate Cortex Using Magnetic Resonance Imaging. Cerebral Cortex, 2, 417-424.
 Barbier, E., Marrett, S., Danek, A., Vortmeyer, A., van Gelderen, P., Duyn, J., Bandettini, P., Grafman, J., & Koretsky, A. (2002). Imaging Cortical Anatomy by High-Resolution MR at 3.0T: Detection of the Stripe of Gennari in Visual Area 17. Magnetic Resonance in Medicine, 48, 735-738.
 Clarke, S., & Miklossy, J. (1990). Occipital cortex in man: Organization of callosal connections, related myelo- and cytoarchitecture, and putative boundaries of functional visual areas. The Journal of Comparative Neurology, 298, 188-214.
 Zeki, S. M. (1969). Representation of central visual fields in prestriate cortex of monkey. Brain Research, 14, 271-291.
 Cragg, B. G., & Ainsworth, A. (1969). The topography of the afferent projections in the circumstriate visual cortex of the monkey studied by the Nauta method. VIsion Research, 9, 733-747.
 Goldstein, K., & Gelb, A. (1918). Psychologische Analysen hirnpathologischer Falle auf Grund von Untersuchungen Hirnverletzter. I. Abhandlung. Zur Psychologie des optischen Wahrnehmungs-und Erkennungsvorganges. Zeitschrift fur die gesamte Neurologie und Psychiatrie, 41, 1-
 Potzl, O., & Redlich, E. (1911). Demonstration eines Falles von bilateraler Affektion beider Occipitallappen. Wiener Klinische Wochenschrift, 24, 517-518.
 Zihl, J., Von Cramon, D., & Mai, N. (1983). Selective disturbance of movement vision after bilateral brain damage. Brain, 106, 313-340.
 Riddoch, G. (1917). Dissociation of visual perceptions due to occipital injuries, with especial reference to appreciation of movement. Brain, 40, 15-57.
 Beckers, G., & Homberg, V. (1992). Cerebral visual motion blindness: transitory akinetopsia induced by transcranial magnetic stimulation of human area V5. Proceedings of the Royal Society of London B: Biological Sciences, 249, 173-178.
 Beckers, G., & Zeki, S. (1995). The consequences of inactivating areas V1 and V5 on visual motion perception. Brain, 118, 49-60.
 Schmitt, F., Grosu, D., Mohr, C., Purdy, D., Salem, K., Scott, K. T., & Stoeckel, B. (2004). 3 Tesla-MRT: Der Erfolg hoherer Feldstarken. Radiologe, 44, 31-48.
 Bernstein, M. A., King, K. F., & Zhou, Z. J. (2004). Handbook of MRI pulse sequences. Amsterdam: Elsevier: Academic Press.
 de Zwart, J. A., Ledden, P. J., van Gelderen, P., Bodurka, J., Chu, R., & Duyn, J. H. (2004). Signal-to-noise ratio and parallel imaging performance of a 16-channel receive-only brain coil array at 3.0 Tesla. Magnetic Resonance in Medicine, 51, 22-26.
 Holmes, C. J., Hoge, R., Collins, L., Woods, R., Toga, A. W., & Evans, A. C. (1998). Enhancement of MR images using registration for signal averaging. Journal of Computer Assisted Tomography, 22, 324-333.
 Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E. J., Johansen-Berg, H., Bannister, P. R., De Luca, M., Drobnjak, I., Flitney, D. E., Niazy, R., Saunders, J., Vickers, J., Zhang, Y., De Stefano, N., Brady, J. M., & Matthews, P. M. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage, 23, 208-219.
 Woolrich, M. W., Jbabdi, S., Patenaude, B., Chappell, M., Makni, S., Behrens, T., Beckmann, C., Jenkinson, M., & Smith, S. M. (2009). Bayesian analysis of neuroimaging data in FSL. Neuroimage, 45, S173-186.
 Rorden, C., & Brett, M. (2000). Stereotaxic display of brain lesions. Behavioural Neurology, 12, 191-200.
 Smith, S. M. (2002). Fast robust automated brain extraction. Human Brain Mapping, 17, 143-155.
 Jenkinson, M., & Smith, S. M. (2001). A global optimisation method for robust affine registration of brain images. Medical Image Analysis, 5, 143-156.
 Abramoff, M. D., Magalhaes, P. J., & Ram, S. J. (2004). Image processing with ImageJ. Biophotonics International, 11, 36-42.
 Friston, K. J., Holmes, A. P., & Worsley, K. J. (1999). How many subjects constitute a study? Neuroimage, 10, 1-5.
 Friston, K. J., Holmes, A. P., Worsley, K. J., Poline, J.-P., Frith, C. D., & Frackowiak, R. S. J. (1995). Statistical parametric maps in functional imaging: a general linear approach. Human Brain Mapping, 2, 189-210.
 Fatterpekar, G. M., Naidich, T. P., Delman, B. N., Aguinaldo, J. G., Gultekin, S. H., Sherwood, C. C., Hof, P. R., Drayer, B. P., & Fayad, Z. A. (2002). Cytoarchitecture of the human cerebral cortex: MR microscopy of excised specimens at 9.4 Tesla. American Journal of
Neuroradiology, 23, 1313-1321.
 Flechsig, P. (1920). Anatomie des menschlichen Gehirns und Ruckenmarks auf myelogenetischer Grundlage. Leipzig: Thieme.
 Turner, R. (2012). Myeloarchitecture – a window for MRI. In Organization for Human Brain Mapping. Beijing, China.
 Benabid, A.L., Chabardes, S., Mitrofanis, J. & Pollak, P. (2009). Deep brain stimulation of the subthalamic nucleus for the treatment of Parkinson’s disease. Lancet Neurology, 8(1): 67-81.