Skip to main content

Novel candidate genes for ECT response prediction—a pilot study analyzing the DNA methylome of depressed patients receiving electroconvulsive therapy



Major depressive disorder (MDD) represents a serious global health concern. The urge for efficient MDD treatment strategies is presently hindered by the incomplete knowledge of its underlying pathomechanism. Despite recent progress (highlighting both genetics and the environment, and thus DNA methylation, to be relevant for its development), 30–50% of MDD patients still fail to reach remission with standard treatment approaches. Electroconvulsive therapy (ECT) is one of the most powerful options for the treatment of pharmacoresistant depression; nevertheless, ECT remission rates barely reach 50% in large-scale naturalistic population-based studies. To optimize MDD treatment strategies and enable personalized medicine in the long- term, prospective indicators of ECT response are thus in great need. Because recent target-driven analyses revealed DNA methylation baseline differences between ECT responder groups, we analyzed the DNA methylome of depressed ECT patients using next-generation sequencing. In this pilot study, we did not only aim to find novel targets for ECT response prediction but also to get a deeper insight into its possible mechanism of action.


Longitudinal DNA methylation analysis of peripheral blood mononuclear cells isolated from a cohort of treatment-resistant MDD patients (n = 12; time points: before and after 1st and last ECT, respectively) using a TruSeq-Methyl Capture EPIC Kit for library preparation, led to the following results: (1) The global DNA methylation differed neither between the four measured time points nor between ECT responders (n = 8) and non-responders (n = 4). (2) Analyzing the DNA methylation variance for every probe (=1476812 single CpG sites) revealed eight novel candidate genes to be implicated in ECT response (protein-coding genes: RNF175, RNF213, TBC1D14, TMC5, WSCD1; genes encoding for putative long non-coding RNA transcripts: AC018685.2, AC098617.1, CLCN3P1). (3) In addition, DNA methylation of two CpG sites (located within AQP10 and TRERF1) was found to change during the treatment course.


We suggest ten novel candidate genes to be implicated in either ECT response or its possible mechanism. Because of the small sample size of our pilot study, our findings must be regarded as preliminary.


The World Health Organization [1] states major depressive disorder (MDD) to be one of the most prevalent mental diseases worldwide. Due to the high number of affected individuals (> 322 million), efficient treatment strategies are required. This need is being challenged by the insufficient knowledge of MDD’s underlying pathophysiology.

Research from recent decades reports nature and nurture to both be relevant for disease development: One’s individual genetic constitution provides a baseline for the vulnerability to certain diseases, but additional environmental factors are often mandatory to provoke their onset [2, 3]. This phenomenon is mediated by epigenetics, i.e., molecular mechanisms (such as DNA methylation (DNAm) and histone modifications) that modulate gene transcription without interfering with the DNA sequence itself [4,5,6,7]. In the case of MDD, animal experiments found the stress reactivity of rodent pups to be associated with their mother’s postnatal grooming behavior. In this context, hippocampal brain cells of neglected animals (in comparison with the pups being intensively cared for) showed a higher DNAm in gene regions encoding the glucocorticoid receptor [8, 9]. As a part of the hypothalamic-pituitary-adrenal (HPA) axis (our central stress response system) disturbances of the latter protein (together with other irregularities) have been suggested to be a cause for the lowered stress resilience found in depressed patients [10,11,12]. The importance of epigenetics for MDD is further underlined by Fuchikami et al., who distinguished depressed subjects from healthy controls simply by analyzing the DNAm of brain-derived neurotrophic factor (BDNF) [13]. BDNF is a neurotrophin shown to be implicated in various neuropsychiatric disorders, including MDD [14,15,16].

Despite this growing body of knowledge, treatment approaches for depression leave much to be desired: The proportion of MDD patients that fail to achieve full remission upon standard medication (30–50%) is still unsatisfyingly high [17]. Electroconvulsive therapy (ECT) has proven superior efficacy and is, therefore, considered to be one of the most powerful options for the treatment of pharmacoresistant depression [18, 19]. However, in naturalistic population-based community-setting studies, ECT remission rates barely reach 50% [20, 21]. To prevent medication failure at baseline and ensure patient-tailored treatment in the long-term, biomarkers predicting ECT response are thus of compelling need. Merely a few clinical characteristics (like age or psychotic symptoms, for instance [21,22,23,24]) serve as a guide for treatment-decision making, but due to MDD’s heterogeneity (compromising various subgroups and a broad spectrum of symptoms), a whole set of biomarkers will be required [25]. In this context, a few biological markers have been recently proposed, as the catechol-O-methyltransferase (COMT) Val158Met (rs4680) [26,27,28] or the dopamine receptor D2 (DRD2) C957T (rs6277) polymorphisms [29, 30]. In the field of epigenetics, our group recently found the DNAm of p11’s promoter (a protein implicated in BDNF production [31]) to reliably predict ECT response in two cohorts of MDD patients [32]. However, none of these targets is likely to reach sufficient sensitivity and specificity to act as an accurate predictor of ECT response alone and other studies on DNAm and ECT are missing [25, 33].

To find further indicators for ECT response prediction, we investigated the methylome of peripheral blood mononuclear cells (PBMCs) isolated from depressed patients undergoing a course of ECT. Most studies conduct their experiments in a target-driven manner and only investigate processes already known to be implicated in MDD. These hypothesis-based analyses provide an essential contribution to the field of biomarker research, though more data-driven approaches are still required so as to not overlook substantial ECT-related information. To address this issue, we used an Illumina EPIC Kit for our study, allowing an analysis of > 3.3 million CpGs located within regions known to be generally implicated in epigenetic mechanisms (such as CpG islands, promoter regions, and open chromatin). We thereby aimed (1) to identify novel targets for the prediction of ECT response and (2) to get a deeper insight into ECT’s general mechanisms.


Patients’ clinical baseline characteristics

Patients’ clinical baseline characteristics are depicted in Table 1. After treatment completion, 10 patients (out of 17) responded to ECT. Four patients had minimally heightened levels of leukocytes (11.2–12.4 × 103/μl), but no signs of infection (i.e., elevated C-reactive protein measures (CRP)). Patients were under medication while receiving ECT, but none were treated with immunomodulatory drugs. During ECT, patients were anesthetized with methohexital (mean = 128.2(± 53.3) mg, minimum = 90 mg, maximum = 250 mg) and remifentanil (89.7(± 49.8) mg, 30 mg, 200 mg) and received succinylcholine for muscle relaxation (114.1(± 45.0) mg, 60 mg, 200 mg). Responders (R) and non-responders (NR) differed only in their body mass index (BMI) (t test, p = 0.015, T = − 2.736, R = 28.4 ± 4.8, NR = 22.4 ± 4.8).

Table 1 Patients’ clinical baseline characteristics (n = 17)

As described in the methods section, patients were excluded from the analysis if the DNAm values of at least one time point were missing (respectively). The clinical baseline characteristics of these patients (n = 12) are reported in Supplementary Table S1. In this subgroup, ECT responders and non-responders differed in their number of total leukocytes (t test, p = 0.048, T = 2.249, R = 8.4 ± 2.2, NR = 5.7 ± 1.2) and their current episode duration (t test, p = 0.026, T = 2.948, R = 24.8 ± 15.2, NR = 60.0 ± 11.3).

ECT and DNA methylation

