NeuroImage

Clinical brain MRIs from pediatric health systems represent a viable resource for modeling early neurodevelopmental trajectories and studying neurodevelopmental risk in real-world populations. However, a limitation to date has been the performance of existing segmentation tools for measuring various brain phenotypes in clinical scans. In particular, many tools underperform in infant scans due to morphological and physical changes such as rapid myelination. Here, we introduce ClinSeg: a robust segmentation approach tailored to early-life clinical MRIs with variable orientation, resolution, and contrast. We leverage existing registration and synthetic data generation tools to construct a training corpus for a 3d U-Net spanning anatomical and contrast diversity, including scans with morphological abnormalities from a pediatric hospital. Validated against manual segmentations, ClinSeg outperforms existing models in infancy while matching them in childhood and adolescence. Finally, ClinSeg enables the construction of reference brain growth trajectories in 11,699 individuals from 0-21 years of age, leading to the detection of more nuanced age-related findings in clinical groups.

Bayesian borrowing of external information can improve trial efficiency, particularly in pediatric and rare disease settings where patient populations are limited, but may introduce bias and inflate the Type~I error rate when the trial differs from external studies. Recent U.S. Food and Drug Administration (FDA) draft Bayesian guidance emphasizes careful evaluation of external information, prior specification, and assessment of operating characteristics. This paper compares three meta-analytic-predictive (MAP)-based methods for Bayesian borrowing: the MAP prior, robust MAP (RMAP) prior, and self-adapting mixture (SAM) prior. An adaptive platform trial design in Rett syndrome is used as a case study. Simulation studies evaluate frequentist operating characteristics under varying prior--data conflict, between-study heterogeneity, treatment effects, and clinically significant differences (CSDs) for the SAM prior. The MAP prior achieved the greatest efficiency when external and current data were compatible but exhibited the largest bias under substantial prior--data conflict. The RMAP priors improved robustness through fixed robust-component weights, whereas the SAM prior adaptively adjusted borrowing and was less sensitive to prior--data conflict while retaining efficiency gains when the data were compatible. Although the CSD influenced the degree of adaptive borrowing, as reflected by effective sample size, it had only a modest impact on frequentist operating characteristics. Sensitivity analyses using a skeptical robust component yielded similar qualitative conclusions, while accentuating the differences between the MAP and RMAP priors. These findings provide guidance for evaluating and selecting MAP-based borrowing strategies before trial implementation, particularly in rare disease settings, consistent with current FDA recommendations.

