Restoring statistical validity in group analyses of motion-corrupted MRI data

The Neurobiology of Life Course Socioeconomic Conditions and Associated Cognitive Performance in Middle to Late Adulthood

Despite major advances, our understanding of the neurobiology of life course socioeconomic conditions is still scarce. This study aimed to provide insight into the pathways linking socioeconomic exposures-household income, last known occupational position, and life course socioeconomic trajectories-with brain microstructure and cognitive performance in middle to late adulthood. We assessed socioeconomic conditions alongside quantitative relaxometry and diffusion-weighted magnetic resonance imaging indicators of brain tissue microstructure and cognitive performance in a sample of community-dwelling men and women (N = 751, aged 50-91 years). We adjusted the applied regression analyses and structural equation models for the linear and nonlinear effects of age, sex, education, cardiovascular risk factors, and the presence of depression, anxiety, and substance use disorders. Individuals from lower-income households showed signs of advanced brain white matter (WM) aging with greater mean diffusivity (MD), lower neurite density, lower myelination, and lower iron content. The association between household income and MD was mediated by neurite density (B = 0.084, p = 0.003) and myelination (B = 0.019, p = 0.009); MD partially mediated the association between household income and cognitive performance (B = 0.017, p < 0.05). Household income moderated the relation between WM microstructure and cognitive performance, such that greater MD, lower myelination, or lower neurite density was only associated with poorer cognitive performance among individuals from lower-income households. Individuals from higher-income households showed preserved cognitive performance even with greater MD, lower myelination, or lower neurite density. These findings provide novel mechanistic insights into the associations between socioeconomic conditions, brain anatomy, and cognitive performance in middle to late adulthood.

A characterization of cardiac-induced noise in R2 * maps of the brain

PURPOSE

Cardiac pulsation increases the noise level in brain maps of the transverse relaxation rate R2 *. Cardiac-induced noise is challenging to mitigate during the acquisition of R2 * mapping data because its characteristics are unknown. In this work, we aim to characterize cardiac-induced noise in brain maps of the MRI parameter R2 *.

METHODS

We designed a sampling strategy to acquire multi-echo 3D data in 12 intervals of the cardiac cycle, monitored with a fingertip pulse-oximeter. We measured the amplitude of cardiac-induced noise in this data and assessed the effect of cardiac pulsation on R2 * maps computed across echoes. The area of k-space that contains most of the cardiac-induced noise in R2 * maps was then identified. Based on these characteristics, we introduced a tentative sampling strategy that aims to mitigate cardiac-induced noise in R2 * maps of the brain.

RESULTS

In inferior brain regions, cardiac pulsation accounts for R2 * variations of up to 3 s-1 across the cardiac cycle (i.e., ∼35% of the overall variability). Cardiac-induced fluctuations occur throughout the cardiac cycle, with a reduced intensity during the first quarter of the cycle. A total of 50% to 60% of the overall cardiac-induced noise is localized near the k-space center (k < 0.074 mm-1 ). The tentative cardiac noise mitigation strategy reduced the variability of R2 * maps across repetitions by 11% in the brainstem and 6% across the whole brain.

CONCLUSION

We provide a characterization of cardiac-induced noise in brain R2 * maps that can be used as a basis for the design of mitigation strategies during data acquisition.

Statistical analyses of motion-corrupted MRI relaxometry data computed from multiple scans

BACKGROUND

Consistent noise variance across data points (i.e. homoscedasticity) is required to ensure the validity of statistical analyses of MRI data conducted using linear regression methods. However, head motion leads to degradation of image quality, introducing noise heteroscedasticity into ordinary-least square analyses.

NEW METHOD

The recently introduced QUIQI method restores noise homoscedasticity by means of weighted least square analyses in which the weights, specific for each dataset of an analysis, are computed from an index of motion-induced image quality degradation. QUIQI was first demonstrated in the context of brain maps of the MRI parameter R2 * , which were computed from a single set of images with variable echo time. Here, we extend this framework to quantitative maps of the MRI parameters R1, R2 * , and MTsat, computed from multiple sets of images.

RESULTS

QUIQI restores homoscedasticity in analyses of quantitative MRI data computed from multiple scans. QUIQI allows for optimization of the noise model by using metrics quantifying heteroscedasticity and free energy.

COMPARISON WITH EXISTING METHODS

QUIQI restores homoscedasticity more effectively than insertion of an image quality index in the analysis design and yields higher sensitivity than simply removing the datasets most corrupted by head motion from the analysis.

CONCLUSION

