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.

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.

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.

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.

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.

Cingulo-opercular control network and disused motor circuits joined in standby mode.

Whole-brain resting-state functional MRI (rs-fMRI) during 2 wk of upper-limb casting revealed that disused motor regions became more strongly connected to the cingulo-opercular network (CON), an executive control network that includes regions of the dorsal anterior cingulate cortex (dACC) and insula. Disuse-driven increases in functional connectivity (FC) were specific to the CON and somatomotor networks and did not involve any other networks, such as the salience, frontoparietal, or default mode networks. Censoring and modeling analyses showed that FC increases during casting were mediated by large, spontaneous activity pulses that appeared in the disused motor regions and CON control regions. During limb constraint, disused motor circuits appear to enter a standby mode characterized by spontaneous activity pulses and strengthened connectivity to CON executive control regions.

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.

Default-mode network streams for coupling to language and control systems.

The human brain is organized into large-scale networks identifiable using resting-state functional connectivity (RSFC). These functional networks correspond with broad cognitive domains; for example, the Default-mode network (DMN) is engaged during internally oriented cognition. However, functional networks may contain hierarchical substructures corresponding with more specific cognitive functions. Here, we used individual-specific precision RSFC to test whether network substructures could be identified in 10 healthy human brains. Across all subjects and networks, individualized network subdivisions were more valid-more internally homogeneous and better matching spatial patterns of task activation-than canonical networks. These measures of validity were maximized at a hierarchical scale that contained ∼83 subnetworks across the brain. At this scale, nine DMN subnetworks exhibited topographical similarity across subjects, suggesting that this approach identifies homologous neurobiological circuits across individuals. Some DMN subnetworks matched known features of brain organization corresponding with cognitive functions. Other subnetworks represented separate streams by which DMN couples with other canonical large-scale networks, including language and control networks. Together, this work provides a detailed organizational framework for studying the DMN in individual humans.

Individual-specific functional connectivity of the amygdala: A substrate for precision psychiatry.

The amygdala is central to the pathophysiology of many psychiatric illnesses. An imprecise understanding of how the amygdala fits into the larger network organization of the human brain, however, limits our ability to create models of dysfunction in individual patients to guide personalized treatment. Therefore, we investigated the position of the amygdala and its functional subdivisions within the network organization of the brain in 10 highly sampled individuals (5 h of fMRI data per person). We characterized three functional subdivisions within the amygdala of each individual. We discovered that one subdivision is preferentially correlated with the default mode network; a second is preferentially correlated with the dorsal attention and fronto-parietal networks; and third subdivision does not have any networks to which it is preferentially correlated relative to the other two subdivisions. All three subdivisions are positively correlated with ventral attention and somatomotor networks and negatively correlated with salience and cingulo-opercular networks. These observations were replicated in an independent group dataset of 120 individuals. We also found substantial across-subject variation in the distribution and magnitude of amygdala functional connectivity with the cerebral cortex that related to individual differences in the stereotactic locations both of amygdala subdivisions and of cortical functional brain networks. Finally, using lag analyses, we found consistent temporal ordering of fMRI signals in the cortex relative to amygdala subdivisions. Altogether, this work provides a detailed framework of amygdala-cortical interactions that can be used as a foundation for models relating aberrations in amygdala connectivity to psychiatric symptoms in individual patients.

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.

Trait-like variants in human functional brain networks.

Resting-state functional magnetic resonance imaging (fMRI) has provided converging descriptions of group-level functional brain organization. Recent work has revealed that functional networks identified in individuals contain local features that differ from the group-level description. We define these features as network variants. Building on these studies, we ask whether distributions of network variants reflect stable, trait-like differences in brain organization. Across several datasets of highly-sampled individuals we show that 1) variants are highly stable within individuals, 2) variants are found in characteristic locations and associate with characteristic functional networks across large groups, 3) task-evoked signals in variants demonstrate a link to functional variation, and 4) individuals cluster into subgroups on the basis of variant characteristics that are related to differences in behavior. These results suggest that distributions of network variants may reflect stable, trait-like, functionally relevant individual differences in functional brain organization.

Organization of Propagated Intrinsic Brain Activity in Individual Humans.

