Development of human white matter pathways in utero over the second and third trimester.

During the second and third trimesters of human gestation, rapid neurodevelopment is underpinned by fundamental processes including neuronal migration, cellular organization, cortical layering, and myelination. In this time, white matter growth and maturation lay the foundation for an efficient network of structural connections. Detailed knowledge about this developmental trajectory in the healthy human fetal brain is limited, in part, due to the inherent challenges of acquiring high-quality MRI data from this population. Here, we use state-of-the-art high-resolution multishell motion-corrected diffusion-weighted MRI (dMRI), collected as part of the developing Human Connectome Project (dHCP), to characterize the in utero maturation of white matter microstructure in 113 fetuses aged 22 to 37 wk gestation. We define five major white matter bundles and characterize their microstructural features using both traditional diffusion tensor and multishell multitissue models. We found unique maturational trends in thalamocortical fibers compared with association tracts and identified different maturational trends within specific sections of the corpus callosum. While linear maturational increases in fractional anisotropy were seen in the splenium of the corpus callosum, complex nonlinear trends were seen in the majority of other white matter tracts, with an initial decrease in fractional anisotropy in early gestation followed by a later increase. The latter is of particular interest as it differs markedly from the trends previously described in ex utero preterm infants, suggesting that this normative fetal data can provide significant insights into the abnormalities in connectivity which underlie the neurodevelopmental impairments associated with preterm birth.

Diffusion MRI of the facial-vestibulocochlear nerve complex: a prospective clinical validation study.

OBJECTIVES: Surgical planning of vestibular schwannoma surgery would benefit greatly from a robust method of delineating the facial-vestibulocochlear nerve complex with respect to the tumour. This study aimed to optimise a multi-shell readout-segmented diffusion-weighted imaging (rs-DWI) protocol and develop a novel post-processing pipeline to delineate the facial-vestibulocochlear complex within the skull base region, evaluating its accuracy intraoperatively using neuronavigation and tracked electrophysiological recordings. METHODS: In a prospective study of five healthy volunteers and five patients who underwent vestibular schwannoma surgery, rs-DWI was performed and colour tissue maps (CTM) and probabilistic tractography of the cranial nerves were generated. In patients, the average symmetric surface distance (ASSD) and 95% Hausdorff distance (HD-95) were calculated with reference to the neuroradiologist-approved facial nerve segmentation. The accuracy of patient results was assessed intraoperatively using neuronavigation and tracked electrophysiological recordings. RESULTS: Using CTM alone, the facial-vestibulocochlear complex of healthy volunteer subjects was visualised on 9/10 sides. CTM were generated in all 5 patients with vestibular schwannoma enabling the facial nerve to be accurately identified preoperatively. The mean ASSD between the annotators' two segmentations was 1.11 mm (SD 0.40) and the mean HD-95 was 4.62 mm (SD 1.78). The median distance from the nerve segmentation to a positive stimulation point was 1.21 mm (IQR 0.81-3.27 mm) and 2.03 mm (IQR 0.99-3.84 mm) for the two annotators, respectively. CONCLUSIONS: rs-DWI may be used to acquire dMRI data of the cranial nerves within the posterior fossa. CLINICAL RELEVANCE STATEMENT: Readout-segmented diffusion-weighted imaging and colour tissue mapping provide 1-2 mm spatially accurate imaging of the facial-vestibulocochlear nerve complex, enabling accurate preoperative localisation of the facial nerve. This study evaluated the technique in 5 healthy volunteers and 5 patients with vestibular schwannoma. KEY POINTS: • Readout-segmented diffusion-weighted imaging (rs-DWI) with colour tissue mapping (CTM) visualised the facial-vestibulocochlear nerve complex on 9/10 sides in 5 healthy volunteer subjects. • Using rs-DWI and CTM, the facial nerve was visualised in all 5 patients with vestibular schwannoma and within 1.21-2.03 mm of the nerve's true intraoperative location. • Reproducible results were obtained on different scanners.

An MR fingerprinting approach for quantitative inhomogeneous magnetization transfer imaging.

