Categories
Nevin Manimala Statistics

Precision brain morphometry using cluster scanning

Imaging Neurosci (Camb). 2024 May 20;2:imag-2-00175. doi: 10.1162/imag_a_00175. eCollection 2024.

ABSTRACT

Measurement error limits the statistical power to detect group differences and longitudinal change in structural MRI morphometric measures (e.g., hippocampal volume, prefrontal cortical thickness). Recent advances in scan acceleration enable extremely fast T1-weighted scans (~1 minute) that achieve morphometric errors that are close to the errors in longer traditional scans. As acceleration allows multiple scans to be acquired in rapid succession, it becomes possible to pool estimates to increase measurement precision, a strategy known as “cluster scanning.” Here, we explored brain morphometry using cluster scanning in a test-retest study of 40 individuals (12 younger adults, 18 cognitively unimpaired older adults, and 10 adults diagnosed with mild cognitive impairment or Alzheimer’s Dementia). Morphometric errors from a single compressed sensing (CS) 1.0 mm scan (CS) were, on average, 12% larger than a traditional scan using the Alzheimer’s Disease Neuroimaging Initiative (ADNI) protocol. Pooled estimates from four clustered CS acquisitions led to errors that were 34% smaller than ADNI despite having a shorter total acquisition time. Given a fixed amount of time, a gain in measurement precision can thus be achieved by acquiring multiple rapid scans instead of a single traditional scan. Errors were further reduced when estimates were pooled from eight CS scans (51% smaller than ADNI). Neither pooling across a break nor pooling across multiple scans of different spatial resolutions boosted this benefit. We discuss the potential of cluster scanning to improve morphometric precision, boost statistical power, and produce more sensitive disease progression biomarkers.

PMID:40800474 | PMC:PMC12247576 | DOI:10.1162/imag_a_00175

Categories
Nevin Manimala Statistics

Balance of power: The choice between trial and participant numbers to optimise the detection of phase-dependent effects

Imaging Neurosci (Camb). 2024 Nov 5;2:imag-2-00345. doi: 10.1162/imag_a_00345. eCollection 2024.

ABSTRACT

The fields of neuroscience and psychology are currently in the midst of aso-called reproducibility crisis, with growing concerns regarding a history ofweak effect sizes and low statistical power in much of the research published inthese fields over the last few decades. Whilst the traditional approach foraddressing this criticism has been to increaseparticipantsample sizes, there are many research contexts in which the number oftrialsper participant may be of equal importance. Thepresent study aimed to compare the relative importance of participants andtrials in the detection of phase-dependent phenomena, which are measured acrossa range of neuroscientific contexts (e.g., neural oscillations, non-invasivebrain stimulation). This was achievable within a simulated environment where onecan manipulate the strength of this phase dependency in two types of outcomevariables: one with normally distributed residuals (idealistic) and onecomparable with motor-evoked potentials (an MEP-like variable). We compared thestatistical power across thousands of experiments with the same number ofsessions per experiment but with different proportions of participants andnumber of sessions per participant (30 participants × 1 session, 15participants × 2 sessions, and 10 participants × 3 sessions), withthe trials being pooled across sessions for each participant. These simulationswere performed for both outcome variables (idealistic and MEP-like) and fourdifferent effect sizes (0.075-“weak,”0.1-“moderate,” 0.125-“strong,”0.15-“very strong”), as well as separate control scenarioswith no true effect. Across all scenarios with (true) discoverable effects, andfor both outcome types, there was a statistical benefit for experimentsmaximising the number of trials rather than the number of participants (i.e., itwas always beneficial to recruit fewer participants but have them complete moretrials). These findings emphasise the importance of obtaining sufficientindividual-level data rather than simply increasing number of participants.

PMID:40800470 | PMC:PMC12290742 | DOI:10.1162/imag_a_00345

Categories
Nevin Manimala Statistics

Spatial-extent inference for testing variance components in reliability and heritability studies

Imaging Neurosci (Camb). 2024 Jan 9;2:imag-2-00058. doi: 10.1162/imag_a_00058. eCollection 2024.

ABSTRACT

