Table of Contents

Interaction of BDNF and COMT Polymorphisms on Paired-Associative Stimulation-Induced Cortical Plasticity

  1. Agnes Flöel1,2,3
  1. 1Department of Neurology,
  2. 2NeuroCure Cluster of Excellence, and
  3. 3Center for Stroke Research Berlin, Charité-Universitätsmedizin Berlin, 10117 Berlin, Germany,
  4. 4Department of Neurology, and
  5. 5Internal Medicine D, Department of Nephrology, Hypertension, and Rheumatology, University of Münster, 48149 Münster, Germany, and
  6. 6Department of Psychiatry and Psychotherapy, Philipps-University Marburg, 35039 Marburg, Germany
  1. Author contributions: A.F. designed research; J.K., S.J., A.S., and J.S. performed research; A.V.W., J.K., A.S., J.S., and A.F. analyzed data; A.V.W., J.K., and A.F. wrote the paper.

Abstract

The common single-nucleotide polymorphism (SNP) brain-derived neurotrophic factor (BDNF) valine-to-methionine substitution at codon 66 (Val66Met) has been associated with differences in memory functions and cortical plasticity following brain stimulation. Other studies could not confirm these results, though, and potential interactions of BDNF carrier status with other learning-relevant SNPs are largely unknown. The present study aimed to evaluate the effects of BDNF Val66Met genotype on paired associative stimulation (PAS)-induced motor cortex plasticity, while additionally taking catechol-O-methyltransferase (COMT) Val158Met and kidney and brain (KIBRA) rs17070145 carrier status into account. Therefore, a cohort of 2 × 16 age- and education-matched healthy young females underwent transcranial magnetic stimulation using an excitatory PAS25 protocol to induce cortical plasticity. Cognitive performance was assessed using implicit grammar- and motor-learning tasks and a detailed neuropsychological test battery. While BDNF carrier status alone did not significantly influence PAS-induced cortical plasticity, we found a significant BDNF × COMT interaction, showing higher plasticity immediately following the PAS25 protocol for the BDNF Val/Val vs Met genotype in COMT Met homozygotes only (ANOVA, p = 0.027). A similar advantage for this group was noted for implicit grammar learning (ANOVA, p = 0.021). Accounting for KIBRA rs17070145 did not explain significant variance. Our findings for the first time demonstrate an interaction of BDNF by COMT on human cortical plasticity. Moreover, they show that genotype-related differences in neurophysiology translate into behavioral differences. These findings might contribute to a better understanding of the mechanisms of interindividual differences in cognition.

Introduction

Common polymorphisms in genes that encode neuromodulatory proteins, including brain-derived neurotrophic factor (BDNF), catechol-O-methyltransferase (COMT), and kidney and brain (KIBRA), may contribute to interindividual differences in cognitive performance (Egan et al., 2001, 2003; Chen et al., 2006; Papassotiropoulos et al., 2006).

One of the most extensively studied polymorphisms, a single-nucleotide polymorphism (SNP) at nucleotide 196 (G/A) in the BDNF gene has been associated with decreased activity-dependent secretion of BDNF, due to a valine-to-methionine substitution at codon 66 (Val66Met) (Egan et al., 2003). Thus, Met allele carriers are thought to exhibit lower BDNF levels than Val/Val homozygotes, and animal models provided good evidence that lower levels of the neurotrophin BDNF are linked with impairments in hippocampal long-term potentiation (LTP), synaptic integrity, and hippocampus-dependent learning (Patterson et al., 1996; Mizuno et al., 2000; Poo, 2001).

In humans, cortical plasticity has been assessed for BDNF Val66Met in several studies. Using noninvasive brain stimulation protocols, reduced plasticity in BDNF Met allele carriers was seen for most, but not all, protocols and techniques (Kleim et al., 2006) (for review, see Cheeran et al., 2009). Cheeran et al. (2008) reported a trend for lower cortical plasticity, as determined using peripheral electric stimulation and subsequent transcranial magnetic stimulation (TMS) (Stefan et al., 2000), in BDNF Met carriers compared with Val/Val homozygotes. Likewise, others reported higher plasticity-induction in Val/Val homozygotes using short-term training and cortical mapping (Kleim et al., 2006; McHughen et al., 2010) or intermittent theta-burst stimulation (iTBS) (Antal et al., 2010). On the other hand, Antal et al. (2010) found larger transcranial direct current stimulation (tDCS)-induced modulation of motor cortex excitability in Met carriers. Others failed to detect significant differences between BDNF Val66Met genotype groups for repetitive TMS- or iTBS-induced plasticity (Li Voti et al., 2011).

