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

Neuroimaging of structural and functional connectivity in preterm infants with intraventricular hemorrhage

Preterm infants with intraventricular hemorrhage (IVH) are known to have some of the worst neurodevelopmental outcomes in all of neonatal medicine, with a growing body of evidence relating these outcomes to underlying disruptions in brain structure and function. This review begins by summarizing state-of-the-art neuroimaging techniques delineating structural and functional connectivity (diffusion and resting state functional MRI) and their application in infants with IVH, including unique technical challenges and emerging methods. We then review studies of altered structural and functional connectivity, highlighting the role of IVH severity and location. We subsequently detail investigations linking structural and functional findings in infancy to later outcomes in early childhood. We conclude with future directions including methodologic considerations for prospective and potentially interventional studies designed to mitigate disruptions to underlying structural and functional connections and improve neurodevelopmental outcomes in this high-risk population.

Diffusion tensor estimation with transformer neural networks

Diffusion tensor imaging (DTI) is a widely used method for studying brain white matter development and degeneration. However, standard DTI estimation methods depend on a large number of high-quality measurements. This would require long scan times and can be particularly difficult to achieve with certain patient populations such as neonates. Here, we propose a method that can accurately estimate the diffusion tensor from only six diffusion-weighted measurements. Our method achieves this by learning to exploit the relationships between the diffusion signals and tensors in neighboring voxels. Our model is based on transformer networks, which represent the state of the art in modeling the relationship between signals in a sequence. In particular, our model consists of two such networks. The first network estimates the diffusion tensor based on the diffusion signals in a neighborhood of voxels. The second network provides more accurate tensor estimations by learning the relationships between the diffusion signals as well as the tensors estimated by the first network in neighboring voxels. Our experiments with three datasets show that our proposed method achieves highly accurate estimations of the diffusion tensor and is significantly superior to three competing methods. Estimations produced by our method with six diffusion-weighted measurements are comparable with those of standard estimation methods with 30-88 diffusion-weighted measurements. Hence, our method promises shorter scan times and more reliable assessment of brain white matter, particularly in non-cooperative patients such as neonates and infants.

Developmental pattern of association fibers and their interaction with associated cortical microstructures in 0-5-month-old infants

Association fibers connect the cortical regions and experience rapid development involving myelination and axonal growth during infancy. Yet, the spatiotemporal patterns of microstructural changes along these tracts, as well as the developmental interaction between the white matter (WM) tracts and the cortical gray matter (cGM) connected to them, are mostly unknown during infancy. In this study, we performed a diffusion MRI-based tractography and microstructure study in a cohort of 89 healthy preterm-born infants with gestational age at birth between 28.1∼36.4 weeks and postmenstrual age at scan between 39.9∼59.9 weeks. Results revealed that several C-shaped fibers, such as the arcuate fasciculus, cingulum, and uncinate fasciculus, demonstrated symmetrical along-tract profiles; and the horizontally oriented running fibers, including the inferior fronto-occipital fasciculus and the inferior longitudinal fasciculus, demonstrated an anterior-posterior developmental gradient. This study characterized the along-tract profiles using fixel-based analysis and revealed that the fiber cross-section (FC) of all five association fibers demonstrated a fluctuating increase with age, while the fiber density (FD) monotonically increase with age. NODDI was utilized to analyze the microstructural development of cGM and indicated cGM connected to the anterior end of the association fibers developed faster than that of the posterior end during 0-5 months. Notably, a mediation analysis was used to explore the relation between the development of WM and associated cGM, and demonstrated a partial mediation effect of FD in WM on the development of intracellular volume (ICV) in cGM and a full mediation effect of ICV on the growth of FD in most fibers, suggesting a predominant mediation of cGM on the WM development. Furthermore, for assessing whether those results were biased by prematurity, we compared preterm- and term-born neonates with matched scan age, gender, and multiple births from the developing human connectome project (dHCP) dataset to assess the effect of preterm-birth, and the results indicated a similar developmental pattern of the association fibers and their attached cGM. These findings presented a comprehensive picture of the major association fibers during early infancy and deciphered the developmental interaction between WM and cGM in this period.

Maternal free fatty acid concentration during pregnancy is associated with newborn hypothalamic microstructure in humans

OBJECTIVE

This study tested the hypothesis, in a prospective cohort study design, that maternal saturated free fatty acid (sFFA) concentration during pregnancy is prospectively associated with offspring (newborn) hypothalamic (HTH) microstructure and to explore the functional relevance of this association with respect to early-childhood body fat percentage (BF%).

METHODS

In N = 94 healthy newborns (born mean 39.3 [SD 1.5] weeks gestation), diffusion-weighted magnetic resonance imaging was performed shortly after birth (25.3 [12.5] postnatal days), and a subgroup (n = 37) underwent a dual-energy x-ray absorptiometry scan in early childhood (4.7 [SD 0.7] years). Maternal sFFA concentration during pregnancy was quantified in fasting blood samples via liquid chromatography-mass spectrometry. Infant HTH microstructural integrity was characterized using mean diffusivity (MD). Multiple linear regression was used to test the association between maternal sFFA and HTH MD, accounting for newborn sex, age at scan, mean white matter MD, and image quality. Multiple linear regression models also tested the association between HTH MD and early-childhood BF%, accounting for breastfeeding status.

RESULTS

Maternal sFFA during pregnancy accounted for 8.3% of the variation in newborn HTH MD (β-std = 0.25; p = 0.006). Furthermore, newborn HTH MD prospectively accounted for 15% of the variation in early-childhood BF% (β-std = 0.32; p = 0.019).

CONCLUSIONS

These findings suggest that maternal overnutrition during pregnancy may influence the development of the fetal hypothalamus, which, in turn, may have clinical relevance for childhood obesity risk.

The Developing Human Connectome Project Neonatal Data Release

The Developing Human Connectome Project has created a large open science resource which provides researchers with data for investigating typical and atypical brain development across the perinatal period. It has collected 1228 multimodal magnetic resonance imaging (MRI) brain datasets from 1173 fetal and/or neonatal participants, together with collateral demographic, clinical, family, neurocognitive and genomic data from 1173 participants, together with collateral demographic, clinical, family, neurocognitive and genomic data. All subjects were studied in utero and/or soon after birth on a single MRI scanner using specially developed scanning sequences which included novel motion-tolerant imaging methods. Imaging data are complemented by rich demographic, clinical, neurodevelopmental, and genomic information. The project is now releasing a large set of neonatal data; fetal data will be described and released separately. This release includes scans from 783 infants of whom: 583 were healthy infants born at term; as well as preterm infants; and infants at high risk of atypical neurocognitive development. Many infants were imaged more than once to provide longitudinal data, and the total number of datasets being released is 887. We now describe the dHCP image acquisition and processing protocols, summarize the available imaging and collateral data, and provide information on how the data can be accessed.

Predicting age and clinical risk from the neonatal connectome

