MRI economics: Balancing sample size and scan duration in brain wide association studies.

A pervasive dilemma in neuroimaging is whether to prioritize sample size or scan duration given fixed resources. Here, we systematically investigate this trade-off in the context of brain-wide association studies (BWAS) using resting-state functional magnetic resonance imaging (fMRI). We find that total scan duration (sample size × scan duration per participant) robustly explains individual-level phenotypic prediction accuracy via a logarithmic model, suggesting that sample size and scan duration are broadly interchangeable. The returns of scan duration eventually diminish relative to sample size, which we explain with principled theoretical derivations. When accounting for fixed costs associated with each participant (e.g., recruitment, non-imaging measures), we find that prediction accuracy in small-scale BWAS might benefit from much longer scan durations (>50 min) than typically assumed. Most existing large-scale studies might also have benefited from smaller sample sizes with longer scan durations. Both logarithmic and theoretical models of the relationships among sample size, scan duration and prediction accuracy explain well-predicted phenotypes better than poorly-predicted phenotypes. The logarithmic and theoretical models are also undermined by individual differences in brain states. These results replicate across phenotypic domains (e.g., cognition and mental health) from two large-scale datasets with different algorithms and metrics. Overall, our study emphasizes the importance of scan time, which is ignored in standard power calculations. Standard power calculations inevitably maximize sample size at the expense of scan duration. The resulting prediction accuracies are likely lower than would be produced with alternate designs, thus impeding scientific discovery. Our empirically informed reference is available for future study design: WEB_APPLICATION_LINK.

Functional parcellation of the neonatal cortical surface.

The cerebral cortex is organized into distinct but interconnected cortical areas, which can be defined by abrupt differences in patterns of resting state functional connectivity (FC) across the cortical surface. Such parcellations of the cortex have been derived in adults and older infants, but there is no widely used surface parcellation available for the neonatal brain. Here, we first demonstrate that existing parcellations, including surface-based parcels derived from older samples as well as volume-based neonatal parcels, are a poor fit for neonatal surface data. We next derive a set of 283 cortical surface parcels from a sample of n = 261 neonates. These parcels have highly homogenous FC patterns and are validated using three external neonatal datasets. The Infomap algorithm is used to assign functional network identities to each parcel, and derived networks are consistent with prior work in neonates. The proposed parcellation may represent neonatal cortical areas and provides a powerful tool for neonatal neuroimaging studies.

The brainstem's red nucleus was evolutionarily upgraded to support goal-directed action.

The red nucleus is a large brainstem structure that coordinates limb movement for locomotion in quadrupedal animals (Basile et al., 2021). The humans red nucleus has a different pattern of anatomical connectivity compared to quadrupeds, suggesting a unique purpose (Hatschek, 1907). Previously the function of the human red nucleus remained unclear at least partly due to methodological limitations with brainstem functional neuroimaging (Sclocco et al., 2018). Here, we used our most advanced resting-state functional connectivity (RSFC) based precision functional mapping (PFM) in highly sampled individuals (n = 5) and large group-averaged datasets (combined N ~ 45,000), to precisely examine red nucleus functional connectivity. Notably, red nucleus functional connectivity to motor-effector networks (somatomotor hand, foot, and mouth) was minimal. Instead, red nucleus functional connectivity along the central sulcus was specific to regions of the recently discovered somato-cognitive action network (SCAN; (Gordon et al., 2023)). Outside of primary motor cortex, red nucleus connectivity was strongest to the cingulo-opercular (CON) and salience networks, involved in action/cognitive control (Dosenbach et al., 2007; Newbold et al., 2021) and reward/motivated behavior (Seeley, 2019), respectively. Functional connectivity to these two networks was organized into discrete dorsal-medial and ventral-lateral zones. Red nucleus functional connectivity to the thalamus recapitulated known structural connectivity of the dento-rubral thalamic tract (DRTT) and could prove clinically useful in functionally targeting the ventral intermediate (VIM) nucleus. In total, our results indicate that far from being a 'motor' structure, the red nucleus is better understood as a brainstem nucleus for implementing goal-directed behavior, integrating behavioral valence and action plans.

Disuse-driven plasticity in the human thalamus and putamen.