QUIQI provides an optimal approach to group-wise analyses of a range of quantitative MRI parameter maps that is robust to inherent homoscedasticity.

[Usefulness of Voxel-Based Quantification (VBQ) Smoothing in Relaxation Time Mapping]

PURPOSE

Voxel-based quantification (VBQ) smoothing is a technique used to smooth quantitative parametric maps in the Montreal Neurological Institute standard space. Although VBQ smoothing could suppress changes in quantitative values at tissue boundaries, its effectiveness on relaxation time (T1 and T2 values and proton density PD) maps has not been investigated. The purpose of this study was to clarify the usefulness of VBQ smoothing in relaxation time mapping.

METHOD

T1 and T2 values and PD maps of the brains of 20 healthy participants were obtained using a two-dimensional multi-dynamic multi-echo sequence. VBQ and Gaussian smoothing were applied to the relaxation time maps by varying the kernel size by 1 mm from 1 to 6 mm. Changes in relaxation time before and after VBQ and Gaussian smoothing for the putamen, caudate nucleus, substantia nigra, and corpus callosum on the relaxation time maps were evaluated.

RESULT

The changes in relaxation time after VBQ smoothing application were smaller than those in that after Gaussian smoothing application. Although the differences in the relaxation time for all tissues before and after VBQ and Gaussian smoothing applications increased with increasing kernel size for all relaxation times for both methods, the changes in the relaxation time for VBQ smoothing were smaller than those in that for Gaussian smoothing.

CONCLUSION

VBQ smoothing can suppress the change in the relaxation time on the boundary of the tissue and is thus a useful smoothing technique in relaxation time mapping.

Arterial hypertension and β-amyloid accumulation have spatially overlapping effects on posterior white matter hyperintensity volume: a cross-sectional study

BACKGROUND

White matter hyperintensities (WMH) in subjects across the Alzheimer's disease (AD) spectrum with minimal vascular pathology suggests that amyloid pathology-not just arterial hypertension-impacts WMH, which in turn adversely influences cognition. Here we seek to determine the effect of both hypertension and Aβ positivity on WMH, and their impact on cognition.

METHODS

We analysed data from subjects with a low vascular profile and normal cognition (NC), subjective cognitive decline (SCD), and amnestic mild cognitive impairment (MCI) enrolled in the ongoing observational multicentre DZNE Longitudinal Cognitive Impairment and Dementia Study (n = 375, median age 70.0 [IQR 66.0, 74.4] years; 178 female; NC/SCD/MCI 127/162/86). All subjects underwent a rich neuropsychological assessment. We focused on baseline memory and executive function-derived from multiple neuropsychological tests using confirmatory factor analysis-, baseline preclinical Alzheimer's cognitive composite 5 (PACC5) scores, and changes in PACC5 scores over the course of three years (ΔPACC5).

RESULTS

Subjects with hypertension or Aβ positivity presented the largest WMH volumes (p_FDR < 0.05), with spatial overlap in the frontal (hypertension: 0.42 ± 0.17; Aβ: 0.46 ± 0.18), occipital (hypertension: 0.50 ± 0.16; Aβ: 0.50 ± 0.16), parietal lobes (hypertension: 0.57 ± 0.18; Aβ: 0.56 ± 0.20), corona radiata (hypertension: 0.45 ± 0.17; Aβ: 0.40 ± 0.13), optic radiation (hypertension: 0.39 ± 0.18; Aβ: 0.74 ± 0.19), and splenium of the corpus callosum (hypertension: 0.36 ± 0.12; Aβ: 0.28 ± 0.12). Elevated global and regional WMH volumes coincided with worse cognitive performance at baseline and over 3 years (p_FDR < 0.05). Aβ positivity was negatively associated with cognitive performance (direct effect-memory: - 0.33 ± 0.08, p_FDR < 0.001; executive: - 0.21 ± 0.08, p_FDR < 0.001; PACC5: - 0.29 ± 0.09, p_FDR = 0.006; ΔPACC5: - 0.34 ± 0.04, p_FDR < 0.05). Splenial WMH mediated the relationship between hypertension and cognitive performance (indirect-only effect-memory: - 0.05 ± 0.02, p_FDR = 0.029; executive: - 0.04 ± 0.02, p_FDR = 0.067; PACC5: - 0.05 ± 0.02, p_FDR = 0.030; ΔPACC5: - 0.09 ± 0.03, p_FDR = 0.043) and WMH in the optic radiation partially mediated that between Aβ positivity and memory (indirect effect-memory: - 0.05 ± 0.02, p_FDR = 0.029).

CONCLUSIONS

