Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy.

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

Diffusion MRI head motion correction methods are highly accurate but impacted by denoising and sampling scheme.

Head motion correction is particularly challenging in diffusion-weighted MRI (dMRI) scans due to the dramatic changes in image contrast at different gradient strengths and directions. Head motion correction is typically performed using a Gaussian Process model implemented in FSL's Eddy. Recently, the 3dSHORE-based SHORELine method was introduced that does not require shell-based acquisitions, but it has not been previously benchmarked. Here we perform a comprehensive evaluation of both methods on realistic simulations of a software fiber phantom that provides known ground-truth head motion. We demonstrate that both methods perform remarkably well, but that performance can be impacted by sampling scheme and the extent of head motion and the denoising strategy applied before head motion correction. Furthermore, we find Eddy benefits from denoising the data first with MP-PCA. In sum, we provide the most extensive known benchmarking of dMRI head motion correction, together with extensive simulation data and a reproducible workflow. PRACTITIONER POINTS: Both Eddy and SHORELine head motion correction methods performed quite well on a large variety of simulated data. Denoising with MP-PCA can improve head motion correction performance when Eddy is used. SHORELine effectively corrects motion in non-shelled diffusion spectrum imaging data.

Development of white matter fiber covariance networks supports executive function in youth.

During adolescence, the brain undergoes extensive changes in white matter structure that support cognition. Data-driven approaches applied to cortical surface properties have led the field to understand brain development as a spatially and temporally coordinated mechanism that follows hierarchically organized gradients of change. Although white matter development also appears asynchronous, previous studies have relied largely on anatomical tract-based atlases, precluding a direct assessment of how white matter structure is spatially and temporally coordinated. Harnessing advances in diffusion modeling and machine learning, we identified 14 data-driven patterns of covarying white matter structure in a large sample of youth. Fiber covariance networks aligned with known major tracts, while also capturing distinct patterns of spatial covariance across distributed white matter locations. Most networks showed age-related increases in fiber network properties, which were also related to developmental changes in executive function. This study delineates data-driven patterns of white matter development that support cognition.

ModelArray: An R package for statistical analysis of fixel-wise data.

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Development of White Matter Fiber Covariance Networks Supports Executive Function in Youth.

The white matter architecture of the human brain undergoes substantial development throughout childhood and adolescence, allowing for more efficient signaling between brain regions that support executive function. Increasingly, the field understands grey matter development as a spatially and temporally coordinated mechanism that follows hierarchically organized gradients of change. While white matter development also appears asynchronous, previous studies have largely relied on anatomical atlases to characterize white matter tracts, precluding a direct assessment of how white matter structure is spatially and temporally coordinated. Here, we leveraged advances in diffusion modeling and unsupervised machine learning to delineate white matter fiber covariance networks comprised of structurally similar areas of white matter in a cross-sectional sample of 939 youth aged 8-22 years. We then evaluated associations between fiber covariance network structural properties with both age and executive function using generalized additive models. The identified fiber covariance networks aligned with the known architecture of white matter while simultaneously capturing novel spatial patterns of coordinated maturation. Fiber covariance networks showed heterochronous increases in fiber density and cross section that generally followed hierarchically organized temporal patterns of cortical development, with the greatest increases in unimodal sensorimotor networks and the most prolonged increases in superior and anterior transmodal networks. Notably, we found that executive function was associated with structural features of limbic and association networks. Taken together, this study delineates data-driven patterns of white matter network development that support cognition and align with major axes of brain maturation.

Development of top-down cortical propagations in youth.

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Spatially-enhanced clusterwise inference for testing and localizing intermodal correspondence.

