TemplateFlow: FAIR-sharing of multi-scale, multi-species brain models.

Reference anatomies of the brain ('templates') and corresponding atlases are the foundation for reporting standardized neuroimaging results. Currently, there is no registry of templates and atlases; therefore, the redistribution of these resources occurs either bundled within existing software or in ad hoc ways such as downloads from institutional sites and general-purpose data repositories. We introduce TemplateFlow as a publicly available framework for human and non-human brain models. The framework combines an open database with software for access, management, and vetting, allowing scientists to share their resources under FAIR-findable, accessible, interoperable, and reusable-principles. TemplateFlow enables multifaceted insights into brains across species, and supports multiverse analyses testing whether results generalize across standard references, scales, and in the long term, species.

Development of structure-function coupling in human brain networks during youth.

The protracted development of structural and functional brain connectivity within distributed association networks coincides with improvements in higher-order cognitive processes such as executive function. However, it remains unclear how white-matter architecture develops during youth to directly support coordinated neural activity. Here, we characterize the development of structure-function coupling using diffusion-weighted imaging and n-back functional MRI data in a sample of 727 individuals (ages 8 to 23 y). We found that spatial variability in structure-function coupling aligned with cortical hierarchies of functional specialization and evolutionary expansion. Furthermore, hierarchy-dependent age effects on structure-function coupling localized to transmodal cortex in both cross-sectional data and a subset of participants with longitudinal data (n = 294). Moreover, structure-function coupling in rostrolateral prefrontal cortex was associated with executive performance and partially mediated age-related improvements in executive function. Together, these findings delineate a critical dimension of adolescent brain development, whereby the coupling between structural and functional connectivity remodels to support functional specialization and cognition.

Unifying the Notions of Modularity and Core-Periphery Structure in Functional Brain Networks during Youth.

At rest, human brain functional networks display striking modular architecture in which coherent clusters of brain regions are activated. The modular account of brain function is pervasive, reliable, and reproducible. Yet, a complementary perspective posits a core-periphery or rich-club account of brain function, where hubs are densely interconnected with one another, allowing for integrative processing. Unifying these two perspectives has remained difficult due to the fact that the methodological tools to identify modules are entirely distinct from the methodological tools to identify core-periphery structure. Here, we leverage a recently-developed model-based approach-the weighted stochastic block model-that simultaneously uncovers modular and core-periphery structure, and we apply it to functional magnetic resonance imaging data acquired at rest in 872 youth of the Philadelphia Neurodevelopmental Cohort. We demonstrate that functional brain networks display rich mesoscale organization beyond that sought by modularity maximization techniques. Moreover, we show that this mesoscale organization changes appreciably over the course of neurodevelopment, and that individual differences in this organization predict individual differences in cognition more accurately than module organization alone. Broadly, our study provides a unified assessment of modular and core-periphery structure in functional brain networks, offering novel insights into their development and implications for behavior.

Associations between Neighborhood SES and Functional Brain Network Development.

Higher socioeconomic status (SES) in childhood is associated with stronger cognitive abilities, higher academic achievement, and lower incidence of mental illness later in development. While prior work has mapped the associations between neighborhood SES and brain structure, little is known about the relationship between SES and intrinsic neural dynamics. Here, we capitalize upon a large cross-sectional community-based sample (Philadelphia Neurodevelopmental Cohort, ages 8-22 years, n = 1012) to examine associations between age, SES, and functional brain network topology. We characterize this topology using a local measure of network segregation known as the clustering coefficient and find that it accounts for a greater degree of SES-associated variance than mesoscale segregation captured by modularity. High-SES youth displayed stronger positive associations between age and clustering than low-SES youth, and this effect was most pronounced for regions in the limbic, somatomotor, and ventral attention systems. The moderating effect of SES on positive associations between age and clustering was strongest for connections of intermediate length and was consistent with a stronger negative relationship between age and local connectivity in these regions in low-SES youth. Our findings suggest that, in late childhood and adolescence, neighborhood SES is associated with variation in the development of functional network structure in the human brain.

