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.

An executive-functions-based reading training enhances sensory-motor systems integration during reading fluency in children with dyslexia.

The Simple View of Reading model suggests that intact language processing and word decoding lead to proficient reading comprehension, with recent studies pointing at executive functions as an important component contributing to reading proficiency. Here, we aimed to determine the underlying mechanism(s) for these changes. Participants include 120 8- to 12-year-old children (n = 55 with dyslexia, n = 65 typical readers) trained on an executive functions-based reading program, including pre/postfunctional MRI and behavioral data collection. Across groups, improved word reading was related to stronger functional connections within executive functions and sensory networks. In children with dyslexia, faster and more accurate word reading was related to stronger functional connections within and between sensory networks. These results suggest greater synchronization of brain systems after the intervention, consistent with the "neural noise" hypothesis in children with dyslexia and support the consideration of including executive functions as part of the Simple View of Reading model.

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.

Maturation of large-scale brain systems over the first month of life.

The period immediately after birth is a critical developmental window, capturing rapid maturation of brain structure and a child's earliest experiences. Large-scale brain systems are present at delivery, but how these brain systems mature during this narrow window (i.e. first weeks of life) marked by heightened neuroplasticity remains uncharted. Using multivariate pattern classification techniques and functional connectivity magnetic resonance imaging, we detected robust differences in brain systems related to age in newborns (n = 262; R2 = 0.51). Development over the first month of life occurred brain-wide, but differed and was more pronounced in brain systems previously characterized as developing early (i.e. sensorimotor networks) than in those characterized as developing late (i.e. association networks). The cingulo-opercular network was the only exception to this organizing principle, illuminating its early role in brain development. This study represents a step towards a normative brain "growth curve" that could be used to identify atypical brain maturation in infancy.

Network-specific selectivity of functional connections in the neonatal brain.

The adult human brain is organized into functional brain networks, groups of functionally connected segregated brain regions. A key feature of adult functional networks is long-range selectivity, the property that spatially distant regions from the same network have higher functional connectivity than spatially distant regions from different networks. Although it is critical to establish the status of functional networks and long-range selectivity during the neonatal period as a foundation for typical and atypical brain development, prior work in this area has been mixed. Although some studies report distributed adult-like networks, other studies suggest that neonatal networks are immature and consist primarily of spatially isolated regions. Using a large sample of neonates (n = 262), we demonstrate that neonates have long-range selective functional connections for the default mode, fronto-parietal, and dorsal attention networks. An adult-like pattern of functional brain networks is evident in neonates when network-detection algorithms are tuned to these long-range connections, when using surface-based registration (versus volume-based registration), and as per-subject data quantity increases. These results help clarify factors that have led to prior mixed results, establish that key adult-like functional network features are evident in neonates, and provide a foundation for studies of typical and atypical brain development.

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.

Specialized late cingulo-opercular network activation elucidates the mechanisms underlying decisions about ambiguity.

Cortical task control networks, including the cingulo-opercular (CO) network play a key role in decision-making across a variety of functional domains. In particular, the CO network functions in a performance reporting capacity that supports successful task performance, especially in response to errors and ambiguity. In two studies testing the contribution of the CO network to ambiguity processing, we presented a valence bias task in which masked clearly and ambiguously valenced emotional expressions were slowly revealed over several seconds. This slow reveal task design provides a window into the decision-making mechanisms as they unfold over the course of a trial. In the main study, the slow reveal task was administered to 32 young adults in the fMRI environment and BOLD time courses were extracted from regions of interest in three control networks. In a follow-up study, the task was administered to a larger, online sample (n = 81) using a more extended slow reveal design with additional unmasking frames. Positive judgments of surprised faces were uniquely accompanied by slower response times and strong, late activation in the CO network. These results support the initial negativity hypothesis, which posits that the default response to ambiguity is negative and positive judgments are associated with a more effortful controlled process, and additionally suggest that this controlled process is mediated by the CO network. Moreover, ambiguous trials were characterized by a second CO response at the end of the trial, firmly placing CO function late in the decision-making process.

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.

BOLD cofluctuation 'events' are predicted from static functional connectivity.

Recent work identified single time points ("events") of high regional cofluctuation in functional Magnetic Resonance Imaging (fMRI) which contain more large-scale brain network information than other, low cofluctuation time points. This suggested that events might be a discrete, temporally sparse signal which drives functional connectivity (FC) over the timeseries. However, a different, not yet explored possibility is that network information differences between time points are driven by sampling variability on a constant, static, noisy signal. Using a combination of real and simulated data, we examined the relationship between cofluctuation and network structure and asked if this relationship was unique, or if it could arise from sampling variability alone. First, we show that events are not discrete - there is a gradually increasing relationship between network structure and cofluctuation; ∼50% of samples show very strong network structure. Second, using simulations we show that this relationship is predicted from sampling variability on static FC. Finally, we show that randomly selected points can capture network structure about as well as events, largely because of their temporal spacing. Together, these results suggest that, while events exhibit particularly strong representations of static FC, there is little evidence that events are unique timepoints that drive FC structure. Instead, a parsimonious explanation for the data is that events arise from a single static, but noisy, FC structure.

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.

Probabilistic mapping of human functional brain networks identifies regions of high group consensus.

