Participants
The study sample included 20 HCs (individuals without any history of psychiatric disorders according to the DSM-5 criteria and no lifetime history of suicide attempt) and 42 patients with MDD [21 SAs with history of suicide attempt within the last 8 days, and 21 affective controls (ACs) without lifetime history of suicide attempt].
HCs were recruited through advertisement. Patients (SAs and ACs) were recruited among the inpatients of the Department of Emergency Psychiatry and Acute Care, Montpellier Academic Hospital, France. All participants were right-handed, Caucasian adults aged between 18 and 70 years. All patients (ACs and SAs) had a current depressive episode according to the DSM-5 criteria. Exclusion criteria were: current neurologic or inflammatory disease; CRP >10mg/L; current antibiotic or anti-inflammatory treatment; history of head trauma with loss of consciousness; current substance misuse, excluding tobacco; lifetime history of manic or hypomanic episode, schizophrenia, or schizoaffective disorder according to the DSM-5; pregnancy and breastfeeding; contraindication to MRI.
The study was conducted in accordance with the CONSORT ethical guidelines. All participants gave their written informed consent. The study was approved by Montpellier University Hospital (CPP Sud Méditerranée IV) ethics committee (N°ID-RCB : 2016-01267-44 (CPP Sud Méditerranée IV) and registered as NCT03052855 in clinicaltrial.gov. Participants in the study received a financial compensation of 50 euros.
Clinical assessment
Participants underwent a standardized interview led by a psychiatrist to determine the depression score and to assess psychopathology according to the DSM-5 criteria using the Mini-International Neuropsychiatric Interview, version 7.0.0. Depressive symptomatology was assessed using the clinician-rated Inventory of Depressive Symptomatology (IDSC-30) 29. Current psychological pain was assessed using a 10-point Likert scale 30. Suicidal ideation severity and intensity in the week before inclusion were assessed with the Columbia–Suicide Severity Rating Scale (C-SSRS) 31, and suicidal intent with the Suicidal Intent Scale (SIS) 32. The number of past suicide attempts and the lethality of the last attempt were recorded. All psychotropic treatments (name, daily dose) in the last 24 hours before MRI were recorded. The medication load was calculated for each patient according to 33.
Laboratory analyses
Blood samples were collected after overnight fasting, between 7:00 and 9:00 AM, the day after admission, in 5 ml EDTA-coated tubes (BD Vacutainer, Franklin Lakes, NJ). High sensitivity CRP in serum was measured with the immunoturbidimetric method and a Cobas8000 modular analyzer (Roche Diagnostic, Meylan, France). Reagents and calibrators were used according to the manufacturer’s guidelines with analytic measuring ranges set at 0.3–350 mg/L.
Plasma was separated by centrifugation at 3000×g at room temperature for 20 minutes, and stored at -80°C for batch analysis by multiplex or ELISA assay of IL-1b, IL-6, IL-2, and tumor necrosis factor alpha (TNF-a). Analytes with >30% missing data were excluded from the analysis (IL-1b, IL-2).
Brain MRI
Brain MRI was performed using a 3T MRI scanner with a 32-channel head coil (MAGNETOM Skyra, Siemens Healthcare, Erlangen, Germany) at the I2FH facility of Montpellier university hospital.
A high-resolution 3DT1 magnetization-prepared rapid gradient-echo (MPRAGE) sequence was acquired with the following parameters: repetition time (TR)/echo time (TE) = 2300/2.900ms, time of inversion = 900 ms, flip angle = 9°, voxel size = 1.2 x1 x1 mm3, and 176 sagittal slices. During rs-fMRI gradient echo-echo planar image acquisition (TR/TE = 3000/30 ms, flip angle = 80°, voxel size= 3.4 x 3.4 x 3.4 mm3, 42 transverse slices, 200 volumes measures, acquisition time = 10 min), participants were instructed to keep the eyes closed, lay still, and not to think about anything in particular.
MRI data preprocessing
Imaging data were preprocessed using MATLAB (TheMathworks Inc., MA) and SPM12 (http://www.fil.ion.ucl.ac.uk/spm/software/spm12/). All images were reoriented using the anterior commissure as reference. Standard DARTEL segmentation was performed on T1 images to extract gray matter, white matter, and cerebrospinal fluid posterior probability distribution. The resulting segmentations were used to create a dedicated DARTEL template that was used to normalize all T1 images (Ashburner, 2007).
Resting state images were processed using a standard processing pipeline that included signal heterogeneity correction using ANTS NI4TK, slice timing, motion correction, co-registration on the T1 images, normalization (based on the T1 images), and spatial smoothing (FWHM=6 mm).
Resting-state metrics
ALFF, ReHo, and Voxel-Mirrored Homotopic Connectivity (VMHC) were estimated to evaluate local alterations of signal fluctuations. ALFF, ReHo and VMHC maps were prepared with the REST toolbox (Resting-State fMRI Data Analysis Toolkit V1.8, REST-Group, http://www.restfmri.net). Functional images were detrended and filtered (0.01-0.08 Hz). Smoothing was performed after ReHo score calculation to avoid artificially increasing the similarity between neighboring voxels. Maps were z-scored. The average gray matter ALFF, ReHo and VMHC values were extracted from 116 regions using the AAL parcellation atlas (AAL 116 toolbox, Tzourio-Mazoyer et al. 2002). Based on our previous fMRI studies 34–36 and literature data 3, additional analyses focused on regions implicated in suicidal vulnerability (Figure 1): 1) orbitofrontal cortex (OFC): orbital parts of the inferior frontal gyrus (IFG), middle frontal gyrus and medial frontal gyrus (corresponding to the Brodmann areas (BA) 11/47); 2) ventrolateral prefrontal cortex (VLPFC): pars opercularis of the IFG (BA44/45); 3) dorsal prefrontal cortex (DPFC): middle frontal gyrus (BA 8/9/46); 4) anterior cingulate cortex (ACC) (BA 24/32); and 5) insula.
Statistical analysis
Statistical analyses were performed with R version 4.1.2 (R Foundation for Statistical Computing, Vienna, Austria). A general linear model was used to model each regional measure as an outcome, while adjusting for age, sex, and education level. To identify the main group effects, a type II analysis of variance was performed with the fitted model. A false discovery rate (FDR) procedure was used for each set of rs-fMRI metrics to correct for multiple comparisons and post hoc analyses were performed for regions with a FDR p <.05. Exploratory correlation analyses were performed to examine possible associations between rs-fMRI metrics with a main group effect and clinical variables (depression, psychological pain intensity and suicidal ideation intensity) in all patients and also suicidal intent and lethality of last suicide attempt in SAs. Exploratory correlation analyses were also performed to identify possible associations between rs-fMRI metrics and peripheral inflammation markers in all patients. All correlation analyses were followed by multiple comparison correction (FDR).