NeuroImage

Presurgical mapping of key white matter (WM) fiber tracts is crucial for intracerebral hemorrhage (ICH) surgery, but it currently relies on tractography from diffusion MRI (dMRI), which has limited applicability in urgent or resource-constrained settings due to long scan times and limited MRI availability. To bridge this gap, we developed a deep learning approach designed to reconstruct critical fiber tracts directly from routinely acquired CT scans, focusing on the corticospinal tract (CST) due to its high clinical relevance. By training a novel network on a curated dataset of 150 paired CT-dMRI scans (101 ICH patients and 49 healthy subjects), we enabled the direct mapping of the CST from CT images alone, bypassing the traditional requirement for dMRI. Our results demonstrate that the CT-derived tracts achieve high anatomical plausibility, with neurosurgeon expert assessments yielding a Likert score of 3.64. Furthermore, the clinical relevance of these reconstructions was validated by a significant correlation between CT-derived tract integrity and patient motor scores (r = 0.726, p = 3.731x 10-7). These findings suggest that in complex clinical scenarios, particularly where dMRI signal quality is compromised by lesion-induced distortion, this CT-based mapping may serve as a useful anatomical reference. Overall, by enabling CT-based WM mapping, our method potentially offers a practical route to gain the key advantages of tractography in resource-limited or time-critical settings

Diffusion MRI tractometry characterizes white matter microstructure along fiber bundles, but standard along-tract profiling collapses measurements across the bundle cross-section, obscuring radial heterogeneity and producing spatially inconsistent units of inference. We present SPECTRA (Spatial Inference for Tractometry), a framework designed to address these limitations through a unified design of parameterization and statistical inference. First, we propose a 2D bundle parameterization that extends along-tract profiling to include a radial dimension defined on the atlas bundle. Second, we develop a two-stage hierarchical false discovery rate (hFDR) procedure for multi-bundle inference, which aggregates evidence at a coarser spatial scale before proceeding to finer-grained inference, with spatial scales derived from a Matern kernel. Across extensive simulation conditions, we found that hFDR improves statistical power and reduces the sample size required to detect effects compared to global FDR correction, while maintaining appropriate error control. We further characterized how sensitivity-specificity tradeoffs depend on sample size, the magnitude, spatial extent, and configurations of effects, thereby providing practical guidance for tractometry study design. In an empirical analysis of mild cognitive impairment and dementia in more than 4,000 subjects across 63 bundles, SPECTRA revealed spatially localized patterns that were absent in 1D profiles. Together, these results demonstrate that spatially resolved parameterization and adaptive error control jointly enable precise mapping of white matter microstructure in large-scale tractometry studies. SPECTRA is openly available as a Python package.

A practical barrier in multimodal neuroimaging is that structural MRI, fMRI, and EEG are often analyzed in modality-specific spaces or reduced to atlas- and sensor-based summaries, limiting the construction of common, interpretable subject-level brain signatures. We evaluate cortical Laplace-Beltrami eigenmode coordinates as a shared geometry-aligned language for structural MRI (sMRI), resting-state fMRI (rs-fMRI), and EEG. In this framework, sMRI morphometric fields are represented by cortical eigenmode coefficients, rs-fMRI by covariance among eigenmode time-series coefficients, and EEG by mode-frequency-condition summaries. Using the Max Planck Institute Leipzig Mind-Brain-Body dataset (MPI-LEMON), we compared unimodal eigenmode-coordinate summaries, multimodal cortical eigenmode-coordinate PCA, conventional atlas/sensor-based PCA and ridge representations, and geometric eigenmode multiview factorization (GEMF). GEMF is a structured decomposition that preserves the modality-native organization of the data objects while separating shared from modality-specific variation. Eigenmode-coordinate representations yielded compact subject-level signatures with strong external validity for chronological age and a secondary cognitive outcome. Multimodal eigenmode-coordinate PCA was among the strongest-performing approaches, reached high age-prediction performance at moderate dimension, and consistently outperformed conventional low-dimensional PCA. GEMF selected an even lower-dimensional shared representation and remained competitive with the benefit of providing interpretable shared and modality-specific factors. These findings support cortical eigenmode coordinates as a practical foundation for compact, interpretable, and multimodally aligned subject-level brain signatures.

