Using brain structural neuroimaging measures to predict psychosis onset for individuals at clinical high-risk.

Machine learning approaches using structural magnetic resonance imaging (sMRI) can be informative for disease classification, although their ability to predict psychosis is largely unknown. We created a model with individuals at CHR who developed psychosis later (CHR-PS+) from healthy controls (HCs) that can differentiate each other. We also evaluated whether we could distinguish CHR-PS+ individuals from those who did not develop psychosis later (CHR-PS-) and those with uncertain follow-up status (CHR-UNK). T1-weighted structural brain MRI scans from 1165 individuals at CHR (CHR-PS+, n = 144; CHR-PS-, n = 793; and CHR-UNK, n = 228), and 1029 HCs, were obtained from 21 sites. We used ComBat to harmonize measures of subcortical volume, cortical thickness and surface area data and corrected for non-linear effects of age and sex using a general additive model. CHR-PS+ (n = 120) and HC (n = 799) data from 20 sites served as a training dataset, which we used to build a classifier. The remaining samples were used external validation datasets to evaluate classifier performance (test, independent confirmatory, and independent group [CHR-PS- and CHR-UNK] datasets). The accuracy of the classifier on the training and independent confirmatory datasets was 85% and 73% respectively. Regional cortical surface area measures-including those from the right superior frontal, right superior temporal, and bilateral insular cortices strongly contributed to classifying CHR-PS+ from HC. CHR-PS- and CHR-UNK individuals were more likely to be classified as HC compared to CHR-PS+ (classification rate to HC: CHR-PS+, 30%; CHR-PS-, 73%; CHR-UNK, 80%). We used multisite sMRI to train a classifier to predict psychosis onset in CHR individuals, and it showed promise predicting CHR-PS+ in an independent sample. The results suggest that when considering adolescent brain development, baseline MRI scans for CHR individuals may be helpful to identify their prognosis. Future prospective studies are required about whether the classifier could be actually helpful in the clinical settings.

Dissecting depression symptoms: Multi-omics clustering uncovers immune-related subgroups and cell-type specific dysregulation.

In a subset of patients with mental disorders, such as depression, low-grade inflammation and altered immune marker concentrations are observed. However, these immune alterations are often assessed by only one data type and small marker panels. Here, we used a transdiagnostic approach and combined data from two cohorts to define subgroups of depression symptoms across the diagnostic spectrum through a large-scale multi-omics clustering approach in 237 individuals. The method incorporated age, body mass index (BMI), 43 plasma immune markers and RNA-seq data from peripheral mononuclear blood cells (PBMCs). Our initial clustering revealed four clusters, including two immune-related depression symptom clusters characterized by elevated BMI, higher depression severity and elevated levels of immune markers such as interleukin-1 receptor antagonist (IL-1RA), C-reactive protein (CRP) and C-C motif chemokine 2 (CCL2 or MCP-1). In contrast, the RNA-seq data mostly differentiated a cluster with low depression severity, enriched in brain related gene sets. This cluster was also distinguished by electrocardiography data, while structural imaging data revealed differences in ventricle volumes across the clusters. Incorporating predicted cell type proportions into the clustering resulted in three clusters, with one showing elevated immune marker concentrations. The cell type proportion and genes related to cell types were most pronounced in an intermediate depression symptoms cluster, suggesting that RNA-seq and immune markers measure different aspects of immune dysregulation. Lastly, we found a dysregulation of the SERPINF1/VEGF-A pathway that was specific to dendritic cells by integrating immune marker and RNA-seq data. This shows the advantages of combining different data modalities and highlights possible markers for further stratification research of depression symptoms.

Functional connectivity signatures of major depressive disorder: machine learning analysis of two multicenter neuroimaging studies.