With the increasing availability of neuroimaging data from multiple modalities-each providing a different lens through which to study brain structure or function-new techniques for comparing, integrating, and interpreting information within and across modalities have emerged. Recent developments include hypothesis tests of associations between neuroimaging modalities, which can be used to determine the statistical significance of intermodal associations either throughout the entire brain or within anatomical subregions or functional networks. While these methods provide a crucial foundation for inference on intermodal relationships, they cannot be used to answer questions about where in the brain these associations are most pronounced. In this paper, we introduce a new method, called CLEAN-R, that can be used both to test intermodal correspondence throughout the brain and also to localize this correspondence. Our method involves first adjusting for the underlying spatial autocorrelation structure within each modality before aggregating information within small clusters to construct a map of enhanced test statistics. Using structural and functional magnetic resonance imaging data from a subsample of children and adolescents from the Philadelphia Neurodevelopmental Cohort, we conduct simulations and data analyses where we illustrate the high statistical power and nominal type I error levels of our method. By constructing an interpretable map of group-level correspondence using spatially-enhanced test statistics, our method offers insights beyond those provided by earlier methods.

Curation of BIDS (CuBIDS): A workflow and software package for streamlining reproducible curation of large BIDS datasets.

The Brain Imaging Data Structure (BIDS) is a specification accompanied by a software ecosystem that was designed to create reproducible and automated workflows for processing neuroimaging data. BIDS Apps flexibly build workflows based on the metadata detected in a dataset. However, even BIDS valid metadata can include incorrect values or omissions that result in inconsistent processing across sessions. Additionally, in large-scale, heterogeneous neuroimaging datasets, hidden variability in metadata is difficult to detect and classify. To address these challenges, we created a Python-based software package titled "Curation of BIDS" (CuBIDS), which provides an intuitive workflow that helps users validate and manage the curation of their neuroimaging datasets. CuBIDS includes a robust implementation of BIDS validation that scales to large samples and incorporates DataLad--a version control software package for data--as an optional dependency to ensure reproducibility and provenance tracking throughout the entire curation process. CuBIDS provides tools to help users perform quality control on their images' metadata and identify unique combinations of imaging parameters. Users can then execute BIDS Apps on a subset of participants that represent the full range of acquisition parameters that are present, accelerating pipeline testing on large datasets.

ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion.

Arterial spin labeled (ASL) magnetic resonance imaging (MRI) is the primary method for noninvasively measuring regional brain perfusion in humans. We introduce ASLPrep, a suite of software pipelines that ensure the reproducible and generalizable processing of ASL MRI data.

Mobile footprinting: linking individual distinctiveness in mobility patterns to mood, sleep, and brain functional connectivity.

Mapping individual differences in behavior is fundamental to personalized neuroscience, but quantifying complex behavior in real world settings remains a challenge. While mobility patterns captured by smartphones have increasingly been linked to a range of psychiatric symptoms, existing research has not specifically examined whether individuals have person-specific mobility patterns. We collected over 3000 days of mobility data from a sample of 41 adolescents and young adults (age 17-30 years, 28 female) with affective instability. We extracted summary mobility metrics from GPS and accelerometer data and used their covariance structures to identify individuals and calculated the individual identification accuracy-i.e., their "footprint distinctiveness". We found that statistical patterns of smartphone-based mobility features represented unique "footprints" that allow individual identification (p < 0.001). Critically, mobility footprints exhibited varying levels of person-specific distinctiveness (4-99%), which was associated with age and sex. Furthermore, reduced individual footprint distinctiveness was associated with instability in affect (p < 0.05) and circadian patterns (p < 0.05) as measured by environmental momentary assessment. Finally, brain functional connectivity, especially those in the somatomotor network, was linked to individual differences in mobility patterns (p < 0.05). Together, these results suggest that real-world mobility patterns may provide individual-specific signatures relevant for studies of development, sleep, and psychopathology.

Developmental coupling of cerebral blood flow and fMRI fluctuations in youth.

