Concurrent mapping of brain ontogeny and phylogeny within a common space: Standardized tractography and applications.

Developmental and evolutionary effects on brain organization are complex, yet linked, as evidenced by the correspondence in cortical area expansion across these vastly different time scales. However, it is still not possible to study concurrently the ontogeny and phylogeny of cortical areal connections, which is arguably more relevant to brain function than allometric measurements. Here, we propose a novel framework that allows the integration of structural connectivity maps from humans (adults and neonates) and nonhuman primates (macaques) onto a common space. We use white matter bundles to anchor the common space and use the uniqueness of cortical connection patterns to these bundles to probe area specialization. This enabled us to quantitatively study divergences and similarities in connectivity over evolutionary and developmental scales, to reveal brain maturation trajectories, including the effect of premature birth, and to translate cortical atlases between diverse brains. Our findings open new avenues for an integrative approach to imaging neuroanatomy.

Coupling cognitive and brainstem dysfunction in multiple sclerosis-related chronic neuropathic limb pain.

Chronic pain in multiple sclerosis is common and difficult to treat. Its mechanisms remain incompletely understood. Dysfunction of the descending pain modulatory system is known to contribute to human chronic pain conditions. However, it is not clear how alterations in executive function influence this network, despite healthy volunteer studies linking function of the descending pain modulatory system, to cognition. In adults with multiple sclerosis-associated chronic neuropathic limb pain, compared to those without pain, we hypothesized altered functional connectivity of the descending pain modulatory system, coupled to executive dysfunction. Specifically we hypothesized reduced mental flexibility, because of potential importance in stimulus reappraisal. To investigate these hypotheses, we conducted a case-control cross-sectional study of 47 adults with relapsing remitting multiple sclerosis (31 with chronic neuropathic limb pain, 16 without pain), employing clinical, neuropsychological, structural, and functional MRI measures. We measured brain lesions and atrophy affecting descending pain modulatory system structures. Both cognitive and affective dysfunctions were confirmed in the chronic neuropathic limb pain group, including reduced mental flexibility (Delis Kaplan Executive Function System card sorting tests P < 0.001). Functional connectivity of rostral anterior cingulate and ventrolateral periaqueductal gray, key structures of the descending pain modulatory system, was significantly lower in the group experiencing chronic neuropathic pain. There was no significant between-group difference in whole-brain grey matter or lesion volumes, nor lesion volume affecting white matter tracts between rostral anterior cingulate and periaqueductal gray. Brainstem-specific lesion volume was higher in the chronic neuropathic limb pain group (P = 0.0017). Differential functional connectivity remained after correction for brainstem-specific lesion volume. Gabapentinoid medications were more frequently used in the chronic pain group. We describe executive dysfunction in people with multiple sclerosis affected by chronic neuropathic pain, along with functional and structural MRI evidence compatible with dysfunction of the descending pain modulatory system. These findings extend understanding of close inter-relationships between cognition, function of the descending pain modulatory system, and chronic pain, both in multiple sclerosis and more generally in human chronic pain conditions. These findings could support application of pharmacological and cognitive interventions in chronic neuropathic pain associated with multiple sclerosis.

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.

The association between inadequate sleep and accelerated brain ageing.

Numerous studies indicate large heterogeneity in brain ageing, which can be attributed to modifiable lifestyle factors, including sleep. Inadequate sleep has been previously linked to gray (GM) and white (WM) matter changes. However, the reported findings are highly inconsistent. By contrast to previous research independently characterizing patterns of either GM or WM changes, we used here linked independent component analysis (FLICA) to examine covariation in GM, and WM in a group of older adults (n = 50). Next, we employed a novel technique to estimate the brain age delta (difference between chronological and brain age assessed using neuroimaging data) and study its associations with sleep quality and sleep fragmentation, hypothesizing that inadequate sleep accelerates brain ageing. FLICA revealed a number of multimodal components, associated with age, sleep quality, and sleep fragmentation. Subsequently, we show significant associations between brain age delta and inadequate sleep, suggesting 2 years deviation above the chronological age. Our findings indicate sensitivity of multimodal approaches and brain age delta in detecting link between inadequate sleep and accelerated brain ageing.

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, B1 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.

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.