Analysis of the global DNAm considering ECT response showed no significant effects for time (F(3, 30) = 2.37, p = 0.09), response (F(1, 10) = 0.05, p > 0.1) and the interaction between time and response (F(3, 30) = 0.14, p > 0.1). The analysis of variance for DNAm with respect to response/non-response for every probe (DMP) showed 13 significant probes located in ten different genes (seven protein-coding and three non-protein-coding (pseudo) genes (which encode for putative long non-coding RNA transcripts instead)) that met the previously established criteria of significance. A detailed presentation of the results regarding the 13 significant probes is shown in Table 2 and depicted in Figs. 1, 2, 3 and 4. Information regarding the genetic loci of our significant probes (and the genetic variants possibly affecting our CpGs of interest) was investigated using Ensembl [34], an internal JBrowse [35], GeneCards® [36], and the NHGRI-EBI GWAS Catalog [37].

Table 2 DNA methylation analysis—results
Fig. 1

Comparison of DNA methylation values between ECT responder groups. ECT responders (blue, n = 8) and non-responders (red, n = 4) differed in their baseline DNA methylation at four single CpG sites (indicated by arrows) located within the ring finger protein 213 gene (RNF213). The horizontal rows represent distinct time points from each patient. Their DNA methylation values were clustered based on their similarity. T1: before the 1st ECT, T2: after the 1st ECT, T3: before the last ECT, T4: after the last ECT

Fig. 2

DNA methylation differences between ECT responder groups: protein-coding genes. DNA methylation of ECT responders (n = 8) and non-responders (n = 4) differed at eight CpG sites located within five different protein-coding genes: TBC1D14 (=TBC1 domain family member 14), RNF175 (=ring finger protein 175), TMC5 (=transmembrane channel-like 5), WSCD1 (=WSC domain containing 1), and RNF213 (=ring finger protein 213). Time point 1: before the 1st ECT, 2: after the 1st ECT, 3: before the last ECT, 4: after the last ECT; error bars: ± SD

Fig. 3

DNA methylation differences between ECT responder groups: long non-coding RNA transcripts. DNA methylation of three CpGs located within gene regions encoding for long non-coding RNA transcripts differed between ECT responders (n = 8) and non-responders (n = 4): AC018685.2, AC098617.1, and CLCN3P1 (=chloride channel voltage-sensitive 3 pseudogene 1). Time point 1: before the 1st ECT, 2: after the 1st ECT, 3: before the last ECT, 4: after the last ECT; error bars: ± SD

Fig. 4

DNA methylation changes during the time course of ECT. DNA methylation of two CpG dinucleotides was found to change during the treatment course. These CpGs were located within the following genes: AQP10 (=aquaporin 10) and TRERF1 (=transcriptional regulating factor 1). Time point 1: before the 1st ECT, 2: after the 1st ECT, 3: before the last ECT, 4: after the last ECT; error bars: ± SD

As the duration of the current depressive episode differed between ECT responders (n = 8) and non-responders (n = 4) in the subgroup of patients analyzed (with ECT non-responders suffering from episodes more than twice as long), we conducted an additional correlation analysis by using a Spearman rank-order correlation test. None of our CpG sites reached statistical significance, though a tendency at the gene locus chr16:19488803 (TMC5) was present (rho = − 0.53, S = 437.53, p = 0.076); Figure S1).


DNAm analysis of 1476812 single CpG sites revealed five novel (protein-coding) candidate genes to be implicated in ECT response (RNF175, RNF213, TBC1D14, TMC5, and WSCD1). Further differences between ECT responder groups were found within gene regions encoding for long non-coding RNA transcripts (AC018685.2, AC098617.1, and CLCN3P1). In all cases (except one: AC098617.1), DNAm differed already at baseline and remained stable throughout the time course. Analyzing ECT’s effect irrespective of clinical outcome, DNAm of merely two CpG sites (located within AQP10 and TRERF1) was found to change during the treatment. Intriguingly, all significant CpGs, but one (chr6:42344977, TRERF1), are known to overlap with a single-nucleotide polymorphism (SNP) directly located within these particular dinucleotides, generating or removing these CpGs and thus DNAm as a consequence. Due to the small group size of the current study, the results must be interpreted with caution, particularly due to the differences between ECT responder and non-responder groups. Nevertheless, the identified genes could be important candidates for therapeutic outcome prediction in future studies.

In this regard, the most striking difference in DNAm between ECT responder groups was present at four CpG sites located within the RNF213 gene. Importantly, RNF213’s DNAm has previously been reported to differ between MDD subjects and healthy controls, though in a much lower magnitude than in our cohort and without being comprehensibly corrected for multiple testing [38, 39]. The gene encodes for a homonymous 591-kDa protein (ring finger protein 213) that contains a RING finger domain mediating protein-protein interaction [40, 41]. Together with its postulated AAA+ ATPase and E3 ligase activity, RNF213 is enabled to unfold and link proteins to ubiquitin [41], a small 8.6-kDa protein whose linkage can lead to diverse outcomes depending on the particular amino acid it is bound to. Among ubiquitin’s various roles, its implication in the proteasome protein degradation system is one of the most pronounced [42,43,44]. By these means, RNF213 contributes to the clearance of two proteins involved in vascular remodeling via the Wnt signaling pathway [45]. Its striking role in vascular development is further supported by clinical studies revealing a particular RNF213 mutant (p.R4859K, caused by a SNP of c.14576G>A) to be strongly associated with Moyamoya disease (MMD)—an occlusive cerebrovascular disorder that is marked by progressive stenosis, a concomitant formation of collateral vessels, and transient seizures [46, 47]. Intriguingly, ample evidence links angiogenesis to either MDD or its treatment. In this context, elevated vascular endothelial growth factor (VEGF) mRNA has been found in depressed subjects [48]. Further support for this notion stems from animal experiments, showing hippocampal angiogenesis to be boosted following electroconvulsive stimulation (ECS) [49]. Moreover, clinical neuroimaging studies report a particular SNP (rs699947, 2578C/A; located within the promoter region of VEGF) to be associated with hippocampal volume changes after ECT treatment [50].

As ECT has been demonstrated to have immunomodulatory properties [51,52,53,54] and to (partially) reverse the immunological irregularities found in MDD patients (or at least in a subgroup thereof) [55, 56], the immune system seems to serve as another link between the strong implication of RNF213’s DNAm and the clinical response to ECT. In this context, RNF213 mRNA was found to be predominantly expressed in immunological tissue [46], and its expression to be enhanced upon pro-inflammatory stimulation [57]. In addition, RNF213 has been reported to affect the number of T regulatory cells [58], an immune cell population shown to be reduced in depressed subjects [59]. Finally, another connection between RNF213 and depression is formed by the let-7 family of miRNAs, i.e., short RNA sequences that were found to suppress the common variant of the RNF213 gene. Within this context, particularly let-7c was shown to be either increased or diminished in MMD and MDD patients [60, 61].

