NRN1 epistasis with BDNF and CACNA1C: mediation effects on symptom severity through neuroanatomical changes in schizophrenia.

The expression of Neuritin-1 (NRN1), a neurotrophic factor crucial for neurodevelopment and synaptic plasticity, is enhanced by the Brain Derived Neurotrophic Factor (BDNF). Although the receptor of NRN1 remains unclear, it is suggested that NRN1's activation of the insulin receptor (IR) pathway promotes the transcription of the calcium voltage-gated channel subunit alpha1 C (CACNA1C). These three genes have been independently associated with schizophrenia (SZ) risk, symptomatology, and brain differences. However, research on how they synergistically modulate these phenotypes is scarce. We aimed to study whether the genetic epistasis between these genes affects the risk and clinical presentation of the disorder via its effect on brain structure. First, we tested the epistatic effect of NRN1 and BDNF or CACNA1C on (i) the risk for SZ, (ii) clinical symptoms severity and functionality (onset, PANSS, CGI and GAF), and (iii) brain cortical structure (thickness, surface area and volume measures estimated using FreeSurfer) in a sample of 86 SZ patients and 89 healthy subjects. Second, we explored whether those brain clusters influenced by epistatic effects mediate the clinical profiles. Although we did not find a direct epistatic impact on the risk, our data unveiled significant effects on the disorder's clinical presentation. Specifically, the NRN1-rs10484320 x BDNF-rs6265 interplay influenced PANSS general psychopathology, and the NRN1-rs4960155 x CACNA1C-rs1006737 interaction affected GAF scores. Moreover, several interactions between NRN1 SNPs and BDNF-rs6265 significantly influenced the surface area and cortical volume of the frontal, parietal, and temporal brain regions within patients. The NRN1-rs10484320 x BDNF-rs6265 epistasis in the left lateral orbitofrontal cortex fully mediated the effect on PANSS general psychopathology. Our study not only adds clinical significance to the well-described molecular relationship between NRN1 and BDNF but also underscores the utility of deconstructing SZ into biologically validated brain-imaging markers to explore their mediation role in the path from genetics to complex clinical manifestation.

Improved estimation of the risk of manic relapse by combining clinical and brain scan data.

INTRODUCTION: Estimating the risk of manic relapse could help the psychiatrist individually adjust the treatment to the risk. Some authors have attempted to estimate this risk from baseline clinical data. Still, no studies have assessed whether the estimation could improve by adding structural magnetic resonance imaging (MRI) data. We aimed to evaluate it. MATERIAL AND METHODS: We followed a cohort of 78 patients with a manic episode without mixed symptoms (bipolar type I or schizoaffective disorder) at 2-4-6-9-12-15-18 months and up to 10 years. Within a cross-validation scheme, we created and evaluated a Cox lasso model to estimate the risk of manic relapse using both clinical and MRI data. RESULTS: The model successfully estimated the risk of manic relapse (Cox regression of the time to relapse as a function of the estimated risk: hazard ratio (HR)=2.35, p=0.027; area under the curve (AUC)=0.65, expected calibration error (ECE)<0.2). The most relevant variables included in the model were the diagnosis of schizoaffective disorder, poor impulse control, unusual thought content, and cerebellum volume decrease. The estimations were poorer when we used clinical or MRI data separately. CONCLUSION: Combining clinical and MRI data may improve the risk of manic relapse estimation after a manic episode. We provide a website that estimates the risk according to the model to facilitate replication by independent groups before translation to clinical settings.

Widespread intra-axonal signal fraction abnormalities in bipolar disorder from multicompartment diffusion MRI: Sensitivity to diagnosis, association with clinical features and pharmacologic treatment.

Despite diffusion tensor imaging (DTI) evidence for widespread fractional anisotropy (FA) reductions in the brain white matter of patients with bipolar disorder, questions remain regarding the specificity and sensitivity of FA abnormalities as opposed to other diffusion metrics in the disorder. We conducted a whole-brain voxel-based multicompartment diffusion MRI study on 316 participants (i.e., 158 patients and 158 matched healthy controls) employing four diffusion metrics: the mean diffusivity (MD) and FA estimated from DTI, and the intra-axonal signal fraction (IASF) and microscopic axonal parallel diffusivity (Dpar) derived from the spherical mean technique. Our findings provide novel evidence about widespread abnormalities in other diffusion metrics in BD. An extensive overlap between the FA and IASF results suggests that the lower FA in patients may be caused by a reduced intra-axonal volume fraction or a higher macromolecular content in the intra-axonal water. We also found a diffuse alteration in MD involving white and grey matter tissue and more localised changes in Dpar. A Machine Learning analysis revealed that FA, followed by IASF, were the most helpful metric for the automatic diagnosis of BD patients, reaching an accuracy of 72%. Number of mood episodes, age of onset/duration of illness, psychotic symptoms, and current treatment with lithium, antipsychotics, antidepressants, and antiepileptics were all significantly associated with microstructure abnormalities. Lithium treatment was associated with less microstructure abnormality.

Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI.

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.

