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.

Rapid motion estimation and correction using self-encoded FID navigators in 3D radial MRI

PURPOSE

To develop a self-navigated motion compensation strategy for 3D radial MRI that can compensate for continuous head motion by measuring rigid body motion parameters with high temporal resolution from the central k-space acquisition point (self-encoded FID navigator) in each radial spoke.

METHODS

A forward model was created from low-resolution calibration data to simulate the effect of relative motion between the coil sensitivity profiles and the underlying object on the self-encoded FID navigator signal. Trajectory deviations were included in the model as low spatial-order field variations. Three volunteers were imaged at 3 T using a modified 3D gradient-echo sequence acquired with a Kooshball trajectory while performing abrupt and continuous head motion. Rigid body-motion parameters were estimated from the central k-space signal of each spoke using a least-squares fitting algorithm. The accuracy of self-navigated motion parameters was assessed relative to an established external tracking system. Quantitative image quality metrics were computed for images with and without retrospective correction using external and self-navigated motion measurements.

RESULTS

Self-encoded FID navigators achieved mean absolute errors of 0.69 ± 0.82 mm and 0.73 ± 0.87° relative to external tracking for maximum motion amplitudes of 12 mm and 10°. Retrospective correction of the 3D radial data resulted in substantially improved image quality for both abrupt and continuous motion paradigms, comparable to external tracking results.

CONCLUSIONS

Accurate rigid body motion parameters can be rapidly obtained from self-encoded FID navigator signals in 3D radial MRI to continuously correct for head movements. This approach is suitable for robust neuroanatomical imaging in subjects that exhibit patterns of large and frequent motion.

Divergent growth of the transient brain compartments in fetuses with nonsyndromic isolated clefts involving the primary and secondary palate

Cleft lip/palate is a common orofacial malformation that often leads to speech/language difficulties as well as developmental delays in affected children, despite surgical repair. Our understanding of brain development in these children is limited. This study aimed to analyze prenatal brain development in fetuses with cleft lip/palate and controls. We examined in utero MRIs of 30 controls and 42 cleft lip/palate fetal cases and measured regional brain volumes. Cleft lip/palate was categorized into groups A (cleft lip or alveolus) and B (any combination of clefts involving the primary and secondary palates). Using a repeated-measures regression model with relative brain hemisphere volumes (%), and after adjusting for multiple comparisons, we did not identify significant differences in regional brain growth between group A and controls. Group B clefts had significantly slower weekly cerebellar growth compared with controls. We also observed divergent brain growth in transient brain structures (cortical plate, subplate, ganglionic eminence) within group B clefts, depending on severity (unilateral or bilateral) and defect location (hemisphere ipsilateral or contralateral to the defect). Further research is needed to explore the association between regional fetal brain growth and cleft lip/palate severity, with the potential to inform early neurodevelopmental biomarkers and personalized diagnostics.

IVIM-Morph: Motion-compensated quantitative Intra-voxel Incoherent Motion (IVIM) analysis for functional fetal lung maturity assessment from diffusion-weighted MRI data

Quantitative analysis of pseudo-diffusion in diffusion-weighted magnetic resonance imaging (DWI) data shows potential for assessing fetal lung maturation and generating valuable imaging biomarkers. Yet, the clinical utility of DWI data is hindered by unavoidable fetal motion during acquisition. We present IVIM-morph, a self-supervised deep neural network model for motion-corrected quantitative analysis of DWI data using the Intra-voxel Incoherent Motion (IVIM) model. IVIM-morph combines two sub-networks, a registration sub-network, and an IVIM model fitting sub-network, enabling simultaneous estimation of IVIM model parameters and motion. To promote physically plausible image registration, we introduce a biophysically informed loss function that effectively balances registration and model-fitting quality. We validated the efficacy of IVIM-morph by establishing a correlation between the predicted IVIM model parameters of the lung and gestational age (GA) using fetal DWI data of 39 subjects. Our approach was compared against six baseline methods: 1) no motion compensation, 2) affine registration of all DWI images to the initial image, 3) deformable registration of all DWI images to the initial image, 4) deformable registration of each DWI image to its preceding image in the sequence, 5) iterative deformable motion compensation combined with IVIM model parameter estimation, and 6) self-supervised deep-learning-based deformable registration. IVIM-morph exhibited a notably improved correlation with gestational age (GA) when performing in-vivo quantitative analysis of fetal lung DWI data during the canalicular phase. Specifically, over 2 test groups of cases, it achieved an R f 2 of 0.44 and 0.52, outperforming the values of 0.27 and 0.25, 0.25 and 0.00, 0.00 and 0.00, 0.38 and 0.00, and 0.07 and 0.14 obtained by other methods. IVIM-morph shows potential in developing valuable biomarkers for non-invasive assessment of fetal lung maturity with DWI data. Moreover, its adaptability opens the door to potential applications in other clinical contexts where motion compensation is essential for quantitative DWI analysis. The IVIM-morph code is readily available at:https://github.com/TechnionComputationalMRILab/qDWI-Morph.

Characterizing microstructural development in the fetal brain using diffusion MRI from 23 to 36 weeks of gestation

We utilized motion-corrected diffusion tensor imaging (DTI) to evaluate microstructural changes in healthy fetal brains during the late second and third trimesters. Data were derived from fetal magnetic resonance imaging scans conducted as part of a prospective study spanning from 2013 March to 2019 May. The study included 44 fetuses between the gestational ages (GAs) of 23 and 36 weeks. We reconstructed fetal brain DTI using a motion-tracked slice-to-volume registration framework. Images were segmented into 14 regions of interest (ROIs) through label propagation using a fetal DTI atlas, with expert refinement. Statistical analysis involved assessing changes in fractional anisotropy (FA) and mean diffusivity (MD) throughout gestation using mixed-effects models, and identifying points of change in trajectory for ROIs with nonlinear trends. Results showed significant GA-related changes in FA and MD in all ROIs except in the thalamus' FA and corpus callosum's MD. Hemispheric asymmetries were found in the FA of the periventricular white matter (pvWM), intermediate zone, and subplate and in the MD of the ganglionic eminence and pvWM. This study provides valuable insight into the normal patterns of development of MD and FA in the fetal brain. These changes are closely linked with cytoarchitectonic changes and display indications of early functional specialization.

Improved myelin water fraction mapping with deep neural networks using synthetically generated 3D data

We introduce a generative model for synthesis of large scale 3D datasets for quantitative parameter mapping of myelin water fraction (MWF). Our model combines a MR physics signal decay model with an accurate probabilistic multi-component parametric T2 model. We synthetically generate a wide variety of high quality signals and corresponding parameters from a wide range of naturally occurring prior parameter values. To capture spatial variation, the generative signal decay model is combined with a generative spatial model conditioned on generic tissue segmentations. Synthesized 3D datasets can be used to train any convolutional neural network (CNN) based architecture for MWF estimation. Our source code is available at: https://github.com/quin-med-harvard-edu/synthmap Reduction of acquisition time at the expense of lower SNR, as well as accuracy and repeatability of MWF estimation techniques, are key factors that affect the adoption of MWF mapping in clinical practice. We demonstrate that the synthetically trained CNN provides superior accuracy over the competing methods under the constraints of naturally occurring noise levels as well as on the synthetically generated images at low SNR levels. Normalized root mean squared error (nRMSE) is less than 7% on synthetic data, which is significantly lower than competing methods. Additionally, the proposed method yields a coefficient of variation (CoV) that is at least 4x better than the competing method on intra-session test-retest reference dataset.

Population Atlas Analysis of Emerging Brain Structural Connections in the Human Fetus

BACKGROUND

A lack of in utero imaging data hampers our understanding of the connections in the human fetal brain. Generalizing observations from postmortem subjects and premature newborns is inaccurate due to technical and biological differences.

PURPOSE

To evaluate changes in fetal brain structural connectivity between 23 and 35 weeks postconceptional age using a spatiotemporal atlas of diffusion tensor imaging (DTI).

STUDY TYPE

Retrospective.

POPULATION

Publicly available diffusion atlases, based on 60 healthy women (age 18-45 years) with normal prenatal care, from 23 and 35 weeks of gestation.

FIELD STRENGTH/SEQUENCE

3.0 Tesla/DTI acquired with diffusion-weighted echo planar imaging (EPI).

ASSESSMENT

We performed whole-brain fiber tractography from DTI images. The cortical plate of each diffusion atlas was segmented and parcellated into 78 regions derived from the Edinburgh Neonatal Atlas (ENA33). Connectivity matrices were computed, representing normalized fiber connections between nodes. We examined the relationship between global efficiency (GE), local efficiency (LE), small-worldness (SW), nodal efficiency (NE), and betweenness centrality (BC) with gestational age (GA) and with laterality.

STATISTICAL TESTS

Linear regression was used to analyze changes in GE, LE, NE, and BC throughout gestation, and to assess changes in laterality. The t-tests were used to assess SW. P-values were corrected using Holm-Bonferroni method. A corrected P-value <0.05 was considered statistically significant.

RESULTS