Diminished Cortical Thickness Is Associated with Impulsive Choice in Adolescence.

Adolescence is characterized by both maturation of brain structure and increased risk of negative outcomes from behaviors associated with impulsive decision-making. One important index of impulsive choice is delay discounting (DD), which measures the tendency to prefer smaller rewards available soon over larger rewards delivered after a delay. However, it remains largely unknown how individual differences in structural brain development may be associated with impulsive choice during adolescence. Leveraging a unique large sample of 427 human youths (208 males and 219 females) imaged as part of the Philadelphia Neurodevelopmental Cohort, we examined associations between delay discounting and cortical thickness within structural covariance networks. These structural networks were derived using non-negative matrix factorization, an advanced multivariate technique for dimensionality reduction, and analyzed using generalized additive models with penalized splines to capture both linear and nonlinear developmental effects. We found that impulsive choice, as measured by greater discounting, was most strongly associated with diminished cortical thickness in structural brain networks that encompassed the ventromedial prefrontal cortex, orbitofrontal cortex, temporal pole, and temporoparietal junction. Furthermore, structural brain networks predicted DD above and beyond cognitive performance. Together, these results suggest that reduced cortical thickness in regions known to be involved in value-based decision-making is a marker of impulsive choice during the critical period of adolescence.SIGNIFICANCE STATEMENT Risky behaviors during adolescence, such as initiation of substance use or reckless driving, are a major source of morbidity and mortality. In this study, we present evidence from a large sample of youths that diminished cortical thickness in specific structural brain networks is associated with impulsive choice. Notably, the strongest association between impulsive choice and brain structure was seen in regions implicated in value-based decision-making; namely, the ventromedial prefrontal and orbitofrontal cortices. Moving forward, such neuroanatomical markers of impulsivity may aid in the development of personalized interventions targeted to reduce risk of negative outcomes resulting from impulsivity during adolescence.

System-level matching of structural and functional connectomes in the human brain.

The brain can be considered as an information processing network, where complex behavior manifests as a result of communication between large-scale functional systems such as visual and default mode networks. As the communication between brain regions occurs through underlying anatomical pathways, it is important to define a "traffic pattern" that properly describes how the regions exchange information. Empirically, the choice of the traffic pattern can be made based on how well the functional connectivity between regions matches the structural pathways equipped with that traffic pattern. In this paper, we present a multimodal connectomics paradigm utilizing graph matching to measure similarity between structural and functional connectomes (derived from dMRI and fMRI data) at node, system, and connectome level. Through an investigation of the brain's structure-function relationship over a large cohort of 641 healthy developmental participants aged 8-22 years, we demonstrate that communicability as the traffic pattern describes the functional connectivity of the brain best, with large-scale systems having significant agreement between their structural and functional connectivity patterns. Notably, matching between structural and functional connectivity for the functionally specialized modular systems such as visual and motor networks are higher as compared to other more integrated systems. Additionally, we show that the negative functional connectivity between the default mode network (DMN) and motor, frontoparietal, attention, and visual networks is significantly associated with its underlying structural connectivity, highlighting the counterbalance between functional activation patterns of DMN and other systems. Finally, we investigated sex difference and developmental changes in brain and observed that similarity between structure and function changes with development.

Motion artifact in studies of functional connectivity: Characteristics and mitigation strategies.

Motion artifacts are now recognized as a major methodological challenge for studies of functional connectivity. As in-scanner motion is frequently correlated with variables of interest such as age, clinical status, cognitive ability, and symptom severity, in-scanner motion has the potential to introduce systematic bias. In this article, we describe how motion-related artifacts influence measures of functional connectivity and discuss the relative strengths and weaknesses of commonly used denoising strategies. Furthermore, we illustrate how motion can bias inference, using a study of brain development as an example. Finally, we highlight directions of ongoing and future research, and provide recommendations for investigators in the field. Hum Brain Mapp, 40:2033-2051, 2019. © 2017 Wiley Periodicals, Inc.