The development of perinatal brain connectivity underpins motor, cognitive and behavioural abilities in later life. Diffusion MRI allows the characterisation of subtle inter-individual differences in structural brain connectivity. Individual brain connectivity maps (connectomes) are by nature high in dimensionality and complex to interpret. Machine learning methods are a powerful tool to uncover properties of the connectome which are not readily visible and can give us clues as to how and why individual developmental trajectories differ. In this manuscript we used Deep Neural Networks and Random Forests to predict demographic and neurodevelopmental characteristics from neonatal structural connectomes in a large sample of babies (n = 524) from the developing Human Connectome Project. We achieved an accurate prediction of post menstrual age (PMA) at scan in term-born infants (mean absolute error (MAE) = 0.72 weeks, r = 0.83 and p<0.001). We also achieved good accuracy when predicting gestational age at birth in a cohort of term and preterm babies scanned at term equivalent age (MAE = 2.21 weeks, r = 0.82, p<0.001). We subsequently used sensitivity analysis to obtain feature relevance from our prediction models, with the most important connections for prediction of PMA and GA found to predominantly involve frontal and temporal regions, thalami, and basal ganglia. From our models of PMA at scan for infants born at term, we computed a brain maturation index (predicted age minus actual age) of individual preterm neonates and found a significant correlation between this index and motor outcome at 18 months corrected age. Our results demonstrate the applicability of machine learning techniques in analyses of the neonatal connectome and suggest that a neural substrate of brain maturation with implications for future neurodevelopment is detectable at term equivalent age from the neonatal connectome.

Through-Plane Super-Resolution With Autoencoders in Diffusion Magnetic Resonance Imaging of the Developing Human Brain

Fetal brain diffusion magnetic resonance images (MRI) are often acquired with a lower through-plane than in-plane resolution. This anisotropy is often overcome by classical upsampling methods such as linear or cubic interpolation. In this work, we employ an unsupervised learning algorithm using an autoencoder neural network for single-image through-plane super-resolution by leveraging a large amount of data. Our framework, which can also be used for slice outliers replacement, overperformed conventional interpolations quantitatively and qualitatively on pre-term newborns of the developing Human Connectome Project. The evaluation was performed on both the original diffusion-weighted signal and the estimated diffusion tensor maps. A byproduct of our autoencoder was its ability to act as a denoiser. The network was able to generalize fetal data with different levels of motions and we qualitatively showed its consistency, hence supporting the relevance of pre-term datasets to improve the processing of fetal brain images.

An ode to fetal, infant, and toddler neuroimaging: Chronicling early clinical to research applications with MRI, and an introduction to an academic society connecting the field

Fetal, infant, and toddler neuroimaging is commonly thought of as a development of modern times (last two decades). Yet, this field mobilized shortly after the discovery and implementation of MRI technology. Here, we provide a review of the parallel advancements in the fields of fetal, infant, and toddler neuroimaging, noting the shifts from clinical to research use, and the ongoing challenges in this fast-growing field. We chronicle the pioneering science of fetal, infant, and toddler neuroimaging, highlighting the early studies that set the stage for modern advances in imaging during this developmental period, and the large-scale multi-site efforts which ultimately led to the explosion of interest in the field today. Lastly, we consider the growing pains of the community and the need for an academic society that bridges expertise in developmental neuroscience, clinical science, as well as computational and biomedical engineering, to ensure special consideration of the vulnerable mother-offspring dyad (especially during pregnancy), data quality, and image processing tools that are created, rather than adapted, for the young brain.

Fiber tracing and microstructural characterization among audiovisual integration brain regions in neonates compared with young adults

Audiovisual integration has been related with cognitive-processing and behavioral advantages, as well as with various socio-cognitive disorders. While some studies have identified brain regions instantiating this ability shortly after birth, little is known about the structural pathways connecting them. The goal of the present study was to reconstruct fiber tracts linking AVI regions in the newborn in-vivo brain and assess their adult-likeness by comparing them with analogous fiber tracts of young adults. We performed probabilistic tractography and compared connective probabilities between a sample of term-born neonates (N = 311; the Developing Human Connectome Project (dHCP, http://www.developingconnectome.org) and young adults (N = 311 The Human Connectome Project; https://www.humanconnectome.org/) by means of a classification algorithm. Furthermore, we computed Dice coefficients to assess between-group spatial similarity of the reconstructed fibers and used diffusion metrics to characterize neonates' AVI brain network in terms of microstructural properties, interhemispheric differences and the association with perinatal covariates and biological sex. Overall, our results indicate that the AVI fiber bundles were successfully reconstructed in a vast majority of neonates, similarly to adults. Connective probability distributional similarities and spatial overlaps of AVI fibers between the two groups differed across the reconstructed fibers. There was a rank-order correspondence of the fibers' connective strengths across the groups. Additionally, the study revealed patterns of diffusion metrics in line with early white matter developmental trajectories and a developmental advantage for females. Altogether, these findings deliver evidence of meaningful structural connections among AVI regions in the newborn in-vivo brain.

Medical Image Segmentation Using Transformer Networks

Deep learning models represent the state of the art in medical image segmentation. Most of these models are fully-convolutional networks (FCNs), namely each layer processes the output of the preceding layer with convolution operations. The convolution operation enjoys several important properties such as sparse interactions, parameter sharing, and translation equivariance. Because of these properties, FCNs possess a strong and useful inductive bias for image modeling and analysis. However, they also have certain important shortcomings, such as performing a fixed and pre-determined operation on a test image regardless of its content and difficulty in modeling long-range interactions. In this work we show that a different deep neural network architecture, based entirely on self-attention between neighboring image patches and without any convolution operations, can achieve more accurate segmentations than FCNs. Our proposed model is based directly on the transformer network architecture. Given a 3D image block, our network divides it into non-overlapping 3D patches and computes a 1D embedding for each patch. The network predicts the segmentation map for the block based on the self-attention between these patch embeddings. Furthermore, in order to address the common problem of scarcity of labeled medical images, we propose methods for pre-training this model on large corpora of unlabeled images. Our experiments show that the proposed model can achieve segmentation accuracies that are better than several state of the art FCN architectures on two datasets. Our proposed network can be trained using only tens of labeled images. Moreover, with the proposed pre-training strategies, our network outperforms FCNs when labeled training data is small.

White matter myelination during early infancy is linked to spatial gradients and myelin content at birth

Development of myelin, a fatty sheath that insulates nerve fibers, is critical for brain function. Myelination during infancy has been studied with histology, but postmortem data cannot evaluate the longitudinal trajectory of white matter development. Here, we obtained longitudinal diffusion MRI and quantitative MRI measures of longitudinal relaxation rate (R1) of white matter in 0, 3 and 6 months-old human infants, and developed an automated method to identify white matter bundles and quantify their properties in each infant's brain. We find that R1 increases from newborns to 6-months-olds in all bundles. R1 development is nonuniform: there is faster development in white matter that is less mature in newborns, and development rate increases along inferior-to-superior as well as anterior-to-posterior spatial gradients. As R1 is linearly related to myelin fraction in white matter bundles, these findings open new avenues to elucidate typical and atypical white matter myelination in early infancy.

Dear reviewers: Responses to common reviewer critiques about infant neuroimaging studies

The field of adult neuroimaging relies on well-established principles in research design, imaging sequences, processing pipelines, as well as safety and data collection protocols. The field of infant magnetic resonance imaging, by comparison, is a young field with tremendous scientific potential but continuously evolving standards. The present article aims to initiate a constructive dialog between researchers who grapple with the challenges and inherent limitations of a nascent field and reviewers who evaluate their work. We address 20 questions that researchers commonly receive from research ethics boards, grant, and manuscript reviewers related to infant neuroimaging data collection, safety protocols, study planning, imaging sequences, decisions related to software and hardware, and data processing and sharing, while acknowledging both the accomplishments of the field and areas of much needed future advancements. This article reflects the cumulative knowledge of experts in the FIT’NG community and can act as a resource for both researchers and reviewers alike seeking a deeper understanding of the standards and tradeoffs involved in infant neuroimaging.

•

The field of infant MRI is young with evolving standards.

•

We address 20 questions that researchers commonly receive reviewers.

•

These come from research ethics boards, grant, and manuscript reviewers.

•

This article reflects the cumulative knowledge of experts in the FIT’NG community.

What's New and What's Next in Diffusion MRI Preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B₁ bias fields, and spatial normalization. The focus will be on “what’s new” since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on “Mapping the Connectome” in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on “what’s next” in dMRI preprocessing.

The Human Connectome Project: A retrospective

The Human Connectome Project (HCP) was launched in 2010 as an ambitious effort to accelerate advances in human neuroimaging, particularly for measures of brain connectivity; apply these advances to study a large number of healthy young adults; and freely share the data and tools with the scientific community. NIH awarded grants to two consortia; this retrospective focuses on the "WU-Minn-Ox" HCP consortium centered at Washington University, the University of Minnesota, and University of Oxford. In just over 6 years, the WU-Minn-Ox consortium succeeded in its core objectives by: 1) improving MR scanner hardware, pulse sequence design, and image reconstruction methods, 2) acquiring and analyzing multimodal MRI and MEG data of unprecedented quality together with behavioral measures from more than 1100 HCP participants, and 3) freely sharing the data (via the ConnectomeDB database) and associated analysis and visualization tools. To date, more than 27 Petabytes of data have been shared, and 1538 papers acknowledging HCP data use have been published. The "HCP-style" neuroimaging paradigm has emerged as a set of best-practice strategies for optimizing data acquisition and analysis. This article reviews the history of the HCP, including comments on key events and decisions associated with major project components. We discuss several scientific advances using HCP data, including improved cortical parcellations, analyses of connectivity based on functional and diffusion MRI, and analyses of brain-behavior relationships. We also touch upon our efforts to develop and share a variety of associated data processing and analysis tools along with detailed documentation, tutorials, and an educational course to train the next generation of neuroimagers. We conclude with a look forward at opportunities and challenges facing the human neuroimaging field from the perspective of the HCP consortium.