Motor adaptation in cortico-striato-thalamo-cortical loops has been studied mainly in animals using invasive electrophysiology. Here, we leverage functional neuroimaging in humans to study motor circuit plasticity in the human subcortex. We employed an experimental paradigm that combined two weeks of upper-extremity immobilization with daily resting-state and motor task fMRI before, during, and after the casting period. We previously showed that limb disuse leads to decreased functional connectivity (FC) of the contralateral somatomotor cortex (SM1) with the ipsilateral somatomotor cortex, increased FC with the cingulo-opercular network (CON) as well as the emergence of high amplitude, fMRI signal pulses localized in the contralateral SM1, supplementary motor area and the cerebellum. From our prior observations, it remains unclear whether the disuse plasticity affects the thalamus and striatum. We extended our analysis to include these subcortical regions and found that both exhibit strengthened cortical FC and spontaneous fMRI signal pulses induced by limb disuse. The dorsal posterior putamen and the central thalamus, mainly CM, VLP and VIM nuclei, showed disuse pulses and FC changes that lined up with fmri task activations from the Human connectome project motor system localizer, acquired before casting for each participant. Our findings provide a novel understanding of the role of the cortico-striato-thalamo-cortical loops in human motor plasticity and a potential link with the physiology of sleep regulation. Additionally, similarities with FC observation from Parkinson Disease (PD) questions a pathophysiological link with limb disuse.

Framewise multi-echo distortion correction for superior functional MRI.

Functional MRI (fMRI) data are severely distorted by magnetic field (B0) inhomogeneities which currently must be corrected using separately acquired field map data. However, changes in the head position of a scanning participant across fMRI frames can cause changes in the B0 field, preventing accurate correction of geometric distortions. Additionally, field maps can be corrupted by movement during their acquisition, preventing distortion correction altogether. In this study, we use phase information from multi-echo (ME) fMRI data to dynamically sample distortion due to fluctuating B0 field inhomogeneity across frames by acquiring multiple echoes during a single EPI readout. Our distortion correction approach, MEDIC (Multi-Echo DIstortion Correction), accurately estimates B0 related distortions for each frame of multi-echo fMRI data. Here, we demonstrate that MEDIC's framewise distortion correction produces improved alignment to anatomy and decreases the impact of head motion on resting-state functional connectivity (RSFC) maps, in higher motion data, when compared to the prior gold standard approach (i.e., TOPUP). Enhanced framewise distortion correction with MEDIC, without the requirement for field map collection, furthers the advantage of multi-echo over single-echo fMRI.

Multi-echo Acquisition and Thermal Denoising Advances Infant Precision Functional Imaging.

The characterization of individual functional brain organization with Precision Functional Mapping has provided important insights in recent years in adults. However, little is known about the ontogeny of inter-individual differences in brain functional organization during human development, but precise characterization of systems organization during periods of high plasticity might be most influential towards discoveries promoting lifelong health. Collecting and analyzing precision fMRI data during early development has unique challenges and emphasizes the importance of novel methods to improve data acquisition, processing, and analysis strategies in infant samples. Here, we investigate the applicability of two such methods from adult MRI research, multi-echo (ME) data acquisition and thermal noise removal with Noise reduction with distribution corrected principal component analysis (NORDIC), in precision fMRI data from three newborn infants. Compared to an adult example subject, T2* relaxation times calculated from ME data in infants were longer and more variable across the brain, pointing towards ME acquisition being a promising tool for optimizing developmental fMRI. The application of thermal denoising via NORDIC increased tSNR and the overall strength of functional connections as well as the split-half reliability of functional connectivity matrices in infant ME data. While our findings related to NORDIC denoising are coherent with the adult literature and ME data acquisition showed high promise, its application in developmental samples needs further investigation. The present work reveals gaps in our understanding of the best techniques for developmental brain imaging and highlights the need for further developmentally-specific methodological advances and optimizations, towards precision functional imaging in infants.

Functional parcellation of the neonatal brain.

The cerebral cortex is organized into distinct but interconnected cortical areas, which can be defined by abrupt differences in patterns of resting state functional connectivity (FC) across the cortical surface. Such parcellations of the cortex have been derived in adults and older infants, but there is no widely used surface parcellation available for the neonatal brain. Here, we first demonstrate that adult- and older infant-derived parcels are a poor fit with neonatal data, emphasizing the need for neonatal-specific parcels. We next derive a set of 283 cortical surface parcels from a sample of n=261 neonates. These parcels have highly homogenous FC patterns and are validated using three external neonatal datasets. The Infomap algorithm is used to assign functional network identities to each parcel, and derived networks are consistent with prior work in neonates. The proposed parcellation may represent neonatal cortical areas and provides a powerful tool for neonatal neuroimaging studies.