Faster family-wise error control for neuroimaging with a parametric bootstrap.

In neuroimaging, hundreds to hundreds of thousands of tests are performed across a set of brain regions or all locations in an image. Recent studies have shown that the most common family-wise error (FWE) controlling procedures in imaging, which rely on classical mathematical inequalities or Gaussian random field theory, yield FWE rates (FWER) that are far from the nominal level. Depending on the approach used, the FWER can be exceedingly small or grossly inflated. Given the widespread use of neuroimaging as a tool for understanding neurological and psychiatric disorders, it is imperative that reliable multiple testing procedures are available. To our knowledge, only permutation joint testing procedures have been shown to reliably control the FWER at the nominal level. However, these procedures are computationally intensive due to the increasingly available large sample sizes and dimensionality of the images, and analyses can take days to complete. Here, we develop a parametric bootstrap joint testing procedure. The parametric bootstrap procedure works directly with the test statistics, which leads to much faster estimation of adjusted p-values than resampling-based procedures while reliably controlling the FWER in sample sizes available in many neuroimaging studies. We demonstrate that the procedure controls the FWER in finite samples using simulations, and present region- and voxel-wise analyses to test for sex differences in developmental trajectories of cerebral blood flow.

Age-Related Effects and Sex Differences in Gray Matter Density, Volume, Mass, and Cortical Thickness from Childhood to Young Adulthood.

Developmental structural neuroimaging studies in humans have long described decreases in gray matter volume (GMV) and cortical thickness (CT) during adolescence. Gray matter density (GMD), a measure often assumed to be highly related to volume, has not been systematically investigated in development. We used T1 imaging data collected on the Philadelphia Neurodevelopmental Cohort to study age-related effects and sex differences in four regional gray matter measures in 1189 youths ranging in age from 8 to 23 years. Custom T1 segmentation and a novel high-resolution gray matter parcellation were used to extract GMD, GMV, gray matter mass (GMM; defined as GMD × GMV), and CT from 1625 brain regions. Nonlinear models revealed that each modality exhibits unique age-related effects and sex differences. While GMV and CT generally decrease with age, GMD increases and shows the strongest age-related effects, while GMM shows a slight decline overall. Females have lower GMV but higher GMD than males throughout the brain. Our findings suggest that GMD is a prime phenotype for the assessment of brain development and likely cognition and that periadolescent gray matter loss may be less pronounced than previously thought. This work highlights the need for combined quantitative histological MRI studies.SIGNIFICANCE STATEMENT This study demonstrates that different MRI-derived gray matter measures show distinct age and sex effects and should not be considered equivalent but complementary. It is shown for the first time that gray matter density increases from childhood to young adulthood, in contrast with gray matter volume and cortical thickness, and that females, who are known to have lower gray matter volume than males, have higher density throughout the brain. A custom preprocessing pipeline and a novel high-resolution parcellation were created to analyze brain scans of 1189 youths collected as part of the Philadelphia Neurodevelopmental Cohort. A clear understanding of normal structural brain development is essential for the examination of brain-behavior relationships, the study of brain disease, and, ultimately, clinical applications of neuroimaging.

Brain state expression and transitions are related to complex executive cognition in normative neurodevelopment.

Adolescence is marked by rapid development of executive function. Mounting evidence suggests that executive function in adults may be driven by dynamic control of neurophysiological processes. Yet, how these dynamics evolve over adolescence and contribute to cognitive development is unknown. In a sample of 780 youth aged 8-22 yr (42.7% male) from the Philadelphia Neurodevelopment Cohort, we use a dynamic graph approach to extract activation states in BOLD fMRI data from 264 brain regions. We construct a graph in which each observation in time is a node and the similarity in brain states at two different times is an edge. Using this graphical approach, we identify two primary brain states reminiscent of intrinsic and task-evoked systems. We show that time spent in these two states is higher in older adolescents, as is the flexibility with which the brain switches between them. Increasing time spent in primary states and flexibility among states relates to increases in a complex executive accuracy factor score over adolescence. Flexibility is more positively associated with accuracy toward early adulthood. These findings suggest that brain state dynamics are associated with complex executive function across a critical period of adolescence.