Spontaneous infra-slow (<0.1 Hz) fluctuations in functional magnetic resonance imaging (fMRI) signals are temporally correlated within large-scale functional brain networks, motivating their use for mapping systems-level brain organization. However, recent electrophysiological and hemodynamic evidence suggest state-dependent propagation of infra-slow fluctuations, implying a functional role for ongoing infra-slow activity. Crucially, the study of infra-slow temporal lag structure has thus far been limited to large groups, as analyzing propagation delays requires extensive data averaging to overcome sampling variability. Here, we use resting-state fMRI data from 11 extensively-sampled individuals to characterize lag structure at the individual level. In addition to stable individual-specific features, we find spatiotemporal topographies in each subject similar to the group average. Notably, we find a set of early regions that are common to all individuals, are preferentially positioned proximal to multiple functional networks, and overlap with brain regions known to respond to diverse behavioral tasks-altogether consistent with a hypothesized ability to broadly influence cortical excitability. Our findings suggest that, like correlation structure, temporal lag structure is a fundamental organizational property of resting-state infra-slow activity.

Removal of high frequency contamination from motion estimates in single-band fMRI saves data without biasing functional connectivity.

Denoising fMRI data requires assessment of frame-to-frame head motion and removal of the biases motion introduces. This is usually done through analysis of the parameters calculated during retrospective head motion correction (i.e., 'motion' parameters). However, it is increasingly recognized that respiration introduces factitious head motion via perturbations of the main (B0) field. This effect appears as higher-frequency fluctuations in the motion parameters (>0.1 â€‹Hz, here referred to as 'HF-motion'), primarily in the phase-encoding direction. This periodicity can sometimes be obscured in standard single-band fMRI (TR 2.0-2.5 â€‹s) due to aliasing. Here we examined (1) how prevalent HF-motion effects are in seven single-band datasets with TR from 2.0 to 2.5 â€‹s and (2) how HF-motion affects functional connectivity. We demonstrate that HF-motion is more common in older adults, those with higher body mass index, and those with lower cardiorespiratory fitness. We propose a low-pass filtering approach to remove the contamination of high frequency effects from motion summary measures, such as framewise displacement (FD). We demonstrate that in most datasets this filtering approach saves a substantial amount of data from FD-based frame censoring, while at the same time reducing motion biases in functional connectivity measures. These findings suggest that filtering motion parameters is an effective way to improve the fidelity of head motion estimates, even in single band datasets. Particularly large data savings may accrue in datasets acquired in older and less fit participants.

Correction of respiratory artifacts in MRI head motion estimates.

Head motion represents one of the greatest technical obstacles in magnetic resonance imaging (MRI) of the human brain. Accurate detection of artifacts induced by head motion requires precise estimation of movement. However, head motion estimates may be corrupted by artifacts due to magnetic main field fluctuations generated by body motion. In the current report, we examine head motion estimation in multiband resting state functional connectivity MRI (rs-fcMRI) data from the Adolescent Brain and Cognitive Development (ABCD) Study and comparison 'single-shot' datasets. We show that respirations contaminate movement estimates in functional MRI and that respiration generates apparent head motion not associated with functional MRI quality reductions. We have developed a novel approach using a band-stop filter that accurately removes these respiratory effects from motion estimates. Subsequently, we demonstrate that utilizing a band-stop filter improves post-processing fMRI data quality. Lastly, we demonstrate the real-time implementation of motion estimate filtering in our FIRMM (Framewise Integrated Real-Time MRI Monitoring) software package.

High-fidelity mapping of repetition-related changes in the parietal memory network.

fMRI studies of human memory have identified a "parietal memory network" (PMN) that displays distinct responses to novel and familiar stimuli, typically deactivating during initial encoding but robustly activating during retrieval. The small size of PMN regions, combined with their proximity to the neighboring default mode network, makes a targeted assessment of their responses in highly sampled subjects important for understanding information processing within the network. Here, we describe an experiment in which participants made semantic decisions about repeatedly-presented stimuli, assessing PMN BOLD responses as items transitioned from experimentally novel to repeated. Data are from the highly-sampled subjects in the Midnight Scan Club dataset, enabling a characterization of BOLD responses at both the group and single-subject level. Across all analyses, PMN regions deactivated in response to novel stimuli and displayed changes in BOLD activity across presentations, but did not significantly activate to repeated items. Results support only a portion of initially hypothesized effects, in particular suggesting that novelty-related deactivations may be less susceptible to attentional/task manipulations than are repetition-related activations within the network. This in turn suggests that novelty and familiarity may be processed as separable entities within the PMN.

Reward-related regions form a preferentially coupled system at rest.