PURPOSE: Magnetization transfer (MT) and inhomogeneous MT (ihMT) contrasts are used in MRI to provide information about macromolecular tissue content. In particular, MT is sensitive to macromolecules, and ihMT appears to be specific to myelinated tissue. This study proposes a technique to characterize MT and ihMT properties from a single acquisition, producing both semiquantitative contrast ratios and quantitative parameter maps. THEORY AND METHODS: Building on previous work that uses multiband RF pulses to efficiently generate ihMT contrast, we propose a cyclic steady-state approach that cycles between multiband and single-band pulses to boost the achieved contrast. Resultant time-variable signals are reminiscent of an MR fingerprinting acquisition, except that the signal fluctuations are entirely mediated by MT effects. A dictionary-based low-rank inversion method is used to reconstruct the resulting images and to produce both semiquantitative MT ratio and ihMT ratio maps, as well as quantitative parameter estimates corresponding to an ihMT tissue model. RESULTS: Phantom and in vivo brain data acquired at 1.5 Tesla demonstrate the expected contrast trends, with ihMT ratio maps showing contrast more specific to white matter, as has been reported by others. Quantitative estimation of semisolid fraction and dipolar T1 was also possible and yielded measurements consistent with literature values in the brain. CONCLUSION: By cycling between multiband and single-band pulses, an entirely MT-mediated fingerprinting method was demonstrated. This proof-of-concept approach can be used to generate semiquantitative maps and quantitatively estimate some macromolecular-specific tissue parameters.

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 â€‹mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.

A data-driven approach to optimising the encoding for multi-shell diffusion MRI with application to neonatal imaging.

Diffusion MRI has the potential to provide important information about the connectivity and microstructure of the human brain during normal and abnormal development, noninvasively and in vivo. Recent developments in MRI hardware and reconstruction methods now permit the acquisition of large amounts of data within relatively short scan times. This makes it possible to acquire more informative multi-shell data, with diffusion sensitisation applied along many directions over multiple b-value shells. Such schemes are characterised by the number of shells acquired, and the specific b-value and number of directions sampled for each shell. However, there is currently no clear consensus as to how to optimise these parameters. In this work, we propose a means of optimising multi-shell acquisition schemes by estimating the information content of the diffusion MRI signal, and optimising the acquisition parameters for sensitivity to the observed effects, in a manner agnostic to any particular diffusion analysis method that might subsequently be applied to the data. This method was used to design the acquisition scheme for the neonatal diffusion MRI sequence used in the developing Human Connectome Project (dHCP), which aims to acquire high quality data and make it freely available to the research community. The final protocol selected by the algorithm, and currently in use within the dHCP, consists of 20 b=0 images and diffusion-weighted images at b = 400, 1000 and 2600 s/mm2 with 64, 88 and 128 directions per shell, respectively.

Inherent and unpredictable bias in multi-component DESPOT myelin water fraction estimation.

Multicomponent driven equilibrium steady-state observation of T1 and T2 (mcDESPOT) aims to quantify the Myelin Water Fraction (MWF) using a two-pool microstructural model. The MWF has been used to track neurodevelopment and neurodegeneration and has been histologically correlated to myelin content. mcDESPOT has a clinically feasible acquisition time and high signal-to-noise ratio (SNR) relative to other MWF techniques. However, disagreement exists in the literature between experimental studies that show MWF maps with plausible grey matter-white matter (GM-WM) contrast and theoretical work that questions the accuracy and precision of mcDESPOT. We demonstrate that mcDESPOT parameter estimation is inaccurate and imprecise if intercompartmental exchange is included in the microstructural model, but that significant bias results if exchange is neglected. The source of apparent MWF contrast is likely due to the complex convergence behaviour of the Stochastic Region Contraction (SRC) method commonly used to fit the mcDESPOT model. mcDESPOT-derived parameter estimates are hence not directly relatable to the underlying microstructural model and are only comparable to others using similar acquisition schemes and fitting constraints.

Automated processing pipeline for neonatal diffusion MRI in the developing Human Connectome Project.

