Skip to main content

BDNF exon IV promoter methylation and antidepressant action: a complex interplay



BDNF exon IV promoter methylation is a potential biomarker for treatment response to antidepressants in MDD. We have previously shown CpG-87 methylation as a successful biomarker for the prediction of non-response to monoaminergic antidepressants like the SSRI Fluoxetine or the SNRI Venlafaxine. This study aimed to dissect the biological evidence and mechanisms for the functionality of CpG-87 methylation in a cell culture model.


We observed a significant interaction between methylation and antidepressant-mediated transcriptional activity in BDNF exon IV promoter. In addition, antidepressant treatment increased the promoter methylation in a concentration-dependent manner. Further single CpG methylation of -87 did not change the promoter activity, but methylation of CREB domain CpG-39 increased the transcriptional activity in an antidepressant-dependent manner. Interestingly, DNMT3a overexpression also increases the BDNF exon IV transcription and more so in Venlafaxine-treated cells.


The study strengthens the previously reported association between antidepressant treatment and BDNF exon IV promoter methylation as well as hints toward the mechanism of action. We argue that potential CpG methylation biomarkers display a complex synergy with the molecular changes at the neighboring CpG positions, thus highlighting the importance of epiallele analyses.


Major depressive disorder (MDD) is a multifaceted disease disorder that imposes a measurable burden on society. Over the past decade, it has been widely accepted that symptoms and manifestations of depression between individuals cannot be empirically explained by dysfunction of any particular biological process. Although matching the pathology to a certain extent with the molecules and pathways involved can draw a subjective classification, more evidence is needed to extrapolate findings to the therapeutic regimens [1]. The major challenge in treating the patients with MDD is the inherent heterogeneity in response rates toward the commonly used antidepressants, the serotonin reuptake inhibitors (SSRI) and serotonin–norepinephrine reuptake inhibitors (SNRI) [2].

Therefore, it is imperative to work systematically toward understanding the underlying causative and predictive biological differences in the patients.

Brain-derived neurotrophic factor (BDNF)-mediated signaling cascades are implicated in neuronal differentiation, synaptic plasticity and neurotransmission, claiming BDNF as a central molecule for various interventions for depression as well as other psychiatric disorders. Compelling evidence is available toward the disruption of BDNF gene in manifesting the depressive behaviors from both animal and patient studies [3]. BDNF upregulation is one of the key functional phenomena to downstream pathways of antidepressant drugs [4, 5].

Transcription of BDNF is an intricate example of temporal and spatial regulation. The BDNF gene contains 11 exons with multiple splicing sites [6]. The last exon codes for BDNF protein and contains two poly (A) stretches resulting in a protein with a short or long 3′ UTR. The remaining 10 exons contain independent internal promoters that allow the transcription of 13 different transcripts each containing a unique 5′ UTR proximate to the coding exon [7]. Transcription factors, epigenetic modifiers and RNA-binding proteins work in synchrony to regulate, modulate and fine-tune the expression of specific BDNF transcript in response to developmental cues and also various environmental stimuli [8].

Neuronal activity-dependent calcium influx is important to process and store the information in response to environmental stimuli. BDNF exon I and IV transcription is regulated by calcium-dependent transcription factors, cAMP response element-binding protein (CREB) and upstream stimulating factor (USF) [9]. CREB is a histone deacetylase and recruits other transcriptional activators, upregulating the transcription of BDNF transcripts [10]. In addition, histone modification like the H3K4me3 mark has been associated with active BDNF promoter [11, 12].

Interestingly, transcriptional regulation of BDNF, in absence of any stimuli occurs via CpG methylation or H3K27 trimethylation. Methylation-dependent regulation has been reported for BDNF exon I, IV, VI and IX [13,14,15].

MDD has a mere heredity of 37%, and therefore the onset, development and progression of disease have been shown to be primarily mediated via epigenetic mechanisms [16]. Social factors like early life stress, childhood maltreatment and maternal separation change BDNF promoter methylation [17]. Serum BDNF levels before treatment have been associated with SSRI response in depression [18]. CpG-specific hyper- and hypomethylation of exon I and IV promoters have been reported in various psychiatric disorders including MDD. Although it remains inconclusive whether antidepressants alter BDNF methylation, differences in the CpG methylation of BDNF internal promoters have been presented as a biomarker for therapy response [19, 20].

We have previously shown that the hypermethylation at CpG-87 (from the transcriptional start site) can successfully predict the non-response to antidepressant therapy [21]. This finding has been replicated in a larger cohort of moderate to severely depressed patients, while it could not be replicated for the patients with mild or minor depression [22]. Similar findings have been published in a Chinese cohort [23]. In our first report, we also provided first in vitro data on the transcriptional activity of the BDNF exon IV promoter. Antidepressant incubation with Fluoxetine (SSRI) or Venlafaxine (SNRI) led to reduced activity in a dual-luciferase reporter assay. This down-regulation was specific to unmethylated promoter fragment.

Here, we aim to provide a more detailed insight into the complex interplay of single CpG methylation and antidepressant effects on the BDNF exon IV promoter activity and BDNF expression (Fig. 1).

