Pushing spatial and temporal resolution for functional and diffusion MRI in the Human Connectome Project

The Duke Mouse Brain Atlas: MRI and light sheet microscopy stereotaxic atlas of the mouse brain

Atlases of the brain are critical resources that make it possible to share data in a common reference frame. Unexpectedly, there is no three-dimensional (3D) stereotaxic atlas of the mouse brain that provides whole brain coverage at macro to single-cell levels. Diffusion tensor images from five perfusion-fixed (in skull) specimens were acquired at 15 micrometers, the highest resolution ever reported. Diffusion tensor imaging yields multiple 3D volumes, each of which highlights unique cytoarchitecture. The averages were mapped into micro-computed tomography of the mouse skull to create external landmarks (bregma and lambda). Light sheet images of the same brains were coregistered, providing cell maps in the same stereotaxic space. The Allen Reference Atlas was registered to the volume to correct the geometric distortion in that atlas and bring it into the stereotaxic space. The resulting multiscalar (13 terabytes) atlas provides a common spatial framework to anneal data across molecular, structural, and functional studies of mice.

A closer look: pediatric neuroimaging at 7T

Although 7 Tesla (7T) field strength MR imaging offers higher signal-to-noise ratio and spatial resolution and improves certain types of tissue contrast, the incorporation of these systems into clinical pediatric neuroradiology has been relatively limited. Following a discussion of available hardware, current regulations, and pediatric specific safety considerations, this article briefly reviews the underlying principles behind the improved image quality attainable with certain techniques at 7T. Subsequently, specific high-performance sequences and techniques are highlighted including MP2RAGE, T2-weighted, and T2*-weighted sequences as well as MR angiography, all with sample images and comparison with standard field strengths. Finally, current clinical neuroradiological applications of 7T are explored with particular focus on focal epilepsy, multiple sclerosis, vascular diseases, and cerebral microbleeds. Ongoing and future innovations in hardware design and sequence development promise continued advancement in 7T neuroimaging and further applications to pediatric neuroradiology.

Towards whole brain mapping of the haemodynamic response function

Functional magnetic resonance imaging time-series are conventionally processed by linear modelling the evoked response as the convolution of the experimental conditions with a stereotyped haemodynamic response function (HRF). However, the neural signal in response to a stimulus can vary according to task, brain region, and subject-specific conditions. Moreover, HRF shape has been suggested to carry physiological information. The BOLD signal across a range of sensorial and cognitive tasks was fitted using a sine series expansion, and modelled signals were deconvolved, thus giving rise to a task-specific deconvolved HRF (dHRF), which was characterized in terms of amplitude, latency, time-to-peak and full-width at half maximum for each task. We found that the BOLD response shape changes not only across activated regions and tasks, but also across subjects despite the age homogeneity of the cohort. Largest variabilities were observed in mean amplitude and latency across tasks and regions, while time-to-peak and full width at half maximum were relatively more consistent. Additionally, the dHRF was found to deviate from canonicity in several brain regions. Our results suggest that the choice of a standard, uniform HRF may be not optimal for all fMRI analyses and may lead to model misspecifications and statistical bias.

The MR neuroimaging protocol for the Accelerating Medicines Partnership® Schizophrenia Program

Neuroimaging with MRI has been a frequent component of studies of individuals at clinical high risk (CHR) for developing psychosis, with goals of understanding potential brain regions and systems impacted in the CHR state and identifying prognostic or predictive biomarkers that can enhance our ability to forecast clinical outcomes. To date, most studies involving MRI in CHR are likely not sufficiently powered to generate robust and generalizable neuroimaging results. Here, we describe the prospective, advanced, and modern neuroimaging protocol that was implemented in a complex multi-site, multi-vendor environment, as part of the large-scale Accelerating Medicines Partnership® Schizophrenia Program (AMP® SCZ), including the rationale for various choices. This protocol includes T1- and T2-weighted structural scans, resting-state fMRI, and diffusion-weighted imaging collected at two time points, approximately 2 months apart. We also present preliminary variance component analyses of several measures, such as signal- and contrast-to-noise ratio (SNR/CNR) and spatial smoothness, to provide quantitative data on the relative percentages of participant, site, and platform (i.e., scanner model) variance. Site-related variance is generally small (typically <10%). For the SNR/CNR measures from the structural and fMRI scans, participant variance is the largest component (as desired; 40-76%). However, for SNR/CNR in the diffusion scans, there is substantial platform-related variance (>55%) due to differences in the diffusion imaging hardware capabilities of the different scanners. Also, spatial smoothness generally has a large platform-related variance due to inherent, difficult to control, differences between vendors in their acquisitions and reconstructions. These results illustrate some of the factors that will need to be considered in analyses of the AMP SCZ neuroimaging data, which will be the largest CHR cohort to date.Watch Dr. Harms discuss this article at https://vimeo.com/1059777228?share=copy#t=0 .

Multivariate Pattern Analysis of fMRI Reveals Striato-Cortical Network Changes in Myoclonus-Dystonia

BACKGROUND

Currently, the pathophysiology of myoclonus-dystonia (M-D) remains insufficiently understood. This study addresses this gap by adding innovative multivariate pattern analysis (MVPA) to traditional univariate analysis of functional magnetic resonance imaging (fMRI) data.

METHODS

Data from 18 M-D patients and 18 age-matched healthy volunteers who performed a finger tapping fMRI task were analyzed. Whole-brain univariate and searchlight (MVPA) analysis with varying hemodynamic response function (HRF) delays were employed to examine brain responses associated with the visually guided motor task.

RESULTS

Distinguishing response patterns between M-D patients and healthy volunteers revealed significant response reductions in the putamen, insula, and visual cortex. Compared to univariate analysis, searchlight analysis was more sensitive for brain activity patterns associated with finger tapping in both M-D patients and healthy volunteers. At short HRF delays, increased (pre)motor cortical responses were evident in M-D patients, whereas such responses emerged at a later HRF delay in healthy volunteers.

CONCLUSION

The task-related effects observed in M-D patients support the involvement of the basal ganglia-thalamo-cortical network. Notably, cerebellar involvement was not strongly implicated in our study. We postulate that inherent deficits in the putamen trigger either premature or downstream compensatory (motor) cortical effects. The potential involvement of the visual cortex in the M-D pathophysiology is new, but its role has been suggested by a previous study investigating visual sensory processing in SGCE gene-positive M-D patients. Our findings, including the innovative searchlight method, pave the way for further studies investigating the complex interplay between brain regions and networks and their role in M-D pathogenesis.

Intrinsic brain mapping of cognitive abilities: A multiple-dataset study on intelligence and its components

This study investigates how functional brain network features contribute to general intelligence and its cognitive components by analyzing three independent cohorts of healthy participants. Cognitive scores were derived from 1) the Wechsler Adult Intelligence Scale (WAIS-IV), 2) the Raven Standard Progressive Matrices (RPM), and 3) the NIH and Penn cognitive batteries from the Human Connectome Project. Factor analysis on the NIH and Penn cognitive batteries yielded latent variables that closely resembled the content of the WAIS-IV indices and RPM. We employed graph theory and a multi-resolution network analysis by varying the modularity parameter (γ) to investigate hierarchical brain-behavior relationships across different scales of brain organization. Brain-behavior associations were quantified using multi-level robust regression analyses to accommodate variability and confounds at the subject-level, node-level, and resolution-level. Our findings reveal consistent brain-behavior relationships across the datasets. Nodal efficiency in fronto-parietal sensorimotor regions consistently played a pivotal role in fluid reasoning, whereas efficiency in visual networks was linked to executive functions and memory. A broad, low-resolution 'task-positive' network emerged as predictive of full-scale IQ scores, indicating a hierarchical brain-behavior coding. Conversely, increased cross-network connections involving default mode and subcortical-limbic networks were associated with reductions in both general and specific cognitive performance. These outcomes highlight the relevance of network efficiency and integration, as well as of the hierarchical organization in supporting specific aspects of intelligence, while recognizing the inherent complexity of these relationships. Our multi-resolution network approach offers new insights into the interplay between multilayer network properties and the structure of cognitive abilities, advancing the understanding of the neural substrates of the intelligence construct.

Multimodal precision MRI of the individual human brain at ultra-high fields

Multimodal neuroimaging, in particular magnetic resonance imaging (MRI), allows for non-invasive examination of human brain structure and function across multiple scales. Precision neuroimaging builds upon this foundation, enabling the mapping of brain structure, function, and connectivity patterns with high fidelity in single individuals. Highfield MRI, operating at magnetic field strengths of 7 Tesla (T) or higher, increases signal-to-noise ratio and opens up possibilities for gains spatial resolution. Here, we share a multimodal Precision Neuroimaging and Connectomics (PNI) 7 T MRI dataset. Ten healthy individuals underwent a comprehensive MRI protocol, including T1 relaxometry, magnetization transfer imaging, T2*-weighted imaging, diffusion MRI, and multi-state functional MRI paradigms, aggregated across three imaging sessions. Alongside anonymized raw MRI data, we release cortex-wide connectomes from different modalities across multiple parcellation scales, and supply "gradients" that compactly characterize spatial patterning of cortical organization. Our precision MRI dataset will advance our understanding of structure-function relationships in the individual human brain and is publicly available via the Open Science Framework.