Multimodal Imaging Brain Markers in Early Adolescence Are Linked with a Physically Active Lifestyle.

The World Health Organization promotes physical exercise and a healthy lifestyle as means to improve youth development. However, relationships between physical lifestyle and human brain development are not fully understood. Here, we asked whether a human brain-physical latent mode of covariation underpins the relationship between physical activity, fitness, and physical health measures with multimodal neuroimaging markers. In 50 12-year old school pupils (26 females), we acquired multimodal whole-brain MRI, characterizing brain structure, microstructure, function, myelin content, and blood perfusion. We also acquired physical variables measuring objective fitness levels, 7 d physical activity, body mass index, heart rate, and blood pressure. Using canonical correlation analysis, we unravel a latent mode of brain-physical covariation, independent of demographics, school, or socioeconomic status. We show that MRI metrics with greater involvement in this mode also showed spatially extended patterns across the brain. Specifically, global patterns of greater gray matter perfusion, volume, cortical surface area, greater white matter extra-neurite density, and resting state networks activity covaried positively with measures reflecting a physically active phenotype (high fit, low sedentary individuals). Showing that a physically active lifestyle is linked with systems-level brain MRI metrics, these results suggest widespread associations relating to several biological processes. These results support the notion of close brain-body relationships and underline the importance of investigating modifiable lifestyle factors not only for physical health but also for brain health early in adolescence.SIGNIFICANCE STATEMENT An active lifestyle is key for healthy development. In this work, we answer the following question: How do brain neuroimaging markers relate with young adolescents' level of physical activity, fitness, and physical health? Combining advanced whole-brain multimodal MRI metrics with computational approaches, we show a robust relationship between physically active lifestyles and spatially extended, multimodal brain imaging-derived phenotypes. Suggesting a wider effect on brain neuroimaging metrics than previously thought, this work underlies the importance of studying physical lifestyle, as well as other brain-body relationships in an effort to foster brain health at this crucial stage in development.

Modelling white matter in gyral blades as a continuous vector field.

Many brain imaging studies aim to measure structural connectivity with diffusion tractography. However, biases in tractography data, particularly near the boundary between white matter and cortical grey matter can limit the accuracy of such studies. When seeding from the white matter, streamlines tend to travel parallel to the convoluted cortical surface, largely avoiding sulcal fundi and terminating preferentially on gyral crowns. When seeding from the cortical grey matter, streamlines generally run near the cortical surface until reaching deep white matter. These so-called "gyral biases" limit the accuracy and effective resolution of cortical structural connectivity profiles estimated by tractography algorithms, and they do not reflect the expected distributions of axonal densities seen in invasive tracer studies or stains of myelinated fibres. We propose an algorithm that concurrently models fibre density and orientation using a divergence-free vector field within gyral blades to encourage an anatomically-justified streamline density distribution along the cortical white/grey-matter boundary while maintaining alignment with the diffusion MRI estimated fibre orientations. Using in vivo data from the Human Connectome Project, we show that this algorithm reduces tractography biases. We compare the structural connectomes to functional connectomes from resting-state fMRI, showing that our model improves cross-modal agreement. Finally, we find that after parcellation the changes in the structural connectome are very minor with slightly improved interhemispheric connections (i.e, more homotopic connectivity) and slightly worse intrahemispheric connections when compared to tracers.

Confound modelling in UK Biobank brain imaging.

