A sensorimotor-association axis of thalamocortical connection development.

Human cortical development follows a sensorimotor-to-association sequence during childhood and adolescence1-6. The brain's capacity to enact this sequence over decades indicates that it relies on intrinsic mechanisms to regulate inter-regional differences in the timing of cortical maturation, yet regulators of human developmental chronology are not well understood. Given evidence from animal models that thalamic axons modulate windows of cortical plasticity7-12, here we evaluate the overarching hypothesis that structural connections between the thalamus and cortex help to coordinate cortical maturational heterochronicity during youth. We first introduce, cortically annotate, and anatomically validate a new atlas of human thalamocortical connections using diffusion tractography. By applying this atlas to three independent youth datasets (ages 8-23 years; total N = 2,676), we reproducibly demonstrate that thalamocortical connections develop along a maturational gradient that aligns with the cortex's sensorimotor-association axis. Associative cortical regions with thalamic connections that take longest to mature exhibit protracted expression of neurochemical, structural, and functional markers indicative of higher circuit plasticity as well as heightened environmental sensitivity. This work highlights a central role for the thalamus in the orchestration of hierarchically organized and environmentally sensitive windows of cortical developmental malleability.

ANTsX neuroimaging-derived structural phenotypes of UK Biobank.

UK Biobank is a large-scale epidemiological resource for investigating prospective correlations between various lifestyle, environmental, and genetic factors with health and disease progression. In addition to individual subject information obtained through surveys and physical examinations, a comprehensive neuroimaging battery consisting of multiple modalities provides imaging-derived phenotypes (IDPs) that can serve as biomarkers in neuroscience research. In this study, we augment the existing set of UK Biobank neuroimaging structural IDPs, obtained from well-established software libraries such as FSL and FreeSurfer, with related measurements acquired through the Advanced Normalization Tools Ecosystem. This includes previously established cortical and subcortical measurements defined, in part, based on the Desikan-Killiany-Tourville atlas. Also included are morphological measurements from two recent developments: medial temporal lobe parcellation of hippocampal and extra-hippocampal regions in addition to cerebellum parcellation and thickness based on the Schmahmann anatomical labeling. Through predictive modeling, we assess the clinical utility of these IDP measurements, individually and in combination, using commonly studied phenotypic correlates including age, fluid intelligence, numeric memory, and several other sociodemographic variables. The predictive accuracy of these IDP-based models, in terms of root-mean-squared-error or area-under-the-curve for continuous and categorical variables, respectively, provides comparative insights between software libraries as well as potential clinical interpretability. Results demonstrate varied performance between package-based IDP sets and their combination, emphasizing the need for careful consideration in their selection and utilization.

Greater white matter degeneration and lower structural connectivity in non-amnestic vs. amnestic Alzheimer's disease.

Introduction: Multimodal evidence indicates Alzheimer's disease (AD) is characterized by early white matter (WM) changes that precede overt cognitive impairment. WM changes have overwhelmingly been investigated in typical, amnestic mild cognitive impairment and AD; fewer studies have addressed WM change in atypical, non-amnestic syndromes. We hypothesized each non-amnestic AD syndrome would exhibit WM differences from amnestic and other non-amnestic syndromes. Materials and methods: Participants included 45 cognitively normal (CN) individuals; 41 amnestic AD patients; and 67 patients with non-amnestic AD syndromes including logopenic-variant primary progressive aphasia (lvPPA, n = 32), posterior cortical atrophy (PCA, n = 17), behavioral variant AD (bvAD, n = 10), and corticobasal syndrome (CBS, n = 8). All had T1-weighted MRI and 30-direction diffusion-weighted imaging (DWI). We performed whole-brain deterministic tractography between 148 cortical and subcortical regions; connection strength was quantified by tractwise mean generalized fractional anisotropy. Regression models assessed effects of group and phenotype as well as associations with grey matter volume. Topological analyses assessed differences in persistent homology (numbers of graph components and cycles). Additionally, we tested associations of topological metrics with global cognition, disease duration, and DWI microstructural metrics. Results: Both amnestic and non-amnestic patients exhibited lower WM connection strength than CN participants in corpus callosum, cingulum, and inferior and superior longitudinal fasciculi. Overall, non-amnestic patients had more WM disease than amnestic patients. LvPPA patients had left-lateralized WM degeneration; PCA patients had reductions in connections to bilateral posterior parietal, occipital, and temporal areas. Topological analysis showed the non-amnestic but not the amnestic group had more connected components than controls, indicating persistently lower connectivity. Longer disease duration and cognitive impairment were associated with more connected components and fewer cycles in individuals' brain graphs. Discussion: We have previously reported syndromic differences in GM degeneration and tau accumulation between AD syndromes; here we find corresponding differences in WM tracts connecting syndrome-specific epicenters. Determining the reasons for selective WM degeneration in non-amnestic AD is a research priority that will require integration of knowledge from neuroimaging, biomarker, autopsy, and functional genetic studies. Furthermore, longitudinal studies to determine the chronology of WM vs. GM degeneration will be key to assessing evidence for WM-mediated tau spread.