Inter-individual variability poses a significant challenge in decoding brain activity across subjects. While standard anatomical registration procedures reduce morphological differences, they fail to capture functional variability between subjects. Functional alignment methods address this issue by establishing voxel-to-voxel correspondences between pairs of individuals, thereby constructing a shared functional space. This shared space can be extended at the group level by generating a functional template. However, despite the availability of toolboxes, functional templates remain underused in fMRI analysis. Adopting this approach is currently difficult due to the diversity of existing methods and the lack of clear guidelines. Comprehensive evaluations of functional templates are limited to movie-watching paradigms. Here, we extensively compare functional alignment methods (Optimal Transport, Procrustes, Ridge regression, and Shared Response Model) and template construction strategies (in-sample, out-of-sample, pairwise) within the more general framework of task decoding. In this framework, decoding accuracy measures how well individual activation patterns align. Across multiple tasks and datasets, we demonstrate that population templates built using Optimal Transport (a) yield the highest decoding accuracy, (b) are not significantly biased by individual subjects, which facilitates generalization to new subjects, and (c) preserve the cortical signal topography.

Diffusion-weighted magnetic resonance imaging (dMRI)-based tractometry enables the quantification of white matter tissue properties in living humans while preserving anatomical specificity. Although tractometry is highly reproducible when the same scanner and acquisition protocol are used, its generalizability across scanners and protocols remains unclear. To address this gap, we performed a traveling-head experiment involving five subjects to evaluate tractometry across progressively different acquisition conditions, including multiple scanners, different scanner models, and two distinct protocols. Tractometry was performed for 20 major white matter tracts using diffusion tensor imaging metrics, neurite orientation dispersion and density imaging (NODDI) metrics, and a semi-quantitative ratio metric (T1w/b0). Generalizability across dataset pairs was quantified using the intraclass correlation coefficient (ICC). Tractometry showed consistently high ICCs when the scanner and protocol were identical; however, ICCs declined as differences in scanner model and acquisition protocol increased. Fractional anisotropy and orientation dispersion index retained relatively high ICCs across these comparisons, whereas other metrics showed marked declines when scanners or protocols differed. ComBat harmonization partially mitigated these declines, but ICCs did not reach the levels observed for datasets acquired using identical scanners and protocols. Finally, the minimum detectable change (MDC) for tractometry in datasets pooled across scanners and protocols varied by tract; for example, the optic radiation showed a lower MDC than the cingulum hippocampus. These findings highlight both the strengths and limitations of tractometry in multisite studies and highlight the importance of quantifying scanner- and protocol-dependent effects for specific metrics and tracts when interpreting measurements from heterogeneous datasets.

There is growing interest in identifying brain function underlying adolescent cognition, personality, and psychopathology. One promising approach is Precision Functional Mapping (PFM) of MRI functional connectivity, a data-intensive method for characterizing individualized brain networks. Foundational studies suggest that PFM can detect stable, task-responsive, and clinically relevant networks. Studies demonstrate that both functional connectivity reliability and network stability improve with increasing data quantity, although benchmark estimates vary across populations, preprocessing pipelines, and MRI acquisition approaches. Accordingly, it is important to understand how PFM performs in adolescent populations and with multi-echo fMRI acquisition. In a case study of eight youth (ages 10-17), we applied PFM to 80-minutes of combined resting-state and task-based fMRI. The resulting networks were highly modular, consistent with adult templates, and without evidence of structural registration artifacts. Functional connectivity reliability compared favorably to prior single-echo studies, with multivariate similarity and ICC estimates showing early stabilization around 10-15 minutes despite continued improvement with additional data. Trait-like stability increased gradually with acquisition time and a Bayesian algorithm (MS-HBM) demonstrated higher stability than Infomap. Across algorithms, stability was greatest in sensory networks (somatomotor, auditory, visual). Furthermore, when evaluating task-based responses to threat and attention paradigms, only the auditory network consistently benefited from individualized mapping over group template networks. These findings suggest that, with constrained scanning time, PFM is especially effective for characterizing sensory and perceptual networks in adolescents. Bridging the methodological divide between deeply sampled individual cases and large-scale developmental studies will require further innovation and validation.