Fig. 1
figure 1

BDNF exon IV promoter schematics. Sequence of the BDNF exon IV promoter analyzed for methylation. Primers (blue) and the CpGs (red) are highlighted. CpG-39 and -35 encompasses the cAMP response element (CRE). CpGs are numbered with respect to BDNF exon IV promoter


Transcriptional activity of BDNF exon IV correlates with promoter methylation

We cloned BDNF exon IV promoter region in a CpG-free vector upstream of the luciferase gene and analyzed the effect of antidepressants on the transcriptional activity. This experiment replicates the findings from our previous paper [21]. We observed a decrease in the BDNF exon IV transcription of unmethylated promoter (F(1, 50) = 11.35; P = 0.0015) in SH-SY5Y neuroblastoma cell line after treatment with antidepressants. The observed reduction was significant for Venlafaxine, Lithium carbonate, Fluoxetine, and Mirtazapine (P = 0.05). Interestingly, the antidepressant-induced reduction in BDNF exon IV transcription was not seen in methylated promoter, thereby explaining a significant interaction between methylation and antidepressant on the promoter activity (F(4, 50) = 4.059; P = 0.0063). An overall reduction in promoter activity was observed in methylated promoter construct (main effect derived from the two-way ANOVA (F(4, 50) = 5.615; P = 0.0008)) (Fig. 2). In the post hoc analysis, we observed a strong and significant reduction in activity in the methylated promoter fragment without antidepressant incubation (control) as compared to the unmethylated construct (post hoc: t(50) = 3.948; P < 0.01).

Fig. 2
figure 2

a Effect of antidepressants on BDNF exon IV promoter activity. The antidepressants at the given concentration decrease the promoter activity: Venlafaxine (50 µM, P < 0.001), Fluoxetine (10 µM, P < 0.05), Mirtazapine (50 µM, P < 0.05) and Lithium carbonate (3 mM, P < 0.01). Overall, methylation reduces the promoter activity (ANOVA, main effect of methylation: P = 0.0008). There is no significant reduction in promoter activity on AD treatment in methylated fragment

Effect of single CpG methylation on BDNF exon IV transcriptional regulation

Clinical data from our previous studies highlighted the importance of CpG-87 methylation for the efficacy of ADs. Patients with methylated CpG-87 showed an increase in BDNF plasma levels after 14 days of treatment. To model the clinical effect of CpG-87, we mutated the cytosines (-87, -39 and -35) to methylated cytosines in BDNFIV-pCpGL construct and compared it to control construct. Additionally, we incubated transfected cells with Venlafaxine and Fluoxetine at the given concentrations for 48 h and harvested cells for luciferase assay. We observed a significant overall effect of methylation (F(3, 151) = 49.41; P < 0.0001). No significant effect of ADs was observed (F(4, 151) = 2.23; P = 0.06).

We observed an interaction between CpG methylation and AD treatment (F(12, 151) = 1.86; P = 0.043). This effect was mainly driven by CpG-39 methylation, as the methylated construct showed significant differences in promoter activity in comparison with BDNFIV unmethylated control only upon AD treatment (Venlafaxine 10 µM, P < 0.01; Venlafaxine 50 µM, Fluoxetine 5 µM and 10 µM, P < 0.0001) (Fig. 3).

Fig. 3
figure 3

Effect of single CpG methylation on BDNF exon IV promoter activity. Methylation at CpG position -35, -87 and -39 does not show any significant change in promoter activity in control condition. In the cells treated with ADs, CpG-39 methylation significantly increases the promoter activity (P < 0.0001). Methylation at CpG-35 and -87 had no significant effect

Effect of antidepressants on BDNF exon IV promoter methylation levels

To test the effect of Venlafaxine and Fluoxetine on BDNF exon IV promoter methylation in a drug-naïve environment, we incubated SHSY-5Y cells at two different concentrations (Venlafaxine 10 µM and 50 µM; Fluoxetine 5 µM and 10 µM) of drugs and harvested cells at four different time points (0.5 h, 2 h, 24 h and 48 h).

We used a linear mixed model computing methylation as a dependent variable and CpG position, drugs, concentration and time point as factors. We observed a significant effect of CpG position (F(10, 1301) = 69.482; P < 0.0001). A significant effect of drugs was found on the BDNFIV methylation (F(1, 1301) = 5.44; P = 0.02). Venlafaxine mainly drove this effect (Fig. 4). The drug concentration had a significant effect on the methylation for both Venlafaxine (F(1, 241) = 13.90; P < 0.0001) and Fluoxetine (F(1, 276) = 5.706; P = 0.025). No overall effect of time point was observed, but a significant change in methylation occurred over time in Fluoxetine-treated cells (F(3, 276) = 2.789; P = 0.041) (data not shown).

Fig. 4
figure 4

