Optimization of tractography of the optic radiations

A detailed spatio-temporal atlas of the white matter tracts for the fetal brain

This study presents the construction of a comprehensive spatiotemporal atlas detailing the development of white matter tracts in the fetal brain using diffusion magnetic resonance imaging (dMRI). Our research leverages data collected from fetal MRI scans conducted between 22 and 37 weeks of gestation, capturing the dynamic changes in the brain's microstructure during this critical period. The atlas includes 60 distinct white matter tracts, including commissural, projection, and association fibers. We employed advanced fetal dMRI processing techniques and tractography to map and characterize the developmental trajectories of these tracts. Our findings reveal that the development of these tracts is characterized by complex patterns of fractional anisotropy (FA) and mean diffusivity (MD), reflecting key neurodevelopmental processes such as axonal growth, involution of the radial-glial scaffolding, and synaptic pruning. This atlas can serve as a useful resource for neuroscience research and clinical practice, improving our understanding of the fetal brain and potentially aiding in the early diagnosis of neurodevelopmental disorders. By detailing the normal progression of white matter tract development, the atlas can be used as a benchmark for identifying deviations that may indicate neurological anomalies or predispositions to disorders.

Parvalbumin as a neurochemical marker of the primate optic radiation

Parvalbumin (PV) is a calcium-binding protein that labels neuronal cell bodies in the magno and parvocellular layers of the primate lateral geniculate nucleus (LGN). Here we demonstrate that PV immunohistochemistry can also be used to trace the optic radiation (OR) of the marmoset monkey (Callithrix jacchus) from its LGN origin to its destinations in the primary visual cortex (V1), thus providing a high-resolution method for identification of the OR with single axon resolution. The emergence of fibers from LGN, their entire course and even the entry points to V1 were clearly defined in coronal, parasagittal, and horizontal sections of marmoset brain. In all cases, the trajectory revealed by PV staining paralleled that defined by high-resolution diffusion tensor imaging (DTI). We found that V1 was the exclusive target for the PV-containing fibers, with abrupt transitions in staining observed in the white matter at the border with area V2, and no evidence of PV-labeled axons feeding into other visual areas. Changes in the pattern of PV staining in the OR were detected following V1 lesions, demonstrating that this method can be used to assess the progress of retrograde degeneration of geniculocortical projections. These results suggest a technically simple approach to advance our understanding of a major white matter structure, which provides a cellular resolution suitable for the detection of microstructural variations during development, health and disease. Understanding the relationship between PV staining and DTI in non-human primates may also offer clues for improving the specificity and sensitivity of OR tractography for clinical purposes.

Tensor- and high-resolution fiber tractography for the delineation of the optic radiation and corticospinal tract in the proximity of intracerebral lesions: a reproducibility and repeatability study

PURPOSE

Fiber tracking (FT) is used in neurosurgical planning for the resection of lesions in proximity to fiber pathways, as it contributes to a substantial amelioration of postoperative neurological impairments. Currently, diffusion-tensor imaging (DTI)-based FT is the most frequently used technique; however, sophisticated techniques such as Q-ball (QBI) for high-resolution FT (HRFT) have suggested favorable results. Little is known about the reproducibility of both techniques in the clinical setting. Therefore, this study aimed to examine the intra- and interrater agreement for the depiction of white matter pathways such as the corticospinal tract (CST) and the optic radiation (OR).

METHODS

Nineteen patients with eloquent lesions in the proximity of the OR or CST were prospectively enrolled. Two different raters independently reconstructed the fiber bundles by applying probabilistic DTI- and QBI-FT. Interrater agreement was evaluated from the comparison between results obtained by the two raters on the same data set acquired in two independent iterations at different timepoints using the Dice Similarity Coefficient (DSC) and the Jaccard Coefficient (JC). Likewise, intrarater agreement was determined for each rater comparing individual results.

RESULTS

DSC values showed substantial intrarater agreement based on DTI-FT (rater 1: mean 0.77 (0.68-0.85); rater 2: mean 0.75 (0.64-0.81); p = 0.673); while an excellent agreement was observed after the deployment of QBI-based FT (rater 1: mean 0.86 (0.78-0.98); rater 2: mean 0.80 (0.72-0.91); p = 0.693). In contrast, fair agreement was observed between both measures for the repeatability of the OR of each rater based on DTI-FT (rater 1: mean 0.36 (0.26-0.77); rater 2: mean 0.40 (0.27-0.79), p = 0.546). A substantial agreement between the measures was noted by applying QBI-FT (rater 1: mean 0.67 (0.44-0.78); rater 2: mean 0.62 (0.32-0.70), 0.665). The interrater agreement was moderate for the reproducibility of the CST and OR for both DSC and JC based on DTI-FT (DSC and JC ≥ 0.40); while a substantial interrater agreement was noted for DSC after applying QBI-based FT for the delineation of both fiber tracts (DSC > 0.6).

CONCLUSIONS

Our findings suggest that QBI-based FT might be a more robust tool for the visualization of the OR and CST adjacent to intracerebral lesions compared with the common standard DTI-FT. For neurosurgical planning during the daily workflow, QBI appears to be feasible and less operator-dependent.

The connectivity-based parcellation of the angular gyrus: fiber dissection and MR tractography study

The angular gyrus (AG) wraps the posterior end of the superior temporal sulcus (STS), so it is considered a continuation of the superior temporal gyrus (STG)/ middle temporal gyrus (MTG) and forms the inferior parietal lobule (IPL) with the supramarginal gyrus (SMG). The AG was functionally divided in the literature, but there is no fiber dissection study in this context. This study divided AG into superior (sAG) and inferior (iAG) parts by focusing on STS. Red, blue silicone-injected eight and four non-silicone-injected human cadaveric cerebrums were dissected via the Klingler method focusing on the AG. White matter (WM) tracts identified during dissection were then reconstructed on the Human Connectome Project 1065 individual template for validation. According to this study, superior longitudinal fasciculus (SLF) II and middle longitudinal fasciculus (MdLF) are associated with sAG; the anterior commissure (AC), optic radiation (OR) with iAG; the arcuate fasciculus (AF), inferior frontooccipital fasciculus (IFOF), and tapetum (Tp) with both parts. In cortical parcellation of AG based on STS, sAG and iAG were associated with different fiber tracts. Although it has been shown in previous studies that there are functionally different subunits with AG parcellation, here, for the first time, other functions of the subunits have been revealed with cadaveric dissection and tractography images.

Associations of body composition with regional brain volumes and white matter microstructure in very preterm infants

OBJECTIVE

To determine associations between body composition and concurrent measures of brain development including (1) Tissue-specific brain volumes and (2) White matter microstructure, among very preterm infants at term equivalent age.

DESIGN

Prospective observational study.

SETTING

Single-centre academic level III neonatal intensive care unit.

PATIENTS

We studied 85 infants born <33 weeks' gestation.

METHODS

At term equivalent age, infants underwent air displacement plethysmography to determine body composition, and brain MRI from which we quantified tissue-specific brain volumes and fractional anisotropy (FA) of white matter tracts. We estimated associations of fat and lean mass Z-scores with each brain outcome, using linear mixed models adjusted for intrafamilial correlation among twins and potential confounding variables.

RESULTS

Median gestational age was 29 weeks (range 23.4-32.9). One unit greater lean mass Z-score was associated with larger total brain volume (10.5 cc, 95% CI 3.8 to 17.2); larger volumes of the cerebellum (1.2 cc, 95% CI 0.5 to 1.9) and white matter (4.5 cc, 95% CI 0.7 to 8.3); and greater FA in the left cingulum (0.3%, 95% CI 0.1% to 0.6%), right uncinate fasciculus (0.2%, 95% CI 0.0% to 0.5%), and right posterior limb of the internal capsule (0.3%, 95% CI 0.03% to 0.6%). Fat Z-scores were not associated with any outcome.

CONCLUSIONS