Machine-learning models are increasingly used to predict cognitive and clinical outcomes from neuroimaging data, yet challenges in fairness and generalizability remain. Large-scale datasets are often racially and ethnically imbalanced, leading to systematic performance disparities, with models typically achieving higher accuracy for majority populations represented in the training data. In this study, we evaluated whether supervised domain adaptation methods--including balanced weighting, two-stage TrAdaBoost, feature augmentation with SrcOnly prediction, and linear interpolation--can mitigate these biases. Using the ABCD dataset, we assessed whether models trained on 80 MRI measures from White American participants could generalize more effectively to African American participants. All domain adaptation methods reduced prediction error for African American participants, particularly for MRI modalities with large baseline disparities (e.g., structural MRI), while offering limited improvements where initial gaps were smaller (e.g., functional connectivity). Among the approaches, balanced weighting performed best and remained stable and beneficial even when only 10 African American participants were used to adapt the original model trained exclusively on White American participants. These findings suggest that simple, low-cost strategies can effectively reduce cross-ethnic performance gaps and improve equity in predictive neuroimaging, offering a practical path forward for future neuroimaging predictive biomarkers. Significant StatementLarge-scale neuroimaging datasets increasingly enable machine-learning models to predict cognitive and clinical outcomes; however, these datasets are often ethnically/racially imbalanced. As a result, predictive models tend to generalize poorly to underrepresented populations. We demonstrate that, across 80 MRI phenotypes, a class of machine-learning approaches collectively known as supervised domain adaptation can substantially reduce cross-ethnicity disparities in neuroimaging-based cognitive prediction, even when only limited data from underrepresented groups are available. Among the methods evaluated, balanced weighting achieved the best performance while maintaining low computational cost. Together, these findings provide a practical and scalable framework for improving fairness and generalizability in neuroimaging-based machine learning under realistic conditions of ethnic/racial imbalance.

Most of what we currently know about brain oscillations is derived from Fourier-based spectral power methods. While widely adopted, these procedures introduce methodological limitations and inherently overlook informative neurophysiological features that often remain unreported. In this study, we characterized resting-state oscillatory activity from human magnetoencephalography (MEG) recordings (N = 128) by integrating two complementary approaches. First, oscillatory episodes were detected at the source level with sBOSC. Subsequently, the ByCycle algorithm was applied to these episodes to extract individual cycle features. Results revealed that the brain engages in oscillatory activity for only [~]25% of the recording time, with occipital and parietal regions accounting for the highest temporal prevalence across canonical frequency bands. Furthermore, oscillatory episodes lasted an average of 4.6 cycles, reinforcing the view of neural oscillations as transient bursts. Region-specific duration and power measures revealed distinct anatomical organizations offering complementary physiological information. Finally, by extracting the instantaneous amplitude, period, and waveform asymmetry of individual cycles, we successfully dissociated sinusoidal occipital alpha waves from the asymmetric sensorimotor mu rhythm. By moving beyond traditional power-centric analyses, this approach provides a comprehensive characterization of spontaneous oscillatory activity, thereby offering new insights into the spatial, temporal, and spectral structure of human brain oscillations.

Focusing on an organisms task at hand is instrumental for intelligent and goal-driven behavior. However, humans and other animals often fail to pay sustained attention across long time intervals. Failing to stay on-task may cause one to miss crucial task-relevant signals, leading to impaired performance, which can have serious consequences. Therefore, it is important to understand the neural basis of attentional lapses. One promising neural marker of attentional lapses is the frontal P3 (fP3) EEG component, which has been suggested to reflect the susceptibility to incoming sensory input. Following this, we hypothesized that the fP3 1) predicts imminent lapses of attention, and 2) that it should predict upcoming lapses of attention across modalities. In two experiments, we found that the fP3 reliably tracked lapses of attention of sustained attention already seconds preceding the crucial visual signal. We further extended this to the auditory domain: Already 1.5s ahead of the incoming auditory target, the fP3 revealed whether that target was detected or not. Detailed topographic analyses did, however, reveal a slight dissociation between modalities in underlying intracranial source configurations. In sum, this work revealed a supramodal neural signature of susceptibility, which tracks lapses of sustained attention seconds ahead of the critical incoming sensory input.

