Skip to main content

Umbrella menu

  • SfN.org
  • eNeuro
  • The Journal of Neuroscience
  • Neuronline
  • BrainFacts.org

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
  • ALERTS
  • FOR AUTHORS
    • Preparing a Manuscript
    • Submission Guidelines
    • Fees
    • Journal Club
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
  • SfN.org
  • eNeuro
  • The Journal of Neuroscience
  • Neuronline
  • BrainFacts.org

User menu

  • Log in
  • Subscribe
  • My alerts
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • Subscribe
  • My alerts
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
  • ALERTS
  • FOR AUTHORS
    • Preparing a Manuscript
    • Submission Guidelines
    • Fees
    • Journal Club
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Articles, Behavioral/Systems/Cognitive

Stress-Related Methylation of the Catechol-O-Methyltransferase Val158 Allele Predicts Human Prefrontal Cognition and Activity

Gianluca Ursini, Valentina Bollati, Leonardo Fazio, Annamaria Porcelli, Luisa Iacovelli, Assia Catalani, Lorenzo Sinibaldi, Barbara Gelao, Raffaella Romano, Antonio Rampino, Paolo Taurisano, Marina Mancini, Annabella Di Giorgio, Teresa Popolizio, Andrea Baccarelli, Antonio De Blasi, Giuseppe Blasi and Alessandro Bertolino
Journal of Neuroscience 4 May 2011, 31 (18) 6692-6698; DOI: https://doi.org/10.1523/JNEUROSCI.6631-10.2011
Gianluca Ursini
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Valentina Bollati
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Leonardo Fazio
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Annamaria Porcelli
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Luisa Iacovelli
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Assia Catalani
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lorenzo Sinibaldi
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Barbara Gelao
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Raffaella Romano
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Antonio Rampino
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paolo Taurisano
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marina Mancini
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Annabella Di Giorgio
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Teresa Popolizio
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea Baccarelli
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Antonio De Blasi
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Giuseppe Blasi
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alessandro Bertolino
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

DNA methylation at CpG dinucleotides is associated with gene silencing, stress, and memory. The catechol-O-methyltransferase (COMT) Val158 allele in rs4680 is associated with differential enzyme activity, stress responsivity, and prefrontal activity during working memory (WM), and it creates a CpG dinucleotide. We report that methylation of the Val158 allele measured from peripheral blood mononuclear cells (PBMCs) of Val/Val humans is associated negatively with lifetime stress and positively with WM performance; it interacts with stress to modulate prefrontal activity during WM, such that greater stress and lower methylation are related to reduced cortical efficiency; and it is inversely related to mRNA expression and protein levels, potentially explaining the in vivo effects. Finally, methylation of COMT in prefrontal cortex and that in PBMCs of rats are correlated. The relationship of methylation of the COMT Val158 allele with stress, gene expression, WM performance, and related brain activity suggests that stress-related methylation is associated with silencing of the gene, which partially compensates the physiological role of the high-activity Val allele in prefrontal cognition and activity. Moreover, these results demonstrate how stress-related DNA methylation of specific functional alleles impacts directly on human brain physiology beyond sequence variation.

Introduction

Single nucleotide polymorphisms (SNPs) can create or abolish CpG dinucleotides, which constitute a site for DNA methylation, an epigenetic mechanism for regulation of gene expression. Earlier studies in rodents have indicated that DNA methylation is sensitive to environmental stressful exposures early in development and later in life (Jirtle and Skinner, 2007; Miller and Sweatt, 2007; Szyf et al., 2008; McGowan et al., 2009; Murgatroyd et al., 2009; Miller et al., 2010). Moreover, other studies in mice have reported that memory formation and maintenance are sensitive to epigenetic programming by early life stress and, more specifically, to DNA methylation of specific neurons (Miller and Sweatt, 2007; Murgatroyd et al., 2009; Miller et al., 2010). Experiments in human monozygotic twins have suggested a relationship between dynamic epigenetic mechanisms and environment (Fraga et al., 2005; Wong et al., 2010). However, to date, there is no evidence in humans of the relationship among stress, DNA methylation, and brain activity during memory.

In humans, prefrontal dopamine (DA) and stress are tightly related with working memory (WM) behavior and related prefrontal cortex (PFC) activity, suggesting that lower DA and greater stress are associated with reduced prefrontal efficiency (Hains and Arnsten, 2008). Catechol-O-methyltransferase (COMT), a key enzyme for inactivation of prefrontal DA (Gogos et al., 1998), contains a SNP (rs4680, G→A, Val→Met), in which the ancestral Val allele is associated with greater enzyme activity, blunted stress responses, and greater prefrontal activity during WM, that is, reduced efficiency (Egan et al., 2001; Chen et al., 2004; Papaleo et al., 2008; van Winkel et al., 2008; Mier et al., 2010; Walder et al., 2010). Notably, rs4680 in COMT exon 4 also abolishes a CpG site so that each Val allele has one CpG site and the Met allele has none.

The objective of the present study was to evaluate the relationship among life stress, COMT methylation at the rs4680 site (C2) measured in peripheral blood mononuclear cells (PBMCs), WM performance, and related prefrontal activity measured by functional magnetic resonance imaging (fMRI) in healthy humans. Moreover, we addressed the effects of DNA methylation on COMT expression in PBMCs to determine a molecular mechanism for the in vivo findings. Finally, we addressed in a sample of rats the question of whether COMT methylation patterns in PBMCs and in brain are correlated.

Materials and Methods

Subjects.

Healthy subjects (n = 84; 32 males; mean age ± SD, 25.90 ± 5.58 years) entered the study. All participants were unrelated, and all were Caucasian; all provided written informed consent. Protocols and procedures were approved by the local Institutional Review Board. Subjects underwent the Structured Clinical Interview for DSM-IV to exclude any Axis I psychiatric disorder. None had a history of significant drug or alcohol abuse (no active drug use in the past 6 months), head trauma with loss of consciousness, or significant medical illness.

Methylation analysis and genotyping.

DNA was extracted from PBMCs using standard procedures and bisulfite treated as described previously (Bollati et al., 2007). Methylation analysis was performed with pyrosequencing (primer sequences available on request) and focused on CpG methylation sites at COMT rs4680 and within the gene promoter. A consensus LINE-1 sequence was also analyzed to estimate global DNA methylation (Bollati et al., 2007). In the COMT promoter, we analyzed methylation of three CpGs, two of which had been already identified as differentially methylated (Murphy et al., 2005; Dempster et al., 2006; Mill et al., 2006). Methylation was expressed as percentage of methylated cytosines divided by the sum of methylated and unmethylated cytosines (%5mC) (Tost and Gut, 2007). SNP rs4680 was genotyped using the pyrosequencing assay designed to interrogate percentage of methylation in this region. In addition, genotypes were double-checked by DNA direct sequencing (primer sequences available on request).

Stress assessment.

Because DNA methylation is sensitive to early and later environmental stimuli (Jirtle and Skinner, 2007; McGowan et al., 2009; Murgatroyd et al., 2009; Miller et al., 2010), we set out to determine a score of lifetime exposure to stress capturing both early and later stressors. To determine this score, we used obstetric complications (OCs) and stressful life events (SLEs) questionnaires as follows. On the basis of interviews administered to mothers using a standard questionnaire developed from other published reports (McIntosh et al., 2002), OCs are referred to as “somatic complications and conditions occurring during pregnancy, labor-delivery and the neonatal period” experienced as an offspring with special focus on the CNS (McNeil et al., 1994) and were determined to obtain a measure of early life stress of the subjects. OC data were rated using the McNeil–Sjöström scale for obstetric complications (McNeil et al., 1994). Consistent with earlier studies (Preti et al., 2000; Verdoux et al., 2002), we used this scale to assign each subject a binary OC score equal to 0 or 1, depending, respectively, on absence or presence of at least one serious OC (severity score ≥ 4 on a scale of 1–6). SLEs were assessed with the aid of a life-history calendar (Caspi et al., 2003). The number of stressful events was normalized to age to obtain an SLE score, which was then converted to a binary measure so that it could be summed to the OC score. Finally, SLE and OC scores were summed to obtain a categorical lifetime stress score, equal to 0 (lower stress), 1, or 2 (greater stress).

WM task.