Clusterwise inference is a popular approach in neuroimaging to increase sensitivity, but most existing methods are currently restricted to the General Linear Model (GLM) for testing mean parameters. Statistical methods for testing variance components, which are critical in neuroimaging studies that involve estimation of narrow-sense heritability or test-retest reliability, are underdeveloped due to methodological and computational challenges, which would potentially lead to low power. We propose a fast and powerful test for variance components called CLEAN-V (CLEAN for testing Variance components). CLEAN-V models the global spatial dependence structure of imaging data and computes a locally powerful variance component test statistic by data-adaptively pooling neighborhood information. Correction for multiple comparisons is achieved by permutations to control family-wise error rate (FWER). Through analysis of task-functional magnetic resonance imaging (fMRI) data from the Human Connectome Project across five tasks and comprehensive data-driven simulations, we show that CLEAN-V outperforms existing methods in detecting test-retest reliability and narrow-sense heritability with significantly improved power, with the detected areas aligning with activation maps. The computational efficiency of CLEAN-V also speaks of its practical utility, and it is available as an R package.

PMID:40800456 | PMC:PMC12224426 | DOI:10.1162/imag_a_00058

Categories
Nevin Manimala Statistics

Superpixel-ComBat modeling: A joint approach for harmonization and characterization of inter-scanner variability in T1-weighted images

Imaging Neurosci (Camb). 2024 Oct 3;2:imag-2-00306. doi: 10.1162/imag_a_00306. eCollection 2024.

ABSTRACT

T1-weighted imaging holds wide applications in clinical and research settings; however, the challenge of inter-scanner variability arises when combining data across scanners, which impedes multi-site research. To address this, post-acquisition harmonization methods such as statistical or deep learning approaches have been proposed to unify cross-scanner images. Nevertheless, how inter-scanner variability manifests in images and derived measures, and how to harmonize it in an interpretable manner, remains underexplored. To broaden our knowledge of inter-scanner variability and leverage it to develop a new harmonization strategy, we devised a pipeline to assess the interpretable inter-scanner variability in matched T1-weighted images across four 3T MRI scanners. The pipeline incorporates ComBat modeling with 3D superpixel parcellation algorithm (namely SP-ComBat), which estimates location and scale effects to quantify the shift and spread in relative signal distributions, respectively, concerning brain tissues in the image domain. The estimated parametric maps revealed significant contrast deviations compared to the joint signal distribution across scanners (p< 0.001), and the identified deviations in signal intensities may relate to differences in the inversion time acquisition parameter. To reduce the inter-scanner variability, we implemented a harmonization strategy involving proper image preprocessing and site effect removal by ComBat-derived parameters, achieving substantial improvement in image quality and significant reduction in variation of volumetric measures of brain tissues (p< 0.001). We also applied SP-ComBat to evaluate and characterize the performance of various image harmonization techniques, demonstrating a new way to assess image harmonization. In addition, we reported various metrics of T1-weighted images to quantify the impact of inter-scanner variation, including signal-to-noise ratio, contrast-to-noise ratio, signal inhomogeneity index, and structural similarity index. This study demonstrates a pipeline that extends the implementation of statistical ComBat method to the image domain in a practical manner for characterizing and harmonizing the inter-scanner variability in T1-weighted images, providing further insight for the studies focusing on the development of image harmonization methodologies and their applications.

PMID:40800451 | PMC:PMC12290534 | DOI:10.1162/imag_a_00306

Categories
Nevin Manimala Statistics

PET imaging of the serotonin 1A receptor in major depressive disorder: Hierarchical multivariate analysis of [ 11 C]WAY100635 overcomes outcome measure discrepancies

Imaging Neurosci (Camb). 2024 Oct 25;2:imag-2-00328. doi: 10.1162/imag_a_00328. eCollection 2024.

ABSTRACT

The serotonin 1A receptor has been linked to both the pathophysiology of major depressive disorder (MDD) and the antidepressant action of serotonin reuptake inhibitors. Most PET studies of the serotonin 1A receptor in MDD used the receptor antagonist radioligand, [carbonyl- C 11 ]WAY100635; however, the interpretation of the combined results has been contentious owing to reports of higher or lower binding in MDD with different outcome measures. The reasons for these divergent results originate from several sources, including properties of the radiotracer itself, which complicate its quantification and interpretation; as well as from previously reported differences between MDD and healthy volunteers in both reference tissue binding and plasma-free fraction, which are typically assumed not to differ. Recently, we have developed two novel hierarchical multivariate methods which we validated for the quantification and analysis of [ C 11 ]WAY100635, which show better accuracy and inferential efficiency compared to standard analysis approaches. Importantly, these new methods should theoretically be more resilient to many of the factors thought to have caused the discrepancies observed in previous studies. We sought to apply these methods in the largest [ C 11 ]WAY100635 sample to date, consisting of 160 individuals, including 103 MDD patients, of whom 50 were not-recently-medicated and 53 were antidepressant-exposed, as well as 57 healthy volunteers. While the outcome measure discrepancies were substantial using conventional univariate analysis, our multivariate analysis techniques instead yielded highly consistent results across PET outcome measures and across pharmacokinetic models, with all approaches showing higher serotonin 1A autoreceptor binding potential in the raphe nuclei of not-recently-medicated MDD patients relative to both healthy volunteers and antidepressant-exposed MDD patients. Moreover, with the additional precision of estimates afforded by this approach, we can show that while binding is also higher in projection areas in this group, these group differences are approximately half of those in the raphe nuclei, which are statistically distinguishable from one another. These results are consistent with the biological role of the serotonin 1A autoreceptor in the raphe nuclei in regulating serotonin neuron firing and release, and with preclinical and clinical evidence of deficient serotonin activity in MDD due to over-expression of autoreceptors resulting from genetic and/or epigenetic effects. These results are also consistent with downregulation of autoreceptors as a mechanism of action of selective serotonin reuptake inhibitors. In summary, the results using multivariate analysis approaches, therefore, demonstrate both face and convergent validity, and may serve to provide a resolution and consensus interpretation for the disparate results of previous studies examining the serotonin 1A receptor in MDD.