BackgroundNeuronavigation based on the standard MNI template (MNI-protocol) offers a cost-effective alternative to the gold-standard individualized T1-weighted MRI approach (T1-protocol). However, it remains unclear whether the reduced anatomical precision of the MNI-protocol compromises its functional efficacy, creating a critical need to verify protocol interchangeability. ObjectiveWe aimed to determine whether the MNI-protocol yields targeting efficacy comparable to the T1-protocol by specifically testing their functional and biophysical equivalence. MethodsWe employed a novel tri-level within-subject framework. The behavioral level assessed functional efficacy via the size congruity effect (SCE) during TMS to the right intraparietal sulcus (IPS). Anatomical accuracy (coil-to-cortex distances) and electromagnetic efficacy (E-field simulations) were evaluated across three distinct regions (right IPS, left dorsolateral prefrontal cortex, and left primary motor cortex) to assess regional generalizability. ResultsThe MNI-protocol demonstrated functional similarity to the T1-protocol, yielding behavioral outcomes that were statistically indistinguishable. This functional equivalence was corroborated by electromagnetic analyses, which revealed nearly identical induced E-field magnitudes and spatial distributions across all three target regions. Although the T1-protocol achieved significantly shorter coil-to-cortex distances, this anatomical advantage did not confer any measurable functional benefit. ConclusionThe MNI-protocol produced behavioral and electromagnetic outcomes equivalent to the T1-protocol. These findings validate the MNI-protocol as a scientifically sound and scalable alternative to individualized MRI-guided targeting, supporting its broader application in diverse research and clinical settings. HighlightsO_LIFunctional equivalence of MNI-vs. T1-guided TMS was systematically tested. C_LIO_LIA novel tri-level framework compared behavioral, anatomical, and E-field metrics. C_LIO_LIMNI- and T1-guided targeting yielded comparable behavioral and E-field outcomes. C_LIO_LIAnatomical proximity does not ensure better behavior or stronger E-field strength. C_LIO_LIMNI-guided targeting offers a robust, practical alternative to individual MRI. C_LI

The thalamus is essential for learning, dynamically engaging with other subcortical and cerebral cortex regions throughout the learning process. Here, the thalamus serves as a critical connector hub and synchroniser within the thalamocortical system of the brain. However, whilst higher order thalamic nuclei are known to be particularly important for this process, the exact contributions of individual higher order and first order thalamic nuclei, alongside their individual involvement with cortical networks and subcortical regions, remains unexplored within the initial phase of learning. In light of this, we analysed fMRI data obtained within a paradigm which is designed to examine initial learning processes within feedback-driven stimulus-response learning, in order to explore thalamic contributions. We investigated dynamic learning-related functional connectivity alterations between various thalamic nuclei with other subcortical regions and cortical networks. Our results show that the initial phase of learning was associated with: (1) decreasing functional connectivity between thalamic nuclei and frontoparietal and cingulo-opercular networks, (2) increasing functional connectivity between thalamic nuclei with default mode and salience networks, (3) decreasing functional connectivity between thalamic nuclei and the putamen, and (4) decreasing functional connectivity amongst higher order thalamic nuclei. Furthermore (5) these dynamic alterations were associated primarily by mediodorsal thalamus. Altogether, these results indicate that higher order thalamic nuclei play a crucial role within initial learning and in the generation of novel goal-directed behaviour. This was demonstrated through enhanced functional connectivity with selected cortical networks which drive goal-directed behaviour, alongside decreased functional connectivity with striatal regions which drive motor selectivity.