Network analysis revealed a significant weekly increase in GE (5.83%/week, 95% CI 4.32-7.37), LE (5.43%/week, 95% CI 3.63-7.25), and presence of SW across GA. No significant hemisphere differences were found in GE (P = 0.971) or LE (P = 0.458). Increasing GA was significantly associated with increasing NE in 41 nodes, increasing BC in 3 nodes, and decreasing BC in 2 nodes.

DATA CONCLUSION

Extensive network development and refinement occur in the second and third trimesters, marked by a rapid increase in global integration and local segregation.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 2.

Atypical fetal brain development in fetuses with non-syndromic isolated musculoskeletal birth defects (niMSBDs)

Non-syndromic, isolated musculoskeletal birth defects (niMSBDs) are among the leading causes of pediatric hospitalization. However, little is known about brain development in niMSBDs. Our study aimed to characterize prenatal brain development in fetuses with niMSBDs and identify altered brain regions compared to controls. We retrospectively analyzed in vivo structural T2-weighted MRIs of 99 fetuses (48 controls and 51 niMSBDs cases). For each group (19-31 and >31 gestational weeks (GW)), we conducted repeated-measures regression analysis with relative regional volume (% brain hemisphere) as a dependent variable (adjusted for age, side, and interactions). Between 19 and 31GW, fetuses with niMSBDs had a significantly (P < 0.001) smaller relative volume of the intermediate zone (-22.9 ± 3.2%) and cerebellum (-16.1 ± 3.5%,) and a larger relative volume of proliferative zones (38.3 ± 7.2%), the ganglionic eminence (34.8 ± 7.3%), and the ventricles (35.8 ± 8.0%). Between 32 and 37 GW, compared to the controls, niMSBDs showed significantly smaller volumes of central regions (-9.1 ± 2.1%) and larger volumes of the cortical plate. Our results suggest there is altered brain development in fetuses with niMSBDs compared to controls (13.1 ± 4.2%). Further basic and translational neuroscience research is needed to better visualize these differences and to characterize the altered development in fetuses with specific niMSBDs.

SlicerCBM: automatic framework for biomechanical analysis of the brain

PURPOSE

Brain shift that occurs during neurosurgery disturbs the brain's anatomy. Prediction of the brain shift is essential for accurate localisation of the surgical target. Biomechanical models have been envisaged as a possible tool for such predictions. In this study, we created a framework to automate the workflow for predicting intra-operative brain deformations.

METHODS

We created our framework by uniquely combining our meshless total Lagrangian explicit dynamics (MTLED) algorithm for computing soft tissue deformations, open-source software libraries and built-in functions within 3D Slicer, an open-source software package widely used for medical research. Our framework generates the biomechanical brain model from the pre-operative MRI, computes brain deformation using MTLED and outputs results in the form of predicted warped intra-operative MRI.

RESULTS

Our framework is used to solve three different neurosurgical brain shift scenarios: craniotomy, tumour resection and electrode placement. We evaluated our framework using nine patients. The average time to construct a patient-specific brain biomechanical model was 3 min, and that to compute deformations ranged from 13 to 23 min. We performed a qualitative evaluation by comparing our predicted intra-operative MRI with the actual intra-operative MRI. For quantitative evaluation, we computed Hausdorff distances between predicted and actual intra-operative ventricle surfaces. For patients with craniotomy and tumour resection, approximately 95% of the nodes on the ventricle surfaces are within two times the original in-plane resolution of the actual surface determined from the intra-operative MRI.

CONCLUSION

Our framework provides a broader application of existing solution methods not only in research but also in clinics. We successfully demonstrated the application of our framework by predicting intra-operative deformations in nine patients undergoing neurosurgical procedures.

Epilepsy Severity Is Associated With Head Circumference and Growth Rate in Infants With Tuberous Sclerosis Complex

BACKGROUND

Abnormal brain growth in tuberous sclerosis complex (TSC) reflects abnormalities in cellular proliferation and differentiation and results in epilepsy and other neurological manifestations. Head circumference (HC) as a proxy for brain volume may provide an easily tracked clinical measure of brain overgrowth and neurological disease burden. This study investigated the relationship between HC and epilepsy severity in infants with TSC.

METHODS

Prospective multicenter observational study of children from birth to three years with TSC. Epilepsy data were collected from clinical history, and HC was collected at study visits at age three, six, nine, 12, 18, 24, and 36 months. Epilepsy severity was classified as no epilepsy, low epilepsy severity (one seizure type and one or two antiepileptic drugs [AEDs]), moderate epilepsy severity (either two to three seizure types and one to two AEDs or one seizure type and more than three AEDs), or high epilepsy severity (two to three seizure types and more than three AEDs).

RESULTS

As a group, children with TSC had HCs approximately 1 S.D. above the mean World Health Organization (WHO) reference by age one year and demonstrated more rapid growth than the normal population reference. Males with epilepsy had larger HCs than those without. Compared with the WHO reference population, infants with TSC and no epilepsy or low or moderate epilepsy had an increased early HC growth rate, whereas those with severe epilepsy had an early larger HC but did not have a faster growth rate.

CONCLUSIONS

Infants and young children with TSC have larger HCs than typical growth norms and have differing rates of head growth depending on the severity of epilepsy.

Brain growth in fetuses with congenital diaphragmatic hernia

BACKGROUND AND PURPOSE

To perform a volumetric evaluation of the brain in fetuses with right or left congenital diaphragmatic hernia (CDH), and to compare brain growth trajectories to normal fetuses.

METHODS

We identified fetal MRIs performed between 2015 and 2020 in fetuses with a diagnosis of CDH. Gestational age (GA) range was 19-40 weeks. Control subjects consisted of normally developing fetuses between 19 and 40 weeks recruited for a separate prospective study. All images were acquired at 3 Tesla and were processed with retrospective motion correction and slice-to-volume reconstruction to generate super-resolution 3-dimensional volumes. These volumes were registered to a common atlas space and segmented in 29 anatomic parcellations.

RESULTS

A total of 174 fetal MRIs in 149 fetuses were analyzed (99 controls [mean GA: 29.2 ± 5.2 weeks], 34 fetuses left-sided CDH [mean GA: 28.4 ± 5.3 weeks], and 16 fetuses right-sided CDH [mean GA: 27 ± 5.4 weeks]). In fetuses with left-sided CDH, brain parenchymal volume was -8.0% (95% confidence interval [CI] [-13.1, -2.5]; p = .005) lower than normal controls. Differences ranged from -11.4% (95% CI [-18, -4.3]; p < .001) in the corpus callosum to -4.6% (95% CI [-8.9, -0.1]; p = .044) in the hippocampus. In fetuses with right-sided CDH, brain parenchymal volume was -10.1% (95% CI [-16.8, -2.7]; p = .008) lower than controls. Differences ranged from -14.1% (95% CI [-21, -6.5]; p < .001) in the ventricular zone to -5.6% (95% CI [-9.3, -1.8]; p = .025) in the brainstem.

CONCLUSION

Left and right CDH are associated with lower fetal brain volumes.

Heterogeneous growth of the insula shapes the human brain

The human cerebrum consists of a precise and stereotyped arrangement of lobes, primary gyri, and connectivity that underlies human cognition [P. Rakic, Nat. Rev. Neurosci. 10, 724-735 (2009)]. The development of this arrangement is less clear. Current models explain individual primary gyrification but largely do not account for the global configuration of the cerebral lobes [T. Tallinen, J. Y. Chung, J. S. Biggins, L. Mahadevan, Proc. Natl. Acad. Sci. U.S.A. 111, 12667-12672 (2014) and D. C. Van Essen, Nature 385, 313-318 (1997)]. The insula, buried in the depths of the Sylvian fissure, is unique in terms of gyral anatomy and size. Here, we quantitatively show that the insula has unique morphology and location in the cerebrum and that these key differences emerge during fetal development. Finally, we identify quantitative differences in developmental migration patterns to the insula that may underlie these differences. We calculated morphologic data in the insula and other lobes in adults (N = 107) and in an in utero fetal brain atlas (N = 81 healthy fetuses). In utero, the insula grows an order of magnitude slower than the other lobes and demonstrates shallower sulci, less curvature, and less surface complexity both in adults and progressively throughout fetal development. Spherical projection analysis demonstrates that the lenticular nuclei obstruct 60 to 70% of radial pathways from the ventricular zone (VZ) to the insula, forcing a curved migration to the insula in contrast to a direct radial pathway. Using fetal diffusion tractography, we identify radial glial fascicles that originate from the VZ and curve around the lenticular nuclei to form the insula. These results confirm existing models of radial migration to the cortex and illustrate findings that suggest differential insular and cerebral development, laying the groundwork to understand cerebral malformations and insular function and pathologies.

Comparison of fMRI language laterality with and without sedation in pediatric epilepsy