The promise of machine learning has fueled the hope for developing diagnostic tools for psychiatry. Initial studies showed high accuracy for the identification of major depressive disorder (MDD) with resting-state connectivity, but progress has been hampered by the absence of large datasets. Here we used regular machine learning and advanced deep learning algorithms to differentiate patients with MDD from healthy controls and identify neurophysiological signatures of depression in two of the largest resting-state datasets for MDD. We obtained resting-state functional magnetic resonance imaging data from the REST-meta-MDD (N = 2338) and PsyMRI (N = 1039) consortia. Classification of functional connectivity matrices was done using support vector machines (SVM) and graph convolutional neural networks (GCN), and performance was evaluated using 5-fold cross-validation. Features were visualized using GCN-Explainer, an ablation study and univariate t-testing. The results showed a mean classification accuracy of 61% for MDD versus controls. Mean accuracy for classifying (non-)medicated subgroups was 62%. Sex classification accuracy was substantially better across datasets (73-81%). Visualization of the results showed that classifications were driven by stronger thalamic connections in both datasets, while nearly all other connections were weaker with small univariate effect sizes. These results suggest that whole brain resting-state connectivity is a reliable though poor biomarker for MDD, presumably due to disease heterogeneity as further supported by the higher accuracy for sex classification using the same methods. Deep learning revealed thalamic hyperconnectivity as a prominent neurophysiological signature of depression in both multicenter studies, which may guide the development of biomarkers in future studies.

Early spinal cord pseudoatrophy in interferon-beta-treated multiple sclerosis.

BACKGROUND AND PURPOSE: Brain pseudoatrophy has been shown to play a pivotal role in the interpretation of brain atrophy measures during the first year of disease-modifying therapy in multiple sclerosis. Whether pseudoatrophy also affects the spinal cord remains unclear. The aim of this study was to analyze the extent of pseudoatrophy in the upper spinal cord during the first 2 years after therapy initiation and compare this to the brain. METHODS: A total of 129 patients from a prospective longitudinal multicentric national cohort study for whom magnetic resonance imaging scans at baseline, 12 months, and 24 months were available were selected for brain and spinal cord volume quantification. Annual percentage brain volume and cord area change were calculated using SIENA (Structural Image Evaluation of Normalized Atrophy) and NeuroQLab, respectively. Linear mixed model analyses were performed to compare patients on interferon-beta therapy (n = 84) and untreated patients (n = 45). RESULTS: Patients treated with interferon-beta demonstrated accelerated annual percentage brain volume and cervical cord area change in the first year after treatment initiation, whereas atrophy rates stabilized to a similar and not significantly different level compared to untreated patients during the second year. CONCLUSIONS: These results suggest that pseudoatrophy occurs not only in the brain, but also in the spinal cord during the first year of interferon-beta treatment.

Pupillometry tracks cognitive load and salience network activity in a working memory functional magnetic resonance imaging task.

The diameter of the human pupil tracks working memory processing and is associated with activity in the frontoparietal network. At the same time, recent neuroimaging research has linked human pupil fluctuations to activity in the salience network. In this combined functional magnetic resonance imaging (fMRI)/pupillometry study, we recorded the pupil size of healthy human participants while they performed a blockwise organized working memory task (N-back) inside an MRI scanner in order to monitor the pupil fluctuations associated neural activity during working memory processing. We first confirmed that mean pupil size closely followed working memory load. Combining this with fMRI data, we focused on blood oxygen level dependent (BOLD) correlates of mean pupil size modeled onto the task blocks as a parametric modulation. Interrogating this modulated task regressor, we were able to retrieve the frontoparietal network. Next, to fully exploit the within-block dynamics, we divided the blocks into 1 s time bins and filled these with corresponding pupil change values (first-order derivative of pupil size). We found that pupil change within N-back blocks was positively correlated with BOLD amplitudes in the areas of the salience network (namely bilateral insula, and anterior cingulate cortex). Taken together, fMRI with simultaneous measurement of pupil parameters constitutes a valuable tool to dissect working memory subprocesses related to both working memory load and salience of the presented stimuli.

FreeSurfer-based segmentation of hippocampal subfields: A review of methods and applications, with a novel quality control procedure for ENIGMA studies and other collaborative efforts.