Clinical adoption of new biomedical techniques depends on establishing reference values against which individual patients can be compared. In resting-state functional MRI (rsfMRI), most biomarker research has relied on the case-control paradigm, whose underlying assumptions are often invalid as diseases are frequently heterogeneous, limiting biomarker generalizability. Normative modeling offers a complementary alternative by characterizing individual deviations against a reference population. However, in rsfMRI, normative modeling has been applied almost exclusively to functional connectivity, with limited attention to age trajectories and sex effects. We address these gaps by developing a spatial normative model of four rsfMRI metrics that capture complementary features of the blood-oxygen-level-dependent (BOLD) signal across age and sex. Five publicly available datasets were aggregated to form a sample of 1,978 participants aged 10-30 years. Four metrics were computed for each of 110 grey matter regions: amplitude of low-frequency fluctuations (ALFF), fractional amplitude of low-frequency fluctuations (fALFF), regional homogeneity (ReHo), and Hurst exponent. A machine-learning model based on hierarchical Bayesian regression with a non-Gaussian likelihood was fitted per metric, modeling non-linear age effects, sex, and multi-site acquisition. Models were well calibrated across all four metrics, with fALFF showing the strongest predictive performance and Hurst exponent the weakest. Normative trajectories varied across brain regions for each metric, but on average, the median of each distribution remained bounded across regions, while the spread was more regionally variable. All four metrics showed predominantly negative slopes with age, indicating a decrease in each metric over the age window. This work provides a normative reference across four rsfMRI metrics that capture distinct features of the BOLD signal, complementing the case-control paradigm and supporting individual-level inference.

We propose a rigorous and reproducible methodology for analyzing functional MRI data, aimed at: (1) demonstrate their efficiency in classifying task-induced brain states with a limited amount of data, (2) present a methodology to identify brain regions critical for classification and reveal their uniqueness across different states, and (3) show, using strong mathematical methods, that the discriminative power of these regions depends not only on their spatial localization but also on their coordinated temporal activity. Through correlation and temporal structure analyses, we demonstrated that top-ranked regions exhibit stronger, more structured, and richer dependencies than low-ranked regions, underscoring the critical role of temporal dynamics in shaping distinct cognitive brain states. Our work addresses the need for a transparent, accessible, and interpretable framework for studying cognitive processes through neuroimaging data. We analyzed fMRI data from 587 healthy participants from the Human Connectome Project across seven cognitive tasks. Finally, we perform a detailed analysis of the identified brain regions to support further neuroscientific interpretation and discussion. Key PointsO_LIClassical machine learning methods effectively classify task-induced brain states from fMRI data with high accuracy (up to 99% for some tasks), demonstrating that simple, interpretable algorithms can successfully decode complex neuroimaging data without requiring advanced deep learning approaches. C_LIO_LIHigh-accuracy brain states require relatively few significant regions suggesting focal neural signatures, while lower-accuracy states involve more distributed activations across multiple brain areas, revealing different levels of neural organization complexity underlying various cognitive processes. C_LIO_LIThe identified brain regions align with established neuroscientific knowledge, with motor tasks activating contralateral sensorimotor areas, language processing engaging left-hemisphere networks, and social cognition recruiting visual motion processing regions, validating the neurobiological relevance of our machine learning approach. C_LIO_LIRigorous mathematical analyses of temporal dynamics demonstrated that the discriminative power of significant brain regions depends not only on spatial localization but also on their coordinated temporal activity. Correlation, temporal structure analyses consistently showed that top-ranked regions exhibit stronger, more structured, and richer dependencies than low-ranked regions, underscoring the critical role of temporal dynamics in shaping distinct cognitive brain states. C_LI

The relationship between cortical morphology and intelligence during adolescence has been widely studied, with existing literature reporting varying degrees of association across different modeling approaches. This study provides a comprehensive comparison of model performance in investigating the association between crystallized intelligence and cortical surface area using data from 11,351 subjects in the Adolescent Brain Cognitive Development (ABCD) study. We evaluate ten widely used models ranging from linear regression to graph convolutional networks across three covariate adjustment formulations: full (no adjustment), partial (age and sex adjusted), and total surface area (TSA) partial (age, sex, and TSA adjusted). Using bootstrap resampling with 50 iterations, we estimate the fraction of variance explained (FVE) for each model. Our results suggest that more complex models do not lead to higher FVE, with LASSO having the highest FVE of 15.9% (full formulation), Ridge at 10.5% (partial formulation), and Principal Component Regression (PCR) with 102 PCs at 2.5% (TSA partial formulation). Our results also reveal that the relationship between cortical surface area and crystallized intelligence is predominantly driven by global factors age, sex, and TSA, rather than by localized cortical surface area.