The developing Human Connectome Project is set to create and make available to the scientific community a 4-dimensional map of functional and structural cerebral connectivity from 20 to 44 weeks post-menstrual age, to allow exploration of the genetic and environmental influences on brain development, and the relation between connectivity and neurocognitive function. A large set of multi-modal MRI data from fetuses and newborn infants is currently being acquired, along with genetic, clinical and developmental information. In this overview, we describe the neonatal diffusion MRI (dMRI) image processing pipeline and the structural connectivity aspect of the project. Neonatal dMRI data poses specific challenges, and standard analysis techniques used for adult data are not directly applicable. We have developed a processing pipeline that deals directly with neonatal-specific issues, such as severe motion and motion-related artefacts, small brain sizes, high brain water content and reduced anisotropy. This pipeline allows automated analysis of in-vivo dMRI data, probes tissue microstructure, reconstructs a number of major white matter tracts, and includes an automated quality control framework that identifies processing issues or inconsistencies. We here describe the pipeline and present an exemplar analysis of data from 140 infants imaged at 38-44 weeks post-menstrual age.

The developing human connectome project: A minimal processing pipeline for neonatal cortical surface reconstruction.

The Developing Human Connectome Project (dHCP) seeks to create the first 4-dimensional connectome of early life. Understanding this connectome in detail may provide insights into normal as well as abnormal patterns of brain development. Following established best practices adopted by the WU-MINN Human Connectome Project (HCP), and pioneered by FreeSurfer, the project utilises cortical surface-based processing pipelines. In this paper, we propose a fully automated processing pipeline for the structural Magnetic Resonance Imaging (MRI) of the developing neonatal brain. This proposed pipeline consists of a refined framework for cortical and sub-cortical volume segmentation, cortical surface extraction, and cortical surface inflation, which has been specifically designed to address considerable differences between adult and neonatal brains, as imaged using MRI. Using the proposed pipeline our results demonstrate that images collected from 465 subjects ranging from 28 to 45 weeks post-menstrual age (PMA) can be processed fully automatically; generating cortical surface models that are topologically correct, and correspond well with manual evaluations of tissue boundaries in 85% of cases. Results improve on state-of-the-art neonatal tissue segmentation models and significant errors were found in only 2% of cases, where these corresponded to subjects with high motion. Downstream, these surfaces will enhance comparisons of functional and diffusion MRI datasets, supporting the modelling of emerging patterns of brain connectivity.

The role of whole-brain diffusion MRI as a tool for studying human in vivo cortical segregation based on a measure of neurite density.

PURPOSE: To investigate whether diffusion MRI can be used to study cortical segregation based on a contrast related to neurite density, thus providing a complementary tool to myelin-based MRI techniques used for myeloarchitecture. METHODS: Several myelin-sensitive MRI methods (e.g., based on T1 , T2 , and T2*) have been proposed to parcellate cortical areas based on their myeloarchitecture. Recent improvements in hardware, acquisition, and analysis methods have opened the possibility of achieving a more robust characterization of cortical microstructure using diffusion MRI. High-quality diffusion MRI data from the Human Connectome Project was combined with recent advances in fiber orientation modeling. The orientational average of the fiber orientation distribution was used as a summary parameter, which was displayed as inflated brain surface views. RESULTS: Diffusion MRI identifies cortical patterns consistent with those previously seen by MRI methods used for studying myeloarchitecture, which have shown patterns of high myelination in the sensorimotor strip, visual cortex, and auditory areas and low myelination in frontal and anterior temporal areas. CONCLUSION: In vivo human diffusion MRI provides a useful complementary noninvasive approach to myelin-based methods used to study whole-brain cortical parcellation, by exploiting a contrast based on tissue microstructure related to neurite density, rather than myelin itself. Magn Reson Med 79:2738-2744, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Tract-specific atrophy in focal epilepsy: Disease, genetics, or seizures?