Association between brain structural network efficiency at term-equivalent age and early development of cerebral palsy in very preterm infants

Very preterm infants (born at less than 32 weeks gestational age) are at high risk for serious motor impairments, including cerebral palsy (CP). The brain network changes that antecede the early development of CP in infants are not well characterized, and a better understanding may suggest new strategies for risk-stratification at term, which could lead to earlier access to therapies. Graph theoretical methods applied to diffusion MRI-derived brain connectomes may help quantify the organization and information transfer capacity of the preterm brain with greater nuance than overt structural or regional microstructural changes. Our aim was to shed light on the pathophysiology of early CP development, before the occurrence of early intervention therapies and other environmental confounders, to help identify the best early biomarkers of CP risk in VPT infants. In a cohort of 395 very preterm infants, we extracted cortical morphometrics and brain volumes from structural MRI and also applied graph theoretical methods to diffusion MRI connectomes, both acquired at term-equivalent age. Metrics from graph network analysis, especially global efficiency, strength values of the major sensorimotor tracts, and local efficiency of the motor nodes and novel non-motor regions were strongly inversely related to early CP diagnosis. These measures remained significantly associated with CP after correction for common risk factors of motor development, suggesting that metrics of brain network efficiency at term may be sensitive biomarkers for early CP detection. We demonstrate for the first time that in VPT infants, early CP diagnosis is anteceded by decreased brain network segregation in numerous nodes, including motor regions commonly-associated with CP and also novel regions that may partially explain the high rate of cognitive impairments concomitant with CP diagnosis. These advanced MRI biomarkers may help identify the highest risk infants by term-equivalent age, facilitating earlier interventions that are informed by early pathophysiological changes.

Motor function and white matter connectivity in children cooled for neonatal encephalopathy

Therapeutic hypothermia reduces the incidence of severe motor disability, such as cerebral palsy, following neonatal hypoxic-ischaemic encephalopathy. However, cooled children without cerebral palsy at school-age demonstrate motor deficits and altered white matter connectivity. In this study, we used diffusion-weighted imaging to investigate the relationship between white matter connectivity and motor performance, measured using the Movement Assessment Battery for Children-2, in children aged 6-8 years treated with therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy at birth, who did not develop cerebral palsy (cases), and matched typically developing controls. Correlations between total motor scores and diffusion properties in major white matter tracts were assessed in 33 cases and 36 controls. In cases, significant correlations (FDR-corrected P < 0.05) were found in the anterior thalamic radiation bilaterally (left: r = 0.513; right: r = 0.488), the cingulate gyrus part of the left cingulum (r = 0.588), the hippocampal part of the left cingulum (r = 0.541), and the inferior fronto-occipital fasciculus bilaterally (left: r = 0.445; right: r = 0.494). No significant correlations were found in controls. We then constructed structural connectivity networks, for 22 cases and 32 controls, in which nodes represent brain regions and edges were determined by probabilistic tractography and weighted by fractional anisotropy. Analysis of whole-brain network metrics revealed correlations (FDR-corrected P < 0.05), in cases, between total motor scores and average node strength (r = 0.571), local efficiency (r = 0.664), global efficiency (r = 0.677), clustering coefficient (r = 0.608), and characteristic path length (r = -0.652). No significant correlations were found in controls. We then investigated edge-level association with motor function using the network-based statistic. This revealed subnetworks which exhibited group differences in the association between motor outcome and edge weights, for total motor scores (P = 0.0109) as well as for balance (P = 0.0245) and manual dexterity (P = 0.0233) domain scores. All three of these subnetworks comprised numerous frontal lobe regions known to be associated with motor function, including the superior frontal gyrus and middle frontal gyrus. The subnetwork associated with total motor scores was highly left-lateralised. These findings demonstrate an association between impaired motor function and brain organisation in school-age children treated with therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy.

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.

An Age-Specific Atlas for Delineation of White Matter Pathways in Children Aged 6-8 Years

Introduction

Diffusion MRI allows non-invasive assessment of white matter connectivity in typical development and of changes due to brain injury or pathology. Probabilistic white matter atlases allow diffusion metrics to be measured in specific white matter pathways, and are a critical component in spatial normalisation for group analysis. However, given the known developmental changes in white matter it may be sub-optimal to use an adult template when assessing data acquired from children.

Methods

By averaging subject-specific fibre bundles from 28 children aged from 6 to 8 years, we created an age-specific probabilistic white matter atlas for 12 major white matter tracts. Using both the newly developed and Johns Hopkins adult atlases, we compared the atlas to subject-specific fibre bundles in two independent validation cohorts, assessing accuracy in terms of volumetric overlap and measured diffusion metrics.

Results

Our age-specific atlas gave better overall performance than the adult atlas, achieving higher volumetric overlap with subject-specific fibre tracking and higher correlation of FA measurements with those measured from subject-specific fibre bundles. Specifically, estimates of FA values for cortico-spinal tract, uncinate fasciculus, forceps minor, cingulate gyrus part of the cingulum and anterior thalamic radiation were all significantly more accurate when estimated with an age-specific atlas.

Discussion

The age-specific atlas allows delineation of white matter tracts in children aged 6-8 years, without the need for tractography, more accurately than when normalising to an adult atlas. To our knowledge, this is the first publicly available probabilistic atlas of white matter tracts for this age group.

Diffusion MRI of the infant brain reveals unique asymmetry patterns during the first-half-year of development