Apparent Diffusion Coefficient fMRI shines light on white matter resting-state connectivity compared to BOLD

Resting-state functional magnetic resonance imaging (fMRI) is used to derive functional connectivity (FC) between brain regions. Typically, blood oxygen level-dependent (BOLD) contrast is used. However, BOLD's reliance on neurovascular coupling poses challenges in reflecting brain activity accurately, leading to reduced sensitivity in white matter (WM). WM BOLD signals have long been considered physiological noise, although recent evidence shows that both stimulus-evoked and resting-state WM BOLD signals resemble those in gray matter (GM), albeit smaller in amplitude. We introduce apparent diffusion coefficient fMRI (ADC-fMRI) as a promising functional contrast for GM and WM FC, capturing activity-driven neuromorphological fluctuations. Our study compares BOLD-fMRI and ADC-fMRI FC in GM and WM, showing that ADC-fMRI mirrors BOLD-fMRI connectivity in GM, while capturing more robust FC in WM. ADC-fMRI displays higher average clustering and average node strength in WM, and higher inter-subject similarity, compared to BOLD. Taken together, this suggests that ADC-fMRI is a reliable tool for exploring FC that incorporates gray and white matter nodes in a novel way.

Single-shot echo planar time-resolved imaging for multi-echo functional MRI and distortion-free diffusion imaging

PURPOSE

To develop a single-shot SNR-efficient distortion-free multi-echo imaging technique for dynamic imaging applications.

METHODS

Echo planar time-resolved imaging (EPTI) was first introduced as a multi-shot technique for distortion-free multi-echo imaging. This work aims to develop single-shot EPTI (ss-EPTI) to achieve improved robustness to motion/physiological noise, increased temporal resolution, and higher SNR efficiency. A new spatiotemporal encoding that enables reduced phase-encoding blips and minimized echo spacing under the single-shot regime was developed, which improves sampling efficiency and enhances spatiotemporal correlation in the k-TE space for improved reconstruction. A continuous readout with minimized deadtime was employed to optimize SNR efficiency. Moreover, k-TE partial Fourier and simultaneous multi-slice acquisition were integrated for further acceleration.

RESULTS

ss-EPTI provided distortion-free imaging with densely sampled multi-echo images at standard resolutions (e.g., ˜1.25 to 3 mm) in a single-shot. Improved SNR efficiency was observed in ss-EPTI due to improved motion/physiological-noise robustness and efficient continuous readout. Its ability to eliminate dynamic distortions-geometric changes across dynamics due to field changes induced by physiological variations or eddy currents-further improved the data's temporal stability. For multi-echo fMRI, ss-EPTI's multi-echo images recovered signal dropout in short-T 2 * $$ {\mathrm{T}}_2^{\ast } $$ regions and provided TE-dependent functional information to distinguish non-BOLD noise for further tSNR improvement. For diffusion MRI, it achieved shortened TEs for improved SNR and provided images free from both B0-induced and diffusion-encoding-dependent eddy-current-induced distortions with multi-TE diffusion metrics.

CONCLUSION

ss-EPTI provides SNR-efficient distortion-free multi-echo imaging with comparable temporal resolutions to ss-EPI, offering a new acquisition tool for dynamic imaging.

Exploring the transmission of cognitive task information through optimal brain pathways

Understanding the large-scale information processing that underlies complex human cognition is the central goal of cognitive neuroscience. While emerging activity flow models demonstrate that cognitive task information is transferred by interregional functional or structural connectivity, graph-theory-based models typically assume that neural communication occurs via the shortest path of brain networks. However, whether the shortest path is the optimal route for empirical cognitive information transmission remains unclear. Based on a large-scale activity flow mapping framework, we found that the performance of activity flow prediction with the shortest path was significantly lower than that with the direct path. The shortest path routing was superior to other network communication strategies, including search information, path ensembles, and navigation. Intriguingly, the shortest path outperformed the direct path in activity flow prediction when the physical distance constraint and asymmetric routing contribution were simultaneously considered. This study not only challenges the shortest path assumption through empirical network models but also suggests that cognitive task information routing is constrained by the spatial and functional embedding of the brain network.

Dynamics of infant white matter maturation from birth to 6 months

The first months after a baby's birth encompass the most rapid period of postnatal change in the human lifespan, but longitudinal trajectories of white matter maturation in this period remain uncharted. Using densely sampled diffusion tensor images collected longitudinally at a mean rate of 1 scan per 1.55 days, we measured non-linear growth and growth rate trajectories of major white matter tracts from birth to 6 months. Growth rates at birth were 6 to 11 times faster than at 6 months, with tracts less mature at birth developing fastest. When matched on chronological age, shorter gestation infants had less mature white matter at birth but faster growth rates than their longer gestation peers; however, growth trajectories were highly similar when corrected for gestational age. This is the first study to estimate white matter trajectories using dense sampling in the first 6 post-natal months, which can inform the study of neurodevelopmental disorders beginning in infancy.

Stable Cortical Body Maps Before and After Arm Amputation

The adult brain's capacity for cortical reorganization remains debated. Using longitudinal neuroimaging in three adults, followed up to five years before and after arm amputation, we compared cortical activity elicited by movement of the hand (pre-amputation) versus phantom hand (post-amputation) and lips (pre/post-amputation). We observed stable representations of both hand and lips. By directly quantifying activity changes across amputation, we overturn decades of animal and human research, demonstrating amputation does not trigger large-scale cortical reorganization.

Tale of Two n-Backs: Diverging Associations of Dorsolateral Prefrontal Cortex Activation With n-Back Task Performance

In studying the neural correlates of working memory (WM) ability via functional magnetic resonance imaging (fMRI) in health and disease, it is relatively uncommon for investigators to report associations between brain activation and measures of task performance. Additionally, how the choice of WM task impacts observed activation-performance relationships is poorly understood. We sought to illustrate the impact of WM task on brain-behavior correlations using two large, publicly available datasets. We conducted between-participants analyses of task-based fMRI data from two publicly available datasets: The Human Connectome Project (HCP; n = 866) and the Queensland Twin Imaging (QTIM) Study (n = 459). Participants performed two distinct variations of the n-back WM task with different stimuli, timings, and response paradigms. Associations between brain activation ([2-back - 0-back] contrast) and task performance (2-back % correct) were investigated separately in each dataset, as well as across datasets, within the dorsolateral prefrontal cortex (dlPFC), medial prefrontal cortex, and whole cortex. Global patterns of activation to task were similar in both datasets. However, opposite associations between activation and task performance were observed in bilateral pre-supplementary motor area and left middle frontal gyrus. Within the dlPFC, HCP participants exhibited a significantly greater activation-performance relationship in bilateral middle frontal gyrus relative to QTIM Study participants. The observation of diverging activation-performance relationships between two large datasets performing variations of the n-back task serves as a critical reminder for investigators to exercise caution when selecting WM tasks and interpreting neural activation in response to a WM task.

RF coil design strategies for improving SNR at the ultrahigh magnetic field of 10.5T

PURPOSE

Toward pushing the boundaries of ultrahigh fields for human brain imaging, we wish to evaluate experimentally achievable SNR relative to ultimate intrinsic SNR (uiSNR) at 10.5T, develop design strategies toward approaching the latter, quantify magnetic field-dependent SNR gains, and demonstrate the feasibility of whole-brain, high-resolution human brain imaging at this uniquely high field strength.

METHODS

A dual row 16-channel self-decoupled transmit (Tx) and receive (Rx) array was developed for 10.5T using custom Tx/Rx switches. A 64-channel receive-only array was built to fit into the 16-channel Tx/Rx array. Electromagnetic modeling and experiments were used to define safe operational power limits. Experimental SNR was evaluated relative to uiSNR at 10.5T and 7T.

RESULTS

The 64-channel Rx array alone captured approximately 50% of the central uiSNR at 10.5T, while an identical array developed for 7T captured about 76% of uiSNR at 7T. The 16-channel Tx/80-channel Rx configuration brought the fraction of uiSNR captured at 10.5T to levels comparable to the 64-channel Rx array at 7T. SNR data displayed an approximateB 0 2 $$ {\mathrm{B}}_0^2 $$ dependence over a large central region when evaluated in the context of uiSNR. Whole-brain, high-resolutionT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted and T1-weighted anatomical and gradient-recalled-echo BOLD-EPI functional MRI images were obtained at 10.5T for the first time with such an advanced array.

CONCLUSION

We demonstrated the ability to approach the uiSNR at 10.5T over the human brain, achieving large SNR gains over 7T, currently the most commonly used ultrahigh-field platform. Whole-brain, high-resolution anatomical and EPI-based functional MRI data were obtained at 10.5T, illustrating the promise of greater than 10T fields in studying the human brain.

qPRF: A system to accelerate population receptive field modeling