Posterior white matter is susceptible to hypertension and Aβ accumulation. Posterior WMH mediate the association between these pathologies and cognitive dysfunction, making them a promising target to tackle the downstream damage related to the potentially interacting and potentiating effects of the two pathologies.

TRIAL REGISTRATION

German Clinical Trials Register (DRKS00007966, 04/05/2015).

Parkinson's disease may disrupt overlapping subthalamic nucleus and pallidal motor networks

There is an ongoing debate about differential clinical outcome and associated adverse effects of deep brain stimulation (DBS) in Parkinson's disease (PD) targeting the subthalamic nucleus (STN) or the globus pallidus pars interna (GPi). Given that functional connectivity profiles suggest beneficial DBS effects within a common network, the empirical evidence about the underlying anatomical circuitry is still scarce. Therefore, we investigate the STN and GPi-associated structural covariance brain patterns in PD patients and healthy controls. We estimate GPi's and STN's whole-brain structural covariance from magnetic resonance imaging (MRI) in a normative mid- to old-age community-dwelling cohort (n = 1184) across maps of grey matter volume, magnetization transfer (MT) saturation, longitudinal relaxation rate (R1), effective transversal relaxation rate (R2*) and effective proton density (PD*). We compare these with the structural covariance estimates in patients with idiopathic PD (n = 32) followed by validation using a reduced size controls' cohort (n = 32). In the normative data set, we observed overlapping spatially distributed cortical and subcortical covariance patterns across maps confined to basal ganglia, thalamus, motor, and premotor cortical areas. Only the subcortical and midline motor cortical areas were confirmed in the reduced size cohort. These findings contrasted with the absence of structural covariance with cortical areas in the PD cohort. We interpret with caution the differential covariance maps of overlapping STN and GPi networks in patients with PD and healthy controls as correlates of motor network disruption. Our study provides face validity to the proposed extension of the currently existing structural covariance methods based on morphometry features to multiparameter MRI sensitive to brain tissue microstructure.

Topography of associations between cardiovascular risk factors and myelin loss in the ageing human brain

Our knowledge of the mechanisms underlying the vulnerability of the brain's white matter microstructure to cardiovascular risk factors (CVRFs) is still limited. We used a quantitative magnetic resonance imaging (MRI) protocol in a single centre setting to investigate the cross-sectional association between CVRFs and brain tissue properties of white matter tracts in a large community-dwelling cohort (n = 1104, age range 46-87 years). Arterial hypertension was associated with lower myelin and axonal density MRI indices, paralleled by higher extracellular water content. Obesity showed similar associations, though with myelin difference only in male participants. Associations between CVRFs and white matter microstructure were observed predominantly in limbic and prefrontal tracts. Additional genetic, lifestyle and psychiatric factors did not modulate these results, but moderate-to-vigorous physical activity was linked to higher myelin content independently of CVRFs. Our findings complement previously described CVRF-related changes in brain water diffusion properties pointing towards myelin loss and neuroinflammation rather than neurodegeneration.

Abnormal brain iron accumulation in obstructive sleep apnea: A quantitative MRI study in the HypnoLaus cohort

Obstructive sleep apnea syndrome (OSA) may be a risk factor for Alzheimer's disease. One of the hallmarks of Alzheimer's disease is disturbed iron homeostasis leading to abnormal iron deposition in brain tissue. To date, there is no empirical evidence to support the hypothesis of altered brain iron homeostasis in patients with obstructive sleep apnea as well. Data were analysed from 773 participants in the HypnoLaus study (mean age 55.9 ± 10.3 years) who underwent polysomnography and brain MRI. Cross-sectional associations were tested between OSA parameters and the MRI effective transverse relaxation rate (R2*) - indicative of iron content - in 68 grey matter regions, after adjustment for confounders. The group with severe OSA (apnea-hypopnea index ≥30/h) had higher iron levels in the left superior frontal gyrus (F_3,760  = 4.79, p = 0.003), left orbital gyri (F_3,760  = 5.13, p = 0.002), right and left middle temporal gyrus (F_3,760  = 4.41, p = 0.004 and F_3,760  = 13.08, p < 0.001, respectively), left angular gyrus (F_3,760  = 6.29, p = 0.001), left supramarginal gyrus (F_3,760  = 4.98, p = 0.003), and right cuneus (F_3,760  = 7.09, p < 0.001). The parameters of nocturnal hypoxaemia were all consistently associated with higher iron levels. Measures of sleep fragmentation had less consistent associations with iron content. This study provides the first evidence of increased brain iron levels in obstructive sleep apnea. The observed iron changes could reflect underlying neuropathological processes that appear to be driven primarily by hypoxaemic mechanisms.