In neural evoked responses, small variations in the timing or duration of responses can be observed when the same functional response is recorded in different trials, different experimental conditions or by different sensors. These variations limit the ability of correlation-based methods to detect similarities between signals. Mutual information (MI) provides an alternative similarity measure, capable of capturing both linear and non-linear dependencies, yet its practical use is hindered by lack of consensus on estimators for continuous data and the limited understanding of the behavior of the estimators on realistic signals. In this work, we investigate how to estimate the similarity of neural evoked responses by systematically comparing sample Pearson correlation with three of the most common MI estimators. We describe their behavior using both simulated signals and real magnetoencephalographic data. In the simulations, the estimators are tested against a set of transformations that depict realistic changes in neural evoked responses. Subsequently, we propose guidelines for defining adaptive lower bounds on the similarity estimates and analyzing the similarity rankings induced by the different estimators. Our findings reveal trade-offs between measures sensitivity and different signal properties. We confirm that Pearson correlation is reliable in describing linear relationships for low-noise signals, and we identify parameter settings that stabilize MI estimators, enabling them to capture complex signal dependencies. Together, these results introduce practical parameter choices and thresholding strategies for mutual information and provide guidance for selecting and interpreting similarity measures in the analysis of neural evoked time series.

Big neuroimaging data enable researchers to study subtle structural and functional brain changes and relationships between brain characteristics and genetics, lifestyle, and disease factors. However, substantial effort is needed to minimise technical, non-biological differences between data batches to avoid incorrect inferences. In this study, we address a previously identified bias in UK Biobank FreeSurfer IDPs derived from only the T1 image compared to those using both T1 and T2-FLAIR by treating the bias as a batch effect and using harmonisation approaches. We investigate and characterise this bias through direct within-participant comparison at the image and IDP level, comparing the results with those seen in the wider UKB sample. We then assess different methods of addressing the effect of missing T2-FLAIR, starting from simple linear regression before moving to ComBat, a widely used harmonisation method, testing different approaches for applying ComBat and showing its similarity to simple linear regression. Finally, we examine how ComBat estimates vary with batch and sample size. Our results show clear benefits in using both T1 and T2-FLAIR data in FreeSurfer, as opposed to just the T1, which is more common, with the pial surface fitting being less likely to fail and showing greater biologically plausible inter-subject variability. This is particularly important for cortical thickness IDPs, where T2-FLAIR omission leads to reduced true variability and systematic underestimation, as shown through within-participant repeat testing. We demonstrate that ComBat can address this bias, with its standard use (i.e., applied separately on different IDP categories) showing the best improvement in cortical thickness measures where the bias is strongest, and we find that it is important not to pool ComBat priors across different classes of IDPs. Our proposed version of ComBat with a reference batch (i.e., estimating mean and variance only from data with T2-FLAIR available) performed best in recovering both mean and variance differences between batches across different IDP classes and offers a promising approach for cases where a reference batch is clearly identifiable. While ComBat reliably corrects mean (additive) batch effects with relatively small sample sizes ({approx}30 subjects per batch), we show that its variance (multiplicative) correction is substantially less stable, requiring much larger sample sizes and becoming unreliable when batches are small or imbalanced, or when there is a large variance difference between them.

Generalizable neuroimaging biomarkers that detect cerebral cortical changes after traumatic brain injury (TBI) and predict patient outcomes are needed to improve care and to develop targeted therapies. We used morphometric inverse divergence (MIND) analysis of structural MRI to investigate cortical gray matter morphological networks cross-sectionally and longitudinally after TBI and correlate these with symptoms, disability and cognition six months after injury. Our findings support the Triple Network Model from functional MRI of post-traumatic alterations in the relationship between task-positive, default mode and salience networks. However, the strongest associations between early cortical similarity metrics and long-term patient outcomes involved the dorsal attention network and the limbic network as well as similarity metrics across Mesulam's hierarchy of laminar differentiation. Since MIND mapping of cortical gray matter networks only requires data that is a routine part of standard clinical MRI protocols and does not need image harmonization across different scanners, this work reports a promising new tool that is immediately available for advancing research and clinical care in TBI.