Substructure of the brain's Cingulo-Opercular network.

The Cingulo-Opercular network (CON) is an executive network of the human brain that regulates actions. CON is composed of many widely distributed cortical regions that are involved in top-down control over both lower-level (i.e., motor) and higher-level (i.e., cognitive) functions, as well as in processing of painful stimuli. Given the topographical and functional heterogeneity of the CON, we investigated whether subnetworks within the CON support separable aspects of action control. Using precision functional mapping (PFM) in 15 participants with > 5 hours of resting state functional connectivity (RSFC) and task data, we identified three anatomically and functionally distinct CON subnetworks within each individual. These three distinct subnetworks were linked to Decisions, Actions, and Feedback (including pain processing), respectively, in convergence with a meta-analytic task database. These Decision, Action and Feedback subnetworks represent pathways by which the brain establishes top-down goals, transforms those goals into actions, implemented as movements, and processes critical action feedback such as pain.

Structure-Function Coupling in Highly Sampled Individual Brains.

Structural connections (SC) between distant regions of the brain support synchronized function known as functional connectivity (FC) and give rise to the large-scale brain networks that enable cognition and behavior. Understanding how SC enables FC is important to understand how injuries to structural connections may alter brain function and cognition. Previous work evaluating whole-brain SC-FC relationships showed that SC explained FC well in unimodal visual and motor areas, but only weakly in association areas, suggesting a unimodal-heteromodal gradient organization of SC-FC coupling. However, this work was conducted in group-averaged SC/FC data. Thus, it could not account for inter-individual variability in the locations of cortical areas and white matter tracts. We evaluated the correspondence of SC and FC within three highly sampled healthy participants. For each participant, we collected 78 minutes of diffusion-weighted MRI for SC and 360 minutes of resting state fMRI for FC. We found that FC was best explained by SC in visual and motor systems, as well as in anterior and posterior cingulate regions. A unimodal-to-heteromodal gradient could not fully explain SC-FC coupling. We conclude that the SC-FC coupling of the anterior-posterior cingulate circuit is more similar to unimodal areas than to heteromodal areas. SIGNIFICANCE STATEMENT: Structural connections between distant regions of the human brain support networked function that enables cognition and behavior. Improving our understanding of how structure enables function could allow better insight into how brain disconnection injuries impair brain function.Previous work using neuroimaging suggested that structure-function relationships vary systematically across the brain, with structure better explaining function in basic visual/motor areas than in higher-order areas. However, this work was conducted in group-averaged data, which may obscure details of individual-specific structure-function relationships.Using individual-specific densely sampled neuroimaging data, we found that in addition to visual/motor regions, structure strongly predicts function in specific circuits of the higher-order cingulate gyrus. The cingulate's structure-function relationship suggests that its organization may be unique among higher-order cortical regions.

Psilocybin desynchronizes brain networks.

1The relationship between the acute effects of psychedelics and their persisting neurobiological and psychological effects is poorly understood. Here, we tracked brain changes with longitudinal precision functional mapping in healthy adults before, during, and for up to 3 weeks after oral psilocybin and methylphenidate (17 MRI visits per participant) and again 6+ months later. Psilocybin disrupted connectivity across cortical networks and subcortical structures, producing more than 3-fold greater acute changes in functional networks than methylphenidate. These changes were driven by desynchronization of brain activity across spatial scales (area, network, whole brain). Psilocybin-driven desynchronization was observed across association cortex but strongest in the default mode network (DMN), which is connected to the anterior hippocampus and thought to create our sense of self. Performing a perceptual task reduced psilocybin-induced network changes, suggesting a neurobiological basis for grounding, connecting with physical reality during psychedelic therapy. The acute brain effects of psilocybin are consistent with distortions of space-time and the self. Psilocybin induced persistent decrease in functional connectivity between the anterior hippocampus and cortex (and DMN in particular), lasting for weeks but normalizing after 6 months. Persistent suppression of hippocampal-DMN connectivity represents a candidate neuroanatomical and mechanistic correlate for psilocybin's pro-plasticity and anti-depressant effects.

Precision Neuroimaging for Localization-Related Psychiatry.

This Viewpoint describes how precision functional mapping may be helpful for associating neuroanatomical regions with specific psychiatric disorders.