Structural hippocampal abnormalities are common in many neurological and psychiatric disorders, and variation in hippocampal measures is related to cognitive performance and other complex phenotypes such as stress sensitivity. Hippocampal subregions are increasingly studied, as automated algorithms have become available for mapping and volume quantification. In the context of the Enhancing Neuro Imaging Genetics through Meta Analysis Consortium, several Disease Working Groups are using the FreeSurfer software to analyze hippocampal subregion (subfield) volumes in patients with neurological and psychiatric conditions along with data from matched controls. In this overview, we explain the algorithm's principles, summarize measurement reliability studies, and demonstrate two additional aspects (subfield autocorrelation and volume/reliability correlation) with illustrative data. We then explain the rationale for a standardized hippocampal subfield segmentation quality control (QC) procedure for improved pipeline harmonization. To guide researchers to make optimal use of the algorithm, we discuss how global size and age effects can be modeled, how QC steps can be incorporated and how subfields may be aggregated into composite volumes. This discussion is based on a synopsis of 162 published neuroimaging studies (01/2013-12/2019) that applied the FreeSurfer hippocampal subfield segmentation in a broad range of domains including cognition and healthy aging, brain development and neurodegeneration, affective disorders, psychosis, stress regulation, neurotoxicity, epilepsy, inflammatory disease, childhood adversity and posttraumatic stress disorder, and candidate and whole genome (epi-)genetics. Finally, we highlight points where FreeSurfer-based hippocampal subfield studies may be optimized.

Subcortical shape alterations in major depressive disorder: Findings from the ENIGMA major depressive disorder working group.

Alterations in regional subcortical brain volumes have been investigated as part of the efforts of an international consortium, ENIGMA, to identify reliable neural correlates of major depressive disorder (MDD). Given that subcortical structures are comprised of distinct subfields, we sought to build significantly from prior work by precisely mapping localized MDD-related differences in subcortical regions using shape analysis. In this meta-analysis of subcortical shape from the ENIGMA-MDD working group, we compared 1,781 patients with MDD and 2,953 healthy controls (CTL) on individual measures of shape metrics (thickness and surface area) on the surface of seven bilateral subcortical structures: nucleus accumbens, amygdala, caudate, hippocampus, pallidum, putamen, and thalamus. Harmonized data processing and statistical analyses were conducted locally at each site, and findings were aggregated by meta-analysis. Relative to CTL, patients with adolescent-onset MDD (≤ 21 years) had lower thickness and surface area of the subiculum, cornu ammonis (CA) 1 of the hippocampus and basolateral amygdala (Cohen's d = -0.164 to -0.180). Relative to first-episode MDD, recurrent MDD patients had lower thickness and surface area in the CA1 of the hippocampus and the basolateral amygdala (Cohen's d = -0.173 to -0.184). Our results suggest that previously reported MDD-associated volumetric differences may be localized to specific subfields of these structures that have been shown to be sensitive to the effects of stress, with important implications for mapping treatments to patients based on specific neural targets and key clinical features.

Brain structural abnormalities in obesity: relation to age, genetic risk, and common psychiatric disorders : Evidence through univariate and multivariate mega-analysis including 6420 participants from the ENIGMA MDD working group.

Emerging evidence suggests that obesity impacts brain physiology at multiple levels. Here we aimed to clarify the relationship between obesity and brain structure using structural MRI (n = 6420) and genetic data (n = 3907) from the ENIGMA Major Depressive Disorder (MDD) working group. Obesity (BMI > 30) was significantly associated with cortical and subcortical abnormalities in both mass-univariate and multivariate pattern recognition analyses independent of MDD diagnosis. The most pronounced effects were found for associations between obesity and lower temporo-frontal cortical thickness (maximum Cohen´s d (left fusiform gyrus) = -0.33). The observed regional distribution and effect size of cortical thickness reductions in obesity revealed considerable similarities with corresponding patterns of lower cortical thickness in previously published studies of neuropsychiatric disorders. A higher polygenic risk score for obesity significantly correlated with lower occipital surface area. In addition, a significant age-by-obesity interaction on cortical thickness emerged driven by lower thickness in older participants. Our findings suggest a neurobiological interaction between obesity and brain structure under physiological and pathological brain conditions.

White matter disturbances in major depressive disorder: a coordinated analysis across 20 international cohorts in the ENIGMA MDD working group.