Functional MRI is an essential component of presurgical language mapping. In clinical settings, young children may be sedated for the MRI with the functional stimuli presented passively. Research has found that sedation changes language activation in healthy adults and children. However, there is limited research comparing sedated and unsedated functional MRI in pediatric epilepsy patients. We compared language activation patterns in children with epilepsy who received sedation for functional MRI to the ones who did not. We retrospectively identified the patients with focal epilepsy who underwent presurgical functional MRI including Auditory Descriptive Decision Task at Boston Children's Hospital from 2014 to 2022. Patients were divided into sedated and awake groups, based on their sedation status during functional MRI. Auditory Descriptive Decision Task stimuli were presented passively to the sedated group per clinical protocol. We extracted language activation maps contrasted against a control task (reverse speech) in the Frontal and Temporal language regions and calculated separate language laterality indexes for each region. We considered positive laterality indexes as left dominant, negative laterality indexes as right dominant, and absolute laterality indexes <0.2 as bilateral. We defined 2 language patterns: typical (i.e., primarily left-sided) and atypical. Typical pattern required at least one left dominant region (either frontal or temporal) and no right dominant region. We then compared the language patterns between the sedated and awake groups. Seventy patients met the inclusion criteria, 25 sedated, and 45 awake. Using the Auditory Descriptive Decision Task paradigm, when adjusted for age, handedness, gender, and laterality of lesion in a weighted logistic regression model, the odds of the atypical pattern were 13.2 times higher in the sedated group compared to the awake group (Confidence Interval: 2.55-68.41, p-value < 0.01). Sedation may alter language activation patterns in pediatric epilepsy patients. Language patterns on sedated functional MRI with passive tasks may not represent language networks during wakefulness, sedation may differentially suppress some networks, or require a different task or method of analysis to capture the awake language network. Given the critical surgical implication of these findings, additional studies are needed to better understand how sedation impacts the functional MRI blood oxygenation level-dependent signal. Consistent with current practice, sedated functional MRI should be interpreted with greater caution and requires additional validation as well as research on post-surgical language outcomes.

Image data and computational grids for computing brain shift and solving the electrocorticography forward problem

This article describes the dataset applied in the research reported in NeuroImage article "Patient-specific solution of the electrocorticography forward problem in deforming brain" [1] that is available for download from the Zenodo data repository (https://zenodo.org/record/7687631) [2]. Preoperative structural and diffusion-weighted magnetic resonance (MR) and postoperative computed tomography (CT) images of a 12-year-old female epilepsy patient under evaluation for surgical intervention were obtained retrospectively from Boston Children's Hospital. We used these images to conduct the analysis at The University of Western Australia's Intelligent Systems for Medicine Laboratory using SlicerCBM [3], our open-source software extension for the 3D Slicer medical imaging platform. As part of the analysis, we processed the images to extract the patient-specific brain geometry; created computational grids, including a tetrahedral grid for the meshless solution of the biomechanical model and a regular hexahedral grid for the finite element solution of the electrocorticography forward problem; predicted the postoperative MRI and DTI that correspond to the brain configuration deformed by the placement of subdural electrodes using biomechanics-based image warping; and solved the patient-specific electrocorticography forward problem to compute the electric potential distribution within the patient's head using the original preoperative and predicted postoperative image data. The well-established and open-source file formats used in this dataset, including Nearly Raw Raster Data (NRRD) files for images, STL files for surface geometry, and Visualization Toolkit (VTK) files for computational grids, allow other research groups to easily reuse the data presented herein to solve the electrocorticography forward problem accounting for the brain shift caused by implantation of subdural grid electrodes.

Morphological Features of Language Regions in Individuals with Tuberous Sclerosis Complex

A significant number of individuals with tuberous sclerosis complex (TSC) exhibit language difficulties. Here, we examined the language-related brain morphometry in 59 participants (7 participants with TSC and comorbid autism spectrum disorder (ASD) (TSC + ASD), 13 with TSC but no ASD (TSC-ASD), 10 with ASD-only (ASD), and 29 typically developing (TD) controls). A hemispheric asymmetry was noted in surface area and gray matter volume of several cortical language areas in TD, ASD, and TSC-ASD groups, but not in TSC + ASD group. TSC + ASD group demonstrated increased cortical thickness and curvature values in multiple language regions for both hemispheres, compared to other groups. After controlling for tuber load in the TSC groups, within-group differences stayed the same but the differences between TSC-ASD and TSC + ASD were no longer statistically significant. These preliminary findings suggest that comorbid ASD in TSC as well as tuber load in TSC is associated with changes in the morphometry of language regions. Future studies with larger sample sizes will be needed to confirm these findings.

Tubers Affecting the Fusiform Face Area Are Associated with Autism Diagnosis

OBJECTIVE

Tuberous sclerosis complex (TSC) is associated with focal brain "tubers" and a high incidence of autism spectrum disorder (ASD). The location of brain tubers associated with autism may provide insight into the neuroanatomical substrate of ASD symptoms.

METHODS

We delineated tuber locations for 115 TSC participants with ASD (n = 31) and without ASD (n = 84) from the Tuberous Sclerosis Complex Autism Center of Excellence Research Network. We tested for associations between ASD diagnosis and tuber burden within the whole brain, specific lobes, and at 8 regions of interest derived from the ASD neuroimaging literature, including the anterior cingulate, orbitofrontal and posterior parietal cortices, inferior frontal and fusiform gyri, superior temporal sulcus, amygdala, and supplemental motor area. Next, we performed an unbiased data-driven voxelwise lesion symptom mapping (VLSM) analysis. Finally, we calculated the risk of ASD associated with positive findings from the above analyses.

RESULTS

There were no significant ASD-related differences in tuber burden across the whole brain, within specific lobes, or within a priori regions derived from the ASD literature. However, using VLSM analysis, we found that tubers involving the right fusiform face area (FFA) were associated with a 3.7-fold increased risk of developing ASD.

INTERPRETATION

Although TSC is a rare cause of ASD, there is a strong association between tuber involvement of the right FFA and ASD diagnosis. This highlights a potentially causative mechanism for developing autism in TSC that may guide research into ASD symptoms more generally. ANN NEUROL 2023;93:577-590.

Detailed anatomic segmentations of a fetal brain diffusion tensor imaging atlas between 23 and 30 weeks of gestation

This work presents detailed anatomic labels for a spatiotemporal atlas of fetal brain Diffusion Tensor Imaging (DTI) between 23 and 30 weeks of post-conceptional age. Additionally, we examined developmental trajectories in fractional anisotropy (FA) and mean diffusivity (MD) across gestational ages (GA). We performed manual segmentations on a fetal brain DTI atlas. We labeled 14 regions of interest (ROIs): cortical plate (CP), subplate (SP), Intermediate zone-subventricular zone-ventricular zone (IZ/SVZ/VZ), Ganglionic Eminence (GE), anterior and posterior limbs of the internal capsule (ALIC, PLIC), genu (GCC), body (BCC), and splenium (SCC) of the corpus callosum (CC), hippocampus, lentiform Nucleus, thalamus, brainstem, and cerebellum. A series of linear regressions were used to assess GA as a predictor of FA and MD for each ROI. The combination of MD and FA allowed the identification of all ROIs. Increasing GA was significantly associated with decreasing FA in the CP, SP, IZ/SVZ/IZ, GE, ALIC, hippocampus, and BCC (p < .03, for all), and with increasing FA in the PLIC and SCC (p < .002, for both). Increasing GA was significantly associated with increasing MD in the CP, SP, IZ/SVZ/IZ, GE, ALIC, and CC (p < .03, for all). We developed a set of expert-annotated labels for a DTI spatiotemporal atlas of the fetal brain and presented a pilot analysis of developmental changes in cerebral microstructure between 23 and 30 weeks of GA.

Integrated response analysis of pediatric low-grade gliomas during and after targeted therapy treatment

Background

Pediatric low-grade gliomas (pLGGs) are the most common central nervous system tumor in children, characterized by RAS/MAPK pathway driver alterations. Genomic advances have facilitated the use of molecular targeted therapies, however, their long-term impact on tumor behavior remains critically unanswered.

Methods

We performed an IRB-approved, retrospective chart and imaging review of pLGGs treated with off-label targeted therapy at Dana-Farber/Boston Children's from 2010 to 2020. Response analysis was performed for BRAFV600E and BRAF fusion/duplication-driven pLGG subsets.

Results

Fifty-five patients were identified (dabrafenib n = 15, everolimus n = 26, trametinib n = 11, and vemurafenib n = 3). Median duration of targeted therapy was 9.48 months (0.12-58.44). The 1-year, 3-year, and 5-year EFS from targeted therapy initiation were 62.1%, 38.2%, and 31.8%, respectively. Mean volumetric change for BRAFV600E mutated pLGG on BRAF inhibitors was -54.11%; median time to best volumetric response was 8.28 months with 9 of 12 (75%) objective RAPNO responses. Median time to largest volume post-treatment was 2.86 months (+13.49%); mean volume by the last follow-up was -14.02%. Mean volumetric change for BRAF fusion/duplication pLGG on trametinib was +7.34%; median time to best volumetric response was 6.71 months with 3 of 7 (43%) objective RAPNO responses. Median time to largest volume post-treatment was 2.38 months (+71.86%); mean volume by the last follow-up was +39.41%.

Conclusions

Our integrated analysis suggests variability in response by pLGG molecular subgroup and targeted therapy, as well as the transience of some tumor growth following targeted therapy cessation.

Real-time shimming with FID navigators

PURPOSE

To implement a method for real-time field control using rapid FID navigator (FIDnav) measurements and evaluate the efficacy of the proposed approach for mitigating dynamic field perturbations and improvingT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality.

METHODS

FIDnavs were embedded in a gradient echo sequence and a subject-specific linear calibration model was generated on the scanner to facilitate rapid shim updates in response to measured FIDnav signals. To confirm the accuracy of FID-navigated field updates, phantom and volunteer scans were performed with online updates of the scanner B0 shim settings. To evaluate improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality with real-time shimming, 10 volunteers were scanned at 3T while performing deep-breathing and nose-touching tasks designed to modulate the B0 field. Quantitative image quality metrics were compared with and without FID-navigated field control. An additional volunteer was scanned at 7T to evaluate performance at ultra-high field.

RESULTS

Applying measured FIDnav shim updates successfully compensated for applied global and linear field offsets in phantoms and across all volunteers. FID-navigated real-time shimming led to a substantial reduction in field fluctuations and a consequent improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality in volunteers performing deep-breathing and nose-touching tasks, with 7.57% ± 6.01% and 8.21% ± 10.90% improvement in peak SNR and structural similarity, respectively.

CONCLUSION

FIDnavs facilitate rapid measurement and application of field coefficients for slice-wise B0 shimming. The proposed approach can successfully counteract spatiotemporal field perturbations and substantially improvesT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality, which is important for a variety of clinical and research applications, particularly at ultra-high field.

Diffusion-derived parameters in lesions, peri-lesion and normal-appearing white matter in multiple sclerosis using tensor, kurtosis and fixel-based analysis

Neuronal damage is the primary cause of long-term disability of multiple sclerosis (MS) patients. Assessment of axonal integrity from diffusion MRI parameters might enable better disease characterisation. 16 diffusion derived measurements from diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), and fixel-based analysis (FBA) in lesions, peri-lesion and normal appearing white matter were investigated. Diffusion MRI scans of 11 MS patients were processed to generate DTI, DKI, and FBA images. Fractional anisotropy (FA) and fibre density (FD) were used to assess axonal integrity across brain regions. Subsequently, 359 lesions were identified, and lesion and peri-lesion segmentation was performed using structural T₁w, T₂w, T₂w-FLAIR, and T₁w post-contrast MRI. The segmentations were then used to extract 16 diffusion MRI parameters from lesion, peri-lesion, and contralateral normal appearing white matter (NAWM). The measurements for axonal integrity, DTI-FA, DKI-FA, FBA-FD, produced similar results. All diffusion MRI parameters were affected in lesions as compared to NAWM (p < 0.001), confirming loss of axonal integrity in lesions. In peri-lesions, most parameters, except FBA-FD, were also significantly different from NAWM, although the effect size was smaller than in lesions. The reduction in axonal integrity in peri-lesions, despite unaffected fibre density estimates, suggests an effect of Wallerian degeneration.