Lean mass-but not fat-at term was associated with larger brain volume and white matter microstructure differences that suggest improved maturation. Lean mass accrual may index brain growth and development.

Reproducible protocol to obtain and measure first-order relay human thalamic white-matter tracts

The "primary" or "first-order relay" nuclei of the thalamus feed the cerebral cortex with information about ongoing activity in the environment or the subcortical motor systems. Because of the small size of these nuclei and the high specificity of their input and output pathways, new imaging protocols are required to investigate thalamocortical interactions in human perception, cognition and language. The goal of the present study was twofold: I) to develop a reconstruction protocol based on in vivo diffusion MRI to extract and measure the axonal fiber tracts that originate or terminate specifically in individual first-order relay nuclei; and, II) to test the reliability of this reconstruction protocol. In left and right hemispheres, we investigated the thalamocortical/corticothalamic axon bundles linking each of the first-order relay nuclei and their main cortical target areas, namely, the lateral geniculate nucleus (optic radiation), the medial geniculate nucleus (acoustic radiation), the ventral posterior nucleus (somatosensory radiation) and the ventral lateral nucleus (motor radiation). In addition, we examined the main subcortical input pathway to the ventral lateral posterior nucleus, which originates in the dentate nucleus of the cerebellum. Our protocol comprised three components: defining regions-of-interest; preprocessing diffusion data; and modeling white-matter tracts and tractometry. We then used computation and test-retest methods to check whether our protocol could reliably reconstruct these tracts of interest and their profiles. Our results demonstrated that the protocol had nearly perfect computational reproducibility and good-to-excellent test-retest reproducibility. This new protocol may be of interest for both basic human brain neuroscience and clinical studies and has been made publicly available to the scientific community.

Associations of Macronutrient Intake Determined by Point-of-Care Human Milk Analysis with Brain Development among very Preterm Infants

Point-of-care human milk analysis is now feasible in the neonatal intensive care unit (NICU) and allows accurate measurement of macronutrient delivery. Higher macronutrient intakes over this period may promote brain growth and development. In a prospective, observational study of 55 infants born at <32 weeks’ gestation, we used a mid-infrared spectroscopy-based human milk analyzer to measure the macronutrient content in repeated samples of human milk over the NICU hospitalization. We calculated daily nutrient intakes from unfortified milk and assigned infants to quintiles based on median intakes over the hospitalization. Infants underwent brain magnetic resonance imaging at term equivalent age to quantify total and regional brain volumes and fractional anisotropy of white matter tracts. Infants in the highest quintile of energy intake from milk, as compared with the lower four quintiles, had larger total brain volume (31 cc, 95% confidence interval [CI]: 5, 56), cortical gray matter (15 cc, 95%CI: 1, 30), and white matter volume (23 cc, 95%CI: 12, 33). Higher protein intake was associated with larger total brain (36 cc, 95%CI: 7, 65), cortical gray matter (22 cc, 95%CI: 6, 38) and deep gray matter (1 cc, 95%CI: 0.1, 3) volumes. These findings suggest innovative strategies to close nutrient delivery gaps in the NICU may promote brain growth for preterm infants.

Fully automated delineation of the optic radiation for surgical planning using clinically feasible sequences

Quadrantanopia caused by inadvertent severing of Meyer's Loop of the optic radiation is a well-recognised complication of temporal lobectomy for conditions such as epilepsy. Dissection studies indicate that the anterior extent of Meyer's Loop varies considerably between individuals. Quantifying this for individual patients is thus an important step to improve the safety profile of temporal lobectomies. Previous attempts to delineate Meyer's Loop using diffusion MRI tractography have had difficulty estimating its full anterior extent, required manual ROI placement, and/or relied on advanced diffusion sequences that cannot be acquired routinely in most clinics. Here we present CONSULT: a pipeline that can delineate the optic radiation from raw DICOM data in a completely automated way via a combination of robust pre-processing, segmentation, and alignment stages, plus simple improvements that bolster the efficiency and reliability of standard tractography. We tested CONSULT on 696 scans of predominantly healthy participants (539 unique brains), including both advanced acquisitions and simpler acquisitions that could be acquired in clinically acceptable timeframes. Delineations completed without error in 99.4% of the scans. The distance between Meyer's Loop and the temporal pole closely matched both averages and ranges reported in dissection studies for all tested sequences. Median scan-rescan error of this distance was 1 mm. When tested on two participants with considerable pathology, delineations were successful and realistic. Through this, we demonstrate not only how to identify Meyer's Loop with clinically feasible sequences, but also that this can be achieved without fundamental changes to tractography algorithms or complex post-processing methods.

Current and Future Challenges of Functional MRI and Diffusion Tractography in the Surgical Setting: From Eloquent Brain Mapping to Neural Plasticity

Decades ago, Spetzler (1986) and Sawaya (1998) provided a rough brain segmentation of the eloquent areas of the brain, aimed to help surgical decisions in cases of vascular malformations and tumors, respectively. Currently in clinical use, their criteria are in need of revision. Defining functions (eg, sensorimotor, language and visual) that should be preserved during surgery seems a straightforward task. In practice, locating the specific areas that could cause a permanent vs transient deficit is not an easy task. This is particularly true for the associative cortex and cognitive domains such as language. The old model, with Broca's and Wernicke's areas at the forefront, has been superseded by a dual-stream model of parallel language processing; named ventral and dorsal pathways. This complicated network of cortical hubs and subcortical white matter pathways needing preservation during surgery is a work in progress. Preserving not only cortical regions but most importantly preserving the connections, or white matter fiber bundles, of core regions in the brain is the new paradigm. For instance, the arcuate fascicululs and inferior fronto-occipital fasciculus are key components of the dorsal and ventral language pathways, respectively; and their damage result in permanent language deficits. Interestedly, the damage of the temporal portions of these bundles -where there is a crossroad with other multiple bundles-, appears to be more important (permanent) than the damage of the frontal portions - where plasticity and contralateral activation could help. Although intraoperative direct cortical and subcortical stimulation have contributed largely, advanced MR techniques such as functional MRI (fMRI) and diffusion tractography (DT), are at the epi-center of our current understanding. Nevertheless, these techniques posse important challenges: such as neurovascular uncoupling or venous bias on fMRI; and appropriate anatomical validation or accurate representation of crossing fibers on DT. These limitations should be well understood and taken into account in clinical practice. Unifying multidisciplinary research and clinical efforts is desirable, so these techniques could contribute more efficiently not only to locate eloquent areas but to improve outcomes and our understanding of neural plasticity. Finally, although there are constant anatomical and functional regions at the individual level, there is a known variability at the inter-individual level. This concept should strengthen the importance of a personalized approach when evaluating these regions on fMRI and DT. It should strengthen the importance of personalized treatments as well, aimed to meet tailored needs and expectations.

Microsurgical and Tractographic Anatomical Study of Transtemporal-Transchoroidal Fissure Approaches to the Ambient Cistern

BACKGROUND

Approaching ambient cistern lesions is still a challenge because of deep location and related white matter tracts (WMTs) and neural structures.

OBJECTIVE

To investigate the white matter anatomy in the course of 3 types of transtemporal-transchoroidal fissure approaches (TTcFA) to ambient cistern by using fiber dissection technique with translumination and magnetic resonance imaging fiber tractography.

METHODS

Eight formalin-fixed cerebral hemispheres were dissected on surgical corridor from the temporal cortex to the ambient cistern by using Klingler's method. The trans-middle temporal gyrus, trans-inferior temporal sulcus (TITS), and trans-inferior temporal gyrus (TITG) approaches were evaluated. WMTs that were identified during dissection were then reconstructed on the Human Connectome Project 1021 individual template for validation.

RESULTS