The functions of the human brain are metabolically expensive and reliant on coupling between cerebral blood flow (CBF) and neural activity, yet how this coupling evolves over development remains unexplored. Here, we examine the relationship between CBF, measured by arterial spin labeling, and the amplitude of low-frequency fluctuations (ALFF) from resting-state magnetic resonance imaging across a sample of 831 children (478 females, aged 8-22 years) from the Philadelphia Neurodevelopmental Cohort. We first use locally weighted regressions on the cortical surface to quantify CBF-ALFF coupling. We relate coupling to age, sex, and executive functioning with generalized additive models and assess network enrichment via spin testing. We demonstrate regionally specific changes in coupling over age and show that variations in coupling are related to biological sex and executive function. Our results highlight the importance of CBF-ALFF coupling throughout development; we discuss its potential as a future target for the study of neuropsychiatric diseases.

A simple permutation-based test of intermodal correspondence.

Many key findings in neuroimaging studies involve similarities between brain maps, but statistical methods used to measure these findings have varied. Current state-of-the-art methods involve comparing observed group-level brain maps (after averaging intensities at each image location across multiple subjects) against spatial null models of these group-level maps. However, these methods typically make strong and potentially unrealistic statistical assumptions, such as covariance stationarity. To address these issues, in this article we propose using subject-level data and a classical permutation testing framework to test and assess similarities between brain maps. Our method is comparable to traditional permutation tests in that it involves randomly permuting subjects to generate a null distribution of intermodal correspondence statistics, which we compare to an observed statistic to estimate a p-value. We apply and compare our method in simulated and real neuroimaging data from the Philadelphia Neurodevelopmental Cohort. We show that our method performs well for detecting relationships between modalities known to be strongly related (cortical thickness and sulcal depth), and it is conservative when an association would not be expected (cortical thickness and activation on the n-back working memory task). Notably, our method is the most flexible and reliable for localizing intermodal relationships within subregions of the brain and allows for generalizable statistical inference.

FlywheelTools: Data Curation and Manipulation on the Flywheel Platform.

The recent and growing focus on reproducibility in neuroimaging studies has led many major academic centers to use cloud-based imaging databases for storing, analyzing, and sharing complex imaging data. Flywheel is one such database platform that offers easily accessible, large-scale data management, along with a framework for reproducible analyses through containerized pipelines. The Brain Imaging Data Structure (BIDS) is the de facto standard for neuroimaging data, but curating neuroimaging data into BIDS can be a challenging and time-consuming task. In particular, standard solutions for BIDS curation are limited on Flywheel. To address these challenges, we developed "FlywheelTools," a software toolbox for reproducible data curation and manipulation on Flywheel. FlywheelTools includes two elements: fw-heudiconv, for heuristic-driven curation of data into BIDS, and flaudit, which audits and inventories projects on Flywheel. Together, these tools accelerate reproducible neuroscience research on the widely used Flywheel platform.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Momentary changes in heart rate variability can detect risk for emotional eating episodes.

Emotion dysregulation is a known risk factor for a variety of maladaptive eating behaviors, including emotional eating (Crockett, Myhre, & Rokke, 2015; Lavender et al., 2015). New passive sensing technologies offer the prospect of detecting emotion dysregulation in real-time through measurement of heart rate variability (HRV), a transdiagnostic bio-signal of emotion regulation, which may in turn signal risk of engaging in a maladaptive eating behavior. In the current study, our primary aim was to test the hypothesis that momentary changes in HRV can be used to detect risk of experiencing an emotional eating episode in an ecologically valid setting using a wrist worn sensor with acceptable classification accuracy. Participants were 21 adults with clinically significant emotional eating behaviors. Participants wore the Empatica E4 wrist-sensor and tracked all emotional eating episodes using ecological momentary assessment for four weeks. Time and frequency domain features of HRV were extracted in the 30-min period preceding emotional eating episodes and control cases (defined as the 30 min prior to an EMA survey that did not contain an emotional eating episode). Support vector machine (SVM) learning models were implemented using time domain and frequency domain features. SVM models using frequency domain features achieved the highest classification accuracy (77.99%), sensitivity (78.75%), and specificity (75.00%), consistent with standards deemed acceptable for the prediction of event-level health behavior. SVM models using time domain features still performed above chance, though were less accurate at classifying episodes (accuracy 63.48%, sensitivity 62.68%, and specificity 70.00%) and did not meet acceptable classification accuracy. Wearable sensors that assess HRV show promise as a tool for capturing risk of engaging in emotional eating episodes.

Autism spectrum disorder (ASD) symptom profiles of children with comorbid Down syndrome (DS) and ASD: A comparison with children with DS-only and ASD-only.

BACKGROUND: Down syndrome (DS) is associated with increased rates of autism spectrum disorder (ASD), characterized by social-communicative impairments (SOC-COM) and repetitive behaviors and interests (RBI). However, little is known about the ASD symptom presentation in children with DS + ASD. AIMS: The current study sought to describe parent-report of SOC-COM and RBI symptoms on the Autism Diagnostic Interview -Revised (ADI-R) in children with DS (n = 22), DS + ASD (n = 11), and ASD (n = 66). METHOD: SOC-COM and RBI scores from the ADI-R were compared across the groups whose autism status was ascertained using the Autism Diagnostic Observation Schedule. RESULTS: Differences in SOC-COM and RBI symptom severity was observed. The general pattern of findings was ASD > DS+ASD > DS. Dissimilar ASD symptom profiles were observed across groups. In ASD, SOC-COM scores were higher than RBI scores; in DS + ASD, similar SOC-COM and RBI scores were observed. Lastly, SOC-COM impairments were highly related to verbal cognition in youth with DS + ASD but not in those with DS or ASD. CONCLUSIONS AND IMPLICATIONS: These findings suggest that children with DS + ASD have a distinct profile of ASD symptoms that differs from peers with either disorder in isolation. Thus, care should be taken in evaluating and designing treatments for this group.

Application of a new dietary pattern analysis method in nutritional epidemiology.

BACKGROUND: Diet plays an important role in chronic disease, and the use of dietary pattern analysis has grown rapidly as a way of deconstructing the complexity of nutritional intake and its relation to health. Pattern analysis methods, such as principal component analysis (PCA), have been used to investigate various dimensions of diet. Existing analytic methods, however, do not fully utilize the predictive potential of dietary assessment data. In particular, these methods are often suboptimal at predicting clinically important variables. METHODS: We propose a new dietary pattern analysis method using the advanced LASSO (Least Absolute Shrinkage and Selection Operator) model to improve the prediction of disease-related risk factors. Despite the potential advantages of LASSO, this is the first time that the model has been adapted for dietary pattern analysis. Hence, the systematic evaluation of the LASSO model as applied to dietary data and health outcomes is highly innovative and novel. Using Food Frequency Questionnaire data from NHANES 2005-2006, we apply PCA and LASSO to identify dietary patterns related to cardiovascular disease risk factors in healthy US adults (n = 2609) after controlling for confounding variables (e.g., age and BMI). Both analyses account for the sampling weights. Model performance in terms of prediction accuracy is evaluated using an independent test set. RESULTS: PCA yields 10 principal components (PCs) that together account for 65% of the variation in the data set and represent distinct dietary patterns. These PCs are then used as predictors in a regression model to predict cardiovascular disease risk factors. We find that LASSO better predicts levels of triglycerides, LDL cholesterol, HDL cholesterol, and total cholesterol (adjusted R2 = 0.861, 0.899, 0.890, and 0.935 respectively) than does the traditional, linear-regression-based, dietary pattern analysis method (adjusted R2 = 0.163, 0.005, 0.235, and 0.024 respectively) when the latter is applied to components derived from PCA. CONCLUSIONS: The proposed method is shown to be an appropriate and promising statistical means of deriving dietary patterns predictive of cardiovascular disease risk. Future studies, involving different diseases and risk factors, will be necessary before LASSO's broader usefulness in nutritional epidemiology can be established.