According to our analysis, another ring finger protein (RNF175) has been linked to ECT response, though its function is less well characterized. Current studies suggest a SNP located within RNF175 (rs981844) to be associated with the response to statins [62], i.e., a group of pharmaceuticals with suggested antidepressant properties [63]. However, despite this sparsity of literature, one thing is clear: RNF213 and RNF175 share their E3 ubiquitin-ligase activity [64], moving ubiquitin again into the spotlight of ECT responsiveness. Its outstanding role is further supported by several studies suggesting the DNAm of other ring finger proteins (as RNF138, RNF130 [65], and RNF2 [66]) to differ between the postmortem brain samples of depressed subjects and healthy controls, although these results do not reach statistical significance after type 1 error correction. However, besides its involvement in the proteasome degradation system, ubiquitin further leads to the removal of cellular structures by mediating selective autophagy [67].

Similarly, DNAm of the TBC1 domain family member 14 (TBC1D14) gene, a negative regulator of starvation-induced autophagy [68], was also found to differ in relation to clinical outcome, indicating ECT responsiveness to not be dependent on the selective type of autophagy alone. Intriguingly, a profound body of evidence illustrates autophagy to be crucially involved in depression [69]. Alcocer-Gómez et al., for instance, reported the expression of autophagy proteins to be upregulated in MDD patients [70]. Moreover, expression of autophagy proteins (i.e., Beclin-1 and light chain 3-II/I (LC3-II/I)) was shown to be elevated in the rat hippocampus following ECS treatment [71], leading Gassen and Rein to hypothesize that autophagic mechanisms (although already heightened at baseline) might still be insufficient in some disease cases [69]. In the context of ECT, boosted autophagic turnover—ensuring efficient recycling of amino and fatty acids and thus the production of urgently needed proteins [72, 73]—would fit well into the picture as an increase in glucose metabolism has been observed in several brain regions (like the hippocampus, for instance) following ECS treatment [74], indicating enhanced metabolic activity thus a higher demand for nutrients. However, the role of autophagy in MDD is still controversial, but the involvement of RNF213, RNF175, and TBC1D14 in cellular degradation hints at a role for these processes in ECT responsiveness, nevertheless.

The connection of the other genes found in our study (being differentially methylated in ECT responder groups) is less clear, yet no less interesting: An intergenic polymorphism (rs75213074) near the WSC domain containing 1 (WSCD1) gene has been previously associated with migraine [75], a neurological disorder sharing several biological abnormalities with depression [76]. Regarding transmembrane channel-like 5 (TMC5), a SNP (rs4780805) located − 17 kbp upstream of its gene was reported to correlate with sleeping duration [77], forming a link to MDD as sleeping patterns are often disturbed in depressed subjects [78].

With respect to ECT’s general effects (whether directly induced, a coincidence, or merely of secondary origin), DNAm of two CpG sites (located within aquaporin 10 (AQP10) and transcriptional regulating factor 1 (TRERF1)) were altered during the treatment course yet had no effect on clinical outcome. The former gene encodes for a water-permeable channel [79] that has not yet to our knowledge been linked to neuropsychiatric disorders. The same holds true for TRERF1, which regulates the expression of a mitochondrial enzyme (cytochrome p450 11A1–CYP11A1) that catalyzes the synthesis of pregnenolone, i.e., the substrate for all known steroids [80, 81]. Irregularities in glucocorticoids (a subtype of steroids) have been repeatedly reported in MDD patients [82], being modulated by ECT as robust rises in cortisol have been found following a single ECT session [83,84,85].

Despite this extensive body of evidence linking the latter mentioned genes to either MDD or ECT, several questions remain. In fact, the role of these particular genes within immune cells (the sample type, we obtained our measures from) is largely unclear: Regarding TRERF1, glucocorticoids are well known for their immunoregulatory properties, mediating diverse effects depending on their concentration [86, 87]. Ubiquitin and autophagy affect immunological processes at multiple points: Ubiquitin regulates signaling cascades involved in the activation of NF-кB, and thus, the subsequent production of pro-inflammatory cytokines [88], whereas autophagy has been reported to mediate anti-inflammatory functions by clearing accumulating proteins, apoptotic bodies, and pathogens [73]. Nevertheless, despite being generally involved in ubiquitin-linkage (RNF213 and RNF175) and the negative regulation of autophagy (TBC1D14), it is largely uncertain whether these particular genes directly contribute to the latter mechanisms in immune cells as well, although indicated by some sources [89, 90].

A clear interpretation of the data is further hindered by the insufficient knowledge of the genetic regulation of the proposed candidate genes. In fact, the consequence of the DNAm differences or changes in the expression of the proposed genes is unknown and can only be assumed on the basis of the current literature. In this context, most of the significant CpG sites are located within introns, i.e., regions which were found to be generally low in CpGs. If prevalent, DNAm of these CpGs was suggested to modulate alternative splicing [91, 92], to suppress transposable elements [93], or to regulate the usage of alternative promoters [91, 94]. Hence, DNAm of these specific loci does not ultimately indicate the inhibition of gene expression (as it has been suggested for promoter regions [91]) but might exert various roles. Some significant CpG sites were also found outside gene bodies, i.e., within promoter flanking regions, but the effect on gene expression is also rather unpredictable at these loci.

We furthermore cannot estimate the contribution of the SNPs located at our candidate loci that are either generating or removing the CpGs and thus DNAm as a consequence. Importantly, recent studies propose DNAm to interact with its underlying genotype (even if the respective CpG sites and SNPs are far apart) [95, 96] and the interaction of genetic and environmental factors to be particularly relevant for disease risk [3, 33, 95]. In the case of our SNPs, their minor allele frequency (MAF) values (see Table 2), indicating the second most common variant at a defined locus, are generally relatively small and the removal of CpGs at these particular sites therefore rather unlikely, but not entirely out of question.

The interpretation of the biological significance of our results is further restricted by the sparsity of studies using the TruSeq Methyl Capture EPIC Library Kit for DNAm analyses. Instead, most researchers use either the Illumina Infinium Human Methylation 450 K or the Infinium MethylationEPIC 850 K BeadChip for their experiments. Because both microarrays do not cover our proposed candidate CpG sites, a comparison to the DNAm of other MDD cohorts or healthy subjects is, unfortunately, unfeasible at these particular loci. Consequently, we are also unable to assess the influence of the patients’ clinical characteristics on DNAm other than in our cohort. In this regard, we found a tendency for the DNAm of chr16:19488803 (TMC5) to be influenced by the patient’s respective current episode duration. Since smoking behavior has also differed greatly between ECT responders (five smokers) and non-responders (one smoker), evaluating its influence on DNAm levels would have been of interest but was not feasible due to the small sample size. Studying the literature regarding its general influence on DNAm (not looking at our particular loci of interest only), we found TRERF1’s DNAm (chr6:42219847) to be affected by tobacco intake [97]. Altogether, we thus cannot entirely rule out a possible impact of the patients’ clinical characteristics on our DNAm analysis.

Further, as MDD is primarily thought to be caused by a malfunctioning of neuronal processes, analyzing the correlation between the DNAm in the brain and the periphery would have been of great value. Unfortunately, the web tools frequently used for this investigation (as BECon [98] or BloodBrain [99]) are based on data obtained from the Illumina 450 K array and are thus not applicable to our results. However, as all of our candidate CpGs (except one) do overlap with a known SNP, the observed DNAm difference between ECT responders and non-responders might be a result of their genotype rather than being a pure epigenetic difference. The extent of their DNAm difference and the stability of its pattern during the treatment course supports this idea. In this case, a high correlation between blood and brain DNAm would be rather likely.