A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging.

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

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.

Novel data-driven subtypes and stages of brain atrophy in the ALS-FTD spectrum.

BACKGROUND: TDP-43 proteinopathies represent a spectrum of neurological disorders, anchored clinically on either end by amyotrophic lateral sclerosis (ALS) and frontotemporal degeneration (FTD). The ALS-FTD spectrum exhibits a diverse range of clinical presentations with overlapping phenotypes, highlighting its heterogeneity. This study was aimed to use disease progression modeling to identify novel data-driven spatial and temporal subtypes of brain atrophy and its progression in the ALS-FTD spectrum. METHODS: We used a data-driven procedure to identify 13 anatomic clusters of brain volume for 57 behavioral variant FTD (bvFTD; with either autopsy-confirmed TDP-43 or TDP-43 proteinopathy-associated genetic variants), 103 ALS, and 47 ALS-FTD patients with likely TDP-43. A Subtype and Stage Inference (SuStaIn) model was trained to identify subtypes of individuals along the ALS-FTD spectrum with distinct brain atrophy patterns, and we related subtypes and stages to clinical, genetic, and neuropathological features of disease. RESULTS: SuStaIn identified three novel subtypes: two disease subtypes with predominant brain atrophy in either prefrontal/somatomotor regions or limbic-related regions, and a normal-appearing group without obvious brain atrophy. The limbic-predominant subtype tended to present with more impaired cognition, higher frequencies of pathogenic variants in TBK1 and TARDBP genes, and a higher proportion of TDP-43 types B, E and C. In contrast, the prefrontal/somatomotor-predominant subtype had higher frequencies of pathogenic variants in C9orf72 and GRN genes and higher proportion of TDP-43 type A. The normal-appearing brain group showed higher frequency of ALS relative to ALS-FTD and bvFTD patients, higher cognitive capacity, higher proportion of lower motor neuron onset, milder motor symptoms, and lower frequencies of genetic pathogenic variants. The overall SuStaIn stages also correlated with evidence for clinical progression including longer disease duration, higher King's stage, and cognitive decline. Additionally, SuStaIn stages differed across clinical phenotypes, genotypes and types of TDP-43 pathology. CONCLUSIONS: Our findings suggest distinct neurodegenerative subtypes of disease along the ALS-FTD spectrum that can be identified in vivo, each with distinct brain atrophy, clinical, genetic and pathological patterns.

ANTsX neuroimaging-derived structural phenotypes of UK Biobank.

UK Biobank is a large-scale epidemiological resource for investigating prospective correlations between various lifestyle, environmental, and genetic factors with health and disease progression. In addition to individual subject information obtained through surveys and physical examinations, a comprehensive neuroimaging battery consisting of multiple modalities provides imaging-derived phenotypes (IDPs) that can serve as biomarkers in neuroscience research. In this study, we augment the existing set of UK Biobank neuroimaging structural IDPs, obtained from well-established software libraries such as FSL and FreeSurfer, with related measurements acquired through the Advanced Normalization Tools Ecosystem. This includes previously established cortical and subcortical measurements defined, in part, based on the Desikan-Killiany-Tourville atlas. Also included are morphological measurements from two recent developments: medial temporal lobe parcellation of hippocampal and extra-hippocampal regions in addition to cerebellum parcellation and thickness based on the Schmahmann anatomical labeling. Through predictive modeling, we assess the clinical utility of these IDP measurements, individually and in combination, using commonly studied phenotypic correlates including age, fluid intelligence, numeric memory, and several other sociodemographic variables. The predictive accuracy of these IDP-based models, in terms of root-mean-squared-error or area-under-the-curve for continuous and categorical variables, respectively, provides comparative insights between software libraries as well as potential clinical interpretability. Results demonstrate varied performance between package-based IDP sets and their combination, emphasizing the need for careful consideration in their selection and utilization.

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.