Alterations in white matter (WM) microstructure have been implicated in the pathophysiology of major depressive disorder (MDD). However, previous findings have been inconsistent, partially due to low statistical power and the heterogeneity of depression. In the largest multi-site study to date, we examined WM anisotropy and diffusivity in 1305 MDD patients and 1602 healthy controls (age range 12-88 years) from 20 samples worldwide, which included both adults and adolescents, within the MDD Working Group of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium. Processing of diffusion tensor imaging (DTI) data and statistical analyses were harmonized across sites and effects were meta-analyzed across studies. We observed subtle, but widespread, lower fractional anisotropy (FA) in adult MDD patients compared with controls in 16 out of 25 WM tracts of interest (Cohen's d between 0.12 and 0.26). The largest differences were observed in the corpus callosum and corona radiata. Widespread higher radial diffusivity (RD) was also observed (all Cohen's d between 0.12 and 0.18). Findings appeared to be driven by patients with recurrent MDD and an adult age of onset of depression. White matter microstructural differences in a smaller sample of adolescent MDD patients and controls did not survive correction for multiple testing. In this coordinated and harmonized multisite DTI study, we showed subtle, but widespread differences in WM microstructure in adult MDD, which may suggest structural disconnectivity in MDD.

Common genetic variants influence human subcortical brain structures.

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.

The brain's hemodynamic response function rapidly changes under acute psychosocial stress in association with genetic and endocrine stress response markers.

Ample evidence links dysregulation of the stress response to the risk for psychiatric disorders. However, we lack an integrated understanding of mechanisms that are adaptive during the acute stress response but potentially pathogenic when dysregulated. One mechanistic link emerging from rodent studies is the interaction between stress effectors and neurovascular coupling, a process that adjusts cerebral blood flow according to local metabolic demands. Here, using task-related fMRI, we show that acute psychosocial stress rapidly impacts the peak latency of the hemodynamic response function (HRF-PL) in temporal, insular, and prefrontal regions in two independent cohorts of healthy humans. These latency effects occurred in the absence of amplitude effects and were moderated by regulatory genetic variants of KCNJ2, a known mediator of the effect of stress on vascular responsivity. Further, hippocampal HRF-PL correlated with both cortisol response and genetic variants that influence the transcriptional response to stress hormones and are associated with risk for major depression. We conclude that acute stress modulates hemodynamic response properties as part of the physiological stress response and suggest that HRF indices could serve as endophenotype of stress-related disorders.

Psychosocial stress reactivity habituates following acute physiological stress.

Acute and chronic stress are important factors in the development of mental disorders. Reliable measurement of stress reactivity is therefore pivotal. Critically, experimental induction of stress often involves multiple "hits" and it is an open question whether individual differences in responses to an earlier stressor lead to habituation, sensitization, or simple additive effects on following events. Here, we investigated the effect of the individual cortisol response to intravenous catheter placement (IVP) on subsequent neural, psychological, endocrine, and autonomous stress reactivity. We used an established psychosocial stress paradigm to measure the acute stress response (Stress) and recovery (PostStress) in 65 participants. Higher IVP-induced cortisol responses were associated with lower pulse rate increases during stress recovery (b = -4.8 bpm, p = .0008) and lower increases in negative affect after the task (b = -4.2, p = .040). While the cortisol response to IVP was not associated with subsequent specific stress-induced neural activation patterns, the similarity of brain responses Pre- and PostStress was higher IVP-cortisol responders (t[64] = 2.35, p = .022) indicating faster recovery. In conclusion, preparatory stress induced by IVP reduced reactivity in a subsequent stress task by modulating the latency of stress recovery. Thus, an individually stronger preceding release of cortisol may attenuate a second physiological response and perceived stress suggesting that relative changes, not absolute levels are crucial for stress attribution. Our study highlights that considering the entire trajectory of stress induction during an experiment is important to develop reliable individual biomarkers.

Looping Star fMRI in Cognitive Tasks and Resting State.