BOLD response can be fitted using the population receptive field (PRF) model to reveal how visual input is represented on the cortex (Dumoulin and Wandell, 2008). Fitting the PRF model costs considerable time, often requiring days to analyze BOLD signals for a small cohort of subjects. We introduce the qPRF ("quick PRF"), a system for accelerated PRF modeling that reduced the computation time by a factor >1,000 without losing goodness-of-fit when compared to another widely available PRF modeling package (Kay et al., 2013) on a benchmark of data from the Human Connectome Project (HCP; Van Essen et al. (2013). The system achieves this level of acceleration by pre-computing a tree-like data structure, which it rapidly searches during the fitting step for an optimal parameter combination. We tested the method on a constrained four-parameter version of the PRF model (Strategy 1 herein) and an unconstrained five-parameter PRF model, which the qPRF fitted at comparable speed (Strategy 2). We show how an additional search step can guarantee optimality of qPRF solutions with little additional time cost (Strategy 3). To assess the quality of qPRF solutions, we compared our Strategy 1 solutions to those provided by Benson et al. (2018) who performed a similar four-parameter fit. Both hemispheres of the 181 subjects in the HCP dataset (a total of 10,753,572 vertices, each with a unique BOLD time series of 1800 frames) were analyzed by qPRF in 12.82 h on an ordinary CPU. The absolute difference in R2 achieved by the qPRF compared to Benson et al. (2018) was negligible, with a median of 0.025% (R2 units being between 0% and 100%). In general, the qPRF yielded a slightly better fitting solution, achieving a greater R2 on 70.2% of vertices. We also assess the qPRF method's model-recovery ability using a simulated dataset. The qPRF may facilitate the development and use of more elaborate models based on the PRF framework and may pave the way for novel clinical applications.

Investigating time-independent and time-dependent diffusion phenomena using steady-state diffusion MRI

Diffusion MRI is a leading method to non-invasively characterise brain tissue microstructure across multiple domains and scales. Diffusion-weighted steady-state free precession (DW-SSFP) is an established imaging sequence for post-mortem MRI, addressing the challenging imaging environment of fixed tissue with short T2 and low diffusivities. However, a current limitation of DW-SSFP is signal interpretation: it is not clear what diffusion 'regime' the sequence probes and therefore its potential to characterise tissue microstructure. Building on Extended Phase Graphs (EPG), I establish two alternative representations of the DW-SSFP signal in terms of (1) conventional b-values (time-independent diffusion) and (2) encoding power-spectra (time-dependent diffusion). The proposed representations provide insights into how different parameter regimes and gradient waveforms impact the diffusion sensitivity of DW-SSFP. I subsequently introduce an approach to incorporate existing biophysical models into DW-SSFP without the requirement of extensive derivations, with time dependence estimated via a Gaussian phase approximation representation of the DW-SSFP signal. Investigations incorporating free-diffusion and tissue-relevant microscopic restrictions (cylinder of varying radius) give excellent agreement to complementary analytical models and Monte Carlo simulations. Experimentally, the time-independent representation is used to derive Tensor and proof-of-principle NODDI estimates in a whole human post-mortem brain. A final SNR-efficiency investigation demonstrates the theoretical potential of DW-SSFP for ultra-high field microstructural imaging.

A double dissociation between semantic and spatial cognition in visual to default network pathways

Processing pathways between sensory and default mode network (DMN) regions support recognition, navigation, and memory but their organisation is not well understood. We show that functional subdivisions of visual cortex and DMN sit at opposing ends of parallel streams of information processing that support visually mediated semantic and spatial cognition, providing convergent evidence from univariate and multivariate task responses, intrinsic functional and structural connectivity. Participants learned virtual environments consisting of buildings populated with objects, drawn from either a single semantic category or multiple categories. Later, they made semantic and spatial context decisions about these objects and buildings during functional magnetic resonance imaging. A lateral ventral occipital to fronto-temporal DMN pathway was primarily engaged by semantic judgements, while a medial visual to medial temporal DMN pathway supported spatial context judgements. These pathways had distinctive locations in functional connectivity space: the semantic pathway was both further from unimodal systems and more balanced between visual and auditory-motor regions compared with the spatial pathway. When semantic and spatial context information could be integrated (in buildings containing objects from a single category), regions at the intersection of these pathways responded, suggesting that parallel processing streams interact at multiple levels of the cortical hierarchy to produce coherent memory-guided cognition.

Multi-study fMRI outlooks on subcortical BOLD responses in the stop-signal paradigm

This study investigates the functional network underlying response inhibition in the human brain, particularly the role of the basal ganglia in successful action cancellation. Functional magnetic resonance imaging (fMRI) approaches have frequently used the stop-signal task to examine this network. We merge five such datasets, using a novel aggregatory method allowing the unification of raw fMRI data across sites. This meta-analysis, along with other recent aggregatory fMRI studies, does not find evidence for the innervation of the hyperdirect or indirect cortico-basal-ganglia pathways in successful response inhibition. What we do find, is large subcortical activity profiles for failed stop trials. We discuss possible explanations for the mismatch of findings between the fMRI results presented here and results from other research modalities that have implicated nodes of the basal ganglia in successful inhibition. We also highlight the substantial effect smoothing can have on the conclusions drawn from task-specific general linear models. First and foremost, this study presents a proof of concept for meta-analytical methods that enable the merging of extensive, unprocessed, or unreduced datasets. It demonstrates the significant potential that open-access data sharing can offer to the research community. With an increasing number of datasets being shared publicly, researchers will have the ability to conduct meta-analyses on more than just summary data.

Category-biased patches encircle core domain-general regions in the human lateral prefrontal cortex

The fine-grained functional organization of the human lateral prefrontal cortex (PFC) remains poorly understood. Previous fMRI studies delineated focal domain-general, or multiple-demand (MD), PFC areas that co-activate during diverse cognitively demanding tasks. While there is some evidence for category-selective (face and scene) patches, in human and non-human primate PFC, these have not been systematically assessed. Recent precision fMRI studies have also revealed sensory-biased PFC patches adjacent to MD regions. To investigate if this topographic arrangement extends to other domains, we analysed two independent fMRI datasets (n=449 and n=37) utilizing the high-resolution multimodal MRI approaches of the Human Connectome Project (HCP). Both datasets included cognitive control tasks and stimuli spanning different categories: faces, places, tools and body parts. Contrasting each stimulus category against the remaining ones revealed focal interdigitated patches of activity located adjacent to core MD regions. The results were robust, replicating across different executive tasks, experimental designs (block and event-related) and at the single subject level. Our results paint a refined view of the fine-grained functional organization of the PFC, revealing a recurring motif of interdigitated domain-specific and domain-general circuits. This organization offers new constraints for models of cognitive control, cortical specialization and development.

A generative model of the connectome with dynamic axon growth

Connectome generative models, otherwise known as generative network models, provide insight into the wiring principles underpinning brain network organization. While these models can approximate numerous statistical properties of empirical networks, they typically fail to explicitly characterize an important contributor to brain organization-axonal growth. Emulating the chemoaffinity-guided axonal growth, we provide a novel generative model in which axons dynamically steer the direction of propagation based on distance-dependent chemoattractive forces acting on their growth cones. This simple dynamic growth mechanism, despite being solely geometry-dependent, is shown to generate axonal fiber bundles with brain-like geometry and features of complex network architecture consistent with the human brain, including lognormally distributed connectivity weights, scale-free nodal degrees, small-worldness, and modularity. We demonstrate that our model parameters can be fitted to individual connectomes, enabling connectome dimensionality reduction and comparison of parameters between groups. Our work offers an opportunity to bridge studies of axon guidance and connectome development, providing new avenues for understanding neural development from a computational perspective.

Closed-loop fMRI at the mesoscopic scale of columns and layers: Can we do it and why would we want to?

Technological advances in fMRI including ultra-high magnetic fields (≥ 7 T) and acquisition methods that increase spatial specificity have paved the way for studies of the human cortex at the scale of layers and columns. This mesoscopic scale promises an improved mechanistic understanding of human cortical function so far only accessible to invasive animal neurophysiology. In recent years, an increasing number of studies have applied such methods to better understand the cortical function in perception and cognition. This future perspective article asks whether closed-loop fMRI studies could equally benefit from these methods to achieve layer and columnar specificity. We outline potential applications and discuss the conceptual and concrete challenges, including data acquisition and volitional control of mesoscopic brain activity. We anticipate an important role of fMRI with mesoscopic resolution for closed-loop fMRI and neurofeedback, yielding new insights into brain function and potentially clinical applications.This article is part of the theme issue 'Neurofeedback: new territories and neurocognitive mechanisms of endogenous neuromodulation'.

Accelerating Heritability, Genetic Correlation, and Genome-Wide Association Imaging Genetic Analyses in Complex Pedigrees

National and international biobanking efforts led to the collection of large and inclusive imaging genetics datasets that enable examination of the contribution of genetic and environmental factors to human brains in illness and health. High-resolution neuroimaging (~104-6 voxels) and genetic (106-8 single nucleotide polymorphic [SNP] variants) data are available in statistically powerful (N = 103-5) epidemiological and disorder-focused samples. Performing imaging genetics analyses at full resolution afforded in these datasets is a formidable computational task even under the assumption of unrelatedness among the subjects. The computational complexity rises as ~N2-3 (where N is the sample size), when accounting for relatedness among subjects. We describe fast, non-iterative simplifications to accelerate classical variance component (VC) methods including heritability, genetic correlation, and genome-wide association in dense and complex empirical pedigrees. These approaches linearize (from N2-3 to N~1) computational effort while maintaining fidelity (r ~ 0.95) with the VC results and take advantage of parallel computing provided by central and graphics processing units (CPU and GPU). We show that the new approaches lead to a 104- to 106-fold reduction in computational complexity-making voxel-wise heritability, genetic correlation, and genome-wide association studies (GWAS) analysis practical for large and complex samples such as those provided by the Amish and Human Connectome Projects (N = 406 and 1052 subjects, respectively) and UK Biobank (N = 31,681). These developments are shared in open-source, SOLAR-Eclipse software.

Track-Weighted Dynamic Functional Connectivity Profiles and Topographic Organization of the Human Pulvinar

The human pulvinar is considered a prototypical associative thalamic nucleus as it represents a key node in several cortico-subcortical networks. Through this extensive connectivity to widespread brain areas, it has been suggested that the pulvinar may play a central role in modulating cortical oscillatory dynamics of complex cognitive and executive functions. Additionally, derangements of pulvinar activity are involved in different neuropsychiatric conditions including Lewy-body disease, Alzheimer's disease, and schizophrenia. Anatomical investigations in nonhuman primates have demonstrated a topographical organization of cortico-pulvinar connectivity along its dorsoventral and rostrocaudal axes; this specific organization shows only partial overlap with the traditional subdivision into subnuclei (anterior, lateral, medial, and inferior) and is thought to coordinate information processing within specific brain networks. However, despite its relevance in mediating higher-order cognitive functions, such a structural and functional organization of the pulvinar in the human brain remains poorly understood. Track-weighted dynamic functional connectivity (tw-dFC) is a recently developed technique that combines structural and dynamic functional connectivity, allowing the identification of white matter pathways underlying the fluctuations observed in functional connectivity between brain regions over time. Herein, we applied a data-driven parcellation approach to reveal topographically organized connectivity clusters within the human pulvinar complex, in two large cohorts of healthy human subjects. Unsupervised clustering of tw-dFC time series within the pulvinar complex revealed dorsomedial, dorsolateral, ventral anterior, and ventral posterior connectivity clusters. Each of these clusters shows functional coupling to specific, widespread cortico-subcortical white matter brain networks. Altogether, our findings represent a relevant step towards a better understanding of pulvinar anatomy and function, and a detailed characterization of his role in healthy and pathological conditions.

Accelerated High-Resolution Deep Learning Reconstruction Turbo Spin Echo MRI of the Knee at 7 T

OBJECTIVES

The aim of this study was to compare the image quality of 7 T turbo spin echo (TSE) knee images acquired with varying factors of parallel-imaging acceleration reconstructed with deep learning (DL)-based and conventional algorithms.

MATERIALS AND METHODS

This was a prospective single-center study. Twenty-three healthy volunteers underwent 7 T knee magnetic resonance imaging. Two-, 3-, and 4-fold accelerated high-resolution fat-signal-suppressing proton density (PD-fs) and T1-weighted coronal 2D TSE acquisitions with an encoded voxel volume of 0.31 × 0.31 × 1.5 mm 3 were acquired. Each set of raw data was reconstructed with a DL-based and a conventional Generalized Autocalibrating Partially Parallel Acquisition (GRAPPA) algorithm. Three readers rated image contrast, sharpness, artifacts, noise, and overall quality. Friedman analysis of variance and the Wilcoxon signed rank test were used for comparison of image quality criteria.

RESULTS

The mean age of the participants was 32.0 ± 8.1 years (15 male, 8 female). Acquisition times at 4-fold acceleration were 4 minutes 15 seconds (PD-fs, Supplemental Video is available at http://links.lww.com/RLI/A938 ) and 3 minutes 9 seconds (T1, Supplemental Video available at http://links.lww.com/RLI/A939 ). At 4-fold acceleration, image contrast, sharpness, noise, and overall quality of images reconstructed with the DL-based algorithm were significantly better rated than the corresponding GRAPPA reconstructions ( P < 0.001). Four-fold accelerated DL-reconstructed images scored significantly better than 2- to 3-fold GRAPPA-reconstructed images with regards to image contrast, sharpness, noise, and overall quality ( P ≤ 0.031). Image contrast of PD-fs images at 2-fold acceleration ( P = 0.087), image noise of T1-weighted images at 2-fold acceleration ( P = 0.180), and image artifacts for both sequences at 2- and 3-fold acceleration ( P ≥ 0.102) of GRAPPA reconstructions were not rated differently than those of 4-fold accelerated DL-reconstructed images. Furthermore, no significant difference was observed for all image quality measures among 2-fold, 3-fold, and 4-fold accelerated DL reconstructions ( P ≥ 0.082).

CONCLUSIONS

This study explored the technical potential of DL-based image reconstruction in accelerated 2D TSE acquisitions of the knee at 7 T. DL reconstruction significantly improved a variety of image quality measures of high-resolution TSE images acquired with a 4-fold parallel-imaging acceleration compared with a conventional reconstruction algorithm.

A review of diffusion MRI in mood disorders: mechanisms and predictors of treatment response

By measuring the molecular diffusion of water molecules in brain tissue, diffusion MRI (dMRI) provides unique insight into the microstructure and structural connections of the brain in living subjects. Since its inception, the application of dMRI in clinical research has expanded our understanding of the possible biological bases of psychiatric disorders and successful responses to different therapeutic interventions. Here, we review the past decade of diffusion imaging-based investigations with a specific focus on studies examining the mechanisms and predictors of therapeutic response in people with mood disorders. We present a brief overview of the general application of dMRI and key methodological developments in the field that afford increasingly detailed information concerning the macro- and micro-structural properties and connectivity patterns of white matter (WM) pathways and their perturbation over time in patients followed prospectively while undergoing treatment. This is followed by a more in-depth summary of particular studies using dMRI approaches to examine mechanisms and predictors of clinical outcomes in patients with unipolar or bipolar depression receiving pharmacological, neurostimulation, or behavioral treatments. Limitations associated with dMRI research in general and with treatment studies in mood disorders specifically are discussed, as are directions for future research. Despite limitations and the associated discrepancies in findings across individual studies, evolving research strongly indicates that the field is on the precipice of identifying and validating dMRI biomarkers that could lead to more successful personalized treatment approaches and could serve as targets for evaluating the neural effects of novel treatments.

Characterization and Mitigation of a Simultaneous Multi-Slice fMRI Artifact: Multiband Artifact Regression in Simultaneous Slices

Simultaneous multi-slice (multiband) acceleration in fMRI has become widespread, but may be affected by novel forms of signal artifact. Here, we demonstrate a previously unreported artifact manifesting as a shared signal between simultaneously acquired slices in all resting-state and task-based multiband fMRI datasets we investigated, including publicly available consortium data from the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) Study. We propose Multiband Artifact Regression in Simultaneous Slices (MARSS), a regression-based detection and correction technique that successfully mitigates this shared signal in unprocessed data. We demonstrate that the signal isolated by MARSS correction is likely nonneural, appearing stronger in neurovasculature than gray matter. Additionally, we evaluate MARSS both against and in tandem with sICA+FIX denoising, which is implemented in HCP resting-state data, to show that MARSS mitigates residual artifact signal that is not modeled by sICA+FIX. MARSS correction leads to study-wide increases in signal-to-noise ratio, decreases in cortical coefficient of variation, and mitigation of systematic artefactual spatial patterns in participant-level task betas. Finally, MARSS correction has substantive effects on second-level t-statistics in analyses of task-evoked activation. We recommend that investigators apply MARSS to multiband fMRI datasets with moderate or higher acceleration factors, in combination with established denoising methods.

Radiofrequency Enhancer to Recover Signal Dropouts in 7 Tesla Diffusion MRI

Diffusion magnetic resonance imaging (dMRI) allows for a non-invasive visualization and quantitative assessment of white matter architecture in the brain by characterizing restrictions on the random motion of water molecules. Ultra-high field MRI scanners, such as those operating at 7 Tesla (7T) or higher, can boost the signal-to-noise ratio (SNR) to improve dMRI compared with what is attainable at conventional field strengths such as 3T or 1.5T. However, wavelength effects at 7T cause reduced transmit magnetic field efficiency in the human brain, mainly in the posterior fossa, manifesting as signal dropouts in this region. Recently, we reported a simple approach of using a wireless radiofrequency (RF) surface array to improve transmit efficiency and signal sensitivity at 7T. In this study, we demonstrate the feasibility and effectiveness of the RF enhancer in improving in vivo dMRI at 7T. The electromagnetic simulation results demonstrated a 2.1-fold increase in transmit efficiency with the use of the RF enhancer. The experimental results similarly showed a 1.9-fold improvement in transmit efficiency and a 1.4-fold increase in normalized SNR. These improvements effectively mitigated signal dropouts in regions with inherently lower SNR, such as the cerebellum, resulting in a better depiction of principal fiber orientations and an enhanced visualization of extended tracts.

Functional Localization of the Human Auditory and Visual Thalamus Using a Thalamic Localizer Functional Magnetic Resonance Imaging Task

Functional magnetic resonance imaging (fMRI) of the auditory and visual sensory systems of the human brain is an active area of investigation in the study of human health and disease. The medial geniculate nucleus (MGN) and lateral geniculate nucleus (LGN) are key thalamic nuclei involved in the processing and relay of auditory and visual information, respectively, and are the subject of blood-oxygen-level-dependent (BOLD) fMRI studies of neural activation and functional connectivity in human participants. However, localization of BOLD fMRI signal originating from neural activity in MGN and LGN remains a technical challenge, due in part to the poor definition of boundaries of these thalamic nuclei in standard T1-weighted and T2-weighted magnetic resonance imaging sequences. Here, we report the development and evaluation of an auditory and visual sensory thalamic localizer (TL) fMRI task that produces participant-specific functionally-defined regions of interest (fROIs) of both MGN and LGN, using 3 Tesla multiband fMRI and a clustered-sparse temporal acquisition sequence, in less than 16 minutes of scan time. We demonstrate the use of MGN and LGN fROIs obtained from the TL fMRI task in standard resting-state functional connectivity (RSFC) fMRI analyses in the same participants. In RSFC analyses, we validated the specificity of MGN and LGN fROIs for signals obtained from primary auditory and visual cortex, respectively, and benchmark their performance against alternative atlas- and segmentation-based localization methods. The TL fMRI task and analysis code (written in Presentation and MATLAB, respectively) have been made freely available to the wider research community.

A Virtual In Vivo Dissection and Analysis of Socioaffective Symptoms Related to Cerebellum-Midbrain Reward Circuitry in Humans

Emerging research in nonhuman animals implicates cerebellar projections to the ventral tegmental area (VTA) in appetitive behaviors, but these circuits have not been characterized in humans. Here, we mapped cerebello-VTA white matter connectivity in a cohort of men and women using probabilistic tractography on diffusion imaging data from the Human Connectome Project. We uncovered the topographical organization of these connections by separately tracking from parcels of cerebellar lobule VI, crus I/II, vermis, paravermis, and cerebrocerebellum. Results revealed that connections between the cerebellum and VTA predominantly originate in the right cerebellar hemisphere, interposed nucleus, and paravermal cortex and terminate mostly ipsilaterally. Paravermal crus I sends the most connections to the VTA compared with other lobules. We discovered a mediolateral gradient of connectivity, such that the medial cerebellum has the highest connectivity with the VTA. Individual differences in microstructure were associated with measures of negative affect and social functioning. By splitting the tracts into quarters, we found that the socioaffective effects were driven by the third quarter of the tract, corresponding to the point at which the fibers leave the deep nuclei. Taken together, we produced detailed maps of cerebello-VTA structural connectivity for the first time in humans and established their relevance for trait differences in socioaffective regulation.

The use of 7T MRI in multiple sclerosis: review and consensus statement from the North American Imaging in Multiple Sclerosis Cooperative

The use of ultra-high-field 7-Tesla (7T) MRI in multiple sclerosis (MS) research has grown significantly over the past two decades. With recent regulatory approvals of 7T scanners for clinical use in 2017 and 2020, the use of this technology for routine care is poised to continue to increase in the coming years. In this context, the North American Imaging in MS Cooperative (NAIMS) convened a workshop in February 2023 to review the previous and current use of 7T technology for MS research and potential future research and clinical applications. In this workshop, experts were tasked with reviewing the current literature and proposing a series of consensus statements, which were reviewed and approved by the NAIMS. In this review and consensus paper, we provide background on the use of 7T MRI in MS research, highlighting this technology's promise for identification and quantification of aspects of MS pathology that are more difficult to visualize with lower-field MRI, such as grey matter lesions, paramagnetic rim lesions, leptomeningeal enhancement and the central vein sign. We also review the promise of 7T MRI to study metabolic and functional changes to the brain in MS. The NAIMS provides a series of consensus statements regarding what is currently known about the use of 7T MRI in MS, and additional statements intended to provide guidance as to what work is necessary going forward to accelerate 7T MRI research in MS and translate this technology for use in clinical practice and clinical trials. This includes guidance on technical development, proposals for a universal acquisition protocol and suggestions for research geared towards assessing the utility of 7T MRI to improve MS diagnostics, prognostics and therapeutic efficacy monitoring. The NAIMS expects that this article will provide a roadmap for future use of 7T MRI in MS.

Heterogenous brain activations across individuals localize to a common network

Task functional magnetic resonance imaging research has generally shielded away from studying individuals due to the low reproducibility. Here, we propose that heterogeneous brain activations across individuals localize to a common network. To test this hypothesis, we use working memory (WM) as our example. First, we showed that discrete-brain-based reproducibility of brain activation during WM across individuals was low. Then, we used activation network mapping (ANM) technique to identify each individual's brain network of WM and found that network-based reproducibility was rather high. Prediction analyses using machine learning algorithms indicated that individual WM networks identified via ANM can predict WM behavioral performance. This predictive ability even outperformed that of brain activations. Our study provides a new explanation on the low reproducibility of brain activations across individuals. The results suggest that ANM can be used to identify individual brain networks of cognitive processes, thus promising broad potential applications.

Comparing single- and multi-post labeling delays for the measurements of resting cerebral and hippocampal blood flow for cerebrovascular testing in midlife adults

Objectives

To assess the reliability and validity of measuring resting cerebral blood flow (CBF) and hippocampal CBF using a single-post-labeling delay (PLD) and a multi-PLD pseudo-continuous arterial spin labeling (pCASL) protocol for cerebrovascular reactivity (CVR) testing.

Methods

25 healthy, midlife adults (57 ± 4 years old) were imaged in a Siemens Prisma 3T magnetic resonance imaging (MRI) scanner. Resting CBF and hippocampal CBF were assessed using two pCASL protocols, our modified single-PLD protocol (pCASL-MOD) to accommodate the needs for CVR testing and the multi-PLD Human Connectome Project (HCP) Lifespan protocol to serve as the reference control (pCASL-HCP). During pCASL-MOD, CVR was calculated as the change in CBF from rest to hypercapnia (+9 mmHg increase in end-tidal partial pressure of carbon dioxide [PETCO2]) and then normalized for PETCO2. The reliability and validity in resting gray matter (GM) CBF, white matter (WM) CBF, and hippocampal CBF between pCASL-MOD and pCASL-HCP protocols were examined using correlation analyses, paired t-tests, and Bland Altman plots.

Results

The pCASL-MOD and pCASL-HCP protocols were significantly correlated for resting GM CBF [r = 0.72; F (1, 23) = 25.24, p < 0.0001], WM CBF [r = 0.57; F (1, 23) = 10.83, p = 0.003], and hippocampal CBF [r = 0.77; F (1, 23) = 32.65, p < 0.0001]. However, pCASL-MOD underestimated resting GM CBF (pCASL-MOD: 53.7 ± 11.1 v. pCASL-HCP: 69.1 ± 13.1 mL/100 g/min; p < 0.0001), WM CBF (pCASL-MOD: 32.4 ± 4.8 v. pCASL-HCP: 35.5 ± 6.9 mL/100 g/min; p = 0.01), and hippocampal CBF (pCASL-MOD: 50.5 ± 9.0 v. pCASL-HCP: 68.1 ± 12.5 mL/100 g/min; p < 0.0001). PETCO2 increased by 8.0 ± 0.7 mmHg to induce CVR (GM CBF: 4.8% ± 2.6%; WM CBF 2.9% ± 2.5%; and hippocampal CBF: 3.4% ± 3.8%).

Conclusion

Our single-PLD pCASL-MOD protocol reliably measured CBF and hippocampal CBF at rest given the significant correlation with the multi-PLD pCASL-HCP protocol. Despite the lower magnitude relative to pCASL-HCP, we recommend using our pCASL-MOD protocol for CVR testing in which an exact estimate of CBF is not required such as the assessment of relative change in CBF to hypercapnia.

Overview of ADNI MRI

The magnetic resonance imaging (MRI) Core has been operating since Alzheimer's Disease Neuroimaging Initiative's (ADNI) inception, providing 20 years of data including reliable, multi-platform standardized protocols, carefully curated image data, and quantitative measures provided by expert investigators. The overarching purposes of the MRI Core include: (1) optimizing and standardizing MRI acquisition methods, which have been adopted by many multicenter studies and trials worldwide and (2) providing curated images and numeric summary values from relevant MRI sequences/contrasts to the scientific community. Over time, ADNI MRI has become increasingly complex. To remain technically current, the ADNI MRI protocol has changed substantially over the past two decades. The ADNI 4 protocol contains nine different imaging types (e.g., three dimensional [3D] T1-weighted and fluid-attenuated inversion recovery [FLAIR]). Our view is that the ADNI MRI data are a greatly underutilized resource. The purpose of this paper is to educate the scientific community on ADNI MRI methods and content to promote greater awareness, accessibility, and use. HIGHLIGHTS: The MRI Core provides multi-platform standardized protocols, carefully curated image data, and quantitative analysis by expert groups. The ADNI MRI protocol has undergone major changes over the past two decades to remain technically current. As of April 25, 2024, the following numbers of image series are available: 17,141 3D T1w; 6877 FLAIR; 3140 T2/PD; 6623 GRE; 3237 dMRI; 2846 ASL; 2968 TF-fMRI; and 2861 HighResHippo (see Table 1 for abbreviations). As of April 25, 2024, the following numbers of quantitative analyses are available: FreeSurfer 10,997; BSI 6120; tensor based morphometry (TBM) and TBM-SYN 12,019; WMH 9944; dMRI 1913; ASL 925; TF-fMRI NFQ 2992; and medial temporal subregion volumes 2726 (see Table 4 for abbreviations). ADNI MRI is an underutilized resource that could be more useful to the research community.

The degenerate coding of psychometric profiles through functional connectivity archetypes

Introduction

Degeneracy in the brain-behavior code refers to the brain's ability to utilize different neural configurations to support similar functions, reflecting its adaptability and robustness. This study aims to explore degeneracy by investigating the non-linear associations between psychometric profiles and resting-state functional connectivity (RSFC).

Methods

The study analyzed RSFC data from 500 subjects to uncover the underlying neural configurations associated with various psychometric outcomes. Self-organized maps (SOM), a type of unsupervised machine learning algorithm, were employed to cluster the RSFC data. And identify distinct archetypal connectivity profiles characterized by unique within- and between-network connectivity patterns.

Results

The clustering analysis using SOM revealed several distinct archetypal connectivity profiles within the RSFC data. Each archetype exhibited unique connectivity patterns that correlated with various cognitive, physical, and socioemotional outcomes. Notably, the interaction between different SOM dimensions was significantly associated with specific psychometric profiles.

Discussion

This study underscores the complexity of brain-behavior interactions and the brain's capacity for degeneracy, where different neural configurations can lead to similar behavioral outcomes. These findings highlight the existence of multiple brain architectures capable of producing similar behavioral outcomes, illustrating the concept of neural degeneracy, and advance our understanding of neural degeneracy and its implications for cognitive and emotional health.

Stable individual differences from dynamic patterns of function: brain network flexibility predicts openness/intellect, intelligence, and psychoticism

A growing understanding of the nature of brain function has led to increased interest in interpreting the properties of large-scale brain networks. Methodological advances in network neuroscience provide means to decompose these networks into smaller functional communities and measure how they reconfigure over time as an index of their dynamic and flexible properties. Recent evidence has identified associations between flexibility and a variety of traits pertaining to complex cognition including creativity and working memory. The present study used measures of dynamic resting-state functional connectivity in data from the Human Connectome Project (n = 994) to test associations with Openness/Intellect, general intelligence, and psychoticism, three traits that involve flexible cognition. Using a machine-learning cross-validation approach, we identified reliable associations of intelligence with cohesive flexibility of parcels in large communities across the cortex, of psychoticism with disjoint flexibility, and of Openness/Intellect with overall flexibility among parcels in smaller communities. These findings are reasonably consistent with previous theories of the neural correlates of these traits and help to expand on previous associations of behavior with dynamic functional connectivity, in the context of broad personality dimensions.

In vivo diffusion MRI of the human heart using a 300 mT/m gradient system

PURPOSE

This work reports for the first time on the implementation and application of cardiac diffusion-weighted MRI on a Connectom MR scanner with a maximum gradient strength of 300 mT/m. It evaluates the benefits of the increased gradient performance for the investigation of the myocardial microstructure.

METHODS

Cardiac diffusion-weighted imaging (DWI) experiments were performed on 10 healthy volunteers using a spin-echo sequence with up to second- and third-order motion compensation (M 2 $$ {M}_2 $$ andM 3 $$ {M}_3 $$ ) andb = 100 ,   450 $$ b=100,450 $$ , and 1000s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ (twice theb max $$ {b}_{\mathrm{max}} $$ commonly used on clinical scanners). Mean diffusivity (MD), fractional anisotropy (FA), helix angle (HA), and secondary eigenvector angle (E2A) were calculated for b = [100, 450]s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ and b = [100, 1000]s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ for bothM 2 $$ {M}_2 $$ andM 3 $$ {M}_3 $$ .

RESULTS

The MD values withM 3 $$ {M}_3 $$ are slightly higher than withM 2 $$ {M}_2 $$ withΔ MD = 0 . 05 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 5 ) $$ \Delta \mathrm{MD}=0.05\pm 0.05\kern0.3em \left[\times 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-5\right) $$ forb max = 450 s / mm 2 $$ {b}_{\mathrm{max}}=450\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ andΔ MD = 0 . 03 ± 0 . 03 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 4 ) $$ \Delta \mathrm{MD}=0.03\pm 0.03\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-4\right) $$ forb max = 1000 s / mm 2 $$ {b}_{\mathrm{max}}=1000\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ . A reduction in MD is observed by increasing theb max $$ {b}_{\mathrm{max}} $$ from 450 to 1000s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ (Δ MD = 0 . 06 ± 0 . 04 [ × 1 0 - 3 mm 2 / s ] ( p = 1 . 6 e - 9 ) $$ \Delta \mathrm{MD}=0.06\pm 0.04\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1.6e-9\right) $$ forM 2 $$ {M}_2 $$ andΔ MD = 0 . 08 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 1 e - 9 ) $$ \Delta \mathrm{MD}=0.08\pm 0.05\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1e-9\right) $$ forM 3 $$ {M}_3 $$ ). The difference between FA, E2A, and HA was not significant in different schemes (p > 0 . 05 $$ p>0.05 $$ ).