OBJECTIVE: To investigate whether genetics, underlying pathology, or repeated seizures contribute to atrophy in specific white matter tracts. METHODS: Medically refractory unilateral temporal lobe epilepsy (TLE) with hippocampal sclerosis (HS-TLE, n = 26) was studied as an archetype of focal epilepsy, using fixel-based analysis of diffusion-weighted imaging. A genetic effect was assessed in first-degree relatives of HS-TLE subjects who did not have epilepsy themselves (HS-1°Rel; n = 26). The role of disease process was uncovered by comparing HS-TLE to unilateral TLE with normal clinical magnetic resonance imaging (MRI-neg TLE; n = 26, matched for seizure severity). The effect of focal seizures was inferred from lateralized atrophy common to both HS-TLE and MRI-neg TLE, in comparison to healthy controls (n = 76). RESULTS: HS-1 °Rel had bilaterally small hippocampi, but no focal white matter atrophy was detected, indicating a limited effect of genetics. HS-TLE subjects had lateralized atrophy of most temporal lobe tracts, and hippocampal volumes in HS-TLE correlated with parahippocampal cingulum and anterior commissure atrophy, indicating an effect of the underlying pathology. Ipsilateral atrophy of the tapetum, uncinate, and inferior fronto-occipital fasciculus was found in both HS-TLE and MRI-neg TLE, suggesting a common lateralized effect of focal seizures. Both epilepsy groups had bilateral atrophy of the dorsal cingulum and corpus callosum fibers, which we interpret as a consequence of bilateral insults (potentially generalized seizures and/or medications). INTERPRETATION: Underlying pathology, repeated focal seizures, and global insults each contribute to atrophy in specific tracts. Genetic factors make less of a contribution in this cohort. A multifactorial model of white matter atrophy in focal epilepsy is proposed. Ann Neurol 2017;81:240-250.

Predicting hand function after hemidisconnection.

Hemidisconnections (i.e. hemispherectomies or hemispherotomies) invariably lead to contralateral hemiparesis. Many patients with a pre-existing hemiparesis, however, experience no deterioration in motor functions, and some can still grasp with their paretic hand after hemidisconnection. The scope of our study was to predict this phenomenon. Hypothesizing that preserved contralateral grasping ability after hemidisconnection can only occur in patients controlling their paretic hands via ipsilateral corticospinal projections already in the preoperative situation, we analysed the asymmetries of the brainstem (by manual magnetic resonance imaging volumetry) and of the structural connectivity of the corticospinal tracts within the brainstem (by magnetic resonance imaging diffusion tractography), assuming that marked hypoplasia or Wallerian degeneration on the lesioned side in patients who can grasp with their paretic hands indicate ipsilateral control. One hundred and two patients who underwent hemidisconnections between 0.8 and 36 years of age were included. Before the operation, contralateral hand function was normal in 3/102 patients, 47/102 patients showed hemiparetic grasping ability and 52/102 patients could not grasp with their paretic hands. After hemidisconnection, 20/102 patients showed a preserved grasping ability, and 5/102 patients began to grasp with their paretic hands only after the operation. All these 25 patients suffered from pre- or perinatal brain lesions. Thirty of 102 patients lost their grasping ability. This group included all seven patients with a post-neonatally acquired or progressive brain lesion who could grasp before the operation, and also all three patients with a preoperatively normal hand function. The remaining 52/102 patients were unable to grasp pre- and postoperatively. On magnetic resonance imaging, the patients with preserved grasping showed significantly more asymmetric brainstem volumes than the patients who lost their grasping ability. Similarly, these patients showed striking asymmetries in the structural connectivity of the corticospinal tracts. In summary, normal preoperative hand function and a post-neonatally acquired or progressive lesion predict a loss of grasping ability after hemidisconnection. A postoperatively preserved grasping ability is possible in hemiparetic patients with pre- or perinatal lesions, and this is highly likely when the brainstem is asymmetric and especially when the structural connectivity of the corticospinal tracts within the brainstem is asymmetric.

The effects of SIFT on the reproducibility and biological accuracy of the structural connectome.