The trans-middle gyrus approach interrupted the U fibers, arcuate fasciculus (AF), the ventral segment of inferior frontoocipital fasciculus (IFOF), the temporal extensions of the anterior commissure (AC) posterior crura, the tapetum (Tp) fibers, and the anterior loop of the optic radiation (OR). The TITS approach interrupted U fibers, inferior longitudinal fasciculus (ILF), IFOF, and OR. The TITG approach interrupted the U fibers, ILF, and OR. The middle longitudinal fasciculus, ILF, and uncinate fasciculus (UF) were not interrupted in the trans-middle gyrus approach and the AF, UF, AC, and Tp fibers were not interrupted in the TITS/gyrus approaches.

CONCLUSION

Surgical planning of the ambient cistern lesions requires detailed knowledge about WMTs. Fiber dissection and tractography techniques improve the orientation during surgery and may help decrease surgical complications.

Comparison of multiple tractography methods for reconstruction of the retinogeniculate visual pathway using diffusion MRI

The retinogeniculate visual pathway (RGVP) conveys visual information from the retina to the lateral geniculate nucleus. The RGVP has four subdivisions, including two decussating and two nondecussating pathways that cannot be identified on conventional structural magnetic resonance imaging (MRI). Diffusion MRI tractography has the potential to trace these subdivisions and is increasingly used to study the RGVP. However, it is not yet known which fiber tracking strategy is most suitable for RGVP reconstruction. In this study, four tractography methods are compared, including constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD‐Stream) methods, and multi‐fiber (UKF‐2T) and single‐fiber (UKF‐1T) unscented Kalman filter (UKF) methods. Experiments use diffusion MRI data from 57 subjects in the Human Connectome Project. The RGVP is identified using regions of interest created by two clinical experts. Quantitative anatomical measurements and expert anatomical judgment are used to assess the advantages and limitations of the four tractography methods. Overall, we conclude that UKF‐2T and iFOD1 produce the best RGVP reconstruction results. The iFOD1 method can better quantitatively estimate the percentage of decussating fibers, while the UKF‐2T method produces reconstructed RGVPs that are judged to better correspond to the known anatomy and have the highest spatial overlap across subjects. Overall, we find that it is challenging for current tractography methods to both accurately track RGVP fibers that correspond to known anatomy and produce an approximately correct percentage of decussating fibers. We suggest that future algorithm development for RGVP tractography should take consideration of both of these two points.

Quantitative visual pathway abnormalities predict visual field defects in patients with pituitary adenomas: a diffusion spectrum imaging study

OBJECTIVES

This study was to investigate clinical applicability of diffusion spectrum imaging (DSI) for quantitative detection of visual pathway abnormalities to predict the degree of visual field defects (VFD) in patients with pituitary adenomas.

METHODS

Sixty-five patients with pituitary adenomas and 33 healthy controls underwent conventional MRI and DSI scanning that allowed high-angular-resolution fiber tracking. Optic chiasmal compression and VFD were confirmed in all patients via radiological and neuro-ophthalmological examinations. Quantitative assessments of chiasmal lift, VFD, and DSI parameters from the optic nerve, optic tract, and optic radiation were performed. Group comparisons and correlation analyses were conducted in patients and controls. Using the 5-fold cross-validation method, the support vector machine classifiers were constructed to predict the degree of visual defects.

RESULTS

The mean values of quantitative anisotropy and generalized fractional anisotropy in optic nerve and optic tract showed significant differences between patients and controls (p < 0.05). These parameters were also significantly correlated with the chiasmal lift distance and degree of visual defects (p < 0.05). All patients were divided into mild (n = 42) and severe (n = 23) VFD groups, using the mean deviation value of -8 dB as the threshold. The classifiers achieved an accuracy of 0.83, sensitivity of 0.78, and specificity of 0.86 to discriminate patients with mild and severe visual defects.

CONCLUSIONS

Using high-angular-resolution fiber tracking, DSI may provide quantitative information to detect visual pathway abnormalities and be a potential diagnostic tool for determining the degree of visual field defects in pituitary adenomas.

KEY POINTS

• Abnormal QA and GFA values of optic nerve and optic tract in adenoma patients • Close relationship between DSI parameters and VFD degree in adenoma patients • The classifiers achieved an accuracy of 0.83, sensitivity of 0.78, and specificity of 0.86 to discriminate patients with mild and severe VFD.

Carpal Bone Segmentation Using Fully Convolutional Neural Network

BACKGROUND

Bone Age Assessment (BAA) refers to a clinical procedure that aims to identify a discrepancy between biological and chronological age of an individual by assessing the bone age growth. Currently, there are two main methods of executing BAA which are known as Greulich-Pyle and Tanner-Whitehouse techniques. Both techniques involve a manual and qualitative assessment of hand and wrist radiographs, resulting in intra and inter-operator variability accuracy and time-consuming. An automatic segmentation can be applied to the radiographs, providing the physician with more accurate delineation of the carpal bone and accurate quantitative analysis.

METHODS

In this study, we proposed an image feature extraction technique based on image segmentation with the fully convolutional neural network with eight stride pixel (FCN-8). A total of 290 radiographic images including both female and the male subject of age ranging from 0 to 18 were manually segmented and trained using FCN-8.

RESULTS AND CONCLUSION

The results exhibit a high training accuracy value of 99.68% and a loss rate of 0.008619 for 50 epochs of training. The experiments compared 58 images against the gold standard ground truth images. The accuracy of our fully automated segmentation technique is 0.78 ± 0.06, 1.56 ±0.30 mm and 98.02% in terms of Dice Coefficient, Hausdorff Distance, and overall qualitative carpal recognition accuracy, respectively.

Quantitative assessment of secondary white matter injury in the visual pathway by pituitary adenomas: a multimodal study at 7-Tesla MRI

OBJECTIVE

The objective of this study was to investigate microstructural damage caused by pituitary macroadenomas by performing probabilistic tractography of the optic tracts and radiations using 7-T diffusion-weighted MRI (DWI). These imaging findings were correlated with neuro-ophthalmological results to assess the utility of ultra-high-field MRI for objective evaluation of damage to the anterior and posterior visual pathways.

METHODS

Probabilistic tractography employing 7-T DWI was used to reconstruct the optic tracts and radiations in 18 patients with adenomas and in 16 healthy volunteers. Optic chiasm compression was found in 66.7% of the patients and visual defects in 61.1%. Diffusion indices were calculated along the projections and correlated with tumor volumes and results from neuro-ophthalmological examinations. Primary visual cortical thicknesses were also assessed.

RESULTS

Fractional anisotropy was reduced by 21.9% in the optic tracts (p < 0.001) and 17.7% in the optic radiations (p < 0.001) in patients with adenomas. Patients showed an 8.5% increase in mean diffusivity of optic radiations compared with healthy controls (p < 0.001). Primary visual cortical thickness was reduced in adenoma patients. Diffusion indices of the visual pathway showed significant correlations with neuro-ophthalmological examination findings.

CONCLUSIONS

Imaging-based quantification of secondary neuronal damage from adenomas strongly correlated with neuro-ophthalmological findings. Diffusion characteristics enabled by ultra-high-field DWI may allow preoperative characterization of visual pathway damage in patients with chiasmatic compression and may inform prognosis for vision recoverability.

Implementation of clinical tractography for pre-surgical planning of space occupying lesions: An investigation of common acquisition and post-processing methods compared to dissection studies

Background and purpose

There is limited standardization of acquisition and processing methods in diffusion tractography for pre-surgical planning, leading to a range of approaches. In this study, a number of representative acquisition variants and post processing methods are considered, to assess their importance when implementing a clinical tractography program.

Methods

Diffusion MRI was undertaken in ten healthy volunteers, using protocols typical of clinical and research scanning: a 32-direction diffusion acquisition with and without peripheral gating, and a non-gated 64 diffusion direction acquisition. All datasets were post-processed using diffusion tensor reconstruction with streamline tractography, and with constrained spherical deconvolution (CSD) with both streamline and probabilistic tractography, to delineate the cortico-spinal tract (CST) and optic radiation (OR). The accuracy of tractography results was assessed against a histological atlas using a novel probabilistic Dice overlap technique, together with direct comparison to tract volumes and distance of Meyer’s loop to temporal pole (ML-TP) from dissections studies. Three clinical case studies of patients with space occupying lesions were also investigated.

