Creation of a novel trigeminal tractography atlas for automated trigeminal nerve identification

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) Study® has collected data from over 10,000 children across 21 sites, providing insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a dataset of harmonized and processed ABCD dMRI data (from release 3) has been created, comprising quality-controlled imaging data from 9,345 subjects, focusing exclusively on the baseline session, i.e., the first time point of the study. This resource required substantial computational time (approx. 50,000 CPU hours) for harmonization, whole-brain tractography, and white matter parcellation. The dataset includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts in full and low resolution, and 804 different dMRI-derived measures per subject (72.3 TB total size). Accessible via the NIMH Data Archive, it offers a large-scale dMRI dataset for studying structural connectivity in child and adolescent neurodevelopment. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Tractography-based automated identification of the retinogeniculate visual pathway with novel microstructure-informed supervised contrastive learning

The retinogeniculate visual pathway (RGVP) is responsible for carrying visual information from the retina to the lateral geniculate nucleus. Identification and visualization of the RGVP are important in studying the anatomy of the visual system and can inform the treatment of related brain diseases. Diffusion MRI (dMRI) tractography is an advanced imaging method that uniquely enables in vivo mapping of the 3D trajectory of the RGVP. Currently, identification of the RGVP from tractography data relies on expert (manual) selection of tractography streamlines, which is time-consuming, has high clinical and expert labor costs, and is affected by inter-observer variability. In this paper, we present a novel deep learning framework, DeepRGVP , to enable fast and accurate identification of the RGVP from dMRI tractography data. We design a novel microstructure-informed supervised contrastive learning method that leverages both streamline label and tissue microstructure information to determine positive and negative pairs. We propose a simple and successful streamline-level data augmentation method to address highly imbalanced training data, where the number of RGVP streamlines is much lower than that of non-RGVP streamlines. We perform comparisons with several state-of-the-art deep learning methods that were designed for tractography parcellation, and we show superior RGVP identification results using DeepRGVP. In addition, we demonstrate a good generalizability of DeepRGVP to dMRI tractography data from neurosurgical patients with pituitary tumors and we show DeepRGVP can successfully identify RGVPs despite the effect of lesions affecting the RGVPs. Overall, our study shows the high potential of using deep learning to automatically identify the RGVP.

Reconstructing the somatotopic organization of the corticospinal tract remains a challenge for modern tractography methods

The corticospinal tract (CST) is a critically important white matter fiber tract in the human brain that enables control of voluntary movements of the body. The CST exhibits a somatotopic organization, which means that the motor neurons that control specific body parts are arranged in order within the CST. Diffusion magnetic resonance imaging (MRI) tractography is increasingly used to study the anatomy of the CST. However, despite many advances in tractography algorithms over the past decade, modern, state-of-the-art methods still face challenges. In this study, we compare the performance of six widely used tractography methods for reconstructing the CST and its somatotopic organization. These methods include constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, unscented Kalman filter (UKF) tractography methods including multi-fiber (UKF2T) and single-fiber (UKF1T) models, the generalized q-sampling imaging (GQI) based deterministic tractography method, and the TractSeg method. We investigate CST somatotopy by dividing the CST into four subdivisions per hemisphere that originate in the leg, trunk, hand, and face areas of the primary motor cortex. A quantitative and visual comparison is performed using diffusion MRI data (N = 100 subjects) from the Human Connectome Project. Quantitative evaluations include the reconstruction rate of the eight anatomical subdivisions, the percentage of streamlines in each subdivision, and the coverage of the white matter-gray matter (WM-GM) interface. CST somatotopy is further evaluated by comparing the percentage of streamlines in each subdivision to the cortical volumes for the leg, trunk, hand, and face areas. Overall, UKF2T has the highest reconstruction rate and cortical coverage. It is the only method with a significant positive correlation between the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex. However, our experimental results show that all compared tractography methods are biased toward generating many trunk streamlines (ranging from 35.10% to 71.66% of total streamlines across methods). Furthermore, the coverage of the WM-GM interface in the largest motor area (face) is generally low (under 40%) for all compared tractography methods. Different tractography methods give conflicting results regarding the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex, indicating that there is generally no clear relationship, and that reconstruction of CST somatotopy is still a large challenge. Overall, we conclude that while current tractography methods have made progress toward the well-known challenge of improving the reconstruction of the lateral projections of the CST, the overall problem of performing a comprehensive CST reconstruction, including clinically important projections in the lateral (hand and face areas) and medial portions (leg area), remains an important challenge for diffusion MRI tractography.

A unified global tractography framework for automatic visual pathway reconstruction

