Regression and alignment for functional data and network topology.

In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of preprocessing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.

ComBatLS: A location- and scale-preserving method for multi-site image harmonization.

Recent work has leveraged massive datasets and advanced harmonization methods to construct normative models of neuroanatomical features and benchmark individuals' morphology. However, current harmonization tools do not preserve the effects of biological covariates including sex and age on features' variances; this failure may induce error in normative scores, particularly when such factors are distributed unequally across sites. Here, we introduce a new extension of the popular ComBat harmonization method, ComBatLS, that preserves biological variance in features' locations and scales. We use UK Biobank data to show that ComBatLS robustly replicates individuals' normative scores better than other ComBat methods when subjects are assigned to sex-imbalanced synthetic "sites". Additionally, we demonstrate that ComBatLS significantly reduces sex biases in normative scores compared to traditional methods. Finally, we show that ComBatLS successfully harmonizes consortium data collected across over 50 studies. R implementation of ComBatLS is available at https://github.com/andy1764/ComBatFamily.

Brain aging patterns in a large and diverse cohort of 49,482 individuals.

Brain aging process is influenced by various lifestyle, environmental and genetic factors, as well as by age-related and often coexisting pathologies. Magnetic resonance imaging and artificial intelligence methods have been instrumental in understanding neuroanatomical changes that occur during aging. Large, diverse population studies enable identifying comprehensive and representative brain change patterns resulting from distinct but overlapping pathological and biological factors, revealing intersections and heterogeneity in affected brain regions and clinical phenotypes. Herein, we leverage a state-of-the-art deep-representation learning method, Surreal-GAN, and present methodological advances and extensive experimental results elucidating brain aging heterogeneity in a cohort of 49,482 individuals from 11 studies. Five dominant patterns of brain atrophy were identified and quantified for each individual by respective measures, R-indices. Their associations with biomedical, lifestyle and genetic factors provide insights into the etiology of observed variances, suggesting their potential as brain endophenotypes for genetic and lifestyle risks. Furthermore, baseline R-indices predict disease progression and mortality, capturing early changes as supplementary prognostic markers. These R-indices establish a dimensional approach to measuring aging trajectories and related brain changes. They hold promise for precise diagnostics, especially at preclinical stages, facilitating personalized patient management and targeted clinical trial recruitment based on specific brain endophenotypic expression and prognosis.

Structure-function coupling in macroscale human brain networks.

Precisely how the anatomical structure of the brain gives rise to a repertoire of complex functions remains incompletely understood. A promising manifestation of this mapping from structure to function is the dependency of the functional activity of a brain region on the underlying white matter architecture. Here, we review the literature examining the macroscale coupling between structural and functional connectivity, and we establish how this structure-function coupling (SFC) can provide more information about the underlying workings of the brain than either feature alone. We begin by defining SFC and describing the computational methods used to quantify it. We then review empirical studies that examine the heterogeneous expression of SFC across different brain regions, among individuals, in the context of the cognitive task being performed, and over time, as well as its role in fostering flexible cognition. Last, we investigate how the coupling between structure and function is affected in neurological and psychiatric conditions, and we report how aberrant SFC is associated with disease duration and disease-specific cognitive impairment. By elucidating how the dynamic relationship between the structure and function of the brain is altered in the presence of neurological and psychiatric conditions, we aim to not only further our understanding of their aetiology but also establish SFC as a new and sensitive marker of disease symptomatology and cognitive performance. Overall, this Review collates the current knowledge regarding the regional interdependency between the macroscale structure and function of the human brain in both neurotypical and neuroatypical individuals.

PARE: A framework for removal of confounding effects from any distance-based dimension reduction method.

Dimension reduction tools preserving similarity and graph structure such as t-SNE and UMAP can capture complex biological patterns in high-dimensional data. However, these tools typically are not designed to separate effects of interest from unwanted effects due to confounders. We introduce the partial embedding (PARE) framework, which enables removal of confounders from any distance-based dimension reduction method. We then develop partial t-SNE and partial UMAP and apply these methods to genomic and neuroimaging data. For lower-dimensional visualization, our results show that the PARE framework can remove batch effects in single-cell sequencing data as well as separate clinical and technical variability in neuroimaging measures. We demonstrate that the PARE framework extends dimension reduction methods to highlight biological patterns of interest while effectively removing confounding effects.