PMID:40800447 | PMC:PMC12290790 | DOI:10.1162/imag_a_00328

Categories
Nevin Manimala Statistics

Unconstrained quantitative magnetization transfer imaging: Disentangling T 1 of the free and semi-solid spin pools

Imaging Neurosci (Camb). 2024 May 20;2:imag-2-00177. doi: 10.1162/imag_a_00177. eCollection 2024.

ABSTRACT

Since the inception of magnetization transfer (MT) imaging, it has been widely assumed that Henkelman’s two spin pools have similar longitudinal relaxation times, which motivated many researchers to constrain them to each other. However, several recent publications reported a T 1 s of thesemi-solid spin poolthat is much shorter than T 1 f of thefree pool. While these studies tailored experiments for robust proofs-of-concept, we here aim to quantify the disentangled relaxation processes on a voxel-by-voxel basis in a clinical imaging setting, that is, with an effective resolution of 1.24mm isotropic and full brain coverage in 12min. To this end, we optimized ahybrid-statepulse sequence for mapping the parameters of an unconstrained MT model. We scanned four people with relapsing-remitting multiple sclerosis (MS) and four healthy controls with this pulse sequence and estimated T 1 f 1.84 s and T 1 s 0.34 s in healthy white matter. Our results confirm the reports that T 1 s T 1 f and we argue that this finding identifies MT as an inherent driver of longitudinal relaxation in brain tissue. Moreover, we estimated a fractional size of the semi-solid spin pool of m 0 s 0.212 , which is larger than previously assumed. An analysis of T 1 f in normal-appearing white matter revealed statistically significant differences between individuals with MS and controls.

PMID:40800438 | PMC:PMC12247553 | DOI:10.1162/imag_a_00177

Categories
Nevin Manimala Statistics

Rescuing missing data in connectome-based predictive modeling

Imaging Neurosci (Camb). 2024 Feb 2;2:imag-2-00071. doi: 10.1162/imag_a_00071. eCollection 2024.

ABSTRACT

Recent evidence suggests brain-phenotype predictions may require very large sample sizes. However, as the sample size increases, missing data also increase. Conventional methods, like complete-case analysis, discard useful information and shrink the sample size. To address the missing data problem, we investigated rescuing these missing data through imputation. Imputation is substituting estimated values for missing data to be used in downstream analyses. We integrated imputation methods into the Connectome-based Predictive Modeling (CPM) framework. Utilizing four open-source datasets-the Human Connectome Project, the Philadelphia Neurodevelopmental Cohort, the UCLA Consortium for Neuropsychiatric Phenomics, and the Healthy Brain Network (HBN)-we validated and compared our framework with different imputation methods against complete-case analysis for both missing connectomes and missing phenotypic measures scenarios. Imputing connectomes exhibited superior prediction performance on real and simulated missing data compared to complete-case analysis. In addition, we found that imputation accuracy was a good indicator for choosing an imputation method for missing phenotypic measures but not informative for missing connectomes. In a real-world example predicting cognition using the HBN, we rescued 628 individuals through imputation, doubling the complete case sample size and increasing the variance explained by the predicted value by 45%. In conclusion, our study is a benchmark for state-of-the-art imputation techniques when dealing with missing connectome and phenotypic data in predictive modeling scenarios. Our results suggest that improving prediction performance can be achieved by strategically addressing missing data through effective imputation methods rather than resorting to the outright exclusion of participants. Our results suggest that rescuing data with imputation, instead of discarding participants with missing information, improves prediction performance.

PMID:40800425 | PMC:PMC12224408 | DOI:10.1162/imag_a_00071

