Crystal Cooper

Establishing moderators and biosignatures of antidepressant response in clinical care (EMBARC): Rationale and design

UNLABELLED Remission rates for Major Depressive Disorder (MDD) are low and unpredictable for any given antidepressant. No biological or clinical marker has demonstrated sufficient ability to match individuals to efficacious treatment. Biosignatures developed from the systematic exploration of multiple biological markers, which optimize treatment selection for individuals (moderators) and provide early indication of ultimate treatment response (mediators) are needed. The rationale and design of a multi-site, placebo-controlled randomized clinical trial of sertraline examining moderators and mediators of treatment response is described. The target sample is 300 participants with early onset (≤30 years) recurrent MDD. Non-responders to an 8-week trial are switched double blind to either bupropion (for sertraline non-responders) or sertraline (for placebo non-responders) for an additional 8 weeks. Clinical moderators include anxious depression, early trauma, gender, melancholic and atypical depression, anger attacks, Axis II disorder, hypersomnia/fatigue, and chronicity of depression. Biological moderator and mediators include cerebral cortical thickness, task-based fMRI (reward and emotion conflict), resting connectivity, diffusion tensor imaging (DTI), arterial spin labeling (ASL), electroencephalograpy (EEG), cortical evoked potentials, and behavioral/cognitive tasks evaluated at baseline and week 1, except DTI, assessed only at baseline. The study is designed to standardize assessment of biomarkers across multiple sites as well as institute replicable quality control methods, and to use advanced data analytic methods to integrate these markers. A Differential Depression Treatment Response Index (DTRI) will be developed. The data, including biological samples (DNA, RNA, and plasma collected before and during treatment), will become available in a public scientific repository. CLINICAL TRIAL REGISTRATION Establishing Moderators and Biosignatures of Antidepressant Response for Clinical Care for Depression (EMBARC). Identifier: NCT01407094. URL: http://clinicaltrials.gov/show/NCT01407094.

Test-retest reliability of freesurfer measurements within and between sites: Effects of visual approval process

In the last decade, many studies have used automated processes to analyze magnetic resonance imaging (MRI) data such as cortical thickness, which is one indicator of neuronal health. Due to the convenience of image processing software (e.g., FreeSurfer), standard practice is to rely on automated results without performing visual inspection of intermediate processing. In this work, structural MRIs of 40 healthy controls who were scanned twice were used to determine the test–retest reliability of FreeSurfer‐derived cortical measures in four groups of subjects—those 25 that passed visual inspection (approved), those 15 that failed visual inspection (disapproved), a combined group, and a subset of 10 subjects (Travel) whose test and retest scans occurred at different sites. Test–retest correlation (TRC), intraclass correlation coefficient (ICC), and percent difference (PD) were used to measure the reliability in the Destrieux and Desikan–Killiany (DK) atlases. In the approved subjects, reliability of cortical thickness/surface area/volume (DK atlas only) were: TRC (0.82/0.88/0.88), ICC (0.81/0.87/0.88), PD (0.86/1.19/1.39), which represent a significant improvement over these measures when disapproved subjects are included. Travel subjects’ results show that cortical thickness reliability is more sensitive to site differences than the cortical surface area and volume. To determine the effect of visual inspection on sample size required for studies of MRI‐derived cortical thickness, the number of subjects required to show group differences was calculated. Significant differences observed across imaging sites, between visually approved/disapproved subjects, and across regions with different sizes suggest that these measures should be used with caution. Hum Brain Mapp 36:3472–3485, 2015.

Neural Correlates of Three Promising Endophenotypes of Depression: Evidence from the EMBARC Study