Thalamo-cortical and cerebello-cortical functional connectivity in development.

The thalamus is a critical relay center for neural pathways involving sensory, motor, and cognitive functions, including cortico-striato-thalamo-cortical and cortico-ponto-cerebello-thalamo-cortical loops. Despite the importance of these circuits, their development has been understudied. One way to investigate these pathways in human development in vivo is with functional connectivity MRI, yet few studies have examined thalamo-cortical and cerebello-cortical functional connectivity in development. Here, we used resting-state functional connectivity to measure functional connectivity in the thalamus and cerebellum with previously defined cortical functional networks in 2 separate data sets of children (7-12 years old) and adults (19-40 years old). In both data sets, we found stronger functional connectivity between the ventral thalamus and the somatomotor face cortical functional network in children compared with adults, extending previous cortico-striatal functional connectivity findings. In addition, there was more cortical network integration (i.e. strongest functional connectivity with multiple networks) in the thalamus in children than in adults. We found no developmental differences in cerebello-cortical functional connectivity. Together, these results suggest different maturation patterns in cortico-striato-thalamo-cortical and cortico-ponto-cerebellar-thalamo-cortical pathways.

A somato-cognitive action network alternates with effector regions in motor cortex.

Motor cortex (M1) has been thought to form a continuous somatotopic homunculus extending down the precentral gyrus from foot to face representations1,2, despite evidence for concentric functional zones3 and maps of complex actions4. Here, using precision functional magnetic resonance imaging (fMRI) methods, we find that the classic homunculus is interrupted by regions with distinct connectivity, structure and function, alternating with effector-specific (foot, hand and mouth) areas. These inter-effector regions exhibit decreased cortical thickness and strong functional connectivity to each other, as well as to the cingulo-opercular network (CON), critical for action5 and physiological control6, arousal7, errors8 and pain9. This interdigitation of action control-linked and motor effector regions was verified in the three largest fMRI datasets. Macaque and pediatric (newborn, infant and child) precision fMRI suggested cross-species homologues and developmental precursors of the inter-effector system. A battery of motor and action fMRI tasks documented concentric effector somatotopies, separated by the CON-linked inter-effector regions. The inter-effectors lacked movement specificity and co-activated during action planning (coordination of hands and feet) and axial body movement (such as of the abdomen or eyebrows). These results, together with previous studies demonstrating stimulation-evoked complex actions4 and connectivity to internal organs10 such as the adrenal medulla, suggest that M1 is punctuated by a system for whole-body action planning, the somato-cognitive action network (SCAN). In M1, two parallel systems intertwine, forming an integrate-isolate pattern: effector-specific regions (foot, hand and mouth) for isolating fine motor control and the SCAN for integrating goals, physiology and body movement.

Using synthetic MR images for distortion correction.

Functional MRI (fMRI) data acquired using echo-planar imaging (EPI) are highly distorted by magnetic field inhomogeneities. Distortion and differences in image contrast between EPI and T1-weighted and T2-weighted (T1w/T2w) images makes their alignment a challenge. Typically, field map data are used to correct EPI distortions. Alignments achieved with field maps can vary greatly and depends on the quality of field map data. However, many public datasets lack field map data entirely. Additionally, reliable field map data is often difficult to acquire in high-motion pediatric or developmental cohorts. To address this, we developed Synth, a software package for distortion correction and cross-modal image registration that does not require field map data. Synth combines information from T1w and T2w anatomical images to construct an idealized undistorted synthetic image with similar contrast properties to EPI data. This synthetic image acts as an effective reference for individual-specific distortion correction. Using pediatric (ABCD: Adolescent Brain Cognitive Development) and adult (MSC: Midnight Scan Club; HCP: Human Connectome Project) data, we demonstrate that Synth performs comparably to field map distortion correction approaches, and often outperforms them. Field map-less distortion correction with Synth allows accurate and precise registration of fMRI data with missing or corrupted field map information.

Accuracy and reliability of diffusion imaging models.