The human brain demonstrates anatomical and functional lateralization/asymmetry between the left and right hemispheres, and such asymmetry is known to start from the early age of life. However, how the asymmetry changes with brain development during infancy remained unknown. In this study, we aimed to systematically investigate the spatiotemporal pattern of brain asymmetry in healthy preterm-born infants during the first-half-year of development, using high angular resolution diffusion MRI. Sixty-five healthy preterm-born infants (gestational age between 25.3-36.6 weeks) were scanned with postmenstrual age (PMA) ranging from term-equivalent age (TEA) to 6-months. At the regional level, we performed a region-of-interest-based analysis by segmenting the brain into 63 symmetrical pairs of regions, based on which the laterality index was assessed and correlated with PMA. At the voxel level, we performed a fixel-based analysis of each fiber component between the native and left-right flipped data, separately in TEA-1 month, 1-3 months, and 3-6 months groups. The infant brains demonstrated extensive regions with structural asymmetry during their first half-of-year of life. A distinct central-peripheral asymmetry pattern was observed in mean diffusivity, namely, leftward lateralization in the neocortex and rightward asymmetry in the deep brain regions. Besides, the posterior brain demonstrated a higher lateralization index compared with the anterior brain in all metrics, which is congruent with the brain developmental pattern from caudal to rostral. Regionally, language processing regions showed a rightward asymmetry, while visuospatial processing regions exhibited leftward lateralization in fractional anisotropy, fibre density, and fibre cross-section measurements, and most white matter regions were lateralized to the left in these measurements. The laterality index of several regions (12 out 63) demonstrated significant developmental changes in mean diffusivity. At the fixel level, the fiber cross-section of inferior fronto-occipital fasciculus showed significant leftward asymmetry and the extent of asymmetry increased with PMA. In summary, the results revealed unique spatiotemporal patterns of macro- and micro-structural asymmetry in early life, which dynamically changed with age. These findings may contribute to the understanding of brain development during infancy.

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.

Multi-Channel 4D Parametrized Atlas of Macro- and Microstructural Neonatal Brain Development

Structural (also known as anatomical) and diffusion MRI provide complimentary anatomical and microstructural characterization of early brain maturation. However, the existing models of the developing brain in time include only either structural or diffusion MRI channels. Furthermore, there is a lack of tools for combined analysis of structural and diffusion MRI in the same reference space. In this work, we propose a methodology to generate a multi-channel (MC) continuous spatio-temporal parametrized atlas of the brain development that combines multiple MRI-derived parameters in the same anatomical space during 37-44 weeks of postmenstrual age range. We co-align structural and diffusion MRI of 170 normal term subjects from the developing Human Connectomme Project using MC registration driven by both T2-weighted and orientation distribution functions channels and fit the Gompertz model to the signals and spatial transformations in time. The resulting atlas consists of 14 spatio-temporal microstructural indices and two parcellation maps delineating white matter tracts and neonatal transient structures. In order to demonstrate applicability of the atlas for quantitative region-specific studies, a comparison analysis of 140 term and 40 preterm subjects scanned at the term-equivalent age is performed using different MRI-derived microstructural indices in the atlas reference space for multiple white matter regions, including the transient compartments. The atlas and software will be available after publication of the article.

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.

Multi-scale structural rich-club organization of the brain in full-term newborns: a combined DWI and fMRI study

Objective.Our understanding of early brain development is limited due to rapid changes in white matter pathways after birth. In this study, we introduced a multi-scale cross-modal approach to investigate the rich club (RC) organization and topology of the structural brain networks in 40 healthy neonates using diffusion-weighted imaging and resting-state fMRI data.Approach.A group independent component analysis was first performed to identify eight resting state networks (RSNs) used as functional modules. A groupwise whole-brain functional parcellation was also performed at five scales comprising 100-900 parcels. The distribution of RC nodes was then investigated within and between the RSNs. We further assessed the distribution of short and long-range RC, feeder and local connections across different parcellation scales.Main results.Sharing the scale-free characteristic of small-worldness, the neonatal structural brain networks exhibited an RC organization at different nodal scales (NSs). The subcortical, sensory-motor and default mode networks were found to be strongly involved in the RC organization of the structural brain networks, especially in the zones where the RSNs overlapped, with an average cross-scale proportion of 45.9%, 28.5% and 10.5%, respectively. A large proportion of the connector hubs were found to be RC members for the coarsest (73%) to finest (92%) NSs. Our results revealed a prominent involvement of cortico-subcortical and cortico-cerebellar white matter pathways in the RC organization of the neonatal brain. Regardless of the NS, the majority (more than 65.2%) of the inter-RSN connections were long distance RC or feeder with an average physical connection of 105.5 and 97.4 mm, respectively. Several key RC regions were identified, including the insula and cingulate gyri, middle and superior temporal gyri, hippocampus and parahippocampus, fusiform gyrus, precuneus, superior frontal and precentral gyri, calcarine fissure and lingual gyrus.Significance.Our results emphasize the importance of the multi-scale connectivity analysis in assessing the cross-scale reproducibility of the connectivity results concerning the global and local topological properties of the brain networks. Our findings may improve our understanding of the early brain development.

Functional and diffusion MRI reveal the neurophysiological basis of neonates' noxious-stimulus evoked brain activity

Understanding the neurophysiology underlying neonatal responses to noxious stimulation is central to improving early life pain management. In this neonatal multimodal MRI study, we use resting-state and diffusion MRI to investigate inter-individual variability in noxious-stimulus evoked brain activity. We observe that cerebral haemodynamic responses to experimental noxious stimulation can be predicted from separately acquired resting-state brain activity (n = 18). Applying this prediction model to independent Developing Human Connectome Project data (n = 215), we identify negative associations between predicted noxious-stimulus evoked responses and white matter mean diffusivity. These associations are subsequently confirmed in the original noxious stimulation paradigm dataset, validating the prediction model. Here, we observe that noxious-stimulus evoked brain activity in healthy neonates is coupled to resting-state activity and white matter microstructure, that neural features can be used to predict responses to noxious stimulation, and that the dHCP dataset could be utilised for future exploratory research of early life pain system neurophysiology.

Aquaporin 4 expression in the hippocampus in sudden infant death syndrome and sudden unexplained death in childhood

Aquaporin 4 (AQP4) is the main membrane water channel in the brain involved in regulating water homeostasis. The water distribution in neural tissue is often dysregulated after hypoxic neural injury. Previous research has indicated that victims of sudden infant death syndrome (SIDS) and sudden unexplained death in childhood (SUDC) have an underlying brain dysfunction that impairs their critical arousal response to hypoxic stress during sleep. The aim of this study was to determine the expression levels of AQP4 in the hippocampus in SIDS/SUDC cases and controls, and compare the findings with AQP4 genotypes that previously have been shown to be associated with SIDS. Immunochemical staining and morphometry were used to evaluate the density of AQP4-positive astrocytes in 30 SIDS/SUDC cases and 26 controls. AQP4-positive cells were counted in grids covering three layers in the hippocampus, which revealed that their count in any of the layers did not differ significantly between cases and controls. A decline in AQP4 expression was observed for infants older than 12 weeks. The AQP4 expression was lower in infants and children with the rs2075575 CT/TT genotype than in those with the CC genotype. This study indicates that AQP4 expression may be influenced by both age and genotype in infants. The role of AQP4 in the pathogenesis of SIDS remains to be elucidated.

MRI of the Neonatal Brain: A Review of Methodological Challenges and Neuroscientific Advances