Novel data-driven subtypes and stages of brain atrophy in the ALS-FTD spectrum.

Background: TDP-43 proteinopathies represents a spectrum of neurological disorders, anchored clinically on either end by amyotrophic lateral sclerosis (ALS) and frontotemporal degeneration (FTD). The ALS-FTD spectrum exhibits a diverse range of clinical presentations with overlapping phenotypes, highlighting its heterogeneity. This study aimed to use disease progression modeling to identify novel data-driven spatial and temporal subtypes of brain atrophy and its progression in the ALS-FTD spectrum. Methods: We used a data-driven procedure to identify 13 anatomic clusters of brain volumes for 57 behavioral variant FTD (bvFTD; with either autopsy-confirmed TDP-43 or TDP-43 proteinopathy-associated genetic variants), 103 ALS, and 47 ALS-FTD patients with likely TDP-43. A Subtype and Stage Inference (SuStaIn) model was trained to identify subtypes of individuals along the ALS-FTD spectrum with distinct brain atrophy patterns, and we related subtypes and stages to clinical, genetic, and neuropathological features of disease. Results: SuStaIn identified three novel subtypes: two disease subtypes with predominant brain atrophy either in prefrontal/somatomotor regions or limbic-related regions, and a normal-appearing group without obvious brain atrophy. The Limbic-predominant subtype tended to present with more impaired cognition, higher frequencies of pathogenic variants in TBK1 and TARDBP genes, and a higher proportion of TDP-43 type B, E and C. In contrast, the Prefrontal/Somatomotor-predominant subtype had higher frequencies of pathogenic variants in C9orf72 and GRN genes and higher proportion of TDP-43 type A. The normal-appearing brain group showed higher frequency of ALS relative to ALS-FTD and bvFTD patients, higher cognitive capacity, higher proportion of lower motor neuron onset, milder motor symptoms, and lower frequencies of genetic pathogenic variants. Overall SuStaIn stages also correlated with evidence for clinical progression including longer disease duration, higher King's stage, and cognitive decline. Additionally, SuStaIn stages differed across clinical phenotypes, genotypes and types of TDP-43 pathology. Conclusions: Our findings suggest distinct neurodegenerative subtypes of disease along the ALS-FTD spectrum that can be identified in vivo, each with distinct brain atrophy, clinical, genetic and pathological patterns.

Establishing the Validity of Compressed Sensing Diffusion Spectrum Imaging.

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of twenty-six participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n=20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

Rates of longitudinal change in 18 F-flortaucipir PET vary by brain region, cognitive impairment, and age in atypical Alzheimer's disease.

INTRODUCTION: Longitudinal positron emission tomography (PET) studies of tau accumulation in Alzheimer's disease (AD) have noted reduced increases or frank decreases in tau signal. We investigated how such reductions related to analytical confounds and disease progression markers in atypical AD. METHODS: We assessed regional and interindividual variation in longitudinal change on 18 F-flortaucipir PET imaging in 24 amyloid beta (Aß)+ patients with atypical, early-onset amnestic or non-amnestic AD plus 62 Aß- and 132 Aß+ Alzheimer's Disease Neuroimaging Initiative (ADNI) participants. RESULTS: In atypical AD, 18 F-flortaucipir uptake slowed or declined over time in areas with high baseline signal and older, more impaired individuals. ADNI participants had reduced longitudinal change in early Braak stage regions relative to late-stage areas. DISCUSSION: Results suggested radioligand uptake plateaus or declines in advanced neurodegeneration. Further research should investigate whether results generalize to other radioligands and whether they relate to changes of the radioligand binding site structure or accessibility.

Mitigating site effects in covariance for machine learning in neuroimaging data.