Results

Tracts produced by CSD with probabilistic tractography provided the greatest overlap with the histological atlas (overlap scores of 44% and 52% for the CST and OR, respectively) and best matched tract volume and ML-TP distance from dissection studies. The acquisition protocols investigated had limited impact on the accuracy of the tractography. In all patients, the CSD based probabilistic tractography created tracts with greatest anatomical plausibility, although in one case anatomically plausible pathways could not be reconstructed without reducing the probabilistic threshold, leading to an increase in false positive tracts.

Conclusions

Advanced post processing techniques such as CSD with probabilistic tractography are vital for pre-surgical planning. However, overall accuracy relative to dissection studies remains limited.

Brain structure in children with congenital visual disorders and visual impairment

AIM

To examine if congenital visual impairment is associated with differences in brain anatomy in children.

METHOD

Ten children (8-12y) with congenital disorders of the peripheral visual system with severe visual impairment (SVI; >0.8 logMAR) or mild-to-moderate visual impairment (MVI; 0.6-0.8 logMAR) were compared to 21 typically sighted comparison (TSC) children. Thalamus volume, grey matter density, white matter microstructure, and integrity of visual tracts were investigated in SVI, MVI, and TSC groups with anatomical and diffusion-weighted magnetic resonance imaging.

RESULTS

Compared to the TSC group, the SVI group had lower white matter integrity in tracts of the visual system (optic radiations: SVI 0.35±0.015, TSC 0.39±0.007 [p=0.022]; posterior corpus callosum: SVI 0.37±0.019; TSC 0.42±0.009 [p=0.033]) and lower left thalamus volume (SVI 4.37±0.087; TSC 4.99±0.339 [p=0.015]). Neuroanatomical differences were greater in the SVI group, while no consistent differences between the MVI and TSC group were observed.

INTERPRETATION

Posterior tracts of the visual system are compromised in children with congenital visual impairment versus those who are typically sighted. The severity of visual input appears to have affected neuroanatomical development as significant reductions were only found in the SVI group.

WHAT THIS PAPER ADDS

Severe visual impairment in mid-childhood is associated with reduced integrity of visual pathways and reduced thalamus volume.

Early white matter development is abnormal in tuberous sclerosis complex patients who develop autism spectrum disorder

Background

Autism spectrum disorder (ASD) is prevalent in tuberous sclerosis complex (TSC), occurring in approximately 50% of patients, and is hypothesized to be caused by disruption of neural circuits early in life. Tubers, or benign hamartomas distributed stochastically throughout the brain, are the most conspicuous of TSC neuropathology, but have not been consistently associated with ASD. Widespread neuropathology of the white matter, including deficits in myelination, neuronal migration, and axon formation, exist and may underlie ASD in TSC. We sought to identify the neural circuits associated with ASD in TSC by identifying white matter microstructural deficits in a prospectively recruited, longitudinally studied cohort of TSC infants.

Methods

TSC infants were recruited within their first year of life and longitudinally imaged at time of recruitment, 12 months of age, and at 24 months of age. Autism was diagnosed at 24 months of age with the ADOS-2. There were 108 subjects (62 TSC-ASD, 55% male; 46 TSC+ASD, 52% male) with at least one MRI and a 24-month ADOS, for a total of 187 MRI scans analyzed (109 TSC-ASD; 78 TSC+ASD). Diffusion tensor imaging properties of multiple white matter fiber bundles were sampled using a region of interest approach. Linear mixed effects modeling was performed to test the hypothesis that infants who develop ASD exhibit poor white matter microstructural integrity over the first 2 years of life compared to those who do not develop ASD.

Results

Subjects with TSC and ASD exhibited reduced fractional anisotropy in 9 of 17 white matter regions, sampled from the arcuate fasciculus, cingulum, corpus callosum, anterior limbs of the internal capsule, and the sagittal stratum, over the first 2 years of life compared to TSC subjects without ASD. Mean diffusivity trajectories did not differ between groups.

Conclusions

Underconnectivity across multiple white matter fiber bundles develops over the first 2 years of life in subjects with TSC and ASD. Future studies examining brain-behavior relationships are needed to determine how variation in the brain structure is associated with ASD symptoms.

Optic Radiation Tractography in Pediatric Brain Surgery Applications: A Reliability and Agreement Assessment of the Tractography Method

Background

Optic radiation (OR) tractography may help predict and reduce post-neurosurgical visual field deficits. OR tractography methods currently lack pediatric and surgical focus.

Purpose

We propose a clinically feasible OR tractography strategy in a pediatric neurosurgery setting and examine its intra-rater and inter-rater reliability/agreements.

Methods

Preoperative and intraoperative MRI data were obtained from six epilepsy and two brain tumor patients on 3 Tesla MRI scanners. Four raters with different clinical experience followed the proposed strategy to perform probabilistic OR tractography with manually drawing anatomical landmarks to reconstruct the OR pathway, based on fiber orientation distributions estimated from high angular resolution diffusion imaging data. Intra- and inter-rater reliabilities/agreements of tractography results were assessed using intraclass correlation coefficient (ICC) and dice similarity coefficient (DSC) across various tractography and OR morphological metrics, including the lateral geniculate body positions, tract volumes, and Meyer's loop position from temporal anatomical landmarks.

Results

Good to excellent intra- and inter-rater reproducibility was demonstrated for the majority of OR reconstructions (ICC = 0.70-0.99; DSC = 0.84-0.89). ICC was higher for non-lesional (0.82-0.99) than lesional OR (0.70-0.99). The non-lesional OR's mean volume was 22.66 cm³; the mean Meyer's loop position was 29.4 mm from the temporal pole, 5.89 mm behind of and 10.26 mm in front of the temporal ventricular horn. The greatest variations (± 1.00-3.00 mm) were observed near pathology, at the tract edges or at cortical endpoints. The OR tractography were used to assist surgical planning and guide lesion resection in all cases, no patient had new visual field deficits postoperatively.

Conclusion

The proposed tractography strategy generates reliable and reproducible OR tractography images that can be reliably implemented in the routine, non-emergency pediatric neurosurgical setting.

Comparative neuroanatomy: Integrating classic and modern methods to understand association fibers connecting dorsal and ventral visual cortex

Comparative neuroanatomy studies improve understanding of brain structure and function and provide insight regarding brain development, evolution, and also what features of the brain are uniquely human. With modern methods such as diffusion MRI (dMRI) and quantitative MRI (qMRI), we are able to measure structural features of the brain with the same methods across human and non-human primates. In this review article, we discuss how recent dMRI measurements of vertical occipital connections in humans and macaques can be compared with previous findings from invasive anatomical studies that examined connectivity, including relatively forgotten classic strychnine neuronography studies. We then review recent progress in understanding the neuroanatomy of vertical connections within the occipitotemporal cortex by combining modern quantitative MRI and classical histological measurements in human and macaque. Finally, we a) discuss current limitations of dMRI and tractography and b) consider potential paths for future investigations using dMRI and tractography for comparative neuroanatomical studies of white matter tracts between species. While we focus on vertical association connections in visual cortex in the present paper, this same approach can be applied to other white matter tracts. Similar efforts are likely to continue to advance our understanding of the neuroanatomical features of the brain that are shared across species, as well as to distinguish the features that are uniquely human.

Reproducibility of Structural and Diffusion Tensor Imaging in the TACERN Multi-Center Study

BACKGROUND

Multi-site MRI studies are often necessary for recruiting sufficiently sized samples when studying rare conditions. However, they require pooling data from multiple scanners into a single data set, and therefore it is critical to evaluate the variability of quantitative MRI measures within and across scanners used in multi-site studies. The aim of this study was to evaluate the reproducibility of structural and diffusion weighted (DW) MRI measurements acquired on seven scanners at five medical centers as part of the Tuberous Sclerosis Complex Autism Center of Excellence Research Network (TACERN) multisite study.