In recent years, exploration of the developing brain has become a major focus for researchers and clinicians in an attempt to understand what allows children to acquire amazing and unique abilities, as well as the impact of early disruptions (eg, prematurity, neonatal insults) that can lead to a wide range of neurodevelopmental disorders. Noninvasive neuroimaging methods such as MRI are essential to establish links between the brain and behavioral changes in newborns and infants. In this review article, we aim to highlight recent and representative studies using the various techniques available: anatomical MRI, quantitative MRI (relaxometry, diffusion MRI), multiparametric approaches, and functional MRI. Today, protocols use 1.5 or 3T MRI scanners, and specialized methodologies have been put in place for data acquisition and processing to address the methodological challenges specific to this population, such as sensitivity to motion. MR sequences must be adapted to the brains of newborns and infants to obtain relevant good soft-tissue contrast, given the small size of the cerebral structures and the incomplete maturation of tissues. The use of age-specific image postprocessing tools is also essential, as signal and contrast differ from the adult brain. Appropriate methodologies then make it possible to explore multiple neurodevelopmental mechanisms in a precise way, and assess changes with age or differences between groups of subjects, particularly through large-scale projects. Although MRI measurements only indirectly reflect the complex series of dynamic processes observed throughout development at the molecular and cellular levels, this technique can provide information on brain morphology, structural connectivity, microstructural properties of gray and white matter, and on the functional architecture. Finally, MRI measures related to clinical, behavioral, and electrophysiological markers have a key role to play from a diagnostic and prognostic perspective in the implementation of early interventions to avoid long-term disabilities in children. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY STAGE: 1.

Improved reproducibility of diffusion MRI of the human brain with a four-way blip-up and down phase-encoding acquisition approach

PURPOSE

To assess the effects of blip-up and -down echo planar imaging (EPI) acquisition designs, with different choices of phase-encoding directions (PEDs) on the reproducibility of diffusion MRI (dMRI)-derived metrics in the human brain.

METHODS

Diffusion MRI data in seven subjects were acquired five times, each with five different protocols. The base design included 64 diffusion directions acquired with anterior-posterior (AP) PED, the first and second protocols added reverse phase-encoded b = 0 s / mm 2 posterior-anterior (PA) PED images. The third one included 32 directions all with PED acquisitions with opposite polarity (AP and PA). The fourth protocol, also with 32 unique directions used four PEDs (AP, PA, right-left (RL), and left-right (LR)). The scan time was virtually identical for all protocols. The variability of diffusion MRI metrics for each subject and each protocol was computed across the different sessions.

RESULTS

The highest reproducibility for all dMRI metrics was obtained with protocol four (AP/PA-RL/LR, ie, four-way PED). Protocols that used only b = 0 s / mm 2 for distortion correction, which are the most widely used designs, had the lowest reproducibility.

CONCLUSIONS

An acquisition design with four PEDs, including all DWIs in addition to b = 0 s / mm 2 images should be used to achieve high reproducibility in diffusion MRI studies.

Minimum TR radiofrequency-pulse design for rapid gradient echo sequences

PURPOSE

A framework to design radiofrequency (RF) pulses specifically to minimize the TR of gradient echo sequences is presented, subject to hardware and physiological constraints.

METHODS

Single-band and multiband (MB) RF pulses can be reduced in duration using variable-rate selective excitation (VERSE) VERSE for a range of flip angles; however, minimum-duration pulses do not guarantee minimum TR because these can lead to a high specific absorption rate (SAR). The optimal RF pulse is found by meeting spatial encoding, peripheral nerve stimulation (PNS) and SAR constraints. A TR reduction for a range of designs is achieved and an application of this in an MB cardiac balanced steady-state free-precession (bSSFP) experiment is presented. Gradient imperfections and their imaging effects are also considered.

RESULTS

Sequence TR with low-time bandwidth product (TBP) pulses, as used in bSSFP, was reduced up to 14%, and the TR when using high TBP pulses, as used in slab-selective imaging, was reduced by up to 72%. A breath-hold cardiac exam was reduced by 46% using both MB and the TR-optimal framework. The importance of RF-based correction of gradient imperfections is demonstrated. PNS was not a practical limitation.

CONCLUSION

The TR-optimal framework designs RF pulses for a range of pulse parameters, specifically to minimize sequence TR.

Disrupted brain connectivity in children treated with therapeutic hypothermia for neonatal encephalopathy

Therapeutic hypothermia following neonatal encephalopathy due to birth asphyxia reduces death and cerebral palsy. However, school-age children without cerebral palsy treated with therapeutic hypothermia for neonatal encephalopathy still have reduced performance on cognitive and motor tests, attention difficulties, slower reaction times and reduced visuo-spatial processing abilities compared to typically developing controls. We acquired diffusion-weighted imaging data from school-age children without cerebral palsy treated with therapeutic hypothermia for neonatal encephalopathy at birth, and a matched control group. Voxelwise analysis (33 cases, 36 controls) confirmed reduced fractional anisotropy in widespread areas of white matter in cases, particularly in the fornix, corpus callosum, anterior and posterior limbs of the internal capsule bilaterally and cingulum bilaterally. In structural brain networks constructed using probabilistic tractography (22 cases, 32 controls), graph-theoretic measures of strength, local and global efficiency, clustering coefficient and characteristic path length were found to correlate with IQ in cases but not controls. Network-based statistic analysis implicated brain regions involved in visuo-spatial processing and attention, aligning with previous behavioural findings. These included the precuneus, thalamus, left superior parietal gyrus and left inferior temporal gyrus. Our findings demonstrate that, despite the manifest successes of therapeutic hypothermia, brain development is impaired in these children.

Imaging the rapidly developing brain: Current challenges for MRI studies in the first five years of life

Rapid and widespread changes in brain anatomy and physiology in the first five years of life present substantial challenges for developmental structural, functional, and diffusion MRI studies. One persistent challenge is that methods best suited to earlier developmental stages are suboptimal for later stages, which engenders a trade-off between using different, but age-appropriate, methods for different developmental stages or identical methods across stages. Both options have potential benefits, but also biases, as pipelines for each developmental stage can be matched on methods or the age-appropriateness of methods, but not both. This review describes the data acquisition, processing, and analysis challenges that introduce these potential biases and attempts to elucidate decisions and make recommendations that would optimize developmental comparisons.

Scattered slice SHARD reconstruction for motion correction in multi-shell diffusion MRI

•

Subject motion in dMRI leads to a set of scattered slices with unique contrast.

•

We introduce a slice-to-volume reconstruction framework for multi-shell HARDI data

•

Based on a data-driven representation as spherical harmonics and radial decomposition (SHARD).

•

The method is evaluated in test-retest scans and in the neonatal dHCP cohort.

•

Results show robust reconstruction in severely motion-corrupted scans.

Diffusion MRI offers a unique probe into neural microstructure and connectivity in the developing brain. However, analysis of neonatal brain imaging data is complicated by inevitable subject motion, leading to a series of scattered slices that need to be aligned within and across diffusion-weighted contrasts. Here, we develop a reconstruction method for scattered slice multi-shell high angular resolution diffusion imaging (HARDI) data, jointly estimating an uncorrupted data representation and motion parameters at the slice or multiband excitation level. The reconstruction relies on data-driven representation of multi-shell HARDI data using a bespoke spherical harmonics and radial decomposition (SHARD), which avoids imposing model assumptions, thus facilitating to compare various microstructure imaging methods in the reconstructed output. Furthermore, the proposed framework integrates slice-level outlier rejection, distortion correction, and slice profile correction. We evaluate the method in the neonatal cohort of the developing Human Connectome Project (650 scans). Validation experiments demonstrate accurate slice-level motion correction across the age range and across the range of motion in the population. Results in the neonatal data show successful reconstruction even in severely motion-corrupted subjects. In addition, we illustrate how local tissue modelling can extract advanced microstructure features such as orientation distribution functions from the motion-corrected reconstructions.