Diffusion MRI streamlines tractography is increasingly being used to characterise and assess the structural connectome of the human brain. However, issues pertaining to quantification of structural connectivity using streamlines reconstructions are well-established in the field, and therefore the validity of any conclusions that may be drawn from these analyses remains ambiguous. We recently proposed a post-processing method entitled "SIFT: Spherical-deconvolution Informed Filtering of Tractograms" as a mechanism for reducing the biases in quantitative measures of connectivity introduced by the streamlines reconstruction method. Here, we demonstrate the advantage of this approach in the context of connectomics in three steps. Firstly, we carefully consider the model imposed by the SIFT method, and the implications this has for connectivity quantification. Secondly, we investigate the effects of SIFT on the reproducibility of structural connectome construction. Thirdly, we compare quantitative measures extracted from structural connectomes derived from streamlines tractography, with and without the application of SIFT, to published estimates drawn from post-mortem brain dissection. The combination of these sources of evidence demonstrates the important role the SIFT methodology has for the robust quantification of structural connectivity of the brain using diffusion MRI.

SIFT2: Enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography.

Diffusion MRI streamlines tractography allows for the investigation of the brain white matter pathways non-invasively. However a fundamental limitation of this technology is its non-quantitative nature, i.e. the density of reconstructed connections is not reflective of the density of underlying white matter fibres. As a solution to this problem, we have previously published the "spherical-deconvolution informed filtering of tractograms (SIFT)" method, which determines a subset of the streamlines reconstruction such that the streamlines densities throughout the white matter are as close as possible to fibre densities estimated using the spherical deconvolution diffusion model; this permits the use of streamline count as a valid biological marker of connection density. Particular aspects of its performance may have however limited its uptake in the diffusion MRI research community. Here we present an alternative to this method, entitled SIFT2, which provides a more logically direct and computationally efficient solution to the streamlines connectivity quantification problem: by determining an appropriate cross-sectional area multiplier for each streamline rather than removing streamlines altogether, biologically accurate measures of fibre connectivity are obtained whilst making use of the complete streamlines reconstruction.

Quantification of voxel-wise total fibre density: Investigating the problems associated with track-count mapping.

A biological parameter that would be valuable to be able to extract from diffusion MRI data is the local white matter axonal density. Track-density imaging (TDI) has been used as if it could provide such a measure; however, this has been the subject of controversy, primarily due to the fact that track-count quantitation is highly sensitive to tracking biases and errors. The spherical-deconvolution informed filtering of tractograms (SIFT) post-processing method was recently introduced to minimise tractography biases, and thus provides a more biologically meaningful measure that could be used in track-count mapping (i.e. TDI following SIFT). The TDI intensity following SIFT ideally corresponds to the orientational average of the fibre orientation distribution (FOD), which corresponds to the total Apparent Fibre Density (AFDtotal) within the AFD framework; in fact, AFDtotal provides a direct measure of local fibre density at native resolution that does not rely on fibre-tracking. In this study, we demonstrate problems associated with quantitative TDI investigations, which can be avoided by using SIFT processing or directly by using AFDtotal maps. We also characterise the intra- and inter-subject reproducibility of TDI maps (with and without SIFT pre-processing) and AFDtotal maps. It is shown that SIFT improves the quantitative characteristics of TDI, but is still vastly inferior to the properties of the AFDtotal parameter itself, because the latter does not require tracking. While standard TDI might be preferable in applications when high anatomical contrast is required, particularly when combined with super-resolution, for voxel-wise quantitation of total tract density (i.e. without tract orientation information) at native resolution, the total AFD maps are preferable to TDI or other related track-count maps. Regardless of the track-count measure, it should be noted that all of these voxel-averaged approaches discard important information that is retained in fibre-specific approaches such as AFD.

Fourier Tract Sampling (FouTS): A framework for improved inference of white matter tracts from diffusion MRI by explicitly modelling tract volume.