Text as signal. A tutorial with case studies focusing on social media (Twitter).

Sentiment analysis is the automated coding of emotions expressed in text. Sentiment analysis and other types of analyses focusing on the automatic coding of textual documents are increasingly popular in psychology and computer science. However, the potential of treating automatically coded text collected with regular sampling intervals as a signal is currently overlooked. We use the phrase "text as signal" to refer to the application of signal processing techniques to coded textual documents sampled with regularity. In order to illustrate the potential of treating text as signal, we introduce the reader to a variety of such techniques in a tutorial with two case studies in the realm of social media analysis. First, we apply finite response impulse filtering to emotion-coded tweets posted during the US Election Week of 2020 and discuss the visualization of the resulting variation in the filtered signal. We use changepoint detection to highlight the important changes in the emotional signals. Then we examine data interpolation, analysis of periodicity via the fast Fourier transform (FFT), and FFT filtering to personal value-coded tweets from November 2019 to October 2020 and link the variation in the filtered signal to some of the epoch-defining events occurring during this period. Finally, we use block bootstrapping to estimate the variability/uncertainty in the resulting filtered signals. After working through the tutorial, the readers will understand the basics of signal processing to analyze regularly sampled coded text.

Axonal T2 estimation using the spherical variance of the strongly diffusion-weighted MRI signal.

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T2 that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T2, as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T2 based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T2 estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system.

Data-driven myelin water imaging based on T1 and T2 relaxometry.

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T1 and T2 relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T1 /T2 data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

Revisiting the T2 spectrum imaging inverse problem: Bayesian regularized non-negative least squares.

Multi-echo T2 magnetic resonance images contain information about the distribution of T2 relaxation times of compartmentalized water, from which we can estimate relevant brain tissue properties such as the myelin water fraction (MWF). Regularized non-negative least squares (NNLS) is the tool of choice for estimating non-parametric T2 spectra. However, the estimation is ill-conditioned, sensitive to noise, and highly affected by the employed regularization weight. The purpose of this study is threefold: first, we want to underline that the apparently innocuous use of two alternative parameterizations for solving the inverse problem, which we called the standard and alternative regularization forms, leads to different solutions; second, to assess the performance of both parameterizations; and third, to propose a new Bayesian regularized NNLS method (BayesReg). The performance of BayesReg was compared with that of two conventional approaches (L-curve and Chi-square (X2) fitting) using both regularization forms. We generated a large dataset of synthetic data, acquired in vivo human brain data in healthy participants for conducting a scan-rescan analysis, and correlated the myelin content derived from histology with the MWF estimated from ex vivo data. Results from synthetic data indicate that BayesReg provides accurate MWF estimates, comparable to those from L-curve and X2, and with better overall stability across a wider signal-to-noise range. Notably, we obtained superior results by using the alternative regularization form. The correlations reported in this study are higher than those reported in previous studies employing the same ex vivo and histological data. In human brain data, the estimated maps from L-curve and BayesReg were more reproducible. However, the T2 spectra produced by BayesReg were less affected by over-smoothing than those from L-curve. These findings suggest that BayesReg is a good alternative for estimating T2 distributions and MWF maps.

Age- and gender-related differences in brain tissue microstructure revealed by multi-component T2 relaxometry.