METHODS

The American College of Radiology (ACR) phantom was imaged monthly to measure reproducibility of signal intensity and uniformity within and across seven 3T scanners from General Electric, Philips, and Siemens vendors. One healthy adult male volunteer was imaged repeatedly on all seven scanners under the TACERN structural and DW protocol (5 b = 0 s/mm2 and 30 b = 1000 s/mm2) over a period of 5 years (age 22-27 years). Reproducibility of inter- and intra-scanner brain segmentation volumes and diffusion tensor imaging metrics fractional anisotropy (FA) and mean diffusivity (MD) within white matter regions was quantified with coefficient of variation.

RESULTS

The American College of Radiology Phantom signal intensity and uniformity were similar across scanners and changed little over time, with a mean intra-scanner coefficient of variation of 3.6 and 1.8%, respectively. The mean inter- and intra-scanner coefficients of variation of brain structure volumes derived from T1-weighted (T1w) images of the human phantom were 3.3 and 1.1%, respectively. The mean inter- and intra-scanner coefficients of variation of FA in white matter regions were 4.5 and 2.5%, while the mean inter- and intra-scanner coefficients of variation of MD in white matter regions were 5.4 and 1.5%.

CONCLUSION

Our results suggest that volumetric and diffusion tensor imaging (DTI) measurements are highly reproducible between and within scanners and provide typical variation amplitudes that can be used as references to interpret future findings in the TACERN network.

The development of white matter structural changes during the process of deterioration of the visual field

Emerging evidence suggests that white matter plasticity in the adult brain is preserved after sensory and behavioral modifications. However, little is known about the progression of structural changes during the process of decline in visual input. Here we studied two groups of patients suffering from advanced retinitis pigmentosa with specific deterioration of the visual field: patients who had lost their peripheral visual field, retaining only central (“tunnel”) vision, and blind patients with complete visual field loss. Testing of these homogeneous groups made it possible to assess the extent to which the white matter is affected by loss of partial visual input and whether partially preserved visual input suffices to sustain stability in tracts beyond the primary visual system. Our results showed gradual changes in diffusivity that are indicative of degenerative processes in the primary visual pathway comprising the optic tract and the optic radiation. Interestingly, changes were also found in tracts of the ventral stream and the corticospinal fasciculus, depicting a gradual reorganisation of these tracts consequentially to the gradual loss of visual field coverage (from intact perception to partial vision to complete blindness). This reorganisation may point to microstructural plasticity underlying adaptive behavior and cross-modal integration after partial visual deprivation.

Meyer's loop tractography for image-guided surgery depends on imaging protocol and hardware

Introduction

Surgical resection is an effective treatment for temporal lobe epilepsy but can result in visual field defects. This could be minimized if surgeons knew the exact location of the anterior part of the optic radiation (OR), the Meyer's loop. To this end, there is increasing prevalence of image-guided surgery using diffusion MRI tractography. Despite considerable effort in developing analysis methods, a wide discrepancy in Meyer's loop reconstructions is observed in the literature. Moreover, the impact of differences in image acquisition on Meyer's loop tractography remains unclear. Here, while employing the same state-of-the-art analysis protocol, we explored the extent to which variance in data acquisition leads to variance in OR reconstruction.

Methods

Diffusion MRI data were acquired for the same thirteen healthy subjects using standard and state-of-the-art protocols on three scanners with different maximum gradient amplitudes (MGA): Siemens Connectom (MGA = 300 mT/m); Siemens Prisma (MGA = 80 mT/m) and GE Excite-HD (MGA = 40 mT/m). Meyer's loop was reconstructed on all subjects and its distance to the temporal pole (ML-TP) was compared across protocols.

Results

A significant effect of data acquisition on the ML-TP distance was observed between protocols (p < .01 to 0.0001). The biggest inter-acquisition discrepancy for the same subject across different protocols was 16.5 mm (mean: 9.4 mm, range: 3.7-16.5 mm).

Conclusion

We showed that variance in data acquisition leads to substantive variance in OR tractography. This has direct implications for neurosurgical planning, where part of the OR is at risk due to an under-estimation of its location using conventional acquisition protocols.

Visualization of Cranial Nerves Using High-Definition Fiber Tractography

BACKGROUND

Recent studies have demonstrated diffusion tensor imaging tractography of cranial nerves (CNs). Spatial and angular resolution, however, is limited with this modality. A substantial improvement in image resolution can be achieved with high-angle diffusion magnetic resonance imaging and atlas-based fiber tracking to provide detailed trajectories of CNs.

OBJECTIVE

To use high-definition fiber tractography to identify CNs in healthy subjects and patients with brain tumors.

METHODS

Five neurologically healthy adults and 3 patients with brain tumors were scanned with diffusion spectrum imaging that allowed high-angular-resolution fiber tracking. In addition, a 488-subject diffusion magnetic resonance imaging template constructed from the Human Connectome Project data was used to conduct atlas space fiber tracking of CNs.

RESULTS

The cisternal portions of most CNs were tracked and visualized in each healthy subject and in atlas fiber tracking. The entire optic radiation, medial longitudinal fasciculus, spinal trigeminal nucleus/tract, petroclival portion of the abducens nerve, and intrabrainstem portion of the facial nerve from the root exit zone to the adjacent abducens nucleus were identified. This suggested that the high-angular-resolution fiber tracking was able to distinguish the facial nerve from the vestibulocochlear nerve complex. The tractography clearly visualized CNs displaced by brain tumors. These tractography findings were confirmed intraoperatively.

CONCLUSION

Using high-angular-resolution fiber tracking and atlas-based fiber tracking, we were able to identify all CNs in unprecedented detail. This implies its potential in localization of CNs during surgical planning.

ABBREVIATIONS

CN, cranial nerveDSI, diffusion spectrum imagingDTI, diffusion tensor imagingHCP, Human Connectome ProjectHDFT, high-definition fiber tractographyMLF, medial longitudinal fasciculusODF, orientation distribution functionROI, region of interest.

Stability metrics for optic radiation tractography: Towards damage prediction after resective surgery

•

The alignment of streamlines is quantified by fiber-to-bundle coherence measures.

•

Reliable ML-TP distance measurement by removal of spurious (deviating) streamlines.

•

Parameter estimation to remove spurious streamlines and to retain the Meyer's loop.

•

The validity of ML-TP distance is estimated by pre and postoperative OR comparisons.

•

The stability metrics are promising to relate OR damage to a visual field deficit.

Background

An accurate delineation of the optic radiation (OR) using diffusion MR tractography may reduce the risk of a visual field deficit after temporal lobe resection. However, tractography is prone to generate spurious streamlines, which deviate strongly from neighboring streamlines and hinder a reliable distance measurement between the temporal pole and the Meyer's loop (ML-TP distance).

New method

Stability metrics are introduced for the automated removal of spurious streamlines near the Meyer's loop. Firstly, fiber-to-bundle coherence (FBC) measures can identify spurious streamlines by estimating their alignment with the surrounding streamline bundle. Secondly, robust threshold selection removes spurious streamlines while preventing an underestimation of the extent of the Meyer's loop. Standardized parameter selection is realized through test–retest evaluation of the variability in ML-TP distance.

Results

The variability in ML-TP distance after parameter selection was below 2 mm for each of the healthy volunteers studied (N = 8). The importance of the stability metrics is illustrated for epilepsy surgery candidates (N = 3) for whom the damage to the Meyer's loop was evaluated by comparing the pre- and post-operative OR reconstruction. The difference between predicted and observed damage is in the order of a few millimeters, which is the error in measured ML-TP distance.

Comparison with existing method(s)

The stability metrics are a novel method for the robust estimate of the ML-TP distance.

Conclusions