Diffusion MRI tractography algorithm development is increasingly moving towards global techniques to incorporate "downstream" information and conditional probabilities between neighbouring tracts. Such approaches also enable white matter to be represented more tangibly than the abstract lines generated by the most common approaches to fibre tracking. However, previously proposed algorithms still use fibre-like models of white matter corresponding to thin strands of white matter tracts rather than the tracts themselves, and therefore require many components for accurate representations, which leads to poorly constrained inverse problems. We propose a novel tract-based model of white matter, the 'Fourier tract', which is able to represent rich tract shapes with a relatively low number of parameters, and explicitly decouples the spatial extent of the modelled tract from its 'Apparent Connection Strength (ACS)'. The Fourier tract model is placed within a novel Bayesian framework, which relates the tract parameters directly to the observed signal, enabling a wide range of acquisition schemes to be used. The posterior distribution of the Bayesian framework is characterised via Markov-chain Monte-Carlo sampling to infer probable values of the ACS and spatial extent of the imaged white matter tracts, providing measures that can be directly applied to many research and clinical studies. The robustness of the proposed tractography algorithm is demonstrated on simulated basic tract configurations, such as curving, twisting, crossing and kissing tracts, and sections of more complex numerical phantoms. As an illustration of the approach in vivo, fibre tracking is performed on a central section of the brain in three subjects from 60 direction HARDI datasets.

Multi-tissue constrained spherical deconvolution for improved analysis of multi-shell diffusion MRI data.

Constrained spherical deconvolution (CSD) has become one of the most widely used methods to extract white matter (WM) fibre orientation information from diffusion-weighted MRI (DW-MRI) data, overcoming the crossing fibre limitations inherent in the diffusion tensor model. It is routinely used to obtain high quality fibre orientation distribution function (fODF) estimates and fibre tractograms and is increasingly used to obtain apparent fibre density (AFD) measures. Unfortunately, CSD typically only supports data acquired on a single shell in q-space. With multi-shell data becoming more and more prevalent, there is a growing need for CSD to fully support such data. Furthermore, CSD can only provide high quality fODF estimates in voxels containing WM only. In voxels containing other tissue types such as grey matter (GM) and cerebrospinal fluid (CSF), the WM response function may no longer be appropriate and spherical deconvolution produces unreliable, noisy fODF estimates. The aim of this study is to incorporate support for multi-shell data into the CSD approach as well as to exploit the unique b-value dependencies of the different tissue types to estimate a multi-tissue ODF. The resulting approach is dubbed multi-shell, multi-tissue CSD (MSMT-CSD) and is compared to the state-of-the-art single-shell, single-tissue CSD (SSST-CSD) approach. Using both simulations and real data, we show that MSMT-CSD can produce reliable WM/GM/CSF volume fraction maps, directly from the DW data, whereas SSST-CSD has a tendency to overestimate the WM volume in voxels containing GM and/or CSF. In addition, compared to SSST-CSD, MSMT-CSD can substantially increase the precision of the fODF fibre orientations and reduce the presence of spurious fODF peaks in voxels containing GM and/or CSF. Both effects translate into more reliable AFD measures and tractography results with MSMT-CSD compared to SSST-CSD.

Mapping somatosensory connectivity in adult mice using diffusion MRI tractography and super-resolution track density imaging.

In this study we combined ultra-high field diffusion MRI fiber tracking and super-resolution track density imaging (TDI) to map the relay locations and connectivity of the somatosensory pathway in paraformaldehyde fixed, C57Bl/6J mouse brains. Super-resolution TDI was used to achieve 20 µm isotropic resolution to inform the 3D topography of the relay locations including thalamic barreloids and brainstem barrelettes, not described previously using MRI methodology. TDI-guided mapping results for thalamo-cortical connectivity were consistent with thalamo-cortical projections labeled using virus mediated fluorescent protein expression. Trigemino-thalamic TDI connectivity maps were concordant with results obtained using anterograde dye tracing from brainstem to thalamus. Importantly, TDI mapping overcame the constraint of tissue distortion observed in mechanically sectioned tissue, enabling 3D reconstruction and long-range connectivity data. In conclusion, our results showed that diffusion micro-imaging at ultra-high field MRI revealed the stereotypical pattern of somatosensory connectivity and is a valuable tool to complement histologic methods, achieving 3D spatial preservation of whole brain networks for characterization in mouse models of human disease.

Dynamic changes in subplate and cortical plate microstructure at the onset of cortical folding in vivo.