Rib-cage morphology is a determinant of thoracic biomechanics, ventilation, and injury response, yet statistical shape models (SSMs) of the rib cage have relied on small cohorts (~100s of individuals) imaged by clinical computed tomography, which over-represents injury and disease. We constructed a surface-based SSM of the complete 24-rib cage from 26,275 standardised whole-body magnetic resonance imaging (MRI) scans of adults aged 19-74 years from the population-based German National Cohort (NAKO). Ribs were segmented with a deep-learning pipeline (a rib-extended SPINEPS model), reconstructed as per-rib surface meshes, and brought into dense vertex-wise correspondence by Gaussian-process morphable registration in Scalismo; the aligned ensemble was summarised by generalised Procrustes analysis and principal component analysis (PCA). Fourteen per-rib geometric descriptors provided a quantitative cross-walk between the abstract PCA modes and named shape features, and associations with sex, age, body size and composition (including body-fat percentage), and smoking exposure were estimated by multivariable regression with Benjamini-Hochberg false-discovery-rate control. Shape variation was strongly concentrated: 28 modes captured 95% of the total variance, and the first three alone accounted for 69.4% (PC1, 42.6%; PC2, 16.3%; PC3, 10.5%) and admitted consistent anatomical readings - a sexually dimorphic axis (PC1), a slender-versus-stout body-habitus contrast (PC2), and a free-rib-size axis at ribs 11-12 (PC3). The sexes were nearly fully separated along PC1 (Cohen's d = 2.52). Body mass and body-fat percentage were the dominant modifiable correlates of rib-cage shape, whereas the association with cumulative smoking exposure was comparatively small. The model is released as a population-representative geometric reference for benchmarking and morphing donor-derived finite-element human-body models and for further large-cohort shape analysis.

Automated lesion segmentation is increasingly central to PSMA PET/CT interpretation, supporting staging, treatment planning, and response assessment at a scale that outpaces available nuclear-medicine expertise. However, automated PSMA-PET/CT whole-body lesion segmentation models are trained on images alone, with no knowledge of where in the body prostate metastases actually tend to occur. Radiologists use clinical domain knowledge of metastatic spread, but its absence in machine learning models produces false positives in anatomically implausible locations and missed lesions in high-risk sites such as the liver. In this work, we explore whether population-level spatial knowledge of metastatic spread can be used to augment deep learning segmentation predictions, and how such a prior should be fused with a network's output, without additional training. We build a data-driven metastasis atlas from 375 expert-annotated whole-body PSMA PET/CT scans and investigate its fusion with a trained segmentation network under a Bayesian framework, in which prediction probabilities from an nnU-Net-based lesion segmentation model serve as the likelihood and the data-driven atlas as the prior. Because metastases occupy only a small fraction of whole-body voxels, the atlas's peak probability is too low, and standard power-scaled or naive Bayesian pooling references lack the tools to deal with this shortcoming. This causes these standard fusion strategies to fail and, in the naive Bayesian case, to sharply degrade performance. We instead derive a calibrated, background-referenced log-odds fusion, one of many possible approaches to combine a population atlas with a deep learning model's predictions, distinct from classical multi-atlas label fusion in that it fuses a single population prior with a trained network's softmax rather than combining several registered atlases. Furthermore, this approach is neutral outside atlas support by construction, reduces exactly to the baseline network when unweighted, and requires no retraining. This atlas fusion significantly improved mean Dice over the baseline nnU-Net on a disjoint internal test set ($+0.011$, Holm-adjusted $p=0.021$) and on an independent external cohort ($+0.0129$, Holm-adjusted $p=3.8\times10^{-16}$), with lesion sensitivity improving from 0.849 to 0.861 internally and Dice improving over baseline in every stratified anatomic region, including the rare, high-risk sites motivating this work, while naive Bayesian pooling degrades performance sharply and power-scaled pooling underperforms it throughout. Our findings suggest that population-level spatial priors can meaningfully augment deep learning predictions in whole-body oncologic segmentation, provided the fusion rule is calibrated to where the prior actually carries signal.

Background: Serial glial fibrillary acidic protein (GFAP) trajectories have become an important framework for contextualizing evolving secondary-injury pathophysiology after moderate-to-severe traumatic brain injury (msTBI). However, total GFAP pools release and clearance signals that may be less useful for longitudinal bedside decisions than a proteoform-resolved assay. We compared total GFAP with neoGFAP, defined here as calpain-generated GFAP proteoforms intended to index active astroglial proteolysis during the subacute phase. Methods: We analyzed 651 serial serum samples from 95 msTBI patients from a previously described single-site cohort. Total GFAP and neoGFAP were measured on the same MSD platform from 6 to 240 hours after injury. Early (6 to 72 h) and late (96 to 240 h) windows, data-derived tertiles, and serial trajectory summaries were calculated directly from serial samples. Models were benchmarked against age plus admission post-resuscitation Glasgow Coma Scale (GCS) and the admission IMPACT extended risk score using five-fold stratified cross-validation. Outcomes were unfavorable outcome (GOSE 1 to 4), less-than-good recovery (GOSE 1 to 6), Disability Rating Scale (DRS) [≥]15, mortality, and neuroimaging worsening at 6 months. Results: The cohort contributed 95 serial biomarker profiles, with 90 participants evaluable for 6-month GOSE and 89 for DRS. Unfavorable outcome occurred in 57/90 (63.3%), and less-than-good recovery in 79/90 (87.8%). For unfavorable outcome, IMPACT plus early neoGFAP reached AUROC 0.85 versus 0.84 for IMPACT plus early total GFAP and 0.81 for IMPACT alone. For less-than-good recovery, IMPACT plus late neoGFAP achieved AUROC 0.90 versus 0.84 for late total GFAP and 0.82 for IMPACT alone. Secondary analyses for DRS, mortality, and neuroimaging worsening showed smaller differences. Conclusions: In this retrospective analysis, neoGFAP provided clearer incremental value than total GFAP for recovery-oriented monitoring, especially when late-window reassessment of patients who remained at risk for less-than-good recovery was required. Results support prospective testing of neoGFAP as a pathophysiology-informed adjunct to serial bedside decision making, repeat-assessment thresholds, and recovery stratification.

Dynamic positron emission tomography (PET) with long axial field-of-view (LAFOV) scanners enables multi-organ imaging and kinetic quantification beyond static (late-phase) imaging; however, the long times typically required for dynamic acquisitions remain clinically impractical. This study evaluates a deep learning (DL) framework to enable abbreviated dynamic PET acquisitions, comparing single-time-window (STW, early dynamic data only) and dual-time-window (DTW, early dynamic data plus a late 5-min static frame) protocols with early dynamic scan durations of 5-30 min and dose levels ranging from 360 MBq to 18 MBq. Seventeen 60-min dynamic [18F]FDG datasets were first motion-corrected using a staggered FALCON pipeline and then used to train and test a spatiotemporal DL model for autoregressive frame prediction. Performance was assessed across the full quantitative workflow, from DL-predicted frames and time-activity curves to organ-based kinetic modeling and voxel-wise parametric imaging in multiple tissues and two patient cohorts. DTW protocols consistently outperformed STW, better preserving late-phase kinetics. For a 15-min early dynamic scan, adding a late 5-min scan reduced mean absolute Ki difference from 23% (STW) to 17% (DTW) in the liver and from 26% to 15% in the thalamus. DTW + DL further reduced errors to [≤]10% in the liver, thalamus, and breast lesion, and 16% in muscle. Our recommended protocol, 15-min early dynamic scan plus a 5-min late scan with DL, remained robust to up to a 5-fold dose reduction (~74 MBq). Overall, these findings support DL-enabled abbreviated, low-dose dynamic LAFOV PET as a clinically feasible approach for accurate kinetic quantification

Diffusion models achieve high reconstruction quality in low-dose computed tomography (LDCT), but their iterative sampling trajectories impose substantial computational costs. Unlike unconditional generation, paired LDCT reconstruction starts from an image that already contains the anatomy and spatial structure of the standard-dose CT (SDCT) target; reconstruction primarily requires correcting dose-related noise and artifacts. We therefore introduce Residual Endpoint Flow Matching (REFM), an LDCT reconstruction method that learns to transport an LDCT image directly toward its paired SDCT endpoint rather than defining a noise-to-image trajectory. REFM predicts the residual velocity along linear interpolations between both images and supports single-step and multi-step reconstruction using the same trained network. We evaluate five model capacities using 1 to 50 Euler steps against deterministic U-Net and diffusion-based baselines. Across all REFM variants, one-step inference consistently provides the highest reconstruction quality. On the TCIA validation set, REFM Base achieves 50.98 dB PSNR and 0.9865 SSIM at 94.54 fps, compared with 50.92 dB, 0.9847, and 9.26 fps for DDPM-10. REFM Small retains 50.71 dB while increasing throughput to 198.56 fps. Without fine-tuning, REFM Base also matches the 25-step DDPM baseline on the external Mayo Clinic dataset, although DDPM remains stronger on synthetically degraded CRLM images. Thus, our results show that exploiting paired anatomical correspondence enables diffusion-level LDCT reconstruction with a single step reconstruction.

Abstract: Objective: Post-traumatic epilepsy (PTE) is a common sequela of traumatic brain injury (TBI). Research indicates that individuals with PTE tend to experience greater cognitive difficulties compared to those with TBI alone. However, it is plausible that a distinct cognitive profile exists that distinguishes between TBI cases with and without PTE. We aimed to identify longitudinal changes in cognitive measures among TBI patients to better assess the changes associated with developing PTE. Setting: Outpatient. Participants: Prospective subjects who had suffered TBI within 6 months post-injury (TBI-6M, n=32), retrospective subjects with pre-existing PTE diagnoses (PTE, n=20), and healthy control subjects (HC, n=41). Design: We examined cognitive performance for TBI patients within 6 months post-injury, then again within 12 months (TBI-12M, n=26), and within 18-months (TBI-18M, n=25), and compared this with cognitive performance among HC and PTE. Main Measures: Cognitive tests administered yielded 15 test components for analysis. We utilized linear mixed effects modeling to examine cohort-level differences cognitive function. Results: 11/15 tests showed a significant performance deficit in the PTE subjects compared to HC. TBI-6M was not significantly different from the PTE subjects; with time, 9/15 tests showed some degree of recovery in TBI subjects. Tests for information processing speed/working memory and executive function showed strong recovery (TBI-6M vs. TBI-18M, SDMT written: p<0.0001, SDMT oral and COWAT: p<0.001). Tests for visual attention/working memory also showed a smaller but significant recovery (TBI-18M vs. PTE, p<0.05). By contrast, tests for verbal memory [HVLT-R Delayed Recall] showed chronic impairment in TBI (TBI-18M vs HC, p<0.0001). TBI subjects generally trend towards recovery in cognitive performance post-TBI. Conclusions: Information processing speed/working memory are strong indicators for TBI recovery, while auditory learning/memory shows chronic impairment. The stagnation of recovery in cognitive domains typically characterized by robust recovery may correlate with an elevated risk of developing PTE.

Objective: Published P300-speller fusion schemes fix prior trust regardless of trial reliability; we tested whether a reliability estimate improves on it. Methods: We reanalyzed 3,373 archived P300-speller selections from 47 people with ALS (BigP3BCI). A fair, matched-search-space comparison, tuning both a fixed weight and an adaptive policy out-of-fold, was evaluated across 22 evaluable language-model priors up to 46.7B parameters. Two representative priors, GPT-2 and a classical 5-gram, additionally received detailed naive and mechanistic analyses. Results: No prior's 95% CI favored adaptive fusion under the fair comparison, despite unexploited oracle headroom at every scale. Under GPT-2, the naive comparison was significantly worse for adaptive fusion; both anchors converged to a degenerate or near-degenerate fair-comparison solution. For the representative anchors, three further controllers failed to convert that headroom into benefit; the fixed-fused posterior's output probability outperformed the best controller for flagging errors (2.8- to 3.8-fold enrichment). Conclusion: A tuned fixed weight is a difficult-to-beat default across the tested scale range; reliability estimation gave no deployable adaptive advantage. Significance: Adaptive weighting should be validated against a fairly tuned baseline across model families and scales; in this dataset, the fused output's confidence identified high-risk selections better than the tested purpose-built ranker.

Background: The Centiloid (CL) scale standardizes global amyloid PET quantification and is widely used to define amyloid positivity. As a global summary measure, however, CL may not fully reflect the regional distribution of amyloid deposition, which can carry additional prognostic information about the rate of cognitive decline. Objective: To develop and externally validate a fixed, regional amyloid PET composite score that complements CL for predicting cognitive decline in Alzheimer's disease. Methods: The Regional Amyloid PET Score (RAPS) was derived from 82 FreeSurfer regions using machine learning with bootstrap stability selection to predict the rate of change in CDR-Sum of Boxes (CDR-SB) in 433 amyloid-positive ADNI [18F]florbetapir participants. The fixed nine-region weights were applied without retraining in a cross-tracer ADNI [18F]florbetaben subset (N = 71; largely overlapping the discovery participants) and two external validation cohorts, NACC SCAN (N = 1531; four tracers) and OASIS-3 (N = 428). Results: RAPS comprised nine regions. In ADNI, RAPS correlated more strongly with CDR-SB slope than CL and showed higher discrimination of rapid decliners (AUC 0.813 vs 0.713). Performance was directionally consistent across validation cohorts; in NACC SCAN, RAPS and CL independently predicted clinical progression. Cross-cohort meta-analysis of the three independent cohorts supported incremental discrimination beyond CL (pooled {Delta}AUC +0.066; I2 = 0%). Conclusions: RAPS, a fixed regional amyloid PET-derived score, may complement CL for prognostic stratification in Alzheimer's disease research.

Background Understanding the mechanisms underlying brain aging and age-related pathological changes is essential for advancing brain health research. Our group previously developed a mechanistic mathematical model of healthy brain, Chamberland et al. (2024) that integrates key biological processes involved in normal aging, from which Alzheimer's disease (AD) related changes may emerge naturally. Objectives To characterize and validate this brain model by evaluating its sensitivity, calibrating its parameters, and assessing generalizability in independent populations. Methods The model represents the evolution of key biological processes associated with brain aging, including amyloid beta (A{beta}), tau pathologies, neuroinflammation, and neuronal death. After identifying the 30 most influential parameters, we calibrated the model using cognitively normal (CN) participants from the AD Neuroimaging Initiative (ADNI) database (n = 211) by minimizing a loss function composed of three outcomes (AB) plaques, tau tangles, and neuronal density). The calibrated model was then applied to the UK Biobank cohort (n = 35,899) of normal controls (aged 44-82 years). The effects of sex and APOE were evaluated using stratified simulations. Results Parameter calibration significantly reduced the prediction errors for A{beta} and tau. Neuronal density predictions showed strong agreement in the UK Biobank cohort. The variance decomposition identified APOE status as a major contributor to variability in A{beta}. Conclusion Our validated brain health model links mechanistic pathways with population data and reproduces neuronal density patterns in an independent cohort. These findings support its use as a framework for studying brain aging and investigating how Alzheimer's disease related pathological changes may emerge with aging.

Background: Hippocampal sclerosis (HS) is a common cause of drug-resistant focal epilepsy (DRFE) and amenable to neurosurgical treatment. Detection relies on MRI but can be challenging. 7 Tesla (T) ultra-high field MRI and automated MRI post-processing tools have independently been shown to improve radiological diagnosis of HS. However, combining these approaches remains underexplored. This study evaluated whether AID-HS, a tool for HS detection developed using 3T MRI, generalises to 7T MRI data. Methods: We collated a dataset of paired 3T and 7T T1-weighted MRI from four epilepsy centres, including 23 patients with HS, 39 healthy controls, and 23 individuals with focal cortical dysplasia as disease controls. Histopathology served as the gold standard for defining HS where available (n=7), otherwise radiological findings (n=16). AID-HS was applied to images acquired at both field strengths, and sensitivity and specificity for detection and lateralisation of HS were compared. Additionally, agreement of hippocampal features across 3T and 7T was evaluated. Results: We found no evidence of a difference in performance of AID-HS between 3T and 7T. Sensitivity for detection of unilateral HS was 63% (12/19) at 3T and 68% (13/19) at 7T (McNemar's exact test p=1.0). Specificity in controls was 97% (60/62) at 3T and 100% (62/62) at 7T (p=0.5). Bilateral HS was correctly flagged in 3 of 4 cases using feature-based criteria, with high specificity in controls. Quantitative hippocampal features showed moderate to good agreement across field strengths (ICC 0.70 to 0.98), with small differences observed for volume and thickness estimates. Conclusion: AID-HS provides robust detection and lateralisation of HS across multiple 7T MRI centres, highlighting its potential to enhance lesion detection. Future work is needed to investigate whether models trained on 7T data can leverage the improved image quality for further gains in HS detection performance.

Despite extensive knowledge of the progressive sequence of cognitive and functional deficits in Alzheimer's Disease (AD), the impact of neurodegeneration on the conscious experience of patients remains largely unexplored. Understanding how the content of consciousness, particularly perceptual awareness, changes with the progression of AD is crucial to enable meaningful person-centred care. This is especially important in severe AD when impairments in language and other cognitive domains mean people are unable to report their experiences. We investigated whether electrophysiological (EEG) and fMRI signatures of perceptual awareness described in healthy older people are present in people with mild-moderate and severe AD using two "no-report" paradigms. Firstly, a visual masking paradigm examined visual awareness negativity (VAN) and late positive (LP) electrophysiological responses and activation in visual cortex and fronto-parietal regions that are characteristically associated with conscious perception of faces; and second, a complex audio-visual (movie) task examined activation in fronto-parietal networks previously associated with perceptual awareness. In healthy older controls we found cortical responses characteristic of awareness in both EEG and fMRI modalities, with VAN and LP markers and widespread occipital, fusiform face area and fronto-parietal activation. In people with mild-moderate AD, there were significant reductions in VAN and LP markers and reduced fronto-parietal activation. In participants with severe AD, who were behaviourally minimally responsive, there was only limited evidence of presence of frontoparietal markers of perceptual awareness, however this may reflect attentional and task insensitivity in people with advanced dementia. These results demonstrate that the brain mechanisms associated with perceptual awareness become increasingly impaired with progression of AD. Specifically, involvement of frontoparietal networks is reduced in AD, which may reflect reduced higher-level awareness. This suggests AD should be considered a disorder of consciousness and should motivate further investigation into the dimensions of awareness affected by the disorder with implications for treatment and management of people with dementia.

Background: Storage symptoms after stroke-isolated urgency, urgency with frequency, and isolated frequency are common but traditionally attributed to a single overactive bladder mechanism via suprapontine disinhibition. However, clinical heterogeneity in symptom presentation suggests distinct underlying mechanisms. We aimed to characterize the neural substrates of three storage symptom subtypes after stroke using comprehensive lesion-symptom mapping. Methods: We prospectively evaluated 1,498 consecutive subacute stroke patients admitted for inpatient rehabilitation (1,105 men, 73.8%; median age 61 years). Storage symptoms were classified into three subtypes: isolated urgency (n=109), urgency with frequency (n=32), and isolated frequency (n=19). Multivariable logistic regression models with Bonferroni correction identified independent predictors across demographic, clinical, white matter hyperintensity (WMH), brain atrophy, and lesion location variables. Results: The three subtypes demonstrated largely distinct sets of independent predictors. The left genu of the corpus callosum (aOR=20.06, 95% CI 7.78-51.74, P<0.001) and the inferior frontal gyrus (aOR=3.48, 95% CI 1.81-6.67, P<0.001) were independently associated with isolated urgency and survived Bonferroni correction, together with a right IFG-insula synergistic effect (OR=21.46, 95% CI 10.49-43.88, P<0.001). Urgency with frequency was associated with a broad fronto-cingulate network-the IFG (aOR=11.45, 95% CI 3.10-42.33, P<0.001, surviving Bonferroni correction) and the ACC (aOR=11.53, 95% CI 2.40-55.49, P=0.002) with diffuse right-hemisphere dominance, older age and brain atrophy. Isolated frequency was associated with anterior corona radiata involvement (aOR=5.46, 95% CI 1.92-15.54, P=0.002) and male sex (aOR=10.62, 95% CI 1.36-82.98, P=0.024), though none reached the strict Bonferroni threshold. Conclusions: These findings identify three mechanistically distinct post-stroke storage symptom subtypes with separable neural substrates, lateralization profiles, and clinical determinants. The triple dissociation across subtypes supports a discrete pathway model over the traditional unitary OAB framework, providing a neuroanatomically grounded basis for subtype-stratified treatment Keywords: storage symptoms; subacute stroke; hemispheric lateralization; structural synergy; lesion-syndrome mapping

Ultrathin endoscopy is highly attractive for real-time tissue imaging in narrow and hard-to-reach regions of the body. A single multimode fibre (MMF) is an attractive probe because of its small diameter, flexibility, and diffraction-limited spatial resolution enabled by the large number of transverse modes guided within a single core. Because the distal fibre tip is inaccessible during endoscopy, reflection-mode imaging, in which the same fibre delivers illumination and collects backscattered light, is more practical than transmission-mode imaging. However, image recovery from the resulting speckle pattern is challenging because light undergoes double-pass propagation through the MMF, with mode coupling and dispersion; the backscattered signal is weak, and the camera records intensity only, without phase information. Here, we propose a single-shot reflection-mode MMF imaging framework that combines a reflected real-valued intensity transmission matrix (reflected-RVITM) with an image restoration network. The reflected-RVITM is calibrated using intensity-only measurements, without interferometry or phase retrieval, and provides a physics-guided initial reconstruction from a single backscattered speckle frame. A restoration network then refines this initial reconstruction instead of inverting the raw speckle. Four restoration backbones are evaluated: HPM-Attention-UNet, GAM, MambaIRv2, and CICPNet. On matched datasets, hybrid models outperformed corresponding networks trained to map raw speckle directly to images. For example, HPM-Attention-UNet on MNIST improved mean PCC from 0.572 to 0.944 (+65.1%). Under domain shift, with training only on Fashion-MNIST and tested on unseen CIFAR scenes, hybrid models achieved mean PCC of 0.61-0.65, compared with 0.36-0.50 for direct learning. This framework is further demonstrated using physical objects at the distal fibre tip. These results demonstrate that a reflected-RVITM physics prior combined with a restoration network enables single-shot image recovery after intensity-only calibration, offering a phase-retrieval-free and generalisable route towards minimally invasive reflection-mode MMF endoscopy.

Autism Spectrum Disorder (ASD) is a heterogeneous neurodevelopmental condition defined by differences in social communication and restricted, repetitive behaviours. As diagnostic criteria have broadened, ASD is now recognised across a wider range of individuals, raising key questions about its structure: does ASD have discrete sub-types, or is it better conceptualised as a continuous, possibly multidimensional, condition? We aim to explore whether a multidimensional continuum model more accurately captures the variability within ASD. We analysed a large SPARK phenotypic dataset of medical history and diagnostic surveys (background history, SCQ, RBS-R; n=36,710 individuals). We apply and compare two traditional statistical approaches, Factor Analysis and Gaussian Mixture Models, with a modern machine learning technique, the Variational Autoencoder (VAE). VAEs reconstructed unseen test data with ~4-fold better accuracy than Factor Analysis, and ~8-fold better accuracy than Gaussian Mixture Models. We identified four stable latent factors across 100 independently trained VAEs. These four dimensions provide an individual behavioural profile that can be visualized using radar-plots, offering a compact way to compare profiles at the person level. Through further analysis, we found evidence for 3 overlapping clusters or subtypes of ASD identified within the 4D latent space. This work aims to inform new ways of modelling ASD using a VAE that will be able to discern between a continuum or a clustered output and that go beyond binary diagnosis, instead reflecting the complex range of trait profiles, with implications for personalised diagnosis and intervention.

Stroke remains a leading cause of mortality and long-term disability worldwide, yet rapid diagnosis is often limited by the shortage of trained radiologists, particularly in resource-constrained settings. Automated analysis of CT imaging offers a potential solution, but existing methods often struggle to achieve clinically generalisable performance while jointly addressing multiple diagnostic tasks. Here we present the Intelligent Integrated Stroke Diagnosis System IISDS, an end-to-end deep learning framework built upon StrokeGNN, a graph-based architecture that integrates 3D contextual feature extraction with U-Net-based 2D lesion segmentation to enable comprehensive stroke analysis from non-contrast CT scans. IISDS performs stroke subtype classification, lesion segmentation and lesion volume estimation within a unified pipeline. To develop and validate the system, we collected and curated BGD-ISD through a collaboration between AI researchers, neurologists, radiologists and clinicians, resulting in a large multi-centre dataset comprising 1,507 CT scans from 597 stroke cases acquired across six hospitals and medical centres in Bangladesh. Across BGD-ISD and multiple publicly available datasets, IISDS achieves state-of-the-art performance on all tasks, improving segmentation accuracy by [&ge;]0.011 Dice score, reducing lesion volume estimation error by [&ge;]0.3 average symmetric surface distance (ASSD), and increasing classification performance by [&ge;]0.018 area under the receiver operating characteristic curve (AUC) compared with existing approaches. These results demonstrate the potential of graph-based deep learning to enable clinically generalisable, automated and scalable stroke diagnosis from CT imaging, supporting rapid clinical decision-making, particularly in healthcare environments with limited access to expert radiological interpretation.

Major depressive disorder (MDD) has been associated with accelerated structural brain aging, yet whether this reflects a pre-existing neurobiological vulnerability, a dynamic acute state effect, or an accumulating biological residual remains unresolved. Across two longitudinal cohorts (N=3220), including a unique sample of 78 initially healthy individuals who transitioned into their first depressive episode during the study course, we systematically tested all three hypotheses. Patients with diagnosed MDD showed elevated MRI-derived brain age relative to healthy controls (1.4 and 2.5 years across cohorts). For the vulnerability hypothesis, individuals scanned prior to their first episode showed no baseline elevation, despite already demonstrating subclinical elevations in self-reported symptom severity, indicating that advanced brain age does not precede illness onset. For the state hypothesis, we found no acceleration of brain aging following the first depressive episode, and longitudinal brain age trajectories were independent of acute clinical symptom severity. Finally, neither episode duration nor recurrence scaled with brain age. Accelerated brain aging in depression is therefore neither an antecedent vulnerability nor an acute state marker of the first episode, but rather a stable biological feature of a long term illness course.

Purpose: The primary objective of the Stanford Knee Osteoarthritis PET/MRI Evaluation (SKOPE) study is to develop and evaluate a multimodal, dynamic [18F]NaF PET-MRI framework for characterizing whole-joint physiology and its relationship to osteoarthritis (OA) risk, pain, and disease progression. Specifically, we aim to integrate dynamic PET with quantitative and anatomical MRI, to characterize structural, compositional, and metabolic features across the knee and surrounding musculoskeletal system, evaluate acute tissue responses to exercise, and identify imaging biomarkers associated with OA risk, pain, and disease progression. Methods: The SKOPE study includes multimodal PET-MRI of the knee and surrounding musculoskeletal tissues, with imaging of the knee, tibia, ankle, thigh, hip, pelvis, and lumbosacral spine. Dynamic [18F]NaF PET is combined with conventional anatomical MRI and quantitative MRI techniques, including quantitative double-echo steady-state (qDESS) T2 mapping of cartilage, Dixon fat-fraction imaging, ultrashort echo time (UTE) T2* mapping of short-T2 tissues, UTE imaging of tibial bone, and zero echo time (ZTE) imaging for bone morphology and pseudo-CT generation. Additional MRI sequences characterize muscle composition, bone and joint anatomy, intervertebral discs, and regional vascular anatomy. Selected scans are acquired before and after a standardized exercise protocol to assess the acute physiological response of the joint. Automated segmentation is used to generate subject-specific masks of muscles, bones, vertebrae, and intervertebral discs. A subset of the MRI protocol is repeated at 1- and 2-year follow-up to assess longitudinal changes. Expected Impact: By combining dynamic bone metabolic imaging with quantitative measures of cartilage, menisci, muscle, bone, fat, vascular structures, and the spine and hip, the SKOPE protocol provides a whole-joint and multijoint framework for studying the structural, metabolic, and physiological processes associated with OA and pain. Exercise and longitudinal imaging further enable assessment of acute tissue responses and changes over time, supporting the development of quantitative imaging biomarkers for OA risk, pain, and disease progression.

Background: Higher Body Mass Index (BMI) is an established risk factor of sleep disturbance. It is not known if the effect is homogeneous across the lifecourse or if there is a particular time point in life that might be best to target. Methods: Two-sample Mendelian randomization (MR) was used to investigated the effect of childhood adiposity (adjusting on adulthood adiposity and obstructive sleep apnea (OSA)) on insomnia, morning chronotype, sleep duration, daytime sleepiness and daytime napping. Similarly, total, and direct effect of adulthood adiposity on these outcomes was explored. We used summary statistics from a genome-wide association study (GWAS) of UK Biobank for childhood and adulthood adiposity (n=453,169) and large-scale consortia of OSA (Million Veteran Program) (n=410,268), insomnia, and chronotype (23andMe) (n=1,978,022 and n=248,1000, respectively). Results: Two-sample univariable MR analysis provided no evidence of an effect of genetically predicted childhood adiposity on later life insomnia (Odds ratio (OR)= 0.94, 95% Confidence interval (CI)= 0.87, 1.03). Whereas, multivariable MR (adjusted for adulthood adiposity) analysis provide strong evidence of direct protective effect of genetically predicted childhood adiposity on later life insomnia (OR= 0.70, CI= 0.64, 0.77). Further, both in univariable and multivariable MR, a strong positive effect of increased childhood body size on morning chronotype was observed (OR= 1.16, CI= 1.01, 1.33 and OR= 1.36, CI= 1.15, 1.62, respectively) after accounting for adulthood body size. In both analysis the estimate did not change considerably after aditionally adjusting for OSA. However, childhood and adulthood adiposity found to be associated with OSA and OSA with insomnia. In both univariable and multivariable analysis, increased body size in adulthood increased the risk of having insomnia and a morning chronotype. Conclusions: The findings suggest that higher body size in childhood is not a risk factor for later life insomnia, whereas higher body size in adulthood was. Further, if healthy body size is maintained in adulthood, high childhood adiposity may decrease the risk of insomnia and increase the risk of being a morning person in later life. Keywords: childhood, adulthood, obesity, insomnia, morning chronotype, medelian randomization