In spite of extensive work, inconsistent findings and lack of specificity in most neuroimaging techniques used to examine age- and gender-related patterns in brain tissue microstructure indicate the need for additional research. Here, we performed the largest Multi-component T2 relaxometry cross-sectional study to date in healthy adults (N = 145, 18-60 years). Five quantitative microstructure parameters derived from various segments of the estimated T2 spectra were evaluated, allowing a more specific interpretation of results in terms of tissue microstructure. We found similar age-related myelin water fraction (MWF) patterns in men and women but we also observed differential male related results including increased MWF content in a few white matter tracts, a faster decline with age of the intra- and extra-cellular water fraction and its T2 relaxation time (i.e. steeper age related negative slopes) and a faster increase in the free and quasi-free water fraction, spanning the whole grey matter. Such results point to a sexual dimorphism in brain tissue microstructure and suggest a lesser vulnerability to age-related changes in women.

Comparison of non-parametric T2 relaxometry methods for myelin water quantification.

Multi-component T2 relaxometry allows probing tissue microstructure by assessing compartment-specific T2 relaxation times and water fractions, including the myelin water fraction. Non-negative least squares (NNLS) with zero-order Tikhonov regularization is the conventional method for estimating smooth T2 distributions. Despite the improved estimation provided by this method compared to non-regularized NNLS, the solution is still sensitive to the underlying noise and the regularization weight. This is especially relevant for clinically achievable signal-to-noise ratios. In the literature of inverse problems, various well-established approaches to promote smooth solutions, including first-order and second-order Tikhonov regularization, and different criteria for estimating the regularization weight have been proposed, such as L-curve, Generalized Cross-Validation, and Chi-square residual fitting. However, quantitative comparisons between the available reconstruction methods for computing the T2 distribution, and between different approaches for selecting the optimal regularization weight, are lacking. In this study, we implemented and evaluated ten reconstruction algorithms, resulting from the individual combinations of three penalty terms with three criteria to estimate the regularization weight, plus non-regularized NNLS. Their performance was evaluated both in simulated data and real brain MRI data acquired from healthy volunteers through a scan-rescan repeatability analysis. Our findings demonstrate the need for regularization. As a result of this work, we provide a list of recommendations for selecting the optimal reconstruction algorithms based on the acquired data. Moreover, the implemented methods were packaged in a freely distributed toolbox to promote reproducible research, and to facilitate further research and the use of this promising quantitative technique in clinical practice.

Model-informed machine learning for multi-component T2 relaxometry.

Recovering the T2 distribution from multi-echo T2 magnetic resonance (MR) signals is challenging but has high potential as it provides biomarkers characterizing the tissue micro-structure, such as the myelin water fraction (MWF). In this work, we propose to combine machine learning and aspects of parametric (fitting from the MRI signal using biophysical models) and non-parametric (model-free fitting of the T2 distribution from the signal) approaches to T2 relaxometry in brain tissue by using a multi-layer perceptron (MLP) for the distribution reconstruction. For training our network, we construct an extensive synthetic dataset derived from biophysical models in order to constrain the outputs with a priori knowledge of in vivo distributions. The proposed approach, called Model-Informed Machine Learning (MIML), takes as input the MR signal and directly outputs the associated T2 distribution. We evaluate MIML in comparison to a Gaussian Mixture Fitting (parametric) and Regularized Non-Negative Least Squares algorithms (non-parametric) on synthetic data, an ex vivo scan, and high-resolution scans of healthy subjects and a subject with Multiple Sclerosis. In synthetic data, MIML provides more accurate and noise-robust distributions. In real data, MWF maps derived from MIML exhibit the greatest conformity to anatomical scans, have the highest correlation to a histological map of myelin volume, and the best unambiguous lesion visualization and localization, with superior contrast between lesions and normal appearing tissue. In whole-brain analysis, MIML is 22 to 4980 times faster than the non-parametric and parametric methods, respectively.

Fast and high-resolution myelin water imaging: Accelerating multi-echo GRASE with CAIPIRINHA.