The human visual pathway starts from the retina, passes through the retinogeniculate visual pathway, the optic radiation, and finally connects to the primary visual cortex. Diffusion MRI tractography is the only technology that can noninvasively reconstruct the visual pathway. However, complete and accurate visual pathway reconstruction is challenging because of the skull base environment and complex fiber geometries. Specifically, the optic nerve within the complex skull base environment can cause abnormal diffusion signals. The crossing and fanning fibers at the optic chiasm, and a sharp turn of Meyer's loop at the optic radiation, contribute to complex fiber geometries of the visual pathway. A fiber trajectory distribution (FTD) function-based tractography method of our previous work and several high sensitivity tractography methods can reveal these complex fiber geometries, but are accompanied by false-positive fibers. Thus, the related studies of the visual pathway mostly applied the expert region of interest selection strategy. However, interobserver variability is an issue in reconstructing an accurate visual pathway. In this paper, we propose a unified global tractography framework to automatically reconstruct the visual pathway. We first extend the FTD function to a high-order streamline differential equation for global trajectory estimation. At the global level, the tractography process is simplified as the estimation of global trajectory distribution coefficients by minimizing the cost between trajectory distribution and the selected directions under the prior guidance by introducing the tractography template as anatomic priors. Furthermore, we use a deep learning-based method and tractography template prior information to automatically generate the mask for tractography. The experimental results demonstrate that our proposed method can successfully reconstruct the visual pathway with high accuracy.

Probabilistic tractography of the extracranial branches of the trigeminal nerve using diffusion tensor imaging

PURPOSE

The peripheral course of the trigeminal nerves is complex and spans multiple bony foramen and tissue compartments throughout the face. Diffusion tensor imaging of these nerves is difficult due to the complex tissue interfaces and relatively low MR signal. The purpose of this work is to develop a method for reliable diffusion tensor imaging-based fiber tracking of the peripheral branches of the trigeminal nerve.

METHODS

We prospectively acquired imaging data from six healthy adult participants with a 3.0-Tesla system, including T2-weighted short tau inversion recovery with variable flip angle (T2-STIR-SPACE) and readout segmented echo planar diffusion weighted imaging sequences. Probabilistic tractography of the ophthalmic, infraorbital, lingual, and inferior alveolar nerves was performed manually and assessed by two observers who determined whether the fiber tracts reached defined anatomical landmarks using the T2-STIR-SPACE volume.

RESULTS

All nerves in all subjects were tracked beyond the trigeminal ganglion. Tracts in the inferior alveolar and ophthalmic nerve exhibited the strongest signal and most consistently reached the most distal landmark (58% and 67%, respectively). All tracts of the inferior alveolar and ophthalmic nerve extended beyond their respective third benchmarks. Tracts of the infraorbital nerve and lingual nerve were comparably lower-signal and did not consistently reach the furthest benchmarks (9% and 17%, respectively).

CONCLUSION

This work demonstrates a method for consistently identifying and tracking the major nerve branches of the trigeminal nerve with diffusion tensor imaging.

Predictive value of preoperative magnetic resonance imaging structural and diffusion indices for the results of trigeminal neuralgia microvascular decompression surgery

PURPOSE

To explore the predictive value of preoperative magnetic resonance imaging structural and diffusion indices of the spinal trigeminal tract (SpTV) on the results of microvascular decompression (MVD) in patients with trigeminal neuralgia (TN).

METHODS

This retrospective study included patients diagnosed with TN and treated with MVD in the Jining First People's Hospital between January 2020 and January 2021. The patients were divided into good and poor results groups according to postoperative pain relief. Logistic regression analysis was performed to explore independent risk factors for poor results of MVD, and their predictive value was examined using receiver operating characteristic (ROC) curves.

RESULTS

A total of 97 TN cases were included, 24 cases with a poor result and 73 with a good result. They were comparable in demographic characteristics. Fractional anisotropy (FA) was lower (P < 0.001), and radial diffusivity (RD) was higher (P < 0.001) in the poor result group compared to the good result group. Patients in the good result group showed a higher proportion of grade 3 neurovascular contact (NVC) (39.7% vs. 16.7%, P = 0.001) and a lower RD (P < 0.001). The multivariate analysis showed that the RD of SpTV (OR = 0.000016, 95% CI: 0.000-0.004, P < 0.001) and NVC (OR = 8.07, 95% CI: 1.67-38.93, P = 0.009) were independently associated with poor results. The area under the curve (AUC) of RD and NVC were 0.848 and 0.710, and their combination achieved an AUC of 0.880.

CONCLUSION

NVC and RD of SpTV are independent risk factors for poor results after MVD surgery, and combining the NVC and RD might achieve relatively high predictive value for poor results.