Neuroimaging studies have implicated a set of striatal and orbitofrontal cortex (OFC) regions that are commonly activated during reward processing tasks. Resting-state functional connectivity (RSFC) studies have demonstrated that the human brain is organized into several functional systems that show strong temporal coherence in the absence of goal-directed tasks. Here we use seed-based and graph-theory RSFC approaches to characterize the systems-level organization of putative reward regions of at rest. Peaks of connectivity from seed-based RSFC patterns for the nucleus accumbens (NAcc) and orbitofrontal cortex (OFC) were used to identify candidate reward regions which were merged with a previously used set of regions (Power et al., 2011). Graph-theory was then used to determine system-level membership for all regions. Several regions previously implicated in reward-processing (NAcc, lateral and medial OFC, and ventromedial prefrontal cortex) comprised a distinct, preferentially coupled system. This RSFC system is stable across a range of connectivity thresholds and shares strong overlap with meta-analyses of task-based reward studies. This reward system shares between-system connectivity with systems implicated in cognitive control and self-regulation, including the fronto-parietal, cingulo-opercular, and default systems. Differences may exist in the pathways through which control systems interact with reward system components. Whereas NAcc is functionally connected to cingulo-opercular and default systems, OFC regions show stronger connectivity with the fronto-parietal system. We propose that future work may be able to interrogate group or individual differences in connectivity profiles using the regions delineated in this work to explore potential relationships to appetitive behaviors, self-regulation failure, and addiction.

Individualized Connectome-Targeted Transcranial Magnetic Stimulation for Neuropsychiatric Sequelae of Repetitive Traumatic Brain Injury in a Retired NFL Player.

OBJECTIVE: The recent advent of individualized resting-state network mapping (RSNM) has revealed substantial interindividual variability in anatomical localization of brain networks identified by using resting-state functional MRI (rsfMRI). RSNM enables personalized targeting of focal neuromodulation techniques such as repetitive transcranial magnetic stimulation (rTMS). rTMS is believed to exert antidepressant efficacy by modulating connectivity between the stimulation site, the default mode network (DMN), and the subgenual anterior cingulate cortex (sgACC). Personalized rTMS may be particularly useful after repetitive traumatic brain injury (TBI), which is associated with neurodegenerative tauopathy in medial temporal limbic structures. These degenerative changes are believed to be related to treatment-resistant neurobehavioral disturbances observed in many retired athletes. METHODS: The authors describe a case in which RSNM was successfully used to target rTMS to treat these neuropsychiatric disturbances in a retired NFL defensive lineman whose symptoms were not responsive to conventional treatments. RSNM was used to identify left-right dorsolateral prefrontal rTMS targets with maximal difference between dorsal attention network and DMN correlations. These targets were spatially distinct from those identified by prior methods. Twenty sessions of left-sided excitatory and right-sided inhibitory rTMS were administered at these targets. RESULTS: Treatment led to improvement in Montgomery-Åsberg Depression Rating Scale (72%), cognitive testing, and headache scales scores. Compared with healthy individuals and subjects with TBI-associated depression, baseline rsfMRI revealed substantially elevated DMN connectivity with the medial temporal lobe (MTL). Serial rsfMRI scans revealed gradual improvement in MTL-DMN connectivity and stimulation site connectivity with sgACC. CONCLUSIONS: These results highlight the possibility of individualized neuromodulation and biomarker-based monitoring for neuropsychiatric sequelae of repetitive TBI.

Local-Global Parcellation of the Human Cerebral Cortex from Intrinsic Functional Connectivity MRI.

A central goal in systems neuroscience is the parcellation of the cerebral cortex into discrete neurobiological "atoms". Resting-state functional magnetic resonance imaging (rs-fMRI) offers the possibility of in vivo human cortical parcellation. Almost all previous parcellations relied on 1 of 2 approaches. The local gradient approach detects abrupt transitions in functional connectivity patterns. These transitions potentially reflect cortical areal boundaries defined by histology or visuotopic fMRI. By contrast, the global similarity approach clusters similar functional connectivity patterns regardless of spatial proximity, resulting in parcels with homogeneous (similar) rs-fMRI signals. Here, we propose a gradient-weighted Markov Random Field (gwMRF) model integrating local gradient and global similarity approaches. Using task-fMRI and rs-fMRI across diverse acquisition protocols, we found gwMRF parcellations to be more homogeneous than 4 previously published parcellations. Furthermore, gwMRF parcellations agreed with the boundaries of certain cortical areas defined using histology and visuotopic fMRI. Some parcels captured subareal (somatotopic and visuotopic) features that likely reflect distinct computational units within known cortical areas. These results suggest that gwMRF parcellations reveal neurobiologically meaningful features of brain organization and are potentially useful for future applications requiring dimensionality reduction of voxel-wise fMRI data. Multiresolution parcellations generated from 1489 participants are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Schaefer2018_LocalGlobal).