PURPOSE: Although several MRI methods have been explored to achieve in vivo myelin quantification, imaging the whole brain in clinically acceptable times and sufficiently high resolution remains challenging. To address this problem, this work investigates the acceleration of multi-echo T2 acquisitions based on the multi-echo gradient and spin echo (GRASE) sequence using CAIPIRINHA undersampling and adapted k-space reordering patterns. METHODS: A prototype multi-echo GRASE sequence supporting CAIPIRINHA parallel imaging was implemented. Multi-echo T2 data were acquired from 12 volunteers using the implemented sequence (1.6 × 1.6 × 1.6 mm3 , 84 slices, acquisition time [TA] = 10:30 min) and a multi-echo spin echo (MESE) sequence as reference (1.6 × 1.6 × 3.2 mm3 , single-slice, TA = 5:41 min). Myelin water fraction (MWF) maps derived from both acquisitions were compared via correlation and Bland-Altman analyses. In addition, scan-rescan datasets were acquired to evaluate the repeatability of the derived maps. RESULTS: Resulting maps from the MESE and multi-echo GRASE sequences were found to be correlated (r = 0.83). The Bland-Altman analysis revealed a mean bias of -0.2% (P = .24) with the limits of agreement ranging from -3.7% to 3.3%. The Pearson's correlation coefficient among MWF values obtained from the scan-rescan datasets was found to be 0.95 and the mean bias equal to 0.11% (P = .32), indicating good repeatability of the retrieved maps. CONCLUSION: By combining a 3D multi-echo GRASE sequence with CAIPIRINHA sampling, whole-brain MWF maps were obtained in 10:30 min with 1.6 mm isotropic resolution. The good correlation with conventional MESE-based maps demonstrates that the implemented sequence may be a promising alternative to time-consuming MESE acquisitions.

Multivariate Brain Functional Connectivity Through Regularized Estimators.

Functional connectivity analyses are typically based on matrices containing bivariate measures of covariability, such as correlations. Although this has been a fruitful approach, it may not be the optimal strategy to fully explore the complex associations underlying brain activity. Here, we propose extending connectivity to multivariate functions relating to the temporal dynamics of a region with the rest of the brain. The main technical challenges of such an approach are multidimensionality and its associated risk of overfitting or even the non-uniqueness of model solutions. To minimize these risks, and as an alternative to the more common dimensionality reduction methods, we propose using two regularized multivariate connectivity models. On the one hand, simple linear functions of all brain nodes were fitted with ridge regression. On the other hand, a more flexible approach to avoid linearity and additivity assumptions was implemented through random forest regression. Similarities and differences between both methods and with simple averages of bivariate correlations (i.e., weighted global brain connectivity) were evaluated on a resting state sample of N = 173 healthy subjects. Results revealed distinct connectivity patterns from the two proposed methods, which were especially relevant in the age-related analyses where both ridge and random forest regressions showed significant patterns of age-related disconnection, almost completely absent from the much less sensitive global brain connectivity maps. On the other hand, the greater flexibility provided by the random forest algorithm allowed detecting sex-specific differences. The generic framework of multivariate connectivity implemented here may be easily extended to other types of regularized models.

Robust Monte-Carlo Simulations in Diffusion-MRI: Effect of the Substrate Complexity and Parameter Choice on the Reproducibility of Results.

Monte-Carlo Diffusion Simulations (MCDS) have been used extensively as a ground truth tool for the validation of microstructure models for Diffusion-Weighted MRI. However, methodological pitfalls in the design of the biomimicking geometrical configurations and the simulation parameters can lead to approximation biases. Such pitfalls affect the reliability of the estimated signal, as well as its validity and reproducibility as ground truth data. In this work, we first present a set of experiments in order to study three critical pitfalls encountered in the design of MCDS in the literature, namely, the number of simulated particles and time steps, simplifications in the intra-axonal substrate representation, and the impact of the substrate's size on the signal stemming from the extra-axonal space. The results obtained show important changes in the simulated signals and the recovered microstructure features when changes in those parameters are introduced. Thereupon, driven by our findings from the first studies, we outline a general framework able to generate complex substrates. We show the framework's capability to overcome the aforementioned simplifications by generating a complex crossing substrate, which preserves the volume in the crossing area and achieves a high packing density. The results presented in this work, along with the simulator developed, pave the way toward more realistic and reproducible Monte-Carlo simulations for Diffusion-Weighted MRI.

Discoidin domain receptor 1 gene variants are associated with decreased white matter fractional anisotropy and decreased processing speed in schizophrenia.