Toward evaluation of multiresolution cortical thickness estimation with FreeSurfer, MaCRUISE, and BrainSuite

Advances in Magnetic Resonance Imaging hardware and methodologies allow for promoting the cortical morphometry with submillimeter spatial resolution. In this paper, we generated 3D self-enhanced high-resolution (HR) MRI imaging, by adapting 1 deep learning architecture, and 3 standard pipelines, FreeSurfer, MaCRUISE, and BrainSuite, have been collectively employed to evaluate the cortical thickness. We systematically investigated the differences in cortical thickness estimation for MRI sequences at multiresolution homologously originated from the native image. It has been revealed that there systematically exhibited the preferences in determining both inner and outer cortical surfaces at higher resolution, yielding most deeper cortical surface placements toward GM/WM or GM/CSF boundaries, which directs a consistent reduction tendency of mean cortical thickness estimation; on the contrary, the lower resolution data will most probably provide a more coarse and rough evaluation in cortical surface reconstruction, resulting in a relatively thicker estimation. Although the differences of cortical thickness estimation at the diverse spatial resolution varied with one another, almost all led to roughly one-sixth to one-fifth significant reduction across the entire brain at the HR, independent to the pipelines we applied, which emphasizes on generally coherent improved accuracy in a data-independent manner and endeavors to cost-efficiency with quantitative opportunities.

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.

Evaluation of white matter microstructure in pediatric onset multiple sclerosis with diffusion compartment imaging

BACKGROUND AND PURPOSE

Pediatric-onset multiple sclerosis (POMS) shows earlier axonal involvement and greater axonal loss than in adults. We aim to characterize the white matter (WM) microstructural changes in POMS using a diffusion compartment imaging (DCI) model and compare it to standard diffusion tensor imaging (DTI).

METHODS

Eleven patients (2 males, mean age 18.8 ± 3.9 years) with a diagnosis of relapsing and remitting POMS (mean age at disease onset 13.8 ± 2.9 years, mean duration 5.1 ± 1.9 years) and healthy controls (8 males, mean age 26.4 ± 6.5 years) were recruited and imaged at 3 T. A 90-gradient set Cube and Sphere acquisition and a novel DCI model known as DIstribution of Anisotropic MicrOstructural eNvironments with Diffusion-weighted imaging (DIAMOND) were used to calculate a single anisotropic compartment, an isotropic compartment, and a free diffusion compartment. Lesions and contralateral normal-appearing white matter (NAWM) in patients and whole brain WM for controls were labeled.

RESULTS

Eleven patients and 11 controls were recruited. When comparing lesions and contralateral NAWM in patients using DCI, compartmental axial diffusivity, radial diffusivity (cRD), and mean diffusivity (cMD) were higher in lesions. Conversely, compartmental fractional anisotropy (cFA) and heterogeneity index were lower in lesions. An analysis of DTI equivalents showed the same trends. In whole-brain NAWM of patients compared to controls, cRD and cMD were higher and cFA was lower in patients.

CONCLUSION

Lesions in POMS can be accurately characterized by a DCI model. Incipient changes in NAWM seen in DCI may not be readily observable by DTI.

Limited utility of structural MRI to identify the epileptogenic zone in young children with tuberous sclerosis

BACKGROUND AND PURPOSE

The success of epilepsy surgery in children with tuberous sclerosis complex (TSC) hinges on identification of the epileptogenic zone (EZ). We studied structural MRI markers of epileptogenic lesions in young children with TSC.

METHODS

We included 26 children with TSC who underwent epilepsy surgery before the age of 3 years at five sites, with 12 months or more follow-up. Two neuroradiologists, blinded to surgical outcome data, reviewed 10 candidate lesions on preoperative MRI for characteristics of the tuber (large affected area, calcification, cyst-like properties) and of focal cortical dysplasia (FCD) features (cortical malformation, gray-white matter junction blurring, transmantle sign). They selected lesions suspect for the EZ based on structural MRI, and reselected after unblinding to seizure onset location on electroencephalography (EEG).

RESULTS

None of the tuber characteristics and FCD features were distinctive for the EZ, indicated by resected lesions in seizure-free children. With structural MRI alone, the EZ was identified out of 10 lesions in 31%, and with addition of EEG data, this increased to 48%. However, rates of identification of resected lesions in non-seizure-free children were similar. Across 251 lesions, interrater agreement was moderate for large size (κ = .60), and fair (κ = .24) for all other features.

CONCLUSIONS

In young children with TSC, the utility of structural MRI features is limited in the identification of the epileptogenic tuber, but improves when combined with EEG data.

LGG-32. Integrated biologic, radiologic and clinical analysis of pediatric low-grade gliomas during and after targeted therapy treatment

BACKGROUND: Pediatric low grade gliomas (pLGGs) are the most common central nervous system tumor in children, characterized by driver alterations in the RAS and MAPK pathways. Genomic advances have facilitated use of molecular targeted therapies, however their long-term impact on tumor behavior remains critically unanswered. METHODS: We performed an IRB-approved, retrospective chart and imaging review of pLGGs treated with off-label targeted therapy at Dana-Farber/Boston Children’s Cancer and Blood Disorders Center from 2010 to 2020. Volumetric analysis was performed for BRAFV600E and BRAF fusion/duplication driven pLGG subsets. RESULTS: Fifty-five patients were identified (dabrafenib n = 15, everolimus n = 26, trametinib n = 11, and vemurafenib n = 3). Targeted agent was used as first or second-line therapy for 58% (32/55). Median duration of targeted therapy was 0.79 years (0.01 – 4.87), and overall median follow-up was 2.50 years (0.01 – 7.39). The 1-year, 3-year, and 5-year EFS from targeted therapy initiation were 62.1%, 38.2%, and 31.8%, respectively. Mean volumetric change for BRAFV600E mutated pLGG on BRAF inhibitors was -54.11%, and median time to best volumetric response was 8.28 months (n = 12). Median time to largest volume post-treatment was 2.86 months. Mean volumetric change for BRAF fusion/duplication pLGG on MEK inhibitors was +7.34% with median time to best volumetric response of 6.71 months (n = 7). Median time to largest volume post-treatment was 2.38 months. CONCLUSIONS: Our integrated clinical and volumetric data suggest the majority of patients receiving BRAF inhibitors or trametinib achieve reduction in tumor volume while on therapy and that tumor stability can be achieved following targeted therapy cessation. Moreover, volumetric analysis shows promise as a tool to assess targeted therapeutic response in pLGGs.