BACKGROUND: Conventional T2 *-weighted functional magnetic resonance imaging (fMRI) is performed with echo-planar imaging (EPI) sequences that create substantial acoustic noise. The loud acoustic noise not only affects the activation of the auditory cortex, but may also interfere with resting state and task fMRI experiments. PURPOSE: To demonstrate the feasibility of a novel, quiet, T2 *, whole-brain blood oxygenation level-dependent (BOLD)-fMRI method, termed Looping Star, compared to conventional multislice gradient-echo EPI. STUDY TYPE: Prospective. PHANTOM/SUBJECTS: Glover stability QA phantom; 10 healthy volunteers. FIELD STRENGTH/SEQUENCE: 3.0T: gradient echo (GE)-EPI and T2 * Looping Star fMRI. ASSESSMENT: Looping Star fMRI was presented and compared to GE-EPI with a working memory (WM) task and resting state (RS) experiments. Temporal stability and acoustic measurements were obtained for both methods. Functional maps and activation accuracy were compared to evaluate the performance of the novel sequence. STATISTICAL TESTS: Mean and standard deviation values were analyzed for temporal stability and acoustic noise tests. Activation maps were assessed with one-sample t-tests and contrast estimates (CE). Paired t-tests and receiver operator characteristic (ROC) were used to compare fMRI sensitivity and performance. RESULTS: Looping Star presented a 98% reduction in sound pressure compared with GE-EPI, with stable temporal stability (0.09% percent fluctuation), but reduced temporal signal-to-noise ratio (tSNR) (mean difference = 15.9%). The novel method yielded consistent activations for RS and WM (83.4% and 69.5% relative BOLD sensitivity), which increased with task difficulty (mean CE 2-back = 0.56 vs. 0-back = 0.08, P < 0.05). A few differences in spatial activations were found between sequences, leading to a 4-8% lower activation accuracy with Looping Star. DATA CONCLUSION: Looping Star provides a suitable approach for whole-brain coverage with sufficient spatiotemporal resolution and BOLD sensitivity, with only 0.5 dB above ambient noise. From the comparison with GE-EPI, further developments of Looping Star fMRI should target increased sensitivity and spatial specificity for both RS and task experiments. LEVEL OF EVIDENCE: 2. TECHNICAL EFFICACY STAGE: 1 J. Magn. Reson. Imaging 2020;52:739-751.

Schema therapy versus cognitive behavioral therapy versus individual supportive therapy for depression in an inpatient and day clinic setting: study protocol of the OPTIMA-RCT.

BACKGROUND: Major depressive disorder represents (MDD) a major cause of disability and disease burden. Beside antidepressant medication, psychotherapy is a key approach of treatment. Schema therapy has been shown to be effective in the treatment of psychiatric disorders, especially personality disorders, in a variety of settings and patient groups. Nevertheless, there is no evidence on its effectiveness for MDD in an inpatient nor day clinic setting and little is known about the factors that drive treatment response in such a target group. METHODS: In the current protocol, we outline OPTIMA (OPtimized Treatment Identification at the MAx Planck Institute): a single-center randomized controlled trial of schema therapy as a treatment approach for MDD in an inpatient and day clinic setting. Over the course of 7 weeks, we compare schema therapy with cognitive behavioral therapy and individual supportive therapy, conducted in individual and group sessions and with no restrictions regarding concurrent antidepressant medication, thus approximating real-life treatment conditions. N = 300 depressed patients are included. All study therapists undergo a specific training and supervision and therapy adherence is assessed. Primary outcome is depressive symptom severity as self-assessment (Beck Depression Inventory-II) and secondary outcomes are clinical ratings of MDD (Montgomery-Asberg Depression Rating Scale), recovery rates after 7 weeks according to the Munich-Composite International Diagnostic Interview, general psychopathology (Brief Symptom Inventory), global functioning (World Health Organization Disability Assessment Schedule), and clinical parameters such as dropout rates. Further parameters on a behavioral, cognitive, psychophysiological, and biological level are measured before, during and after treatment and in 2 follow-up assessments after 6 and 24 months after end of treatment. DISCUSSION: To our knowledge, the OPTIMA-Trial is the first to investigate the effectiveness of schema therapy as a treatment approach of MDD, to investigate mechanisms of change, and explore predictors of treatment response in an inpatient and day clinic setting by using such a wide range of parameters. Insights from OPTIMA will allow more integrative approaches of psychotherapy of MDD. Especially, the identification of intervention-specific markers of treatment response can improve evidence-based clinical decision for individualizing treatment. TRIAL REGISTRATION: Identifier on clinicaltrials.gov : NCT03287362 ; September, 12, 2017.