DDR1 has been linked to schizophrenia (SZ) and myelination. Here, we tested whether DDR1 variants in people at risk for SZ influence white matter (WM) structural variations and cognitive processing speed (PS). First, following a case-control design (Study 1), SZ patients (N = 1193) and controls (N = 1839) were genotyped for rs1264323 and rs2267641 at DDR1, and the frequencies were compared. We replicated the association between DDR1 and SZ (rs1264323, adjusted P = 0.015). Carriers of the rs1264323AA combined with the rs2267641AC or CC genotype are at risk to develop SZ compared to the other genotype combinations. Second, SZ patients (Study 2, N = 194) underwent an evaluation of PS using the Trail Making Test (TMT) and DDR1 genotyping. To compare PS between DDR1 genotype groups, we conducted an analysis of covariance (including rs1264323 as a covariate) and found that SZ patients with the rs2267641CC genotype had decreased PS compared to patients with the AA and AC genotypes. Third, 54 patients (Study 3) from Study 2 were selected based on rs1264323 genotype to undergo reevaluation, including a DTI-MRI brain scan. To test for associations between PS, WM microstructure and DDR1 genotype, we first localized those WM regions where fractional anisotropy (FA) was correlated with PS and tested whether FA showed differences between the rs1264323 genotypes. SZ patients with the rs1264323AA genotype showed decreased FA in WM regions associated with decreased PS. We conclude that DDR1 variants may confer a risk of SZ through WM microstructural alterations leading to cognitive dysfunction.

Sparse wars: A survey and comparative study of spherical deconvolution algorithms for diffusion MRI.

Spherical deconvolution methods are widely used to estimate the brain's white-matter fiber orientations from diffusion MRI data. In this study, eight spherical deconvolution algorithms were implemented and evaluated. These included two model selection techniques based on the extended Bayesian information criterion (i.e., best subset selection and the least absolute shrinkage and selection operator), iteratively reweighted l2- and l1-norm approaches to approximate the l0-norm, sparse Bayesian learning, Cauchy deconvolution, and two accelerated Richardson-Lucy algorithms. Results from our exhaustive evaluation show that there is no single optimal method for all different fiber configurations, suggesting that further studies should be conducted to find the optimal way of combining solutions from different methods. We found l0-norm regularization algorithms to resolve more accurately fiber crossings with small inter-fiber angles. However, in voxels with very dominant fibers, algorithms promoting more sparsity are less accurate in detecting smaller fibers. In most cases, the best algorithm to reconstruct fiber crossings with two fibers did not perform optimally in voxels with one or three fibers. Therefore, simplified validation systems as employed in a number of previous studies, where only two fibers with similar volume fractions were tested, should be avoided as they provide incomplete information. Future studies proposing new reconstruction methods based on high angular resolution diffusion imaging data should validate their results by considering, at least, voxels with one, two, and three fibers, as well as voxels with dominant fibers and different diffusion anisotropies.

Converging Medial Frontal Resting State and Diffusion-Based Abnormalities in Borderline Personality Disorder.

BACKGROUND: The psychological profile of patients with borderline personality disorder (BPD), with impulsivity and emotional dysregulation as core symptoms, has guided the search for abnormalities in specific brain areas such as the hippocampal-amygdala complex and the frontomedial cortex. However, whole-brain imaging studies so far have delivered highly heterogeneous results involving different brain locations. METHODS: Functional resting-state and diffusion magnetic resonance imaging data were acquired in patients with BPD and in an equal number of matched control subjects (n = 60 for resting and n = 43 for diffusion). While mean diffusivity and fractional anisotropy brain images were generated from diffusion data, amplitude of low-frequency fluctuations and global brain connectivity images were used for the first time to evaluate BPD-related brain abnormalities from resting functional acquisitions. RESULTS: Whole-brain analyses using a p = .05 corrected threshold showed a convergence of alterations in BPD patients in genual and perigenual structures, with frontal white matter fractional anisotropy abnormalities partially encircling areas of increased mean diffusivity and global brain connectivity. Additionally, a cluster of enlarged amplitude of low-frequency fluctuations (high resting activity) was found involving part of the left hippocampus and amygdala. In turn, this cluster showed increased resting functional connectivity with the anterior cingulate. CONCLUSIONS: With a multimodal approach and without using a priori selected regions, we prove that structural and functional abnormality in BPD involves both temporolimbic and frontomedial structures as well as their connectivity. These structures have been previously related to behavioral and clinical symptoms in patients with BPD.

Anisotropic kernels for coordinate-based meta-analyses of neuroimaging studies.