A network control theory pipeline for studying the dynamics of the structural connectome.

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Psilocybin desynchronizes the human brain.

A single dose of psilocybin, a psychedelic that acutely causes distortions of space-time perception and ego dissolution, produces rapid and persistent therapeutic effects in human clinical trials1-4. In animal models, psilocybin induces neuroplasticity in cortex and hippocampus5-8. It remains unclear how human brain network changes relate to subjective and lasting effects of psychedelics. Here we tracked individual-specific brain changes with longitudinal precision functional mapping (roughly 18 magnetic resonance imaging visits per participant). Healthy adults were tracked before, during and for 3 weeks after high-dose psilocybin (25 mg) and methylphenidate (40 mg), and brought back for an additional psilocybin dose 6-12 months later. Psilocybin massively disrupted functional connectivity (FC) in cortex and subcortex, acutely causing more than threefold greater change than methylphenidate. These FC changes were driven by brain desynchronization across spatial scales (areal, global), which dissolved network distinctions by reducing correlations within and anticorrelations between networks. Psilocybin-driven FC changes were strongest in the default mode network, which is connected to the anterior hippocampus and is thought to create our sense of space, time and self. Individual differences in FC changes were strongly linked to the subjective psychedelic experience. Performing a perceptual task reduced psilocybin-driven FC changes. Psilocybin caused persistent decrease in FC between the anterior hippocampus and default mode network, lasting for weeks. Persistent reduction of hippocampal-default mode network connectivity may represent a neuroanatomical and mechanistic correlate of the proplasticity and therapeutic effects of psychedelics.

Harmonization of Longitudinal Diffusion Tensor Imaging Data of the Pediatric Cervical and Thoracic Spinal Cord at 3T Using Longitudinal ComBat.

Diffusion tensor imaging (DTI) of the spinal cord has been extensively used to identify biomarkers for spinal cord pathology. Previously, the longitudinal ComBat (longComBat) technique was examined to reduce scanner effects in multi-site, multi-scanner spinal cord DTI data. This study aimed to assess its effectiveness on longitudinal scans using a single-scanner pediatric dataset, including healthy and spinal cord injury (SCI) subjects. Two identical datasets were collected from 42 healthy and 27 SCI subjects with a 2-hour interval between scans on a 3T Siemens MRI scanner. Axial DTI images of the entire cervical and thoracic spinal cord were obtained, and various average diffusion tensor metrics (FA, MD, RD, & AD) were measured at each vertebral level. Pearson correlation and intraclass correlation coefficients were used to evaluate inter- and intra-subject agreement pre- and post-harmonization. Minimal improvement in agreement was observed with the mean square residual (MSR) model, while the restricted maximum likelihood estimator (REML) model showed reduced intra-subject agreement in all the tensor metrics. The significant variability between longitudinal DTI scans within a single scanner was likely due to physiological motion rather than scanner effects. Post-harmonization using the longComBat MSR model showed limited improvement in agreement.

Neural Responses to Auditory Food Stimuli Following Cognitive Behavioral Therapy for Binge-Eating Disorder.

OBJECTIVE: Adults with binge-eating disorder (BED), compared with those without BED, demonstrate higher blood-oxygen-level-dependent (BOLD) response to food cues in reward-related regions of the brain. It is not known whether cognitive behavioral therapy (CBT) can reverse this reward system hyperactivation. This randomized controlled trial (RCT) assessed changes in BOLD response to binge-eating cues following CBT versus wait-list control (WLC). METHOD: Females with BED (N = 40) were randomized to CBT or WLC. Participants completed assessments at baseline and 16 weeks including measures of eating and appetite and functional magnetic resonance imaging (fMRI) to measure BOLD response while listening to personalized scripts of binge-eating and neutral-relaxing cues. Data were analyzed using general linear models with mixed effects. RESULTS: Overall retention rate was 87.5%. CBT achieved significantly greater reductions in binge-eating episodes than WLC (mean ± standard error decline of 14.6 ± 2.7 vs. 5.7 ± 2.8 episodes in the past 28 days, respectively; p = 0.03). CBT and WLC did not differ significantly in changes in neural responses to binge-eating stimuli during the fMRI sessions. Compared with WLC, CBT had significantly greater improvements in reward-based eating drive, disinhibition, and hunger as assessed by questionnaires (ps < 0.05). DISCUSSION: CBT was effective in reducing binge eating, but, contrary to our hypothesis, CBT did not improve BOLD response to auditory binge-eating stimuli in reward regions of the brain. Further studies are needed to assess mechanisms underlying improvements with CBT for BED. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT03604172.

DeepComBat: A statistically motivated, hyperparameter-robust, deep learning approach to harmonization of neuroimaging data.

Neuroimaging data acquired using multiple scanners or protocols are increasingly available. However, such data exhibit technical artifacts across batches which introduce confounding and decrease reproducibility. This is especially true when multi-batch data are analyzed using complex downstream models which are more likely to pick up on and implicitly incorporate batch-related information. Previously proposed image harmonization methods have sought to remove these batch effects; however, batch effects remain detectable in the data after applying these methods. We present DeepComBat, a deep learning harmonization method based on a conditional variational autoencoder and the ComBat method. DeepComBat combines the strengths of statistical and deep learning methods in order to account for the multivariate relationships between features while simultaneously relaxing strong assumptions made by previous deep learning harmonization methods. As a result, DeepComBat can perform multivariate harmonization while preserving data structure and avoiding the introduction of synthetic artifacts. We apply this method to cortical thickness measurements from a cognitive-aging cohort and show DeepComBat qualitatively and quantitatively outperforms existing methods in removing batch effects while preserving biological heterogeneity. Additionally, DeepComBat provides a new perspective for statistically motivated deep learning harmonization methods.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Computer-vision analysis of craniofacial dysmorphology in 22q11.2 deletion syndrome and psychosis spectrum disorders.

BACKGROUND: Minor physical anomalies (MPAs) are congenital morphological abnormalities linked to disruptions of fetal development. MPAs are common in 22q11.2 deletion syndrome (22q11DS) and psychosis spectrum disorders (PS) and likely represent a disruption of early embryologic development that may help identify overlapping mechanisms linked to psychosis in these disorders. METHODS: Here, 2D digital photographs were collected from 22q11DS (n = 150), PS (n = 55), and typically developing (TD; n = 93) individuals. Photographs were analyzed using two computer-vision techniques: (1) DeepGestalt algorithm (Face2Gene (F2G)) technology to identify the presence of genetically mediated facial disorders, and (2) Emotrics-a semi-automated machine learning technique that localizes and measures facial features. RESULTS: F2G reliably identified patients with 22q11DS; faces of PS patients were matched to several genetic conditions including FragileX and 22q11DS. PCA-derived factor loadings of all F2G scores indicated unique and overlapping facial patterns that were related to both 22q11DS and PS. Regional facial measurements of the eyes and nose were smaller in 22q11DS as compared to TD, while PS showed intermediate measurements. CONCLUSIONS: The extent to which craniofacial dysmorphology 22q11DS and PS overlapping and evident before the impairment or distress of sub-psychotic symptoms may allow us to identify at-risk youths more reliably and at an earlier stage of development.

Multivariate testing and effect size measures for batch effect evaluation in radiomic features.

While precision medicine applications of radiomics analysis are promising, differences in image acquisition can cause "batch effects" that reduce reproducibility and affect downstream predictive analyses. Harmonization methods such as ComBat have been developed to correct these effects, but evaluation methods for quantifying batch effects are inconsistent. In this study, we propose the use of the multivariate statistical test PERMANOVA and the Robust Effect Size Index (RESI) to better quantify and characterize batch effects in radiomics data. We evaluate these methods in both simulated and real radiomics features extracted from full-field digital mammography (FFDM) data. PERMANOVA demonstrated higher power than standard univariate statistical testing, and RESI was able to interpretably quantify the effect size of site at extremely large sample sizes. These methods show promise as more powerful and interpretable methods for the detection and quantification of batch effects in radiomics studies.

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.

Assessing the Stability of Photon-Counting CT: Insights from a Two-Year Longitudinal Study.

Background: Among the advancements in computed tomography (CT) technology, photon-counting computed tomography (PCCT) stands out as a significant innovation, providing superior spectral imaging capabilities while simultaneously reducing radiation exposure. Its long-term stability is important for clinical care, especially longitudinal studies, but is currently unknown. Purpose: This study sets out to comprehensively analyze the long-term stability of a first-generation clinical PCCT scanner. Materials and Methods: Over a two-year period, from November 2021 to November 2023, we conducted weekly identical experiments utilizing the same multi-energy CT protocol. These experiments included various tissue-mimicking inserts to rigorously assess the stability of Hounsfield Units (HU) and image noise in Virtual Monochromatic Images (VMIs) and iodine density maps. Throughout this period, notable software and hardware modifications were meticulously recorded. Each week, VMIs and iodine density maps were reconstructed and analyzed to evaluate quantitative stability over time. Results: Spectral results consistently demonstrated the quantitative stability of PCCT. VMIs exhibited stable HU values, such as variation in relative error for VMI 70 keV measuring 0.11% and 0.30% for single-source and dual-source modes, respectively. Similarly, noise levels remained stable with slight fluctuations linked to software changes for VMI 40 and 70 keV that corresponded to changes of 8 and 1 HU, respectively. Furthermore, iodine density quantification maintained stability and showed significant improvement with software and hardware changes, especially in dual-source mode with nominal errors decreasing from 1.44 to 0.03 mg/mL. Conclusion: This study provides the first long-term reproducibility assessment of quantitative PCCT imaging, highlighting its potential for the clinical arena. This study indicates its long-term utility in diagnostic radiology, especially for longitudinal studies.

pNet: A toolbox for personalized functional networks modeling.

Personalized functional networks (FNs) derived from functional magnetic resonance imaging (fMRI) data are useful for characterizing individual variations in the brain functional topography associated with the brain development, aging, and disorders. To facilitate applications of the personalized FNs with enhanced reliability and reproducibility, we develop an open-source toolbox that is user-friendly, extendable, and includes rigorous quality control (QC), featuring multiple user interfaces (graphics, command line, and a step-by-step guideline) and job-scheduling for high performance computing (HPC) clusters. Particularly, the toolbox, named personalized functional network modeling (pNet), takes fMRI inputs in either volumetric or surface type, ensuring compatibility with multiple fMRI data formats, and computes personalized FNs using two distinct modeling methods: one method optimizes the functional coherence of FNs, while the other enhances their independence. Additionally, the toolbox provides HTML-based reports for QC and visualization of personalized FNs. The toolbox is developed in both MATLAB and Python platforms with a modular design to facilitate extension and modification by users familiar with either programming language. We have evaluated the toolbox on two fMRI datasets and demonstrated its effectiveness and user-friendliness with interactive and scripting examples. pNet is publicly available at https://github.com/MLDataAnalytics/pNet.

Diagnosis concealment behaviors and disclosure beliefs are associated with health and quality of life in people with multiple sclerosis.