Major depressive disorder (MDD) is clinically, and likely pathophysiologically, heterogeneous. A potentially fruitful approach to parsing this heterogeneity is to focus on promising endophenotypes. Guided by the NIMH Research Domain Criteria initiative, we used source localization of scalp-recorded EEG resting data to examine the neural correlates of three emerging endophenotypes of depression: neuroticism, blunted reward learning, and cognitive control deficits. Data were drawn from the ongoing multi-site EMBARC study. We estimated intracranial current density for standard EEG frequency bands in 82 unmedicated adults with MDD, using Low-Resolution Brain Electromagnetic Tomography. Region-of-interest and whole-brain analyses tested associations between resting state EEG current density and endophenotypes of interest. Neuroticism was associated with increased resting gamma (36.5–44 Hz) current density in the ventral (subgenual) anterior cingulate cortex (ACC) and orbitofrontal cortex (OFC). In contrast, reduced cognitive control correlated with decreased gamma activity in the left dorsolateral prefrontal cortex (dlPFC), decreased theta (6.5–8 Hz) and alpha2 (10.5–12 Hz) activity in the dorsal ACC, and increased alpha2 activity in the right dlPFC. Finally, blunted reward learning correlated with lower OFC and left dlPFC gamma activity. Computational modeling of trial-by-trial reinforcement learning further indicated that lower OFC gamma activity was linked to reduced reward sensitivity. Three putative endophenotypes of depression were found to have partially dissociable resting intracranial EEG correlates, reflecting different underlying neural dysfunctions. Overall, these findings highlight the need to parse the heterogeneity of MDD by focusing on promising endophenotypes linked to specific pathophysiological abnormalities.

Test-retest reliability of freesurfer measurements within and between sites

In the last decade, many studies have used automated processes to analyze magnetic resonance imaging (MRI) data such as cortical thickness, which is one indicator of neuronal health. Due to the convenience of image processing software (e.g., FreeSurfer), standard practice is to rely on automated results without performing visual inspection of intermediate processing. In this work, structural MRIs of 40 healthy controls who were scanned twice were used to determine the test–retest reliability of FreeSurfer‐derived cortical measures in four groups of subjects—those 25 that passed visual inspection (approved), those 15 that failed visual inspection (disapproved), a combined group, and a subset of 10 subjects (Travel) whose test and retest scans occurred at different sites. Test–retest correlation (TRC), intraclass correlation coefficient (ICC), and percent difference (PD) were used to measure the reliability in the Destrieux and Desikan–Killiany (DK) atlases. In the approved subjects, reliability of cortical thickness/surface area/volume (DK atlas only) were: TRC (0.82/0.88/0.88), ICC (0.81/0.87/0.88), PD (0.86/1.19/1.39), which represent a significant improvement over these measures when disapproved subjects are included. Travel subjects’ results show that cortical thickness reliability is more sensitive to site differences than the cortical surface area and volume. To determine the effect of visual inspection on sample size required for studies of MRI‐derived cortical thickness, the number of subjects required to show group differences was calculated. Significant differences observed across imaging sites, between visually approved/disapproved subjects, and across regions with different sizes suggest that these measures should be used with caution. Hum Brain Mapp 36:3472–3485, 2015.

Demonstrating test-retest reliability of electrophysiological measures for healthy adults in a multisite study of biomarkers of antidepressant treatment response

Growing evidence suggests that loudness dependency of auditory evoked potentials (LDAEP) and resting EEG alpha and theta may be biological markers for predicting response to antidepressants. In spite of this promise, little is known about the joint reliability of these markers, and thus their clinical applicability. New standardized procedures were developed to improve the compatibility of data acquired with different EEG platforms, and used to examine test-retest reliability for the three electrophysiological measures selected for a multisite project-Establishing Moderators and Biosignatures of Antidepressant Response for Clinical Care (EMBARC). Thirty-nine healthy controls across four clinical research sites were tested in two sessions separated by about 1 week. Resting EEG (eyes-open and eyes-closed conditions) was recorded and LDAEP measured using binaural tones (1000 Hz, 40 ms) at five intensities (60-100 dB SPL). Principal components analysis of current source density waveforms reduced volume conduction and provided reference-free measures of resting EEG alpha and N1 dipole activity to tones from auditory cortex. Low-resolution electromagnetic tomography (LORETA) extracted resting theta current density measures corresponding to rostral anterior cingulate (rACC), which has been implicated in treatment response. There were no significant differences in posterior alpha, N1 dipole, or rACC theta across sessions. Test-retest reliability was .84 for alpha, .87 for N1 dipole, and .70 for theta rACC current density. The demonstration of good-to-excellent reliability for these measures provides a template for future EEG/ERP studies from multiple testing sites, and an important step for evaluating them as biomarkers for predicting treatment response.