To acquire larger samples for answering complex questions in neuroscience, researchers have increasingly turned to multi-site neuroimaging studies. However, these studies are hindered by differences in images acquired across multiple sites. These effects have been shown to bias comparison between sites, mask biologically meaningful associations, and even introduce spurious associations. To address this, the field has focused on harmonizing data by removing site-related effects in the mean and variance of measurements. Contemporaneously with the increase in popularity of multi-center imaging, the use of machine learning (ML) in neuroimaging has also become commonplace. These approaches have been shown to provide improved sensitivity, specificity, and power due to their modeling the joint relationship across measurements in the brain. In this work, we demonstrate that methods for removing site effects in mean and variance may not be sufficient for ML. This stems from the fact that such methods fail to address how correlations between measurements can vary across sites. Data from the Alzheimer's Disease Neuroimaging Initiative is used to show that considerable differences in covariance exist across sites and that popular harmonization techniques do not address this issue. We then propose a novel harmonization method called Correcting Covariance Batch Effects (CovBat) that removes site effects in mean, variance, and covariance. We apply CovBat and show that within-site correlation matrices are successfully harmonized. Furthermore, we find that ML methods are unable to distinguish scanner manufacturer after our proposed harmonization is applied, and that the CovBat-harmonized data retain accurate prediction of disease group.

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.

The ANTsX ecosystem for quantitative biological and medical imaging.

The Advanced Normalizations Tools ecosystem, known as ANTsX, consists of multiple open-source software libraries which house top-performing algorithms used worldwide by scientific and research communities for processing and analyzing biological and medical imaging data. The base software library, ANTs, is built upon, and contributes to, the NIH-sponsored Insight Toolkit. Founded in 2008 with the highly regarded Symmetric Normalization image registration framework, the ANTs library has since grown to include additional functionality. Recent enhancements include statistical, visualization, and deep learning capabilities through interfacing with both the R statistical project (ANTsR) and Python (ANTsPy). Additionally, the corresponding deep learning extensions ANTsRNet and ANTsPyNet (built on the popular TensorFlow/Keras libraries) contain several popular network architectures and trained models for specific applications. One such comprehensive application is a deep learning analog for generating cortical thickness data from structural T1-weighted brain MRI, both cross-sectionally and longitudinally. These pipelines significantly improve computational efficiency and provide comparable-to-superior accuracy over multiple criteria relative to the existing ANTs workflows and simultaneously illustrate the importance of the comprehensive ANTsX approach as a framework for medical image analysis.

Transdiagnostic dimensions of psychopathology explain individuals' unique deviations from normative neurodevelopment in brain structure.

Psychopathology is rooted in neurodevelopment. However, clinical and biological heterogeneity, together with a focus on case-control approaches, have made it difficult to link dimensions of psychopathology to abnormalities of neurodevelopment. Here, using the Philadelphia Neurodevelopmental Cohort, we built normative models of cortical volume and tested whether deviations from these models better predicted psychiatric symptoms compared to raw cortical volume. Specifically, drawing on the p-factor hypothesis, we distilled 117 clinical symptom measures into six orthogonal psychopathology dimensions: overall psychopathology, anxious-misery, externalizing disorders, fear, positive psychosis symptoms, and negative psychosis symptoms. We found that multivariate patterns of deviations yielded improved out-of-sample prediction of psychopathology dimensions compared to multivariate patterns of raw cortical volume. We also found that correlations between overall psychopathology and deviations in ventromedial prefrontal, inferior temporal, and dorsal anterior cingulate cortices were stronger than those observed for specific dimensions of psychopathology (e.g., anxious-misery). Notably, these same regions are consistently implicated in a range of putatively distinct disorders. Finally, we performed conventional case-control comparisons of deviations in a group of individuals with depression and a group with attention-deficit hyperactivity disorder (ADHD). We observed spatially overlapping effects between these groups that diminished when controlling for overall psychopathology. Together, our results suggest that modeling cortical brain features as deviations from normative neurodevelopment improves prediction of psychiatric symptoms in out-of-sample testing, and that p-factor models of psychopathology may assist in separating biomarkers that are disorder-general from those that are disorder-specific.