CONCLUSION

This work demonstrates cardiac DWI in vivo with higher b-value and higher order of motion compensated diffusion gradient waveforms than is commonly used. Increasing the motion compensation order fromM 2 $$ {M}_2 $$ toM 3 $$ {M}_3 $$ and the maximum b-value from 450 to 1000  s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ affected the MD values but FA and the angular metrics (HA and E2A) remained unchanged. Our work paves the way for cardiac DWI on the next-generation MR scanners with high-performance gradient systems.

Sleep Quality Moderates the Associations between Cardiorespiratory Fitness and Hippocampal and Entorhinal Volume in Middle-Aged and Older Adults

INTRODUCTION/PURPOSE

As individuals age, the entorhinal cortex (ERC) and hippocampus-crucial structures for memory-tend to atrophy, with related cognitive decline. Simultaneously, lifestyle factors that can be modified, such as exercise and sleep, have been separately linked to slowing of brain atrophy and functional decline. However, the synergistic impact of fitness and sleep on susceptible brain structures in aging adults remains uncertain.

METHODS

We examined both independent and interactive associations of fitness and subjective sleep quality with regard to ERC thickness and hippocampal volume in 598 middle-aged and older adults from the Human Connectome Lifespan Aging Project. Cardiorespiratory fitness was assessed using the 2-min walk test, whereas subjective sleep quality was measured with the continuous Pittsburgh Sleep Quality Index global score. High-resolution structural magnetic resonance imaging was used to examine mean ERC thickness and bilateral hippocampal volume. Through multiple linear regression analyses, we investigated the moderating effects of subjective sleep quality on the association between fitness and brain structure, accounting for age, sex, education, body mass index, gait speed, and subjective physical activity.