Accounting for Dynamic Fluctuations across Time when Examining fMRI Test-Retest Reliability: Analysis of a Reward Paradigm in the EMBARC Study.

Longitudinal investigation of the neural correlates of reward processing in depression may represent an important step in defining effective biomarkers for antidepressant treatment outcome prediction, but the reliability of reward-related activation is not well understood. Thirty-seven healthy control participants were scanned using fMRI while performing a reward-related guessing task on two occasions, approximately one week apart. Two main contrasts were examined: right ventral striatum (VS) activation fMRI BOLD signal related to signed prediction errors (PE) and reward expectancy (RE). We also examined bilateral visual cortex activation coupled to outcome anticipation. Significant VS PE-related activity was observed at the first testing session, but at the second testing session, VS PE-related activation was significantly reduced. Conversely, significant VS RE-related activity was observed at time 2 but not time 1. Increases in VS RE-related activity from time 1 to time 2 were significantly associated with decreases in VS PE-related activity from time 1 to time 2 across participants. Intraclass correlations (ICCs) in VS were very low. By contrast, visual cortex activation had much larger ICCs, particularly in individuals with high quality data. Dynamic changes in brain activation are widely predicted, and failure to account for these changes could lead to inaccurate evaluations of the reliability of functional MRI signals. Conventional measures of reliability cannot distinguish between changes specified by algorithmic models of neural function and noisy signal. Here, we provide evidence for the former possibility: reward-related VS activations follow the pattern predicted by temporal difference models of reward learning but have low ICCs.

Harmonization of cortical thickness measurements across scanners and sites

&NA; With the proliferation of multi‐site neuroimaging studies, there is a greater need for handling non‐biological variance introduced by differences in MRI scanners and acquisition protocols. Such unwanted sources of variation, which we refer to as “scanner effects”, can hinder the detection of imaging features associated with clinical covariates of interest and cause spurious findings. In this paper, we investigate scanner effects in two large multi‐site studies on cortical thickness measurements across a total of 11 scanners. We propose a set of tools for visualizing and identifying scanner effects that are generalizable to other modalities. We then propose to use ComBat, a technique adopted from the genomics literature and recently applied to diffusion tensor imaging data, to combine and harmonize cortical thickness values across scanners. We show that ComBat removes unwanted sources of scan variability while simultaneously increasing the power and reproducibility of subsequent statistical analyses. We also show that ComBat is useful for combining imaging data with the goal of studying life‐span trajectories in the brain. HighlightsCortical thickness (CT) measurements are highly scanner specific.Identifying scanner effects is crucial for inference and biomarker development.We propose to use ComBat to harmonize cortical thickness values across scanners.

A comparison of structural connectivity in anxious depression versus non-anxious depression

BACKGROUND Major depressive disorder (MDD) and anxiety disorders are highly co-morbid. Research has shown conflicting evidence for white matter alteration and amygdala volume reduction in mood and anxiety disorders. To date, no studies have examined differences in structural connectivity between anxious depressed and non-anxious depressed individuals. This study compared fractional anisotropy (FA) and density of selected white matter tracts and amygdala volume between anxious depressed and non-anxious depressed individuals. METHODS 64- direction DTI and T1 scans were collected from 110 unmedicated subjects with MDD, 39 of whom had a co-morbid anxiety disorder diagnosis. Region of interest (ROI) and tractography methods were performed to calculate amygdala volume and FA in the uncinate fasciculus, respectively. Diffusion connectometry was performed to identify whole brain group differences in white matter health. Correlations were computed between biological and clinical measures. RESULTS Tractography and ROI analyses showed no significant differences between bilateral FA values or bilateral amygdala volumes when comparing the anxious depressed and non-anxious depressed groups. The diffusion connectometry analysis showed no significant differences in anisotropy between the groups. Furthermore, there were no significant relationships between MRI-based and clinical measures. CONCLUSION The lack of group differences could indicate that structural connectivity and amygdalae volumes of those with anxious-depression are not significantly altered by a co-morbid anxiety disorder. Improving understanding of anxiety co-morbid with MDD would facilitate development of treatments that more accurately target the underlying networks.