Dealing with confounds is an essential step in large cohort studies to address problems such as unexplained variance and spurious correlations. UK Biobank is a powerful resource for studying associations between imaging and non-imaging measures such as lifestyle factors and health outcomes, in part because of the large subject numbers. However, the resulting high statistical power also raises the sensitivity to confound effects, which therefore have to be carefully considered. In this work we describe a set of possible confounds (including non-linear effects and interactions that researchers may wish to consider for their studies using such data). We include descriptions of how we can estimate the confounds, and study the extent to which each of these confounds affects the data, and the spurious correlations that may arise if they are not controlled. Finally, we discuss several issues that future studies should consider when dealing with confounds.

The developing Human Connectome Project (dHCP) automated resting-state functional processing framework for newborn infants.

The developing Human Connectome Project (dHCP) aims to create a detailed 4-dimensional connectome of early life spanning 20-45 weeks post-menstrual age. This is being achieved through the acquisition of multi-modal MRI data from over 1000 in- and ex-utero subjects combined with the development of optimised pre-processing pipelines. In this paper we present an automated and robust pipeline to minimally pre-process highly confounded neonatal resting-state fMRI data, robustly, with low failure rates and high quality-assurance. The pipeline has been designed to specifically address the challenges that neonatal data presents including low and variable contrast and high levels of head motion. We provide a detailed description and evaluation of the pipeline which includes integrated slice-to-volume motion correction and dynamic susceptibility distortion correction, a robust multimodal registration approach, bespoke ICA-based denoising, and an automated QC framework. We assess these components on a large cohort of dHCP subjects and demonstrate that processing refinements integrated into the pipeline provide substantial reduction in movement related distortions, resulting in significant improvements in SNR, and detection of high quality RSNs from neonates.

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.

Improved fibre dispersion estimation using b-tensor encoding.

Measuring fibre dispersion in white matter with diffusion magnetic resonance imaging (MRI) is limited by an inherent degeneracy between fibre dispersion and microscopic diffusion anisotropy (i.e., the diffusion anisotropy expected for a single fibre orientation). This means that estimates of fibre dispersion rely on strong assumptions, such as constant microscopic anisotropy throughout the white matter or specific biophysical models. Here we present a simple approach for resolving this degeneracy using measurements that combine linear (conventional) and spherical tensor diffusion encoding. To test the accuracy of the fibre dispersion when our microstructural model is only an approximation of the true tissue structure, we simulate multi-compartment data and fit this with a single-compartment model. For such overly simplistic tissue assumptions, we show that the bias in fibre dispersion is greatly reduced (~5x) for single-shell linear and spherical tensor encoding data compared with single-shell or multi-shell conventional data. In in-vivo data we find a consistent estimate of fibre dispersion as we reduce the b-value from 3 to 1.5 ms/µm2, increase the repetition time, increase the echo time, or increase the diffusion time. We conclude that the addition of spherical tensor encoded data to conventional linear tensor encoding data greatly reduces the sensitivity of the estimated fibre dispersion to the model assumptions of the tissue microstructure.

Towards HCP-Style macaque connectomes: 24-Channel 3T multi-array coil, MRI sequences and preprocessing.

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

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.

Joint modelling of diffusion MRI and microscopy.

The combination of diffusion MRI (dMRI) with microscopy provides unique opportunities to study microstructural features of tissue, particularly when acquired in the same sample. Microscopy is frequently used to validate dMRI microstructure models, addressing the indirect nature of dMRI signals. Typically, these modalities are analysed separately, and microscopy is taken as a gold standard against which dMRI-derived parameters are validated. Here we propose an alternative approach in which we combine dMRI and microscopy data obtained from the same tissue sample to drive a single, joint model. This simultaneous analysis allows us to take advantage of the breadth of information provided by complementary data acquired from different modalities. By applying this framework to a spherical-deconvolution analysis, we are able to overcome a known degeneracy between fibre dispersion and radial diffusion. Spherical-deconvolution based approaches typically estimate a global fibre response function to determine the fibre orientation distribution in each voxel. However, the assumption of a 'brain-wide' fibre response function may be challenged if the diffusion characteristics of white matter vary across the brain. Using a generative joint dMRI-histology model, we demonstrate that the fibre response function is dependent on local anatomy, and that current spherical-deconvolution based models may be overestimating dispersion and underestimating the number of distinct fibre populations per voxel.