During fMRI, all subjects completed a blocked paradigm of the n-back task (Bertolino et al., 2010). Briefly, n-back refers to how far back in the sequence of stimuli the subject has to recall. The stimuli consisted of numbers (1–4) shown in random sequence and displayed at the points of a diamond-shaped box. There was a visually paced motor task that also served as a non-memory-guided control condition (0-back) that simply required subjects to identify the stimulus currently seen. In the WM condition, the task required recollection of a stimulus seen two stimuli previously (2-back) while continuing to encode additionally incoming stimuli.

fMRI data acquisition.

Blood oxygen level-dependent (BOLD) fMRI was performed on a GE Signa 3T scanner (General Electric), equipped with a standard quadrature head coil as previously described (Bertolino et al., 2010). A gradient-echo planar imaging sequence (repetition time, 2000 ms; echo time, 28 ms; 20 interleaved axial slices; thickness, 4 mm; gap, 1 mm; voxel size, 3.75 mm (isotropic); flip angle, 90°; field of view, 24 cm; matrix, 64 × 64) was used to acquire 120 volumes while subjects performed the WM task. The first four scans were discarded to allow for a T1 equilibration effect.

fMRI data analysis.

Analysis of the fMRI data was completed using statistical parametric mapping (SPM5; http://www.fil.ion.ucl.ac.uk/spm). All preprocessing was performed as described previously (Bertolino et al., 2010). A boxcar model convolved with the hemodynamic response function at each voxel was generated. In the first-level analysis, linear contrasts were computed, producing a t statistical map at each voxel for the 2-back condition, assuming the 0-back condition as a baseline. All individual contrast images were entered in a second-level random-effects analysis. A one-sample t test (p < 0.005; k = 5), performed on the whole group, was used to assess the functional network involved in performance of the WM task. A multiple regression model, entering percentage of methylation of COMT C2 and lifetime stress score as predictors, was used to evaluate the main effects as well as the interaction between C2 methylation and stress in Val/Val subjects. Because of our strong hypothesis about prefrontal behavior and activity, the statistical threshold of p < 0.05 after familywise error (FWE) small-volume correction for multiple comparisons was used to identify significant responses in bilateral PFC, using three 10 mm radius spheres centered on coordinates published in previous WM studies (Wager and Smith, 2003; Habeck et al., 2005; Tan et al., 2007). Activation outside these regions of interest was FWE whole-brain corrected (p < 0.05; k = 5). Similar analyses were also performed for the other genotype groups and for the whole sample when evaluating the potential association with methylation of the COMT promoter or of LINE-1. For anatomical localization, we converted statistical maxima of activation to conform to the standard space of Talairach and Tournoux.

COMT m-RNA levels by quantitative reverse transcription-PCR.

We analyzed COMT mRNA expression in 46 healthy individuals (14 Val/Val, 18 Val/Met, 14 Met/Met). MB-COMT (long transcript) and total-COMT (both long and short transcripts) mRNA levels were assessed using the comparative CT method with β-actin as reference (control) gene (Gao et al., 2010) (primers available on request). Real-time quantitative PCR was performed using ABI PRISM 7700 Sequence Detection System in combination with continuous SYBR Green detection (Applied Biosystems).

Western blotting.

We further evaluated the relationship between C2 methylation and COMT expression (both MB- and S-COMT proteins) in 21 Val-carrier healthy individuals (7 Val/Val, 14 Val/Met), given that only the Val allele has a CpG dinucleotide at C2. PBMCs were isolated as previously described (Castri et al., 2007). Aliquots (2 μl) of samples were used for the protein determination by Bio-Rad Protein Assay kit (Bio-Rad); 80 μg of total proteins were separated by SDS-PAGE, using 15% polyacrylamide gels, blotted onto nitrocellulose and probed overnight using commercial anti-COMT antibody (1:7000, Millipore Bioscience Research Reagents, Cat. AB5873). Membranes were then washed and reprobed with anti-tubulin (1:5000, Millipore, Cat. 05-661) to verify equal protein loading. Antibody against COMT was able to recognize both S- and MB-COMT. Blots were then incubated in horseradish peroxidase-conjugated secondary antibodies and target proteins visualized by ECL detection (Pierce), followed by quantification by Quantity One software (Bio-Rad). Optical density values were normalized to tubulin for variation in loading and transfer.

Methylation analysis in rat lymphocytes and PFC.

Adult male Wistar rats (Charles River) of 4 months of age, weighing 491 ± 11 g, were housed two per cage in a temperature-controlled room at 21 ± 1°C, with 50 ± 5% of humidity and a 12 h light–dark period (lights on 7:00 A.M. to 7:00 P.M.). Food (Standard Diet 4RF21; Charles River) and tap water were provided ad libitum. Animals were killed by decapitation between 9:00 and 11:00 A.M., trunk blood (5 ml) was collected in tubes containing EDTA (2 mg/ml blood), brains were removed from the skull, and the PFC was rapidly dissected using a standard procedure (Gearhart et al., 2006), frozen on dry ice, and kept at −80°C until processed. Similarly, the hippocampus and striatum were also dissected from the same rats as control regions. Both anatomical regions are sensitive to stress like the PFC, but they have a different relationship with rs4680 genotype in response to several cognitive paradigms (Smolka et al., 2005; Bertolino et al., 2006; Yacubian et al., 2007; Krugel et al., 2009). Consistent with these functional imaging data, both regions differentially express COMT when compared with PFC, thus potentially implicating different mechanisms for gene expression (Sesack et al., 1998; Matsumoto et al., 2003). Analysis of rat LINE-1 and COMT methylation was performed with pyrosequencing (primer sequences available on request).

Results

The relationship among COMT methylation, rs4680 genotype, and stress

As expected, methylation of the C2 rs4680 site was affected by rs4680 genotype (ANOVA: F(2,81) = 7336.8, p < 0.0001; Val/Val > Val/Met > Met/Met) given that only the Val allele has a CpG site in this position. Therefore, separate analyses were conducted in the three genotype groups. ANOVA demonstrated that Val/Val subjects with greater stress scores have significantly reduced methylation in C2 compared with subjects with lower scores (F(2,20) = 4.1, p = 0.03; post hoc with Fisher's least significant difference: “2” < “0,” p = 0.01; “2” < “1,” p = 0.06) (Fig. 1A). Consistent with earlier studies (Dempster et al., 2006), we also found an effect of rs4680 genotype on methylation of the COMT promoter sites, indicating that subjects carrying two Met alleles have greater methylation than Val/Met subjects, who, in turn, have greater methylation than Val/Val subjects (ANOVA: F(2,81) = 11.2, p < 0.01; data not shown). However, no significant relationship was found between methylation of sites in the COMT promoter or in LINE-1 and stress in the whole sample (ANOVAs: F values ≤ 0.1, p values > 0.9) and in any genotype group (ANOVAs: F values < 2.0, p values > 0.1).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

C2 methylation in Val/Val subjects: relationship with stress and WM behavioral performance. A, Mean ± 0.95 confidence interval of the effect of stress on methylation of COMT C2 rs4680 in Val/Val homozygotes—subjects with greater stress scores have reduced methylation. B, Scatter plot of the correlation between methylation of COMT C2 rs4680 and WM accuracy—Val/Val subjects with greater methylation have greater accuracy (see Results for statistics).

The relationship among COMT methylation, stress, and prefrontal cognition and activity

As expected given that DNA methylation may reduce gene expression, C2 methylation in Val/Val subjects was positively correlated with 2-back WM accuracy (n = 19, Spearman's ρ = 0.59, p = 0.007), with greater methylation predicting greater percentage of correct responses (Fig. 1B). No statistically significant relationship was found between WM accuracy and stress (F(2,16) = 2.1, p > 0.1). Consistent with earlier fMRI studies demonstrating an inverse relationship between Val alleles and PFC activity during WM, multiple regression in SPM5 demonstrated the following in Val/Val subjects: a negative correlation between C2 methylation and bilateral PFC activity during WM [x = −38, y = 23, z = 9, Brodmann area 45/13 (BA45/13), k = 11, Z = 3.38, pFWE-corrected = 0.01; x = 45, y = 15, z = 3, BA47, k = 7, Z = 3.54, pFWE-corrected = 0.009]; a positive correlation between stress scores and activity in left PFC (x = −45, y = 48, z = −3, BA10/47, k = 6, Z = 2.9, pFWE-corrected = 0.04); and an interaction between C2 methylation and stress scores on bilateral prefrontal activity during WM, suggesting that greater stress and lower methylation are related to greater activity (i.e., inefficient; x = −45, y = 48, z = −3, BA10/47, k = 11, Z = 3.26, pFWE-corrected = 0.02; x = 49, y = 11, z = 6, BA44/45, k = 8, Z = 3.18, pFWE-corrected = 0.03) (Fig. 2). On the other hand, no significant relationship was found between methylation of sites in the COMT promoter or in LINE-1 and stress, WM performance, and related PFC activity in any other genotype group or in the whole sample (p values > 0.1). Lack of a correlation among COMT promoter methylation, stress, and prefrontal cognition or activity suggests that methylation of the rs4680 CpG site and of the promoter respond to different mechanisms, at least in the context of the relationship with stress and WM. Finally, no significant relationship between activity in other brain regions, including parietal cortex, and methylation of COMT or stress was evident at the statistical threshold identified a priori (data not shown).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Interaction between C2 methylation, stress scores, and prefrontal activity during WM in Val/Val subjects. A, Coronal section of the interaction between stress scores and C2 methylation on BOLD fMRI response in PFC. B, Scatter plot of the interaction showing that subjects with greater stress scores and lower methylation have greater activity in PFC, i.e., less efficiency (see Results for statistics).

The relationship between methylation of Val158 and COMT expression

To determine a molecular mechanism for the in vivo findings, we evaluated the effects of methylation on COMT expression in PBMCs by quantitative real-time reverse transcription-PCR and Western blotting analysis. As predicted, C2 methylation in Val/Val subjects was inversely correlated with MB-COMT (the isoform more expressed in brain) mRNA expression (n = 14, ρ = −0.57, p = 0.03) (Fig. 3) and with levels of MB-COMT protein (n = 7, ρ = −0.78, p = 0.03), whereas no relationship was found with S-COMT protein levels (n = 7, ρ = −0.14, p = 0.7), further suggesting the specificity of this relationship. In Val/Met subjects, we found no relationship between C2 methylation and COMT mRNA or protein levels (p values > 0.1).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Scatter plot of the correlation between C2 methylation and t scores of MB-COMT expression levels (transformed −ΔΔCt) in PBMCs of Val/Val subjects; greater methylation is correlated with lower expression (see Results for statistics).

Correlation of COMT methylation in brain and PBMCs

Because at present it is impossible to determine methylation in the brain of living healthy humans and we have evaluated DNA methylation from PBMCs, we addressed in a sample of rats the question of whether COMT methylation patterns in PBMCs and in brain are correlated. Moreover, because rats do not have the Val158Met mutation, we analyzed methylation on a CpG (location: chromosome 11: 84,577,576-7) in the COMT region (Leu151) corresponding to the region where the human rs4680 site is found. This analysis revealed a direct correlation between methylation in PBMCs and in PFC (n = 16, ρ = 0.49, p = 0.05), indicating that COMT methylation in PBMCs can be used as a proxy for COMT methylation in PFC (Fig. 4), as also suggested by others (Murphy et al., 2005). However, methylation in PBMCs was not significantly correlated with methylation in hippocampus (n = 16, ρ = 0.33, p = 0.2) or in striatum (n = 15, ρ = −0.13, p = 0.6). Finally, in the same rats, LINE-1 methylation in PBMCs and PFC were not significantly correlated, further suggesting a relatively specific relationship for COMT.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Scatter plot of the correlation between methylation of COMT in PBMCs and in PFC of a group of 16 rats—methylation analysis on a CpG (chromosome 11: 84,577,576-7) in the COMT region (Leu151) corresponding to the region where the human rs4680 site is found shows a direct correlation between methylation in PBMCs and in PFC (see Results for statistics).

Discussion

This study demonstrates for the first time in living humans that environment-related methylation of alleles within functional SNPs is an important mechanism for regulation of gene expression and behavior-related brain activity. Earlier studies have indicated that DNA methylation, prefrontal DA levels, and prefrontal activity are modulated by stress (Arnsten and Goldman-Rakic, 1998; Pani et al., 2000; Arnsten, 2009). Moreover, a series of earlier studies have demonstrated that COMT is a key factor in prefrontal DA inactivation (Gogos et al., 1998). Furthermore, the enzyme containing the Val158 allele has higher activity (Lotta et al., 1995), presumably leading to lower prefrontal DA. Consistently, the Val allele has been repeatedly associated with inefficient prefrontal activity, that is, greater activity for a given level of behavioral performance (Egan et al., 2001; Bertolino et al., 2004; Mier et al., 2010). Other studies have also indicated that the Val158 allele is a protective factor against stress, suggesting specific behavioral correlates of pleiotropic effects (Mattay et al., 2003; Papaleo et al., 2008; van Winkel et al., 2008; Walder et al., 2010). Consistent with this earlier body of work, here we have demonstrated that methylation of the Val allele in homozygous subjects is sensitive to stress; is correlated with WM behavioral performance; and interacts with stress, modulating PFC responses during WM performance such that greater stress and lower methylation are correlated with more inefficient prefrontal activity. These in vivo results are consistent with the role of DNA methylation in gene silencing and alternative transcription, even when outside of CpG islands (Jirtle and Skinner, 2007; Szyf et al., 2008; Irizarry et al., 2009; Maunakea et al., 2010). In fact, methylation of the Val allele in rs4680 is inversely related to MB-COMT expression and protein levels in Val/Val subjects, suggesting that silencing of the gene partially compensates the physiological role of the high-activity Val allele in prefrontal DA inactivation. This interpretation is also consistent with previous findings indicating that Val behaves as a high-activity allele when found on a high genetic expression background, but not otherwise (Meyer-Lindenberg et al., 2006; Nackley et al., 2006). In other words, the in vivo effects on memory performance and brain activity of the high-activity ancestral allele can be partially compensated by greater methylation, which is related to reduced expression of the COMT isoform (MB-) mainly responsible for inactivation of prefrontal DA.

It is interesting to note that no significant relationship among COMT rs4680 methylation, stress, prefrontal cognition, and COMT expression was found in Val/Met subjects. These negative results suggest that C2 methylation is not sufficient to modulate COMT function in heterozygous subjects, probably because methylation of both alleles is necessary to modify transcription and translation. In fact, COMT function in heterozygote subjects is associated with expression and activity of both alleles (Chen et al., 2004). The presence of one Met allele in heterozygotes may compensate for the effects related to methylation of the Val allele on expression and/or activity of COMT, preventing demonstration of a significant correlation. Furthermore, we cannot exclude the possibility that other CpGs modulate function of the chromosome containing a Met allele and might be sensitive to stress. For this reason, future studies of COMT methylation in Val/Met subjects should take into account the possibility of identifying and characterizing other CpGs present on chromosomes carrying the Met allele in addition to C2, which is present only on the Val allele.

The mechanisms correlating methylation of the Val allele in rs4680 with COMT transcription are not known. However, some speculations are possible. Cytosine methylation is most commonly associated with transcriptional repression mediated by direct interference with binding of transcription factors and/or through recruitment of methylated DNA-binding proteins and repressive complexes (Meaney and Ferguson-Smith, 2010). Consistently, bioinformatic tools to perform gene regulation analyses and search of complex regulatory patterns (http://www.genomatix.de/) suggest that rs4680 (G→A, Val→Met) leads to loss of a binding site for AHRARNT.01, a transcription factor belonging to the AHR/ARNT heterodimers and AHR-related factors family. Therefore, methylation of rs4680 CpG (C2) on the Val allele might affect gene expression through interference with binding of this transcription factor, potentially involved in several functions, including response to environmental toxicants and diverse developmental processes, such as neuronal differentiation (Oesch-Bartlomowicz et al., 2005; Sartor et al., 2009). Interestingly, AHR agonists at the lowest effective dose have been reported to induce an increase in MB-COMT expression in rat brain cortex (Desaulniers et al., 2005), and the induction of gene expression by AHR/ARNT has been associated with modification of specific chromatin marks (Schnekenburger et al., 2007). Consistent with these speculations and with earlier studies demonstrating that intragenic methylation has a major role in regulating alternative promoters within gene bodies (Maunakea et al., 2010), SNP rs4680 also seems to be located within a COMT alternative promoter (name: GXP_2244264; locus: COMT/GXL_744584; position: chromosome 22: 19950722–19951322, 601 bp; transcripts: ENST0000428707). Until further experiments are specifically performed, these contentions about the mechanisms by which methylation of rs4680 is correlated with COMT expression must be considered speculative.

Some limitations of our study concerning the assessment of stress and methylation have to be acknowledged. To evaluate the potential relationship between stress and methylation, we used a categorical lifetime stress score. This is indeed a crude measure of stress that does not consider all nuances of experiencing stress. The main reason why we used this crude measure of stress is that methylation is subject to environmental modifications that can take place both early and later in the life of an individual. Despite its simplistic nature, our measure does take into account early and later events that can be correlated with methylation. Moreover, further analyses separating OCs and SLEs also indicate that both scores correlate with methylation in Val homozygotes (results available on request). Despite all limitations intrinsic to such an oversimplification of measuring stress, our results suggest that the relationship with methylation has behavioral, physiological, and molecular implications.

Another limitation is that we have evaluated DNA methylation from PBMCs rather than directly in the brain, which is at present impossible in living healthy humans. Because epigenetic profiles can differ across tissues and cell types, this is an unavoidable limitation of our study. However, we believe that this limitation is tempered by our experiments in rats demonstrating that COMT methylation patterns in PBMCs and in PFC are correlated, suggesting that COMT methylation in PBMCs can be used as a proxy for COMT methylation in PFC. These results are also consistent with another report (Murphy et al., 2005). On the other hand, methylation in PBMCs did not correlate with methylation in hippocampus or striatum. Albeit weak, the correlation in the hippocampus was in the same direction as in PFC. In the striatum, where expression of COMT is very low, there was no relationship whatsoever; if anything, it was in the opposite direction. We believe that these further data suggest that other factors in addition to methylation of this CpG may be involved in regulation of MB-COMT expression in hippocampus and striatum. In other words, COMT methylation may have an important role in gene expression in PBMCs and in PFC, less so in other brain regions where it is expressed. Consistently, Li et al. (2010) have recently demonstrated in mice that B2 SINE, a transposon in the 3′ UTR, is an important source of variation in COMT expression in the hippocampus and striatum but not in PFC, suggesting that specific transpositions can alter mRNA isoform use and modulate behavior.

Even if we cannot exclude rs4680 methylation as being at least in part stochastic or inherited, our results suggest that it might mediate the relationship between stress and WM prefrontal activity in Val/Val subjects, thus providing a molecular mechanism for this gene–environment interaction. It is also plausible that rs4680 methylation could mediate other gene–environment interactions involving COMT (Caspi et al., 2005; Kolassa et al., 2010). More generally, our findings indicate that better understanding of the complex interplay among genome, epigenome, and environment can be obtained by studying methylation of CpG sites created or abolished by SNPs or mutations, which during evolution may originate themselves from “mutagenic” CpG methylation (Holliday and Grigg, 1993; Zemach et al., 2010). In fact, DNA methylation is sensitive to environment, and methylated cytosines (5-methylcytosine) can undergo spontaneous deamination to thymine, causing G/C to A/T transition and thus creating genetic variants and abolishing CpG sites that may be involved in regulation of gene expression. At the same time, these genetic variants creating or abolishing a CpG site could alter response to environmental factors that might act via methylation.

In conclusion, our findings suggest that genetic variation at COMT rs4680 and stress interact with and through methylation on prefrontal activity, thus demonstrating how epigenetics may modulate the effect of genetic variation and create a two-way bridge between genes and environment.

Notes

Supplemental material (primer sequences and supplemental results) is available at http://www.psichiat.uniba.it/png/ursinijneurosci2011SM.pdf. This material has not been peer reviewed.

Footnotes

  • We thank Banca Carime S.p.A. for providing an unrestricted grant used to support this work. We thank Dr. Elizabeth Tunbridge for providing total-COMT mRNA primer sequences and for help in designing MB-COMT expression experiments; Riccarda Lomuscio and Dr. Rita Masellis and Dr. Giuseppe Desimeis for help with data acquisition; and Dr. Giovanna Punzi for helpful discussions.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Dr. Alessandro Bertolino, Dipartimento di Scienze Neurologiche e Psichiatriche, Università degli Studi di Bari, Piazza Giulio Cesare, 11, 70124 Bari, Italy. a.bertolino{at}psichiat.uniba.it

References

  1. ↵
    1. Arnsten AF
    (2009) Stress signalling pathways that impair prefrontal cortex structure and function. Nat Rev Neurosci 10:410–422.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Arnsten AF,
    2. Goldman-Rakic PS
    (1998) Noise stress impairs prefrontal cortical cognitive function in monkeys: evidence for a hyperdopaminergic mechanism. Arch Gen Psychiatry 55:362–368.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Bertolino A,
    2. Caforio G,
    3. Blasi G,
    4. De Candia M,
    5. Latorre V,
    6. Petruzzella V,
    7. Altamura M,
    8. Nappi G,
    9. Papa S,
    10. Callicott JH,
    11. Mattay VS,
    12. Bellomo A,
    13. Scarabino T,
    14. Weinberger DR,
    15. Nardini M
    (2004) Interaction of COMT (Val(108/158)Met) genotype and olanzapine treatment on prefrontal cortical function in patients with schizophrenia. Am J Psychiatry 161:1798–1805.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bertolino A,
    2. Rubino V,
    3. Sambataro F,
    4. Blasi G,
    5. Latorre V,
    6. Fazio L,
    7. Caforio G,
    8. Petruzzella V,
    9. Kolachana B,
    10. Hariri A,
    11. Meyer-Lindenberg A,
    12. Nardini M,
    13. Weinberger DR,
    14. Scarabino T
    (2006) Prefrontal–hippocampal coupling during memory processing is modulated by COMT Val158Met genotype. Biol Psychiatry 60:1250–1258.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bertolino A,
    2. Taurisano P,
    3. Pisciotta NM,
    4. Blasi G,
    5. Fazio L,
    6. Romano R,
    7. Gelao B,
    8. Lo Bianco L,
    9. Lozupone M,
    10. Di Giorgio A,
    11. Caforio G,
    12. Sambataro F,
    13. Niccoli-Asabella A,
    14. Papp A,
    15. Ursini G,
    16. Sinibaldi L,
    17. Popolizio T,
    18. Sadee W,
    19. Rubini G
    (2010) Genetically determined measures of striatal D2 signaling predict prefrontal activity during working memory performance. PLoS One 5:e9348.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Bollati V,
    2. Baccarelli A,
    3. Hou L,
    4. Bonzini M,
    5. Fustinoni S,
    6. Cavallo D,
    7. Byun HM,
    8. Jiang J,
    9. Marinelli B,
    10. Pesatori AC,
    11. Bertazzi PA,
    12. Yang AS
    (2007) Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res 67:876–880.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Caspi A,
    2. Sugden K,
    3. Moffitt TE,
    4. Taylor A,
    5. Craig IW,
    6. Harrington H,
    7. McClay J,
    8. Mill J,
    9. Martin J,
    10. Braithwaite A,
    11. Poulton R
    (2003) Influence of life stress on depression: moderation by a polymorphism in the 5-HTT gene. Science 301:386–389.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Caspi A,
    2. Moffitt TE,
    3. Cannon M,
    4. McClay J,
    5. Murray R,
    6. Harrington H,
    7. Taylor A,
    8. Arseneault L,
    9. Williams B,
    10. Braithwaite A,
    11. Poulton R,
    12. Craig IW
    (2005) Moderation of the effect of adolescent-onset cannabis use on adult psychosis by a functional polymorphism in the catechol-O-methyltransferase gene: longitudinal evidence of a gene X environment interaction. Biol Psychiatry 57:1117–1127.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Castri P,
    2. Iacovelli L,
    3. De Blasi A,
    4. Giubilei F,
    5. Moretti A,
    6. Capone FT,
    7. Nicoletti F,
    8. Orzi F
    (2007) Reduced insulin-induced phosphatidylinositol-3-kinase activation in peripheral blood mononuclear leucocytes from patients with Alzheimer's disease. Eur J Neurosci 26:2469–2472.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Chen J,
    2. Lipska BK,
    3. Halim N,
    4. Ma QD,
    5. Matsumoto M,
    6. Melhem S,
    7. Kolachana BS,
    8. Hyde TM,
    9. Herman MM,
    10. Apud J,
    11. Egan MF,
    12. Kleinman JE,
    13. Weinberger DR
    (2004) Functional analysis of genetic variation in catechol-O-methyltransferase (COMT): effects on mRNA, protein, and enzyme activity in postmortem human brain. Am J Hum Genet 75:807–821.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Dempster EL,
    2. Mill J,
    3. Craig IW,
    4. Collier DA
    (2006) The quantification of COMT mRNA in post mortem cerebellum tissue: diagnosis, genotype, methylation and expression. BMC Med Genet 7:10.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Desaulniers D,
    2. Xiao GH,
    3. Leingartner K,
    4. Chu I,
    5. Musicki B,
    6. Tsang BK
    (2005) Comparisons of brain, uterus, and liver mRNA expression for cytochrome p450s, DNA methyltransferase-1, and catechol-O-methyltransferase in prepubertal female Sprague-Dawley rats exposed to a mixture of aryl hydrocarbon receptor agonists. Toxicol Sci 86:175–184.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Egan MF,
    2. Goldberg TE,
    3. Kolachana BS,
    4. Callicott JH,
    5. Mazzanti CM,
    6. Straub RE,
    7. Goldman D,
    8. Weinberger DR
    (2001) Effect of COMT Val108/158 Met genotype on frontal lobe function and risk for schizophrenia. Proc Natl Acad Sci U S A 98:6917–6922.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Fraga MF,
    2. Ballestar E,
    3. Paz MF,
    4. Ropero S,
    5. Setien F,
    6. Ballestar ML,
    7. Heine-Suñer D,
    8. Cigudosa JC,
    9. Urioste M,
    10. Benitez J,
    11. Boix-Chornet M,
    12. Sanchez-Aguilera A,
    13. Ling C,
    14. Carlsson E,
    15. Poulsen P,
    16. Vaag A,
    17. Stephan Z,
    18. Spector TD,
    19. Wu Y-Z,
    20. Plass C,
    21. Esteller M
    (2005) Epigenetic differences arise during the lifetime of monozygotic twins. Proc Natl Acad Sci U S A 102:10604–10609.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Gao MC,
    2. Bellugi U,
    3. Dai L,
    4. Mills DL,
    5. Sobel EM,
    6. Lange K,
    7. Korenberg JR
    (2010) Intelligence in Williams Syndrome is related to STX1A, which encodes a component of the presynaptic SNARE complex. PLoS One 5:e10292.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Gearhart DA,
    2. Middlemore ML,
    3. Terry AV
    (2006) ELISA methods to measure cholinergic markers and nerve growth factor receptors in cortex, hippocampus, prefrontal cortex, and basal forebrain from rat brain. J Neurosci Methods 150:159–173.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Gogos JA,
    2. Morgan M,
    3. Luine V,
    4. Santha M,
    5. Ogawa S,
    6. Pfaff D,
    7. Karayiorgou M
    (1998) Catechol-O-methyltransferase-deficient mice exhibit sexually dimorphic changes in catecholamine levels and behavior. Proc Natl Acad Sci U S A 95:9991–9996.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Habeck C,
    2. Rakitin BC,
    3. Moeller J,
    4. Scarmeas N,
    5. Zarahn E,
    6. Brown T,
    7. Stern Y
    (2005) An event-related fMRI study of the neural networks underlying the encoding, maintenance, and retrieval phase in a delayed-match-to-sample task. Brain Res Cogn Brain Res 23:207–220.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Hains AB,
    2. Arnsten AF
    (2008) Molecular mechanisms of stress-induced prefrontal cortical impairment: implications for mental illness. Learn Mem 15:551–564.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Holliday R,
    2. Grigg GW
    (1993) DNA methylation and mutation. Mutat Res 285:61–67.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Irizarry RA,
    2. Ladd-Acosta C,
    3. Wen B,
    4. Wu Z,
    5. Montano C,
    6. Onyango P,
    7. Cui H,
    8. Gabo K,
    9. Rongione M,
    10. Webster M,
    11. Ji H,
    12. Potash JB,
    13. Sabunciyan S,
    14. Feinberg AP
    (2009) The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet 41:178–186.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Jirtle RL,
    2. Skinner MK
    (2007) Environmental epigenomics and disease susceptibility. Nat Rev Genet 8:253–262.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Kolassa IT,
    2. Kolassa S,
    3. Ertl V,
    4. Papassotiropoulos A,
    5. De Quervain DJ
    (2010) The risk of posttraumatic stress disorder after trauma depends on traumatic load and the catechol-O-methyltransferase Val(158)Met polymorphism. Biol Psychiatry 67:304–308.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Krugel LK,
    2. Biele G,
    3. Mohr PN,
    4. Li SC,
    5. Heekeren HR
    (2009) Genetic variation in dopaminergic neuromodulation influences the ability to rapidly and flexibly adapt decisions. Proc Natl Acad Sci U S A 106:17951–17956.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Li Z,
    2. Mulligan MK,
    3. Wang X,
    4. Miles MF,
    5. Lu L,
    6. Williams RW
    (2010) A transposon in Comt generates mRNA variants and causes widespread expression and behavioral differences among mice. PLoS One 5:e12181.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Lotta T,
    2. Vidgren J,
    3. Tilgmann C,
    4. Ulmanen I,
    5. Melén K,
    6. Julkunen I,
    7. Taskinen J
    (1995) Kinetics of human soluble and membrane-bound catechol-O-methyltransferase: a revised mechanism and description of the thermolabile variant of the enzyme. Biochemistry 34:4202–4210.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Matsumoto M,
    2. Weickert CS,
    3. Beltaifa S,
    4. Kolachana B,
    5. Chen J,
    6. Hyde TM,
    7. Herman MM,
    8. Weinberger DR,
    9. Kleinman JE
    (2003) Catechol-O-methyltransferase (COMT) mRNA expression in the dorsolateral prefrontal cortex of patients with schizophrenia. Neuropsychopharmacology 28:1521–1530.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Mattay VS,
    2. Goldberg TE,
    3. Fera F,
    4. Hariri AR,
    5. Tessitore A,
    6. Egan MF,
    7. Kolachana B,
    8. Callicott JH,
    9. Weinberger DR
    (2003) Catechol O-methyltransferase val158-met genotype and individual variation in the brain response to amphetamine. Proc Natl Acad Sci U S A 100:6186–6191.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Maunakea AK,
    2. Nagarajan RP,
    3. Bilenky M,
    4. Ballinger TJ,
    5. D'Souza C,
    6. Fouse SD,
    7. Johnson BE,
    8. Hong C,
    9. Nielsen C,
    10. Zhao Y,
    11. Turecki G,
    12. Delaney A,
    13. Varhol R,
    14. Thiessen N,
    15. Shchors K,
    16. Heine VM,
    17. Rowitch DH,
    18. Xing X,
    19. Fiore C,
    20. Schillebeeckx M,
    21. et al.
    (2010) Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature 466:253–257.
    OpenUrlCrossRefPubMed
  30. ↵
    1. McGowan PO,
    2. Sasaki A,
    3. D'Alessio AC,
    4. Dymov S,
    5. Labonté B,
    6. Szyf M,
    7. Turecki G,
    8. Meaney MJ
    (2009) Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat Neurosci 12:342–348.
    OpenUrlCrossRefPubMed
  31. ↵
    1. McIntosh AM,
    2. Holmes S,
    3. Gleeson S,
    4. Burns JK,
    5. Hodges AK,
    6. Byrne MM,
    7. Dobbie R,
    8. Miller P,
    9. Lawrie SM,
    10. Johnstone EC
    (2002) Maternal recall bias, obstetric history and schizophrenia. Br J Psychiatry 181:520–525.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. McNeil TF,
    2. Cantor-Graae E,
    3. Torrey EF,
    4. Sjöström K,
    5. Bowler A,
    6. Taylor E,
    7. Rawlings R,
    8. Higgins ES
    (1994) Obstetric complications in histories of monozygotic twins discordant and concordant for schizophrenia. Acta Psychiatr Scand 89:196–204.
    OpenUrlPubMed
  33. ↵
    1. Meaney MJ,
    2. Ferguson-Smith AC
    (2010) Epigenetic regulation of the neural transcriptome: the meaning of the marks. Nat Neurosci 13:1313–1318.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Meyer-Lindenberg A,
    2. Nichols T,
    3. Callicott JH,
    4. Ding J,
    5. Kolachana B,
    6. Buckholtz J,
    7. Mattay VS,
    8. Egan M,
    9. Weinberger DR
    (2006) Impact of complex genetic variation in COMT on human brain function. Mol Psychiatry 11:867–877.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Mier D,
    2. Kirsch P,
    3. Meyer-Lindenberg A
    (2010) Neural substrates of pleiotropic action of genetic variation in COMT: a meta-analysis. Mol Psychiatry 15:918–927.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Mill J,
    2. Dempster E,
    3. Caspi A,
    4. Williams B,
    5. Moffitt T,
    6. Craig I
    (2006) Evidence for monozygotic twin (MZ) discordance in methylation level at two CpG sites in the promoter region of the catechol-O-methyltransferase (COMT) gene. Am J Med Genet B Neuropsychiatr Genet 141B:421–425.
    OpenUrlPubMed
  37. ↵
    1. Miller CA,
    2. Sweatt JD
    (2007) Covalent modification of DNA regulates memory formation. Neuron 53:857–869.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Miller CA,
    2. Gavin CF,
    3. White JA,
    4. Parrish RR,
    5. Honasoge A,
    6. Yancey CR,
    7. Rivera IM,
    8. Rubio MD,
    9. Rumbaugh G,
    10. Sweatt JD
    (2010) Cortical DNA methylation maintains remote memory. Nat Neurosci 13:664–666.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Murgatroyd C,
    2. Patchev AV,
    3. Wu Y,
    4. Micale V,
    5. Bockmühl Y,
    6. Fischer D,
    7. Holsboer F,
    8. Wotjak CT,
    9. Almeida OF,
    10. Spengler D
    (2009) Dynamic DNA methylation programs persistent adverse effects of early-life stress. Nat Neurosci 12:1559–1566.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Murphy BC,
    2. O'Reilly RL,
    3. Singh SM
    (2005) Site-specific cytosine methylation in S-COMT promoter in 31 brain regions with implications for studies involving schizophrenia. Am J Med Genet B Neuropsychiatr Genet 133B:37–42.
    OpenUrlPubMed
  41. ↵
    1. Nackley AG,
    2. Shabalina SA,
    3. Tchivileva IE,
    4. Satterfield K,
    5. Korchynskyi O,
    6. Makarov SS,
    7. Maixner W,
    8. Diatchenko L
    (2006) Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science 314:1930–1933.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Oesch-Bartlomowicz B,
    2. Huelster A,
    3. Wiss O,
    4. Antoniou-Lipfert P,
    5. Dietrich C,
    6. Arand M,
    7. Weiss C,
    8. Bockamp E,
    9. Oesch F
    (2005) Aryl hydrocarbon receptor activation by cAMP vs. dioxin: divergent signaling pathways. Proc Natl Acad Sci U S A 102:9218–9223.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Pani L,
    2. Porcella A,
    3. Gessa GL
    (2000) The role of stress in the pathophysiology of the dopaminergic system. Mol Psychiatry 5:14–21.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Papaleo F,
    2. Crawley JN,
    3. Song J,
    4. Lipska BK,
    5. Pickel J,
    6. Weinberger DR,
    7. Chen J
    (2008) Genetic dissection of the role of catechol-O-methyltransferase in cognition and stress reactivity in mice. J Neurosci 28:8709–8723.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Preti A,
    2. Cardascia L,
    3. Zen T,
    4. Marchetti M,
    5. Favaretto G,
    6. Miotto P
    (2000) Risk for obstetric complications and schizophrenia. Psychiatry Res 96:127–139.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Sartor MA,
    2. Schnekenburger M,
    3. Marlowe JL,
    4. Reichard JF,
    5. Wang Y,
    6. Fan Y,
    7. Ma C,
    8. Karyala S,
    9. Halbleib D,
    10. Liu X,
    11. Medvedovic M,
    12. Puga A
    (2009) Genomewide analysis of aryl hydrocarbon receptor binding targets reveals an extensive array of gene clusters that control morphogenetic and developmental programs. Environ Health Perspect 117:1139–1146.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Schnekenburger M,
    2. Peng L,
    3. Puga A
    (2007) HDAC1 bound to the Cyp1a1 promoter blocks histone acetylation associated with Ah receptor-mediated trans-activation. Biochim Biophys Acta 1769:569–578.
    OpenUrlPubMed
  48. ↵
    1. Sesack SR,
    2. Hawrylak VA,
    3. Matus C,
    4. Guido MA,
    5. Levey AI
    (1998) Dopamine axon varicosities in the prelimbic division of the rat prefrontal cortex exhibit sparse immunoreactivity for the dopamine transporter. J Neurosci 18:2697–2708.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Smolka MN,
    2. Schumann G,
    3. Wrase J,
    4. Grüsser SM,
    5. Flor H,
    6. Mann K,
    7. Braus DF,
    8. Goldman D,
    9. Büchel C,
    10. Heinz A
    (2005) Catechol-O-methyltransferase val158met genotype affects processing of emotional stimuli in the amygdala and prefrontal cortex. J Neurosci 25:836–842.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Szyf M,
    2. McGowan P,
    3. Meaney MJ
    (2008) The social environment and the epigenome. Environ Mol Mutagen 49:46–60.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Tan HY,
    2. Chen Q,
    3. Goldberg TE,
    4. Mattay VS,
    5. Meyer-Lindenberg A,
    6. Weinberger DR,
    7. Callicott JH
    (2007) Catechol-O-methyltransferase Val158Met modulation of prefrontal–parietal–striatal brain systems during arithmetic and temporal transformations in working memory. J Neurosci 27:13393–13401.
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Tost J,
    2. Gut IG
    (2007) DNA methylation analysis by pyrosequencing. Nat Protoc 2:2265–2275.
    OpenUrlCrossRefPubMed
  53. ↵
    1. van Winkel R,
    2. Henquet C,
    3. Rosa A,
    4. Papiol S,
    5. Fananás L,
    6. De Hert M,
    7. Peuskens J,
    8. van Os J,
    9. Myin-Germeys I
    (2008) Evidence that the COMT(Val158Met) polymorphism moderates sensitivity to stress in psychosis: an experience-sampling study. Am J Med Genet B Neuropsychiatr Genet 147B:10–17.
    OpenUrlPubMed
  54. ↵
    1. Verdoux H,
    2. Sutter AL,
    3. Glatigny-Dallay E,
    4. Minisini A
    (2002) Obstetrical complications and the development of postpartum depressive symptoms: a prospective survey of the MATQUID cohort. Acta Psychiatr Scand 106:212–219.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Wager TD,
    2. Smith EE
    (2003) Neuroimaging studies of working memory: a meta-analysis. Cogn Affect Behav Neurosci 3:255–274.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Walder DJ,
    2. Trotman HD,
    3. Cubells JF,
    4. Brasfield J,
    5. Tang YL,
    6. Walker EF
    (2010) Catechol-O-methyltransferase modulation of cortisol secretion in psychiatrically at-risk and healthy adolescents. Psychiatr Genet 20:166–170.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Wong CC,
    2. Caspi A,
    3. Williams B,
    4. Craig IW,
    5. Houts R,
    6. Ambler A,
    7. Moffitt TE,
    8. Mill J
    (2010) A longitudinal study of epigenetic variation in twins. Epigenetics 5:516–526.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Yacubian J,
    2. Sommer T,
    3. Schroeder K,
    4. Gläscher J,
    5. Kalisch R,
    6. Leuenberger B,
    7. Braus DF,
    8. Büchel C
    (2007) Gene-gene interaction associated with neural reward sensitivity. Proc Natl Acad Sci U S A 104:8125–8130.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Zemach A,
    2. McDaniel IE,
    3. Silva P,
    4. Zilberman D
    (2010) Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science 328:916–919.
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

The Journal of Neuroscience: 31 (18)
Journal of Neuroscience
Vol. 31, Issue 18
4 May 2011
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Stress-Related Methylation of the Catechol-O-Methyltransferase Val158 Allele Predicts Human Prefrontal Cognition and Activity
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Citation Tools
Stress-Related Methylation of the Catechol-O-Methyltransferase Val158 Allele Predicts Human Prefrontal Cognition and Activity
Gianluca Ursini, Valentina Bollati, Leonardo Fazio, Annamaria Porcelli, Luisa Iacovelli, Assia Catalani, Lorenzo Sinibaldi, Barbara Gelao, Raffaella Romano, Antonio Rampino, Paolo Taurisano, Marina Mancini, Annabella Di Giorgio, Teresa Popolizio, Andrea Baccarelli, Antonio De Blasi, Giuseppe Blasi, Alessandro Bertolino
Journal of Neuroscience 4 May 2011, 31 (18) 6692-6698; DOI: 10.1523/JNEUROSCI.6631-10.2011

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Stress-Related Methylation of the Catechol-O-Methyltransferase Val158 Allele Predicts Human Prefrontal Cognition and Activity
Gianluca Ursini, Valentina Bollati, Leonardo Fazio, Annamaria Porcelli, Luisa Iacovelli, Assia Catalani, Lorenzo Sinibaldi, Barbara Gelao, Raffaella Romano, Antonio Rampino, Paolo Taurisano, Marina Mancini, Annabella Di Giorgio, Teresa Popolizio, Andrea Baccarelli, Antonio De Blasi, Giuseppe Blasi, Alessandro Bertolino
Journal of Neuroscience 4 May 2011, 31 (18) 6692-6698; DOI: 10.1523/JNEUROSCI.6631-10.2011
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Notes
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Responses to this article

Respond to this article

Jump to comment:

  • Reply to Rousselet and Pernet: Use of inferential statistics to reject the null hypothesis and the implications of COMT methylation
    Alessandro Bertolino
    Published on: 30 November 2011
  • Robust statistics show no evidence for a relationship between methylation and several variables
    Guillaume A. Rousselet
    Published on: 19 September 2011
  • Published on: (30 November 2011)
    Reply to Rousselet and Pernet: Use of inferential statistics to reject the null hypothesis and the implications of COMT methylation
    • Alessandro Bertolino, Professor
    • Other Contributors:
      • Gianluca Ursini, Giuseppe Blasi

    Using a translational approach, we evaluated in living humans whether stress-related methylation of the COMT Val158 allele is associated with working memory (WM) performance, with related prefrontal cortex (PFC) activity, and with MB-COMT expression and protein levels (Ursini et al., 2011). Our research hypothesis was based on a vast literature on the effects of stress on prefrontal dopamine, on the...

    Show More

    Using a translational approach, we evaluated in living humans whether stress-related methylation of the COMT Val158 allele is associated with working memory (WM) performance, with related prefrontal cortex (PFC) activity, and with MB-COMT expression and protein levels (Ursini et al., 2011). Our research hypothesis was based on a vast literature on the effects of stress on prefrontal dopamine, on the effects of COMT on prefrontal dopamine, WM behavior and related activity, as well as on recent studies suggesting that intragenic methylation is important for gene function. In particular, the findings we reported are biologically predictable, internally consistent, and coherent in demonstrating how the physiological effects of the Val allele (presumably through modulation of dopamine levels) on memory performance and brain activity can be partially modulated by its methylation. We also demonstrated in a sample of rats a positive correlation between COMT methylation in PBMCs and in PFC, further supporting the biological plausibility of the relationship between COMT methylation in PBMCs and brain phenotypes reported in humans. In their comment, Drs. Rousselet and Pernet have not taken into the necessary consideration the strong biological background that led to these experiments, avoiding any biological discussion of our data and overlooking other consistent and important results reported in the paper.

    We believe it is important to underline that inferential statistics make propositions about populations based on a research hypothesis. In frequentist inference, testing should be used to reject a null hypothesis, thus affirming the alternative hypothesis, which is the research hypothesis. Moreover, when data from other independent studies and experiments are available, biological statistical testing should consider not only a single analysis, but also the biological plausibility as suggested by the consistency of the results with other earlier reports and with other results of the actual study. Furthermore, replication in independent studies is always necessary. More specifically, Rousselet and Pernet are concerned about the possibility that some of the results may not be as statistically significant as reported in the paper because of the presence of potential outliers. We firmly believe that the letter by Drs. Rousselet and Pernet may be an important stimulus for a reflection on why and how to identify outliers and, more generally, on the risk of type I and type II errors in biological sciences. We also firmly believe that the quality of a study should not be judged solely on how low a p-value is, a problem too often overstated. For example, when analyzing significant correlations, a more important indicator of biological significance is the proportion of variance explained (coefficient of determination) which provides prediction of future outcomes based on the current model.

    An outlier is defined as an observation that "appears" to be inconsistent with other observations in the data set (Barnett and Lewis, 1985). An outlier has a low probability that it originates from the same statistical distribution as the other observations in the data set. As Drs. Schott and Duzel (Schott and Duzel, 2011) underline in their response to a similar letter by Drs. Rousselet and Pernet (Rousselet and Pernet, 2011), the rationale for outlier exclusion is the assumption that they represent artifactual results from experimental errors, such as sampling or measurement problems. Keeping in mind the purpose of outlier detection, the first rule of thumb to identify potential outliers is to visually inspect the distribution of the data and evaluate whether there is any data point presumably not originating from a similar distribution as the others. Once a potential outlier in the univariate or bivariate space is identified, the second most obvious step is to go back to the data and identify potential sources of experimental error. In fact, we had adopted an experimental procedure that was also designed to reduce the risk of artifactual results. Specifically, the methylation analysis was performed in replicate, twice in separate plates, and mean values were used in the analysis. The coefficient of variation (SD/mean) was very low (1.2%), thus confirming very low variability of these results. Similarly, we performed mRNA expression measurements in triplicate and, as is commonly done, the analysis was repeated for samples in which the standard deviation of theshold cycle between replicates was > 0.25. Therefore, the SD between replicates was <0.25 for all subjects included in the analysis, in which mean values of the replicates were used. Moreover, we applied Grubb's test to WM accuracy data in order to detect outliers using the mean and SD from a large sample of healthy subjects that is more representative of the general population compared with the sample in our paper (see below, N=584; mean=78.3; SD=28.6).

    All of this having been said, in the sample included in our paper there was only one potential univariate outlier for methylation and this subject was not excluded for the following reasons. After the identification of this human subject with lower methylation (76.9%), we went back and checked the data relative to this subject for whom % methylation was 76.8 in the first plate and 76.9 in the second plate. Moreover, lack of experimental errors was also supported by the lack of significant deviations of % methylation from the general distribution of subjects from the same plates. Finally, we have expanded our sample to include 86 Val/Val healthy individuals and this subject does not result as a potential outlier as indicated by Grubb's Test using mean and SD from this larger population (N=86; mean=92.4; SD=5.4). Similarly, the same subject also had fairly low WM performance (61.9 % of correct responses). Considering our database of 584 healthy subjects which is more representative of the general population (N=584; mean % correct responses=78.3; SD=28.6), this subject cannot be considered as an outlier in terms of behavioral performance. Interestingly, the same subject is a potential outlier for MB-COMT mRNA expression. As we had done for the methylation analysis, we went back to our experiments and found no specific technical reason to exclude this subject, including the fact that the triplicate results were very similar. Furthermore, the relationship between COMT methylation, expression, and WM accuracy in this subject was consistent with the distribution of the rest of the data, with the general relationship between DNA methylation and expression known in the literature, and with the predictable relationship between COMT methylation and WM accuracy, thus suggesting that these values result from biological variation. These data points may be located in the right tail of a normal distribution, rather than being outliers.

    In our paper, after screening for potential outliers, we used Spearman correlation to test our hypotheses on the relationships between data variables. This type of correlation uses variable ranks instead of their values, calculating a non-parametric correlation coefficient. It should be noted that Spearman's correlation coefficient is used to estimate a relation between variables independently of their distribution., because distribution-specific information is lost when the data are replaced by their ranks. This coefficient is also less sensitive to outliers, which is important when working with experimental data. Using Spearman correlation, we analyzed whether COMT methylation in PBMCs was related to WM performance, MB-COMT expression and COMT methylation in PFC. After removing what they report as being bivariate outliers, Rousselet and Pernet found that some of the results in figures 1B, 3, and 4 are, in their opinion, not significant. It is important to note that they detected outliers with a single statistical method based on median absolute deviation (MAD). However, MAD is also known to underestimate the standard deviation in relatively small sample sizes (Rousseeuw and Croux, 1993; Motulsky and Brown, 2006). Other statistical methods for detection of outliers exist and they may lead to different conclusions. ROUT (www.graphpad.com) is another such method that combines robust regression and bivariate outlier detection (Motulsky and Brown, 2006). Notably, the application of this alternative method can lead to different results, compared with the analyses by Rousselet and Pernet. For example, application of the ROUT method to our data presented in Figure 1B did not result in detection of any outlier and the results indicated a statistically significant relationship (N=19; R=0.25; p=0.02). Below, we provide more details related to different statistical methods to detect outliers in bivariate space relative to all findings criticized by Rousselet and Pernet. In our opinion, all these additional analyses strongly make the point that outlier detection and removal is ultimately a choice of the investigators which should be based on test-retest of the methods and on hypotheses about the biological plausibility/implausibility of the data-points. In other words, outlier detection based on only one statistical method without carefully considering technical reproducibility or biological diversity may result in somewhat arbitrary decisions on the data. Therefore, we believe that outlier detection should take into account the possibilities that differences in data points may arise from biological diversity and from experimental mistakes. To estimate this possibility it is useful to consider: the experimental design of a study, the distribution of values within a population more representative than the sample studied, and the relationship between the studied variables on the basis of literature and in the context of the research.



    Tab.1. Bivariate outlier detection based on different statistical methods leads to different results. Only some of the used statistical methods identify one potential bivariate outlier in the analyses concerning Figure 1B and 3 (the same subject with 76.9% methylation who is a potential univariate outlier). After removal of this potential outlier, the correlation between COMT methylation and WM accuracy remains significant (N=18; Spearman Rho=0.54; p=0.018), while the correlation between COMT methylation and mRNA expression becomes a weak statistical trend even though around 20% of the variance is still explained (N=13; Spearman Rho= -0.46; p=0.11). However, we believe it would be inappropriate to exclude this subject from the analysis, since this data point is unlikely to be artifactual as detailed above in the text. Moreover, the relationship between methylation and accuracy for this subject is consistent with the distribution of the rest of the data and, notably, is consistent with our fMRI results (figure 2), which were not considered by Rousselet and Pernet. Similarly, the relationship between methylation and expression for this subject is consistent with the distribution of the rest of the data and with the general relationship between DNA methylation and expression known in the literature. Also, to further support and extend the results of the relationship between methylation and mRNA expression, we performed a Western Blotting analysis, which was not considered by Rousselet and Pernet, reporting a biologically consistent and plausible relationship between COMT methylation and MB-COMT protein expression (N=7; Rho= -0.78; p=0.03). Even in the analysis concerning Figure 4, only the Cook's distance method identifies two potential bivariate outliers. Removing these two potential outliers from the analyses (subjects with PBMCs and PFC % methylation respectively equal to 57.5 and 65.2; 62.1 and 74.3), the correlation remains significant and the proportion of variance explained is about 30% (N=14; Rho=0.55; p=0.04). It is interesting to note that even using the skipped correlation approach as applied by Rousselet and Perned the proportion of variance explained is greater than what originally reported in the paper. However, similar to the previous analyses, we believe it would be inappropriate to consider these data points as outliers, given that the methylation analysis is performed in replicate and there are no univariate outliers for PFC and PBMCs methylation. *All the analyses, except ROUT method, were performed with Statistica (www.statsoft.com).

    Conclusions
    In summary, the relationship between stress, methylation and prefrontal cognition that we demonstrated in our paper remains robust after the use of different statistical techniques to account for possible outliers. This relationship is particularly strong if we consider the biologically predictable convergence and internal consistency of our results. With all these considerations, we advocate that robustness of the hypothesis and of the methodological acquisition of data have to be considered first and foremost when judging the biological plausibility of any result.


    REFERENCES

    Barnett V and Lewis T (1985). Outliers in Statistical Data (2nd Ed.). New York, NY: John Wiley & Sons.

    Motulsky HJ, Brown RE (2006) Detecting outliers when fitting data with nonlinear regression - a new method based on robust nonlinear regression and the false discovery rate. BMC Bioinformatics 7:123.

    Rousseeuw PJ, Croux C (1993) Alternatives to the median absolute deviation. J Amer Stat Assoc 88:1273-1283.

    Rousselet GA, Pernet CR (2011) Robust statistics show no evidence for a relationship between fiber density and memory performance. Proc Natl Acad Sci U S A 108:E598.

    Schott BH and Duzel ED (2011) Reply to Rousselet et al.: A robust relationship between fiber density and memory recall. Proc Natl Acad Sci U S A 108:E599.

    Ursini G, Bollati V, Fazio L, Porcelli A, Iacovelli L, Catalani A, Sinibaldi L, Gelao B, Romano R, Rampino A, Taurisano P, Mancini M, Di Giorgio A, Popolizio T, Baccarelli A, De Blasi A, Blasi G, Bertolino A (2011) Stress-related methylation of the catechol-O-methyltransferase Val 158 allele predicts human prefrontal cognition and activity. J Neurosci 31:6692-6698.

    Conflict of Interest:

    None declared

    Show Less
    Competing Interests: None declared.
  • Published on: (19 September 2011)
    Robust statistics show no evidence for a relationship between methylation and several variables
    • Guillaume A. Rousselet, Senior Lecturer
    • Other Contributors:
      • Cyril R. Pernet

    Ursini et al. report significant Spearman correlations in their Figures 1, 3 and 4. Visual inspection of these figures suggests that most of the points are clustered together and show no obvious relationship between variables. To confirm this impression, we extracted data from the figures using the Mac software GraphClick, version 3.0 (Arizona Software, 2008).

    Using these extracted data points, we could replicate the auth...

    Show More

    Ursini et al. report significant Spearman correlations in their Figures 1, 3 and 4. Visual inspection of these figures suggests that most of the points are clustered together and show no obvious relationship between variables. To confirm this impression, we extracted data from the figures using the Mac software GraphClick, version 3.0 (Arizona Software, 2008).

    Using these extracted data points, we could replicate the authors' original correlation values. However, using a robust statistical technique called a skipped correlation, whereby the bivariate outliers are removed and Spearman correlation is applied to the remaining data points with adjustments for pruning the data (Wilcox, 2004, 2005), none of the correlations are significant. Outliers were detected in the bivariate space using the MAD (median absolute deviation) rule on the data points orthogonally projected onto the lines linking these same data point to the center of the cloud that minimized the covariance determinant (i.e., the scatter, Rousseeuw, 1984; Rousseeuw & Van Driessen, 1999).

    Thus, the authors' conclusions that rely on these correlations appear to be unsubstantiated. These results provide a good illustration of a problem too often overlooked: Spearman correlation is sensitive to bivariate outliers.

    Figure 1B:

    Authors: Spearman r=0.59, p=0.007

    Our analyses: Spearman r=0.591, p=0.0076

    Skipped correlation: rp=0.392, t=1.756, tcrit=2.685

    For the skipped correlation, the observed t value (t) must be as large or larger than the critical t value (tcrit) for the effect to be considered significant assuming a type I error rate of 0.05.

    Figure 3:

    Authors: Spearman r=-0.57, p=0.03

    Our analyses: Spearman r=-0.577, p=0.031

    Skipped correlation: rp=-0.326, t=1.196, tcrit=2.816

    Figure 4:

    Authors: Spearman r=0.49, p=0.05

    Our analyses: Spearman r=0.50, p=0.048

    Skipped correlation: rp=0.547, t=2.446, tcrit=2.754

    References

    Rousseeuw, P. J. (1984). Least Median of Squares Regression. Journal of the American Statistical Association, 79(388), 871-880.

    Rousseeuw, P. J., & Van Driessen, K. (1999). A fast algorithm for the minimum covariance determinant estimator. Technometrics, 41(3), 212-223.

    Wilcox, R. (2004). Inferences based on a skipped correlation coefficient. Journal of Applied Statistics, 31(2), 131-143.

    Wilcox, R. R. (2005). Introduction to Robust Estimation and Hypothesis Testing (2nd Ed.). New York, NY: Elsevier Academic Press.

    Conflict of Interest:

    None declared

    Show Less
    Competing Interests: None declared.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Articles

Behavioral/Systems/Cognitive

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Behavioral/Systems/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(JNeurosci logo)
(SfN logo)

Copyright © 2021 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.