Scan-Specific Generative Neural Network for MRI Super-Resolution Reconstruction

The interpretation and analysis of Magnetic resonance imaging (MRI) benefit from high spatial resolution. Unfortunately, direct acquisition of high spatial resolution MRI is time-consuming and costly, which increases the potential for motion artifact, and suffers from reduced signal-to-noise ratio (SNR). Super-resolution reconstruction (SRR) is one of the most widely used methods in MRI since it allows for the trade-off between high spatial resolution, high SNR, and reduced scan times. Deep learning has emerged for improved SRR as compared to conventional methods. However, current deep learning-based SRR methods require large-scale training datasets of high-resolution images, which are practically difficult to obtain at a suitable SNR. We sought to develop a methodology that allows for dataset-free deep learning-based SRR, through which to construct images with higher spatial resolution and of higher SNR than can be practically obtained by direct Fourier encoding. We developed a dataset-free learning method that leverages a generative neural network trained for each specific scan or set of scans, which in turn, allows for SRR tailored to the individual patient. With the SRR from three short duration scans, we achieved high quality brain MRI at an isotropic spatial resolution of 0.125 cubic mm with six minutes of imaging time for T2 contrast and an average increase of 7.2 dB (34.2%) in SNR to these short duration scans. Motion compensation was achieved by aligning the three short duration scans together. We assessed our technique on simulated MRI data and clinical data acquired from 15 subjects. Extensive experimental results demonstrate that our approach achieved superior results to state-of-the-art methods, while in parallel, performed at reduced cost as scans delivered with direct high-resolution acquisition.

Fetal Brain Volume Predicts Neurodevelopment in Congenital Heart Disease

BACKGROUND

Neurodevelopmental impairment is common in children with congenital heart disease (CHD), but postnatal variables explain only 30% of the variance in outcomes. To explore whether the antecedents for neurodevelopmental disabilities might begin in utero, we analyzed whether fetal brain volume predicted subsequent neurodevelopmental outcome in children with CHD.

METHODS

Fetuses with isolated CHD and sociodemographically comparable healthy control fetuses underwent fetal brain magnetic resonance imaging and 2-year neurodevelopmental evaluation with the Bayley Scales of Infant and Toddler Development, Third Edition (Bayley-III) and the Adaptive Behavior Assessment System, Third Edition (ABAS-3). Hierarchical regression evaluated potential predictors of Bayley-III and ABAS-3 outcomes in the CHD group, including fetal total brain volume adjusted for gestational age and sex, sociodemographic characteristics, birth measures, and medical history.

RESULTS

The CHD group (n=52) had lower Bayley-III cognitive, language, and motor scores than the control group (n=26), but fetal brain volumes were similar. Within the CHD group, larger fetal total brain volume correlated with higher Bayley-III cognitive, language, and motor scores and ABAS-3 adaptive functioning scores (r=0.32-0.47; all P<0.05), but this was not noted in the control group. Fetal brain volume predicted 10% to 21% of the variance in neurodevelopmental outcome measures in univariate analyses. Multivariable models that also included social class and postnatal factors explained 18% to 45% of the variance in outcome, depending on developmental domain. Moreover, in final multivariable models, fetal brain volume was the most consistent predictor of neurodevelopmental outcome across domains.

CONCLUSIONS

Small fetal brain volume is a strong independent predictor of 2-year neurodevelopmental outcomes and may be an important imaging biomarker of future neurodevelopmental risk in CHD. Future studies are needed to support this hypothesis. Our findings support inclusion of fetal brain volume in risk stratification models and as a possible outcome in fetal neuroprotective intervention studies.

Normal Growth, Sexual Dimorphism, and Lateral Asymmetries at Fetal Brain MRI

Background Tools in image reconstruction, motion correction, and segmentation have enabled the accurate volumetric characterization of fetal brain growth at MRI. Purpose To evaluate the volumetric growth of intracranial structures in healthy fetuses, accounting for gestational age (GA), sex, and laterality with use of a spatiotemporal MRI atlas of fetal brain development. Materials and Methods T2-weighted 3.0-T half-Fourier acquired single-shot turbo spin-echo sequence MRI was performed in healthy fetuses from prospectively recruited pregnant volunteers from March 2013 to May 2019. A previously validated section-to-volume reconstruction algorithm was used to generate intensity-normalized superresolution three-dimensional volumes that were registered to a fetal brain MRI atlas with 28 anatomic regions of interest. Atlas-based segmentation was performed and manually refined. Labels included the bilateral hippocampus, amygdala, caudate nucleus, lentiform nucleus, thalamus, lateral ventricle, cerebellum, cortical plate, hemispheric white matter, internal capsule, ganglionic eminence, ventricular zone, corpus callosum, brainstem, hippocampal commissure, and extra-axial cerebrospinal fluid. For fetuses younger than 31 weeks of GA, the subplate and intermediate zones were delineated. A linear regression analysis was used to determine weekly age-related change adjusted for sex and laterality. Results The final analytic sample consisted of 122 MRI scans in 98 fetuses (mean GA, 29 weeks ± 5 [range, 20-38 weeks]). All structures had significant volume growth with increasing GA (P < .001). Weekly age-related change for individual structures in the brain parenchyma ranged from 2.0% (95% CI: 0.9, 3.1; P < .001) in the hippocampal commissure to 19.4% (95% CI: 18.7, 20.1; P < .001) in the cerebellum. The largest sex-related differences were 22.1% higher volume in male fetuses for the lateral ventricles (95% CI: 10.9, 34.4; P < .001). There was rightward volumetric asymmetry of 15.6% for the hippocampus (95% CI: 14.2, 17.2; P < .001) and leftward volumetric asymmetry of 8.1% for the lateral ventricles (95% CI: 3.7, 12.2; P < .001). Conclusion With use of a spatiotemporal MRI atlas, volumetric growth of the fetal brain showed complex trajectories dependent on structure, gestational age, sex, and laterality. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Rollins in this issue.

Brain charts for the human lifespan

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

Self-supervised IVIM DWI parameter estimation with a physics based forward model

PURPOSE

To assess the robustness and repeatability of intravoxel incoherent motion model (IVIM) parameter estimation for the diffusion-weighted MRI in the abdominal organs under the constraints of noisy diffusion signal using a novel neural network method.

METHODS

Clinically acquired abdominal scans of Crohn's disease patients were retrospectively analyzed with regions segmented in the kidney cortex, spleen, liver, and bowel. A novel IVIM parameter fitting method based on the principle of a physics guided self-supervised convolutional neural network that does not require reference parameter estimates for training was compared to a conventional non-linear least squares (NNLS) algorithm, and a voxelwise trained artificial neural network (ANN).

RESULTS

Results showed substantial increase in parameter robustness to the noise corrupted signal. In an intra-session repeatability experiment, the proposed method showed reduced coefficient of variation (CoV) over multiple acquisitions in comparison to conventional NLLS method and comparable performance to ANN. The use of D and f estimates from the proposed method led to the smallest misclassification error in linear discriminant analysis for characterization between normal and abnormal Crohn's disease bowel tissue. The fitting of D∗ parameter remains to be challenging.

CONCLUSION

The proposed method yields robust estimates of D and f IVIM parameters under the constraints of noisy diffusion signal. This indicates a potential for the use of the proposed method in conjunction with accelerated DW-MRI acquisition strategies, which would typically result in lower signal to noise ratio.

Reducing the Effects of Motion Artifacts in fMRI: A Structured Matrix Completion Approach

Functional MRI (fMRI) is widely used to study the functional organization of normal and pathological brains. However, the fMRI signal may be contaminated by subject motion artifacts that are only partially mitigated by motion correction strategies. These artifacts lead to distance-dependent biases in the inferred signal correlations. To mitigate these spurious effects, motion-corrupted volumes are censored from fMRI time series. Censoring can result in discontinuities in the fMRI signal, which may lead to substantial alterations in functional connectivity analysis. We propose a new approach to recover the missing entries from censoring based on structured low rank matrix completion. We formulated the artifact-reduction problem as the recovery of a super-resolved matrix from unprocessed fMRI measurements. We enforced a low rank prior on a large structured matrix, formed from the samples of the time series, to recover the missing entries. The recovered time series, in addition to being motion compensated, are also slice-time corrected at a fine temporal resolution. To achieve a fast and memory-efficient solution for our proposed optimization problem, we employed a variable splitting strategy. We validated the algorithm with simulations, data acquired under different motion conditions, and datasets from the ABCD study. Functional connectivity analysis showed that the proposed reconstruction resulted in connectivity matrices with lower errors in pair-wise correlation than non-censored and censored time series based on a standard processing pipeline. In addition, seed-based correlation analyses showed improved delineation of the default mode network. These demonstrate that the method can effectively reduce the adverse effects of motion in fMRI analysis.