The above-described results for BDNF Val66Met may in part be explained by study-specific differences in the stimulation protocols or study populations. However, potential gene–gene interactions, e.g., between BDNF Val66Met and COMT Val158Met (Erickson et al., 2008; Nagel et al., 2008), had not been considered (Kleim et al., 2006; Cheeran et al., 2008). Notably, data from studies on cognition suggest that a significant interactive effect may exist that may turn a nonsignificant association, for example between BDNF and executive function, into a significant association if the interaction between BDNF and COMT Val158Met is taken into account (Nagel et al., 2008). The COMT Met allele has been linked with less enzymatic activity of COMT and thus higher prefrontal dopamine (Lindenberger et al., 2008). Several studies reported better executive function and memory performance in COMT Met carriers; however, others could not confirm these results (for review, see Witte and Flöel, 2011). For an SNP at intron 9 in the KIBRA gene (rs17070145), both advantages (Papassotiropoulos et al., 2006) and disadvantages have been reported for a C-to-T substitution (for review, see Wersching et al., 2011). The KIBRA protein is discussed to be involved in brain development and memory formation as a postsynaptic scaffold protein (Johannsen et al., 2008; Yoshihama et al., 2009). No study so far assessed the impact of KIBRA rs17070145 or COMT Val158Met on cortical plasticity using neurophysiological measurements.

Together, the impact of BDNF Val66Met on neuronal plasticity remains controversial, and interactions of BDNF carrier status with other learning-relevant SNPs are not yet understood. Therefore, the present study aimed to determine the effects of Val66Met on LTP-like synaptic plasticity, while additionally accounting for COMT Val158Met and KIBRA rs17070145 carrier status. Moreover, a comprehensive assessment of implicit learning ability and neuropsychological test performance was taken. To avoid confounding by gender (Wersching et al., 2011), a female-only cohort was studied comprising 16 healthy, young Val homozygotes and 16 Met carriers who were matched for age and education.

Materials and Methods

Subjects

Seventy-six healthy young volunteers were initially recruited via advertisements in local newspapers and university websites. Peripheral blood samples were collected for genetic analyses if meeting the following inclusion criteria: (1) no history of neurological or psychiatric disorders, (2) no drug abuse or CNS medication, (3) right-handedness (>70% on the Edinburgh Handedness Inventory) (Oldfield, 1971), (4) native German speaker, and (5) normal or corrected-to-normal vision. To avoid confounding by gender (Wersching et al., 2011), only female participants were chosen. Of the initial sample, 16 female BDNF Met allele carriers (15 Val/Met, 1 Met/Met; age, 26.9 ± 8.5 years, mean ± SD) were identified and pairwise matched for age and years of education to 16 homozygous BDNF Val/Val carriers. These 32 women were invited to participate in three experimental sessions, assessing cortical plasticity as induced by the excitatory paired associative stimulation (PAS25) protocol, implicit learning in the motor domain [serial reaction time task (SRTT)] and in the verbal domain [artificial grammar learning (AGL)], as well as neuropsychological performance (Fig. 1, experimental outline). All participants were of Caucasian ethnicity (for demographic details, see Table 1).

Figure 1.

Study flow. At the screening session, subjects underwent medical examination and donated blood samples for genotyping analyses. After identification of 16 female BDNF Met carriers and 16 age- and education-matched BDNF Val/Val carriers, these subjects were tested on cognitive performance using standard neuropsychological tests (day 1; for details, see Materials and Methods). In a randomized order on day 2 or 3, subjects underwent assessment of plasticity of the motor cortex, using a PAS25 protocol (adapted from Stefan et al., 2000). On the other day, implicit motor and grammar learning were assessed using a SRTT and an AGL task.

Table 1.

Demographics of BDNF genotype groups

The study was approved by the local ethics committee of the University of Münster and was conducted in accordance with the Declaration of Helsinki. All participants gave written informed consent.

Genotyping