CNTSeg: A multimodal deep-learning-based network for cranial nerves tract segmentation

The segmentation of cranial nerves (CNs) tracts based on diffusion magnetic resonance imaging (dMRI) provides a valuable quantitative tool for the analysis of the morphology and course of individual CNs. Tractography-based approaches can describe and analyze the anatomical area of CNs by selecting the reference streamlines in combination with ROIs-based (regions-of-interests) or clustering-based. However, due to the slender structure of CNs and the complex anatomical environment, single-modality data based on dMRI cannot provide a complete and accurate description, resulting in low accuracy or even failure of current algorithms in performing individualized CNs segmentation. In this work, we propose a novel multimodal deep-learning-based multi-class network for automated cranial nerves tract segmentation without using tractography, ROI placement or clustering, called CNTSeg. Specifically, we introduced T1w images, fractional anisotropy (FA) images, and fiber orientation distribution function (fODF) peaks into the training data set, and design the back-end fusion module which uses the complementary information of the interphase feature fusion to improve the segmentation performance. CNTSeg has achieved the segmentation of 5 pairs of CNs (i.e. optic nerve CN II, oculomotor nerve CN III, trigeminal nerve CN V, and facial-vestibulocochlear nerve CN VII/VIII). Extensive comparisons and ablation experiments show promising results and are anatomically convincing even for difficult tracts. The code will be openly available at https://github.com/IPIS-XieLei/CNTSeg.

Vestibular paroxysmia entails vestibular nerve function, microstructure and endolymphatic space changes linked to root-entry zone neurovascular compression

Combining magnetic resonance imaging (MRI) sequences that permit the determination of vestibular nerve angulation (NA = change of nerve caliber or direction), structural nerve integrity via diffusion tensor imaging (DTI), and exclusion of endolymphatic hydrops (ELH) via delayed gadolinium-enhanced MRI of the inner ear (iMRI) could increase the diagnostic accuracy in patients with vestibular paroxysmia (VP). Thirty-six participants were examined, 18 with VP (52.6 ± 18.1 years) and 18 age-matched with normal vestibulocochlear testing (NP 50.3 ± 16.5 years). This study investigated whether (i) NA, (ii) DTI changes, or (iii) ELH occur in VP, and (iv) to what extent said parameters relate. Methods included vestibulocochlear testing and MRI data analyses for neurovascular compression (NVC) and NA verification, DTI and ELS quantification. As a result, (i) NA increased NVC specificity. (ii) DTI structural integrity was reduced on the side affected by VP (p < 0.05). (iii) 61.1% VP showed mild ELH and higher asymmetry indices than NP (p > 0.05). (iv) "Disease duration" and "total number of attacks" correlated with the decreased structural integrity of the affected nerve in DTI (p < 0.001). NVC distance within the nerve's root-entry zone correlated with nerve function (Roh = 0.72, p < 0.001), nerve integrity loss (Roh = - 0.638, p < 0.001), and ELS volume (Roh = - 0.604, p < 0.001) in VP. In conclusion, this study is the first to link eighth cranial nerve function, microstructure, and ELS changes in VP to clinical features and increased vulnerability of NVC in the root-entry zone. Combined MRI with NVC or NA verification, DTI and ELS quantification increased the diagnostic accuracy at group-level but did not suffice to diagnose VP on a single-subject level due to individual variability and lack of diagnostic specificity.

Case Report: The Imperfect Association Between Craniofacial Lesion Burden and Pain in Fibrous Dysplasia

Patients with fibrous dysplasia (FD) often present with craniofacial lesions that affect the trigeminal nerve system. Debilitating pain, headache, and migraine are frequently experienced by FD patients with poor prognosis, while some individuals with similar bone lesions are asymptomatic. The clinical and biological factors that contribute to the etiopathogenesis of pain in craniofacial FD are largely unknown. We present two adult females with comparable craniofacial FD lesion size and location, as measured by ¹⁸F-sodium fluoride positron emission tomography/computed tomography (PET/CT), yet their respective pain phenotypes differed significantly. Over 4 weeks, the average pain reported by Patient A was 0.4/0–10 scale. Patient B reported average pain of 7.8/0–10 scale distributed across the entire skull and left facial region. Patient B did not experience pain relief from analgesics or more aggressive treatments (denosumab). In both patients, evaluation of trigeminal nerve divisions (V1, V2, and V3) with CT and magnetic resonance imaging (MRI) revealed nerve compression and displacement with more involvement of the left trigeminal branches relative to the right. First-time employment of diffusion MRI and tractography suggested reduced apparent fiber density within the cisternal segment of the trigeminal nerve, particularly for Patient B and in the left hemisphere. These cases highlight heterogeneous clinical presentation and neurobiological properties in craniofacial FD and also, the disconnect between peripheral pathology and pain severity. We hypothesize that a detailed phenotypic characterization of patients that incorporates an advanced imaging approach probing the trigeminal system may provide enhanced insights into the variable experiences with pain in craniofacial FD.