Pretreatment Rostral Anterior Cingulate Cortex Theta Activity in Relation to Symptom Improvement in Depression: A Randomized Clinical Trial

Importance Major depressive disorder (MDD) remains challenging to treat. Although several clinical and demographic variables have been found to predict poor antidepressant response, these markers have not been robustly replicated to warrant implementation in clinical care. Increased pretreatment rostral anterior cingulate cortex (rACC) theta activity has been linked to better antidepressant outcomes. However, no prior study has evaluated whether this marker has incremental predictive validity over clinical and demographic measures. Objective To determine whether increased pretreatment rACC theta activity would predict symptom improvement regardless of randomization arm. Design, Setting, and Participants A multicenter randomized clinical trial enrolled outpatients without psychosis and with chronic or recurrent MDD between July 29, 2011, and December 15, 2015 (Establishing Moderators and Biosignatures of Antidepressant Response for Clinical Care [EMBARC]). Patients were consecutively recruited from 4 university hospitals: 634 patients were screened, 296 were randomized to receive sertraline hydrochloride or placebo, 266 had electroencephalographic (EEG) recordings, and 248 had usable EEG data. Resting EEG data were recorded at baseline and 1 week after trial onset, and rACC theta activity was extracted using source localization. Intent-to-treat analysis was conducted. Data analysis was performed from October 7, 2016, to January 19, 2018. Interventions An 8-week course of sertraline or placebo. Main Outcomes and Measures The 17-item Hamilton Rating Scale for Depression score (assessed at baseline and weeks 1, 2, 3, 4, 6, and 8). Results The 248 participants (160 [64.5%] women, 88 [35.5%] men) with usable EEG data had a mean (SD) age of 36.75 (13.15) years. Higher rACC theta activity at both baseline (b = −1.05; 95% CI, −1.77 to −0.34; P = .004) and week 1 (b = −0.83; 95% CI, −1.60 to −0.06; P < .04) predicted greater depressive symptom improvement, even when controlling for clinical and demographic variables previously linked with treatment outcome. These effects were not moderated by treatment arm. The rACC theta marker, in combination with clinical and demographic variables, accounted for an estimated 39.6% of the variance in symptom change (with 8.5% of the variance uniquely attributable to the rACC theta marker). Conclusions and Relevance Increased pretreatment rACC theta activity represents a nonspecific prognostic marker of treatment outcome. This is the first study to date to demonstrate that rACC theta activity has incremental predictive validity. Trial Registration clinicaltrials.gov Identifier: NCT01407094

Neuroticism and Individual Differences in Neural Function in Unmedicated Major Depression: Findings From the EMBARC Study

BACKGROUND Personality dysfunction represents one of the only predictors of differential response between active treatments for depression to have replicated. In this study, we examine whether depressed patients with higher neuroticism scores, a marker of personality dysfunction, show differences versus depressed patients with lower scores in the functioning of two brain regions associated with treatment response, the anterior cingulate and anterior insula cortices. METHODS Functional magnetic resonance imaging data during an emotional Stroop task were collected from 135 adults diagnosed with major depressive disorder at four academic medical centers participating in the Establishing Moderators and Biosignatures of Antidepressant Response for Clinical Care (EMBARC) study. Secondary analyses were conducted including a sample of 28 healthy individuals. RESULTS In whole-brain analyses, higher neuroticism among depressed adults was associated with increased activity in and connectivity with the right anterior insula cortex to incongruent compared to congruent emotional stimuli (ks>281, ps<0.05 FWE corrected), covarying for concurrent psychiatric distress. We also observed an unanticipated relationship between neuroticism and reduced activity in the precuneus (k=269, p<0.05 FWE corrected). Exploratory analyses including healthy individuals suggested that associations between neuroticism and brain function may be nonlinear over the full range of neuroticism scores. CONCLUSIONS This study provides convergent evidence for the importance of the right anterior insula cortex as a brain-based marker of clinically meaningful individual differences in neuroticism among adults with depression. This is a critical next step in linking personality dysfunction, a replicated clinical predictor of differential antidepressant treatment response, with differences in underlying brain function.