Conventional connectome edge weights, such as number of streamlines (NOS) or diffusion tensor imaging (DTI) metrics, lack specificity to microstructural details which may hold relevance for macroscale brain organisation. Since biophysical diffusion modelling offers greater specificity to microstructure, we investigated whether parameters from the Standard Model of diffusion in white matter provide informative alternatives for connectome weights - namely the intra-axonal signal fraction (f) and perpendicular extra-axonal diffusivity [Formula], as proxies of axonal density and myelination, respectively. Using diffusion MRI data from healthy adults, we constructed structural networks at four parcellation scales, weighted by f, [Formula], NOS, fractional anisotropy (FA) and radial diffusivity (RD). While all weights reproduced expected small-world properties, only [Formula] and normalised NOS captured non-random properties of local organisation across all spatial scales. We then correlated each weighted connectome with resting-state fMRI functional connectivity and intracranial measurements of conduction velocity. At the whole-brain level, although NOS gave strongest coupling with fMRI functional connectivity, only [Formula] exhibited significant structure-function coupling across all spatial scales and modalities. At the regional level, [Formula] and RD gave highest consistency in structure-function coupling across spatial scales. Thus, connectome weights derived from [Formula] capture meaningful aspects of brain network organisation with functional relevance.

Studies in rodents and humans using invasive electrophysiology have established that neural replay is a ubiquitous phenomenon in the brain that is associated with a wide range of cognitive functions, including memory, planning and decision making. Yet, invasively recording in humans remains difficult, and hence knowledge about replay in humans remains scarce. Hence, to comprehensively understand replay in humans, we need reliable approaches that can detect it non-invasively. Several main non-invasive approaches have been proposed, but we lack a full comparative validation against known ground truth signals. In this study, we present FASTIMAGES, a benchmark dataset from seventy participants with parallel fMRI (n = 40, previously published) and MEG (n=30) recordings containing known neural sequences evoked by fast visual stimulation as well as functional localizer trials. The neural sequences were elicited by five different visual stimuli shown in sequences at speeds of 132, 164, 228 and 612 milliseconds onset-to-onset intervals. Using this dataset, we investigate two existing statistical methods for sequence detection, namely Temporally Delayed Linear Modelling (TDLM, developed for MEG by Liu et al., 2021) and Slope Order Dynamic Analysis (SODA, developed for fMRI by Wittkuhn & Schuck, 2021). We examine the underlying assumptions of each method, analyse their resulting strengths and weaknesses in application to MEG and fMRI. We demonstrate that both approaches excel in their native modality (TDLM for MEG and SODA for fMRI), with comparable effect sizes given idealized conditions in this benchmark. Cross-modality transfer remains challenging. Finally, the FASTIMAGES dataset provides data with known and clearly expressed sequences and can be used to benchmark and validate future sequence detection methods under idealized conditions.

Normative modeling quantifies individual deviation from population norms by estimating the conditional mean and variance of brain-derived measures as functions of clinically relevant parameters such as age. The rapid growth of multicenter consortia has created an urgent need for normative models that incorporate batch harmonization. Several harmonization methods based on linear mixed models--ComBat, GAMLSS, HBR, and Generalized Normative Modeling (GNM)--offer explicit formulations of the mean and variance, making them natural candidates for batch-harmonized normative modeling; yet the absence of a unified theoretical framework leaves it unclear whether and how these methods support the computation of batch-harmonized z-scores. We bridge this gap by writing existing harmonization methods as special cases of a single location-scale equation, y = m(x, {Theta})+{sigma}(x, {Theta}){varepsilon} , which we term the unified form of batch harmonization equation for normative modeling. The methods differ only in the functional forms of m and{sigma} , how batch parameters enter{Theta} , and how{Theta} is estimated. This unified form yields both harmonized data y* and site-invariant z-scores from the same model, providing a common theoretical language for harmonized normative modeling. Building on this framework, we evaluate the underlying regression engines (parametric, spline, Gaussian process, kernel, deep learning), sensitivity to outliers, computational scalability, and federated decomposability for privacy-preserving multi-center computation. By clarifying what each method assumes, what it delivers, and where the boundaries of current methodology lie, the unified equation establishes a principled foundation for method selection and charts a path toward reliable, scalable, and privacy-aware normative modeling across multi-center neuroimaging.