Structural brain measures linked to clinical phenotypes in major depression replicate across clinical centres.

Abnormalities in brain structural measures, such as cortical thickness and subcortical volumes, are observed in patients with major depressive disorder (MDD) who also often show heterogeneous clinical features. This study seeks to identify the multivariate associations between structural phenotypes and specific clinical symptoms, a novel area of investigation. T1-weighted magnetic resonance imaging measures were obtained using 3 T scanners for 178 unmedicated depressed patients at four academic medical centres. Cortical thickness and subcortical volumes were determined for the depressed patients and patients' clinical presentation was characterized by 213 item-level clinical measures, which were grouped into several large, homogeneous categories by K-means clustering. The multivariate correlations between structural and cluster-level clinical-feature measures were examined using canonical correlation analysis (CCA) and confirmed with both 5-fold and leave-one-site-out cross-validation. Four broad types of clinical measures were detected based on clustering: an anxious misery composite (composed of item-level depression, anxiety, anhedonia, neuroticism and suicidality scores); positive personality traits (extraversion, openness, agreeableness and conscientiousness); reported history of physical/emotional trauma; and a reported history of sexual abuse. Responses on the item-level anxious misery measures were negatively associated with cortical thickness/subcortical volumes in the limbic system and frontal lobe; reported childhood history of physical/emotional trauma and sexual abuse measures were negatively correlated with entorhinal thickness and left hippocampal volume, respectively. In contrast, the positive traits measures were positively associated with hippocampal and amygdala volumes and cortical thickness of the highly-connected precuneus and cingulate cortex. Our findings suggest that structural brain measures may reflect neurobiological mechanisms underlying MDD features.

Longitudinal ComBat: A method for harmonizing longitudinal multi-scanner imaging data.

While aggregation of neuroimaging datasets from multiple sites and scanners can yield increased statistical power, it also presents challenges due to systematic scanner effects. This unwanted technical variability can introduce noise and bias into estimation of biological variability of interest. We propose a method for harmonizing longitudinal multi-scanner imaging data based on ComBat, a method originally developed for genomics and later adapted to cross-sectional neuroimaging data. Using longitudinal cortical thickness measurements from 663 participants in the Alzheimer's Disease Neuroimaging Initiative (ADNI) study, we demonstrate the presence of additive and multiplicative scanner effects in various brain regions. We compare estimates of the association between diagnosis and change in cortical thickness over time using three versions of the ADNI data: unharmonized data, data harmonized using cross-sectional ComBat, and data harmonized using longitudinal ComBat. In simulation studies, we show that longitudinal ComBat is more powerful for detecting longitudinal change than cross-sectional ComBat and controls the type I error rate better than unharmonized data with scanner included as a covariate. The proposed method would be useful for other types of longitudinal data requiring harmonization, such as genomic data, or neuroimaging studies of neurodevelopment, psychiatric disorders, or other neurological diseases.

Leveraging multi-shell diffusion for studies of brain development in youth and young adulthood.

Diffusion weighted imaging (DWI) has advanced our understanding of brain microstructure evolution over development. Recently, the use of multi-shell diffusion imaging sequences has coincided with advances in modeling the diffusion signal, such as Neurite Orientation Dispersion and Density Imaging (NODDI) and Laplacian-regularized Mean Apparent Propagator MRI (MAPL). However, the relative utility of recently-developed diffusion models for understanding brain maturation remains sparsely investigated. Additionally, despite evidence that motion artifact is a major confound for studies of development, the vulnerability of metrics derived from contemporary models to in-scanner motion has not been described. Accordingly, in a sample of 120 youth and young adults (ages 12-30) we evaluated metrics derived from diffusion tensor imaging (DTI), NODDI, and MAPL for associations with age and in-scanner head motion at multiple scales. Specifically, we examined mean white matter values, white matter tracts, white matter voxels, and connections in structural brain networks. Our results revealed that multi-shell diffusion imaging data can be leveraged to robustly characterize neurodevelopment, and demonstrate stronger age effects than equivalent single-shell data. Additionally, MAPL-derived metrics were less sensitive to the confounding effects of head motion. Our findings suggest that multi-shell imaging data and contemporary modeling techniques confer important advantages for studies of neurodevelopment.

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.