Effect of patient factors, center, and era on Fontan timing: An observational study using the Pediatric Health Information Systems Database

BACKGROUND

There are no consensus guidelines defining optimal timing for the Fontan operation, the last planned surgery in staged palliation for single-ventricle heart disease.

OBJECTIVES

Identify patient-level characteristics, center-level variation, and secular trends driving Fontan timing.

METHODS

A retrospective observational study of subjects who underwent Fontan from 2007 to 2021 at centers in the Pediatric Health Information Systems database was performed using linear mixed-effects modeling in which age at Fontan was regressed on patient characteristics and date of operation with center as random effect.

RESULTS

We included 10,305 subjects (40.4% female, 44% non-white) at 47 centers. Median age at Fontan was 3.4 years (IQR 2.6-4.4). Hypoplastic left heart syndrome (-4.4 months, 95%CI -5.5 to -3.3) and concomitant conditions (-2.6 months, 95%CI -4.1 to -1.1) were associated with younger age at Fontan. Subjects with technology-dependence (+4.6 months, 95%CI 3.1-6.1) were older at Fontan. Black (+4.1 months, 95%CI 2.5-5.7) and Asian (+8.3 months, 95%CI 5.4-11.2) race were associated with older age at Fontan. There was significant variation in Fontan timing between centers. Center accounted for 10% of variation (ICC 0.10, 95%CI 0.07-0.14). Center surgical volume was not associated with Fontan timing (P = .21). Operation year was associated with age at Fontan, with a 3.1 month increase in age for every 5 years (+0.61 months, 95%CI 0.48-0.75).

CONCLUSIONS

After adjusting for patient-level characteristics there remains significant inter-center variation in Fontan timing. Age at Fontan has increased. Future studies addressing optimal Fontan timing are warranted.

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.

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.

Spatial nonstationarity of image noise in widefield optical imaging and its effects on cluster-based inference for resting-state functional connectivity

BACKGROUND

Resting-state functional connectivity (RSFC) analysis with widefield optical imaging (WOI) is a potentially powerful tool to develop imaging biomarkers in mouse models of disease before translating them to human neuroimaging with functional magnetic resonance imaging (fMRI). The delineation of such biomarkers depends on rigorous statistical analysis. However, statistical understanding of WOI data is limited. In particular, cluster-based analysis of neuroimaging data depends on assumptions of spatial stationarity (i.e., that the distribution of cluster sizes under the null is equal at all brain locations). Whether actual data deviate from this assumption has not previously been examined in WOI.

NEW METHOD

In this manuscript, we characterize the effects of spatial nonstationarity in WOI RSFC data and adapt a "two-pass" technique from fMRI to correct cluster sizes and mitigate spatial bias, both parametrically and nonparametrically. These methods are tested on multi-institutional data.

RESULTS AND COMPARISON WITH EXISTING METHODS

We find that spatial nonstationarity has a substantial effect on inference in WOI RSFC data with false positives much more likely at some brain regions than others. This pattern of bias varies between imaging systems, contrasts, and mouse ages, all of which could affect experimental reproducibility if not accounted for.

CONCLUSIONS

Both parametric and nonparametric corrections for nonstationarity result in significant improvements in spatial bias. The proposed methods are simple to implement and will improve the robustness of inference in optical neuroimaging data.

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.

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.

Disparities in seizure outcomes revealed by large language models

OBJECTIVE

Large-language models (LLMs) can potentially revolutionize health care delivery and research, but risk propagating existing biases or introducing new ones. In epilepsy, social determinants of health are associated with disparities in care access, but their impact on seizure outcomes among those with access remains unclear. Here we (1) evaluated our validated, epilepsy-specific LLM for intrinsic bias, and (2) used LLM-extracted seizure outcomes to determine if different demographic groups have different seizure outcomes.

MATERIALS AND METHODS

We tested our LLM for differences and equivalences in prediction accuracy and confidence across demographic groups defined by race, ethnicity, sex, income, and health insurance, using manually annotated notes. Next, we used LLM-classified seizure freedom at each office visit to test for demographic outcome disparities, using univariable and multivariable analyses.

RESULTS

We analyzed 84 675 clinic visits from 25 612 unique patients seen at our epilepsy center. We found little evidence of bias in the prediction accuracy or confidence of outcome classifications across demographic groups. Multivariable analysis indicated worse seizure outcomes for female patients (OR 1.33, P ≤ .001), those with public insurance (OR 1.53, P ≤ .001), and those from lower-income zip codes (OR ≥1.22, P  ≤ .007). Black patients had worse outcomes than White patients in univariable but not multivariable analysis (OR 1.03, P = .66).

CONCLUSION

We found little evidence that our LLM was intrinsically biased against any demographic group. Seizure freedom extracted by LLM revealed disparities in seizure outcomes across several demographic groups. These findings quantify the critical need to reduce disparities in the care of people with epilepsy.

Generalizable Links Between Borderline Personality Traits and Functional Connectivity

BACKGROUND

Symptoms of borderline personality disorder (BPD) often manifest during adolescence, but the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here, we aimed to investigate how multivariate patterns of functional connectivity are associated with borderline personality traits in large samples of young adults and adolescents.

METHODS

We used functional magnetic resonance imaging data from young adults and adolescents from the HCP-YA (Human Connectome Project Young Adult) (n = 870, ages 22-37 years, 457 female) and the HCP-D (Human Connectome Project Development) (n = 223, ages 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five-Factor Inventory. A ridge regression model with cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity.

RESULTS

Multivariate functional connectivity patterns significantly predicted out-of-sample BPD scores in unseen data in young adults (HCP-YA ppermuted = .001) and older adolescents (HCP-D ppermuted = .001). Regional predictive capacity was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD scores aligned with those associated with development in youth.

CONCLUSIONS

Individual differences in functional connectivity in developmentally sensitive regions are associated with borderline personality traits.

Impact of white matter hyperintensities on structural connectivity and cognition in cognitively intact ADNI participants

We used indirect brain mapping with virtual lesion tractography to test the hypothesis that the extent of white matter tract disconnection due to white matter hyperintensities (WMH) is associated with corresponding tract-specific cognitive performance decrements. To estimate tract disconnection, WMH masks were extracted from FLAIR MRI data of 481 cognitively intact participants in the Alzheimer's Disease Neuroimaging Initiative (ADNI) and used as regions of avoidance for fiber tracking in diffusion MRI data from 50 healthy young participants from the Human Connectome Project. Estimated tract disconnection in the right inferior fronto-occipital fasciculus, right frontal aslant tract, and right superior longitudinal fasciculus mediated the effects of WMH volume on executive function. Estimated tract disconnection in the left uncinate fasciculus mediated the effects of WMH volume on memory and in the right frontal aslant tract on language. In a subset of ADNI control participants with amyloid data, positive status increased the probability of periventricular WMH and moderated the relationship between WMH burden and tract disconnection in executive function performance.

Integrating imaging and genomic data for the discovery of distinct glioblastoma subtypes: a joint learning approach

Glioblastoma is a highly heterogeneous disease, with variations observed at both phenotypical and molecular levels. Personalized therapies would be facilitated by non-invasive in vivo approaches for characterizing this heterogeneity. In this study, we developed unsupervised joint machine learning between radiomic and genomic data, thereby identifying distinct glioblastoma subtypes. A retrospective cohort of 571 IDH-wildtype glioblastoma patients were included in the study, and pre-operative multi-parametric MRI scans and targeted next-generation sequencing (NGS) data were collected. L21-norm minimization was used to select a subset of 12 radiomic features from the MRI scans, and 13 key driver genes from the five main signal pathways most affected in glioblastoma were selected from the genomic data. Subtypes were identified using a joint learning approach called Anchor-based Partial Multi-modal Clustering on both radiomic and genomic modalities. Kaplan-Meier analysis identified three distinct glioblastoma subtypes: high-risk, medium-risk, and low-risk, based on overall survival outcome (p < 0.05, log-rank test; Hazard Ratio = 1.64, 95% CI 1.17-2.31, Cox proportional hazard model on high-risk and low-risk subtypes). The three subtypes displayed different phenotypical and molecular characteristics in terms of imaging histogram, co-occurrence of genes, and correlation between the two modalities. Our findings demonstrate the synergistic value of integrated radiomic signatures and molecular characteristics for glioblastoma subtyping. Joint learning on both modalities can aid in better understanding the molecular basis of phenotypical signatures of glioblastoma, and provide insights into the biological underpinnings of tumor formation and progression.

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 pre-processing 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.

Genetic and Clinical Correlates of AI-Based Brain Aging Patterns in Cognitively Unimpaired Individuals

Importance

Brain aging elicits complex neuroanatomical changes influenced by multiple age-related pathologies. Understanding the heterogeneity of structural brain changes in aging may provide insights into preclinical stages of neurodegenerative diseases.

Objective

To derive subgroups with common patterns of variation in participants without diagnosed cognitive impairment (WODCI) in a data-driven manner and relate them to genetics, biomedical measures, and cognitive decline trajectories.

Design, Setting, and Participants

Data acquisition for this cohort study was performed from 1999 to 2020. Data consolidation and harmonization were conducted from July 2017 to July 2021. Age-specific subgroups of structural brain measures were modeled in 4 decade-long intervals spanning ages 45 to 85 years using a deep learning, semisupervised clustering method leveraging generative adversarial networks. Data were analyzed from July 2021 to February 2023 and were drawn from the Imaging-Based Coordinate System for Aging and Neurodegenerative Diseases (iSTAGING) international consortium. Individuals WODCI at baseline spanning ages 45 to 85 years were included, with greater than 50 000 data time points.

Exposures

Individuals WODCI at baseline scan.

Main Outcomes and Measures

Three subgroups, consistent across decades, were identified within the WODCI population. Associations with genetics, cardiovascular risk factors (CVRFs), amyloid β (Aβ), and future cognitive decline were assessed.

Results

In a sample of 27 402 individuals (mean [SD] age, 63.0 [8.3] years; 15 146 female [55%]) WODCI, 3 subgroups were identified in contrast with the reference group: a typical aging subgroup, A1, with a specific pattern of modest atrophy and white matter hyperintensity (WMH) load, and 2 accelerated aging subgroups, A2 and A3, with characteristics that were more distinct at age 65 years and older. A2 was associated with hypertension, WMH, and vascular disease-related genetic variants and was enriched for Aβ positivity (ages ≥65 years) and apolipoprotein E (APOE) ε4 carriers. A3 showed severe, widespread atrophy, moderate presence of CVRFs, and greater cognitive decline. Genetic variants associated with A1 were protective for WMH (rs7209235: mean [SD] B = -0.07 [0.01]; P value = 2.31 × 10-9) and Alzheimer disease (rs72932727: mean [SD] B = 0.1 [0.02]; P value = 6.49 × 10-9), whereas the converse was observed for A2 (rs7209235: mean [SD] B = 0.1 [0.01]; P value = 1.73 × 10-15 and rs72932727: mean [SD] B = -0.09 [0.02]; P value = 4.05 × 10-7, respectively); variants in A3 were associated with regional atrophy (rs167684: mean [SD] B = 0.08 [0.01]; P value = 7.22 × 10-12) and white matter integrity measures (rs1636250: mean [SD] B = 0.06 [0.01]; P value = 4.90 × 10-7).

Conclusions and Relevance

The 3 subgroups showed distinct associations with CVRFs, genetics, and subsequent cognitive decline. These subgroups likely reflect multiple underlying neuropathologic processes and affect susceptibility to Alzheimer disease, paving pathways toward patient stratification at early asymptomatic stages and promoting precision medicine in clinical trials and health care.

Transcriptional cartography integrates multiscale biology of the human cortex

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Imaging chronic active lesions in multiple sclerosis: a consensus statement

Chronic active lesions (CAL) are an important manifestation of chronic inflammation in multiple sclerosis (MS) and have implications for non-relapsing biological progression. In recent years, the discovery of innovative magnetic resonance imaging (MRI) and PET derived biomarkers has made it possible to detect CAL, and to some extent quantify them, in the brain of persons with MS, in vivo. Paramagnetic rim lesions on susceptibility-sensitive MRI sequences, MRI-defined slowly expanding lesions on T1-weighted (T1-w) and T2-w scans, and 18-kDa translocator protein-positive lesions on PET are promising candidate biomarkers of CAL. While partially overlapping, these biomarkers do not have equivalent sensitivity and specificity to histopathological CAL. Standardization in the use of available imaging measures for CAL identification, quantification, and monitoring is lacking. To fast-forward clinical translation of CAL, the North American Imaging in Multiple Sclerosis Cooperative developed a Consensus Statement, which provides guidance for the radiological definition and measurement of CAL. The proposed manuscript presents this Consensus Statement, summarizes the multistep process leading to it, and identifies the remaining major gaps in knowledge.

A multicenter pilot study evaluating simplified central vein assessment for the diagnosis of multiple sclerosis

BACKGROUND

The central vein sign (CVS) is a proposed magnetic resonance imaging (MRI) biomarker for multiple sclerosis (MS); the optimal method for abbreviated CVS scoring is not yet established.

OBJECTIVE

The aim of this study was to evaluate the performance of a simplified approach to CVS assessment in a multicenter study of patients being evaluated for suspected MS.

METHODS

Adults referred for possible MS to 10 sites were recruited. A post-Gd 3D T2*-weighted MRI sequence (FLAIR*) was obtained in each subject. Trained raters at each site identified up to six CVS-positive lesions per FLAIR* scan. Diagnostic performance of CVS was evaluated for a diagnosis of MS which had been confirmed using the 2017 McDonald criteria at thresholds including three positive lesions (Select-3*) and six positive lesions (Select-6*). Inter-rater reliability assessments were performed.

RESULTS

Overall, 78 participants were analyzed; 37 (47%) were diagnosed with MS, and 41 (53%) were not. The mean age of participants was 45 (range: 19-64) years, and most were female (n = 55, 71%). The area under the receiver operating characteristic curve (AUROC) for the simplified counting method was 0.83 (95% CI: 0.73-0.93). Select-3* and Select-6* had sensitivity of 81% and 65% and specificity of 68% and 98%, respectively. Inter-rater agreement was 78% for Select-3* and 83% for Select-6*.

CONCLUSION

A simplified method for CVS assessment in patients referred for suspected MS demonstrated good diagnostic performance and inter-rater agreement.

Alcohol for seizure induction in the epilepsy monitoring unit

RATIONALE

Seizure induction techniques are used in the epilepsy monitoring unit (EMU) to increase diagnostic yield and reduce length of stay. There are insufficient data on the efficacy of alcohol as an induction technique.

METHODS

We performed a retrospective cohort study using six years of EMU data at our institution. We compared cases who received alcohol for seizure induction to matched controls who did not. The groups were matched on the following variables: age, reason for admission, length of stay, number of antiseizure medications (ASM) at admission, whether ASMs were tapered during admission, and presence of interictal epileptiform discharges. We used both propensity score and exact matching strategies. We compared the likelihood of epileptic seizures and nonepileptic events in cases versus controls using Kaplan-Meier time-to-event analysis, as well as odds ratios for these outcomes occurring at any time during the admission.

RESULTS

We analyzed 256 cases who received alcohol (median dose 2.5 standard drinks) and 256 propensity score-matched controls. Cases who received alcohol were no more likely than controls to have an epileptic seizure (X2(1) = 0.01, p = 0.93) or nonepileptic event (X2(1) = 2.1, p = 0.14) in the first 48 h after alcohol administration. For the admission overall, cases were no more likely to have an epileptic seizure (OR 0.89, 95 % CI 0.61-1.28, p = 0.58), nonepileptic event (OR 0.97, CI 0.62-1.53, p = 1.00), nor require rescue benzodiazepine (OR 0.63, CI 0.35-1.12, p = 0.15). Stratified analyses revealed no increased risk of epileptic seizure in any subgroups. Sensitivity analysis using exact matching showed that results were robust to matching strategy.

CONCLUSIONS

Alcohol was not an effective induction technique in the EMU. This finding has implications for counseling patients with epilepsy about the risks of drinking alcohol in moderation in their daily lives.

Five dominant dimensions of brain aging are identified via deep learning: associations with clinical, lifestyle, and genetic measures

Brain aging is a complex process influenced by various lifestyle, environmental, and genetic factors, as well as by age-related and often co-existing pathologies. MRI and, more recently, AI methods have been instrumental in understanding the neuroanatomical changes that occur during aging in large and diverse populations. However, the multiplicity and mutual overlap of both pathologic processes and affected brain regions make it difficult to precisely characterize the underlying neurodegenerative profile of an individual from an MRI scan. Herein, we leverage a state-of-the art deep representation learning method, Surreal-GAN, and present both methodological advances and extensive experimental results that allow us to elucidate the heterogeneity of brain aging in a large and diverse cohort of 49,482 individuals from 11 studies. Five dominant patterns of neurodegeneration were identified and quantified for each individual by their respective (herein referred to as) R-indices. Significant associations between R-indices and distinct biomedical, lifestyle, and genetic factors provide insights into the etiology of observed variances. Furthermore, baseline R-indices showed predictive value for disease progression and mortality. These five R-indices contribute to MRI-based precision diagnostics, prognostication, and may inform stratification into clinical trials.

Multi-contrast high-field quality image synthesis for portable low-field MRI using generative adversarial networks and paired data

Introduction

Portable low-field strength (64mT) MRI scanners promise to increase access to neuroimaging for clinical and research purposes, however these devices produce lower quality images compared to high-field scanners. In this study, we developed and evaluated a deep learning architecture to generate high-field quality brain images from low-field inputs using a paired dataset of multiple sclerosis (MS) patients scanned at 64mT and 3T.

Methods

A total of 49 MS patients were scanned on portable 64mT and standard 3T scanners at Penn (n=25) or the National Institutes of Health (NIH, n=24) with T1-weighted, T2-weighted and FLAIR acquisitions. Using this paired data, we developed a generative adversarial network (GAN) architecture for low- to high-field image translation (LowGAN). We then evaluated synthesized images with respect to image quality, brain morphometry, and white matter lesions.

Results

Synthetic high-field images demonstrated visually superior quality compared to low-field inputs and significantly higher normalized cross-correlation (NCC) to actual high-field images for T1 (p=0.001) and FLAIR (p<0.001) contrasts. LowGAN generally outperformed the current state-of-the-art for low-field volumetrics. For example, thalamic, lateral ventricle, and total cortical volumes in LowGAN outputs did not differ significantly from 3T measurements. Synthetic outputs preserved MS lesions and captured a known inverse relationship between total lesion volume and thalamic volume.

Conclusions

LowGAN generates synthetic high-field images with comparable visual and quantitative quality to actual high-field scans. Enhancing portable MRI image quality could add value and boost clinician confidence, enabling wider adoption of this technology.

Genomic loci influence patterns of structural covariance in the human brain

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

Personalized functional brain network topography is associated with individual differences in youth cognition

Individual differences in cognition during childhood are associated with important social, physical, and mental health outcomes in adolescence and adulthood. Given that cortical surface arealization during development reflects the brain's functional prioritization, quantifying variation in the topography of functional brain networks across the developing cortex may provide insight regarding individual differences in cognition. We test this idea by defining personalized functional networks (PFNs) that account for interindividual heterogeneity in functional brain network topography in 9-10 year olds from the Adolescent Brain Cognitive Development℠ Study. Across matched discovery (n = 3525) and replication (n = 3447) samples, the total cortical representation of fronto-parietal PFNs positively correlates with general cognition. Cross-validated ridge regressions trained on PFN topography predict cognition in unseen data across domains, with prediction accuracy increasing along the cortex's sensorimotor-association organizational axis. These results establish that functional network topography heterogeneity is associated with individual differences in cognition before the critical transition into adolescence.

CoCoA: conditional correlation models with association size

Many scientific questions can be formulated as hypotheses about conditional correlations. For instance, in tests of cognitive and physical performance, the trade-off between speed and accuracy motivates study of the two variables together. A natural question is whether speed-accuracy coupling depends on other variables, such as sustained attention. Classical regression techniques, which posit models in terms of covariates and outcomes, are insufficient to investigate the effect of a third variable on the symmetric relationship between speed and accuracy. In response, we propose a conditional correlation model with association size, a likelihood-based statistical framework to estimate the conditional correlation between speed and accuracy as a function of additional variables. We propose novel measures of the association size, which are analogous to effect sizes on the correlation scale while adjusting for confound variables. In simulation studies, we compare likelihood-based estimators of conditional correlation to semiparametric estimators adapted from genomic studies and find that the former achieves lower bias and variance under both ideal settings and model assumption misspecification. Using neurocognitive data from the Philadelphia Neurodevelopmental Cohort, we demonstrate that greater sustained attention is associated with stronger speed-accuracy coupling in a complex reasoning task while controlling for age. By highlighting conditional correlations as the outcome of interest, our model provides complementary insights to traditional regression modeling and partitioned correlation analyses.

Evaluation of resampling-based inference for topological features of neuroimages

Many recent studies have demonstrated the inflated type 1 error rate of the original Gaussian random field (GRF) methods for inference of neuroimages and identified resampling (permutation and bootstrapping) methods that have better performance. There has been no evaluation of resampling procedures when using robust (sandwich) statistical images with different topological features (TF) used for neuroimaging inference. Here, we consider estimation of distributions TFs of a statistical image and evaluate resampling procedures that can be used when exchangeability is violated. We compare the methods using realistic simulations and study sex differences in life-span age-related changes in gray matter volume in the Nathan Kline Institute Rockland sample. We find that our proposed wild bootstrap and the commonly used permutation procedure perform well in sample sizes above 50 under realistic simulations with heteroskedasticity. The Rademacher wild bootstrap has fewer assumptions than the permutation and performs similarly in samples of 100 or more, so is valid in a broader range of conditions. We also evaluate the GRF-based pTFCE method and show that it has inflated error rates in samples less than 200. Our R package, pbj , is available on Github and allows the user to reproducibly implement various resampling-based group level neuroimage analyses.

Similarity-based multimodal regression

To better understand complex human phenotypes, large-scale studies have increasingly collected multiple data modalities across domains such as imaging, mobile health, and physical activity. The properties of each data type often differ substantially and require either separate analyses or extensive processing to obtain comparable features for a combined analysis. Multimodal data fusion enables certain analyses on matrix-valued and vector-valued data, but it generally cannot integrate modalities of different dimensions and data structures. For a single data modality, multivariate distance matrix regression provides a distance-based framework for regression accommodating a wide range of data types. However, no distance-based method exists to handle multiple complementary types of data. We propose a novel distance-based regression model, which we refer to as Similarity-based Multimodal Regression (SiMMR), that enables simultaneous regression of multiple modalities through their distance profiles. We demonstrate through simulation, imaging studies, and longitudinal mobile health analyses that our proposed method can detect associations between clinical variables and multimodal data of differing properties and dimensionalities, even with modest sample sizes. We perform experiments to evaluate several different test statistics and provide recommendations for applying our method across a broad range of scenarios.

Trends in Ductus Arteriosus Stent Versus Blalock-Taussig-Thomas Shunt Use and Comparison of Cost, Length of Stay, and Short-Term Outcomes in Neonates With Ductal-Dependent Pulmonary Blood Flow: An Observational Study Using the Pediatric Health Information Systems Database

BACKGROUND

The modified Blalock-Taussig-Thomas shunt is the gold standard palliation for securing pulmonary blood flow in infants with ductal-dependent pulmonary blood flow. Recently, the ductus arteriosus stent (DAS) has become a viable alternative.

METHODS AND RESULTS

This was a retrospective multicenter study of neonates ≤30 days undergoing DAS or Blalock-Taussig-Thomas shunt placement between January 1, 2017 and December 31, 2020 at hospitals reporting to the Pediatric Health Information Systems database. We performed generalized linear mixed-effects modeling to evaluate trends in intervention and intercenter variation, propensity score adjustment and inverse probability weighting with linear mixed-effects modeling to analyze length of stay and cost of hospitalization, and generalized linear mixed modeling to analyze differences in 30-day outcomes. There were 1874 subjects (58% male, 61% White) from 45 centers (29% DAS). Odds of DAS increased with time (odds ratio [OR] 1.23, annually, P<0.01 [95% CI, 1.10-1.38]) with significant intercenter variation (median OR, 3.81 [95% CI, 2.74-5.91]). DAS was associated with shorter hospital length of stay (ratio of geometric means, 0.76 [95% CI, 0.63-0.91]), shorter intensive care unit length of stay (ratio of geometric means, 0.77 [95% CI, 0.61-0.97]), and less expensive hospitalization (ratio of geometric means, 0.70 [95% CI, 0.56-0.87]). Intervention was not significantly associated with odds of 30-day transplant-free survival (OR,1.18 [95% CI, 0.70-1.99]) or freedom from catheter reintervention (OR, 1.02 [95% CI, 0.65-1.58]), but DAS was associated with 30-day freedom from composite adverse outcome (OR, 1.51 [95% CI, 1.11-2.05]).

CONCLUSIONS

Use of DAS is increasing, but there is variability across centers. Though odds of transplant-free survival and reintervention were not significantly different after DAS, and DAS was associated with shorter length of stay and lower in-hospital costs.

XCP-D: A Robust Pipeline for the post-processing of fMRI data

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

Mapping the Relationship of White Matter Lesions to Depression in Multiple Sclerosis

BACKGROUND

Multiple sclerosis (MS) is an immune-mediated neurological disorder, and up to 50% of patients experience depression. We investigated how white matter network disruption is related to depression in MS.

METHODS

Using electronic health records, 380 participants with MS were identified. Depressed individuals (MS+Depression group; n = 232) included persons who had an ICD-10 depression diagnosis, had a prescription for antidepressant medication, or screened positive via Patient Health Questionnaire (PHQ)-2 or PHQ-9. Age- and sex-matched nondepressed individuals with MS (MS-Depression group; n = 148) included persons who had no prior depression diagnosis, had no psychiatric medication prescriptions, and were asymptomatic on PHQ-2 or PHQ-9. Research-quality 3T structural magnetic resonance imaging was obtained as part of routine care. We first evaluated whether lesions were preferentially located within the depression network compared with other brain regions. Next, we examined if MS+Depression patients had greater lesion burden and if this was driven by lesions in the depression network. Primary outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain.

RESULTS

MS lesions preferentially affected fascicles within versus outside the depression network (β = 0.09, 95% CI = 0.08 to 0.10, p < .001). MS+Depression patients had more lesion burden (β = 0.06, 95% CI = 0.01 to 0.10, p = .015); this was driven by lesions within the depression network (β = 0.02, 95% CI = 0.003 to 0.040, p = .020).

CONCLUSIONS

We demonstrated that lesion location and burden may contribute to depression comorbidity in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression patients had more disease than MS-Depression patients, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

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 the spatial structure 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 genomics 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 phenotype associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Intersite brain MRI volumetric biases persist even in a harmonized multisubject study of multiple sclerosis

BACKGROUND AND PURPOSE

Multicenter study designs involving a variety of MRI scanners have become increasingly common. However, these present the issue of biases in image-based measures due to scanner or site differences. To assess these biases, we imaged 11 volunteers with multiple sclerosis (MS) with scan and rescan data at four sites.

METHODS

Images were acquired on Siemens or Philips scanners at 3 Tesla. Automated white matter lesion detection and whole-brain, gray and white matter, and thalamic volumetry were performed, as well as expert manual delineations of T1 magnetization-prepared rapid acquisition gradient echo and T2 fluid-attenuated inversion recovery lesions. Random-effect and permutation-based nonparametric modeling was performed to assess differences in estimated volumes within and across sites.

RESULTS

Random-effect modeling demonstrated model assumption violations for most comparisons of interest. Nonparametric modeling indicated that site explained >50% of the variation for most estimated volumes. This expanded to >75% when data from both Siemens and Philips scanners were included. Permutation tests revealed significant differences between average inter- and intrasite differences in most estimated brain volumes (P < .05). The automatic activation of spine coil elements during some acquisitions resulted in a shading artifact in these images. Permutation tests revealed significant differences between thalamic volume measurements from acquisitions with and without this artifact.

CONCLUSION

Differences in brain volumetry persisted across MR scanners despite protocol harmonization. These differences were not well explained by variance component modeling; however, statistical innovations for mitigating intersite differences show promise in reducing biases in multicenter studies of MS.

The overlapping genetic architecture of psychiatric disorders and cortical brain structure

Both psychiatric vulnerability and cortical structure are shaped by the cumulative effect of common genetic variants across the genome. However, the shared genetic underpinnings between psychiatric disorders and brain structural phenotypes, such as thickness and surface area of the cerebral cortex, remains elusive. In this study, we employed pleiotropy-informed conjunctional false discovery rate analysis to investigate shared loci across genome-wide association scans of regional cortical thickness, surface area, and seven psychiatric disorders in approximately 700,000 individuals of European ancestry. Aggregating regional measures, we identified 50 genetic loci shared between psychiatric disorders and surface area, as well as 26 genetic loci shared with cortical thickness. Risk alleles exhibited bidirectional effects on both cortical thickness and surface area, such that some risk alleles for each disorder increased regional brain size while other risk alleles decreased regional brain size. Due to bidirectional effects, in many cases we observed extensive pleiotropy between an imaging phenotype and a psychiatric disorder even in the absence of a significant genetic correlation between them. The impact of genetic risk for psychiatric disorders on regional brain structure did exhibit a consistent pattern across highly comorbid psychiatric disorders, with 80% of the genetic loci shared across multiple disorders displaying consistent directions of effect. Cortical patterning of genetic overlap revealed a hierarchical genetic architecture, with the association cortex and sensorimotor cortex representing two extremes of shared genetic influence on psychiatric disorders and brain structural variation. Integrating multi-scale functional annotations and transcriptomic profiles, we observed that shared genetic loci were enriched in active genomic regions, converged on neurobiological and metabolic pathways, and showed differential expression in postmortem brain tissue from individuals with psychiatric disorders. Cumulatively, these findings provide a significant advance in our understanding of the overlapping polygenic architecture between psychopathology and cortical brain structure.

Brain Growth Charts for Quantitative Analysis of Pediatric Clinical Brain MRI Scans with Limited Imaging Pathology

Background Clinically acquired brain MRI scans represent a valuable but underused resource for investigating neurodevelopment due to their technical heterogeneity and lack of appropriate controls. These barriers have curtailed retrospective studies of clinical brain MRI scans compared with more costly prospectively acquired research-quality brain MRI scans. Purpose To provide a benchmark for neuroanatomic variability in clinically acquired brain MRI scans with limited imaging pathology (SLIPs) and to evaluate if growth charts from curated clinical MRI scans differed from research-quality MRI scans or were influenced by clinical indication for the scan. Materials and Methods In this secondary analysis of preexisting data, clinical brain MRI SLIPs from an urban pediatric health care system (individuals aged ≤22 years) were scanned across nine 3.0-T MRI scanners. The curation process included manual review of signed radiology reports and automated and manual quality review of images without gross pathology. Global and regional volumetric imaging phenotypes were measured using two image segmentation pipelines, and clinical brain growth charts were quantitatively compared with charts derived from a large set of research controls in the same age range by means of Pearson correlation and age at peak volume. Results The curated clinical data set included 532 patients (277 male; median age, 10 years [IQR, 5-14 years]; age range, 28 days after birth to 22 years) scanned between 2005 and 2020. Clinical brain growth charts were highly correlated with growth charts derived from research data sets (22 studies, 8346 individuals [4947 male]; age range, 152 days after birth to 22 years) in terms of normative developmental trajectories predicted by the models (median r = 0.979). Conclusion The clinical indication of the scans did not significantly bias the output of clinical brain charts. Brain growth charts derived from clinical controls with limited imaging pathology were highly correlated with brain charts from research controls, suggesting the potential of curated clinical MRI scans to supplement research data sets. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Ertl-Wagner and Pai in this issue.

Early Magnetic Resonance Imaging Features of New Paramagnetic Rim Lesions in Multiple Sclerosis

OBJECTIVE

To determine early magnetic resonance imaging (MRI) features of new multiple sclerosis (MS) lesions that will develop into paramagnetic rim lesions (PRLs), which have been associated with progressive tissue injury in MS.

METHODS

New contrast-enhancing lesions observed on routine clinical MRI were imaged at 7 T within 4 weeks of observation, and 3 and 6 months later. The 6-month MRI was used to classify PRL status (PRL or non-PRL). The relationship between early lesion characteristics and subsequent PRL status was assessed using generalized linear mixed effects models. Random forest classification was performed to classify early predictors of subsequent PRL status.

RESULTS

From 93 contrast-enhancing lesions in 23 MS patients, 37 lesions developed into a PRL. In lesions that developed into PRLs compared with those that did not, the average lesion T1 on the initial 7 T MRI was 1994 ms compared with 1,670 ms (p-value <0.001), and the average volume was 168.7 mL compared with 44 mL (p-value <0.001) in lesions that did not. These volume differences were also found on 3 T scans (p-value <0.001), and for intensity-normalized T1 -w (p-value = 0.011) and fluid-attenuated inversion recovery (p-value = 0.005). The area under the receiver operating characteristic curve for the random forest classification with leave-one-out cross-validation was found to be 0.86 using initial 7 T features.

INTERPRETATION

New MS lesions that evolve into PRLs can be identified early in lesion evolution. These findings suggest that biological mechanisms underlying PRL development begin early, which has important implications for clinical trials targeting PRLs development and subsequent therapeutics. ANN NEUROL 2023;94:736-744.

Performance of an automated deep learning algorithm to identify hepatic steatosis within noncontrast computed tomography scans among people with and without HIV

PURPOSE

Hepatic steatosis (fatty liver disease) affects 25% of the world's population, particularly people with HIV (PWH). Pharmacoepidemiologic studies to identify medications associated with steatosis have not been conducted because methods to evaluate liver fat within digitized images have not been developed. We determined the accuracy of a deep learning algorithm (automatic liver attenuation region-of-interest-based measurement [ALARM]) to identify steatosis within clinically obtained noncontrast abdominal CT images compared to manual radiologist review and evaluated its performance by HIV status.

METHODS

We performed a cross-sectional study to evaluate the performance of ALARM within noncontrast abdominal CT images from a sample of patients with and without HIV in the US Veterans Health Administration. We evaluated the ability of ALARM to identify moderate-to-severe hepatic steatosis, defined by mean absolute liver attenuation <40 Hounsfield units (HU), compared to manual radiologist assessment.

RESULTS

Among 120 patients (51 PWH) who underwent noncontrast abdominal CT, moderate-to-severe hepatic steatosis was identified in 15 (12.5%) persons via ALARM and 12 (10%) by radiologist assessment. Percent agreement between ALARM and radiologist assessment of absolute liver attenuation <40 HU was 95.8%. Sensitivity, specificity, positive predictive value, and negative predictive value of ALARM were 91.7% (95%CI, 51.5%-99.8%), 96.3% (95%CI, 90.8%-99.0%), 73.3% (95%CI, 44.9%-92.2%), and 99.0% (95%CI, 94.8%-100%), respectively. No differences in performance were observed by HIV status.

CONCLUSIONS

ALARM demonstrated excellent accuracy for moderate-to-severe hepatic steatosis regardless of HIV status. Application of ALARM to radiographic repositories could facilitate real-world studies to evaluate medications associated with steatosis and assess differences by HIV status.

Myelination and excitation-inhibition balance synergistically shape structure-function coupling across the human cortex

Recent work has demonstrated that the relationship between structural and functional connectivity varies regionally across the human brain, with reduced coupling emerging along the sensory-association cortical hierarchy. The biological underpinnings driving this expression, however, remain largely unknown. Here, we postulate that intracortical myelination and excitation-inhibition (EI) balance mediate the heterogeneous expression of structure-function coupling (SFC) and its temporal variance across the cortical hierarchy. We employ atlas- and voxel-based connectivity approaches to analyze neuroimaging data acquired from two groups of healthy participants. Our findings are consistent across six complementary processing pipelines: 1) SFC and its temporal variance respectively decrease and increase across the unimodal-transmodal and granular-agranular gradients; 2) increased myelination and lower EI-ratio are associated with more rigid SFC and restricted moment-to-moment SFC fluctuations; 3) a gradual shift from EI-ratio to myelination as the principal predictor of SFC occurs when traversing from granular to agranular cortical regions. Collectively, our work delivers a framework to conceptualize structure-function relationships in the human brain, paving the way for an improved understanding of how demyelination and/or EI-imbalances induce reorganization in brain disorders.

Automated analysis of low-field brain MRI in cerebral malaria

A central challenge of medical imaging studies is to extract biomarkers that characterize disease pathology or outcomes. Modern automated approaches have found tremendous success in high-resolution, high-quality magnetic resonance images. These methods, however, may not translate to low-resolution images acquired on magnetic resonance imaging (MRI) scanners with lower magnetic field strength. In low-resource settings where low-field scanners are more common and there is a shortage of radiologists to manually interpret MRI scans, it is critical to develop automated methods that can augment or replace manual interpretation, while accommodating reduced image quality. We present a fully automated framework for translating radiological diagnostic criteria into image-based biomarkers, inspired by a project in which children with cerebral malaria (CM) were imaged using low-field 0.35 Tesla MRI. We integrate multiatlas label fusion, which leverages high-resolution images from another sample as prior spatial information, with parametric Gaussian hidden Markov models based on image intensities, to create a robust method for determining ventricular cerebrospinal fluid volume. We also propose normalized image intensity and texture measurements to determine the loss of gray-to-white matter tissue differentiation and sulcal effacement. These integrated biomarkers have excellent classification performance for determining severe brain swelling due to CM.

Characterizing Subcortical Structural Heterogeneity in Autism

Autism presents with significant phenotypic and neuroanatomical heterogeneity, and neuroimaging studies of the thalamus, globus pallidus and striatum in autism have produced inconsistent and contradictory results. These structures are critical mediators of functions known to be atypical in autism, including sensory gating and motor function. We examined both volumetric and fine-grained localized shape differences in autism using a large (n=3145, 1045-1318 after strict quality control), cross-sectional dataset of T1-weighted structural MRI scans from 32 sites, including both males and females (assigned-at-birth). We investigated three potentially important sources of neuroanatomical heterogeneity: sex, age, and intelligence quotient (IQ), using a meta-analytic technique after strict quality control to minimize non-biological sources of variation. We observed no volumetric differences in the thalamus, globus pallidus, or striatum in autism. Rather, we identified a variety of localized shape differences in all three structures. Including age, but not sex or IQ, in the statistical model improved the fit for both the pallidum and striatum, but not for the thalamus. Age-centered shape analysis indicated a variety of age-dependent regional differences. Overall, our findings help confirm that the neurodevelopment of the striatum, globus pallidus and thalamus are atypical in autism, in a subtle location-dependent manner that is not reflected in overall structure volumes, and that is highly non-uniform across the lifespan.

Using network control theory to study the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. 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 dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Generalizable links between symptoms of borderline personality disorder and functional connectivity

Background |

Symptoms of borderline personality disorder (BPD) often manifest in adolescence, yet the underlying relationship between these debilitating symptoms and the development of functional brain networks is not well understood. Here we aimed to investigate how multivariate patterns of functional connectivity are associated with symptoms of BPD in a large sample of young adults and adolescents.

Methods |

We used high-quality functional Magnetic Resonance Imaging (fMRI) data from young adults from the Human Connectome Project: Young Adults (HCP-YA; N = 870, ages 22-37 years, 457 female) and youth from the Human Connectome Project: Development (HCP-D; N = 223, age range 16-21 years, 121 female). A previously validated BPD proxy score was derived from the NEO Five Factor Inventory (NEO-FFI). A ridge regression model with 10-fold cross-validation and nested hyperparameter tuning was trained and tested in HCP-YA to predict BPD scores in unseen data from regional functional connectivity, while controlling for in-scanner motion, age, and sex. The trained model was further tested on data from HCP-D without further tuning. Finally, we tested how the connectivity patterns associated with BPD aligned with age-related changes in connectivity.

Results |

Multivariate functional connectivity patterns significantly predicted out-of-sample BPD proxy scores in unseen data in both young adults (HCP-YA; pperm = 0.001) and older adolescents (HCP-D; pperm = 0.001). Predictive capacity of regions was heterogeneous; the most predictive regions were found in functional systems relevant for emotion regulation and executive function, including the ventral attention network. Finally, regional functional connectivity patterns that predicted BPD proxy scores aligned with those associated with development in youth.

Conclusion |

Individual differences in functional connectivity in developmentally-sensitive regions are associated with the symptoms of BPD.

Penalized decomposition using residuals (PeDecURe) for feature extraction in the presence of nuisance variables

Neuroimaging data are an increasingly important part of etiological studies of neurological and psychiatric disorders. However, mitigating the influence of nuisance variables, including confounders, remains a challenge in image analysis. In studies of Alzheimer's disease, for example, an imbalance in disease rates by age and sex may make it difficult to distinguish between structural patterns in the brain (as measured by neuroimaging scans) attributable to disease progression and those characteristic of typical human aging or sex differences. Concerningly, when not properly accounted for, nuisance variables pose threats to the generalizability and interpretability of findings from these studies. Motivated by this critical issue, in this work, we examine the impact of nuisance variables on feature extraction methods and propose Penalized Decomposition Using Residuals (PeDecURe), a new method for obtaining nuisance variable-adjusted features. PeDecURe estimates primary directions of variation which maximize covariance between partially residualized imaging features and a variable of interest (e.g., Alzheimer's diagnosis) while simultaneously mitigating the influence of nuisance variation through a penalty on the covariance between partially residualized imaging features and those variables. Using features derived using PeDecURe's first direction of variation, we train a highly accurate and generalizable predictive model, as evidenced by its robustness in testing samples with different underlying nuisance variable distributions. We compare PeDecURe to commonly used decomposition methods (principal component analysis (PCA) and partial least squares) as well as a confounder-adjusted variation of PCA. We find that features derived from PeDecURe offer greater accuracy and generalizability and lower correlations with nuisance variables compared with the other methods. While PeDecURe is primarily motivated by challenges that arise in the analysis of neuroimaging data, it is broadly applicable to data sets with highly correlated features, where novel methods to handle nuisance variables are warranted.

Long-term epilepsy outcome dynamics revealed by natural language processing of clinic notes

OBJECTIVE

Electronic medical records allow for retrospective clinical research with large patient cohorts. However, epilepsy outcomes are often contained in free text notes that are difficult to mine. We recently developed and validated novel natural language processing (NLP) algorithms to automatically extract key epilepsy outcome measures from clinic notes. In this study, we assessed the feasibility of extracting these measures to study the natural history of epilepsy at our center.

METHODS

We applied our previously validated NLP algorithms to extract seizure freedom, seizure frequency, and date of most recent seizure from outpatient visits at our epilepsy center from 2010 to 2022. We examined the dynamics of seizure outcomes over time using Markov model-based probability and Kaplan-Meier analyses.

RESULTS

Performance of our algorithms on classifying seizure freedom was comparable to that of human reviewers (algorithm F1 = .88 vs. human annotatorκ = .86). We extracted seizure outcome data from 55 630 clinic notes from 9510 unique patients written by 53 unique authors. Of these, 30% were classified as seizure-free since the last visit, 48% of non-seizure-free visits contained a quantifiable seizure frequency, and 47% of all visits contained the date of most recent seizure occurrence. Among patients with at least five visits, the probabilities of seizure freedom at the next visit ranged from 12% to 80% in patients having seizures or seizure-free at the prior three visits, respectively. Only 25% of patients who were seizure-free for 6 months remained seizure-free after 10 years.

SIGNIFICANCE

Our findings demonstrate that epilepsy outcome measures can be extracted accurately from unstructured clinical note text using NLP. At our tertiary center, the disease course often followed a remitting and relapsing pattern. This method represents a powerful new tool for clinical research with many potential uses and extensions to other clinical questions.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

Mapping the relationship of white matter lesions to depression in multiple sclerosis

Importance

Multiple sclerosis (MS) is an immune-mediated neurological disorder that affects nearly one million people in the United States. Up to 50% of patients with MS experience depression.

Objective

To investigate how white matter network disruption is related to depression in MS.

Design

Retrospective case-control study of participants who received research-quality 3-tesla neuroimaging as part of MS clinical care from 2010-2018. Analyses were performed from May 1 to September 30, 2022.

Setting

Single-center academic medical specialty MS clinic.

Participants

Participants with MS were identified via the electronic health record (EHR). All participants were diagnosed by an MS specialist and completed research-quality MRI at 3T. After excluding participants with poor image quality, 783 were included. Inclusion in the depression group (MS+Depression) required either: 1) ICD-10 depression diagnosis (F32-F34.*); 2) prescription of antidepressant medication; or 3) screening positive via Patient Health Questionnaire-2 (PHQ-2) or -9 (PHQ-9). Age- and sex-matched nondepressed comparators (MS-Depression) included persons with no depression diagnosis, no psychiatric medications, and were asymptomatic on PHQ-2/9.

Exposure

Depression diagnosis.

Main Outcomes and Measures

We first evaluated if lesions were preferentially located within the depression network compared to other brain regions. Next, we examined if MS+Depression patients had greater lesion burden, and if this was driven by lesions specifically in the depression network. Outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain. Secondary measures included between-diagnosis lesion burden, stratified by brain network. Linear mixed-effects models were employed.

Results

Three hundred-eighty participants met inclusion criteria, (232 MS+Depression: age[SD]=49[12], %females=86; 148 MS-Depression: age[SD]=47[13], %females=79). MS lesions preferentially affected fascicles within versus outside the depression network (β=0.09, 95% CI=0.08-0.10, P<0.001). MS+Depression had more white matter lesion burden (β=0.06, 95% CI=0.01-0.10, P=0.015); this was driven by lesions within the depression network (β=0.02, 95% CI 0.003-0.040, P=0.020).

Conclusions and Relevance

We provide new evidence supporting a relationship between white matter lesions and depression in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression had more disease than MS-Depression, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

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.

Resting state functional connectivity demonstrates increased segregation in bilateral temporal lobe epilepsy

OBJECTIVE

Temporal lobe epilepsy (TLE) is the most common type of focal epilepsy. An increasingly identified subset of patients with TLE consists of those who show bilaterally independent temporal lobe seizures. The purpose of this study was to leverage network neuroscience to better understand the interictal whole brain network of bilateral TLE (BiTLE).

METHODS

In this study, using a multicenter resting state functional magnetic resonance imaging (rs-fMRI) data set, we constructed whole-brain functional networks of 19 patients with BiTLE, and compared them to those of 75 patients with unilateral TLE (UTLE). We quantified resting-state, whole-brain topological properties using metrics derived from network theory, including clustering coefficient, global efficiency, participation coefficient, and modularity. For each metric, we computed an average across all brain regions, and iterated this process across network densities. Curves of network density vs each network metric were compared between groups. Finally, we derived a combined metric, which we term the "integration-segregation axis," by combining whole-brain average clustering coefficient and global efficiency curves, and applying principal component analysis (PCA)-based dimensionality reduction.

RESULTS

Compared to UTLE, BiTLE had decreased global efficiency (p = .031), and decreased whole brain average participation coefficient across a range of network densities (p = .019). Modularity maximization yielded a larger number of smaller communities in BiTLE than in UTLE (p = .020). Differences in network properties separate BiTLE and UTLE along the integration-segregation axis, with regions within the axis having a specificity of up to 0.87 for BiTLE. Along the integration-segregation axis, UTLE patients with poor surgical outcomes were distributed in the same regions as BiTLE, and network metrics confirmed similar patterns of increased segregation in both BiTLE and poor outcome UTLE.

SIGNIFICANCE

Increased interictal whole-brain network segregation, as measured by rs-fMRI, is specific to BiTLE, as well as poor surgical outcome UTLE, and may assist in non-invasively identifying this patient population prior to intracranial electroencephalography or device implantation.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

DeepComBat: A Statistically Motivated, Hyperparameter-Robust, Deep Learning Approach to Harmonization of Neuroimaging Data

Neuroimaging data from multiple batches (i.e. acquisition sites, scanner manufacturer, datasets, etc.) are increasingly necessary to gain new insights into the human brain. However, multi-batch data, as well as extracted radiomic features, exhibit pronounced technical artifacts across batches. These batch effects introduce confounding into the data and can obscure biological effects of interest, decreasing the generalizability and reproducibility of findings. This is especially true when multi-batch data is used alongside complex downstream analysis models, such as machine learning methods. Image harmonization methods seeking to remove these batch effects are important for mitigating these issues; however, significant multivariate batch effects remain in the data following harmonization by current state-of-the-art statistical and deep learning methods. We present DeepCombat, a deep learning harmonization method based on a conditional variational autoencoder architecture and the ComBat harmonization model. DeepCombat learns and removes subject-level batch effects by accounting for the multivariate relationships between features. Additionally, DeepComBat relaxes a number of strong assumptions commonly made by previous deep learning harmonization methods and is empirically robust across a wide range of hyperparameter choices. We apply this method to neuroimaging data from a large cognitive-aging cohort and find that DeepCombat outperforms existing methods, as assessed by a battery of machine learning methods, in removing scanner effects from cortical thickness measurements while preserving biological heterogeneity. Additionally, DeepComBat provides a new perspective for statistically-motivated deep learning harmonization methods.

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.

Assessment of Neuroanatomical Endophenotypes of Autism Spectrum Disorder and Association With Characteristics of Individuals With Schizophrenia and the General Population

Importance

Autism spectrum disorder (ASD) is associated with significant clinical, neuroanatomical, and genetic heterogeneity that limits precision diagnostics and treatment.

Objective

To assess distinct neuroanatomical dimensions of ASD using novel semisupervised machine learning methods and to test whether the dimensions can serve as endophenotypes also in non-ASD populations.

Design, Setting, and Participants

This cross-sectional study used imaging data from the publicly available Autism Brain Imaging Data Exchange (ABIDE) repositories as the discovery cohort. The ABIDE sample included individuals diagnosed with ASD aged between 16 and 64 years and age- and sex-match typically developing individuals. Validation cohorts included individuals with schizophrenia from the Psychosis Heterogeneity Evaluated via Dimensional Neuroimaging (PHENOM) consortium and individuals from the UK Biobank to represent the general population. The multisite discovery cohort included 16 internationally distributed imaging sites. Analyses were performed between March 2021 and March 2022.

Main Outcomes and Measures

The trained semisupervised heterogeneity through discriminative analysis models were tested for reproducibility using extensive cross-validations. It was then applied to individuals from the PHENOM and the UK Biobank. It was hypothesized that neuroanatomical dimensions of ASD would display distinct clinical and genetic profiles and would be prominent also in non-ASD populations.

Results

Heterogeneity through discriminative analysis models trained on T1-weighted brain magnetic resonance images of 307 individuals with ASD (mean [SD] age, 25.4 [9.8] years; 273 [88.9%] male) and 362 typically developing control individuals (mean [SD] age, 25.8 [8.9] years; 309 [85.4%] male) revealed that a 3-dimensional scheme was optimal to capture the ASD neuroanatomy. The first dimension (A1: aginglike) was associated with smaller brain volume, lower cognitive function, and aging-related genetic variants (FOXO3; Z = 4.65; P = 1.62 × 10-6). The second dimension (A2: schizophrenialike) was characterized by enlarged subcortical volumes, antipsychotic medication use (Cohen d = 0.65; false discovery rate-adjusted P = .048), partially overlapping genetic, neuroanatomical characteristics to schizophrenia (n = 307), and significant genetic heritability estimates in the general population (n = 14 786; mean [SD] h2, 0.71 [0.04]; P < 1 × 10-4). The third dimension (A3: typical ASD) was distinguished by enlarged cortical volumes, high nonverbal cognitive performance, and biological pathways implicating brain development and abnormal apoptosis (mean [SD] β, 0.83 [0.02]; P = 4.22 × 10-6).

Conclusions and Relevance

This cross-sectional study discovered 3-dimensional endophenotypic representation that may elucidate the heterogeneous neurobiological underpinnings of ASD to support precision diagnostics. The significant correspondence between A2 and schizophrenia indicates a possibility of identifying common biological mechanisms across the 2 mental health diagnoses.

Intrinsic activity development unfolds along a sensorimotor-association cortical axis in youth

Animal studies of neurodevelopment have shown that recordings of intrinsic cortical activity evolve from synchronized and high amplitude to sparse and low amplitude as plasticity declines and the cortex matures. Leveraging resting-state functional MRI (fMRI) data from 1,033 youths (ages 8-23 years), we find that this stereotyped refinement of intrinsic activity occurs during human development and provides evidence for a cortical gradient of neurodevelopmental change. Declines in the amplitude of intrinsic fMRI activity were initiated heterochronously across regions and were coupled to the maturation of intracortical myelin, a developmental plasticity regulator. Spatiotemporal variability in regional developmental trajectories was organized along a hierarchical, sensorimotor-association cortical axis from ages 8 to 18. The sensorimotor-association axis furthermore captured variation in associations between youths' neighborhood environments and intrinsic fMRI activity; associations suggest that the effects of environmental disadvantage on the maturing brain diverge most across this axis during midadolescence. These results uncover a hierarchical neurodevelopmental axis and offer insight into the progression of cortical plasticity in humans.

The etiology and evolution of magnetic resonance imaging-visible perivascular spaces: Systematic review and meta-analysis

Objectives

Perivascular spaces have been involved in neuroinflammatory and neurodegenerative diseases. Upon a certain size, these spaces can become visible on magnetic resonance imaging (MRI), referred to as enlarged perivascular spaces (EPVS) or MRI-visible perivascular spaces (MVPVS). However, the lack of systematic evidence on etiology and temporal dynamics of MVPVS hampers their diagnostic utility as MRI biomarker. Thus, the goal of this systematic review was to summarize potential etiologies and evolution of MVPVS.

Methods

In a comprehensive literature search, out of 1,488 unique publications, 140 records assessing etiopathogenesis and dynamics of MVPVS were eligible for a qualitative summary. 6 records were included in a meta-analysis to assess the association between MVPVS and brain atrophy.

Results

Four overarching and partly overlapping etiologies of MVPVS have been proposed: (1) Impairment of interstitial fluid circulation, (2) Spiral elongation of arteries, (3) Brain atrophy and/or perivascular myelin loss, and (4) Immune cell accumulation in the perivascular space. The meta-analysis in patients with neuroinflammatory diseases did not support an association between MVPVS and brain volume measures [R: −0.15 (95%-CI −0.40–0.11)]. Based on few and mostly small studies in tumefactive MVPVS and in vascular and neuroinflammatory diseases, temporal evolution of MVPVS is slow.

Conclusion

Collectively, this study provides high-grade evidence for MVPVS etiopathogenesis and temporal dynamics. Although several potential etiologies for MVPVS emergence have been proposed, they are only partially supported by data. Advanced MRI methods should be employed to further dissect etiopathogenesis and evolution of MVPVS. This can benefit their implementation as an imaging biomarker.

Systematic review registration

https://www.crd.york.ac.uk/prospero/display_record.php?RecordID=346564, identifier CRD42022346564.

Three-dimensional printing of patient-specific computed tomography lung phantoms: a reader study

In modern clinical decision-support algorithms, heterogeneity in image characteristics due to variations in imaging systems and protocols hinders the development of reproducible quantitative measures including for feature extraction pipelines. With the help of a reader study, we investigate the ability to provide consistent ground-truth targets by using patient-specific 3D-printed lung phantoms. PixelPrint was developed for 3D-printing lifelike computed tomography (CT) lung phantoms by directly translating clinical images into printer instructions that control density on a voxel-by-voxel basis. Data sets of three COVID-19 patients served as input for 3D-printing lung phantoms. Five radiologists rated patient and phantom images for imaging characteristics and diagnostic confidence in a blinded reader study. Effect sizes of evaluating phantom as opposed to patient images were assessed using linear mixed models. Finally, PixelPrint's production reproducibility was evaluated. Images of patients and phantoms had little variation in the estimated mean (0.03-0.29, using a 1-5 scale). When comparing phantom images to patient images, effect size analysis revealed that the difference was within one-third of the inter- and intrareader variabilities. High correspondence between the four phantoms created using the same patient images was demonstrated by PixelPrint's production repeatability tests, with greater similarity scores between high-dose acquisitions of the phantoms than between clinical-dose acquisitions of a single phantom. We demonstrated PixelPrint's ability to produce lifelike CT lung phantoms reliably. These phantoms have the potential to provide ground-truth targets for validating the generalizability of inference-based decision-support algorithms between different health centers and imaging protocols and for optimizing examination protocols with realistic patient-based phantoms. Classification: CT lung phantoms, reader study.