The stability metrics are a promising tool for clinical trial studies, in which the damage to the OR can be related to the visual field deficit that may occur after epilepsy surgery.

Active delineation of Meyer's loop using oriented priors through MAGNEtic tractography (MAGNET)

Streamline tractography algorithms infer connectivity from diffusion MRI (dMRI) by following diffusion directions which are similarly aligned between neighboring voxels. However, not all white matter (WM) fascicles are organized in this manner. For example, Meyer's loop is a highly curved portion of the optic radiation (OR) that exhibits a narrow turn, kissing and crossing pathways, and changes in fascicle dispersion. From a neurosurgical perspective, damage to Meyer's loop carries a potential risk of inducing vision deficits to the patient, especially during temporal lobe resection surgery. To prevent such impairment, achieving an accurate delineation of Meyer's loop with tractography is thus of utmost importance. However, current algorithms tend to under-estimate the full extent of Meyer's loop, mainly attributed to the aforementioned rule for connectivity which requires a direction to be chosen across a field of orientations. In this article, it was demonstrated that MAGNEtic Tractography (MAGNET) can benefit Meyer's loop delineation by incorporating anatomical knowledge of the expected fiber orientation to overcome local ambiguities. A new ROI-mechanism was proposed which supplies additional information to streamline reconstruction algorithms by the means of oriented priors. Their results showed that MAGNET can accurately generate Meyer's loop in all of our 15 child subjects (8 males; mean age 10.2 years ± 3.1). It effectively improved streamline coverage when compared with deterministic tractography, and significantly reduced the distance between the anterior-most portion of Meyer's loop and the temporal pole by 16.7 mm on average, a crucial landmark used for preoperative planning of temporal lobe surgery. Hum Brain Mapp 38:509-527, 2017. © 2016 Wiley Periodicals, Inc.

Vision after trans-sylvian or temporobasal selective amygdalohippocampectomy: a prospective randomised trial

BACKGROUND

Selective amygdalohippocampectomy (SAH) is an accepted surgical procedure for treatment of pharmacoresistant mesial temporal lobe epilepsy, but it may lead to postoperative visual field deficits (VFDs). Here we present a prospective randomised trial comparing the postoperative VFDs after either a trans-sylvian or temporobasal approach for SAH.

METHOD

Forty-eight patients were randomly assigned to trans-sylvian (n = 24) or temporobasal (n = 24) SAH. Postoperative VFD were quantitatively evaluated using automated static and kinetic perimetry. In 24 cases, diffusion tensor imaging-based deterministic fibre-tracking of the optic radiation was performed. The primary endpoint was absence of postoperative VFD. The secondary endpoint was seizure outcome and driving ability.

RESULTS

Three patients (13 %) from the trans-sylvian group showed no VFD, compared to 11 patients (46 %) from the temporobasal group without VFD (p = 0.01, RR = 3.7; CI = 1.2-11.5). Fifteen patients from each group (63 %) became completely seizure-free (ILAE1). Among those seizure-free cases, five trans-sylvian (33 %) and ten temporobasal (66 %) patients could apply for a driving licence (NNT = 3) when VFDs were considered. Although the trans-sylvian group experienced more frequent VFDs, the mean functional visual impairment showed a tendency to be less pronounced compared with the temporobasal group. DTI-based tracking of the optic radiation revealed that a lower distance of optic radiation to the temporal base correlated with increased rate of VFD in the temporobasal group.

CONCLUSIONS

Temporobasal SAH shows significantly fewer VFDs and equal seizure-free rate compared with the trans-sylvian SAH. However, in patients in whom the optic radiation is close to the temporal base, the trans-sylvian approach may be a preferred alternative.

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

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.

The white matter tracts of the cerebrum in ventricular surgery and hydrocephalus

OBJECTIVE The relationship of the white matter tracts to the lateral ventricles is important when planning surgical approaches to the ventricles and in understanding the symptoms of hydrocephalus. The authors' aim was to explore the relationship of the white matter tracts of the cerebrum to the lateral ventricles using fiber dissection technique and MR tractography and to discuss these findings in relation to approaches to ventricular lesions. METHODS Forty adult human formalin-fixed cadaveric hemispheres (20 brains) and 3 whole heads were examined using fiber dissection technique. The dissections were performed from lateral to medial, medial to lateral, superior to inferior, and inferior to superior. MR tractography showing the lateral ventricles aided in the understanding of the 3D relationships of the white matter tracts with the lateral ventricles. RESULTS The relationship between the lateral ventricles and the superior longitudinal I, II, and III, arcuate, vertical occipital, middle longitudinal, inferior longitudinal, inferior frontooccipital, uncinate, sledge runner, and lingular amygdaloidal fasciculi; and the anterior commissure fibers, optic radiations, internal capsule, corona radiata, thalamic radiations, cingulum, corpus callosum, fornix, caudate nucleus, thalamus, stria terminalis, and stria medullaris thalami were defined anatomically and radiologically. These fibers and structures have a consistent relationship to the lateral ventricles. CONCLUSIONS Knowledge of the relationship of the white matter tracts of the cerebrum to the lateral ventricles should aid in planning more accurate surgery for lesions within the lateral ventricles.

Individual Assessment of Brain Tissue Changes in MS and the Effect of Focal Lesions on Short-Term Focal Atrophy Development in MS: A Voxel-Guided Morphometry Study

We performed voxel-guided morphometry (VGM) investigating the mechanisms of brain atrophy in multiple sclerosis (MS) related to focal lesions. VGM maps detect regional brain changes when comparing 2 time points on high resolution T1-weighted (T1w) magnetic resonace imaging (MRI). Two T1w MR datasets from 92 relapsing-remitting MS patients obtained 12 months apart were analysed with VGM. New lesions and volume changes of focal MS lesions as well as in the surrounding tissue were identified by visual inspection on colour coded VGM maps. Lesions were dichotomized in active and inactive lesions. Active lesions, defined by either new lesions (NL) (volume increase > 5% in VGM), chronic enlarging lesions (CEL) (pre-existent T1w lesions with volume increase > 5%), or chronic shrinking lesions (CSL) (pre-existent T1w lesions with volume reduction > 5%) in VGM, were accompanied by tissue shrinkage in surrounding and/or functionally related regions. Volume loss within the corpus callosum was highly correlated with the number of lesions in its close proximity. Volume loss in the lateral geniculate nucleus was correlated with lesions along the optic radiation. VGM analysis provides strong evidence that all active lesion types (NL, CEL, and CSL) contribute to brain volume reduction in the vicinity of lesions and/or in anatomically and functionally related areas of the brain.

Ensemble Tractography

Tractography uses diffusion MRI to estimate the trajectory and cortical projection zones of white matter fascicles in the living human brain. There are many different tractography algorithms and each requires the user to set several parameters, such as curvature threshold. Choosing a single algorithm with specific parameters poses two challenges. First, different algorithms and parameter values produce different results. Second, the optimal choice of algorithm and parameter value may differ between different white matter regions or different fascicles, subjects, and acquisition parameters. We propose using ensemble methods to reduce algorithm and parameter dependencies. To do so we separate the processes of fascicle generation and evaluation. Specifically, we analyze the value of creating optimized connectomes by systematically combining candidate streamlines from an ensemble of algorithms (deterministic and probabilistic) and systematically varying parameters (curvature and stopping criterion). The ensemble approach leads to optimized connectomes that provide better cross-validated prediction error of the diffusion MRI data than optimized connectomes generated using a single-algorithm or parameter set. Furthermore, the ensemble approach produces connectomes that contain both short- and long-range fascicles, whereas single-parameter connectomes are biased towards one or the other. In summary, a systematic ensemble tractography approach can produce connectomes that are superior to standard single parameter estimates both for predicting the diffusion measurements and estimating white matter fascicles.