Many recent developments surrounding the functional network organization of the human brain have focused on data that have been averaged across groups of individuals. While such group-level approaches have shed considerable light on the brain's large-scale distributed systems, they conceal individual differences in network organization, which recent work has demonstrated to be common and widespread. This individual variability produces noise in group analyses, which may average together regions that are part of different functional systems across participants, limiting interpretability. However, cost and feasibility constraints may limit the possibility for individual-level mapping within studies. Here our goal was to leverage information about individual-level brain organization to probabilistically map common functional systems and identify locations of high inter-subject consensus for use in group analyses. We probabilistically mapped 14 functional networks in multiple datasets with relatively high amounts of data. All networks show "core" (high-probability) regions, but differ from one another in the extent of their higher-variability components. These patterns replicate well across four datasets with different participants and scanning parameters. We produced a set of high-probability regions of interest (ROIs) from these probabilistic maps; these and the probabilistic maps are made publicly available, together with a tool for querying the network membership probabilities associated with any given cortical location. These quantitative estimates and public tools may allow researchers to apply information about inter-subject consensus to their own fMRI studies, improving inferences about systems and their functional specializations.

Network variants are similar between task and rest states.

Recent work has demonstrated that individual-specific variations in functional networks (termed "network variants") can be identified in individuals using resting state functional magnetic resonance imaging (fMRI). These network variants exhibit reliability over time, suggesting that they may be trait-like markers of individual differences in brain organization. However, while networks variants are reliable at rest, is is still untested whether they are stable between task and rest states. Here, we use precision data from the Midnight Scan Club (MSC) to demonstrate that (1) task data can be used to identify network variants reliably, (2) these network variants show substantial spatial overlap with those observed in rest, although state-specific effects are present, (3) network variants assign to similar canonical functional networks in task and rest states, and (4) single tasks or a combination of multiple tasks produce similar network variants to rest. Together, these findings further reinforce the trait-like nature of network variants and demonstrate the utility of using task data to define network variants.

A Critical, Event-Related Appraisal of Denoising in Resting-State fMRI Studies.

This article advances two parallel lines of argument about resting-state functional magnetic resonance imaging (fMRI) signals, one empirical and one conceptual. The empirical line creates a four-part organization of the text: (1) head motion and respiration commonly cause distinct, major, unwanted influences (artifacts) in fMRI signals; (2) head motion and respiratory changes are, confoundingly, both related to psychological and clinical and biological variables of interest; (3) many fMRI denoising strategies fail to identify and remove one or the other kind of artifact; and (4) unremoved artifact, due to correlations of artifacts with variables of interest, renders studies susceptible to identifying variance of noninterest as variance of interest. Arising from these empirical observations is a conceptual argument: that an event-related approach to task-free scans, targeting common behaviors during scanning, enables fundamental distinctions among the kinds of signals present in the data, information which is vital to understanding the effects of denoising procedures. This event-related perspective permits statements like "Event X is associated with signals A, B, and C, each with particular spatial, temporal, and signal decay properties". Denoising approaches can then be tailored, via performance in known events, to permit or suppress certain kinds of signals based on their desirability.

Evaluating the Sensitivity of Resting-State BOLD Variability to Age and Cognition after Controlling for Motion and Cardiovascular Influences: A Network-Based Approach.

Recent functional magnetic resonance imaging (fMRI) studies report that moment-to-moment variability in the BOLD signal is related to differences in age and cognition and, thus, may be sensitive to age-dependent decline. However, head motion and/or cardiovascular health (CVH) may contaminate these relationships. We evaluated relationships between resting-state BOLD variability, age, and cognition, after characterizing and controlling for motion-related and cardiovascular influences, including pulse, blood pressure, BMI, and white matter hyperintensities (WMH), in a large (N = 422) resting-state fMRI sample of cognitively normal individuals (age 43-89). We found that resting-state BOLD variability was negatively related to age and positively related to cognition after maximally controlling for head motion. Age relationships also survived correction for CVH, but were greatly reduced when correcting for WMH alone. Our results suggest that network-based machine learning analyses of resting-state BOLD variability might yield reliable, sensitive measures to characterize age-related decline across a broad range of networks. Age-related differences in resting-state BOLD variability may be largely sensitive to processes related to WMH burden.

A set of functionally-defined brain regions with improved representation of the subcortex and cerebellum.

An important aspect of network-based analysis is robust node definition. This issue is critical for functional brain network analyses, as poor node choice can lead to spurious findings and misleading inferences about functional brain organization. Two sets of functional brain nodes from our group are well represented in the literature: (1) 264 volumetric regions of interest (ROIs) reported in Power et al., 2011, and (2) 333 cortical surface parcels reported in Gordon et al., 2016. However, subcortical and cerebellar structures are either incompletely captured or missing from these ROI sets. Therefore, properties of functional network organization involving the subcortex and cerebellum may be underappreciated thus far. Here, we apply a winner-take-all partitioning method to resting-state fMRI data to generate novel functionally-constrained ROIs in the thalamus, basal ganglia, amygdala, hippocampus, and cerebellum. We validate these ROIs in three datasets using several criteria, including agreement with existing literature and anatomical atlases. Further, we demonstrate that combining these ROIs with established cortical ROIs recapitulates and extends previously described functional network organization. This new set of ROIs is made publicly available for general use, including a full list of MNI coordinates and functional network labels.

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.