RESULTS

We found that greater cardiorespiratory fitness, but not subjective sleep quality, was positively associated with bilateral hippocampal volume and ERC thickness. Notably, significant interaction effects suggest that poor subjective sleep quality was associated with a weaker association between fitness and both hippocampal volume and ERC thickness.

CONCLUSIONS

Findings suggest the potential importance of both cardiorespiratory fitness and subjective sleep quality in preserving critical, age-vulnerable brain structures. Interventions targeting brain health should consider potential combined effects of sleep and fitness on brain health.

Contributions of network structure, chemoarchitecture and diagnostic categories to transitions between cognitive topographies

The mechanisms linking the brain's network structure to cognitively relevant activation patterns remain largely unknown. Here, by leveraging principles of network control, we show how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic database. Specifically, we systematically integrated large-scale multimodal neuroimaging data from functional magnetic resonance imaging, diffusion tractography, cortical morphometry and positron emission tomography to simulate how anatomically guided transitions between cognitive states can be reshaped by neurotransmitter engagement or by changes in cortical thickness. Our model incorporates neurotransmitter-receptor density maps (18 receptors and transporters) and maps of cortical thickness pertaining to a wide range of mental health, neurodegenerative, psychiatric and neurodevelopmental diagnostic categories (17,000 patients and 22,000 controls). The results provide a comprehensive look-up table charting how brain network organization and chemoarchitecture interact to manifest different cognitive topographies, and establish a principled foundation for the systematic identification of ways to promote selective transitions between cognitive topographies.