Preterm birth leads to impaired rich-club organization and fronto-paralimbic/limbic structural connectivity in newborns

Prematurity disrupts brain development during a critical period of brain growth and organization and is known to be associated with an increased risk of neurodevelopmental impairments. Investigating whole-brain structural connectivity alterations accompanying preterm birth may provide a better comprehension of the neurobiological mechanisms related to the later neurocognitive deficits observed in this population. Using a connectome approach, we aimed to study the impact of prematurity on neonatal whole-brain structural network organization at term-equivalent age. In this cohort study, twenty-four very preterm infants at term-equivalent age (VPT-TEA) and fourteen full-term (FT) newborns underwent a brain MRI exam at term age, comprising T2-weighted imaging and diffusion MRI, used to reconstruct brain connectomes by applying probabilistic constrained spherical deconvolution whole-brain tractography. The topological properties of brain networks were quantified through a graph-theoretical approach. Furthermore, edge-wise connectivity strength was compared between groups. Overall, VPT-TEA infants' brain networks evidenced increased segregation and decreased integration capacity, revealed by an increased clustering coefficient, increased modularity, increased characteristic path length, decreased global efficiency and diminished rich-club coefficient. Furthermore, in comparison to FT, VPT-TEA infants had decreased connectivity strength in various cortico-cortical, cortico-subcortical and intra-subcortical networks, the majority of them being intra-hemispheric fronto-paralimbic and fronto-limbic. Inter-hemispheric connectivity was also decreased in VPT-TEA infants, namely through connections linking to the left precuneus or left dorsal cingulate gyrus - two regions that were found to be hubs in FT but not in VPT-TEA infants. Moreover, posterior regions from Default-Mode-Network (DMN), namely precuneus and posterior cingulate gyrus, had decreased structural connectivity in VPT-TEA group. Our finding that VPT-TEA infants' brain networks displayed increased modularity, weakened rich-club connectivity and diminished global efficiency compared to FT infants suggests a delayed transition from a local architecture, focused on short-range connections, to a more distributed architecture with efficient long-range connections in those infants. The disruption of connectivity in fronto-paralimbic/limbic and posterior DMN regions might underlie the behavioral and social cognition difficulties previously reported in the preterm population.

Test-retest reliability and long-term stability of three-tissue constrained spherical deconvolution methods for analyzing diffusion MRI data

PURPOSE

Several recent studies have used a three-tissue constrained spherical deconvolution pipeline to obtain quantitative metrics of brain tissue microstructure from diffusion-weighted MRI data. The three tissue compartments, consisting of white matter, gray matter, and CSF-like (free water) signals, are potentially useful in the evaluation of brain microstructure in a range of pathologies. However, the reliability and long-term stability of these metrics have not yet been evaluated.

METHODS

This study examined estimates of whole-brain microstructure for the three tissue compartments, in three separate test-retest cohorts. Each cohort had different lengths of time between baseline and retest, ranging from within the same scanning session in the shortest interval to 3 months in the longest interval. Each cohort was also collected with different acquisition parameters.

RESULTS

The CSF-like compartment displayed the greatest reliability across all cohorts, with intraclass correlation coefficient (ICC) values being above 0.95 in each cohort. White matter-like and gray matter-like compartments both demonstrated very high reliability in the immediate cohort (both ICC > 0.90); however, this declined in the 3-month interval cohort to both compartments having ICC > 0.80. Regional CSF-like signal fraction was examined in bilateral hippocampus and had an ICC > 0.80 in each cohort.

CONCLUSION

The three-tissue constrained spherical deconvolution techniques provide reliable and stable estimates of tissue-microstructure composition, up to 3 months longitudinally in a control population. This forms an important basis for further investigations using three-tissue constrained spherical deconvolution techniques to track changes in microstructure across a variety of brain pathologies.

Infant FreeSurfer: An automated segmentation and surface extraction pipeline for T1-weighted neuroimaging data of infants 0-2 years

The development of automated tools for brain morphometric analysis in infants has lagged significantly behind analogous tools for adults. This gap reflects the greater challenges in this domain due to: 1) a smaller-scaled region of interest, 2) increased motion corruption, 3) regional changes in geometry due to heterochronous growth, and 4) regional variations in contrast properties corresponding to ongoing myelination and other maturation processes. Nevertheless, there is a great need for automated image-processing tools to quantify differences between infant groups and other individuals, because aberrant cortical morphologic measurements (including volume, thickness, surface area, and curvature) have been associated with neuropsychiatric, neurologic, and developmental disorders in children. In this paper we present an automated segmentation and surface extraction pipeline designed to accommodate clinical MRI studies of infant brains in a population 0-2 year-olds. The algorithm relies on a single channel of T1-weighted MR images to achieve automated segmentation of cortical and subcortical brain areas, producing volumes of subcortical structures and surface models of the cerebral cortex. We evaluated the algorithm both qualitatively and quantitatively using manually labeled datasets, relevant comparator software solutions cited in the literature, and expert evaluations. The computational tools and atlases described in this paper will be distributed to the research community as part of the FreeSurfer image analysis package.

Non-negative data-driven mapping of structural connections with application to the neonatal brain

Mapping connections in the neonatal brain can provide insight into the crucial early stages of neurodevelopment that shape brain organisation and lay the foundations for cognition and behaviour. Diffusion MRI and tractography provide unique opportunities for such explorations, through estimation of white matter bundles and brain connectivity. Atlas-based tractography protocols, i.e. a priori defined sets of masks and logical operations in a template space, have been commonly used in the adult brain to drive such explorations. However, rapid growth and maturation of the brain during early development make it challenging to ensure correspondence and validity of such atlas-based tractography approaches in the developing brain. An alternative can be provided by data-driven methods, which do not depend on predefined regions of interest. Here, we develop a novel data-driven framework to extract white matter bundles and their associated grey matter networks from neonatal tractography data, based on non-negative matrix factorisation that is inherently suited to the non-negative nature of structural connectivity data. We also develop a non-negative dual regression framework to map group-level components to individual subjects. Using in-silico simulations, we evaluate the accuracy of our approach in extracting connectivity components and compare with an alternative data-driven method, independent component analysis. We apply non-negative matrix factorisation to whole-brain connectivity obtained from publicly available datasets from the Developing Human Connectome Project, yielding grey matter components and their corresponding white matter bundles. We assess the validity and interpretability of these components against traditional tractography results and grey matter networks obtained from resting-state fMRI in the same subjects. We subsequently use them to generate a parcellation of the neonatal cortex using data from 323 new-born babies and we assess the robustness and reproducibility of this connectivity-driven parcellation.

Rich-club in the brain's macrostructure: Insights from graph theoretical analysis

The brain is a complex network. Growing evidence supports the critical roles of a set of brain regions within the brain network, known as the brain’s cores or hubs. These regions require high energy cost but possess highly efficient neural information transfer in the brain’s network and are termed the rich-club. The rich-club of the brain network is essential as it directly regulates functional integration across multiple segregated regions and helps to optimize cognitive processes. Here, we review the recent advances in rich-club organization to address the fundamental roles of the rich-club in the brain and discuss how these core brain regions affect brain development and disorders. We describe the concepts of the rich-club behind network construction in the brain using graph theoretical analysis. We also highlight novel insights based on animal studies related to the rich-club and illustrate how human studies using neuroimaging techniques for brain development and psychiatric/neurological disorders may be relevant to the rich-club phenomenon in the brain network.