Despite these limitations—with the small group size being the most restricting one—our study might be of great value for future approaches, as data generated by the TruSeq Methyl Capture EPIC Kit is still a rarity, especially in the context of MDD.


DNAm of 13 single CpG sites (located within ten genes encoding for either a protein or a long-coding RNA transcript) was found to differ between ECT responder groups or to alter within the treatment course of ECT. The data of the current work thus provides a deeper insight into ECT-associated effects and suggests novel candidate genes for ECT response prediction. Due to a small sample size, the findings must be regarded as preliminary; a replication in larger cohorts is required.


Study design

Our cohort of depressed ECT patients (n = 17) was acquired at the Department of Psychiatry, Social Psychiatry and Psychotherapy at the Hannover Medical School (Germany). The study complies with the ethical principles of the Declaration of Helsinki (1964, including its later amendments) and was approved by the Ethics Committee of the Hannover Medical School (NEKTOR-Registry: 2842-2015). Written informed consent was signed by all participants prior to study inclusion. As this is a naturalistic long-term observational study of one patient cohort, participants were grouped according to their clinical course (to either ECT responder or non-responder) after treatment completion.


MDD diagnosis was established using the International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD-10) and depression severity assessed via two psychometric questionnaires, namely the Beck Depression Inventory (BDI-II) and the Montgomery-Åsberg Depression Rating Scale (MADRS). However, only the latter test served for the assignment of patients to clinical outcome groups. In this context, a decrease of ≥ 50% in MADRS scores was interpreted as a treatment response. The Mini-Mental State Examination (MMSE) was conducted at the same time points (i.e., at baseline and after the first and last ECT). Patients who had an autoimmune, infectious, or schizophrenic disorder were excluded from our study. Heightened levels of CRP, a prominent leukocytosis, or medication with immunomodulatory drugs were additional exclusion criteria.

Application of ECT and sample collection

During the actual treatment course, ECT was applied three times weekly for up to 4 weeks, followed by maintenance ECTs applied only once a week. Right unilateral electrical stimulation was performed using an ultra-brief impulse device (Thymatron® System IV, Somatics, LLC). The seizure threshold was assessed at the first ECT session (based on an age-dependent method), and the stimulus intensity adjusted according to the recorded motoric and electroencephalographic (EEG) seizure duration. If the patient did not show any improvement of symptoms after two following weeks of treatment, bilateral electrode placement was considered. During ECT, patients were anesthetized with methohexital and remifentanil while muscle relaxation was achieved with succinylcholine. Fasting blood samples were taken at four different time points, namely directly before (i.e., 8 a.m.–10 a.m.) and 15 min after the first and the last ECT session, respectively. Samples were stored at 4 °C until further processing (3 hours maximum).

Sample processing

PBMCs—isolation and thawing

PBMCs were isolated by gradient centrifugation as described elsewhere [100]. Based on the recommendations of Mallone et al. [101], changes have been applied to the latter procedure. After isolation, PBMCs were kept at − 196 °C until thawing. Thawing was performed according to a protocol published by the Helmholtz Institute in Munich [102], though (as in the previous case) adaptions were made to meet our requirements. A detailed description of all steps performed to isolate, freeze, and thaw PBMCs is provided in the supplements.

HeLa—thawing and splitting

HeLa cells (immortal cervical cancer cells) were required as a quality control for sequencing. For this purpose, early passage cells were thawed and split according to a protocol published elsewhere [103].

DNA isolation