Assessment of perivascular space filtering methods using a three-dimensional computational model

Growing interest surrounds the assessment of perivascular spaces (PVS) on magnetic resonance imaging (MRI) and their validation as a clinical biomarker of adverse brain health. Nonetheless, the limits of validity of current state-of-the-art segmentation methods are still unclear. Here, we propose an open-source three-dimensional computational framework comprising 3D digital reference objects and evaluate the performance of three PVS filtering methods under various spatiotemporal imaging considerations (including sampling, motion artefacts, and Rician noise). Specifically, we study the performance of the Frangi, Jerman and RORPO filters in enhancing PVS-like structures to facilitate segmentation. Our findings were three-fold. First, as long as voxels are isotropic, RORPO outperforms the other two filters, regardless of imaging quality. Unlike the Frangi and Jerman filters, RORPO's performance does not deteriorate as PVS volume increases. Second, the performance of all "vesselness" filters is heavily influenced by imaging quality, with sampling and motion artefacts being the most damaging for these types of analyses. Third, none of the filters can distinguish PVS from other hyperintense structures (e.g. white matter hyperintensities, stroke lesions, or lacunes) effectively, the area under precision-recall curve dropped substantially (Frangi: from 94.21 [IQR 91.60, 96.16] to 43.76 [IQR 25.19, 63.38]; Jerman: from 94.51 [IQR 91.90, 95.37] to 58.00 [IQR 35.68, 64.87]; RORPO: from 98.72 [IQR 95.37, 98.96] to 71.87 [IQR 57.21, 76.63] without and with other hyperintense structures, respectively). The use of our computational model enables comparing segmentation methods and identifying their advantages and disadvantages, thereby providing means for testing and optimising pipelines for ongoing and future studies.

GnRH replacement rescues cognition in Down syndrome

At the present time, no viable treatment exists for cognitive and olfactory deficits in Down syndrome (DS). We show in a DS model (Ts65Dn mice) that these progressive nonreproductive neurological symptoms closely parallel a postpubertal decrease in hypothalamic as well as extrahypothalamic expression of a master molecule that controls reproduction-gonadotropin-releasing hormone (GnRH)-and appear related to an imbalance in a microRNA-gene network known to regulate GnRH neuron maturation together with altered hippocampal synaptic transmission. Epigenetic, cellular, chemogenetic, and pharmacological interventions that restore physiological GnRH levels abolish olfactory and cognitive defects in Ts65Dn mice, whereas pulsatile GnRH therapy improves cognition and brain connectivity in adult DS patients. GnRH thus plays a crucial role in olfaction and cognition, and pulsatile GnRH therapy holds promise to improve cognitive deficits in DS.

CYP2C19 expression modulates affective functioning and hippocampal subiculum volume-a large single-center community-dwelling cohort study

Given controversial findings of reduced depressive symptom severity and increased hippocampus volume in CYP2C19 poor metabolizers, we sought to provide empirical evidence from a large-scale single-center longitudinal cohort in the community-dwelling adult population-Colaus|PsyCoLaus in Lausanne, Switzerland (n = 4152). We looked for CYP2C19 genotype-related behavioral and brain anatomy patterns using a comprehensive set of psychometry, water diffusion- and relaxometry-based magnetic resonance imaging (MRI) data (BrainLaus, n = 1187). Our statistical models tested for differential associations between poor metabolizer and other metabolizer status with imaging-derived indices of brain volume and tissue properties that explain individuals' current and lifetime mood characteristics. The observed association between CYP2C19 genotype and lifetime affective status showing higher functioning scores in poor metabolizers, was mainly driven by female participants (ß = 3.9, p = 0.010). There was no difference in total hippocampus volume between poor metabolizer and other metabolizer, though there was higher subiculum volume in the right hippocampus of poor metabolizers (ß = 0.03, p_FDRcorrected = 0.036). Our study supports the notion of association between mood phenotype and CYP2C19 genotype, however, finds no evidence for concomitant hippocampus volume differences, with the exception of the right subiculum.

Error quantification in multi-parameter mapping facilitates robust estimation and enhanced group level sensitivity