Spatiotemporal changes in diffusivity and anisotropy in fetal brain tractography

Population averaged diffusion atlases can be utilized to characterize complex microstructural changes with less bias than data from individual subjects. In this study, a fetal diffusion tensor imaging (DTI) atlas was used to investigate tract-based changes in anisotropy and diffusivity in vivo from 23 to 38 weeks of gestational age (GA). Healthy pregnant volunteers with typically developing fetuses were imaged at 3 T. Acquisition included structural images processed with a super-resolution algorithm and DTI images processed with a motion-tracked slice-to-volume registration algorithm. The DTI from individual subjects were used to generate 16 templates, each specific to a week of GA; this was accomplished by means of a tensor-to-tensor diffeomorphic deformable registration method integrated with kernel regression in age. Deterministic tractography was performed to outline the forceps major, forceps minor, bilateral corticospinal tracts (CST), bilateral inferior fronto-occipital fasciculus (IFOF), bilateral inferior longitudinal fasciculus (ILF), and bilateral uncinate fasciculus (UF). The mean fractional anisotropy (FA) and mean diffusivity (MD) was recorded for all tracts. For a subset of tracts (forceps major, CST, and IFOF) we manually divided the tractograms into anatomy conforming segments to evaluate within-tract changes. We found tract-specific, nonlinear, age related changes in FA and MD. Early in gestation, these trends appear to be dominated by cytoarchitectonic changes in the transient white matter fetal zones while later in gestation, trends conforming to the progression of myelination were observed. We also observed significant (local) heterogeneity in within-tract developmental trajectories for the CST, IFOF, and forceps major.

An MRI method for parcellating the human striatum into matrix and striosome compartments in vivo

The mammalian striatum is comprised of intermingled tissue compartments, matrix and striosome. Though indistinguishable by routine histological techniques, matrix and striosome have distinct embryologic origins, afferent/efferent connections, surface protein expression, intra-striatal location, susceptibilities to injury, and functional roles in a range of animal behaviors. Distinguishing the compartments previously required post-mortem tissue and/or genetic manipulation; we aimed to identify matrix/striosome non-invasively in living humans. We used diffusion MRI (probabilistic tractography) to identify human striatal voxels with connectivity biased towards matrix-favoring or striosome-favoring regions (determined by prior animal tract-tracing studies). Segmented striatal compartments replicated the topological segregation and somatotopic organization identified in animal matrix/striosome studies. Of brain regions mapped in prior studies, our human brain data confirmed 93% of the compartment-selective structural connectivity demonstrated in animals. Test-retest assessment on repeat scans found a voxel classification error rate of 0.14%. Fractional anisotropy was significantly higher in matrix-like voxels, while mean diffusivity did not differ between the compartments. As mapped by the Talairach human brain atlas, 460 regions were significantly biased towards either matrix or striosome. Our method allows the study of striatal compartments in human health and disease, in vivo, for the first time.

Gradient-Guided Isotropic MRI Reconstruction from Anisotropic Acquisitions

The trade-off between image resolution, signal-to-noise ratio (SNR), and scan time in any magnetic resonance imaging (MRI) protocol is inevitable and unavoidable. Super-resolution reconstruction (SRR) has been shown effective in mitigating these factors, and thus, has become an important approach in addressing the current limitations of MRI. In this work, we developed a novel, image-based MRI SRR approach based on anisotropic acquisition schemes, which utilizes a new gradient guidance regularization method that guides the high-resolution (HR) reconstruction via a spatial gradient estimate. Further, we designed an analytical solution to propagate the spatial gradient fields from the low-resolution (LR) images to the HR image space and exploited these gradient fields over multiple scales with a dynamic update scheme for more accurate edge localization in the reconstruction. We also established a forward model of image formation and inverted it along with the proposed gradient guidance. The proposed SRR method allows subject motion between volumes and is able to incorporate various acquisition schemes where the LR images are acquired with arbitrary orientations and displacements, such as orthogonal and through-plane origin-shifted scans. We assessed our proposed approach on simulated data as well as on the data acquired on a Siemens 3T MRI scanner containing 45 MRI scans from 14 subjects. Our experimental results demonstrate that our approach achieved superior reconstructions compared to state-of-the-art methods, both in terms of local spatial smoothness and edge preservation, while, in parallel, at reduced, or at the same cost as scans delivered with direct HR acquisition.

Deep learning-based parameter estimation in fetal diffusion-weighted MRI

Diffusion-weighted magnetic resonance imaging (DW-MRI) of fetal brain is challenged by frequent fetal motion and signal to noise ratio that is much lower than non-fetal imaging. As a result, accurate and robust parameter estimation in fetal DW-MRI remains an open problem. Recently, deep learning techniques have been successfully used for DW-MRI parameter estimation in non-fetal subjects. However, none of those prior works has addressed the fetal brain because obtaining reliable fetal training data is challenging. To address this problem, in this work we propose a novel methodology that utilizes fetal scans as well as scans from prematurely-born infants. High-quality newborn scans are used to estimate accurate maps of the parameter of interest. These parameter maps are then used to generate DW-MRI data that match the measurement scheme and noise distribution that are characteristic of fetal data. In order to demonstrate the effectiveness and reliability of the proposed data generation pipeline, we used the generated data to train a convolutional neural network (CNN) to estimate color fractional anisotropy (CFA). We evaluated the trained CNN on independent sets of fetal data in terms of reconstruction accuracy, precision, and expert assessment of reconstruction quality. Results showed significantly lower reconstruction error (n=100,p<0.001) and higher reconstruction precision (n=20,p<0.001) for the proposed machine learning pipeline compared with standard estimation methods. Expert assessments on 20 fetal test scans showed significantly better overall reconstruction quality (p<0.001) and more accurate reconstruction of 11 regions of interest (p<0.001) with the proposed method.

Learning to estimate the fiber orientation distribution function from diffusion-weighted MRI

Estimation of white matter fiber orientation distribution function (fODF) is the essential first step for reliable brain tractography and connectivity analysis. Most of the existing fODF estimation methods rely on sub-optimal physical models of the diffusion signal or mathematical simplifications, which can impact the estimation accuracy. In this paper, we propose a data-driven method that avoids some of these pitfalls. Our proposed method is based on a multilayer perceptron that learns to map the diffusion-weighted measurements, interpolated onto a fixed spherical grid in the q space, to the target fODF. Importantly, we also propose methods for synthesizing reliable simulated training data. We show that the model can be effectively trained with simulated or real training data. Our phantom experiments show that the proposed method results in more accurate fODF estimation and tractography than several competing methods including the multi-tensor model, Bayesian estimation, spherical deconvolution, and two other machine learning techniques. On real data, we compare our method with other techniques in terms of accuracy of estimating the ground-truth fODF. The results show that our method is more accurate than other methods, and that it performs better than the competing methods when applied to under-sampled diffusion measurements. We also compare our method with the Sparse Fascicle Model in terms of expert ratings of the accuracy of reconstruction of several commissural, projection, association, and cerebellar tracts. The results show that the tracts reconstructed with the proposed method are rated significantly higher by three independent experts. Our study demonstrates the potential of data-driven methods for improving the accuracy and robustness of fODF estimation.

Super-resolution reconstruction of T2-weighted thick-slice neonatal brain MRI scans

Background and Purpose

Super‐resolutionreconstruction (SRR) can be used to reconstruct 3‐dimensional (3D) high‐resolution (HR) volume from several 2‐dimensional (2D) low‐resolution (LR) stacks of MRI slices. The purpose is to compare lengthy 2D T2‐weighted HR image acquisition of neonatal subjects with 3D SRR from several LR stacks in terms of image quality for clinical and morphometric assessments.

Methods

LR brain images were acquired from neonatal subjects to reconstruct isotropic 3D HR volumes by using SRR algorithm. Quality assessments were done by an experienced pediatric radiologist using scoring criteria adapted to newborn anatomical landmarks. The Wilcoxon signed‐rank test was used to compare scoring results between HR and SRR images. For quantitative assessments, morphology‐based segmentation was performed on both HR and SRR images and Dice coefficients between the results were computed. Additionally, simple linear regression was performed to compare the tissue volumes.

Results

No statistical difference was found between HR and SRR structural scores using Wilcoxon signed‐rank test (p = .63, Z = .48). Regarding segmentation results, R² values for the volumes of gray matter, white matter, cerebrospinal fluid, basal ganglia, cerebellum, and total brain volume including brain stem ranged between .95 and .99. Dice coefficients between the segmented regions from HR and SRR ranged between .83 ± .04 and .96 ± .01.

Conclusion

Qualitative and quantitative assessments showed that 3D SRR of several LR images produces images that are of comparable quality to standard 2D HR image acquisition for healthy neonatal imaging without loss of anatomical details with similar edge definition allowing the detection of fine anatomical structures and permitting comparable morphometric measurement.