Genomic DNA (gDNA) of PBMCs and HeLa cells was isolated using the AllPrep DNA/RNA 96 kit (#80311; QIAGEN N.V.). Minor changes have been made to the recommended procedure; a detailed description is to be found in the supplements. After isolation, gDNA of patients was kept at − 80 °C and gDNA of HeLa cells at − 4 °C (being supplemented with 0.5 mM UltraPureTM EDTA (#11568896; Invitrogen AG) in the latter case; storage duration: maximum 3 weeks).

Library generation, quality control, and quantification

Five hundred nanograms of total gDNA per sample was pipetted into microTUBE AFA Fiber Pre-Slit Snap-Cap 6x16mm tubes (#520045; Covaris, Inc.) and subsequently sheared using a Covaris S220 Ultrasonicator. The sheared gDNA was utilized as input for preparing targeted methylseq libraries with the TruSeq-Methyl Capture EPIC Library Kit (#FC-151-1003; Illumina, Inc.), allowing a preparation of up to 48 libraries at four-plex within less than 2 days. All steps were performed as recommended in the Illumina user document 1000000001643 v01 May 2017, though one additional purification step was introduced at the end of the standard procedure, using 1× Agencourt® AMPure® XP Beads (#A63881; Beckman Coulter, Inc.). Four-plex DNA samples were barcoded by a single indexing (6 bp) approach using 12 different DNA indexes. All generated DNA libraries were amplified with 11–13 cycles of PCR using a KAPA HiFi HotStart Uracil+ Ready Mix (2X) enzyme (#KK2801; Kapa Biosystems), which was not included in the kit. Fragment length distribution of individual libraries was monitored using the Bioanalyzer High Sensitivity DNA Assay (#5067-4626; Agilent Technologies). Quantification of libraries was performed by the use of the Qubit® dsDNA HS Assay Kit (#Q32854; ThermoFisher Scientific). Importantly, gDNA of HeLa cells served as a quality control and was thus added to each run of library preparation and subsequent sequencing (1–2 aliquots per run).

Library denaturation and sequencing run

Equimolar amounts of twelve individually barcoded libraries were pooled. Each analyzed library relevant to this project constitutes 8.7% of an overall flow cell capacity. The library pool was denatured with sodium chloride (NaOH, #72082-100 ml; SIGMA-ALDRICH Co.) and was finally diluted to 1.8 pM according to the Denature and Dilute Libraries Guide (Document #15048776 v02; Illumina, Inc.). 1.3 ml of the denatured pool was loaded on an Illumina NextSeq 550 Sequencer using a High Output Flow Cell for paired-end reads (Document #20024907; Illumina, Inc.). Paired-end sequencing was performed with 76 cycles, a 6-base barcode index, and 25% calibration control v3 PhiX library (#FC-110-3001; Illumina, Inc.). This level of PhiX was required as the samples were relatively GC rich. In total, 8 NextSeq runs were performed.

Data processing

Sequence data analysis

Illumina 75 bp paired-end datasets were demultiplexed using bcl2fastq version (Illumina, Inc.). Fastq files were then subjected to quality control with FASTQC and MultiQC [104] and analyzed using the nf-core methylseq pipeline with the Bismark software (version 1.5dev) [105, 106]. The genome reference hg38 from Ensembl was used without decoy sequences. The pipeline was modified to only make methylation calls for sites covered by at least 5 reads, instead of the default 1 read. Coverage and methylation calls were converted into bigwig format and visualized in the JBrowse web application [35]. All data were analyzed on the MHH HPC-seq SLURM research cluster. Where several runs were necessary to achieve sufficient coverage, FASTQ files were combined before analysis.

Quality control

In order to control for technical variability (which was necessary as our experiment included several runs of library preparation and sequencing), a detailed quality control of the measurement was carried out. In this context, isolated gDNA from HeLa cells was added to each sequencing run. The resulting average HeLa-cell-probe correlation (whose calculation included all CpG sites of all measurements with a previously defined minimum coverage of least 5 reads) was Mr = 0.97 and based on 14 measured samples and 2016707 observed CpG probes each, indicating a sufficiently good accuracy of our measurements. Furthermore, the measured samples were examined with density plots and dendrograms. One patient sample showed a conspicuous value distribution pattern and was different from all other samples in the cluster analysis. Therefore, a faulty measurement was assumed, leading to the subsequent exclusion of this sample from further analyses. In order to avoid gender bias, chromosome X and Y were additionally excluded. After all quality control measures, 1476812 CpG probes (per sample) were examined for final statistics.

Statistical analyses

Demographics of patients were normally distributed. T tests and Fisher’s exact tests were used for the analysis of demographic, psychometric, and other clinical baseline differences between ECT responders and non-responders. Regarding the methylation analysis of our cohort, we first checked for ECT-associated changes considering the overall DNAm levels (=mean of all measured CpG sites), being investigated with repeated measures ANOVA modeling approach implemented in the lmer package [107] and separated for “response” as an outcome definition. Second, we performed a detailed analysis for the differences in methylation with respect to response/non-response for every probe (DMP) with a series of repeated measures ANOVA. As a correction for multiple testing, a false discovery rate (FDR) of < 5% in combination with a minimum variance of 0.1 was defined as significant. Only samples with complete data at all four time points (n = 12) were included in the analyses. Statistical analyses of the patients’ clinical baseline characteristics were performed using IBM SPSS Statistics 25.0 for Windows (IBM Corp.), the methylation analyses were conducted within the R (3.6.1) environment on Windows 10.0.18362.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request because they include genetic information of patients.



Aquaporin 10


Beck Depression Inventory


Brain-derived neurotrophic factor


Body mass index




Cytosine-phosphate-guanine dinucleotide


C-reactive protein


Cytochrome p450 11A1


Differences in methylation for every probe


DNA methylation


Dopamine receptor D2


Electroconvulsive stimulation


Electroconvulsive therapy




False discovery rate


Genomic DNA




International Statistical Classification of Diseases and Related Health Problems 10th Revision


Light chain 3-II/I


Montgomery-Åsberg Depression Rating Scale


Minor allele frequency


Major depressive disorder


Moyamoya disease


Mini-Mental State Examination


Sodium chloride


Peripheral blood mononuclear cells


Ring finger protein 175/213


Standard deviation


Single-nucleotide polymorphism


TBC1 domain family member 14


Transmembrane channel-like 5


Transcriptional regulating factor 1


Vascular endothelial growth factor


WSC domain containing 1


  1. 1.

    WHO, World Health Organization. Depression and other common mental disorders: global health estimates: World Health Organization; 2017. p. 1–24.

  2. 2.

    Stepniak B, Papiol S, Hammer C, Ramin A, Everts S, Hennig L, et al. Accumulated environmental risk determining age at schizophrenia onset: a deep phenotyping-based study. Lancet Psychiatry. 2014;1(6):444–53.

    Article  PubMed  Google Scholar 

  3. 3.

    Czamara D, Eraslan G, Page CM, Lahti J, Lahti-Pulkkinen M, Hämäläinen E, et al. Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns. Nat Commun. 2019;10(1):2548.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Jaenisch R, Bird A. Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003;33:245–54.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Richards EJ. Inherited epigenetic variation—revisiting soft inheritance. Nat Rev Genet. 2006;7(5):395–401.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Klengel T, Mehta D, Anacker C, Rex-Haffner M, Pruessner JC, Pariante CM, et al. Allele-specific FKBP5 DNA demethylation mediates gene–childhood trauma interactions. Nat Neurosci. 2013;16:33–41.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Moore DS. Behavioral epigenetics. Wiley Interdiscip Rev Syst Biol Med. 2017;9(1):e1333.

    Article  Google Scholar 

  8. 8.

    Weaver ICG, Cervoni N, Champagne FA, D’Alessio AC, Sharma S, Seckl JR, et al. Epigenetic programming by maternal behavior. Nat Neurosci. 2004;7(8):847–54.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Franklin TB, Russig H, Weiss IC, Grff J, Linder N, Michalon A, et al. Epigenetic transmission of the impact of early stress across generations. Biol Psychiatry. 2010;68(5):408–15.

    Article  PubMed  Google Scholar 

  10. 10.

    Binder EB, Salyakina D, Lichtner P, Wochnik GM, Ising M, Pütz B, et al. Polymorphisms in FKBP5 are associated with increased recurrence of depressive episodes and rapid response to antidepressant treatment. Nat Genet. 2004;36(12):1319–25.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Ising M, Horstmann S, Kloiber S, Lucae S, Binder EB, Kern N, et al. Combined dexamethasone/corticotropin releasing hormone test predicts treatment response in major depression–a potential biomarker? Biol Psychiatry. 2007;62(1):47–54.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Carvalho LA, Bergink V, Sumaski L, Wijkhuijs J, Hoogendijk WJ, Birkenhager TK, et al. Inflammatory activation is associated with a reduced glucocorticoid receptor alpha/beta expression ratio in monocytes of inpatients with melancholic major depressive disorder. Transl Psychiatry. 2014;4(1):e344.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Fuchikami M, Morinobu S, Segawa M, Okamoto Y, Yamawaki S, Ozaki N, et al. DNA methylation profiles of the brain-derived neurotrophic factor (BDNF) gene as a potent diagnostic biomarker in major depression. PLoS One. 2011;6(8):e23881.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Chen B, Dowlatshahi D, MacQueen GM, Wang J-F, Young LT. Increased hippocampal bdnf immunoreactivity in subjects treated with antidepressant medication. Biol Psychiatry. 2001;50(4):260–5.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Molendijk ML, Spinhoven P, Polak M, Bus BAA, Penninx BWJH, Elzinga BM. Serum BDNF concentrations as peripheral manifestations of depression: evidence from a systematic review and meta-analyses on 179 associations (N = 9484). Mol Psychiatry. 2014;19(7):791–800.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Mariga A, Mitre M, Chao MV. Consequences of brain-derived neurotrophic factor withdrawal in CNS neurons and implications in disease. Neurobiol Dis. 2017;97:73–9.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Gaynes BN, Rush AJ, Trivedi MH, Wisniewski SR, Spencer D, Fava M. The STAR*D study: treating depression in the real world. Cleve Clin J Med. 2008;75(1):57–66.

  18. 18.

    Lisanby SH. Electroconvulsive therapy for depression. 2007;357(19):1939–45.

  19. 19.

    Heijnen WT, Birkenhäger TK, Wierdsma AI, van den Broek WW. Antidepressant pharmacotherapy failure and response to subsequent electroconvulsive therapy: a meta-analysis. J Clin Psychopharmacol. 2010;30(5):616–9.

    Article  PubMed  Google Scholar 

  20. 20.

    Prudic J, Olfson M, Marcus SC, Fuller RB, Sackeim HA. Effectiveness of electroconvulsive therapy in community settings. Biol Psychiatry. 2004;55(3):301–12.

    Article  PubMed  Google Scholar 

  21. 21.

    Brus O, Cao Y, Gustafsson E, Hultén M, Landen M, Lundberg J, et al. Self-assessed remission rates after electroconvulsive therapy of depressive disorders. Eur Psychiatry. 2017;45:154–60.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    O’Connor MK, Knapp R, Husain M, Rummans TA, Petrides G, Smith G, et al. The influence of age on the response of major depression to electroconvulsive therapy: a C.O.R.E. report. Am J Geriatr Psychiatry. 2001;9(4):382–90.

    Article  PubMed  Google Scholar 

  23. 23.

    Petrides G, Fink M, Husain MM, Knapp RG, John Rush A, Mueller M, et al. ECT remission rates in psychotic versus nonpsychotic depressed patients: a report from CORE. J ECT. 2001;17(4):244–53.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Van Diermen L, Van Den Ameele S, Kamperman AM, Sabbe BCG, Vermeulen T, Birkenhäger TK. Prediction of electroconvulsive therapy response and remission in major depression: meta-analysis. Br J Psychiatry. 2018;212(2):71–80.

  25. 25.

    Jentsch MC, Van Buel EM, Bosker FJ, Gladkevich AV, Klein HC, Oude Voshaar RC, et al. Biomarker approaches in major depressive disorder evaluated in the context of current hypotheses. Biomark Med. 2015;9(3):277–97.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    SNPedia. Accessed 12 Feb 2020.

  27. 27.

    Anttila S, Huuhka K, Huuhka M, Illi A, Rontu R, Leinonen E, et al. Catechol-O-methyltransferase (COMT) polymorphisms predict treatment response in electroconvulsive therapy. Pharm J. 2008;8(2):113–6.

    CAS  Article  Google Scholar 

  28. 28.

    Lin Z, He H, Zhang C, Wang Z, Jiang M, Li Q, et al. Influence of Val108/158Met COMT gene polymorphism on the efficacy of modified electroconvulsive therapy in patients with treatment resistant depression. Cell Biochem Biophys. 2015;71(3):1387–93.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Huuhka K, Anttila S, Huuhka M, Hietala J, Huhtala H, Mononen N, et al. Dopamine 2 receptor C957T and catechol-o-methyltransferase Val158Met polymorphisms are associated with treatment response in electroconvulsive therapy. Neurosci Lett. 2008;448(1):79–83.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Bousman CA, Katalinic N, Martin DM, Smith DJ, Ingram A, Dowling N, et al. Effects of COMT, DRD2, BDNF, and APOE genotypic variation on treatment efficacy and cognitive side effects of electroconvulsive therapy. J ECT. 2015;31(2):129–35.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Tsai S-J. The P11, tPA/plasminogen system and brain-derived neurotrophic factor: implications for the pathogenesis of major depression and the therapeutic mechanism of antidepressants. Med Hypotheses. 2007;68(1):180–3.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Neyazi A, Theilmann W, Brandt C, Rantamäki T, Matsui N, Rhein M, et al. P11 promoter methylation predicts the antidepressant effect of electroconvulsive therapy. Transl Psychiatry. 2018;8(1):25.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Lohoff FW. Overview of the genetics of major depressive disorder. Curr Psychiatry Rep. 2010;12(6):539–46.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Ensembl. Accessed 03 Dec 2019.

  35. 35.

    Buels R, Yao E, Diesh CM, Hayes RD, Munoz-Torres M, Helt G, et al. JBrowse: a dynamic web platform for genome visualization and analysis. Genome Biol. 2016;17(1).

  36. 36.

    Genecards. Accessed 08 Dec 2019.

  37. 37.

    GWAS Catalog. The NHGRI-EBI Catalog of published genome-wide association studies. Accessed 08 Dec 2019.

  38. 38.

    Byrne EM, Carrillo-Roa T, Henders AK, Bowdler L, McRae AF, Heath AC, et al. Monozygotic twins affected with major depressive disorder have greater variance in methylation than their unaffected co-twin. Transl Psychiatry. 2013;3(6):e269.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Dempster EL, Wong CC, Lester KJ, Burrage J, Gregory AM, Mill J, et al. Genome-wide methylomic analysis of monozygotic twins discordant for adolescent depression. Biol Psychiatry. 2014;76(12):977–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Laity JH, Lee BM, Wright PE. Zinc finger proteins: new insights into structural and functional diversity. Curr Opin Struct Biol. 2001;11(1):39–46.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Morito D, Nishikawa K, Hoseki J, Kitamura A, Kotani Y, Kiso K, et al. Moyamoya disease-associated protein mysterin/RNF213 is a novel AAA+ ATPase, which dynamically changes its oligomeric state. Sci Rep. 2014;4:4442.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Ciechanover A, Elias S, Heller H, Ferber S, Hershko A. Characterization of the heat-stable polypeptide of the ATP-dependent proteolytic system from reticulocytes. J Biol Chem. 1980;255(16):7525–8.

    CAS  PubMed  Google Scholar 

  43. 43.

    Hershko A, Eytan E, Ciechanover A, Haas AL. Immunochemical analysis of the turnover of ubiquitin-protein conjugates in intact cells. Relationship to the breakdown of abnormal proteins. J Biol Chem. 1982;257(23):13964–70.

    CAS  PubMed  Google Scholar 

  44. 44.

    Bach SV, Hegde AN. The proteasome and epigenetics: zooming in on histone modifications. Biomol Concepts. 2016;7(4):215–27.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Scholz B, Korn C, Wojtarowicz J, Mogler C, Augustin I, Boutros M, et al. Endothelial RSPO3 controls vascular stability and pruning through non-canonical WNT/Ca2+/NFAT signaling. Dev Cell. 2016;36(1):79–93.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Kamada F, Aoki Y, Narisawa A, Abe Y, Komatsuzaki S, Kikuchi A, et al. A genome-wide association study identifies RNF213 as the first Moyamoya disease gene. J Hum Genet. 2011;56(1):34–40.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Nakayama C, Mikami T, Ukai R, Chiba R, Enatsu R, Suzuki H, et al. Moyamoya disease with epileptic nystagmus: a case report. J Clin Neurosci. 2019;70:251–4.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Berent D, Macander M, Szemraj J, Orzechowska A, Gałecki P. Vascular endothelial growth factor A gene expression level is higher in patients with major depressive disorder and not affected by cigarette smoking, hyperlipidemia or treatment with statins. Acta Neurobiol Exp. 2014;74:82–90.

    Google Scholar 

  49. 49.

    Hellsten J, West MJ, Arvidsson A, Ekstrand J, Jansson L, Wennström M, et al. Electroconvulsive seizures induce angiogenesis in adult rat hippocampus. Biol Psychiatry. 2005;58(11):871–8.

    Article  PubMed  Google Scholar 

  50. 50.

    Van Den Bossche MJA, Emsell L, Dols A, Vansteelandt K, De Winter F-L, Van den Stock J, et al. Hippocampal volume change following ECT is mediated by rs699947 in the promotor region of VEGF. Transl Psychiatry. 2019;9(1):191.

    CAS  Article  Google Scholar 

  51. 51.

    Fischler B, Bocken R, Schneider I, De Waele M, Thielemans K, Derde MP. Immune changes induced by electroconvulsive therapy (ECT). Ann N Y Acad Sci. 1992;650(1):326–30.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Kronfol Z, Nair MP, Weinberg V, Young EA, Aziz M. Acute effects of electroconvulsive therapy on lymphocyte natural killer cell activity in patients with major depression. J Affect Disord. 2002;71(1):211–5.

    Article  PubMed  Google Scholar 

  53. 53.

    Lehtimäki K, Keränen T, Huuhka M, Palmio J, Hurme M, Leinonen E, et al. Increase in plasma proinflammatory cytokines after electroconvulsive therapy in patients with depressive disorder. J ECT. 2008;24(1):88–91.

    Article  PubMed  Google Scholar 

  54. 54.

    Fluitman SBAHA, Heijnen CJ, Denys DAJP, Nolen WA, Balk FJ, Westenberg HGM. Electroconvulsive therapy has acute immunological and neuroendocrine effects in patients with major depressive disorder. J Affect Disord. 2011;131(1–3):388–92.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Hestad KA, Tønseth S, Støen CD, Ueland T, Aukrust P. Raised plasma levels of tumor necrosis factor α in patients with depression: normalization during electroconvulsive therapy. J ECT. 2003;19(4):183–8.

  56. 56.

    van Buel EM, Patas K, Peters M, Bosker FJ, Eisel ULM, Klein HC. Immune and neurotrophin stimulation by electroconvulsive therapy: is some inflammation needed after all? Transl Psychiatry. 2015;5(7):e609.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Ohkubo K, Sakai Y, Inoue H, Akamine S, Ishizaki Y, Matsushita Y, et al. Moyamoya disease susceptibility gene RNF213 links inflammatory and angiogenic signals in endothelial cells. Sci Rep. 2015;5(1):13191.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Kanoke A, Fujimura M, Niizuma K, Fujimura T, Kakizaki A, Ito A, et al. Temporal profile of magnetic resonance angiography and decreased ratio of regulatory T cells after immunological adjuvant administration to mice lacking RNF213, a susceptibility gene for moyamoya disease. Brain Res. 2016;1642:1–9.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Grosse L, Carvalho LA, Birkenhager TK, Hoogendijk WJ, Kushner SA, Drexhage HA, et al. Circulating cytotoxic T cells and natural killer cells as potential predictors for antidepressant response in melancholic depression. Restoration of T regulatory cell populations after antidepressant therapy. Psychopharmacology. 2016;233(9):1679–88.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Zhao S, Gong Z, Zhang J, Xu X, Liu P, Guan W, et al. Elevated serum microRNA Let-7c in moyamoya disease. J Stroke Cerebrovasc Dis. 2015;24(8):1709–14.

    Article  PubMed  Google Scholar 

  61. 61.

    Gururajan A, Naughton ME, Scott KA, O’Connor RM, Moloney G, Clarke G, et al. MicroRNAs as biomarkers for major depression: a role for let-7b and let-7c. Transl Psychiatry. 2016;6(8):e862.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Postmus I, Trompet S, Deshmukh HA, Barnes MR, Li X, Warren HR, et al. Pharmacogenetic meta-analysis of genome-wide association studies of LDL cholesterol response to statins. Nat Commun. 2014;5(1):5068.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Salagre E, Fernandes BS, Dodd S, Brownstein DJ, Berk M. Statins for the treatment of depression: a meta-analysis of randomized, double-blind, placebo-controlled trials. J Affect Disord. 2016;200:235–42.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Kaneko M, Iwase I, Yamasaki Y, Takai T, Wu Y, Kanemoto S, et al. Genome-wide identification and gene expression profiling of ubiquitin ligases for endoplasmic reticulum protein degradation. Sci Rep. 2016;6(August):1–10.

    CAS  Article  Google Scholar 

  65. 65.

    Sabunciyan S, Aryee MJ, Irizarry RA, Rongione M, Webster MJ, Kaufman WE, et al. Genome-wide DNA methylation scan in major depressive disorder. PLoS One. 2012;7(4):e34451.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Nagy C, Suderman M, Yang J, Szyf M, Mechawar N, Ernst C, et al. Astrocytic abnormalities and global DNA methylation patterns in depression and suicide. Mol Psychiatry. 2015;20(3):320–8.

  67. 67.

    Peng H, Yang J, Li G, You Q, Han W, Li T, et al. Ubiquitylation of p62/sequestosome1 activates its autophagy receptor function and controls selective autophagy upon ubiquitin stress. Cell Res. 2017;27(5):657–74.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Longatti A, Tooze SA. Recycling endosomes contribute to autophagosome formation. Autophagy. 2012;8(11):1682–3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Gassen NC, Rein T. Is there a role of autophagy in depression and antidepressant action? Front Psychiatry. 2019;10:337.

    Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Alcocer-Gómez E, Casas-Barquero N, Núñez-Vasco J, Navarro-Pando J, Bullon P. Psychological status in depressive patients correlates with metabolic gene expression. CNS Neurosci Ther. 2017;23(10):843–5.

    Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Li P, Hao XC, Luo J, Lv F, Wei K, Min S. Propofol mitigates learning and memory impairment after electroconvulsive shock in depressed rats by inhibiting autophagy in the hippocampus. Med Sci Monit. 2016;22:1702–8.

  72. 72.

    Yang Z, Huang J, Geng J, Nair U, Klionsky DJ. Atg22 recycles amino acids to link the degradative and recycling functions of autophagy. Mol Biol Cell. 2006;17(12):5094–104.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  73. 73.

    Ma Y, Galluzzi L, Zitvogel L, Kroemer G. Review autophagy and cellular immune responses. Immunity. 2013;39(2):211–27.

    CAS  Article  PubMed  Google Scholar 

  74. 74.

    McCormick LM, Boles Ponto LL, Pierson RK, Johnson HJ, Magnotta V, Brumm MC. Metabolic correlates of antidepressant and antipsychotic response in patients with psychotic depression undergoing electroconvulsive therapy. J ECT. 2007;23(4):265–73.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Gormley P, Anttila V, Winsvold BS, Palta P, Esko T, Pers TH, et al. Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine. Nat Genet. 2016;48(8):856–66.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  76. 76.

    Yang Y, Ligthart L, Terwindt GM, Boomsma DI, Rodriguez-Acevedo AJ, Nyholt DR. Genetic epidemiology of migraine and depression. Cephalalgia. 2016;36(7):679–91.

    CAS  Article  PubMed  Google Scholar 

  77. 77.

    Byrne EM, Gehrman PR, Medland SE, Nyholt DR, Heath AC, Madden PAF, et al. A genome-wide association study of sleep habits and insomnia. Am J Med Genet Part B Neuropsychiatr Genet. 2013;162(5):439–51.

    CAS  Article  Google Scholar 

  78. 78.

    O’Brien EM, Chelminski I, Young D, Dalrymple K, Hrabosky J, Zimmerman M. Severe insomnia is associated with more severe presentation and greater functional deficits in depression. J Psychiatr Res. 2011;45(8):1101–5.

    Article  PubMed  Google Scholar 

  79. 79.

    Gotfryd K, Mósca AF, Missel JW, Truelsen SF, Wang K, Spulber M, et al. Human adipose glycerol flux is regulated by a pH gate in AQP10. Nat Commun. 2018;9(1):4749.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  80. 80.

    Gizard F, El-Alfy M, Duguay Y, Lavallée B, DeWitte F, Staels B, et al. Function of the transcriptional regulating protein of 132 kDa (TReP-132) on human P450scc gene expression. Endocr Res. 2002;28(4):559–74.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Gizard F, Teissier E, Dufort I, Luc G, Luu-The V, Staels B, et al. The transcriptional regulating protein of 132 kDa (TReP-132) differentially influences steroidogenic pathways in human adrenal NCI-H295 cells. J Mol Endocrinol. 2004;32(2):557–69.

    CAS  Article  PubMed  Google Scholar 

  82. 82.

    Burke HM, Davis MC, Otte C, Mohr DC. Depression and cortisol responses to psychological stress: a meta-analysis. Psychoneuroendocrinology. 2005;30(9):846–56.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Deakin JFW, Ferrier IN, Crow TJ, Johnstoneand EC, Lawler P. Effects of ECT on pituitary hormone release: relationship to seizure, clinical variables and outcome. Brit J Psychiat. 1983;143(6):618–24.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Aperia B, Thorén M, Zettergren M, Wetterberg L. Plasma pattern of adrenocorticotropin and cortisol during electroconvulsive therapy in patients with major depressive illness. Acta Psychiatr Scand. 1984;70(4):361–9.

    CAS  Article  PubMed  Google Scholar 

  85. 85.

    Kronfol Z, Hamdan-Allen G, Goel K, Hill EM. Effects of single and repeated electroconvulsive therapy sessions on plasma ACTH, prolactin, growth hormone and cortisol concentrations. Psychoneuroendocrinology. 1991;16(4):345–52.

    CAS  Article  PubMed  Google Scholar 

  86. 86.

    Talabér G, Jondal M, Okret S. Extra-adrenal glucocorticoid synthesis: Immune regulation and aspects on local organ homeostasis. Mol Cell Endocrinol. 2013;380(1):89–98.

    CAS  Article  PubMed  Google Scholar 

  87. 87.

    Cain DW, Cidlowski JA. Immune regulation by glucocorticoids. Nat Rev Immunol. 2017;17(4):233–47.

    CAS  Article  PubMed  Google Scholar 

  88. 88.

    Hu H, Sun SC. Ubiquitin signaling in immune responses. Cell Res. 2016;26(4):457–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  89. 89.

    Reactome. Accessed 10 Dec 2019.

  90. 90.

    Reactome.,R-HSA-1280218. Accessed 10 Dec 2019.

  91. 91.

    Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.

    CAS  Article  PubMed  Google Scholar 

  92. 92.

    Shayevitch R, Askayo D, Keydar I, Ast G. The importance of DNA methylation of exons on alternative splicing. RNA. 2018;24(10):1351–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  93. 93.

    Yoder JA, Walsh CP, Bestor TH. Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 1997;13(8):335–40.

    CAS  Article  PubMed  Google Scholar 

  94. 94.

    Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D’Souza C, Fouse SD, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466(7303):253–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  95. 95.

    Bell AF, Carter CS, Steer CD, Golding J, Davis JM, Steffen AD, et al. Interaction between oxytocin receptor DNA methylation and genotype is associated with risk of postpartum depression in women without depression in pregnancy. Front Genet. 2015;6.

  96. 96.

    Fisher VA, Wang L, Deng X, Sarnowski C, Cupples LA, Liu CT. Do changes in DNA methylation mediate or interact with SNP variation? A pharmacoepigenetic analysis. BMC Genet. 2018;19(1):70.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  97. 97.

    Zeilinger S, Kühnel B, Klopp N, Baurecht H, Kleinschmidt A, Gieger C, et al. Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS One. 2013;8(5):e63812.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  98. 98.

    Edgar RD, Jones MJ, Meaney MJ, Turecki G, Kobor MS. BECon: a tool for interpreting DNA methylation findings from blood in the context of brain. Transl Psychiatry. 2017;7(8):e1187.

  99. 99.

    Hannon E, Lunnon K, Schalkwyk L, Mill J. Interindividual methylomic variation across blood, cortex, and cerebellum: implications for epigenetic studies of neurological and neuropsychiatric phenotypes. Epigenetics. 2015;10(11):1024–32.

  100. 100.

    Sieve I, Ricke-Hoch M, Kasten M, Battmer K, Stapel B, Falk CS, et al. A positive feedback loop between IL-1β, LPS and NEU1 may promote atherosclerosis by enhancing a pro-inflammatory state in monocytes and macrophages. Vasc Pharmacol. 2018;103–105:16–28.

    CAS  Article  Google Scholar 

  101. 101.

    Mallone R, Mannering SI, Brooks-Worrell BM, Durinovic-Belló I, Cilio CM, Wong FS, et al. Isolation and preservation of peripheral blood mononuclear cells for analysis of islet antigen-reactive T cell responses: position statement of the T-Cell Workshop Committee of the Immunology of Diabetes Society. Clin Exp Immunol. 2011;163(1):33–49.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  102. 102.

    Pohla, H. Helmholtz Zentrum Munich CCG immune monitoring protocol: PBMC isolation, cryopreservation and thawing. Version 1.0. 2007. Accessed 03 Oct 2019.

  103. 103.

    Woods Hole Physiology Course. 2006.,%20HeLa%20and%20CHO%20cells.pdf?sequence=1; Accessed 09 Dec 2019.

  104. 104.

    Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: Summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  105. 105.

    Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  106. 106.

    Ewels P, Peltzer A, Fillinger S, Alneberg J, Patel H, Wilm A, et al. nf-core: community curated bioinformatics pipelines. BioRxiv. 2019;610741.

  107. 107.

    Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.

    Article  Google Scholar 

Download references


The authors are sincerely grateful to Dr. Britta Stapel and Dr. Abdul Qayyum Khan for their scientific advice. Further, we owe our profound gratitude to the patients that participated in our study. Last but not least, we thank Hannover Medical School for financial support.


The research was supported by in-house funding from Hannover Medical School. The grant provider had no role in study design, collection, interpretation, or analysis of data nor in writing of the report or the decision to submit the paper for publication.

Author information




NM: Processed the samples. Assisted with statistical analyses and the graphical illustration. Wrote the manuscript (including the supplementary material section) with support from TZ, MD, CD, AN, and HF. TZ: Analyzed the data statistically (incl. the graphical illustration of the results). KJ + MD: Provided enormous support regarding the sample processing. CD: Processed the data (section “sequence data analysis”). LW + SB: Substantial contribution to the conception of the work. Revised the study for important intellectual content. HBM + FE: Recruitment and treatment of the cohort (provided clinical data and blood samples). AN + HF: Supervised the study. An essential contribution to the conception and design of the work, the acquisition of patient groups, and the analysis and interpretation of data. Drafted the work and revised it for important intellectual content. Strong guidance regarding the writing of the paper, and the presentation of data. The authors have approved the final manuscript.

Corresponding author

Correspondence to Nicole Moschny.

Ethics declarations

Ethics approval and consent to participate

The study complies with the ethical principles of the Declaration of Helsinki (1964, including its later amendments) and was approved by the Ethics Committee of the Hannover Medical School (NEKTOR-Registry: 2842-2015). Written informed consent was signed by all participants prior to study inclusion.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Moschny, N., Zindler, T., Jahn, K. et al. Novel candidate genes for ECT response prediction—a pilot study analyzing the DNA methylome of depressed patients receiving electroconvulsive therapy. Clin Epigenet 12, 114 (2020).

Download citation


  • Depression
  • DNA methylation
  • Electroconvulsive therapy
  • EWAS
  • Personalized medicine
  • Response prediction
  • Single-nucleotide polymorphism
  • RNF213
  • Ubiquitin
  • Autophagy