DNA was extracted from whole blood using a blood mini-kit (Quiagen). We used GenBank sequences (http://genome.ucsc.edu) to identify the common coding variant in the BDNF gene, a G → A polymorphism responsible for the Val66Met change (Sen et al., 2003). Based on a primer database (http://primer3.sourceforge.net/), we designed sequencing primers (BDNF; forward, AGAGTGATGACCATCCTTTTCC; reverse, GCCCGAACTTTCTGGTCCTC) using a Big Dye Terminator Cycle Sequencing Ready Reactio Kit (Applied Biosystems), resulting in a 241 bp PCR product. Amplification reactions were performed using Taq-Polymerase (Quiagen), followed by gel electrophoresis of the PCR products. Sequencing was done on a 3730 DNA-Analyzer (Applied Biosystems). For further details, see Sen et al. (2003).

To assess potential interaction effects between BDNF genotype and other learning-relevant polymorphisms, participants were additionally screened for COMT Val158Met and KIBRA rs17070145 genotype status using standard procedures (for details, see Witte et al., 2010; Wersching et al., 2011). Genotyping revealed 25 COMT Met allele carriers (10 Met/Met, 15 Val/Met), 7 COMT Val/Val carriers, 14 KIBRA rs17070145 T-allele carriers (T → C allele substitution at intron 9; 4 T/T, 10 C/T), and 18 C/C carriers (for group stratification based on BDNF genotype, see Table 2). Distributions of BDNF and COMT genotypes were in Hardy–Weinberg equilibrium (χ2 < 1.65, p > 0.05).

Table 2.

Distribution of BNDF × COMT and BDNF × KIBRA genotype

Neurophysiological measurements

Determination of baseline stimulation parameters.

Single-pulse TMS was performed using a Magstim 200 connected to a figure-eight coil with an outer diameter of 9.5 cm. TMS was delivered to the motor cortex of the left hemisphere. The optimal coil position (“hot spot”) was defined as the location over the motor cortex producing the highest and most stable motor evoked potential (MEP) amplitudes of the abductor digiti minimi (ADM) muscle. During stimulation the coil was held tangentially to the participant's skull at an angle of 45° to the sagittal plane and with the handle pointing backward (Fig. 2). Resting motor threshold (RMT) was defined as the lowest stimulation intensity that could evoke MEPs of 50 μV in 5 of 10 consecutive pulses. The intensity of stimulation was given as the percentage of maximal stimulator output. MEPs of the ADM were recorded via surface EMG activity using disposable Ag-AgCl electrodes in a belly tendon montage. MEP signals were amplified and bandpass filtered between 100 Hz and 10 kHz (Portable Keypoint).

Figure 2.

Cortical plasticity assessed by TMS using the excitatory PAS25 protocol. TMS of the hand field of the contralateral motor cortex (M1) was used to evoke MEPs of the ADM muscle of the right hand. During PAS25, TMS was simultaneously administered with peripheral electrical stimulation of the right N. ulnaris for 90 times (flashes) with an interstimulus interval of 25 ms. Single TMS-induced MEPs were collected at baseline (pre) and at post 0, post 15, post 30, and post 60 after PAS25. An increase in MEP amplitudes after PAS25 indicates induction of synaptic plasticity in M1.

PAS-induced cortical plasticity.

All subjects underwent the PAS protocol that was similar to the one initially described by Stefan et al. (2000) (see also List et al., 2011). PAS consisted of 90 pairs of electrical stimulation of the ulnar nerve (300% of individual sensory threshold) at the wrist combined with suprathreshold TMS impulse over the hot spot of the ADM over the contralateral hemisphere at a rate of 0.05 Hz (total duration of stimulation protocol, 30 min). The stimulator output was set to induce MEP values of 0.5–1 mV. Interstimulus interval was set to 25 ms, which has been shown to induce LTP-like increase in MEP amplitude (Stefan et al., 2000). MEP amplitude measurements (average of 20 responses) were performed at baseline (pre) and at 0 min (post 0), 15 min (post 15), 30 min (post 30), and 60 min (post 60) after PAS25 (Fig. 2). As primary outcome measures, MEP amplitudes of each time point were averaged and normalized to the MEP amplitude at baseline for each subject. Additionally, the grand average of adjusted MEPs at post 0, post 15, post 30, and post 60 were calculated, indicating responders (value >1.0) and nonresponders (value ≤1.0) to the PAS25 protocol (Muller-Dahlhaus et al., 2008). To maintain a standardized level of attention during the PAS intervention, subjects were instructed to stay alert, voluntarily relax the ADM of the dominant arm, and count the number of ulnar nerve stimulations. Muscle relaxation was continuously monitored by visual feedback from the surface EMG.

Implicit learning

Artificial grammar learning.

To assess implicit learning in the language domain, we used an AGL task, originally designed by Reber (1967). Therefore, 200 sequences of the letters M, V, S, R, and X were generated following a finite state grammar by Reber (1967). The Reber grammar defines the possible arrangements of the letters dependent on a grammatical rule [e.g., that the letter S can follow letters M (repeatedly) and V (once), but not X and R; thus, MSVS and MSSVS would be correct, but not MSXS, MSSXS, or MSVSS]. One hundred of the correct letter strings, consisting of 4 to 12 letters, were used as grammatical stimuli in an initial acquisition task. Here, subjects were asked to pay close attention and to retype each string that appeared on a computer screen, without feedback. Subjects were unaware of the underlying concept to ensure implicit learning of the grammar. Subsequently, subjects were told that the strings followed a certain grammatical rule, and asked to discriminate between grammatical and nongrammatical items in a following classification task. As stimuli for this classification phase, another 100 grammatical strings plus 50 nongrammatical strings (created by switching two nonterminal letters) were presented in random order on the screen. Subjects were instructed to respond spontaneously whether they thought the presented stimulus was grammatical or nongrammatical.

Numbers of “hits” and “misses” [for grammatical (G) items], and “false alarms” and “correct rejections” [for nongrammatical (NG) items] of the classification task were assessed. Overall hit rate for G items and overall false alarms for NG items created the “endorsement rate” (percentage of items perceived as grammatical), a frequently used variable in the AGL literature to assess syntactic violation detection (Lieberman et al., 2004; Forkstam et al., 2006; de Vries et al., 2010; Kurten et al., 2010). Additionally, d′ values were calculated, a measure of classification accuracy controlling for response bias (Hochhaus, 1972; Flöel et al., 2009).

Serial reaction time task.

To assess implicit learning in the motor domain, we measured performance of an SRTT (Nissen and Bullemer, 1987; Cleeremans and McClelland, 1991; Jiménez et al., 1996), using a modified version with a probabilistic sequence (85% sequential elements plus 15% random elements). The sequential structure was generated by a finite-state grammar that defined the permissible transitions between successive stimuli. Each stimulus, presented on a computer screen, consisted of an asterisk as substitution of one of five dashes in a row, serving as indicator of the finger 1 (thumb), 2, 3, 4, or 5 (little finger). Subjects had their right hand on a custom-made keypad, which provided a key underlying each fingertip, and were asked to respond to each stimulus by pressing the corresponding key as fast as possible. Subjects had to perform 8 blocks of 250 sequential and random stimuli each. Five hundred milliseconds (ms) after each keypress, the subsequent stimulus appeared, and subjects had a resting interval of 3 min between blocks. No feedback was provided during the task. Reaction times (RTs, in ms) were recorded throughout the experiment and averaged for random and sequential stimuli among blocks. While subjects were familiarized with the task during the first four blocks, mean RT differences between random and sequential stimuli of the last four blocks was taken as outcome measure of procedural motor learning, similar to previous studies (Cleeremans and McClelland, 1991; Rösser et al., 2008).

Neuropychological testing

We used a comprehensive neuropsychological test battery to assess memory functions, executive functions, and general intellectual abilities (for an overview, see Table 3).

Table 3.

Neuropsychological testing

Briefly, in the Auditory Verbal Learning Task (German version) (Helmstaedter and Kurthen, 2001), subjects were asked to learn and recall a list of 15 words, with delayed retrieval and recognition trials after 30 min (delayed episodic memory). In subtests of the Wechsler Memory Scale-Revised (Elwood, 1991), subjects had to recall and modify small sequences of digits and blocks, and to learn and recall verbal and visual pairs, with immediate and delayed trials. In the Rey-Osterrieth Complex Figure Test (Meyers et al., 1996), participants were instructed to copy a detailed figure and to recall it after a delay of 30 min. In the Trail Making Test (Spreen and Strauss, 1998) subjects have to join with a pencil on a sheet different digits (part A) or alternating digits and letters (part B) as fast as possible in the right order. In the Regensburger Verbal Fluency Test (German version) (Aschenbrenner et al., 2000), subjects were asked to produce as many words as possible, according to a specific initial letter or a semantic category, or to switch between two letters or two categories. The short version of the Hamburg Wechsler Intelligence Scale for Adults-revised (German version) (Tewes, 1991) consists of five scales with different demands: completing figures (detecting missing details in pictures), vocabulary test (explaining the meaning of terms), mosaic test (replicating patterns with small cubes), finding similarities between different terms, and mathematical reasoning (mental arithmetic).

Statistical analyses

As primary analysis, a repeated-measures ANOVA (ANOVARM) of changes in adjusted MEP size after the PAS protocol was conducted with repeated factor “TIME” (baseline, post 0, post 15, post 30, post 60) and between-subject factor “GENOTYPE” (BDNF Met carriers, BDNF Val/Val carriers). To account for COMT and KIBRA carrier status, COMT and KIBRA genotype were entered into the model as additional between-subject factors. Student's t tests were run for post hoc comparisons if appropriate. Responders (grand average >1) and nonresponders to the PAS protocol were compared using χ2 tests.

Performance in the implicit learning tasks (AGL, SRTT) and in the neuropsychological test scores were compared using paired t tests (BDNF) or univariate ANOVA (BDNF/COMT/KIBRA interactions).

Correlations between PAS-induced plasticity and implicit learning tasks were assessed using Pearson's or Spearman's correlations, as appropriate.

Statistical analyses were done using SPSS 19 (PASW, SPSS, IBM), and levels of significance were set to p < 0.05.

Results

Neurophysiological baseline measurements

No differences were found between BDNF Val66Met genotype groups with regard to RMT, baseline MEPs, stimulation intensities or perceptual threshold (Table 4).

Table 4.

TMS parameters at baseline dependent on BDNF-genotype

PAS-induced plasticity

One subject (BDNF Met carrier) had to be excluded from further analysis, due to extreme outliers in post-MEP sizes and grand average (data >3 SDs from the mean). The matched subject (BDNF Val/Val carrier) was additionally excluded to counterbalance group sizes, leaving 30 subjects for group analyses.

For BDNF only, ANOVARM did not show a significant effect of BDNF genotype × TIME after PAS (F(3,84) = 2.02, p = 0.117); however, trend-wise larger MEPs at post 0, post 15, and post 60 in BDNF Val/Val carriers were noted compared with Met carriers (Fig. 3A).

Figure 3.

Effects of BDNF Val66Met alone and stratified for COMT Val158Met carrier status on cortical plasticity using PAS25. A, BDNF Val/Val carriers (triangles, dashed line) exhibited larger MEP sizes than BDNF Met carriers (circles, gray line) after LTP induction relative to baseline (pre), most pronounced at post 0 (ANOVARM, p = 0.117). B, In COMT Met/Met carriers only, BDNF Val/Val genotype was significantly associated with an increase in MEP sizes pre vs post 0 compared with BDNF Met carriers, indicating higher PAS-induced plasticity in the combined BDNF Val/Val + COMT Met/Met genotype group (p = 0.01). C, No differences were found within COMT Val carriers. Bars indicate SEM. *p < 0.01 according to ANOVARM post hoc t test.

When adding COMT Val158Met gene carrier status to the model, ANOVARM detected a significant interaction term BDNF × COMT × TIME (F(3,78) = 3.23, p = 0.027) (Fig. 3B,C). Post hoc t tests showed that in the group of homozygous COMT Met/Met carriers only, the BDNF Val/Val genotype was associated with a larger increase in MEP sizes at post 0, compared with BDNF Met carriers (p = 0.01; Fig. 3B). No differences were found within the group of COMT Val carriers (p > 0.05; Fig. 3C).

When adding KIBRA gene carrier status to the model, no significant effects were observed in combination with BDNF genotype, nor with both BDNF and COMT status (ANOVARM, p > 0.05, data not shown).

In the two BDNF groups, distribution of responders (grand average, >1) and nonresponders were not significantly different (χ2 (1) = 1.21, p > 0.05). Likewise, no group differences were found for COMT and KIBRA interactions (p > 0.05).

Implicit learning

Artificial grammar learning

No significant differences between the two BDNF genotype groups (n = 32) for performance in the AGL task were found: both groups had comparable scores (percentage correct) in the classification task, d′ values, and endorsement rates for grammatical and nongrammatical items (Table 5; Fig. 4A) (all p > 0.05).

Table 5.

Artificial grammar learning dependent on BDNF genotype

Figure 4.

Effects of BDNF Val66Met genotype alone and stratified for COMT Val158Met carrier status on artificial grammar learning. A, Endorsement rate, i.e., the percentage of items perceived as grammatical, was comparable between BDNF genotype groups for G and NG items. B, C, When also considering COMT carrier status, a significant BDNF × COMT interaction was found for NG (ANOVARM, p = 0.02, gray boxes), showing that in COMT Met/Met carriers (B), BDNF Val/Val carriers (gray circles) performed better than BDNF Met carriers (black triangles), and vice versa in COMT Val carriers (C) (n.s., all post hoc t tests p < 0.14). Bars indicate SEM.

However, adding COMT Val158Met carrier status revealed distinct behavioral differences, namely a significant BDNF X COMT interaction term for endorsement rate of nongrammatical items (ANOVA, F(28,1) = 5.97, p = 0.021) (Fig. 4B,C). Accordingly, BDNF Val/Val carriers that were also COMT Met/Met carriers showed higher task performance (as indicated by a lower tendency to perceive nongrammatical items as grammatical), compared with BDNF Met carriers of this group (Fig. 4B) (post hoc t test, p = 0.07). In contrast, in the group of COMT Val carriers, BDNF Val/Val carriers exhibited lower task performance (as indicated by higher endorsement rates, thus a higher tendency to perceive nongrammatical items as grammatical), compared with BDNF Met carriers (Fig. 4C) (post hoc t test, p = 0.014, n.s.).

No further interactions with COMT or KIBRA were found (ANOVA, all p > 0.05).

Serial reaction time task

Both BDNF genotype groups showed similar motor learning performance in the SRTT (n = 32), with slightly higher values for BDNF Val/Val carriers (Fig. 5A) (ANOVARM, p > 0.05). In COMT Met/Met carriers, an advantage of BDNF Val/Val carriers compared with BDNF Met carriers was noted (Fig. 5B), similar to results of AGL (see Artificial grammar learning), which failed to reach statistical significance though (ANOVA p > 0.05). No differences were found in COMT-Val carriers.

Figure 5.

Effects of BDNF Val66Met genotype alone and stratified for COMT Val158Met carrier status on motor learning performance. A, Considering BDNF genotype alone, both groups exhibited similar motor learning, i.e., mean reaction time differences between sequential and random stimuli of the SRTT. B, BDNF Val/Val carriers (gray circles) that were also COMT Met/Met carriers showed a trend for better performance compared with BDNF Met carriers (black triangles; ANOVA, n.s.). C, No differences were found in COMT Val carriers.

No interaction effects were found with regard to KIBRA and COMT × KIBRA carrier status (all p > 0.05).

Notably, SRTT performance correlated significantly with PAS-induced plasticity (post 0; Spearman‘s correlation, r = 0.52, p = 0.003) (Fig. 6). Performance in the AGL task did not correlate with PAS-induced plasticity (all p > 0.05).

Figure 6.

Correlation of motor learning performance and PAS-induced cortical plasticity. Mean reaction time differences between sequential and random stimuli (reaction time differences, x-axis) correlated significantly with the relative motor evoked potential (MEP, y-axis) at post 0 following the PAS25 protocol (r = 0.52, p = 0.003). Gray circles, BDNF Met carriers; black triangles, BDNF Val/Val carriers.

Neuropsychological performance

No significant differences were found between BDNF genotype groups on neuropsychological test performance (all t values <1.7, all p values >0.1) (Table 6).

Table 6.

Neuropsychological performance dependent on BDNF genotype

No significant associations with COMT or KIBRA were found (p > 0.05), except a BDNF × COMT interaction for the subtest immediate learning for verbal paired associates (F(1,28) = 8.21, p = 0.008), showing higher scores in BDNF Val/Val carriers compared with BDNF Met carriers in the group of COMT Val carriers (post hoc t test, t = 2.3, p = 0.032).

Discussion

We here present first-time evidence for a selective interaction effect of BDNF Val66Met and COMT Val158Met genotype on PAS-induced cortical plasticity in young females: while BDNF carrier status alone did not account for interindividual differences, plasticity was significantly higher in BDNF Val homozygotes if they were also COMT Met homozygotes, compared with BDNF Met carriers. These genotype-dependent effects were also found for grammar learning.

BDNF Val66Met and brain stimulation

In a key study by Kleim et al. (2006) (replicated by the same group, McHughen et al., 2010), increased MEP sizes and motor map representations after training of motor skills were found in young healthy Val/Val carriers using TMS. Enhanced plasticity in Val/Val carriers has been also observed using other brain stimulation protocols, i.e., TBS (Cheeran et al., 2008; Antal et al., 2010) and pharyngeal excitatory stimulation (Jayasekeran et al., 2011). An “advantage” of the BDNF Val allele with regard to plasticity has been explained by higher activity-dependent secretion of BDNF in Val carriers (Egan et al., 2003), as BDNF is implicated in LTP (Poo, 2001). However, opposite effects were also reported, i.e., increased stimulation-induced plasticity in Met carriers using tDCS (Antal et al., 2010) and in response to inhibitory 1 Hz repetitive (r) TMS (Jayasekeran et al., 2011). Notably, recent studies indicated that an overexpression of BDNF may negatively affect plasticity and learning, due to hyperactivation of inhibitory pathways (for review, see Cunha et al., 2010). In the present study, we did not find significant effects on cortical plasticity induction using the PAS protocol for BDNF genotype alone. This finding is in line with previous studies reporting null findings (quadripulse stimulation; Nakamura et al., 2011) or trends only (PAS; Cheeran et al., 2008) if looking at BDNF genotype and measures of cortical plasticity.

The divergent results might be in part explained by methodological differences in stimulation techniques. For example, as the BDNF Val66Met polymorphism influences secretion of BDNF but does not influence its function (Egan et al., 2003), a longer duration of stimulation may be able to induce sufficient secretion of BDNF even in BDNF Met allele carriers. Thus, a stimulation protocol of 30 min, as used in our own study and in that of Nakamura et al. (2011), might have diminished potential genotype-dependent differences, in contrast to shorter protocols of 190 s or 13 min. Other intrinsic differences in protocols used to induce cortical plasticity in humans, e.g., stimulation mode or pulse frequency (Ziemann, 2004), may additionally contribute to differential findings between studies. We chose the PAS protocol (Stefan et al., 2000) for the present study as a frequently used and sensitive tool to investigate cortical plasticity and its relation to age and cognitive function in healthy subjects and neurological patients (Morgante et al., 2006; Weise et al., 2006; Frantseva et al., 2008; Muller-Dahlhaus et al., 2008; Castel-Lacanal et al., 2009; Fathi et al., 2010; Zeller et al., 2010; List et al., 2011, 2012).

Interaction effects of BDNF Val66Met and COMT Val158Met

Additional possible confounders in this context are interactions of BDNF Val66Met with other SNPs that are suggested to modulate LTP, such as COMT and KIBRA. For example, the supposed BDNF genotype-dependent differences in BDNF concentrations and thus plasticity induction may be masked or evened out by COMT Val158Met-dependent differences in dopaminergic transmission, as dopamine is also known to affect LTP and neuronal plasticity (Wolf et al., 2003; Li et al., 2010). Indeed, after stratifying for COMT Val158Met genotype, we observed significantly increased PAS-induced cortical plasticity in BDNF Val homozygotes if subjects were also COMT Met homozygotes, compared with BDNF Met carriers. The COMT Met allele has been linked to less enzymatic activity of COMT and thus higher prefrontal dopamine (Lindenberger et al., 2008), yet whether COMT Val158Met carrier status influences human cortical plasticity has not been assessed so far. Several studies reported better executive function and memory performance in COMT Met carriers; however, others could not confirm these results (for review, see Witte and Flöel, 2011). Our findings for cortical plasticity indicated that the BDNF Val/Val genotype is beneficial only in combination with the low-activity COMT Met allele. Notably, group differences were most prominent immediately after PAS, yet diminished at 15, 30, and 60 min of follow-up. Bearing in mind the duration of the stimulation protocol (i.e., 30 min in the present study, see BDNF Val66Met and brain stimulation, above), the BDNF by COMT interaction found here might have affected immediate LTP-like plasticity only (for a detailed discussion, see also Stefan et al., 2000). A similar time course was also observed by Cheeran et al. (2008) for rTMS-induced effects, who found a significant BDNF Val66Met group difference immediately after a 15 min rTMS protocol, but not 10 min later. However, other stimulation protocols yielded longer-lasting group differences for BDNF Val66Met [e.g., following iTBS/cTBS (Cheeran et al., 2008) and tDCS (Antal et al., 2010)].

Immediate LTP affects the probability of transmitter release and has been suggested to be crucially important in “everyday memory” (Kandel, 2000; Wang et al., 2010). Thus, higher LTP-like plasticity directly after PAS should translate into behavioral advantages. Accordingly, the BDNF by COMT genotype group with increased LTP showed higher immediate learning success in the grammar task, and, although not reaching statistical significance, in the motor learning task. Moreover, LTP induction directly after PAS correlated with motor learning performance (Fig. 6). This is in line with previous studies that demonstrated higher LTP-like plasticity, induced with a 30 or 15 min PAS protocol, to be correlated with better cognitive functions in patients with cerebral white matter lesions (List et al., 2011) and better motor functions in patients with multiple sclerosis (Zeller et al., 2010). These findings, together with studies on LTP and learning in animal experiments (deToledo-Morrell et al., 1988; Tsien et al., 1996; Rosenzweig and Barnes, 2003), support the hypothesis that higher LTP, or LTP-like plasticity, correlates with better learning success.

One possible explanation for the genotype interaction effect observed in the current study might be that the positive effects of the BDNF-Val and COMT-Met alleles on BDNF concentrations and dopamine levels are small; therefore, subtle differences could only be seen when pooling the two genotypes to a “best of” combination. However, a genuine interaction effect may also be present and has been reported by previous studies for BDNF and COMT (Nagel et al., 2008). This is further supported by the fact that the high-activity BDNF Val homozygotes did not uniformly show better performance; rather, if paired with the “disadvantageous” lower dopamine COMT Val carriers, they showed significantly lower performance in grammar learning and immediate learning of verbal paired associates, compared with BDNF Met carriers of this group. Here, it may be hypothesized that reduced dopaminergic transmission in combination with high levels of BDNF would provoke rather adverse actions with regard to plasticity, e.g., hyperactivation of inhibitory pathways (Cunha et al., 2010). However, the exact underlying mechanisms remain elusive and have to be addressed in future experimental studies, e.g., using advanced molecular-genetic tools to explicitly study the effects of high BDNF concentration and low dopamine levels on neuronal activity. Notably, higher activity of COMT might not always lead to behavioral disadvantages, for example, regarding emotional tasks with higher demands of flexible versus stable neuronal processing within prefronto-limbic circuits (for a detailed discussion, see Witte and Flöel, 2011).

A slightly different interaction effect between BDNF and COMT has been reported for executive functions in older adults (Nagel et al., 2008). Yet in young adults in the Nagel et al. (2008) study, working memory performance for the BDNF by COMT interaction showed a similar pattern compared with that of our subjects in the grammar learning task (compare Nagel et al., 2008, their Fig. 2D, nonsignificant results). In sum, the impact of BDNF and COMT genotype on cortical plasticity and learning ability appears to be highly complex and remains to be fully elucidated, yet interactions between these two SNPs explained significantly more variance than considering BDNF genotype alone (this study; Nagel et al., 2008).

Other learning-relevant SNPs

For rs17070145 in the KIBRA gene, behavioral results have been mixed, ranging from beneficial effects of the T-allele (Papassotiropoulos et al., 2006; Almeida et al., 2008; Schaper et al., 2008; Vassos et al., 2010; Yasuda et al., 2010; Kauppi et al., 2011), to null findings (Need et al., 2008), to even poorer performance in T-allele carriers (Nacmias et al., 2008; Wersching et al., 2011). In line with Need et al. (2008), we did not find significant effects for KIBRA genotype on PAS-induced cortical plasticity, implicit learning, or neuropsychological performance. As previous studies with positive results included sample sizes of >60 subjects, the proposed effects of KIBRA genotype might be too small to yield significant results in our sample, especially with regard to triple interactions including BDNF and COMT genotype. Upcoming studies including larger sample sizes have to address this issue.

Strength and limitations

One limitation is the sample size of n = 32, which might impact the statistical power to detect small genotype-dependent effects (see also Other learning-relevant SNPs, above). However, we used a sensitive measure, PAS-induced cortical plasticity, as primary outcome, and focused on one polymorphism, BDNF. For BDNF, subjects were carefully pairwise matched for age and education, and sample size was similar to previous studies using cortical plasticity measurements (Kleim et al., 2006; Cheeran et al., 2008; Antal et al., 2010; Nakamura et al., 2011). Still, no significant effect of BDNF alone on cortical plasticity or learning ability was detected. Moreover, as gender is a crucial factor that might interact with BDNF, COMT, and KIBRA (Wersching et al., 2011; Witte and Flöel, 2011), we included females only to fully avoid gender effects. Second, including only females entails that results cannot be directly transferred to males. Third, ceiling effects in our sample of healthy young subjects with regard to neuropsychological test scores may mask genotype-driven effects of BDNF, COMT, and KIBRA for these outcome parameters. Note that participants did not reach ceiling levels for cortical plasticity or implicit learning.

Conclusion

Our findings for the first time demonstrate an interaction of BDNF Val66Met and COMT Val158Met genotype on human cortical plasticity, explaining more interindividual variance than BDNF genotype alone. Moreover, they show that genotype-related differences in neurophysiology also translate into behavioral differences. As a consequence, upcoming studies should account for possible interaction effects between different learning-relevant SNPs, to further clarify the ambiguous effects of single genotypes in the past. In sum, our findings might contribute to a better understanding of the mechanisms underlying interindividual differences in cognition, which might help to improve preventive and therapeutic strategies against neurological or psychiatric diseases (Mayor, 2007).

Footnotes

    • Received December 2, 2011.
    • Revision received January 12, 2012.
    • Accepted January 19, 2012.
  • This work was supported by grants from the Deutsche Forschungsgemeinschaft (Fl 379-4/2, Fl 379-8/1, and DFG-Exc 257), the Else-Kröner Fresenius Stiftung (2009-141), the Interdisciplinary Center of Clinical Research Münster (Floe3/004/08), and the Bundesministerium für Bildung und Forschung (FKZ 0315673A and 01EO0801).

  • Correspondence should be addressed to Dr. Agnes Flöel, Universitätsmedizin Charité, CCM, Department of Neurology, Charitéplatz 1, 10117 Berlin, Germany. agnes.floeel{at}charite.de

References

  1. ↵
  2. ↵
  3. ↵
  4. ↵
  5. ↵
  6. ↵
  7. ↵
  8. ↵
  9. ↵
  10. ↵
  11. ↵
  12. ↵
  13. ↵
  14. ↵
  15. ↵
  16. ↵
  17. ↵
  18. ↵
  19. ↵
  20. ↵
  21. ↵
  22. ↵
  23. ↵
  24. ↵
  25. ↵
  26. ↵
  27. ↵
  28. ↵
  29. ↵
  30. ↵
  31. ↵
  32. ↵
  33. ↵
  34. ↵
  35. ↵
  36. ↵
  37. ↵
  38. ↵
  39. ↵
  40. ↵
  41. ↵
  42. ↵
  43. ↵
  44. ↵
  45. ↵
  46. ↵
  47. ↵
  48. ↵
  49. ↵
  50. ↵
  51. ↵
  52. ↵
  53. ↵
  54. ↵
  55. ↵
  56. ↵
  57. ↵
  58. ↵
  59. ↵
  60. ↵
  61. ↵
  62. ↵
  63. ↵
  64. ↵
  65. ↵
  66. ↵
  67. ↵
  68. ↵
  69. ↵