Diffusion MRI and tractography opened a new avenue for studying white matter fascicles and their tissue properties in the living human brain. There are many different tractography methods, and each requires the user to set several parameters. A limitation of tractography is that the results depend on the selection of algorithms and parameters. Here, we analyze an ensemble method, Ensemble Tractography (ET), that reduces the effect of algorithm and parameter selection. ET creates a large set of candidate streamlines using an ensemble of algorithms and parameter values and then selects the streamlines with strong support from the data using a global fascicle evaluation method. Compared to single parameter connectomes, ET connectomes predict diffusion MRI signals better and cover a wider range of white matter volume. Importantly, ET connectomes include both short- and long-association fascicles, which are not typically found together in single-parameter connectomes.

Automated retinofugal visual pathway reconstruction with multi-shell HARDI and FOD-based analysis

Diffusion MRI tractography provides a non-invasive modality to examine the human retinofugal projection, which consists of the optic nerves, optic chiasm, optic tracts, the lateral geniculate nuclei (LGN) and the optic radiations. However, the pathway has several anatomic features that make it particularly challenging to study with tractography, including its location near blood vessels and bone-air interface at the base of the cerebrum, crossing fibers at the chiasm, somewhat-tortuous course around the temporal horn via Meyer's Loop, and multiple closely neighboring fiber bundles. To date, these unique complexities of the visual pathway have impeded the development of a robust and automated reconstruction method using tractography. To overcome these challenges, we develop a novel, fully automated system to reconstruct the retinofugal visual pathway from high-resolution diffusion imaging data. Using multi-shell, high angular resolution diffusion imaging (HARDI) data, we reconstruct precise fiber orientation distributions (FODs) with high order spherical harmonics (SPHARM) to resolve fiber crossings, which allows the tractography algorithm to successfully navigate the complicated anatomy surrounding the retinofugal pathway. We also develop automated algorithms for the identification of ROIs used for fiber bundle reconstruction. In particular, we develop a novel approach to extract the LGN region of interest (ROI) based on intrinsic shape analysis of a fiber bundle computed from a seed region at the optic chiasm to a target at the primary visual cortex. By combining automatically identified ROIs and FOD-based tractography, we obtain a fully automated system to compute the main components of the retinofugal pathway, including the optic tract and the optic radiation. We apply our method to the multi-shell HARDI data of 215 subjects from the Human Connectome Project (HCP). Through comparisons with post-mortem dissection measurements, we demonstrate the retinotopic organization of the optic radiation including a successful reconstruction of Meyer's loop. Then, using the reconstructed optic radiation bundle from the HCP cohort, we construct a probabilistic atlas and demonstrate its consistency with a post-mortem atlas. Finally, we generate a shape-based representation of the optic radiation for morphometry analysis.

Improved Framework for Tractography Reconstruction of the Optic Radiation

The optic radiation (OR) is one of the major components of the visual system and a key structure at risk in white matter diseases such as multiple sclerosis (MS). However, it is challenging to perform track reconstruction of the OR using diffusion MRI due to a sharp change of direction in the Meyer's loop and the presence of kissing and crossing fibers along the pathway. As such, we aimed to provide a highly precise and reproducible framework for tracking the OR from thalamic and visual cortex masks. The framework combined the generation of probabilistic streamlines by high order fiber orientation distributions estimated with constrained spherical deconvolution and an automatic post-processing based on anatomical exclusion criteria (AEC) to compensate for the presence of anatomically implausible streamlines. Specifically, those ending in the contralateral hemisphere, cerebrospinal fluid or grey matter outside the visual cortex were automatically excluded. We applied the framework to two distinct high angular resolution diffusion-weighted imaging (HARDI) acquisition protocols on one cohort, comprised of ten healthy volunteers and five MS patients. The OR was successfully delineated in both HARDI acquisitions in the healthy volunteers and MS patients. Quantitative evaluation of the OR position was done by comparing the results with histological reference data. Compared with histological mask, the OR reconstruction into a template (OR-TCT) was highly precise (percentage of voxels within the OR-TCT correctly defined as OR), ranging from 0.71 to 0.83. The sensitivity (percentage of voxels in histological reference mask correctly defined as OR in OR-TCT) ranged from 0.65 to 0.81 and the accuracy (measured by F1 score) was 0.73 to 0.77 in healthy volunteers. When AEC was not applied the precision and accuracy decreased. The absolute agreement between both HARDI datasets measured by the intraclass correlation coefficient was 0.73. This improved framework allowed us to reconstruct the OR with high reliability and accuracy independently of the acquisition parameters. Moreover, the reconstruction was possible even in the presence of tissue damage due to MS. This framework could also be applied to other tracts with complex configuration.

Tractography of the optic radiation: a repeatability and reproducibility study