The Connectivity Fingerprint of the Fusiform Gyrus Captures the Risk of Developing Autism in Infants with Tuberous Sclerosis Complex

Tuberous sclerosis complex (TSC) is a rare genetic disorder characterized by benign tumors throughout the body; it is generally diagnosed early in life and has a high prevalence of autism spectrum disorder (ASD), making it uniquely valuable in studying the early development of autism, before neuropsychiatric symptoms become apparent. One well-documented deficit in ASD is an impairment in face processing. In this work, we assessed whether anatomical connectivity patterns of the fusiform gyrus, a central structure in face processing, capture the risk of developing autism early in life. We longitudinally imaged TSC patients at 1, 2, and 3 years of age with diffusion compartment imaging. We evaluated whether the anatomical connectivity fingerprint of the fusiform gyrus was associated with the risk of developing autism measured by the Autism Observation Scale for Infants (AOSI). Our findings suggest that the fusiform gyrus connectivity captures the risk of developing autism as early as 1 year of age and provides evidence that abnormal fusiform gyrus connectivity increases with age. Moreover, the identified connections that best capture the risk of developing autism involved the fusiform gyrus and limbic and paralimbic regions that were consistent with the ASD phenotype, involving an increased number of left-lateralized structures with increasing age.

Free-breathing radial stack-of-stars three-dimensional Dixon gradient echo sequence in abdominal magnetic resonance imaging in sedated pediatric patients

BACKGROUND

There is a strong need for improvements in motion robust T1-weighted abdominal imaging sequences in children to enable high-quality, free-breathing imaging.

OBJECTIVE

To compare imaging time and quality of a radial stack-of-stars, free-breathing T1-weighted gradient echo acquisition (volumetric interpolated breath-hold examination [VIBE]) three-dimensional (3-D) Dixon sequence in sedated pediatric patients undergoing abdominal magnetic resonance imaging (MRI) against conventional Cartesian T1-weighed sequences.

MATERIALS AND METHODS

This study was approved by the institutional review board with informed consent obtained from all subjects. Study subjects included 31 pediatric patients (19 male, 12 female; median age: 5 years; interquartile range: 5 years) undergoing abdominal MRI at 3 tesla with a free-breathing T1-weighted radial stack-of-stars 3-D VIBE Dixon prototype sequence, StarVIBE Dixon (radial technique), between October 2018 and June 2019 with previous abdominal MR imaging using conventional Cartesian T1-weighed imaging (traditional technique). MRI component times were recorded as well as the total number of non-contrast T1-weighted sequences. Two radiologists independently rated images for quality using a scale from 1 to 5 according to the following metrics: overall image quality, hepatic edge sharpness, hepatic vessel clarity and respiratory motion robustness. Scores were compared between the groups.

RESULTS

Mean T1-weighted imaging times for all subjects were 3.63 min for radial exams and 8.01 min for traditional exams (P<0.001), and total non-contrast imaging time was 32.7 min vs. 43.9 min (P=0.002). Adjusted mean total MRI time for all subjects was 60.2 min for radial exams and 65.7 min for traditional exams (P=0.387). The mean number of non-contrast T1-weighted sequences performed in radial MRI exams was 1.0 compared to 1.9 (range: 0-6) in traditional exams (P<0.001). StarVIBE Dixon outperformed Cartesian methods in all quality metrics. The mean overall image quality (scale 1-5) was 3.95 for radial exams and 3.31 for traditional exams (P<0.001).

CONCLUSION

Radial stack-of-stars 3-D VIBE Dixon during free-breathing abdominal MRI in pediatric patients offers improved image quality compared to Cartesian T1-weighted imaging techniques with decreased T1-weighted and total non-contrast imaging time. This has important implications for children undergoing sedation for imaging.

A machine learning-based method for estimating the number and orientations of major fascicles in diffusion-weighted magnetic resonance imaging

Accurate modeling of diffusion-weighted magnetic resonance imaging measurements is necessary for accurate brain connectivity analysis. Existing methods for estimating the number and orientations of fascicles in an imaging voxel either depend on non-convex optimization techniques that are sensitive to initialization and measurement noise, or are prone to predicting spurious fascicles. In this paper, we propose a machine learning-based technique that can accurately estimate the number and orientations of fascicles in a voxel. Our method can be trained with either simulated or real diffusion-weighted imaging data. Our method estimates the angle to the closest fascicle for each direction in a set of discrete directions uniformly spread on the unit sphere. This information is then processed to extract the number and orientations of fascicles in a voxel. On realistic simulated phantom data with known ground truth, our method predicts the number and orientations of crossing fascicles more accurately than several classical and machine learning methods. It also leads to more accurate tractography. On real data, our method is better than or compares favorably with other methods in terms of robustness to measurement down-sampling and also in terms of expert quality assessment of tractography results.

Physics-based iterative reconstruction for dual-source and flying focal spot computed tomography

PURPOSE

For single-source helical Computed Tomography (CT), both Filtered-Back Projection (FBP) and statistical iterative reconstruction have been investigated. However, for dual-source CT with flying focal spot (DS-FFS CT), a statistical iterative reconstruction that accurately models the scanner geometry and acquisition physics remains unknown to researchers. Therefore, our purpose is to present a novel physics-based iterative reconstruction method for DS-FFS CT and assess its image quality.

METHODS

Our algorithm uses precise physics models to reconstruct from the native cone-beam geometry and interleaved dual-source helical trajectory of a DS-FFS CT. To do so, we construct a noise physics model to represent data acquisition noise and a prior image model to represent image noise and texture. In addition, we design forward system models to compute the locations of deflected focal spots, the dimension, and sensitivity of voxels and detector units, as well as the length of intersection between x-rays and voxels. The forward system models further represent the coordinated movement between the dual sources by computing their x-ray coverage gaps and overlaps at an arbitrary helical pitch. With the above models, we reconstruct images by an advanced Consensus Equilibrium (CE) numerical method to compute the maximum a posteriori estimate to a joint optimization problem that simultaneously fits all models.

RESULTS

We compared our reconstruction with Siemens ADMIRE, which is the clinical standard hybrid iterative reconstruction (IR) method for DS-FFS CT, in terms of spatial resolution, noise profile, and image artifacts through both phantoms and clinical scan datasets. Experiments show that our reconstruction has a higher spatial resolution, with a Task-Based Modulation Transfer Function (MTF_task ) consistently higher than the clinical standard hybrid IR. In addition, our reconstruction shows a reduced magnitude of image undersampling artifacts than the clinical standard.

CONCLUSIONS

By modeling a precise geometry and avoiding data rebinning or interpolation, our physics-based reconstruction achieves a higher spatial resolution and fewer image artifacts with smaller magnitude than the clinical standard hybrid IR.

Fast and High-Resolution Neonatal Brain MRI Through Super-Resolution Reconstruction From Acquisitions With Variable Slice Selection Direction

The brain of neonates is small in comparison to adults. Imaging at typical resolutions such as one cubic mm incurs more partial voluming artifacts in a neonate than in an adult. The interpretation and analysis of MRI of the neonatal brain benefit from a reduction in partial volume averaging that can be achieved with high spatial resolution. Unfortunately, direct acquisition of high spatial resolution MRI is slow, which increases the potential for motion artifact, and suffers from reduced signal-to-noise ratio. The purpose of this study is thus that using super-resolution reconstruction in conjunction with fast imaging protocols to construct neonatal brain MRI images at a suitable signal-to-noise ratio and with higher spatial resolution than can be practically obtained by direct Fourier encoding. We achieved high quality brain MRI at a spatial resolution of isotropic 0.4 mm with 6 min of imaging time, using super-resolution reconstruction from three short duration scans with variable directions of slice selection. Motion compensation was achieved by aligning the three short duration scans together. We applied this technique to 20 newborns and assessed the quality of the images we reconstructed. Experiments show that our approach to super-resolution reconstruction achieved considerable improvement in spatial resolution and signal-to-noise ratio, while, in parallel, substantially reduced scan times, as compared to direct high-resolution acquisitions. The experimental results demonstrate that our approach allowed for fast and high-quality neonatal brain MRI for both scientific research and clinical studies.

Quantitative In vivo MRI Assessment of Structural Asymmetries and Sexual Dimorphism of Transient Fetal Compartments in the Human Brain

Structural asymmetries and sexual dimorphism of the human cerebral cortex have been identified in newborns, infants, children, adolescents, and adults. Some of these findings were linked with cognitive and neuropsychiatric disorders, which have roots in altered prenatal brain development. However, little is known about structural asymmetries or sexual dimorphism of transient fetal compartments that arise in utero. Thus, we aimed to identify structural asymmetries and sexual dimorphism in the volume of transient fetal compartments (cortical plate [CP] and subplate [SP]) across 22 regions. For this purpose, we used in vivo structural T2-weighted MRIs of 42 healthy fetuses (16.43-36.86 gestational weeks old, 15 females). We found significant leftward asymmetry in the volume of the CP and SP in the inferior frontal gyrus. The orbitofrontal cortex showed significant rightward asymmetry in the volume of CP merged with SP. Males had significantly larger volumes in regions belonging to limbic, occipital, and frontal lobes, which were driven by a significantly larger SP. Lastly, we did not observe sexual dimorphism in the growth trajectories of the CP or SP. In conclusion, these results support the hypothesis that structural asymmetries and sexual dimorphism in relative volumes of cortical regions are present during prenatal brain development.