The biological classification of mental disorders (BeCOME) study: a protocol for an observational deep-phenotyping study for the identification of biological subtypes.

BACKGROUND: A major research finding in the field of Biological Psychiatry is that symptom-based categories of mental disorders map poorly onto dysfunctions in brain circuits or neurobiological pathways. Many of the identified (neuro) biological dysfunctions are "transdiagnostic", meaning that they do not reflect diagnostic boundaries but are shared by different ICD/DSM diagnoses. The compromised biological validity of the current classification system for mental disorders impedes rather than supports the development of treatments that not only target symptoms but also the underlying pathophysiological mechanisms. The Biological Classification of Mental Disorders (BeCOME) study aims to identify biology-based classes of mental disorders that improve the translation of novel biomedical findings into tailored clinical applications. METHODS: BeCOME intends to include at least 1000 individuals with a broad spectrum of affective, anxiety and stress-related mental disorders as well as 500 individuals unaffected by mental disorders. After a screening visit, all participants undergo in-depth phenotyping procedures and omics assessments on two consecutive days. Several validated paradigms (e.g., fear conditioning, reward anticipation, imaging stress test, social reward learning task) are applied to stimulate a response in a basic system of human functioning (e.g., acute threat response, reward processing, stress response or social reward learning) that plays a key role in the development of affective, anxiety and stress-related mental disorders. The response to this stimulation is then read out across multiple levels. Assessments comprise genetic, molecular, cellular, physiological, neuroimaging, neurocognitive, psychophysiological and psychometric measurements. The multilevel information collected in BeCOME will be used to identify data-driven biologically-informed categories of mental disorders using cluster analytical techniques. DISCUSSION: The novelty of BeCOME lies in the dynamic in-depth phenotyping and omics characterization of individuals with mental disorders from the depression and anxiety spectrum of varying severity. We believe that such biology-based subclasses of mental disorders will serve as better treatment targets than purely symptom-based disease entities, and help in tailoring the right treatment to the individual patient suffering from a mental disorder. BeCOME has the potential to contribute to a novel taxonomy of mental disorders that integrates the underlying pathomechanisms into diagnoses. TRIAL REGISTRATION: Retrospectively registered on June 12, 2019 on ClinicalTrials.gov (TRN: NCT03984084).

Association of smoking but not HLA-DRB1*15:01, APOE or body mass index with brain atrophy in early multiple sclerosis.

BACKGROUND: The course of multiple sclerosis (MS) shows substantial inter-individual variability. The underlying determinants of disease severity likely involve genetic and environmental factors. OBJECTIVE: The aim of this study was to assess the impact of APOE and HLA polymorphisms as well as smoking and body mass index (BMI) in the very early MS course. METHODS: Untreated patients ( n = 263) with a recent diagnosis of relapsing-remitting (RR) MS or clinically isolated syndrome underwent standardized magnetic resonance imaging (MRI). Genotyping was performed for single-nucleotide polymorphisms (SNPs) rs3135388 tagging the HLA-DRB1*15:01 haplotype and rs7412 (Ɛ2) and rs429358 (Ɛ4) in APOE. Linear regression analyses were applied based on the three SNPs, smoking and BMI as exposures and MRI surrogate markers for disease severity as outcomes. RESULTS: Current smoking was associated with reduced gray matter fraction, lower brain parenchymal fraction and increased cerebrospinal fluid fraction in comparison to non-smoking, whereas no effect was observed on white matter fraction. BMI and the SNPs in HLA and APOE were not associated with structural MRI parameters. CONCLUSIONS: Smoking may have an unfavorable effect on the gray matter fraction as a potential measure of MS severity already in early MS. These findings may impact patients' counseling upon initial diagnosis of MS.

Disentangling reward anticipation with simultaneous pupillometry / fMRI.

The reward system may provide an interesting intermediate phenotype for anhedonia in affective disorders. Reward anticipation is characterized by an increase in arousal, and previous studies have linked the anterior cingulate cortex (ACC) to arousal responses such as dilation of the pupil. Here, we examined pupil dynamics during a reward anticipation task in forty-six healthy human subjects and evaluated its neural correlates using functional magnetic resonance imaging (fMRI). Pupil size showed a strong increase during monetary reward anticipation, a moderate increase during verbal reward anticipation and a decrease during control trials. For fMRI analyses, average pupil size and pupil change were computed in 1-s time bins during the anticipation phase. Activity in the ventral striatum was inversely related to the pupil size time course, indicating an early onset of activation and a role in reward prediction processing. Pupil dilations were linked to increased activity in the salience network (dorsal ACC and bilateral insula), which likely triggers an increase in arousal to enhance task performance. Finally, increased pupil size preceding the required motor response was associated with activity in the ventral attention network. In sum, pupillometry provides an effective tool for disentangling different phases of reward anticipation, with relevance for affective symptomatology.

A diffusion model-free framework with echo time dependence for free-water elimination and brain tissue microstructure characterization.

PURPOSE: The compartmental nature of brain tissue microstructure is typically studied by diffusion MRI, MR relaxometry or their correlation. Diffusion MRI relies on signal representations or biophysical models, while MR relaxometry and correlation studies are based on regularized inverse Laplace transforms (ILTs). Here we introduce a general framework for characterizing microstructure that does not depend on diffusion modeling and replaces ill-posed ILTs with blind source separation (BSS). This framework yields proton density, relaxation times, volume fractions, and signal disentanglement, allowing for separation of the free-water component. THEORY AND METHODS: Diffusion experiments repeated for several different echo times, contain entangled diffusion and relaxation compartmental information. These can be disentangled by BSS using a physically constrained nonnegative matrix factorization. RESULTS: Computer simulations, phantom studies, together with repeatability and reproducibility experiments demonstrated that BSS is capable of estimating proton density, compartmental volume fractions and transversal relaxations. In vivo results proved its potential to correct for free-water contamination and to estimate tissue parameters. CONCLUSION: Formulation of the diffusion-relaxation dependence as a BSS problem introduces a new framework for studying microstructure compartmentalization, and a novel tool for free-water elimination.

Multi-echo EPI of human fear conditioning reveals improved BOLD detection in ventromedial prefrontal cortex.

Standard T2* weighted functional magnetic resonance imaging (fMRI) performed with echo-planar imaging (EPI) suffers from signal loss in the ventromedial prefrontal cortex (vmPFC) due to macroscopic field inhomogeneity. However, this region is of special interest to affective neuroscience and psychiatry. The Multi-echo EPI (MEPI) approach has several advantages over EPI but its performance against EPI in the vmPFC has not yet been examined in a study with sufficient statistical power using a task specifically eliciting activity in this region. We used a fear conditioning task with MEPI to compare the performance of MEPI and EPI in vmPFC and control regions in 32 healthy young subjects. We analyzed activity associated with short (12ms), standard (29ms) and long (46ms) echo times, and a voxel-wise combination of these three echo times. Behavioral data revealed successful differentiation of the conditioned versus safety stimulus; activity in the vmPFC was shown by the contrast "safety stimulus > conditioned stimulus" as in previous research and proved significantly stronger with the combined MEPI than standard single-echo EPI. Then, we aimed to demonstrate that the additional cluster extent (ventral extension) detected in the vmPFC with MEPI reflects activation in a relevant cluster (i.e., not just non-neuronal noise). To do this, we used resting state data from the same subjects to show that the time-course of this region was both connected to bilateral amygdala and the default mode network. Overall, we demonstrate that MEPI (by means of the weighted sum combination approach) outperforms standard EPI in vmPFC; MEPI performs always at least as good as the best echo time for a given brain region but provides all necessary echo times for an optimal BOLD sensitivity for the whole brain. This is relevant for affective neuroscience and psychiatry given the critical role of the vmPFC in emotion regulation.