The impact of in-scanner head motion on structural connectivity derived from diffusion MRI.

Multiple studies have shown that data quality is a critical confound in the construction of brain networks derived from functional MRI. This problem is particularly relevant for studies of human brain development where important variables (such as participant age) are correlated with data quality. Nevertheless, the impact of head motion on estimates of structural connectivity derived from diffusion tractography methods remains poorly characterized. Here, we evaluated the impact of in-scanner head motion on structural connectivity using a sample of 949 participants (ages 8-23 years old) who passed a rigorous quality assessment protocol for diffusion magnetic resonance imaging (dMRI) acquired as part of the Philadelphia Neurodevelopmental Cohort. Structural brain networks were constructed for each participant using both deterministic and probabilistic tractography. We hypothesized that subtle variation in head motion would systematically bias estimates of structural connectivity and confound developmental inference, as observed in previous studies of functional connectivity. Even following quality assurance and retrospective correction for head motion, eddy currents, and field distortions, in-scanner head motion significantly impacted the strength of structural connectivity in a consistency- and length-dependent manner. Specifically, increased head motion was associated with reduced estimates of structural connectivity for network edges with high inter-subject consistency, which included both short- and long-range connections. In contrast, motion inflated estimates of structural connectivity for low-consistency network edges that were primarily shorter-range. Finally, we demonstrate that age-related differences in head motion can both inflate and obscure developmental inferences on structural connectivity. Taken together, these data delineate the systematic impact of head motion on structural connectivity, and provide a critical context for identifying motion-related confounds in studies of structural brain network development.

Quantitative assessment of structural image quality.

Data quality is increasingly recognized as one of the most important confounding factors in brain imaging research. It is particularly important for studies of brain development, where age is systematically related to in-scanner motion and data quality. Prior work has demonstrated that in-scanner head motion biases estimates of structural neuroimaging measures. However, objective measures of data quality are not available for most structural brain images. Here we sought to identify quantitative measures of data quality for T1-weighted volumes, describe how these measures relate to cortical thickness, and delineate how this in turn may bias inference regarding associations with age in youth. Three highly-trained raters provided manual ratings of 1840 raw T1-weighted volumes. These images included a training set of 1065 images from Philadelphia Neurodevelopmental Cohort (PNC), a test set of 533 images from the PNC, as well as an external test set of 242 adults acquired on a different scanner. Manual ratings were compared to automated quality measures provided by the Preprocessed Connectomes Project's Quality Assurance Protocol (QAP), as well as FreeSurfer's Euler number, which summarizes the topological complexity of the reconstructed cortical surface. Results revealed that the Euler number was consistently correlated with manual ratings across samples. Furthermore, the Euler number could be used to identify images scored "unusable" by human raters with a high degree of accuracy (AUC: 0.98-0.99), and out-performed proxy measures from functional timeseries acquired in the same scanning session. The Euler number also was significantly related to cortical thickness in a regionally heterogeneous pattern that was consistent across datasets and replicated prior results. Finally, data quality both inflated and obscured associations with age during adolescence. Taken together, these results indicate that reliable measures of data quality can be automatically derived from T1-weighted volumes, and that failing to control for data quality can systematically bias the results of studies of brain maturation.

Benchmarking of participant-level confound regression strategies for the control of motion artifact in studies of functional connectivity.