Antidepressants increase BDNF exon IV promoter methylation. Venlafaxine at 50 µM significantly increases BDNF exon IV promoter methylation in comparison with 10 µM (*P < 0.0001) and 10 µM Fluoxetine in comparison with 5 µM (#P = 0.025)

Effect of DNMT overexpression on BDNF exon IV methylation and transcriptional regulation

To get an insight into methylation and BDNF exon IV transcription, we co-transfected the luciferase construct harboring the BDNF promoter with DNMT3a, DNMT3b and DNMT1 mammalian expression vectors and their empty vector controls pHTC (DNMT1) and pCDNA3 (DNMT3a and 3b). We found that DNMT3a increases the BDNF exon IV promoter activity whereas DNMT3b and DNMT1 have no effect.

Further, we treated the transfected cells with ADs to investigate, if overexpression of methyltransferases has any effect on BDNF transcriptional activity in the presence of ADs, we observed a significant effect of the drugs (F(2, 24) = 9.33; P = 0.001). The DNMT3a-mediated increase on promoter activity was replicable (F(3, 24) = 29.03; P < 0.0001) (Fig. 5a, b).

Fig. 5
figure 5

Effect of DNMT overexpression on BDNF exon IV promoter activity. a DNMT 3a overexpression significantly increases BDNF exon IV promoter activity (P < 0.0001). DNMT3b and DNMT1 do not have any significant effect. b Venlafaxine further increases the DNMT3a-mediated increase in promoter activity (P < 0.01). Fluoxetine has no significant effect

Additionally, we sequenced the cells transfected with DNMTs for methylation analysis to address, if DNMT overexpression has a direct effect on BDNF exon IV promoter CpG methylation. We also used nanaomycin as a chemical modulator as it has been shown to inhibit DNMT3b. We observed no significant effect of DNMT overexpression or nanaomycin treatment on BDNF exon IV promoter methylation.


BDNF exon IV promoter activity is regulated by methylation in a CpG-dependent manner in the absence of neuronal activity [8]. BDNF exon IV promoter methylation has been associated with multiple psychiatric disorders along with MDD [19] as well as treatment response to ADs [24] and other available therapy regimens [25]. Previous work from our lab has confirmed the CpG-87 hypomethylation as a biomarker for unsuccessful AD treatment in two different studies [21, 22]. This work aimed to explore the link between AD action and CpG-87 methylation mechanistically.

To understand the effect of ADs, methylation and their interaction on the transcriptional activity of BDNF exon IV, we cloned BDNF exon IV promoter region in a CpG-free luciferase vector [26] and transfected it into a neuroblastoma cell line SH-SY5Y and treated them with ADs 24 h after transfection for 48 h. The ADs decreased BDNF exon IV promoter activity and this effect was nullified by methylation arguing in favor of a possible link between antidepressant-mediated BDNF increase in serum of patients harboring methylation at CpG-87 locus. These findings corroborated with our previous study where we used the luciferase vector (Promega) containing CpGs demonstrating that the reduction in the promoter activity is due to the CpG methylation of BDNF exon IV promoter and not the vector itself. It must be noted that we methylated the promoter fragment using bacterial methylase M.SssI, thereby methylating every CpG in the promoter to approximately 90% that does not correspond to the physiological methylation pattern.

The linear effect of pharmacological treatments on BDNF exon IV promoter methylation has not been clearly elucidated, mainly because of the lack of pre-treatment methylation data from the patients and lack of drug-naïve cohorts [19]. The clinical data from our previous study show that antidepressants do not change the methylation over 6 weeks of treatment [21]. This patient cohort did not contain any drug-naïve patients, and therefore we could have missed potential AD-mediated changes in methylation. To address this in vitro, we incubated SH-SY5Y cells with ADs and analyzed the CpG methylation of BDNF exon IV promoter. Both Venlafaxine and Fluoxetine increased the mean promoter methylation in a concentration-dependent manner. We did not observe any change in the methylation of CpG-87 upon treatment, suggesting that the treatment itself does not directly modulate CpG-87. CpG-147, CpG-66, CpG-39 and CpG-35 (Additional file 1: Fig. S1) drove the differences in the mean methylation. Our data justify the scope for further studies investigating the molecular mechanisms underlying the antidepressant moderated epigenetic changes.

Further, to differentiate the regulatory effect of methylation of individual CpGs from the whole promoter, we mutated the CpGs of interest to methylated CpGs in the BDNFIV-pCpGL basic construct and performed luciferase assays with or without ADs. We included CpG-87 to model the clinical findings and CpG-35 and -39 due to their proximity with CRE (CREB binding element). CREB has been shown to regulate BDNF protein [27,28,29] and is involved in antidepressant response [30]. We observed no difference in the BDNF exon IV promoter activity upon methylating CpG-35, and -87 in comparison with the unmethylated controls, but CpG-39 methylation increased the transcriptional activity of the promoter. Further plasmid harboring methylated CpG-39 showed a significant increase in promoter activity in the Venlafaxine- and Fluoxetine-treated cells, suggesting yet again possible drug–methylation interaction.

Methyl-CpG-binding protein (MeCP2), an important epigenetic regulator, [31] has been implicated in depression and other psychiatric disorders [32]. MeCP2 can bind to both methylated cytosine (5mC) and 5-hydroxymethylcytosine (5hmC), thereby acting as both transcriptional repressor and activator [33]. Studies have shown that MeCP2 binds to methylated CpGs on BDNF exon IV promoter and represses transcription by blocking CRE [8] or recruiting CREB to a repressor complex [34]. Contradictorily, a strong line of evidence suggests MeCP2 physically and functionally interacts with CREB facilitating the interaction with CRE on the BDNF and other active promoters [35]. Therefore, we speculate that the observed increase in transcriptional activity can be due to enhanced MECP2-CREB complex recruitment to the promoter (Fig. 6).

Fig. 6
figure 6

Proposed molecular mechanism for methylation-specific BDNF exon IV transcriptional activation. CREB, a known transcriptional factor binds to cAMP response element (CRE) on BDNF exon IV promoter and upregulate transcription by recruiting co-activators like CREB-binding protein (CBP). MeCP2 binds across genome, preferably to methylated CpGs and functions as chromatin modifier. Site-specific methylation of BDNF exon IV promoter increases MeCP2 occupancy at the methylated CpGs. MeCP2 and CREB form a transcriptional activating complex. No increase in CREB–CRE interaction was observed, if methylated CpG is far from CRE as with CpG-87 (a) or at the 3′ end of CRE as observed with CpG-35 (b). Methylation at CpG-39 allows the recruitment of MeCP2-CREB protein complex and increases the transcriptional activity at BDNF exon IV promoter (c). Illustration was created with

Antidepressants have been shown to phosphorylate MeCP2 at S421 [36], but we did not observe MeCP2 phosphorylation in SH-SY5Y cells upon antidepressant treatment after 48 h (data not shown). Furthermore, no clear evidence exists for dissociation of MeCP2 from DNA upon calcium influx-induced neuronal depolarization and phosphorylation at specific serine 421 [37], which has been speculated to be a potential mechanism of antidepressant action [38]. Compelling evidence arguing against an important role of S421 phosphorylation for the DNA binding capacity of MeCP2 has been reported by Cohen and co-workers. In S421A mice lacking a functional serine phosphorylation site, no significant change in MeCP2-DNA binding profile was observed in cortical neurons [39].

Therefore, we speculate that the antidepressant-mediated increase in BDNF exon IV promoter activity in the methylated CpG-39 promoter construct is due to an increase in promoter methylation at the CpGs that synchronously upregulate the transcription.

We also conclude that CpG-87 methylation does not directly increase BDNF exon IV transcription. A unique population of cells might contribute to the association of positive treatment response to the methylation (1–5%) at this position (measured in blood), which can be identified by emerging techniques like single-cell epigenomics [40]. Due to technical limitation, we could only methylate single CpG at a time, so we did not account for permutation and combination of methylation patterns that probably exist in the patients. We will be sequencing a larger cohort of MDD patients using Oxford Nanopore Technologies that allows single-strand methylation analysis that can dissect complex epi-allelic methylation patterns.

DNA methyltransferases specifically methylate regulatory CpGs and modulate gene transcription. A suggestive link exists between DNMTs expression and MDD in the literature [41]. To address whether the AD-dependent BDNF increase in serum/plasma levels of patients and animals is associated with in vitro methylation changes, we overexpressed DNMTs and treated the cells with ADs. Of the methyltransferases, only DNMT3a increased exon IV activity in control and Venlafaxine-treated cells. Interestingly, of the two de novo methyltransferases, DNMT3a is expressed in adult brain, whereas DNMT3b is only found in neural progenitor cells [17]. When we bisulfite-sequenced the DNMT overexpressed cells, we did not observe any change in the CpG methylation at BDNF exon IV promoter. Taken together, our findings suggest an indirect effect of DNMT3a overexpression on an increase in BDNF exon IV promoter activity. It will be very interesting to get a deeper insight into the expression pattern of DNMTs in the responders and non-responders.

Our system relies on artificially curated methylation marks that do not reflect the complex scenario in vivo. Therefore, it will be interesting to use Cas9-mediated methylation editing by recruiting methyl modifiers (DNMT/TET) to study the interaction between BDNF exon IV promoter methylation and antidepressant action on its regulation. Another limitation of the study is that we used a neuroblastoma cell line, thus lacking the communication between different cell types that is essential for antidepressant action. Therefore, we plan to replicate the findings in a mixed cell culture model consisting neurons, astrocytes and myelinating oligodendrocytes.


Our study provides evidence for AD-dependent methylation changes on the BDNF exon IV promoter and emphasizes the importance of studying such interactions in the drug-naïve cohorts of MDD patients. It also describes the functional relevance of studying CpG-specific methylation and its cross talk with AD-mediated increase in BDNF exon IV promoter activity. Further, we highlight the importance of studying DNMTs, especially DNMT3a in the context of MDD.


Cell culture

Human neuroblastoma cell line SH-SY5Y (ATCC, Wesel) was cultured in DMEM-Ham’s F12 medium 1:1 (Biochrom AG, Berlin) supplemented with 10% FBS Superior (Biochrom AG, Berlin) and 1% penicillin/streptomycin (PAN Biotech GmbH, Aidenbach). Cells were maintained at 37 °C in a humidified atmosphere of 5% CO2 in atmosphere. Cells were incubated with indicated concentration and for the indicated time with the drugs.

Insert amplification and cloning

BDNF exon IV promoter region (GRCh38; chr11: 27701325-27702402) was amplified (primers: Inf-BDNF-pCpGL Fw/Rv) from the genomic DNA using KOD Hot Start Taq DNA Polymerase (Novagen, Merck Chemicals, Darmstadt, Germany). Primers contained pCpGL CpG-free luciferase vector-specific overhangs (red) that were adapted to the PCR product. The 1172 bp fragment with adapter sequences was cloned into HindIII HF-linearized pCpGL firefly luciferase reporter plasmid [26] by recombination using In-Fusion® HD Cloning Kit (Clontech, cat. no. 639648) (Additional file 2: Fig. S2.). Similarly, DNA methyltransferase 1 (DNMT1) cDNA were PCR amplified (primers: Inf-DNMT1-halo Fw and Rv) from a mixture of human cell line cDNA using KOD Hot Start Taq DNA Polymerase and cloned in pHTC HaloTag®CMV-neo Vector (Promega, cat. no. G7711). The plasmids pcDNA3/Myc-DNMT3a and pcDNA3/Myc-DNMT3b were gifts from Arthur Riggs (Addgene, cat. no. 35521 and cat. no. 35522, respectively [42]. Plasmids were transformed in NEB alpha competent cells except pCpGLempty and pCpGL-BDNFIV that were transformed in piR1 Competent E. coli (Invitrogen) with 5 µL recombined plasmid and plated on LB-Agar containing 100 µg/mL Ampicillin: Grown colonies were transferred to 3 mL LB-Amp-Medium and plasmids were isolated using the NucleoSpin Plasmid Kit (Macherey–Nagel, Düren). After verifying the plasmids containing the hBDNF exon IV promoter region by sequencing (primer: pCpGL-seq_Fw), DNA purification was done using endotoxine-free plasmid isolation kit (NucleoBond Xtra Midi EF Kit; Macherey–Nagel, Düren). Expression of DNMT1, DNMT3a and DNMT3b was validated by western blotting (data not shown). Refer to Table 1 for primer sequences. The underlined sequence is gene-specific.

Table 1 Primer list

In vitro methylation of BDNF exon IV promoter construct

The pCpGL-BDNF IV promoter reporter (4 µg) was incubated (4 h at 37 °C) with CpG Methyltransferase M.SssI (20 U; New England Biolabs, Ipswich, MA) in buffer containing 640 µM S-adenosylmethionine (New England Biolabs, Ipswich, MA). Subsequently, the methylated plasmid was purified with NucleoSpin Plasmid Kit (Macherey–Nagel, Düren).

Site-directed mutagenesis

CpGs-87, -39 and -35 were methylated using the Q5 Site-Directed Mutagenesis Kit (NEB, E0554) following the manufacturer’s protocol. Briefly, the pCpGL-BDNFIV promoter reporter construct was reverse-amplified using mut -87, mut -39 and mut -35 or mut-cntrl Fw/Rv primers and Q5 Hot Start High-Fidelity DNA Polymerase. The primers were back-to-back primers and were designed using NEBaseChanger. Cytosine at the desired position was replaced with methylated cytosine. Multiple reactions were pooled to gain a sufficient starting DNA amount. The PCR product was combined with the KLD Mix, containing a blend of kinase, ligase and DpnI enzymes, and was incubated for 15 min. After KLD treatment, the reaction was purified using NucleoSpin Plasmid Kit (Macherey–Nagel) and was used directly for transfection. The efficiency of mutagenesis was assessed by Sanger sequencing. All three CpGs were about 60% methylated. Refer to Table 1 for primer sequences.

Transfection and luciferase-based reporter gene assays

1.5 × 105 SH-SY5Y cells were plated in 96 well plates for 24 h before transfection. 100 ng/well of pCpGL-BDNFIV promoter luciferase reporter plasmid or methylated pCpGL-BDNFIV promoter luciferase reporter plasmid were co-transfected with 2 ng/well pGL4.74 [hRluc/TK] renilla luciferase control plasmid (Promega, Madison, WI) for normalization (3 columns per reporter plasmid) using Lipofectamine 3000 (company). 50 ng of the luciferase constructs was used for transfection, both for control and for mutated CpGs. Cells were harvested 48 h after transfection. For DNMT overexpression experiments, 25 ng of DNA was used for transfection for both DNMTs and empty vectors. Antidepressants were added 12 h post-transfection. After incubation over 48 h, cells were lysed and luminescence was measured using the Dual-Luciferase Reporter Assay System (Promega, Madison) and the GloMax-Multi + Luminometer (Promega, Madison), according to the manufacturer’s protocol.

Luciferase assays

Post-transfection, media was removed and cells were washed with PBS. Assay was performed using Dual-Luciferase Reporter Assay System (Promega), and the luminescence was recorded using a double injector program on GloMax-Multi Detection System. The assay was applied in triplicates and was repeated 3 times. For analysis, firefly readings were normalized by the respective renilla luminescence readings and normalized values were imported to GraphPad Prism 5 (Graph Pad Inc., San Diego, CA).

Bisulfite sequencing of the BDNF exon IV promoter region

Genomic DNA was extracted from cells using NucleoMag Blood 200µL (Macherey–Nagel) according to the manufacturer’s protocol. Afterward, 500 ng of genomic DNA was modified by sodium bisulfite using the EpiTect® 96 Bisulfite Kit (QIAGEN).

Primers were designed to amplify a region covering a fragment of 277 base pairs (bp) from − 200 to + 77 bp relative to the starting bp of the BDNF exon IV promoter (NCBI NC_000011.9: 27723103-27723380) containing 13 CpG sites within the promoter region of BDNF exon IV. To amplify a specific product, a semi-nested polymerase chain reaction (PCR) was performed using 1 µL bisulfite modified DNA or PCR round 1 product with 0.4 µL of each primer and 5 µL HotStarTaq Master Mix (Qiagen) filled to a reaction volume of 10 µL.

Sequencing of the reverse strand was performed using a BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) on an Applied Biosystems® 3500xL Genetic Analyzer (Applied Biosystems; POP-7™ polymer) according to the manufacturer’s instructions. Primer sequences are given in Table 1 (Amplification: bisBDNF_IV_forw1/2 and bisBDNF_IV_rev; Sequencing: BDNF_IV_seq). Refer to Table 1 for primer sequences.

Statistical analysis

The methylation of each of the 11 CpGs in the CpG island was assessed using the Epigenetic Sequencing Methylation analysis Software (ESME) [43] which compares each methylation site to the original sequence of the promoter. To analyze the effect of antidepressant treatment on BDNF exon IV promoter methylation, mixed linear models for repeated measurements were performed including the factors: CpG position, drugs (Venlafaxine/ Fluoxetine), time point (0.5 h, 2 h, 24 h and 48 h) and drug concentration (5 µM and 10 µM for Fluoxetine and 10 µM and 50 µM for Venlafaxine). For the model assessing the effect of DNMTs and nanaomycin, CpG position and methyl modulators were used as factors. Data were analyzed employing IBM SPSS Statistics for Windows, Version 21.0. (Armonk, NY: IBM Corp.) and illustrated using Graph Pad Prism 5 (Graph Pad Inc., San Diego, CA). All the luciferase data were analyzed using two-way ANOVA in GraphPad Prism 5.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.





Brain-derived neurotrophic factor


CAMP response element-binding protein


CREB-binding protein


DNA methyltransferase


Major depressive disorder


Methyl-CpG-binding protein 2


Serotonin reuptake inhibitors


Serotonin–norepinephrine reuptake inhibitors


  1. Kraus C, Kadriu B, Lanzenberger R, Zarate CA Jr, Kasper S. Prognosis and improved outcomes in major depression: a review. Transl Psychiatry. 2019;9(1):127.

    Article  Google Scholar 

  2. Pandarakalam JP. Challenges of treatment-resistant depression. Psychiatr Danub. 2018;30(3):273–84.

    Article  CAS  Google Scholar 

  3. Fuchikami M, Morinobu S, Segawa M, Okamoto Y, Yamawaki S, Ozaki N, Inoue T, Kusumi I, Koyama T, Tsuchiyama K, 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.

    Article  CAS  Google Scholar 

  4. Bjorkholm C, Monteggia LM. BDNF: a key transducer of antidepressant effects. Neuropharmacology. 2016;102:72–9.

    Article  Google Scholar 

  5. Altar CA, Whitehead RE, Chen R, Wortwein G, Madsen TM. Effects of electroconvulsive seizures and antidepressant drugs on brain-derived neurotrophic factor protein in rat brain. Biol Psychiat. 2003;54(7):703–9.

    Article  CAS  Google Scholar 

  6. Pruunsild P, Kazantseva A, Aid T, Palm K, Timmusk T. Dissecting the human BDNF locus: bidirectional transcription, complex splicing, and multiple promoters. Genomics. 2007;90(3):397–406.

    Article  CAS  Google Scholar 

  7. Boulle F, van den Hove DL, Jakob SB, Rutten BP, Hamon M, van Os J, Lesch KP, Lanfumey L, Steinbusch HW, Kenis G. Epigenetic regulation of the BDNF gene: implications for psychiatric disorders. Mol Psychiatry. 2012;17(6):584–96.

    Article  CAS  Google Scholar 

  8. Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, Fan G, Sun YE. DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science. 2003;302(5646):890–3.

    Article  CAS  Google Scholar 

  9. Hing B, Sathyaputri L, Potash JB. A comprehensive review of genetic and epigenetic mechanisms that regulate BDNF expression and function with relevance to major depressive disorder. Am J Med Genet B Neuropsychiatr Genet. 2018;177(2):143–67.

    Article  CAS  Google Scholar 

  10. Finkbeiner S, Tavazoie SF, Maloratsky A, Jacobs KM, Harris KM, Greenberg ME. CREB: a major mediator of neuronal neurotrophin responses. Neuron. 1997;19(5):1031–47.

    Article  CAS  Google Scholar 

  11. Gupta S, Kim SY, Artis S, Molfese DL, Schumacher A, Sweatt JD, Paylor RE, Lubin FD. Histone methylation regulates memory formation. J Neurosci Off J Soc Neurosci. 2010;30(10):3589–99.

    Article  CAS  Google Scholar 

  12. Collins BE, Greer CB, Coleman BC, Sweatt JD. Histone H3 lysine K4 methylation and its role in learning and memory. Epigenetics Chromatin. 2019;12(1):7.

    Article  Google Scholar 

  13. Mitchelmore C, Gede L. Brain derived neurotrophic factor: epigenetic regulation in psychiatric disorders. Brain Res. 2014;1586:162–72.

    Article  CAS  Google Scholar 

  14. Januar V, Saffery R, Ryan J. Epigenetics and depressive disorders: a review of current progress and future directions. Int J Epidemiol. 2015;44(4):1364–87.

    Article  Google Scholar 

  15. Januar V, Ancelin ML, Ritchie K, Saffery R, Ryan J. BDNF promoter methylation and genetic variation in late-life depression. Transl Psychiatry. 2015;5:e619.

    Article  CAS  Google Scholar 

  16. Vialou V, Feng J, Robison AJ, Nestler EJ. Epigenetic mechanisms of depression and antidepressant action. Annu Rev Pharmacol Toxicol. 2013;53:59–87.

    Article  CAS  Google Scholar 

  17. Yu NK, Baek SH, Kaang BK. DNA methylation-mediated control of learning and memory. Mol Brain. 2011;4:5.

    Article  CAS  Google Scholar 

  18. Wolkowitz OM, Wolf J, Shelly W, Rosser R, Burke HM, Lerner GK, Reus VI, Nelson JC, Epel ES, Mellon SH. Serum BDNF levels before treatment predict SSRI response in depression. Prog Neuropsychopharmacol Biol Psychiatry. 2011;35(7):1623–30.

    Article  CAS  Google Scholar 

  19. Zheleznyakova GY, Cao H, Schioth HB. BDNF DNA methylation changes as a biomarker of psychiatric disorders: literature review and open access database analysis. Behavioral Brain Funct BBF. 2016;12(1):17.

    Article  Google Scholar 

  20. Ikegame T, Bundo M, Murata Y, Kasai K, Kato T, Iwamoto K. DNA methylation of the BDNF gene and its relevance to psychiatric disorders. J Hum Genet. 2013;58(7):434–8.

    Article  CAS  Google Scholar 

  21. Tadic A, Muller-Engling L, Schlicht KF, Kotsiari A, Dreimuller N, Kleimann A, Bleich S, Lieb K, Frieling H. Methylation of the promoter of brain-derived neurotrophic factor exon IV and antidepressant response in major depression. Mol Psychiatry. 2014;19(3):281–3.

    Article  CAS  Google Scholar 

  22. Lieb K, Dreimuller N, Wagner S, Schlicht K, Falter T, Neyazi A, Muller-Engling L, Bleich S, Tadic A, Frieling H. BDNF plasma levels and BDNF Exon IV promoter methylation as predictors for antidepressant treatment response. Front Psych. 2018;9:511.

    Article  Google Scholar 

  23. Wang P, Zhang C, Lv Q, Bao C, Sun H, Ma G, Fang Y, Yi Z, Cai W. Association of DNA methylation in BDNF with escitalopram treatment response in depressed Chinese Han patients. Eur J Clin Pharmacol. 2018;74(8):1011–20.

    Article  CAS  Google Scholar 

  24. Webb LM, Phillips KE, Ho MC, Veldic M, Blacker CJ. The Relationship between DNA methylation and antidepressant medications: a systematic review. Int J Mol Sci. 2020;21(3):826.

    Article  CAS  Google Scholar 

  25. Kleimann A, Kotsiari A, Sperling W, Groschl M, Heberlein A, Kahl KG, Hillemacher T, Bleich S, Kornhuber J, Frieling H. BDNF serum levels and promoter methylation of BDNF exon I, IV and VI in depressed patients receiving electroconvulsive therapy. J Neural Transm. 2015;122(6):925–8.

    Article  CAS  Google Scholar 

  26. Klug M, Rehli M. Functional analysis of promoter CpG methylation using a CpG-free luciferase reporter vector. Epigenetics. 2006;1(3):127–30.

    Article  Google Scholar 

  27. Shieh PB, Hu SC, Bobb K, Timmusk T, Ghosh A. Identification of a signaling pathway involved in calcium regulation of BDNF expression. Neuron. 1998;20(4):727–40.

    Article  CAS  Google Scholar 

  28. West AE, Chen WG, Dalva MB, Dolmetsch RE, Kornhauser JM, Shaywitz AJ, Takasu MA, Tao X, Greenberg ME. Calcium regulation of neuronal gene expression. Proc Natl Acad Sci USA. 2001;98(20):11024–31.

    Article  CAS  Google Scholar 

  29. Tao X, Finkbeiner S, Arnold DB, Shaywitz AJ, Greenberg ME. Ca2+ influx regulates BDNF transcription by a CREB family transcription factor-dependent mechanism. Neuron. 1998;20(4):709–26.

    Article  CAS  Google Scholar 

  30. Blendy JA. The role of CREB in depression and antidepressant treatment. Biol Psychiat. 2006;59(12):1144–50.

    Article  CAS  Google Scholar 

  31. Hansen JC, Ghosh RP, Woodcock CL. Binding of the Rett syndrome protein, MeCP2, to methylated and unmethylated DNA and chromatin. IUBMB Life. 2010;62(10):732–8.

    Article  CAS  Google Scholar 

  32. Ausio J. MeCP2 and the enigmatic organization of brain chromatin: implications for depression and cocaine addiction. Clin Epigenetics. 2016;8:58.

    Article  Google Scholar 

  33. Li Y, Wang H, Muffat J, Cheng AW, Orlando DA, Loven J, Kwok SM, Feldman DA, Bateup HS, Gao Q, et al. Global transcriptional and translational repression in human-embryonic-stem-cell-derived Rett syndrome neurons. Cell Stem Cell. 2013;13(4):446–58.

    Article  Google Scholar 

  34. Tai DJ, Liu YC, Hsu WL, Ma YL, Cheng SJ, Liu SY, Lee EH. MeCP2 SUMOylation rescues Mecp2-mutant-induced behavioural deficits in a mouse model of Rett syndrome. Nat Commun. 2016;7:10552.

    Article  CAS  Google Scholar 

  35. Chahrour M, Jung SY, Shaw C, Zhou X, Wong ST, Qin J, Zoghbi HY. MeCP2, a key contributor to neurological disease, activates and represses transcription. Science. 2008;320(5880):1224–9.

    Article  CAS  Google Scholar 

  36. Jin HJ, Pei L, Li YN, Zheng H, Yang S, Wan Y, Mao L, Xia YP, He QW, Li M, et al. Alleviative effects of fluoxetine on depressive-like behaviors by epigenetic regulation of BDNF gene transcription in mouse model of post-stroke depression. Sci Rep. 2017;7(1):14926.

    Article  Google Scholar 

  37. Zhou Z, Hong EJ, Cohen S, Zhao WN, Ho HY, Schmidt L, Chen WG, Lin Y, Savner E, Griffith EC, et al. Brain-specific phosphorylation of MeCP2 regulates activity-dependent Bdnf transcription, dendritic growth, and spine maturation. Neuron. 2006;52(2):255–69.

    Article  CAS  Google Scholar 

  38. Hutchinson AN, Deng JV, Cohen S, West AE. Phosphorylation of MeCP2 at Ser421 contributes to chronic antidepressant action. J Neurosci Off J Soc Neurosci. 2012;32(41):14355–63.

    Article  CAS  Google Scholar 

  39. Cohen S, Gabel HW, Hemberg M, Hutchinson AN, Sadacca LA, Ebert DH, Harmin DA, Greenberg RS, Verdine VK, Zhou Z, et al. Genome-wide activity-dependent MeCP2 phosphorylation regulates nervous system development and function. Neuron. 2011;72(1):72–85.

    Article  CAS  Google Scholar 

  40. Schwartzman O, Tanay A. Single-cell epigenomics: techniques and emerging applications. Nat Rev Genet. 2015;16(12):716–26.

    Article  CAS  Google Scholar 

  41. Duan Z, Lu J. DNA methyltransferases in depression: an update. Front Psych. 2020;11:538683.

    Article  Google Scholar 

  42. Chen ZX, Mann JR, Hsieh CL, Riggs AD, Chedin F. Physical and functional interactions between the human DNMT3L protein and members of the de novo methyltransferase family. J Cell Biochem. 2005;95(5):902–17.

    Article  CAS  Google Scholar 

  43. Lewin J, Schmitt AO, Adorjan P, Hildmann T, Piepenbrock C. Quantitative DNA methylation analysis based on four-dye trace data from direct sequencing of PCR amplificates. Bioinformatics. 2004;20(17):3005–12.

    Article  CAS  Google Scholar 

Download references


Open Access funding enabled and organized by Projekt DEAL. In-house funds.

Author information

Authors and Affiliations



HF and HP contributed to conceptualization and formal analysis; HP, ABo and SG investigated the study; HP contributed to writing—original draft and preparation, and visualization; HP, ABu and HF contributed to writing—review and editing; HF supervised the study. All authors have read and agreed to the submitted version of the manuscript.

Corresponding author

Correspondence to Hansi Pathak.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

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

Additional file1

. Figure S1: Effect of antidepressants on single CpG methylation after 48 h of incubation with Venlafaxine (50 µM, 10 µM) and Fluoxetine (5 µM, 10 µM).

Additional file2

. Figure S2 Plasmid map of pCpGL-BDNFIV construct. BDNF exon IV promoter insert (orange) is cloned upstream of luciferase gene (neon green). Beginning of exon IV is indicated (grey). The PCR-amplified fragment used for methylation analysis (yellow) and the CpGs (pink) are marked. The transcription factor binding sites CARE1 and CRE are highlighted (purple). The plasmid map was made using Benchling [Biology Software]. (2022). Retrieved from

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pathak, H., Borchert, A., Garaali, S. et al. BDNF exon IV promoter methylation and antidepressant action: a complex interplay. Clin Epigenet 14, 187 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: