Microstructural Mapping of Neural Pathways in Alzheimer's Disease using Macrostructure-Informed Normative Tractometry.

Introduction: Diffusion MRI is sensitive to the microstructural properties of brain tissues, and shows great promise in detecting the effects of degenerative diseases. However, many approaches analyze single measures averaged over regions of interest, without considering the underlying fiber geometry. Methods: Here, we propose a novel Macrostructure-Informed Normative Tractometry (MINT) framework, to investigate how white matter microstructure and macrostructure are jointly altered in mild cognitive impairment (MCI) and dementia. We compare MINT-derived metrics with univariate metrics from diffusion tensor imaging (DTI), to examine how fiber geometry may impact interpretation of microstructure. Results: In two multi-site cohorts from North America and India, we find consistent patterns of microstructural and macrostructural anomalies implicated in MCI and dementia; we also rank diffusion metrics' sensitivity to dementia. Discussion: We show that MINT, by jointly modeling tract shape and microstructure, has potential to disentangle and better interpret the effects of degenerative disease on the brain's neural pathways.

Harmonizing three-dimensional MRI using pseudo-warping field guided GAN.

In pursuit of cultivating automated models for magnetic resonance imaging (MRI) to aid in diagnostics, an escalating demand for extensive, multisite, and heterogeneous brain imaging datasets has emerged. This potentially introduces biased outcomes when directly applied for subsequent analysis. Researchers have endeavored to address this issue by pursuing the harmonization of MRIs. However, most existing image-based harmonization methods for MRI are tailored for 2D slices, which may introduce inter-slice variations when they are combined into a 3D volume. In this study, we aim to resolve inconsistencies between slices by introducing a pseudo-warping field. This field is created randomly and utilized to transform a slice into an artificially warped subsequent slice. The objective of this pseudo-warping field is to ensure that generators can consistently harmonize adjacent slices to another domain, without being affected by the varying content present in different slices. Furthermore, we construct unsupervised spatial and recycle loss to enhance the spatial accuracy and slice-wise consistency across the 3D images. The results demonstrate that our model effectively mitigates inter-slice variations and successfully preserves the anatomical details of the images during the harmonization process. Compared to generative harmonization models that employ 3D operators, our model exhibits greater computational efficiency and flexibility.

Subclinical variations on ECG and their associations with structural brain aging networks.

Impaired cardiac function is associated with cognitive impairment and brain imaging features of aging. Cardiac arrhythmias, including atrial fibrillation, are implicated in clinical and subclinical brain injuries. Even in the absence of a clinical diagnosis, subclinical or prodromal substrates of arrhythmias, including an abnormally long or short P-wave duration (PWD), a measure associated with atrial abnormalities, have been associated with stroke and cognitive decline. However, the extent to which PWD has subclinical influences on overall aging patterns of the brain is not clearly understood. Here, using neuroimaging and ECG data from the UK Biobank, we use a novel regional "brain age" method to identify the brain aging networks associated with abnormal PWD. We find that PWD is inversely associated with accelerated brain aging in the sensorimotor, frontoparietal, ventral attention, and dorsal attention networks, even in the absence of overt cardiac diseases. These findings suggest that detrimental aging outcomes may result from subclinically abnormal PWD.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI.

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

White matter diffusion estimates in obsessive-compulsive disorder across 1653 individuals: machine learning findings from the ENIGMA OCD Working Group.

White matter pathways, typically studied with diffusion tensor imaging (DTI), have been implicated in the neurobiology of obsessive-compulsive disorder (OCD). However, due to limited sample sizes and the predominance of single-site studies, the generalizability of OCD classification based on diffusion white matter estimates remains unclear. Here, we tested classification accuracy using the largest OCD DTI dataset to date, involving 1336 adult participants (690 OCD patients and 646 healthy controls) and 317 pediatric participants (175 OCD patients and 142 healthy controls) from 18 international sites within the ENIGMA OCD Working Group. We used an automatic machine learning pipeline (with feature engineering and selection, and model optimization) and examined the cross-site generalizability of the OCD classification models using leave-one-site-out cross-validation. Our models showed low-to-moderate accuracy in classifying (1) "OCD vs. healthy controls" (Adults, receiver operator characteristic-area under the curve = 57.19 ± 3.47 in the replication set; Children, 59.8 ± 7.39), (2) "unmedicated OCD vs. healthy controls" (Adults, 62.67 ± 3.84; Children, 48.51 ± 10.14), and (3) "medicated OCD vs. unmedicated OCD" (Adults, 76.72 ± 3.97; Children, 72.45 ± 8.87). There was significant site variability in model performance (cross-validated ROC AUC ranges 51.6-79.1 in adults; 35.9-63.2 in children). Machine learning interpretation showed that diffusivity measures of the corpus callosum, internal capsule, and posterior thalamic radiation contributed to the classification of OCD from HC. The classification performance appeared greater than the model trained on grey matter morphometry in the prior ENIGMA OCD study (our study includes subsamples from the morphometry study). Taken together, this study points to the meaningful multivariate patterns of white matter features relevant to the neurobiology of OCD, but with low-to-moderate classification accuracy. The OCD classification performance may be constrained by site variability and medication effects on the white matter integrity, indicating room for improvement for future research.

Normative modelling of brain morphometry across the lifespan with CentileBrain: algorithm benchmarking and model optimisation.

The value of normative models in research and clinical practice relies on their robustness and a systematic comparison of different modelling algorithms and parameters; however, this has not been done to date. We aimed to identify the optimal approach for normative modelling of brain morphometric data through systematic empirical benchmarking, by quantifying the accuracy of different algorithms and identifying parameters that optimised model performance. We developed this framework with regional morphometric data from 37 407 healthy individuals (53% female and 47% male; aged 3-90 years) from 87 datasets from Europe, Australia, the USA, South Africa, and east Asia following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The multivariate fractional polynomial regression (MFPR) emerged as the preferred algorithm, optimised with non-linear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3000 study participants. This model can inform about the biological and behavioural implications of deviations from typical age-related neuroanatomical changes and support future study designs. The model and scripts described here are freely available through CentileBrain.

Deep Normative Tractometry for Identifying Joint White Matter Macro- and Micro-structural Abnormalities in Alzheimer's Disease.

This study introduces the Deep Normative Tractometry (DNT) framework, that encodes the joint distribution of both macrostructural and microstructural profiles of the brain white matter tracts through a variational autoencoder (VAE). By training on data from healthy controls, DNT learns the normative distribution of tract data, and can delineate along-tract micro-and macro-structural abnormalities. Leveraging a large sample size via generative pre-training, we assess DNT's generalizability using transfer learning on data from an independent cohort acquired in India. Our findings demonstrate DNT's capacity to detect widespread diffusivity abnormalities along tracts in mild cognitive impairment and Alzheimer's disease, aligning closely with results from the Bundle Analytics (BUAN) tractometry pipeline. By incorporating tract geometry information, DNT may be able to distinguish disease-related abnormalities in anisotropy from tract macrostructure, and shows promise in enhancing fine-scale mapping and detection of white matter alterations in neurodegenerative conditions.

Multi-site benchmark classification of major depressive disorder using machine learning on cortical and subcortical measures.

Machine learning (ML) techniques have gained popularity in the neuroimaging field due to their potential for classifying neuropsychiatric disorders. However, the diagnostic predictive power of the existing algorithms has been limited by small sample sizes, lack of representativeness, data leakage, and/or overfitting. Here, we overcome these limitations with the largest multi-site sample size to date (N = 5365) to provide a generalizable ML classification benchmark of major depressive disorder (MDD) using shallow linear and non-linear models. Leveraging brain measures from standardized ENIGMA analysis pipelines in FreeSurfer, we were able to classify MDD versus healthy controls (HC) with a balanced accuracy of around 62%. But after harmonizing the data, e.g., using ComBat, the balanced accuracy dropped to approximately 52%. Accuracy results close to random chance levels were also observed in stratified groups according to age of onset, antidepressant use, number of episodes and sex. Future studies incorporating higher dimensional brain imaging/phenotype features, and/or using more advanced machine and deep learning methods may yield more encouraging prospects.

Mapping gray and white matter volume abnormalities in early-onset psychosis: an ENIGMA multicenter voxel-based morphometry study.

INTRODUCTION: Regional gray matter (GM) alterations have been reported in early-onset psychosis (EOP, onset before age 18), but previous studies have yielded conflicting results, likely due to small sample sizes and the different brain regions examined. In this study, we conducted a whole brain voxel-based morphometry (VBM) analysis in a large sample of individuals with EOP, using the newly developed ENIGMA-VBM tool. METHODS: 15 independent cohorts from the ENIGMA-EOP working group participated in the study. The overall sample comprised T1-weighted MRI data from 482 individuals with EOP and 469 healthy controls. Each site performed the VBM analysis locally using the standardized ENIGMA-VBM tool. Statistical parametric T-maps were generated from each cohort and meta-analyzed to reveal voxel-wise differences between EOP and healthy controls as well as the individual-based association between GM volume and age of onset, chlorpromazine (CPZ) equivalent dose, and other clinical variables. RESULTS: Compared with healthy controls, individuals with EOP showed widespread lower GM volume encompassing most of the cortex, with the most marked effect in the left median cingulate (Hedges' g = 0.55, p = 0.001 corrected), as well as small clusters of lower white matter (WM), whereas no regional GM or WM volumes were higher in EOP. Lower GM volume in the cerebellum, thalamus and left inferior parietal gyrus was associated with older age of onset. Deficits in GM in the left inferior frontal gyrus, right insula, right precentral gyrus and right superior frontal gyrus were also associated with higher CPZ equivalent doses. CONCLUSION: EOP is associated with widespread reductions in cortical GM volume, while WM is affected to a smaller extent. GM volume alterations are associated with age of onset and CPZ equivalent dose but these effects are small compared to case-control differences. Mapping anatomical abnormalities in EOP may lead to a better understanding of the role of psychosis in brain development during childhood and adolescence.

Smaller total and subregional cerebellar volumes in posttraumatic stress disorder: a mega-analysis by the ENIGMA-PGC PTSD workgroup.

Although the cerebellum contributes to higher-order cognitive and emotional functions relevant to posttraumatic stress disorder (PTSD), prior research on cerebellar volume in PTSD is scant, particularly when considering subregions that differentially map on to motor, cognitive, and affective functions. In a sample of 4215 adults (PTSD n = 1642; Control n = 2573) across 40 sites from the ENIGMA-PGC PTSD working group, we employed a new state-of-the-art deep-learning based approach for automatic cerebellar parcellation to obtain volumetric estimates for the total cerebellum and 28 subregions. Linear mixed effects models controlling for age, gender, intracranial volume, and site were used to compare cerebellum volumes in PTSD compared to healthy controls (88% trauma-exposed). PTSD was associated with significant grey and white matter reductions of the cerebellum. Compared to controls, people with PTSD demonstrated smaller total cerebellum volume, as well as reduced volume in subregions primarily within the posterior lobe (lobule VIIB, crus II), vermis (VI, VIII), flocculonodular lobe (lobule X), and corpus medullare (all p-FDR < 0.05). Effects of PTSD on volume were consistent, and generally more robust, when examining symptom severity rather than diagnostic status. These findings implicate regionally specific cerebellar volumetric differences in the pathophysiology of PTSD. The cerebellum appears to play an important role in higher-order cognitive and emotional processes, far beyond its historical association with vestibulomotor function. Further examination of the cerebellum in trauma-related psychopathology will help to clarify how cerebellar structure and function may disrupt cognitive and affective processes at the center of translational models for PTSD.

Mapping Brain Structure Variability in Chronic Pain: The Role of Widespreadness and Pain Type and Its Mediating Relationship With Suicide Attempt.

BACKGROUND: Chronic pain affects nearly 20% of the U.S. POPULATION: It is a leading cause of disability globally and is associated with a heightened risk for suicide. The role of the central nervous system in the perception and maintenance of chronic pain has recently been accepted, but specific brain circuitries involved have yet to be mapped across pain types in a large-scale study. METHODS: We used data from the UK Biobank (N = 21,968) to investigate brain structural alterations in individuals reporting chronic pain compared with pain-free control participants and their mediating effect on history of suicide attempt. RESULTS: Chronic pain and, more notably, chronic multisite pain was associated with, on average, lower surface area throughout the cortex after adjusting for demographic, clinical, and neuropsychiatric confounds. Only participants with abdominal pain showed lower subcortical volumes, including the amygdala and brainstem, and lower cerebellum volumes. Participants with chronic headaches showed a widespread thicker cortex compared with control participants. Mediation analyses revealed that precuneus thickness mediated the relationship of chronic multisite pain and history of suicide attempt. Mediating effects were also identified specific to localized pain, with the strongest effect being amygdala volume in individuals with chronic abdominal pain. CONCLUSIONS: Results support a widespread effect of chronic pain on brain structure and distinct brain structures underlying chronic musculoskeletal pain, visceral pain, and headaches. Mediation effects of regions in the extended ventromedial prefrontal cortex subsystem suggest that exacerbated negative internal states, negative self-referencing, and impairments in future planning may underlie suicidal behaviors in individuals with chronic pain.

Normative Modeling of Brain Morphometry in Clinical High Risk for Psychosis.

Importance: The lack of robust neuroanatomical markers of psychosis risk has been traditionally attributed to heterogeneity. A complementary hypothesis is that variation in neuroanatomical measures in individuals at psychosis risk may be nested within the range observed in healthy individuals. Objective: To quantify deviations from the normative range of neuroanatomical variation in individuals at clinical high risk for psychosis (CHR-P) and evaluate their overlap with healthy variation and their association with positive symptoms, cognition, and conversion to a psychotic disorder. Design, Setting, and Participants: This case-control study used clinical-, IQ-, and neuroimaging software (FreeSurfer)-derived regional measures of cortical thickness (CT), cortical surface area (SA), and subcortical volume (SV) from 1340 individuals with CHR-P and 1237 healthy individuals pooled from 29 international sites participating in the Enhancing Neuroimaging Genetics Through Meta-analysis (ENIGMA) Clinical High Risk for Psychosis Working Group. Healthy individuals and individuals with CHR-P were matched on age and sex within each recruitment site. Data were analyzed between September 1, 2021, and November 30, 2022. Main Outcomes and Measures: For each regional morphometric measure, deviation scores were computed as z scores indexing the degree of deviation from their normative means from a healthy reference population. Average deviation scores (ADS) were also calculated for regional CT, SA, and SV measures and globally across all measures. Regression analyses quantified the association of deviation scores with clinical severity and cognition, and 2-proportion z tests identified case-control differences in the proportion of individuals with infranormal (z < -1.96) or supranormal (z > 1.96) scores. Results: Among 1340 individuals with CHR-P, 709 (52.91%) were male, and the mean (SD) age was 20.75 (4.74) years. Among 1237 healthy individuals, 684 (55.30%) were male, and the mean (SD) age was 22.32 (4.95) years. Individuals with CHR-P and healthy individuals overlapped in the distributions of the observed values, regional z scores, and all ADS values. For any given region, the proportion of individuals with CHR-P who had infranormal or supranormal values was low (up to 153 individuals [<11.42%]) and similar to that of healthy individuals (<115 individuals [<9.30%]). Individuals with CHR-P who converted to a psychotic disorder had a higher percentage of infranormal values in temporal regions compared with those who did not convert (7.01% vs 1.38%) and healthy individuals (5.10% vs 0.89%). In the CHR-P group, only the ADS SA was associated with positive symptoms (ß = -0.08; 95% CI, -0.13 to -0.02; P = .02 for false discovery rate) and IQ (ß = 0.09; 95% CI, 0.02-0.15; P = .02 for false discovery rate). Conclusions and Relevance: In this case-control study, findings suggest that macroscale neuromorphometric measures may not provide an adequate explanation of psychosis risk.

Examining the association between posttraumatic stress disorder and disruptions in cortical networks identified using data-driven methods.

Posttraumatic stress disorder (PTSD) is associated with lower cortical thickness (CT) in prefrontal, cingulate, and insular cortices in diverse trauma-affected samples. However, some studies have failed to detect differences between PTSD patients and healthy controls or reported that PTSD is associated with greater CT. Using data-driven dimensionality reduction, we sought to conduct a well-powered study to identify vulnerable networks without regard to neuroanatomic boundaries. Moreover, this approach enabled us to avoid the excessive burden of multiple comparison correction that plagues vertex-wise methods. We derived structural covariance networks (SCNs) by applying non-negative matrix factorization (NMF) to CT data from 961 PTSD patients and 1124 trauma-exposed controls without PTSD. We used regression analyses to investigate associations between CT within SCNs and PTSD diagnosis (with and without accounting for the potential confounding effect of trauma type) and symptom severity in the full sample. We performed additional regression analyses in subsets of the data to examine associations between SCNs and comorbid depression, childhood trauma severity, and alcohol abuse. NMF identified 20 unbiased SCNs, which aligned closely with functionally defined brain networks. PTSD diagnosis was most strongly associated with diminished CT in SCNs that encompassed the bilateral superior frontal cortex, motor cortex, insular cortex, orbitofrontal cortex, medial occipital cortex, anterior cingulate cortex, and posterior cingulate cortex. CT in these networks was significantly negatively correlated with PTSD symptom severity. Collectively, these findings suggest that PTSD diagnosis is associated with widespread reductions in CT, particularly within prefrontal regulatory regions and broader emotion and sensory processing cortical regions.

The Influence of Brain MRI Defacing Algorithms on Brain-Age Predictions via 3D Convolutional Neural Networks.

In brain imaging research, it is becoming standard practice to remove the face from the individual's 3D structural MRI scan to ensure data privacy standards are met. Face removal - or 'defacing' - is being advocated for large, multi-site studies where data is transferred across geographically diverse sites. Several methods have been developed to limit the loss of important brain data by accurately and precisely removing non-brain facial tissue. At the same time, deep learning methods such as convolutional neural networks (CNNs) are increasingly being used in medical imaging research for diagnostic classification and prognosis in neurological diseases. These neural networks train predictive models based on patterns in large numbers of images. Because of this, defacing scans could remove informative data. Here, we evaluated 4 popular defacing methods to identify the effects of defacing on 'brain age' prediction - a common benchmarking task of predicting a subject's chronological age from their 3D T1-weighted brain MRI. We compared brain-age calculations using defaced MRIs to those that were directly brain extracted, and those with both brain and face. Significant differences were present when comparing average per-subject error rates between algorithms in both the defaced brain data and the extracted facial tissue. Results also indicated brain age accuracy depends on defacing and the choice of algorithm. In a secondary analysis, we also examined how well comparable CNNs could predict chronological age from the facial region only (the extracted portion of the defaced image), as well as visualize areas of importance in facial tissue for predictive tasks using CNNs. We obtained better performance in age prediction when using the extracted face portion alone than images of the brain, suggesting the need for caution when defacing methods are used in medical image analysis.

Cortical microstructural associations with CSF amyloid and pTau.

Diffusion MRI (dMRI) can be used to probe microstructural properties of brain tissue and holds great promise as a means to non-invasively map Alzheimer's disease (AD) pathology. Few studies have evaluated multi-shell dMRI models such as neurite orientation dispersion and density imaging (NODDI) and mean apparent propagator (MAP)-MRI in cortical gray matter where many of the earliest histopathological changes occur in AD. Here, we investigated the relationship between CSF pTau181 and Aß1-42 burden and regional cortical NODDI and MAP-MRI indices in 46 cognitively unimpaired individuals, 18 with mild cognitive impairment, and two with dementia (mean age: 71.8 ± 6.2 years) from the Alzheimer's Disease Neuroimaging Initiative. We compared findings to more conventional cortical thickness measures. Lower CSF Aß1-42 and higher pTau181 were associated with cortical dMRI measures reflecting less hindered or restricted diffusion and greater diffusivity. Cortical dMRI measures, but not cortical thickness measures, were more widely associated with Aß1-42 than pTau181 and better distinguished Aß+ from Aß- participants than pTau+ from pTau- participants. dMRI associations mediated the relationship between CSF markers and delayed logical memory performance, commonly impaired in early AD. dMRI metrics sensitive to early AD pathogenesis and microstructural damage may be better measures of subtle neurodegeneration in comparison to standard cortical thickness and help to elucidate mechanisms underlying cognitive decline.

A Comprehensive Corpus Callosum Segmentation Tool for Detecting Callosal Abnormalities and Genetic Associations from Multi Contrast MRIs.

Structural alterations of the midsagittal corpus callosum (midCC) have been associated with a wide range of brain disorders. The midCC is visible on most MRI contrasts and in many acquisitions with a limited field-of-view. Here, we present an automated tool for segmenting and assessing the shape of the midCC from T1w, T2w, and FLAIR images. We train a UNet on images from multiple public datasets to obtain midCC segmentations. A quality control algorithm is also built-in, trained on the midCC shape features. We calculate intraclass correlations (ICC) and average Dice scores in a test-retest dataset to assess segmentation reliability. We test our segmentation on poor quality and partial brain scans. We highlight the biological significance of our extracted features using data from over 40,000 individuals from the UK Biobank; we classify clinically defined shape abnormalities and perform genetic analyses.

Evaluating Fiber Orientation Dispersion Measures Computed From Single-Shell Diffusion MRI.

Fiber orientation dispersion is one of the fundamental features that can be estimated from diffusion magnetic resonance imaging (dMRI) of the brain. Several approaches have been proposed to estimate dispersion from single- and multi-shell dMRI acquisitions. Here, we derive solutions to bring these proposed methods to a standard orientation dispersion index (ODI) with the goal of making them comparable across different dMRI acquisitions. To illustrate the utility of the measures in studying brain aging, we further examined the age-dependent trajectory of the different single- and multi-shell ODI estimates in the white matter across the lifespan.Clinical Relevance- This work computes metrics of brain microstructure that can be adapted for large neuroimaging initiatives that aim to study the brain's development and aging, and to identify deviations that may serve as biomarkers of brain disease.

Brain age predicted using graph convolutional neural network explains neurodevelopmental trajectory in preterm neonates.

OBJECTIVES: Dramatic brain morphological changes occur throughout the third trimester of gestation. In this study, we investigated whether the predicted brain age (PBA) derived from graph convolutional network (GCN) that accounts for cortical morphometrics in third trimester is associated with postnatal abnormalities and neurodevelopmental outcome. METHODS: In total, 577 T1 MRI scans of preterm neonates from two different datasets were analyzed; the NEOCIVET pipeline generated cortical surfaces and morphological features, which were then fed to the GCN to predict brain age. The brain age index (BAI; PBA minus chronological age) was used to determine the relationships among preterm birth (i.e., birthweight and birth age), perinatal brain injuries, postnatal events/clinical conditions, BAI at postnatal scan, and neurodevelopmental scores at 30 months. RESULTS: Brain morphology and GCN-based age prediction of preterm neonates without brain lesions (mean absolute error [MAE]: 0.96 weeks) outperformed conventional machine learning methods using no topological information. Structural equation models (SEM) showed that BAI mediated the influence of preterm birth and postnatal clinical factors, but not perinatal brain injuries, on neurodevelopmental outcome at 30 months of age. CONCLUSIONS: Brain morphology may be clinically meaningful in measuring brain age, as it relates to postnatal factors, and predicting neurodevelopmental outcome. CLINICAL RELEVANCE STATEMENT: Understanding the neurodevelopmental trajectory of preterm neonates through the prediction of brain age using a graph convolutional neural network may allow for earlier detection of potential developmental abnormalities and improved interventions, consequently enhancing the prognosis and quality of life in this vulnerable population. KEY POINTS: •Brain age in preterm neonates predicted using a graph convolutional network with brain morphological changes mediates the pre-scan risk factors and post-scan neurodevelopmental outcomes. •Predicted brain age oriented from conventional deep learning approaches, which indicates the neurodevelopmental status in neonates, shows a lack of sensitivity to perinatal risk factors and predicting neurodevelopmental outcomes. •The new brain age index based on brain morphology and graph convolutional network enhances the accuracy and clinical interpretation of predicted brain age for neonates.

Multiomic approach and Mendelian randomization analysis identify causal associations between blood biomarkers and subcortical brain structure volumes.

Alterations in subcortical brain structure volumes have been found to be associated with several neurodegenerative and psychiatric disorders. At the same time, genome-wide association studies (GWAS) have identified numerous common variants associated with brain structure. In this study, we integrate these findings, aiming to identify proteins, metabolites, or microbes that have a putative causal association with subcortical brain structure volumes via a two-sample Mendelian randomization approach. This method uses genetic variants as instrument variables to identify potentially causal associations between an exposure and an outcome. The exposure data that we analyzed comprised genetic associations for 2994 plasma proteins, 237 metabolites, and 103 microbial genera. The outcome data included GWAS data for seven subcortical brain structure volumes including accumbens, amygdala, caudate, hippocampus, pallidum, putamen, and thalamus. Eleven proteins and six metabolites were found to have a significant association with subcortical structure volumes, with nine proteins and five metabolites replicated using independent exposure data. We found causal associations between accumbens volume and plasma protease c1 inhibitor as well as strong association between putamen volume and Agouti signaling protein. Among metabolites, urate had the strongest association with thalamic volume. No significant associations were detected between the microbial genera and subcortical brain structure volumes. We also observed significant enrichment for biological processes such as proteolysis, regulation of the endoplasmic reticulum apoptotic signaling pathway, and negative regulation of DNA binding. Our findings provide insights to the mechanisms through which brain volumes may be affected in the pathogenesis of neurodevelopmental and psychiatric disorders and point to potential treatment targets for disorders that are associated with subcortical brain structure volumes.