Matched neurofeedback during fMRI differentially activates reward-related circuits in active and sham groups

BACKGROUND AND PURPOSE

Functional MRI neurofeedback (fMRI-nf) leverages the brain's ability to self-regulate its own activity. However, self-regulation processes engaged during fMRI-nf are incompletely understood. Here, we used matched feedback in an fMRI-nf experimental protocol to investigate whether brain processes recognize true neurofeedback signals.

METHODS

We implemented an existing fMRI-nf protocol to train lateralized motor activity using a finger-tap task in conjunction with real-time feedback. Twelve healthy, right-handed, adult participants were assigned into age- and sex-matched active and sham study groups. Matched participant pairs received the same visual feedback, based on brain activity of the participant from the active group. We compared group-averaged activation maps before, during, and after neurofeedback, and analyzed changes in lateralized motor activity due to neurofeedback.

RESULTS

Active and sham groups demonstrated different brain activation to the same feedback during neurofeedback. In particular, there was higher activation in visual cortex, secondary somatosensory cortex, and right inferior frontal gyrus in the active group compared to the sham group. Conversely, sham participants demonstrated higher activation in anterior cingulate cortex, left frontal pole, and posterior superior temporal gyrus. Despite differing brain activations during neurofeedback, neither group demonstrated significant improvement in lateralized motor activity from pre to postfeedback scan in the same session. We also observed no significant difference between pre and postfeedback activation maps, suggesting that no significant finger-tap related functional reorganization had occurred.

CONCLUSIONS

These findings suggest that fMRI neurofeedback paradigms that monitor or incorporate activity from regions reported here would provide enhanced efficacy for research investigation and clinical intervention.

A structural brain network of genetic vulnerability to psychiatric illness

Psychiatry is undergoing a paradigm shift from the acceptance of distinct diagnoses to a representation of psychiatric illness that crosses diagnostic boundaries. How this transition is supported by a shared neurobiology remains largely unknown. In this study, we first identify single nucleotide polymorphisms (SNPs) associated with psychiatric disorders based on 136 genome-wide association studies. We then conduct a joint analysis of these SNPs and brain structural connectomes in 678 healthy children in the PING study. We discovered a strong, robust, and transdiagnostic mode of genome-connectome covariation which is positively and specifically correlated with genetic risk for psychiatric illness at the level of individual SNPs. Similarly, this mode is also significantly positively correlated with polygenic risk scores for schizophrenia, alcohol use disorder, major depressive disorder, a combined bipolar disorder-schizophrenia phenotype, and a broader cross-disorder phenotype, and significantly negatively correlated with a polygenic risk score for educational attainment. The resulting "vulnerability network" is shown to mediate the influence of genetic risks onto behaviors related to psychiatric vulnerability (e.g., marijuana, alcohol, and caffeine misuse, perceived stress, and impulsive behavior). Its anatomy overlaps with the default-mode network, with a network of cognitive control, and with the occipital cortex. These findings suggest that the brain vulnerability network represents an endophenotype funneling genetic risks for various psychiatric illnesses through a common neurobiological root. It may form part of the neural underpinning of the well-recognized but poorly explained overlap and comorbidity between psychiatric disorders.

Transfer learning in medical image segmentation: New insights from analysis of the dynamics of model parameters and learned representations

We present a critical assessment of the role of transfer learning in training fully convolutional networks (FCNs) for medical image segmentation. We first show that although transfer learning reduces the training time on the target task, improvements in segmentation accuracy are highly task/data-dependent. Large improvements are observed only when the segmentation task is more challenging and the target training data is smaller. We shed light on these observations by investigating the impact of transfer learning on the evolution of model parameters and learned representations. We observe that convolutional filters change little during training and still look random at convergence. We further show that quite accurate FCNs can be built by freezing the encoder section of the network at random values and only training the decoder section. At least for medical image segmentation, this finding challenges the common belief that the encoder section needs to learn data/task-specific representations. We examine the evolution of FCN representations to gain a deeper insight into the effects of transfer learning on the training dynamics. Our analysis shows that although FCNs trained via transfer learning learn different representations than FCNs trained with random initialization, the variability among FCNs trained via transfer learning can be as high as that among FCNs trained with random initialization. Moreover, feature reuse is not restricted to the early encoder layers; rather, it can be more significant in deeper layers. These findings offer new insights and suggest alternative ways of training FCNs for medical image segmentation.

Tuber Locations Associated with Infantile Spasms Map to a Common Brain Network

OBJECTIVE

Approximately 50% of patients with tuberous sclerosis complex develop infantile spasms, a sudden onset epilepsy syndrome associated with poor neurological outcomes. An increased burden of tubers confers an elevated risk of infantile spasms, but it remains unknown whether some tuber locations confer higher risk than others. Here, we test whether tuber location and connectivity are associated with infantile spasms.

METHODS

We segmented tubers from 123 children with (n = 74) and without (n = 49) infantile spasms from a prospective observational cohort. We used voxelwise lesion symptom mapping to test for an association between spasms and tuber location. We then used lesion network mapping to test for an association between spasms and connectivity with tuber locations. Finally, we tested the discriminability of identified associations with logistic regression and cross-validation as well as statistical mediation.

RESULTS

Tuber locations associated with infantile spasms were heterogenous, and no single location was significantly associated with spasms. However, >95% of tuber locations associated with spasms were functionally connected to the globi pallidi and cerebellar vermis. These connections were specific compared to tubers in patients without spasms. Logistic regression found that globus pallidus connectivity was a stronger predictor of spasms (odds ratio [OR] = 1.96, 95% confidence interval [CI] = 1.10-3.50, p = 0.02) than tuber burden (OR = 1.65, 95% CI = 0.90-3.04, p = 0.11), with a mean receiver operating characteristic area under the curve of 0.73 (±0.1) during repeated cross-validation.

INTERPRETATION

Connectivity between tuber locations and the bilateral globi pallidi is associated with infantile spasms. Our findings lend insight into spasm pathophysiology and may identify patients at risk. ANN NEUROL 2021;89:726-739.

Magnetic Resonance Imaging (MRI) and Spectroscopy in Succinic Semialdehyde Dehydrogenase Deficiency

Succinic semialdehyde dehydrogenase (SSADH) deficiency is an autosomal recessive disorder of γ-aminobutyric acid (GABA) degradation, resulting in elevations of brain GABA and γ-hydroxybutyric acid (GHB). Previous magnetic resonance (MR) spectroscopy studies have shown increased levels of Glx in SSADH deficiency patients. Here in this work, we measure brain GABA in a large cohort of SSADH deficiency patients using advanced MR spectroscopy techniques that allow separation of GABA from overlapping metabolite peaks. We observed significant increases in GABA concentrations in SSADH deficiency patients for all 3 brain regions that were evaluated. Although GABA levels were higher in all 3 regions, each region had different patterns in terms of GABA changes with respect to age. We also report results from structural magnetic resonance imaging (MRI) of the same cohort compared with age-matched controls. We consistently observed signal hyperintensities in globus pallidus and cerebellar dentate nucleus.

Magic DIAMOND: Multi-fascicle diffusion compartment imaging with tensor distribution modeling and tensor-valued diffusion encoding

Diffusion tensor imaging provides increased sensitivity to microstructural tissue changes compared to conventional anatomical imaging but also presents limited specificity. To tackle this problem, the DIAMOND model subdivides the voxel content into diffusion compartments and draws from diffusion-weighted data to estimate compartmental non-central matrix-variate Gamma distributions of diffusion tensors. It models each sub-voxel fascicle separately, resolving crossing white-matter pathways and allowing for a fascicle-element (fixel) based analysis of microstructural features. Alternatively, specific features of the intra-voxel diffusion tensor distribution can be selectively measured using tensor-valued diffusion-weighted acquisition schemes. However, the impact of such schemes on estimating brain microstructural features has only been studied in a handful of parametric single-fascicle models. In this work, we derive a general Laplace transform for the non-central matrix-variate Gamma distribution, which enables the extension of DIAMOND to tensor-valued encoded data. We then evaluate this "Magic DIAMOND" model in silico and in vivo on various combinations of tensor-valued encoded data. Assessing uncertainty on parameter estimation via stratified bootstrap, we investigate both voxel-based and fixel-based metrics by carrying out multi-peak tractography. We demonstrate using in silico evaluations that tensor-valued diffusion encoding significantly improves Magic DIAMOND's accuracy. Most importantly, we show in vivo that our estimated metrics can be robustly mapped along tracks across regions of fiber crossing, which opens new perspectives for tractometry and microstructure mapping along specific white-matter tracts.