qPRF: A system to accelerate population receptive field decoding

Patterns of BOLD response can be decoded using the population receptive field (PRF) model to reveal how visual input is represented on the cortex (Dumoulin and Wandell, 2008). The time cost of evaluating the PRF model is high, often requiring days to decode BOLD signals for a small cohort of subjects. We introduce the qPRF, an efficient method for decoding that reduced the computation time by a factor of 1436 when compared to another widely available PRF decoder (Kay, Winawer, Mezer and Wandell, 2013) on a benchmark of data from the Human Connectome Project (HCP; Van Essen, Smith, Barch, Behrens, Yacoub and Ugurbil, 2013). With a specially designed data structure and an efficient search algorithm, the qPRF optimizes the five PRF model parameters according to a least-squares criterion. To verify the accuracy of the qPRF solutions, we compared them to those provided by Benson, Jamison, Arcaro, Vu, Glasser, Coalson, Van Essen, Yacoub, Ugurbil, Winawer and Kay (2018). Both hemispheres of the 181 subjects in the HCP data set (a total of 10,753,572 vertices, each with a unique BOLD time series of 1800 frames) were decoded by qPRF in 15.2 hours on an ordinary CPU. The absolute difference inR 2 reported by Benson et al. and achieved by the qPRF was negligible, with a median of 0.39% (R 2 units being between 0% and 100%). In general, the qPRF yielded a slightly better fitting solution, achieving a greaterR 2 on 99.7% of vertices. The qPRF may facilitate the development and computation of more elaborate models based on the PRF framework, as well as the exploration of novel clinical applications.