Multi-Parameter Mapping (MPM) is a comprehensive quantitative neuroimaging protocol that enables estimation of four physical parameters (longitudinal and effective transverse relaxation rates R₁ and R₂^*, proton density PD, and magnetization transfer saturation MT_sat) that are sensitive to microstructural tissue properties such as iron and myelin content. Their capability to reveal microstructural brain differences, however, is tightly bound to controlling random noise and artefacts (e.g. caused by head motion) in the signal. Here, we introduced a method to estimate the local error of PD, R₁, and MT_sat maps that captures both noise and artefacts on a routine basis without requiring additional data. To investigate the method's sensitivity to random noise, we calculated the model-based signal-to-noise ratio (mSNR) and showed in measurements and simulations that it correlated linearly with an experimental raw-image-based SNR map. We found that the mSNR varied with MPM protocols, magnetic field strength (3T vs. 7T) and MPM parameters: it halved from PD to R₁ and decreased from PD to MT_sat by a factor of 3-4. Exploring the artefact-sensitivity of the error maps, we generated robust MPM parameters using two successive acquisitions of each contrast and the acquisition-specific errors to down-weight erroneous regions. The resulting robust MPM parameters showed reduced variability at the group level as compared to their single-repeat or averaged counterparts. The error and mSNR maps may better inform power-calculations by accounting for local data quality variations across measurements. Code to compute the mSNR maps and robustly combined MPM maps is available in the open-source hMRI toolbox.

Test-retest reliability of multi-parametric maps (MPM) of brain microstructure

Multiparameter mapping (MPM) is a quantitative MRI protocol that is promising for studying microstructural brain changes in vivo with high specificity. Reliability values are an important prior knowledge for efficient study design and facilitating replicable findings in development, aging and neuroplasticity research. To explore longitudinal reliability of MPM we acquired the protocol in 31 healthy young subjects twice over a rescan interval of 4 weeks. We assessed the within-subject coefficient of variation (WCV), the between-subject coefficient of variation (BCV), and the intraclass correlation coefficient (ICC). Using these metrics, we investigated the reliability of (semi-) quantitative magnetization transfer saturation (MT_sat), proton density (PD), transversal relaxation (R2*) and longitudinal relaxation (R1). To increase relevance for explorative studies in development and training-induced plasticity, we assess reliability both on local voxel- as well as ROI-level. Finally, we disentangle contributions and interplay of within- and between-subject variability to ICC and assess the optimal degree of spatial smoothing applied to the data. We reveal evidence that voxelwise ICC reliability of MPMs is moderate to good with median values in cortex (subcortical GM): MT: 0.789 (0.447) PD: 0.553 (0.264) R1: 0.555 (0.369) R2*: 0.624 (0.477). The Gaussian smoothing kernel of 2 to 4 mm FWHM resulted in optimal reproducibility. We discuss these findings in the context of longitudinal intervention studies and the application to research designs in neuroimaging field.

In vivo Estimation of Axonal Morphology From Magnetic Resonance Imaging and Electroencephalography Data

Purpose

We present a novel approach that allows the estimation of morphological features of axonal fibers from data acquired in vivo in humans. This approach allows the assessment of white matter microscopic properties non-invasively with improved specificity.

Theory

The proposed approach is based on a biophysical model of Magnetic Resonance Imaging (MRI) data and of axonal conduction velocity estimates obtained with Electroencephalography (EEG). In a white matter tract of interest, these data depend on (1) the distribution of axonal radius [P(r)] and (2) the g-ratio of the individual axons that compose this tract [g(r)]. P(r) is assumed to follow a Gamma distribution with mode and scale parameters, M and θ, and g(r) is described by a power law with parameters α and β.

Methods

MRI and EEG data were recorded from 14 healthy volunteers. MRI data were collected with a 3T scanner. MRI-measured g-ratio maps were computed and sampled along the visual transcallosal tract. EEG data were recorded using a 128-lead system with a visual Poffenberg paradigm. The interhemispheric transfer time and axonal conduction velocity were computed from the EEG current density at the group level. Using the MRI and EEG measures and the proposed model, we estimated morphological properties of axons in the visual transcallosal tract.

Results

The estimated interhemispheric transfer time was 11.72 ± 2.87 ms, leading to an average conduction velocity across subjects of 13.22 ± 1.18 m/s. Out of the 4 free parameters of the proposed model, we estimated θ – the width of the right tail of the axonal radius distribution – and β – the scaling factor of the axonal g-ratio, a measure of fiber myelination. Across subjects, the parameter θ was 0.40 ± 0.07 μm and the parameter β was 0.67 ± 0.02 μm^−α.

Conclusion

The estimates of axonal radius and myelination are consistent with histological findings, illustrating the feasibility of this approach. The proposed method allows the measurement of the distribution of axonal radius and myelination within a white matter tract, opening new avenues for the combined study of brain structure and function, and for in vivo histological studies of the human brain.