Peak-based meta-analyses of neuroimaging studies create, for each study, a brain map of effect size or peak likelihood by convolving a kernel with each reported peak. A kernel is a small matrix applied in order that voxels surrounding the peak have a value similar to, but slightly lower than that of the peak. Current kernels are isotropic, i.e., the value of a voxel close to a peak only depends on the Euclidean distance between the voxel and the peak. However, such perfect spheres of effect size or likelihood around the peak are rather implausible: a voxel that correlates with the peak across individuals is more likely to be part of the cluster of significant activation or difference than voxels uncorrelated with the peak. This paper introduces anisotropic kernels, which assign different values to the different neighboring voxels based on the spatial correlation between them. They are specifically developed for effect-size signed differential mapping (ES-SDM), though might be easily implemented in other meta-analysis packages such as activation likelihood estimation (ALE). The paper also describes the creation of the required correlation templates for gray matter/BOLD response, white matter, cerebrospinal fluid, and fractional anisotropy. Finally, the new method is validated by quantifying the accuracy of the recreation of effect size maps from peak information. This empirical validation showed that the optimal degree of anisotropy and full-width at half-maximum (FWHM) might vary largely depending on the specific data meta-analyzed. However, it also showed that the recreation substantially improved and did not depend on the FWHM when full anisotropy was used. Based on these results, we recommend the use of fully anisotropic kernels in ES-SDM and ALE, unless optimal meta-analysis-specific parameters can be estimated based on the recreation of available statistical maps. The new method and templates are freely available at http://www.sdmproject.com/.

Structural abnormalities in bipolar euthymia: a multicontrast molecular diffusion imaging study.

BACKGROUND: Evidence from decades of magnetic resonance imaging (MRI) research in bipolar disorder has been summarized in meta-analyses of various MRI modalities. Notably, although structural MRI studies suggest gray matter reductions are restricted to specific cortical regions, functional MRI has also shown involvement of subcortical dysfunction. Such disparity in results is open to discussion and requires further exploration with additional MRI modalities. METHODS: We applied whole-brain high angular resolution molecular diffusion imaging to compare different properties of the water diffusion process in brain tissues, using different contrasts. Specifically, we looked at fractional anisotropy, mean diffusivity, probability of return to the origin, and generalized fractional anisotropy in a sample of 40 euthymic patients with bipolar disorder and 40 well-matched healthy control subjects. RESULTS: Convergent abnormalities were detected by contrasts in various tissue types. Apart from alterations in white matter (in corpus callosum, cingulum bundle, corona radiata, and superior fronto-occipital fasciculus) and cortical gray matter (in medial frontal cortex, left insula, Heschl's gyrus, and cerebellum), three of the contrasts (fractional anisotropy, mean diffusivity, and generalized fractional anisotropy) revealed abnormalities in subcortical structures, including the hippocampus, the thalamus and the caudate nucleus. CONCLUSIONS: Our findings point to a wider pattern of axonal pathology in bipolar disorder than previously thought. Although findings related to cortical gray matter are consistent with structural meta-analyses, subcortical abnormalities suggest a cytoarchitectonic basis for previously reported subcortical dysfunction. Diffusion results could be interpreted in terms of loss of tissue volume and/or altered membrane permeability, agreeing with both hypotheses of mitochondrial malfunction and neuroinflammation.

Validity of modulation and optimal settings for advanced voxel-based morphometry.

Voxel-based morphometry (VBM) is a widely-used structural neuroimaging technique for comparing meso- and macroscopic regional brain volumes between patients and controls in vivo, but some of its steps, particularly the modulation, lack an experimental validation. The aims of this study were two-fold: a) to assess the effects of modulation to detect mesoscopic (i.e. between microscopic and macroscopic) abnormalities on published, classic VBM; and b) to suggest a set of potentially optimal settings for detecting mesoscopic abnormalities with new, advanced, high-resolution diffeomorphic VBM normalization algorithms. Sensitivity and false positive rate after modulating or not in classic VBM using different software packages and spatial statistics, and after setting a range of different parameters in advanced VBM (ANTS-SyN), were calculated in 10 VBM comparisons of 32 altered vs. 32 unaltered gray matter images from different healthy controls. Simulated brain abnormalities comprised mesoscopic volume differences mainly due to cortical thinning. In classic VBM, modulation was associated with a substantial decrease of the sensitivity to detect mesoscopic abnormalities (p<0.001). Optimal settings for advanced VBM included the omission of modulation, the use of large smoothing kernels, and the application of voxel-based or threshold-free cluster enhancement (TFCE) spatial statistics. The modulation-related decrease in sensitivity was due to an increase in variance, and it was more severe in higher-resolution normalization algorithms. Findings from this study suggest the use of unmodulated VBM to detect mesoscopic abnormalities such as cortical thinning.