The connectivity-based architecture of the human piriform cortex

The anatomy of the human piriform cortex (PC) is poorly understood. We used a bimodal connectivity-based-parcellation approach to investigate subregions of the PC and its connectional differentiation from the amygdala. One hundred (55 % female) genetically unrelated subjects from the Human Connectome Project were included. A region of interest (ROI) was delineated bilaterally covering PC and amygdala, and functional and structural connectivity of this ROI with the whole gray matter was computed. Spectral clustering was performed to obtain bilateral parcellations at granularities of k = 2-10 clusters and combined bimodal parcellations were computed. Validity of parcellations was assessed via their mean individual-to-group similarity per adjusted rand index (ARI). Individual-to-group similarity was higher than chance in both modalities and in all clustering solutions. The amygdala was clearly distinguished from PC in structural parcellations, and olfactory amygdala was connectionally more similar to amygdala than to PC. At higher granularities, an anterior and ventrotemporal and a posterior frontal cluster emerged within PC, as well as an additional temporal cluster at their boundary. Functional parcellations also showed a frontal piriform cluster, and similar temporal clusters were observed with less consistency. Results from bimodal parcellations were similar to the structural parcellations. Consistent results were obtained in a validation cohort. Distinction of the human PC from the amygdala, including its olfactory subregions, is possible based on its structural connectivity alone. The canonical fronto-temporal boundary within PC was reproduced in both modalities and with consistency. All obtained parcellations are freely available.

Main Paths of Brain Fibers in Diffusion Images Mixed with a Noise to Improve Performance of Tractography Algorithm-Evaluation in Phantom

Background

Some voxels may alter the tractography results due to unintentional alteration of noises and other unwanted factors.

Objective

This study aimed to investigate the effect of local phase features on tractography results providing data are mixed by a Gaussian or random distribution noise.

Material and Methods

In this simulation study, a mask was firstly designed based on the local phase features to decrease false-negative and -positive tractography results. The local phase features are calculated according to the local structures of images, which can be zero-dimensional, meaning just one point (equivalent to noise in tractography algorithm), a line (equivalent to a simple fiber), or an edge (equivalent to structures more complex than a simple fiber). A digital phantom evaluated the feasibility current model with the maximum complexities of configurations in fibers, including crossing fibers. In this paper, the diffusion images were mixed separately by a Gaussian or random distribution noise in 2 forms a zero-mean noise and a noise with a mean of data.

Results

The local mask eliminates the pixels of unfitted values with the main structures of images, due to noise or other interferer factors.

Conclusion

The local phase features of diffusion images are an innovative solution to determine principal diffusion directions.

Sex and mental health are related to subcortical brain microstructure

Some mental health problems such as depression and anxiety are more common in females, while others such as autism and attention deficit/hyperactivity (AD/H) are more common in males. However, the neurobiological origins of these sex differences are poorly understood. Animal studies have shown substantial sex differences in neuronal and glial cell structure, while human brain imaging studies have shown only small differences, which largely reflect overall body and brain size. Advanced diffusion MRI techniques can be used to examine intracellular, extracellular, and free water signal contributions and provide unique insights into microscopic cellular structure. However, the extent to which sex differences exist in these metrics of subcortical gray matter structures implicated in psychiatric disorders is not known. Here, we show large sex-related differences in microstructure in subcortical regions, including the hippocampus, thalamus, and nucleus accumbens in a large sample of young adults. Unlike conventional T1-weighted structural imaging, large sex differences remained after adjustment for age and brain volume. Further, diffusion metrics in the thalamus and amygdala were associated with depression, anxiety, AD/H, and antisocial personality problems. Diffusion MRI may provide mechanistic insights into the origin of sex differences in behavior and mental health over the life course and help to bridge the gap between findings from experimental, epidemiological, and clinical mental health research.

Exploring complex and integrated information during sleep

The Integrated Information Theory is a theoretical framework that aims to elucidate the nature of consciousness by postulating that it emerges from the integration of information within a system, and that the degree of consciousness depends on the extent of information integration within the system. When consciousness is lost, the core complex of consciousness proposed by the Integrated Information Theory disintegrates, and Φ measures, which reflect the level of integrated information, are expected to diminish. This study examined the predictions of the Integrated Information Theory using the global brain network acquired via functional magnetic resonance imaging during various tasks and sleep. We discovered that the complex located within the frontoparietal network remained constant regardless of task content, while the regional distribution of the complex collapsed in the initial stages of sleep. Furthermore, Φ measures decreased as sleep progressed under limited analysis conditions. These findings align with predictions made by the Integrated Information Theory and support its postulates.

Multi-hierarchy Network Configuration Can Predict Brain States and Performance

The brain is a hierarchical modular organization that varies across functional states. Network configuration can better reveal network organization patterns. However, the multi-hierarchy network configuration remains unknown. Here, we propose an eigenmodal decomposition approach to detect modules at multi-hierarchy, which can identify higher-layer potential submodules and is consistent with the brain hierarchical structure. We defined three metrics: node configuration matrix, combinability, and separability. Node configuration matrix represents network configuration changes between layers. Separability reflects network configuration from global to local, whereas combinability shows network configuration from local to global. First, we created a random network to verify the feasibility of the method. Results show that separability of real networks is larger than that of random networks, whereas combinability is smaller than random networks. Then, we analyzed a large data set incorporating fMRI data from resting and seven distinct tasking conditions. Experiment results demonstrates the high similarity in node configuration matrices for different task conditions, whereas the tasking states have less separability and greater combinability between modules compared with the resting state. Furthermore, the ability of brain network configuration can predict brain states and cognition performance. Crucially, derived from tasks are highlighted with greater power than resting, showing that task-induced attributes have a greater ability to reveal individual differences. Together, our study provides novel perspectives for analyzing the organization structure of complex brain networks at multi-hierarchy, gives new insights to further unravel the working mechanisms of the brain, and adds new evidence for tasking states to better characterize and predict behavioral traits.

RF coil design strategies for improving SNR at the ultrahigh magnetic field of 10.5 Tesla

Purpose

To develop multichannel transmit and receive arrays towards capturing the ultimate-intrinsic-SNR (uiSNR) at 10.5 Tesla (T) and to demonstrate the feasibility and potential of whole-brain, high-resolution human brain imaging at this high field strength.

Methods

A dual row 16-channel self-decoupled transmit (Tx) array was converted to a 16Tx/Rx transceiver using custom transmit/receive switches. A 64-channel receive-only (64Rx) array was built to fit into the 16Tx/Rx array. Electromagnetic modeling and experiments were employed to define safe operation limits of the resulting 16Tx/80Rx array and obtain FDA approval for human use.

Results

The 64Rx array alone captured approximately 50% of the central uiSNR at 10.5T while the identical 7T 64Rx array captured ∼76% of uiSNR at this lower field strength. The 16Tx/80Rx configuration brought the fraction of uiSNR captured at 10.5T to levels comparable to the performance of the 64Rx array at 7T. SNR data obtained at the two field strengths with these arrays displayed dependent increases over a large central region. Whole-brain high resolution T 2 * and T 1 weighted anatomical and gradient-recalled echo EPI BOLD fMRI images were obtained at 10.5T for the first time with such an advanced array, illustrating the promise of >10T fields in studying the human brain.

Conclusion

We demonstrated the ability to approach the uiSNR at 10.5T over the human brain with a novel, high channel count array, achieving large SNR gains over 7T, currently the most commonly employed ultrahigh field platform, and demonstrate high resolution and high contrast anatomical and functional imaging at 10.5T.

Personal Income Performance Correlates with Brain Structural Network Modularity but Not Intelligence Quotient

Introduction: This study aims to use diffusion tensor imaging (DTI) in conjunction with brain graph techniques to define brain structural connectivity and investigate its association with personal income (PI) in individuals of various ages and intelligence quotients (IQ). Methods: MRI examinations were performed on 55 male subjects (mean age: 40.1 ± 9.4 years). Graph data and metrics were generated, and DTI images were analyzed using tract-based spatial statistics (TBSS). All subjects underwent the Wechsler Adult Intelligence Scale for a reliable estimation of the full-scale IQ (FSIQ), which includes verbal comprehension index, perceptual reasoning index, working memory index, and processing speed index. The performance score was defined as the monthly PI normalized by the age of the subject. Results: The analysis of global graph metrics showed that modularity correlated positively with performance score (p = 0.003) and negatively with FSIQ (p = 0.04) and processing speed index (p = 0.005). No significant correlations were found between IQ indices and performance scores. Regional analysis of graph metrics showed modularity differences between right and left networks in sub-cortical (p = 0.001) and frontal (p = 0.044) networks. TBSS analysis showed greater axial and mean diffusivities in the high-performance group in correlation with their modular brain organization. Conclusion: This study showed that PI performance is strongly correlated with a modular organization of brain structural connectivity, which implies short and rapid networks, providing automatic and unconscious brain processing. Additionally, the lack of correlation between performance and IQ suggests a reduced role of academic reasoning skills in performance to the advantage of high uncertainty decision-making networks.

A tale of two n-backs: Diverging associations of dorsolateral prefrontal cortex activation with n-back task performance

Background

In studying the neural correlates of working memory (WM) ability via functional magnetic resonance imaging (fMRI) in health and disease, it is relatively uncommon for investigators to report associations between brain activation and measures of task performance. Additionally, how the choice of WM task impacts observed activation-performance relationships is poorly understood. We sought to illustrate the impact of WM task on brain-behavior correlations using two large, publicly available datasets.

Methods

We conducted between-participants analyses of task-based fMRI data from two publicly available datasets: the Human Connectome Project (HCP; n = 866) and the Queensland Twin Imaging (QTIM) Study (n = 459). Participants performed two distinct variations of the n-back WM task with different stimuli, timings, and response paradigms. Associations between brain activation ([2-back - 0-back] contrast) and task performance (2-back % correct) were investigated separately in each dataset, as well as across datasets, within the dorsolateral prefrontal cortex (dlPFC), medial prefrontal cortex, and whole cortex.

Results

Global patterns of activation to task were similar in both datasets. However, opposite associations between activation and task performance were observed in bilateral pre-supplementary motor area and left middle frontal gyrus. Within the dlPFC, HCP participants exhibited a significantly greater activation-performance relationship in bilateral middle frontal gyrus relative to QTIM Study participants.

Conclusions

The observation of diverging activation-performance relationships between two large datasets performing variations of the n-back task serves as a critical reminder for investigators to exercise caution when selecting WM tasks and interpreting neural activation in response to a WM task.

An in vivo Dissection, and Analysis of Socio-Affective Symptoms related to Cerebellum-Midbrain Reward Circuitry in Humans

Emerging research in non-human animals implicates cerebellar projections to the ventral tegmental area (VTA) in appetitive behaviors, but these circuits have not been characterized in humans. Here, we mapped cerebello-VTA white-matter connectivity in humans using probabilistic tractography on diffusion imaging data from the Human Connectome Project. We uncovered the topographical organization of these connections by separately tracking from parcels of cerebellar lobule VI, crus I/II, vermis, paravermis, and cerebrocerebellum. Results revealed that connections from the cerebellum to the VTA predominantly originate in the right hemisphere, interposed nucleus, and paravermal cortex, and terminate mostly ipsilaterally. Paravermal crus I sends the most connections to the VTA compared to other lobules. We discovered a medial-to-lateral gradient of connectivity, such that the medial cerebellum has the highest connectivity with the VTA. Individual differences in microstructure were associated with measures of negative affect and social functioning. By splitting the tracts into quarters, we found that the socio-affective effects were driven by the third quarter of the tract, corresponding to the point at which the fibers leave the deep nuclei. Taken together, we produced detailed maps of cerebello-VTA structural connectivity for the first time in humans and established their relevance for trait differences in socio-affective regulation.

Characterization of pulsations in the brain and cerebrospinal fluid using ultra-high field magnetic resonance imaging

The development of innovative non-invasive neuroimaging methods and biomarkers is critical for studying brain disease. Imaging of cerebrospinal fluid (CSF) pulsatility may inform the brain fluid dynamics involved in clearance of cerebral metabolic waste. In this work, we developed a methodology to characterize the frequency and spatial localization of whole brain CSF pulsations in humans. Using 7 Tesla (T) human magnetic resonance imaging (MRI) and ultrafast echo-planar imaging (EPI), in-vivo images were obtained to capture pulsations of the CSF signal. Physiological data were simultaneously collected and compared with the 7 T MR data. The primary components of signal pulsations were identified using spectral analysis, with the most evident frequency bands identified around 0.3, 1.2, and 2.4 Hz. These pulsations were mapped spatially and temporally onto the MR image domain and temporally onto the physiological measures of electrocardiogram and respiration. We identified peaks in CSF pulsations that were distinct from peaks in grey matter and white matter regions. This methodology may provide novel in vivo biomarkers of disrupted brain fluid dynamics.

Differences in Cortical Activation During Dorsiflexion and Plantarflexion in Chronic Ankle Instability: A Task-fMRI Study

BACKGROUND

Chronic ankle instability is a common sports injury that often presents with increased plantarflexion and restricted dorsiflexion. The cumulative effect of peripheral injuries may induce neuroplasticity in the central nervous system. However, the relationship between dorsiflexion or plantarflexion and the central nervous system in patients with chronic ankle instability remains unknown.

QUESTIONS/PURPOSES

(1) Is there a difference in region and voxel (volume pixel) of cortical activation during plantarflexion and dorsiflexion between patients with chronic ankle instability and a control group with normal ankle function? (2) Is there a correlation between activation of sensorimotor-related brain regions and three clinical measurement scales of ankle function and disease severity in patients with chronic ankle instability?

METHODS

Between December 2020 and May 2022, we treated 400 patients who had chronic ankle instability. Ten percent (40 patients; mean ± standard deviation age 29 ± 7 years; 17 male patients) were randomly selected to participate in this study. We recruited 42 volunteers with normal ankle function (mean age 28 ± 5 years; 21 male participants) matched by age and education level. A total of 2.5% (1 of 40) of patients with bilateral chronic ankle instability and 30% (12 of 40) with left-sided chronic ankle injury did not meet our inclusion criteria and were excluded from the study. The control group underwent MRI with good image quality. Finally, 27 patients with chronic ankle instability (mean age 26 ± 5 years; 10 male patients) and 42 participants with normal ankle function were enrolled. Ankle function and disease severity were assessed using three clinical scales: the Cumberland Ankle Instability Tool, Karlsson-Peterson Ankle Function Score, and the American Orthopedic Foot and Ankle Society Score. A uniplanar and nonweightbearing ankle dorsiflexion-plantarflexion paradigm (a recognized model or pattern) was performed using a short-block design during the functional MRI scan. This experimental design included a series of on-off periods consisting of movement and a rest period. From 15° of plantarflexion to 15° of dorsiflexion, the manipulator allowed 30° of ankle rotation. The cerebral excitability patterns between patients with chronic ankle instability and controls were analyzed using t-tests. We retained voxels with p values less than 0.05 in a voxel-level family-wise error correction. Clusters with voxel numbers greater than 10 were retained. The Cohen d coefficient was used to calculate between-group effect sizes. Spearman analysis was performed to explore the correlation between activation regions and the three clinical assessment scales.

RESULTS

In the patient group, cortical activation was greater during plantarflexion than during dorsiflexion, which was different from that in the control group. The between-group comparison showed that patients with chronic ankle instability had reduced activation in the ipsilateral precuneus (cluster size = 35 voxels [95% CI -0.23 to 0.07]; p < 0.001) during dorsiflexion, whereas during plantarflexion, chronic ankle instability caused increased activation in the ipsilateral superior temporal gyrus (cluster size = 90 voxels [95% CI -0.73 to -0.13]; p < 0.001), precuneus (cluster size = 18 voxels [95% CI -0.56 to -0.19]; p < 0.001), supplementary motor area (cluster size = 57 voxels [95% CI -0.31 to 0.00]; p < 0.001), superior frontal gyrus (cluster size = 43 voxels [95% CI -0.82 to -0.29]; p < 0.001), medial part of the superior frontal gyrus (cluster size = 39 voxels [95% CI 0.41 to 0.78]; p < 0.001), and contralateral postcentral gyrus (cluster size = 100 voxels [95% CI -0.32 to 0.02]; p < 0.001). Patients with chronic ankle instability showed a large effect size compared with controls (Cohen d > 0.8). During plantarflexion, the number of activated voxels in the supplementary motor area had a modest, positive correlation with the Karlsson-Peterson Ankle Function Score (r = 0.52; p = 0.01), and the number of activated voxels in the primary motor cortex (M1) and primary sensory cortex (S1) had a weak, positive correlation with the American Orthopedic Foot and Ankle Society Score in patients with chronic ankle instability (M1: r = 0.45; p = 0.02, S1: r = 0.49; p = 0.01).

CONCLUSION

Compared with volunteers with normal ankle function, patients with chronic ankle instability had increased cortical activation during plantarflexion and decreased cortical activation during dorsiflexion. We analyzed the central neural mechanisms of chronic ankle instability in patients with sports injuries and provided a theoretical basis for the development of new central and peripheral interventions in the future.

CLINICAL RELEVANCE

Because there was a positive correlation between the neural activity in sensorimotor-related regions during plantarflexion and clinical severity, clinicians might one day be able to help patients who have chronic ankle instability with neuromuscular rehabilitation by applying electrical stimulation to specific targets (such as S1M1 and the supplementary motor area) or by increasing activation of sensorimotor neurons through ankle movement.