Automatic, Age Consistent Reconstruction of the Corpus Callosum Guided by Coherency From In Utero Diffusion-Weighted MRI

Reconstruction of white matter connectivity in the fetal brain from in utero diffusion-weighted magnetic resonance imaging (MRI) faces many challenges, including subject motion, small anatomical scale, and limited image resolution and signal. These issues are compounded by the need to track significant changes in structural connectivity throughout development. We present an automated method for improved reliability and completeness of tract extraction across a wide range of gestational ages, based on the geometry of coherent patterns in streamline tractography, and apply it to the reconstruction of the corpus callosum. This method, focused specifically at addressing the challenges of fetal brain imaging, avoids depending on a tractography atlas, and handles variations in size, shape, and tissue properties of developing brains, both between subjects and across ages. Although tractography from in utero MRI generally suffers from a significant number of misleading and missing pathways, we demonstrate the feasibility of extracting the coherent bundle of the corpus callosum while avoiding inappropriate diversions into other tracts.

Multiscale Structure-Function Gradients in the Neonatal Connectome

The adult functional connectome is well characterized by a macroscale spatial gradient of connectivity traversing from unimodal toward higher-order transmodal cortices that recapitulates known principles of hierarchical organization and myelination patterns. Despite an emerging literature assessing connectome properties in neonates, the presence of connectome gradients and particularly their correspondence to microstructure remains largely unknown. We derived connectome gradients using unsupervised techniques applied to functional connectivity data from 40 term-born neonates. A series of cortex-wide analysis examined associations to magnetic resonance imaging-derived morphological parameters (cortical thickness, sulcal depth, curvature), measures of tissue microstructure (intracortical T1w/T2w intensity, superficial white matter diffusion parameters), and subcortico-cortical functional connectivity. Our findings indicate that the primary neonatal connectome gradient runs between sensorimotor and visual anchors and captures specific associations to cortical and superficial white matter microstructure as well as thalamo-cortical connectivity. A second gradient indicated an anterior-to-posterior asymmetry in macroscale connectivity alongside an immature differentiation between unimodal and transmodal areas, indicating a connectome-level circuitry en route to an adult-like organization. Our findings reveal an important coordination of structural and functional interactions in the neonatal connectome across spatial scales. Observed associations were replicable across individual neonates, suggesting consistency and generalizability.

Evaluation of Six Phase Encoding Based Susceptibility Distortion Correction Methods for Diffusion MRI

Purpose: Susceptibility distortions impact diffusion MRI data analysis and is typically corrected during preprocessing. Correction strategies involve three classes of methods: registration to a structural image, the use of a fieldmap, or the use of images acquired with opposing phase encoding directions. It has been demonstrated that phase encoding based methods outperform the other two classes, but unfortunately, the choice of which phase encoding based method to use is still an open question due to the absence of any systematic comparisons.

Methods: In this paper we quantitatively evaluated six popular phase encoding based methods for correcting susceptibility distortions in diffusion MRI data. We employed a framework that allows for the simulation of realistic diffusion MRI data with susceptibility distortions. We evaluated the ability for methods to correct distortions by comparing the corrected data with the ground truth. Four diffusion tensor metrics (FA, MD, eigenvalues and eigenvectors) were calculated from the corrected data and compared with the ground truth. We also validated two popular indirect metrics using both simulated data and real data. The two indirect metrics are the difference between the corrected LR and AP data, and the FA standard deviation over the corrected LR, RL, AP, and PA data.

Results: We found that DR-BUDDI and TOPUP offered the most accurate and robust correction compared to the other four methods using both direct and indirect evaluation metrics. EPIC and HySCO performed well in correcting b₀ images but produced poor corrections for diffusion weighted volumes, and also they produced large errors for the four diffusion tensor metrics. We also demonstrate that the indirect metric (the difference between corrected LR and AP data) gives a different ordering of correction quality than the direct metric.

Conclusion: We suggest researchers to use DR-BUDDI or TOPUP for susceptibility distortion correction. The two indirect metrics (the difference between corrected LR and AP data, and the FA standard deviation) should be interpreted together as a measure of distortion correction quality. The performance ranking of the various tools inferred from direct and indirect metrics differs slightly. However, across all tools, the results of direct and indirect metrics are highly correlated indicating that the analysis of indirect metrics may provide a good proxy of the performance of a correction tool if assessment using direct metrics is not feasible.

Comparison of quality control methods for automated diffusion tensor imaging analysis pipelines

The processing of brain diffusion tensor imaging (DTI) data for large cohort studies requires fully automatic pipelines to perform quality control (QC) and artifact/outlier removal procedures on the raw DTI data prior to calculation of diffusion parameters. In this study, three automatic DTI processing pipelines, each complying with the general ENIGMA framework, were designed by uniquely combining multiple image processing software tools. Different QC procedures based on the RESTORE algorithm, the DTIPrep protocol, and a combination of both methods were compared using simulated ground truth and artifact containing DTI datasets modeling eddy current induced distortions, various levels of motion artifacts, and thermal noise. Variability was also examined in 20 DTI datasets acquired in subjects with vascular cognitive impairment (VCI) from the multi-site Ontario Neurodegenerative Disease Research Initiative (ONDRI). The mean fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were calculated in global brain grey matter (GM) and white matter (WM) regions. For the simulated DTI datasets, the measure used to evaluate the performance of the pipelines was the normalized difference between the mean DTI metrics measured in GM and WM regions and the corresponding ground truth DTI value. The performance of the proposed pipelines was very similar, particularly in FA measurements. However, the pipeline based on the RESTORE algorithm was the most accurate when analyzing the artifact containing DTI datasets. The pipeline that combined the DTIPrep protocol and the RESTORE algorithm produced the lowest standard deviation in FA measurements in normal appearing WM across subjects. We concluded that this pipeline was the most robust and is preferred for automated analysis of multisite brain DTI data.

Music enhances structural maturation of emotional processing neural pathways in very preterm infants

Prematurity disrupts brain maturation by exposing the developing brain to different noxious stimuli present in the neonatal intensive care unit (NICU) and depriving it from meaningful sensory inputs during a critical period of brain development, leading to later neurodevelopmental impairments. Musicotherapy in the NICU environment has been proposed to promote sensory stimulation, relevant for activity-dependent brain plasticity, but its impact on brain structural maturation is unknown. Neuroimaging studies have demonstrated that music listening triggers neural substrates implied in socio-emotional processing and, thus, it might influence networks formed early in development and known to be affected by prematurity. Using multi-modal MRI, we aimed to evaluate the impact of a specially composed music intervention during NICU stay on preterm infant's brain structure maturation. 30 preterm newborns (out of which 15 were exposed to music during NICU stay and 15 without music intervention) and 15 full-term newborns underwent an MRI examination at term-equivalent age, comprising diffusion tensor imaging (DTI), used to evaluate white matter maturation using both region-of-interest and seed-based tractography approaches, as well as a T2-weighted image, used to perform amygdala volumetric analysis. Overall, WM microstructural maturity measured through DTI metrics was reduced in preterm infants receiving the standard-of-care in comparison to full-term newborns, whereas preterm infants exposed to the music intervention demonstrated significantly improved white matter maturation in acoustic radiations, external capsule/claustrum/extreme capsule and uncinate fasciculus, as well as larger amygdala volumes, in comparison to preterm infants with standard-of-care. These results suggest a structural maturational effect of the proposed music intervention on premature infants' auditory and emotional processing neural pathways during a key period of brain development.

Quantifying nerve decussation abnormalities in the optic chiasm

•

Diffusion MRI is capable of detecting structural abnormalities of the optic chiasm.

•

Quantification of crossing strength in optic chiasm is of promise for albinism diagnostics.

•

Optic chiasm is a powerful test model for neuroimaging methods resolving crossing fibers.

Objective

The human optic chiasm comprises partially crossing optic nerve fibers. Here we used diffusion MRI (dMRI) for the in-vivo identification of the abnormally high proportion of crossing fibers found in the optic chiasm of people with albinism.

Methods

In 9 individuals with albinism and 8 controls high-resolution 3T dMRI data was acquired and analyzed with a set of methods for signal modeling [Diffusion Tensor (DT) and Constrained Spherical Deconvolution (CSD)], tractography, and streamline filtering (LiFE, COMMIT, and SIFT2). The number of crossing and non-crossing streamlines and their weights after filtering entered ROC-analyses to compare the discriminative power of the methods based on the area under the curve (AUC). The dMRI results were cross-validated with fMRI estimates of misrouting in a subset of 6 albinotic individuals.

Results

We detected significant group differences in chiasmal crossing for both unfiltered DT (p = 0.014) and CSD tractograms (p = 0.0009) also reflected by AUC measures (for DT and CSD: 0.61 and 0.75, respectively), underlining the discriminative power of the approach. Estimates of crossing strengths obtained with dMRI and fMRI were significantly correlated for CSD (R² = 0.83, p = 0.012). The results show that streamline filtering methods in combination with probabilistic tracking, both optimized for the data at hand, can improve the detection of crossing in the human optic chiasm.

Conclusions

Especially CSD-based tractography provides an efficient approach to detect structural abnormalities in the optic chiasm. The most realistic results were obtained with filtering methods with parameters optimized for the data at hand.

Significance

Our findings demonstrate a novel anatomy-driven approach for the individualized diagnostics of optic chiasm abnormalities.

Towards an optimised processing pipeline for diffusion magnetic resonance imaging data: Effects of artefact corrections on diffusion metrics and their age associations in UK Biobank

Increasing interest in the structural and functional organisation of the human brain encourages the acquisition of big data sets comprising multiple neuroimaging modalities, often accompanied by additional information obtained from health records, cognitive tests, biomarkers and genotypes. Diffusion weighted magnetic resonance imaging data enables a range of promising imaging phenotypes probing structural connections as well as macroanatomical and microstructural properties of the brain. The reliability and biological sensitivity and specificity of diffusion data depend on processing pipeline. A state-of-the-art framework for data processing facilitates cross-study harmonisation and reduces pipeline-related variability. Using diffusion magnetic resonance imaging (MRI) data from 218 individuals in the UK Biobank, we evaluate the effects of different processing steps that have been suggested to reduce imaging artefacts and improve reliability of diffusion metrics. In lack of a ground truth, we compared diffusion metric sensitivity to age between pipelines. By comparing distributions and age sensitivity of the resulting diffusion metrics based on different approaches (diffusion tensor imaging, diffusion kurtosis imaging and white matter tract integrity), we evaluate a general pipeline comprising seven postprocessing blocks: noise correction; Gibbs ringing correction; evaluation of field distortions; susceptibility, eddy-current and motion-induced distortion corrections; bias field correction; spatial smoothing and final diffusion metric estimations. Based on this evaluation, we suggest an optimised processing pipeline for diffusion weighted MRI data.

Diffusion MRI: Assessment of the Impact of Acquisition and Preprocessing Methods Using the BrainVISA-Diffuse Toolbox

Diffusion MR images are prone to severe geometric distortions induced by head movement, eddy-current and inhomogeneity of magnetic susceptibility. Various correction methods have been proposed that depend on the choice of the acquisition settings and potentially provide highly different data quality. However, the impact of this choice has not been evaluated in terms of the ratio between scan time and preprocessed data quality. This study aims at investigating the impact of six well-known preprocessing methods, each associated to specific acquisition settings, on the outcome of diffusion analyses. For this purpose, we developed a comprehensive toolbox called Diffuse which automatically guides the user to the best preprocessing pipeline according to the input data. Using MR images of 20 subjects from the HCP dataset, we compared the six pre-processing pipelines regarding the following criteria: the ability to recover brain’s true geometry, the tensor model estimation and derived indices in the white matter, and finally the spatial dispersion of six well known connectivity pathways. As expected the pipeline associated to the longer acquisition fully repeated with reversed phase-encoding (RPE) yielded the higher data quality and was used as a reference to evaluate the other pipelines. In this way, we highlighted several significant aspects of other pre-processing pipelines. Our results first established that eddy-current correction improves the tensor-fitting performance with a localized impact especially in the corpus callosum. Concerning susceptibility distortions, we showed that the use of a field map is not sufficient and involves additional smoothing, yielding to an artificial decrease of tensor-fitting error. Of most importance, our findings demonstrate that, for an equivalent scan time, the acquisition of a b0 volume with RPE ensures a better brain’s geometry reconstruction and local improvement of tensor quality, without any smoothing of the image. This was found to be the best scan time/data quality compromise. To conclude, this study highlights and attempts to quantify the strong dependence of diffusion metrics on acquisition settings and preprocessing methods.

Asymmetry Spectrum Imaging for Baby Diffusion Tractography

Fiber tractography in baby diffusion MRI is challenging due to the low and spatially-varying diffusion anisotropy, causing most tractography algorithms to yield streamlines that fall short of reaching the cortex. In this paper, we introduce a method called asymmetry spectrum imaging (ASI) to improve the estimation of white matter pathways in the baby brain by (i) incorporating an asymmetric fiber orientation model to resolve subvoxel fiber configurations such as fanning and bending, and (ii) explicitly modeling the range (or spectrum) of typical diffusion length scales in the developing brain. We validated ASI using in-vivo baby diffusion MRI data from the Baby Connectome Project (BCP), demonstrating that ASI can characterize complex subvoxel fiber configurations and accurately estimate the fiber orientation distribution function in spite of changes in diffusion patterns. This, in turn, results in significantly better diffusion tractography in the baby brain.

Optimising neonatal fMRI data analysis: Design and validation of an extended dHCP preprocessing pipeline to characterise noxious-evoked brain activity in infants

The infant brain is unlike the adult brain, with considerable differences in morphological, neurodynamic, and haemodynamic features. As the majority of current MRI analysis tools were designed for use in adults, a primary objective of the Developing Human Connectome Project (dHCP) is to develop optimised methodological pipelines for the analysis of neonatal structural, resting state, and diffusion MRI data. Here, in an independent neonatal dataset we have extended and optimised the dHCP fMRI preprocessing pipeline for the analysis of stimulus-response fMRI data. We describe and validate this extended dHCP fMRI preprocessing pipeline to analyse changes in brain activity evoked following an acute noxious stimulus applied to the infant's foot. We compare the results obtained from this extended dHCP pipeline to results obtained from a typical FSL FEAT-based analysis pipeline, evaluating the pipelines' outputs using a wide range of tests. We demonstrate that a substantial increase in spatial specificity and sensitivity to signal can be attained with a bespoke neonatal preprocessing pipeline through optimised motion and distortion correction, ICA-based denoising, and haemodynamic modelling. The improved sensitivity and specificity, made possible with this extended dHCP pipeline, will be paramount in making further progress in our understanding of the development of sensory processing in the infant brain.