Our main objective was to evaluate the repeatability and reproducibility of optic radiation (OR) reconstruction from diffusion MRI (dMRI) data. 14 adults were scanned twice with the same 60-direction dMRI sequence. Peaks in the diffusion profile were estimated with the single tensor (ST), Q-ball (QSH) and persistent angular structure (PAS) methods. Segmentation of the OR was performed by two experimenters with probabilistic tractography based on a manually drawn region-of-interest (ROI) protocol typically employed for OR segmentation, with both standard and extended sets of ROIs. The repeatability and reproducibility were assessed by calculating the intra-class correlation coefficient (ICC) of intra- and inter-rater experiments, respectively. ICCs were calculated for commonly used dMRI metrics (FA, MD, AD, RD) and anatomical dimensions of the optic radiation (distance from Meyer's loop to the temporal pole, ML-TP), as well as the Dice similarity coefficient (DSC) between the raters' OR segmentation. Bland-Altman plots were also calculated to investigate bias and variability in the reproducibility measurements. The OR was successfully reconstructed in all subjects by both raters. The ICC was found to be in the good to excellent range for both repeatability and reproducibility of the dMRI metrics, DSC and ML-TP distance. The Bland-Altman plots did not show any apparent systematic bias for any quantities. Overall, higher ICC values were found for the multi-fiber methods, QSH and PAS, and for the standard set of ROIs. Considering the good to excellent repeatability and reproducibility of all the quantities investigated, these findings support the use of multi-fiber OR reconstruction with a limited number of manually drawn ROIs in clinical applications utilizing either OR microstructure characterization or OR dimensions, as is the case in neurosurgical planning for temporal lobectomy.

Probabilistic MRI tractography of the optic radiation using constrained spherical deconvolution: a feasibility study

Background and Purpose

Imaging the optic radiation (OR) is of considerable interest in studying diseases affecting the visual pathway and for pre-surgical planning of temporal lobe resections. The purpose of this study was to investigate the clinical feasibility of using probabilistic diffusion tractography based on constrained spherical deconvolution (CSD) to image the optic radiation. It was hypothesized that CSD would provide improved tracking of the OR compared with the widely used ball-and-stick model.

Methods

Diffusion weighted MRI (30 directions) was performed on twenty patients with no known visual deficits. Tractography was performed using probabilistic algorithms based on fiber orientation distribution models of local white matter trajectories. The performance of these algorithms was evaluated by comparing computational times and receiver operating characteristic results, and by correlation of anatomical landmark distances to dissection estimates.

Results

The results showed that it was consistently feasible to reconstruct individual optic radiations from clinically practical (4.5 minute acquisition) diffusion weighted imaging data sets using CSD. Tractography based on the CSD model resulted in significantly shorter computational times, improved receiver operating characteristic results, and shorter Meyer’s loop to temporal pole distances (in closer agreement with dissection studies) when compared to the ball-and-stick based algorithm.

Conclusions

Accurate tractography of the optic radiation can be accomplished using diffusion MRI data collected within a clinically practical timeframe. CSD based tractography was faster, more accurate and had better correlation with known anatomical landmarks than ball-and-stick tractography.

Using fMRI non-local means denoising to uncover activation in sub-cortical structures at 1.5 T for guided HARDI tractography

In recent years, there has been ever-increasing interest in combining functional magnetic resonance imaging (fMRI) and diffusion magnetic resonance imaging (dMRI) for better understanding the link between cortical activity and connectivity, respectively. However, it is challenging to detect and validate fMRI activity in key sub-cortical areas such as the thalamus, given that they are prone to susceptibility artifacts due to the partial volume effects (PVE) of surrounding tissues (GM/WM interface). This is especially true on relatively low-field clinical MR systems (e.g., 1.5 T). We propose to overcome this limitation by using a spatial denoising technique used in structural MRI and more recently in diffusion MRI called non-local means (NLM) denoising, which uses a patch-based approach to suppress the noise locally. To test this, we measured fMRI in 20 healthy subjects performing three block-based tasks : eyes-open closed (EOC) and left/right finger tapping (FTL, FTR). Overall, we found that NLM yielded more thalamic activity compared to traditional denoising methods. In order to validate our pipeline, we also investigated known structural connectivity going through the thalamus using HARDI tractography: the optic radiations, related to the EOC task, and the cortico-spinal tract (CST) for FTL and FTR. To do so, we reconstructed the tracts using functionally based thalamic and cortical ROIs to initiates seeds of tractography in a two-level coarse-to-fine fashion. We applied this method at the single subject level, which allowed us to see the structural connections underlying fMRI thalamic activity. In summary, we propose a new fMRI processing pipeline which uses a recent spatial denoising technique (NLM) to successfully detect sub-cortical activity which was validated using an advanced dMRI seeding strategy in single subjects at 1.5 T.

Application of diffusion tensor imaging and tractography of the optic radiation in anterior temporal lobe resection for epilepsy: a systematic review

BACKGROUND

Approximately 50-100% of patients with temporal lobe epilepsy undergoing anterior temporal lobe resection (ATLR) will suffer a postoperative visual field defect (VFD) due to disruption of the optic radiation (OpR).

OBJECTIVE

We conducted a systematic review of the literature to examine the role of DTI and tractography in ATLR and its potential in reducing the incidence of postoperative VFD.

METHODS

We conducted an electronic literature search using PubMed, Embase, Web of Science and BMJ case report databases. Eligibility for study inclusion was determined on abstract screening using the following criteria: the study must have been (1) an original investigation or case report in humans; (2) investigating the OpR with DTI in cases of ATLR in temporal lobe epilepsy; (3) investigating postoperative VFD. All forms of ATLR and ways of assessing VFD were included to reflect clinical practice.

RESULTS

13 studies (four case reports, eight prospective observational studies, one prospective comparative trial) were included in the review, 179 (mean±SD, 13.8±12.6; range, 1-48) subjects were investigated using DTI. The time of postoperative VFD measurement differed between the detected studies, ranging from two weeks to nine years following ATLR. A modest number of studies and insufficient statistical homogeneity precluded meta-analysis. However, DTI methods were consistently accurate at quantifying and predicting postoperative damage to the OpR. These methods revealed a correlation between the extent of OpR damage and the severity of postoperative VFD. The first and only trial with 15 subjects compared to 23 controls reported that using intraoperative tractography in ATLR significantly reduces the occurrence of postoperative VFD on comparison to conventional surgical planning.

CONCLUSIONS

DTI shows potential to be an effective method used in planning ATLR. Findings from a single modest sized study suggest that tractography may be employed as part of intraoperative navigation techniques in order to avoid injury to the OpR. Further research needs to be conducted to ensure the applicability and effectiveness of this technology before implementation in routine clinical practice.

Neural correlates of impaired vision in adolescents born extremely preterm and/or extremely low birthweight

BACKGROUND

Adolescents born extremely preterm (EP; <28 weeks' gestation) and/or extremely low birthweight (ELBW; <1000 g) experience high rates of visual impairments, however the potential neural correlates of visual impairments in EP/ELBW adolescents require further investigation. This study aimed to: 1) compare optic radiation and primary visual cortical structure between EP/ELBW adolescents and normal birthweight controls; 2) investigate associations between perinatal factors and optic radiation and primary visual cortical structure in EP/ELBW adolescents; 3) investigate associations between optic radiation and primary visual cortical structure in EP/ELBW adolescents and the odds of impaired vision.

METHODS

196 EP/ELBW adolescents and 143 controls underwent magnetic resonance imaging at a mean age of 18 years. Optic radiations were delineated using constrained spherical deconvolution based probabilistic tractography. Primary visual cortices were segmented using FreeSurfer software. Diffusion tensor variables and tract volume of the optic radiations, as well as volume, surface area and thickness of the primary visual cortices, were estimated.

RESULTS

Axial, radial and mean diffusivities within the optic radiations, and primary visual cortical thickness, were higher in the EP/ELBW adolescents than controls. Within EP/ELBW adolescents, postnatal corticosteroid exposure was associated with altered optic radiation diffusion values and lower tract volume, while decreasing gestational age at birth was associated with increased primary visual cortical volume, area and thickness. Furthermore, decreasing optic radiation fractional anisotropy and tract volume, and increasing optic radiation diffusivity in EP/ELBW adolescents were associated with increased odds of impaired vision, whereas primary visual cortical measures were not associated with the odds of impaired vision.

CONCLUSIONS

Optic radiation and primary visual cortical structure are altered in EP/ELBW adolescents compared with controls, with the greatest alterations seen in those exposed to postnatal corticosteroids and those born earliest. Structural alterations to the optic radiations may increase the risk of impaired vision in EP/ELBW adolescents.

A mathematical framework for the registration and analysis of multi-fascicle models for population studies of the brain microstructure

Diffusion tensor imaging (DTI) is unable to represent the diffusion signal arising from multiple crossing fascicles and freely diffusing water molecules. Generative models of the diffusion signal, such as multi-fascicle models, overcome this limitation by providing a parametric representation for the signal contribution of each population of water molecules. These models are of great interest in population studies to characterize and compare the brain microstructural properties. Central to population studies is the construction of an atlas and the registration of all subjects to it. However, the appropriate definition of registration and atlasing methods for multi-fascicle models have proven challenging. This paper proposes a mathematical framework to register and analyze multi-fascicle models. Specifically, we define novel operators to achieve interpolation, smoothing and averaging of multi-fascicle models. We also define a novel similarity metric to spatially align multi-fascicle models. Our framework enables simultaneous comparisons of different microstructural properties that are confounded in conventional DTI. The framework is validated on multi-fascicle models from 24 healthy subjects and 38 patients with tuberous sclerosis complex, 10 of whom have autism. We demonstrate the use of the multi-fascicle models registration and analysis framework in a population study of autism spectrum disorder.

Epilepsy surgery, vision, and driving: what has surgery taught us and could modern imaging reduce the risk of visual deficits?

Up to 40% of patients with temporal lobe epilepsy (TLE) are refractory to medication. Surgery is an effective treatment but may cause new neurologic deficits including visual field deficits (VFDs). The ability to drive after surgery is a key goal, but a postoperative VFD precludes driving in 4-50% of patients even if seizure-free. VFDs are a consequence of damage to the most anterior portion of the optic radiation, Meyer's loop. Anatomic dissection reveals that the anterior extent of Meyer's loop is highly variable and may clothe the temporal horn, a key landmark entered during temporal lobe epilepsy surgery. Experience from surgery since the 1940s has shown that VFDs are common (48-100%) and that the degree of resection affects the frequency or severity of the deficit. The pseudowedge shape of the deficit has led to a revised retinotopic model of the organization of the optic radiation. Evidence suggests that the left optic radiation is more anterior and thus at greater risk. Alternative surgical approaches, such as selective amygdalo-hippocampectomy, may reduce this risk, but evidence is conflicting or lacking. The optic radiation can be delineated in vivo using diffusion tensor imaging tractography, which has been shown to be useful in predicting the postoperative VFDs and in surgical planning. These data are now being used for surgical guidance with the aim of reducing the severity of VFDs. Compensation for brain shift occurring during surgery can be performed using intraoperative magnetic resonance imaging (MRI), but the additional utility of this expensive technique remains unproven.