BACKGROUND: People with multiple sclerosis (pwMS) struggle with whether, how, and how much to disclose their diagnosis. They often expend resources to conceal their diagnosis, and hold beliefs that it may negatively affect their personal relationships and/or professional opportunities. To better understand these effects, we developed a measure to quantify concealment behaviors and disclosure beliefs. Our main objective is to evaluate relationships of DISCO-MS responses to health and quality of life in a multinational cohort. METHODS: Survey responses were obtained for DISCO-MS and PROMIS-MS scales: global health, communication, social roles participation, anxiety, depression, emotional / behavioral dyscontrol, fatigue, lower extremity function, positive affect / well-being, social roles satisfaction, sleep, stigma, upper extremity function, cognitive function, bladder control, bowel control, visual function. Simple linear regression assessed associations. RESULTS: 263 pwMS were include. Higher concealment was associated with higher anxiety (beta= 0.15 [0.07, 0.23]), depression (beta = 0.13 [0.05, 0.21]), emotional dyscontrol (beta = 0.12 [0.04, 0.20]), lower affect / well-being (beta = -0.13 [-0.21, - 0.05]). Higher anticipation of negative consequences of disclosure was associated with lower self-reported physical (beta = -0.15) and mental health (beta = -0.14), lower positive affect / well-being, social roles satisfaction, higher anxiety, depression, emotional dyscontrol, sleep disturbance, and higher perceived stigma. DISCUSSION: These results reveal potential consequences of diagnosis concealment for physical and mental health and quality of life. Raising awareness and implementing interventions may mitigate negative repercussions of concealment.

A general exposome factor explains individual differences in functional brain network topography and cognition in youth.

Childhood environments are critical in shaping cognitive neurodevelopment. With the increasing availability of large-scale neuroimaging datasets with deep phenotyping of childhood environments, we can now build upon prior studies that have considered relationships between one or a handful of environmental and neuroimaging features at a time. Here, we characterize the combined effects of hundreds of inter-connected and co-occurring features of a child's environment ("exposome") and investigate associations with each child's unique, multidimensional pattern of functional brain network organization ("functional topography") and cognition. We apply data-driven computational models to measure the exposome and define personalized functional brain networks in pre-registered analyses. Across matched discovery (n=5139, 48.5% female) and replication (n=5137, 47.1% female) samples from the Adolescent Brain Cognitive Development study, the exposome was associated with current (ages 9-10) and future (ages 11-12) cognition. Changes in the exposome were also associated with changes in cognition after accounting for baseline scores. Cross-validated ridge regressions revealed that the exposome is reflected in functional topography and can predict performance across cognitive domains. Importantly, a single measure capturing a child's exposome could more accurately and parsimoniously predict cognition than a wealth of personalized neuroimaging data, highlighting the importance of children's complex, multidimensional environments in cognitive neurodevelopment.

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.

Disparities in Genetic Testing for Neurologic Disorders.

BACKGROUND AND OBJECTIVES: Genetic testing is now the standard of care for many neurologic conditions. Health care disparities are unfortunately widespread in the US health care system, but disparities in the utilization of genetic testing for neurologic conditions have not been studied. We tested the hypothesis that access to and results of genetic testing vary according to race, ethnicity, sex, socioeconomic status, and insurance status for adults with neurologic conditions. METHODS: We analyzed retrospective data from patients who underwent genetic evaluation and testing through our institution's neurogenetics program. We tested for differences between demographic groups in 3 steps of a genetic evaluation pathway: (1) attending a neurogenetic evaluation, (2) completing genetic testing, and (3) receiving a diagnostic result. We compared patients on this genetic evaluation pathway with the population of all neurology outpatients at our institution, using univariate and multivariable logistic regression analyses. RESULTS: Between 2015 and 2022, a total of 128,440 patients were seen in our outpatient neurology clinics and 2,540 patients underwent genetic evaluation. Black patients were less than half as likely as White patients to be evaluated (odds ratio [OR] 0.49, p < 0.001), and this disparity was similar after controlling for other demographic factors in multivariable analysis. Patients from the least wealthy quartile of zip codes were also less likely to be evaluated (OR 0.67, p < 0.001). Among patients who underwent evaluation, there were no disparities in the likelihood of completing genetic testing, nor in the likelihood of a diagnostic result after adjusting for age. Analyses restricted to specific indications for genetic testing supported these findings. DISCUSSION: We observed unequal utilization of our clinical neurogenetics program for patients from marginalized and minoritized demographic groups, especially Black patients. Among patients who do undergo evaluation, all groups benefit similarly from genetic testing when it is indicated. Understanding and removing barriers to accessing genetic testing will be essential to health care equity and optimal care for all patients with neurologic disorders.