Diffusion imaging aims to non-invasively characterize the anatomy and integrity of the brain's white matter fibers. We evaluated the accuracy and reliability of commonly used diffusion imaging methods as a function of data quantity and analysis method, using both simulations and highly sampled individual-specific data (927-1442 diffusion weighted images [DWIs] per individual). Diffusion imaging methods that allow for crossing fibers (FSL's BedpostX [BPX], DSI Studio's Constant Solid Angle Q-Ball Imaging [CSA-QBI], MRtrix3's Constrained Spherical Deconvolution [CSD]) estimated excess fibers when insufficient data were present and/or when the data did not match the model priors. To reduce such overfitting, we developed a novel Bayesian Multi-tensor Model-selection (BaMM) method and applied it to the popular ball-and-stick model used in BedpostX within the FSL software package. BaMM was robust to overfitting and showed high reliability and the relatively best crossing-fiber accuracy with increasing amounts of diffusion data. Thus, sufficient data and an overfitting resistant analysis method enhance precision diffusion imaging. For potential clinical applications of diffusion imaging, such as neurosurgical planning and deep brain stimulation (DBS), the quantities of data required to achieve diffusion imaging reliability are lower than those needed for functional MRI.

Reproducible brain-wide association studies require thousands of individuals.

Magnetic resonance imaging (MRI) has transformed our understanding of the human brain through well-replicated mapping of abilities to specific structures (for example, lesion studies) and functions1-3 (for example, task functional MRI (fMRI)). Mental health research and care have yet to realize similar advances from MRI. A primary challenge has been replicating associations between inter-individual differences in brain structure or function and complex cognitive or mental health phenotypes (brain-wide association studies (BWAS)). Such BWAS have typically relied on sample sizes appropriate for classical brain mapping4 (the median neuroimaging study sample size is about 25), but potentially too small for capturing reproducible brain-behavioural phenotype associations5,6. Here we used three of the largest neuroimaging datasets currently available-with a total sample size of around 50,000 individuals-to quantify BWAS effect sizes and reproducibility as a function of sample size. BWAS associations were smaller than previously thought, resulting in statistically underpowered studies, inflated effect sizes and replication failures at typical sample sizes. As sample sizes grew into the thousands, replication rates began to improve and effect size inflation decreased. More robust BWAS effects were detected for functional MRI (versus structural), cognitive tests (versus mental health questionnaires) and multivariate methods (versus univariate). Smaller than expected brain-phenotype associations and variability across population subsamples can explain widespread BWAS replication failures. In contrast to non-BWAS approaches with larger effects (for example, lesions, interventions and within-person), BWAS reproducibility requires samples with thousands of individuals.

Individualized Functional Subnetworks Connect Human Striatum and Frontal Cortex.

The striatum and cerebral cortex are interconnected via multiple recurrent loops that play a major role in many neuropsychiatric conditions. Primate corticostriatal connections can be precisely mapped using invasive tract-tracing. However, noninvasive human research has not mapped these connections with anatomical precision, limited in part by the practice of averaging neuroimaging data across individuals. Here we utilized highly sampled resting-state functional connectivity MRI for individual-specific precision functional mapping (PFM) of corticostriatal connections. We identified ten individual-specific subnetworks linking cortex-predominately frontal cortex-to striatum, most of which converged with nonhuman primate tract-tracing work. These included separable connections between nucleus accumbens core/shell and orbitofrontal/medial frontal gyrus; between anterior striatum and dorsomedial prefrontal cortex; between dorsal caudate and lateral prefrontal cortex; and between middle/posterior putamen and supplementary motor/primary motor cortex. Two subnetworks that did not converge with nonhuman primates were connected to cortical regions associated with human language function. Thus, precision subnetworks identify detailed, individual-specific, neurobiologically plausible corticostriatal connectivity that includes human-specific language networks.

Parallel hippocampal-parietal circuits for self- and goal-oriented processing.

The hippocampus is critically important for a diverse range of cognitive processes, such as episodic memory, prospective memory, affective processing, and spatial navigation. Using individual-specific precision functional mapping of resting-state functional MRI data, we found the anterior hippocampus (head and body) to be preferentially functionally connected to the default mode network (DMN), as expected. The hippocampal tail, however, was strongly preferentially functionally connected to the parietal memory network (PMN), which supports goal-oriented cognition and stimulus recognition. This anterior-posterior dichotomy of resting-state functional connectivity was well-matched by differences in task deactivations and anatomical segmentations of the hippocampus. Task deactivations were localized to the hippocampal head and body (DMN), relatively sparing the tail (PMN). The functional dichotomization of the hippocampus into anterior DMN-connected and posterior PMN-connected parcels suggests parallel but distinct circuits between the hippocampus and medial parietal cortex for self- versus goal-oriented processing.