Since initial reports regarding the impact of motion artifact on measures of functional connectivity, there has been a proliferation of participant-level confound regression methods to limit its impact. However, many of the most commonly used techniques have not been systematically evaluated using a broad range of outcome measures. Here, we provide a systematic evaluation of 14 participant-level confound regression methods in 393 youths. Specifically, we compare methods according to four benchmarks, including the residual relationship between motion and connectivity, distance-dependent effects of motion on connectivity, network identifiability, and additional degrees of freedom lost in confound regression. Our results delineate two clear trade-offs among methods. First, methods that include global signal regression minimize the relationship between connectivity and motion, but result in distance-dependent artifact. In contrast, censoring methods mitigate both motion artifact and distance-dependence, but use additional degrees of freedom. Importantly, less effective de-noising methods are also unable to identify modular network structure in the connectome. Taken together, these results emphasize the heterogeneous efficacy of existing methods, and suggest that different confound regression strategies may be appropriate in the context of specific scientific goals.

hyve, a compositional visualisation engine for brain imaging data.

Visualisations facilitate the interpretation of geometrically structured data and results. However, heterogeneous geometries-such as volumes, surfaces, and networks-have traditionally mandated different software approaches. We introduce hyve, a Python library that uses a compositional functional framework to enable parametric implementation of custom visualisations for different brain geometries. Under this framework, users compose a reusable visualisation protocol from geometric primitives for representing data geometries, input primitives for common data formats and research objectives, and output primitives for producing interactive displays or configurable snapshots. hyve also writes documentation for user-constructed protocols, automates serial production of multiple visualisations, and includes an API for semantically organising an editable multi-panel figure. Through the seamless composition of input, output, and geometric primitives, hyve supports creating visualisations for a range of neuroimaging research objectives.

Efficient coding in the economics of human brain connectomics.

In systems neuroscience, most models posit that brain regions communicate information under constraints of efficiency. Yet, evidence for efficient communication in structural brain networks characterized by hierarchical organization and highly connected hubs remains sparse. The principle of efficient coding proposes that the brain transmits maximal information in a metabolically economical or compressed form to improve future behavior. To determine how structural connectivity supports efficient coding, we develop a theory specifying minimum rates of message transmission between brain regions to achieve an expected fidelity, and we test five predictions from the theory based on random walk communication dynamics. In doing so, we introduce the metric of compression efficiency, which quantifies the trade-off between lossy compression and transmission fidelity in structural networks. In a large sample of youth (n = 1,042; age 8-23 years), we analyze structural networks derived from diffusion-weighted imaging and metabolic expenditure operationalized using cerebral blood flow. We show that structural networks strike compression efficiency trade-offs consistent with theoretical predictions. We find that compression efficiency prioritizes fidelity with development, heightens when metabolic resources and myelination guide communication, explains advantages of hierarchical organization, links higher input fidelity to disproportionate areal expansion, and shows that hubs integrate information by lossy compression. Lastly, compression efficiency is predictive of behavior-beyond the conventional network efficiency metric-for cognitive domains including executive function, memory, complex reasoning, and social cognition. Our findings elucidate how macroscale connectivity supports efficient coding and serve to foreground communication processes that utilize random walk dynamics constrained by network connectivity.

Temporal sequences of brain activity at rest are constrained by white matter structure and modulated by cognitive demands.

A diverse set of white matter connections supports seamless transitions between cognitive states. However, it remains unclear how these connections guide the temporal progression of large-scale brain activity patterns in different cognitive states. Here, we analyze the brain's trajectories across a set of single time point activity patterns from functional magnetic resonance imaging data acquired during the resting state and an n-back working memory task. We find that specific temporal sequences of brain activity are modulated by cognitive load, associated with age, and related to task performance. Using diffusion-weighted imaging acquired from the same subjects, we apply tools from network control theory to show that linear spread of activity along white matter connections constrains the probabilities of these sequences at rest, while stimulus-driven visual inputs explain the sequences observed during the n-back task. Overall, these results elucidate the structural underpinnings of cognitively and developmentally relevant spatiotemporal brain dynamics.

Optimization of energy state transition trajectory supports the development of executive function during youth.