Characterizing Microstructural Tissue Properties in Multiple Sclerosis with Diffusion MRI at 7 T and 3 T: The Impact of the Experimental Design.

The recent introduction of advanced magnetic resonance (MR) imaging techniques to characterize focal and global degeneration in multiple sclerosis (MS), like the Composite Hindered and Restricted Model of Diffusion, or CHARMED, diffusional kurtosis imaging (DKI) and Neurite Orientation Dispersion and Density Imaging (NODDI) made available new tools to image axonal pathology non-invasively in vivo. These methods already showed greater sensitivity and specificity compared to conventional diffusion tensor-based metrics (e.g., fractional anisotropy), overcoming some of its limitations. While previous studies uncovered global and focal axonal degeneration in MS patients compared to healthy controls, here our aim is to investigate and compare different diffusion MRI acquisition protocols in their ability to highlight microstructural differences between MS and control tissue over several much used models. For comparison, we contrasted the ability of fractional anisotropy measurements to uncover differences between lesion, normal-appearing white matter (WM), gray matter and healthy tissue under the same imaging protocols. We show that: (1) focal and diffuse differences in several microstructural parameters are observed under clinical settings; (2) advanced models (CHARMED, DKI and NODDI) have increased specificity and sensitivity to neurodegeneration when compared to fractional anisotropy measurements; and (3) both high (3 T) and ultra-high fields (7 T) are viable options for imaging tissue change in MS lesions and normal appearing WM, while higher b-values are less beneficial under the tested short-time (10 min acquisition) conditions.

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.

Automated quality control for within and between studies diffusion MRI data using a non-parametric framework for movement and distortion correction.

Diffusion MRI data can be affected by hardware and subject-related artefacts that can adversely affect downstream analyses. Therefore, automated quality control (QC) is of great importance, especially in large population studies where visual QC is not practical. In this work, we introduce an automated diffusion MRI QC framework for single subject and group studies. The QC is based on a comprehensive, non-parametric approach for movement and distortion correction: FSL EDDY, which allows us to extract a rich set of QC metrics that are both sensitive and specific to different types of artefacts. Two different tools are presented: QUAD (QUality Assessment for DMRI), for single subject QC and SQUAD (Study-wise QUality Assessment for DMRI), which is designed to enable group QC and facilitate cross-studies harmonisation efforts.

Extending the Human Connectome Project across ages: Imaging protocols for the Lifespan Development and Aging projects.

The Human Connectome Projects in Development (HCP-D) and Aging (HCP-A) are two large-scale brain imaging studies that will extend the recently completed HCP Young-Adult (HCP-YA) project to nearly the full lifespan, collecting structural, resting-state fMRI, task-fMRI, diffusion, and perfusion MRI in participants from 5 to 100+ years of age. HCP-D is enrolling 1300+ healthy children, adolescents, and young adults (ages 5-21), and HCP-A is enrolling 1200+ healthy adults (ages 36-100+), with each study collecting longitudinal data in a subset of individuals at particular age ranges. The imaging protocols of the HCP-D and HCP-A studies are very similar, differing primarily in the selection of different task-fMRI paradigms. We strove to harmonize the imaging protocol to the greatest extent feasible with the completed HCP-YA (1200+ participants, aged 22-35), but some imaging-related changes were motivated or necessitated by hardware changes, the need to reduce the total amount of scanning per participant, and/or the additional challenges of working with young and elderly populations. Here, we provide an overview of the common HCP-D/A imaging protocol including data and rationales for protocol decisions and changes relative to HCP-YA. The result will be a large, rich, multi-modal, and freely available set of consistently acquired data for use by the scientific community to investigate and define normative developmental and aging related changes in the healthy human brain.