Categories
Nevin Manimala Statistics

Motion-invariant variational autoencoding of brain structural connectomes

Imaging Neurosci (Camb). 2024 Oct 7;2:imag-2-00303. doi: 10.1162/imag_a_00303. eCollection 2024.

ABSTRACT

Mapping of human brain structural connectomes via diffusion magnetic resonance imaging (dMRI) offers a unique opportunity to understand brain structural connectivity and relate it to various human traits, such as cognition. However, head displacement during image acquisition can compromise the accuracy of connectome reconstructions and subsequent inference results. We develop a generative model to learn low-dimensional representations of structural connectomes invariant to motion-induced artifacts, so that we can link brain networks and human traits more accurately, and generate motion-adjusted connectomes. We apply the proposed model to data from the Adolescent Brain Cognitive Development (ABCD) study and the Human Connectome Project (HCP) to investigate how our motion-invariant connectomes facilitate understanding of the brain network and its relationship with cognition. Empirical results demonstrate that the proposed motion-invariant variational autoencoder (inv-VAE) outperforms its competitors in various aspects. In particular, motion-adjusted structural connectomes are more strongly associated with a wide array of cognition-related traits than other approaches without motion adjustment.

PMID:40800413 | PMC:PMC12290590 | DOI:10.1162/imag_a_00303

Categories
Nevin Manimala Statistics

Progression of white matter hyperintensities is related to blood pressure increases and global cognitive decline – A registered report

Imaging Neurosci (Camb). 2024 Jun 24;2:imag-2-00188. doi: 10.1162/imag_a_00188. eCollection 2024.

ABSTRACT

White matter hyperintensities (WMH) reflect cerebral small vessel disease (cSVD), a major brain pathology contributing to cognitive decline and dementia. Vascular risk factors, including higher diastolic blood pressure (DBP), have been associated with the progression of WMH yet longitudinal studies have not comprehensively assessed these effects for abdominal obesity or reported sex/gender-specific effects. In this pre-registered analysis of a longitudinal population-based neuroimaging cohort, we investigated the association of baseline DBP and waist-to-hip ratio with WMH progression in linear mixed models. We also examined the relationship of WMH progression and executive and global cognitive function. We conducted gender interaction and stratified analyses. We included data from 596 individuals (44.1 % females, mean age = 63.2 years) with two MRI scans over approximately 6 years. We did not find a significant association of baseline DBP with WMH progression. WMH progression significantly predicted global cognitive decline but not decline in executive function. In exploratory analyses, increases in DBP as well as baseline and increase in systolic blood pressure were associated with WMH progression, confined to frontal periventricular regions. There was no association of WHR nor any gender-specific associations with WMH progression. Adequate BP control might contribute to limit WMH progression and negative effects on global cognitive function in the middle-aged to older population for men and women.

PMID:40800400 | PMC:PMC12272209 | DOI:10.1162/imag_a_00188

Categories
Nevin Manimala Statistics

VertexWiseR: A package for simplified vertex-wise analyses of whole-brain and hippocampal surfaces in R

Imaging Neurosci (Camb). 2024 Nov 14;2:imag-2-00372. doi: 10.1162/imag_a_00372. eCollection 2024.

ABSTRACT

Currently, whole-brain vertex-wise analyses on brain surfaces commonly require specially configured operating systems/environments to run and are largely inaccessible to R users. As such, these analyses are inconvenient to execute and inaccessible to many aspiring researchers. To address these limitations, we present VertexWiseR, a user-friendly R package, to run cortical and hippocampal surface vertex-wise analyses, in just about any computer, requiring minimal technical expertise and computational resources. The package allows cohort-wise anatomical surface data to be highly compressed into a single, compact, easy-to-share file. Users can then run a range of vertex-wise statistical analyses with that single file without requiring a special operating system/environment and direct access to the preprocessed file directories. This enables the user to easily take the analyses “offline”, which would be highly appropriate and conducive in classroom settings. This R package includes a conventional suite of tools for extracting, manipulating, analyzing, and visualizing vertex-wise data, and is designed to be easy for beginners to use. Furthermore, it also contains novel or advanced functionalities such as hippocampal surface analyses, meta-analytic decoding, threshold-free cluster enhancement, and mixed-effects models that would appeal to experienced researchers as well. In the current report, we showcase these functionalities in the analyses of two publicly accessible datasets. Overall, our R package opens up new frontiers for the R’s user base/community and makes such neuroimaging analyses accessible to the masses.

PMID:40800380 | PMC:PMC12330379 | DOI:10.1162/imag_a_00372