Executive function develops during adolescence, yet it remains unknown how structural brain networks mature to facilitate activation of the fronto-parietal system, which is critical for executive function. In a sample of 946 human youths (ages 8-23y) who completed diffusion imaging, we capitalized upon recent advances in linear dynamical network control theory to calculate the energetic cost necessary to activate the fronto-parietal system through the control of multiple brain regions given existing structural network topology. We found that the energy required to activate the fronto-parietal system declined with development, and the pattern of regional energetic cost predicts unseen individuals' brain maturity. Finally, energetic requirements of the cingulate cortex were negatively correlated with executive performance, and partially mediated the development of executive performance with age. Our results reveal a mechanism by which structural networks develop during adolescence to reduce the theoretical energetic costs of transitions to activation states necessary for executive function.

Adolescents are known for taking risks, from driving too fast to experimenting with drugs and alcohol. Such behaviors tend to decrease as individuals move into adulthood. Most people in their mid-twenties have greater self-control than they did as teenagers. They are also often better at planning, sustaining attention, and inhibiting impulsive behaviors. These skills, which are known as executive functions, develop over the course of adolescence. Executive functions rely upon a series of brain regions distributed across the frontal lobe and the lobe that sits just behind it, the parietal lobe. Fiber tracts connect these regions to form a fronto-parietal network. These fiber tracts are also referred to as white matter due to the whitish fatty material that surrounds and insulates them. Cui et al. now show that changes in white matter networks have implications for teen behavior. Almost 950 healthy young people aged between 8 and 23 years underwent a type of brain scan called diffusion-weighted imaging that visualizes white matter. The scans revealed that white matter networks in the frontal and parietal lobes mature over adolescence. This makes it easier for individuals to activate their fronto-parietal networks by decreasing the amount of energy required. Cui et al. show that a computer model can predict the maturity of a person's brain based on the energy needed to activate their fronto-parietal networks. These changes help explain why executive functions improve during adolescence. This in turn explains why behaviors such as risk-taking tend to decrease with age. That said, adults with various psychiatric disorders, such as ADHD and psychosis, often show impaired executive functions. In the future, it may be possible to reduce these impairments by applying magnetic fields to the scalp to reduce the activity of specific brain regions. The techniques used in the current study could help reveal which brain regions to target with this approach.

Analysis of task-based functional MRI data preprocessed with fMRIPrep.

Functional magnetic resonance imaging (fMRI) is a standard tool to investigate the neural correlates of cognition. fMRI noninvasively measures brain activity, allowing identification of patterns evoked by tasks performed during scanning. Despite the long history of this technique, the idiosyncrasies of each dataset have led to the use of ad-hoc preprocessing protocols customized for nearly every different study. This approach is time consuming, error prone and unsuitable for combining datasets from many sources. Here we showcase fMRIPrep (http://fmriprep.org), a robust tool to prepare human fMRI data for statistical analysis. This software instrument addresses the reproducibility concerns of the established protocols for fMRI preprocessing. By leveraging the Brain Imaging Data Structure to standardize both the input datasets (MRI data as stored by the scanner) and the outputs (data ready for modeling and analysis), fMRIPrep is capable of preprocessing a diversity of datasets without manual intervention. In support of the growing popularity of fMRIPrep, this protocol describes how to integrate the tool in a task-based fMRI investigation workflow.

Evaluation of confound regression strategies for the mitigation of micromovement artifact in studies of dynamic resting-state functional connectivity and multilayer network modularity.

Dynamic functional connectivity reflects the spatiotemporal organization of spontaneous brain activity in health and disease. Dynamic functional connectivity may be susceptible to artifacts induced by participant motion. This report provides a systematic evaluation of 12 commonly used participant-level confound regression strategies designed to mitigate the effects of micromovements in a sample of 393 youths (ages 8-22 years). Each strategy was evaluated according to a number of benchmarks, including (a) the residual association between participant motion and edge dispersion, (b) distance-dependent effects of motion on edge dispersion, (c) the degree to which functional subnetworks could be identified by multilayer modularity maximization, and (d) measures of module reconfiguration, including node flexibility and node promiscuity. Results indicate variability in the effectiveness of the evaluated pipelines across benchmarks. Methods that included global signal regression were the most consistently effective de-noising strategies.