Cortical gyrification takes place predominantly during the second to third trimester, alongside other fundamental developmental processes, such as the development of white matter connections, lamination of the cortex and formation of neural circuits. The mechanistic biology that drives the formation cortical folding patterns remains an open question in neuroscience. In our previous work, we modelled the in utero diffusion signal to quantify the maturation of microstructure in transient fetal compartments, identifying patterns of change in diffusion metrics that reflect critical neurobiological transitions occurring in the second to third trimester. In this work, we apply the same modelling approach to explore whether microstructural maturation of these compartments is correlated with the process of gyrification. We quantify the relationship between sulcal depth and tissue anisotropy within the cortical plate (CP) and underlying subplate (SP), key transient fetal compartments often implicated in mechanistic hypotheses about the onset of gyrification. Using in utero high angular resolution multi-shell diffusion-weighted imaging (HARDI) from the Developing Human Connectome Project (dHCP), our analysis reveals that the anisotropic, tissue component of the diffusion signal in the SP and CP decreases immediately prior to the formation of sulcal pits in the fetal brain. By back-projecting a map of folded brain regions onto the unfolded brain, we find evidence for cytoarchitectural differences between gyral and sulcal areas in the late second trimester, suggesting that regional variation in the microstructure of transient fetal compartments precedes, and thus may have a mechanistic function, in the onset of cortical folding in the developing human brain.

SIFT: Spherical-deconvolution informed filtering of tractograms.

Diffusion MRI allows the structural connectivity of the whole brain (the 'tractogram') to be estimated in vivo non-invasively using streamline tractography. The biological accuracy of these data sets is however limited by the inherent biases associated with the reconstruction method. Here we propose a method to retrospectively improve the accuracy of these reconstructions, by selectively filtering out streamlines from the tractogram in a manner that improves the fit between the streamline reconstruction and the underlying diffusion images. This filtering is guided by the results of spherical deconvolution of the diffusion signal, hence the acronym SIFT: spherical-deconvolution informed filtering of tractograms. Data sets processed by this algorithm show a marked reduction in known reconstruction biases, and improved biological plausibility. Emerging methods in diffusion MRI, particularly those that aim to characterise and compare the structural connectivity of the brain, should benefit from the improved accuracy of the reconstruction.

Track-weighted functional connectivity (TW-FC): a tool for characterizing the structural-functional connections in the brain.

MRI provides a powerful tool for studying the functional and structural connections in the brain non-invasively. The technique of functional connectivity (FC) exploits the intrinsic temporal correlations of slow spontaneous signal fluctuations to characterise brain functional networks. In addition, diffusion MRI fibre-tracking can be used to study the white matter structural connections. In recent years, there has been considerable interest in combining these two techniques to provide an overall structural-functional description of the brain. In this work we applied the recently proposed super-resolution track-weighted imaging (TWI) methodology to demonstrate how whole-brain fibre-tracking data can be combined with FC data to generate a track-weighted (TW) FC map of FC networks. The method was applied to data from 8 healthy volunteers, and illustrated with (i) FC networks obtained using a seeded connectivity-based analysis (seeding in the precuneus/posterior cingulate cortex, PCC, known to be part of the default mode network), and (ii) with FC networks generated using independent component analysis (in particular, the default mode, attention, visual, and sensory-motor networks). TW-FC maps showed high intensity in white matter structures connecting the nodes of the FC networks. For example, the cingulum bundles show the strongest TW-FC values in the PCC seeded-based analysis, due to their major role in the connection between medial frontal cortex and precuneus/posterior cingulate cortex; similarly the superior longitudinal fasciculus was well represented in the attention network, the optic radiations in the visual network, and the corticospinal tract and corpus callosum in the sensory-motor network. The TW-FC maps highlight the white matter connections associated with a given FC network, and their intensity in a given voxel reflects the functional connectivity of the part of the nodes of the network linked by the structural connections traversing that voxel. They therefore contain a different (and novel) image contrast from that of the images used to generate them. The results shown in this study illustrate the potential of the TW-FC approach for the fusion of structural and functional data into a single quantitative image. This technique could therefore have important applications in neuroscience and neurology, such as for voxel-based comparison studies.