Automatic oculomotor nerve identification based on data‐driven fiber clustering

The oculomotor nerve (OCN) is the main motor nerve innervating eye muscles and can be involved in multiple flammatory, compressive, or pathologies. The diffusion magnetic resonance imaging (dMRI) tractography is now widely used to describe the trajectory of the OCN. However, the complex cranial structure leads to difficulties in fiber orientation distribution (FOD) modeling, fiber tracking, and region of interest (ROI) selection. Currently, the identification of OCN relies on expert manual operation, resulting in challenges, such as the carries high clinical, time‐consuming, and labor costs. Thus, we propose a method that can automatically identify OCN from dMRI tractography. First, we choose the multi‐shell multi‐tissue constraint spherical deconvolution (MSMT‐CSD) FOD estimation model and deterministic tractography to describe the 3D trajectory of the OCN. Then, we rely on the well‐established computational pipeline and anatomical expertise to create a data‐driven OCN tractography atlas from 40 HCP data. We identify six clusters belonging to the OCN from the atlas, including the structures of three kinds of positional relationships (pass between, pass through, and go around) with the red nuclei and two kinds of positional relationships with medial longitudinal fasciculus. Finally, we apply the proposed OCN atlas to identify the OCN automatically from 40 new HCP subjects and two patients with brainstem cavernous malformation. In terms of spatial overlap and visualization, experiment results show that the automatically and manually identified OCN fibers are consistent. Our proposed OCN atlas provides an effective tool for identifying OCN by avoiding the traditional selection strategy of ROIs.

Tractography of the ansa lenticularis in the human brain

The aim of this study was to make a thorough investigation of the trajectory of the ansa lenticularis (AL) and its subcomponents using high-resolution fiber-tracking tractography. The subcomponents of the AL were reconstructed from one region of interest (ROI) in the area of the globus pallidus combined with another ROI in the red nucleus, substantia nigra, subthalamic nucleus, or thalamus. This fiber-tracking protocol was tested in an HCP-1065 template, 35 healthy subjects from Massachusetts General Hospital (MGH), and 20 healthy subjects from the human connectome project (HCP) using generalized q-sampling imaging (GQI)-based tractography. Quantitative anisotropy and fractional anisotropy were also computed for the AL subcomponents. The subcomponents of the AL could be reconstructed in the HCP-1065 template, 35 MGH healthy subjects, and 20 HCP healthy subjects. The AL descends from the globus pallidus and joins the ansa peduncularis for a short distance, subdividing later into fibers that continue separately to the red nucleus, substantia nigra, subthalamic nucleus, and thalamus. The study demonstrated the trajectory of the ansa lenticularis and its subcomponents using GQI-based tractography, improving our understanding of the anatomical connectivity between the globus pallidus and the thalamo-subthalamic region in the human brain. One Sentence Summary The investigation of the ansa lenticularis and its subcomponents using high-resolution diffusion images based tractography.

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.

Deep white matter analysis (DeepWMA): Fast and consistent tractography segmentation

White matter tract segmentation, i.e. identifying tractography fibers (streamline trajectories) belonging to anatomically meaningful fiber tracts, is an essential step to enable tract quantification and visualization. In this study, we present a deep learning tractography segmentation method (DeepWMA) that allows fast and consistent identification of 54 major deep white matter fiber tracts from the whole brain. We create a large-scale training tractography dataset of 1 million labeled fiber samples, and we propose a novel 2D multi-channel feature descriptor (FiberMap) that encodes spatial coordinates of points along each fiber. We learn a convolutional neural network (CNN) fiber classification model based on FiberMap and obtain a high fiber classification accuracy of 90.99% on the training tractography data with ground truth fiber labels. Then, the method is evaluated on a test dataset of 597 diffusion MRI scans from six independently acquired populations across genders, the lifespan (1 day - 82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). We perform comparisons with two state-of-the-art tract segmentation methods. Experimental results show that our method obtains a highly consistent tract segmentation result, where on average over 99% of the fiber tracts are successfully identified across all subjects under study, most importantly, including neonates and patients with space-occupying brain tumors. We also demonstrate good generalization of the method to tractography data from multiple different fiber tracking methods. The proposed method leverages deep learning techniques and provides a fast and efficient tool for brain white matter segmentation in large diffusion MRI tractography datasets.